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 );