Genas_system
Eugeniy Mikhailov and Simon Rochester
License GPL.
Solving 4 level atom in 0->1 configuration,
with field propagation along spatial axis Z
no Doppler broadening
Fields:
light field, default circularly polarized
rf field along x
static B-field along z
M=0 upper state sublevel can be shifted independently
We are solving
dE/dz+(1/c)*dE/dt=i*eta*rho_ji, where j level is higher then i.
Note that E is actually a Rabi frequency of electromagnetic field not the EM field
in xmds terms it looks like
dE_dz = i*eta*rhoji - 1/c*L[E], here we moved t dependence to Fourier space
VERY IMPORTANT: all Rabi frequency should be given in [1/s], if you want to
normalize it to something else look drho/dt equation.
No need to renormalizes eta as long as its express through i
the upper level decay rate in the same units as Rabi frequency.
z
Ex Ey
r11 r22 r33 r12 r13 r14 r23 r24 r34 r44
rf
200 200
density_matrix
E_field rfField
Lt
E_field
density_matrix