diff options
| -rw-r--r-- | compass_lin.m | 20 |
1 files changed, 15 insertions, 5 deletions
diff --git a/compass_lin.m b/compass_lin.m index 0ebd20c..18316f8 100644 --- a/compass_lin.m +++ b/compass_lin.m @@ -70,10 +70,14 @@ phi=pi*2/8; theta=0/4; theta=pi/2; %theta=65/180*pi; + %small ellipticity angle psi (0 to pi/2) % 0 will give linear polarization -% pi/2 circular polarization -psi=pi/2; +% pi/4 circular polarization +%psi_el=pi/4; +psi_el=.2*pi/4; +%psi_el=0*pi/4; + @@ -97,11 +101,17 @@ for phi=phis; wm=wd-(wp-wd); light_positive_freq=[ wp, wd, wp-wd]; - % we define light as linearly polarized - % where phi is angle between light polarization and axis x - [E_field_lab_pos_freq.x, E_field_lab_pos_freq.y] = rotXpolarization(phi, E_field_lab_pos_freq.linear); + % we define light as linearly polarized along x + E_field_lab_pos_freq.x = E_field_lab_pos_freq.linear; + E_field_lab_pos_freq.y = 0; + % now we add small elasticity + E_field_lab_pos_freq.y=sin(psi_el)*E_field_lab_pos_freq.x*(1i); % order is important + E_field_lab_pos_freq.x=cos(psi_el)*E_field_lab_pos_freq.x; + % set phi angle between light polarization and axis x + [E_field_lab_pos_freq.x, E_field_lab_pos_freq.y] = rotLinPolarization(phi, E_field_lab_pos_freq.x, E_field_lab_pos_freq.y); E_field_lab_pos_freq.z=E_field_zero; + % now we transfor x,y,z, to x',y', and z' with respect to magnetic field az z' axis E_field_pos_freq=xyz_lin2atomic_axis_polarization(theta, E_field_lab_pos_freq); |