ColdAtoms.jl Documentation
ColdAtoms.release_evolve
— Methodrelease_evolve(tspan, cord, atom_params, trap_params; eps=1e-3)
Simulates evolution of atoms after turning the trap off and recapture results
Input
tspan
– vector specifying the points of time for which output should be displayedcord
– vector of initial coordiantes and velocities $[x,y,z,v_{x},v_{y},v_{z}]$atom_params
– vector [atom mass in a.u., atom temperature in $\mu K$]trap_params
– vector [trap depth $U_{0}$ in $\mu K$, beam waist radius in $\mu m$, beam Rayleigh length in $\mu m$]eps
– (optional, default:1e-3
) cutoff to regularize Metropolis sampler, atoms that have energy over $U_{0}(1-eps)$ are considered to be out of trap
Output
Binary list of length tspan, which contains 1 at idx $i$ if atom is recaptured at tspan[i]
and 0 otherwise
ColdAtoms.release_recapture
— Methodrelease_recapture(tspan, trap_params, atom_params, N; freq=10, skip=1000, eps=1e-3, harmonic=true)
Simulate release and recapture experiment to estimate atom's temperature.
Input
tspan
– vector specifying the points of time for which output should be displayedtrap_params
– vector [trap depth $U_{0}$ in $\mu K$, beam waist radius in $\mu m$, beam Rayleigh length in $\mu m$]atom_params
– vector [atom mass in a.u., atom temperature in $\mu K$]N
– number of Monte-Carlo samples, the same as number of atomsfreq
– (optional, default:10
) number of Metropolis steps skipped between samples to reduce sample dependencyskip
– (optional, default:1000
) number of Metropolis steps skipped before the Markov Chain is considered to reach stationary distributioneps
– (optional, default:1e-3
) cutoff to regularize Metropolis sampler, atoms that have energy over $U_{0}(1-eps)$ are considered to be out of trapharmonic
– (optional, default:true
) uses harmonic approximation of gaussian beam if set totrue
, otherwise uses Metropolis sampler
Output
List of recapture probabilities corresponding to times in tspan
and acceptance rate of Metropolis algorithm. If harmonic
is set to true
, acceptance rate is set to 1.0
ColdAtoms.simulation
— Methodsimulation(
tspan, ψ0,
atom_params,
trap_params,
samples,
f,
red_laser_phase_amplitudes,
blue_laser_phase_amplitudes,
red_laser_params,
blue_laser_params,
detuning_params,
decay_params;
atom_motion=true,
free_motion=true,
laser_noise=true,
spontaneous_decay=true,
parallel=false
)
Simulate two-photon Rydberg excitation of single atom with several sources of decoherence
Input
tspan
– vector specifying the points of time for which output should be displayedψ0
– initial wavefunction vector of normalized complex amplitudes $[c_{g}, c_{p}, c_{r}]$atom_params
– vector [atom mass in a.u., atom temperature in $\mu K$]trap_params
– vector [trap depth $U_{0}$ in $\mu K$, beam waist radius in $\mu m$, beam Rayleigh length in $\mu m$]samples
– Monte-Carlo samples of initial atom coordinates and velocities, can be received using samples_generatef
– frequencies at which laser phase noise is sampledred_laser_phase_amplitudes
– amplitudes of red laser phase noise for correspoding frequenciesf
blue_laser_phase_amplitudes
– amplitudes of blue laser phase noise for correspoding frequenciesf
red_laser_params
– write explanationblue_laser_params
– write explanationdetuning_params
– vector [Δ0, δ0], which sets detuning from intermediate level and Rydberg level correspondinglydecay_params
– write explanationatom_motion
– (optional, default:true
) if set to true, atom motion is includedfree_motion
– (optional, default:true
) if set to true, trap is turned offlaser_noise
– (optional, default:true
) if set to true, laser phase noise is includedspontaneous_decay
– (optional, default:true
) if set to true, spontaneous decay from intermediate level is includedparallel
– (optional, default:false
) parallel implementation is under developmentn
– (optional, default:1
) super-gauss parameter
Output
Outputs Monte-Carlo averaged density matrix and squared density matrix for error calculation with elements ordered in correspondence with order ground, intermediate, Rydberg