Initializing Python Environment...
(Downloading NumPy/SciPy, please wait)
Model
$$ \ddot{Q} + 2\gamma \dot{Q} + \omega_0^2 Q + a_4 Q^3 = E_{\text{ext}}(t) $$ $$ E_{\text{ext}}(t) = E_0 \exp\left(-\frac{t^2}{t_0^2}\right) \sin(\omega_{\text{THz}} t) $$
Parameters
Phonon freq. (THz)
1.0
Damping
0.6
Anharmonicity a4
20
Driving freq. (THz)
1.0
Field strength E0
80
Pulse duration (ps)
2.0
Solving 1D ODE...
Calculate 2D Spectrum
Calculating 2D Spectrum...
(It might take 20 ~ 40 s)