Calculates the universal physics constants (see www.tcaep.co.uk/science/constant) INPUTS: None OUTPUTS: [qe,me,mp,mn,e0,mu0,re,mc2,clum,alpha,kB] - qe : Absolute value of the electron charge (C) - me : Electron rest mass (Kg) - mp : Proton rest mass (Kg) - mn : Neutron rest mass (Kg) - e0 : Free space permeability (F/m) - mu0 : Free space magnetic permeability (H/m) - re : Classical electron radius (m) - mc2 : Electron energy rest mass (keV) - clum : Speed of light (m/s) - alpha : Fine structure constant - kB : Boltzmann's constant (J/K-1) by Joan Decker <jodecker@mit.edu> (MIT/RLE) and Yves Peysson <yves.peysson@cea.fr> (CEA/DRFC)
0001 function [varargout] = pc_dke_yp(gauss) 0002 % 0003 % Calculates the universal physics constants 0004 % (see www.tcaep.co.uk/science/constant) 0005 % 0006 % INPUTS: None 0007 % 0008 % 0009 % OUTPUTS: [qe,me,mp,mn,e0,mu0,re,mc2,clum,alpha,kB] 0010 % 0011 % - qe : Absolute value of the electron charge (C) 0012 % - me : Electron rest mass (Kg) 0013 % - mp : Proton rest mass (Kg) 0014 % - mn : Neutron rest mass (Kg) 0015 % - e0 : Free space permeability (F/m) 0016 % - mu0 : Free space magnetic permeability (H/m) 0017 % - re : Classical electron radius (m) 0018 % - mc2 : Electron energy rest mass (keV) 0019 % - clum : Speed of light (m/s) 0020 % - alpha : Fine structure constant 0021 % - kB : Boltzmann's constant (J/K-1) 0022 % 0023 % by Joan Decker <jodecker@mit.edu> (MIT/RLE) and Yves Peysson <yves.peysson@cea.fr> (CEA/DRFC) 0024 % 0025 qe = 1.602176462e-19;%Absolute value of the electron charge (C) 0026 me = 9.10938188e-31;%Electron rest mass (Kg) 0027 mp = 1.67262158e-27;%Proton rest mass (Kg) 0028 mn = 1.67492716e-27;%Neutron rest mass (Kg) 0029 e0 = 8.854187818e-12;%Free space permittivity (F/m) 0030 mu0 = 12.566370614e-7;%Free space permeability (F/m) 0031 re = 2.817940285e-15;%Classical electron radius (m) 0032 mc2 = 510.998902;%Electron energy rest mass (keV) 0033 clum = 299792458;%Speed of light (m/s) 0034 alpha = 7.297352533e-3;%Fine structure constant (~1/137) 0035 kB = 1.3806504e-23;%Boltzmann's constant (J/K-1) 0036 % 0037 if nargin > 0,%CGS units 0038 qe = qe*clum*1e1;% statcoul 0039 me = me*1e3;% g 0040 mp = mp*1e3;% g 0041 mn = mn*1e3;% g 0042 e0 = NaN;% does not exist in CGS units 0043 mu0 = NaN;% does not exist in CGS units 0044 re = re*1e2;% cm 0045 clum = clum*1e2;% cm/s 0046 mc2 = mc2*qe*1e11/clum;% erg 0047 alpha = alpha; 0048 kB = kB*1e7; % erg/K 0049 end 0050 % 0051 if nargout >= 1, varargout{1} = qe;end 0052 if nargout >= 2, varargout{2} = me;end 0053 if nargout >= 3, varargout{3} = mp;end 0054 if nargout >= 4, varargout{4} = mn;end 0055 if nargout >= 5, varargout{5} = e0;end 0056 if nargout >= 6, varargout{6} = mu0;end 0057 if nargout >= 7, varargout{7} = re;end 0058 if nargout >= 8, varargout{8} = mc2;end 0059 if nargout >= 9, varargout{9} = clum;end 0060 if nargout >= 10, varargout{10} = alpha;end 0061 if nargout >= 11, varargout{11} = kB;end 0062 0063