diff options
Diffstat (limited to 'liouville.m')
| -rw-r--r-- | liouville.m | 23 |
1 files changed, 20 insertions, 3 deletions
diff --git a/liouville.m b/liouville.m index 736f0d2..63badba 100644 --- a/liouville.m +++ b/liouville.m @@ -54,9 +54,26 @@ atom_properties.dipole_elements=dipole_elements; E_field_drive =[0, 0 , Ed, 0 , Edc, 0, 0 ]; E_field_probe =[0, Ep, 0 , Epc, 0 , 0, 0 ]; E_field_zero =[0, 0 , 0 , 0 , 0 , 0, 0 ]; -E_field.linear = E_field_zero + (0.00000+0.00000i)*E_field_probe + (0.00000+0.00000i)*E_field_drive; -E_field.right = E_field_zero + (1.00000+0.00000i)*E_field_probe + (1.00000+0.00000i)*E_field_drive; -E_field.left = E_field_zero + (1.00000+0.00000i)*E_field_probe + (1.00000+0.00000i)*E_field_drive; +E_field_lab.linear = E_field_zero + (1.00000+0.00000i)*E_field_probe + (1.00000+0.00000i)*E_field_drive; +E_field_lab.right = E_field_zero + (0.00000+0.00000i)*E_field_probe + (0.00000+0.00000i)*E_field_drive; +E_field_lad.left = E_field_zero + (0.00000+0.00000i)*E_field_probe + (0.00000+0.00000i)*E_field_drive; + +% phi is angle between linear polarization and axis x +phi=pi/4; +% theta is angle between lab z axis (light porpagation dirextion) and magnetic field axis (z') +theta=pi/2; + +% we define light as linearly polarized +E_field_lab.x=cos(phi)*E_field_lab.linear; +E_field_lab.y=sin(phi)*E_field_lab.linear; +E_field_lab.z=E_field_zero; + +basis_transformation; % load subroutines +coord_transf_m = lin2circ * oldlin2newlin(theta); +E_field.right = coord_transf_m(1,1)*E_field_lab.x + coord_transf_m(1,2)*E_field_lab.y + coord_transf_m(1,3)*E_field_lab.z; +E_field.left = coord_transf_m(2,1)*E_field_lab.x + coord_transf_m(2,2)*E_field_lab.y + coord_transf_m(2,3)*E_field_lab.z; +E_field.linear = coord_transf_m(3,1)*E_field_lab.x + coord_transf_m(3,2)*E_field_lab.y + coord_transf_m(3,3)*E_field_lab.z; + fprintf (stderr, "tuning laser in forloop to set conditions vs detuning\n"); fflush (stderr); |