Interactive Tutorials

# Import required libraries
import numpy as np
import matplotlib.pyplot as plt
from pcm import JMAKModel, PCMDevice

# Initialize the PCM device model
device = PCMDevice()

# Check default parameters
print(f"Amorphous resistance: {device.R_amorphous} Ω")
print(f"Crystalline resistance: {device.R_crystalline} Ω")

# Create a simple voltage pulse
t = np.linspace(0, 100e-9, 1000)  # 100 ns time window
V = np.zeros_like(t)
V[(t > 10e-9) & (t < 60e-9)] = 2.0  # 50 ns pulse at 2V

# Simulate device response
result = device.simulate_voltage_pulse(V, t)

# Plot results
fig, axes = plt.subplots(2, 2, figsize=(12, 8))

axes[0,0].plot(t*1e9, V)
axes[0,0].set_xlabel('Time (ns)')
axes[0,0].set_ylabel('Voltage (V)')

axes[0,1].plot(t*1e9, result['temperature'])
axes[0,1].set_xlabel('Time (ns)')
axes[0,1].set_ylabel('Temperature (K)')

axes[1,0].semilogy(t*1e9, result['resistance'])
axes[1,0].set_xlabel('Time (ns)')
axes[1,0].set_ylabel('Resistance (Ω)')

axes[1,1].plot(t*1e9, result['phase'])
axes[1,1].set_xlabel('Time (ns)')
axes[1,1].set_ylabel('Crystalline Fraction')

plt.tight_layout()
plt.show()

# JMAK equation: X(t) = 1 - exp(-(kt)^n)
# where k = k0 * exp(-Ea/kT)

from pcm import JMAKModel

# Initialize with custom parameters
model = JMAKModel(params={
    'Ea': 1.8,   # Activation energy (eV)
    'k0': 1e16,  # Pre-exponential factor (1/s)
    'n': 3.0     # Avrami exponent
})

# Simulate at different temperatures
temperatures = [400, 450, 500, 550, 600]  # Kelvin
time = np.logspace(-9, -3, 100)  # 1 ns to 1 ms

plt.figure(figsize=(10, 6))

for T in temperatures:
    X = model.simulate_isothermal(T, time)
    plt.semilogx(time * 1e6, X, label=f'{T} K')

plt.xlabel('Time (μs)')
plt.ylabel('Crystalline Fraction')
plt.legend()
plt.grid(True, alpha=0.3)
plt.title('Isothermal Crystallization Kinetics')
plt.show()