diff options
Diffstat (limited to 'faraday/psr_vs_detuning_combo.m')
| -rw-r--r-- | faraday/psr_vs_detuning_combo.m | 149 |
1 files changed, 149 insertions, 0 deletions
diff --git a/faraday/psr_vs_detuning_combo.m b/faraday/psr_vs_detuning_combo.m new file mode 100644 index 0000000..b2ea9d3 --- /dev/null +++ b/faraday/psr_vs_detuning_combo.m @@ -0,0 +1,149 @@ +1; + +data_dir='results/'; +output_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) +ret=ouput_psr_vs_detuning_combo( ... + 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 ... + , detuning_freq, B_field, theta, phi, psi_el + , output_dir ); + +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%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) +ret=ouput_psr_vs_detuning_combo( ... + 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 ... + , detuning_freq, B_field, theta, phi, psi_el + , output_dir ); + +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%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) +ret=ouput_psr_vs_detuning_combo( ... + 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 ... + , detuning_freq, B_field, theta, phi, psi_el ... + , output_dir ); + +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%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) +ret=ouput_psr_vs_detuning_combo( ... + 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 ... + , detuning_freq, B_field, theta, phi, psi_el ... + , output_dir ); + +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%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) +ret=ouput_psr_vs_detuning_combo( ... + 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 ... + , detuning_freq, B_field, theta, phi, psi_el ... + , output_dir ); + +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%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) +ret=ouput_psr_vs_detuning_combo( ... + 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 ... + , detuning_freq, B_field, theta, phi, psi_el ... + , output_dir ); + +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%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) +ret=ouput_psr_vs_detuning_combo( ... + 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 ... + , detuning_freq, B_field, theta, phi, psi_el ... + , output_dir ); |