haeffnerlab / ionsim.jl-examples Goto Github PK

Multi-mode interactions notebook: Laser polarization is not orthogonal to wavevector.

L.k = (x̂ + ẑ)/√2
L.ϵ = ẑ

I had to add these to run the notebook

using Pkg
Pkg.add("StochasticDiffEq")
Pkg.add("DSP")

Currently links back to examples.ionsim.org

ψi_ion = dm(C["S-1/2"])
mode.N = 50
ψi_mode = thermalstate(mode, 10)
ψi = ψi_ion ⊗ ψi_mode
tspan = 0:0.01:50
h = hamiltonian(T, timescale=1e-6, rwa_cutoff=Inf, time_dependent_eta=false, displacement="truncated")
tout, sol = timeevolution.schroedinger_dynamic(tspan, ψi, h);

I guess mode.N should be 50 to match the numerical and analytical plots. mode.N = 100 does not make both plots identical.

tspan was missing. I added tspan.