1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
|
1;
data_dir='results/';
N_detun_steps=1000;
detuning_p_min=-200.0;
%detuning_p_max=-detuning_p_min;
detuning_p_max=1000;
detuning_freq=linspace(detuning_p_min,detuning_p_max,N_detun_steps);
gmg=.7; % gyro magnetic ration for ground level
zeeman_splitting=+0.000;
B_field=zeeman_splitting/gmg;
%[psr_rad]=psr_vs_detuning(Ep, psi_el, B_field, theta, phi)
% phi is angle between linear polarization and axis x
phi=pi/4;
% theta is angle between lab z axis (light propagation direction) and magnetic field axis (z')
theta=0;
% psi_el is the ellipticity parameter (phase difference between left and right polarization)
psi_el=-30/180*pi;
figure(6);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
% zero magnetic field,, 30 degree ellipticity
zeeman_splitting=+0.000;
B_field=zeeman_splitting/gmg;
psi_el=30/180*pi;
%[psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, psr_rad_grEp_pos_el, psr_rad_grEp_neg_el] =make_representative_psr_vs_detuning_for_given_B_and_psi_el(detuning_freq, B_field, psi_el, theta, phi);
fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el)
load(fname)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
% 0.1 G magnetic field,, 30 degree ellipticity
zeeman_splitting=+0.070;
B_field=zeeman_splitting/gmg;
psi_el=30/180*pi;
%[psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, psr_rad_grEp_pos_el, psr_rad_grEp_neg_el] =make_representative_psr_vs_detuning_for_given_B_and_psi_el(detuning_freq, B_field, theta, phi, psi_el);
fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el)
load(fname)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
% 0.0001 G magnetic field,, 30 degree ellipticity
zeeman_splitting=+0.000070;
B_field=zeeman_splitting/gmg;
psi_el=30/180*pi;
%[psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, psr_rad_grEp_pos_el, psr_rad_grEp_neg_el] =make_representative_psr_vs_detuning_for_given_B_and_psi_el(detuning_freq, B_field, theta, phi, psi_el);
fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el)
load(fname)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
% 1.0 G magnetic field,, 30 degree ellipticity
zeeman_splitting=+0.70;
B_field=zeeman_splitting/gmg;
psi_el=30/180*pi;
%[psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, psr_rad_grEp_pos_el, psr_rad_grEp_neg_el] =make_representative_psr_vs_detuning_for_given_B_and_psi_el(detuning_freq, B_field, theta, phi, psi_el);
fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el)
load(fname)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
% -0.1 G magnetic field,, 30 degree ellipticity
zeeman_splitting=-0.070;
B_field=zeeman_splitting/gmg;
psi_el=30/180*pi;
%[psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, psr_rad_grEp_pos_el, psr_rad_grEp_neg_el] =make_representative_psr_vs_detuning_for_given_B_and_psi_el(detuning_freq, B_field, theta, phi, psi_el);
fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el)
load(fname)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
% -0.0001 G magnetic field,, 30 degree ellipticity
zeeman_splitting=-0.000070;
B_field=zeeman_splitting/gmg;
psi_el=30/180*pi;
%[psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, psr_rad_grEp_pos_el, psr_rad_grEp_neg_el] =make_representative_psr_vs_detuning_for_given_B_and_psi_el(detuning_freq, B_field, theta, phi, psi_el);
fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el)
load(fname)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
% -1.0 G magnetic field,, 30 degree ellipticity
zeeman_splitting=-0.70;
B_field=zeeman_splitting/gmg;
psi_el=30/180*pi;
%[psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, psr_rad_grEp_pos_el, psr_rad_grEp_neg_el] =make_representative_psr_vs_detuning_for_given_B_and_psi_el(detuning_freq, B_field, theta, phi, psi_el);
fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el)
load(fname)
|
