% silica_plot_web_d
c = 299792458;
c_j = c^2/2 ;
ps = 1e-12 ;
meter_ps = (1/c) / ps ; % Picosecond meter for 1 sec
c_1000 = c / 1000 ; % light in kilometers
km_ps = (1 / c_1000) / ps ; % Light takes this many picoseconds to cover a kilometer.
% The refractive index multiplied by this number tells you how much longer a
% slower beam of light will take to cross the same distance.
% Manually Interpolated From silicaindexofrefractionversuswavelength.gif
si_data = [600 1.4576;650 1.4565;700 1.4554;750 1.4543;800 1.4533;850 1.4526;900 1.4517;…
950 1.4510;1000 1.4504;1050 1.4497;1100 1.4492;1150 1.4485;…
1200 1.4480;1250 1.4475;1300 1.4470;1350 1.4465;1400 1.4458;1450 1.4453;…
1500 1.4446;1550 1.4440;1600 1.4434;1650 1.4426;1700 1.4418] ;
% Manually Interpolated From Newport Fiber_Basics.pdf
silica_dispersion_data = [600 -290;650 -225;700 -170;750 -130;800 -105;850 -80;900 -64;…
950 -50;1000 -40;1050 -30;1100 -22;1150 -18;…
1200 -12;1250 -7;1300 +2;1350 8;1400 12;1450 15;…
1500 18;1550 20;1600 20;1650 20;1700 20] ;
% A crude calculation for the temporal dilation within silica.A complete
% equation would have to integrate all the temporal spaces in any form of
% dense matter at nearly atomic level scales. It is usually better than 50%
% accurate.
[t_rate_silica] = t_solid_pk(2203) ;
x_map_axis = [] ;
y1_map_axis = [] ;
y2_map_axis = [] ;
y3_map_axis = [] ;
y4_map_axis = [] ;
y5_map_axis = [] ;
y6_map_axis = [] ;
for step=1:1:23
si_dis = silica_dispersion_data(step,2) ; % Dispersion data
% Fresnel Parameters
c_time_rate_f = (si_data(step,2));
v_time_rate_f = 1/c_time_rate_f ;
v_photon_km_f =4.826672457e+06 -(c_time_rate_f * km_ps) ; % 1300nm zero
% Temporal Parameters Synthetic Refraction
photon_test=si_data(step,1);
e_nm = e_photon_nm(photon_test) ;
% Photon energy stress in spacetime
[photon_rate] = t_photon_ver_2(photon_test) ;
% Photon velocity reduction by temporal stress in a solid
photon_solid = t_rate_silica + (photon_rate / t_rate_silica) ;
% Temporal correction factor of velocity within time dilated spacetime
[v_true,v_dilation] = e_tk(photon_solid,1000) ;
c_v_dilation = (c – v_dilation) / c ;
v_photon_km_syn =3.3544e+06 – ((1 / c_v_dilation) * km_ps) ; % 1300nm zero
x_map_axis = [x_map_axis; photon_test] ;
y1_map_axis = [y1_map_axis; c_time_rate_f] ;
y2_map_axis = [y2_map_axis; v_time_rate_f] ;
y3_map_axis = [y3_map_axis; v_photon_km_f] ;
y4_map_axis = [y4_map_axis; si_dis] ;
y5_map_axis = [y5_map_axis; c_v_dilation] ;
y6_map_axis = [y6_map_axis; v_photon_km_syn] ;
end
ref_1300 = y6_map_axis(15,1) ;
end_plot = 6 ;
subplot(end_plot,1,1) ;
plot(x_map_axis,y1_map_axis,’r’) ; %Plot
legend(‘Silica Refractive Index’); %
legend(‘Location’,’northeast’) ;
grid on; %Turn on grid
xlabel(‘Photon NM’)
ylabel(‘Silica Refractive Index’)
subplot(end_plot,1,2) ;
plot(x_map_axis,y2_map_axis,’r’) ; %Plot
legend(‘Internal Dogma Velocity vs Wavelength’); %
legend(‘Location’,’southeast’) ;
grid on; %Turn on grid
xlabel(‘Photon NM’)
ylabel(‘Decimal Velocity c’)
subplot(end_plot,1,3) ;
plot(x_map_axis,y3_map_axis,’r’) ; %Plot
legend(‘Picoseconds Per Kilometer centered @1300nm’); %
legend(‘Location’,’southeast’) ;
grid on; %Turn on grid
xlabel(‘Photon NM’)
ylabel(‘Dogma Dispersion’)
subplot(end_plot,1,4) ;
plot(x_map_axis,y4_map_axis,’b’) ; %Plot
legend(‘Picoseconds Per Kilometer to 1500nm’ ); %
legend(‘Location’,’southeast’) ;
grid on; %Turn on grid
xlabel(‘Photon NM’)
ylabel(‘Silica Dispersion’)
subplot(end_plot,1,5) ;
plot(x_map_axis,y5_map_axis,’b’) ; %Plot
legend(‘Internal Synthetic Velocity vs Wavelength’); %
legend(‘Location’,’southeast’) ;
grid on; %Turn on grid
xlabel(‘Photon NM’)
ylabel(‘Decimal Velocity c’)
subplot(end_plot,1,6) ;
plot(x_map_axis,y6_map_axis,’b’) ; %Plot
legend(‘Picoseconds Per Kilometer to 1500nm’ ); %
legend(‘Location’,’southeast’) ;
grid on; %Turn on grid
xlabel(‘Photon NM’)
ylabel(‘Synthetic Dispersion’)