/home66/gary/public_html/cloudy/c08_branch/source/physconst.h File Reference

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Variables

const double EE = 2.718281828459045235360287
const double EULER = 0.577215664901532860606512090082
const double PI = 3.141592653589793238462643
const double PI2 = 6.283185307179586476925287
const double PI4 = 12.56637061435917295385057
const double PI8 = 25.13274122871834590770115
const double SQRT2 = 1.414213562373095048801689
const double SQRTPI = 1.772453850905516027298167
const double SQRTPIBY2 = 1.253314137315500251207883
const double LN_TWO = 0.6931471805599453094172321
const double LN_TEN = 2.302585092994045684017991
const double LOG10_E = 0.4342944819032518276511289
const double OPTDEP2EXTIN = 1.085736204758129569127822
const double RADIAN = 57.29577951308232087679815
const double SOLAR_MASS = 1.9889e33
const double SOLAR_LUMINOSITY = 3.846e33
const double AU = 1.49597870691e13
const double ATOMIC_MASS_UNIT = 1.660538782e-24
const double ELECTRON_MASS = 9.10938215e-28
const double PROTON_MASS = 1.672621637e-24
const double BOLTZMANN = 1.3806504e-16
const double SPEEDLIGHT = 2.99792458e10
const double HPLANCK = 6.62606896e-27
const double GRAV_CONST = 6.67428e-8
const double ELEM_CHARGE = 1.602176487e-19
const double RYD_INF = 1.0973731568527e5
const double HIONPOT = 0.999466508219
const double PARSEC = AU*RADIAN*3600.
const double H_BAR = HPLANCK/(2.*PI)
const double ELEM_CHARGE_ESU = ELEM_CHARGE*SPEEDLIGHT/10.
const double ELECTRIC_CONST = 1.e11/(PI4*pow2(SPEEDLIGHT))
const double HION_LTE_POP = pow2(HPLANCK)/(PI2*BOLTZMANN*ELECTRON_MASS)
const double SAHA = sqrt(pow3(HION_LTE_POP))
const double ERG1CM = HPLANCK*SPEEDLIGHT
const double T1CM = HPLANCK*SPEEDLIGHT/BOLTZMANN
const double WAVNRYD = 1./RYD_INF
const double RYDLAM = 1.e8/RYD_INF
const double EN1RYD = HPLANCK*SPEEDLIGHT*RYD_INF
const double TE1RYD = HPLANCK*SPEEDLIGHT*RYD_INF/BOLTZMANN
const double EVDEGK = ELEM_CHARGE*1.e7/BOLTZMANN
const double EVRYD = HPLANCK*SPEEDLIGHT*RYD_INF/ELEM_CHARGE*1.e-7
const double EN1EV = EN1RYD/EVRYD
const double FR1RYD = SPEEDLIGHT*RYD_INF
const double HNU3C2 = 2.*HPLANCK*SPEEDLIGHT*pow3(RYD_INF)
const double FR1RYDHYD = SPEEDLIGHT*RYD_INF*HIONPOT
const double HBAReV = H_BAR/EN1EV
const double RYDLAMHYD = RYDLAM/HIONPOT
const double STEFAN_BOLTZ = pow2(PI*pow2(BOLTZMANN))/(60.*pow3(H_BAR)*pow2(SPEEDLIGHT))
const double FREQ_1EV = SPEEDLIGHT*RYD_INF/EVRYD
const double FINE_STRUCTURE = pow2(ELEM_CHARGE_ESU)/SPEEDLIGHT/H_BAR
const double FINE_STRUCTURE2 = pow2(FINE_STRUCTURE)
const double BOHR_RADIUS_CM = FINE_STRUCTURE/(PI4*RYD_INF)
const double TWO_PHOT_CONST = 9.*pow3(FINE_STRUCTURE2)*FR1RYD/2048.
const double COLL_CONST = SAHA*BOLTZMANN/HPLANCK
const double MILNE_CONST = SPEEDLIGHT*sqrt(pow3(FINE_STRUCTURE2)*pow3(TE1RYD)/PI)
const double TRANS_PROB_CONST = PI4*HPLANCK*FINE_STRUCTURE/ELECTRON_MASS


Detailed Description

 * physical constants used by Cloudy, mostly taken from
 * >>refer	phys	const	Mohr P.J., Taylor B.N., & Newell D.B., Codata 2006, http://www.physics.nist.gov/constants
 * <BR><BR>
 * NB - these are all printed with the "print constants" command, 
 * which is in parse_print.cpp, so any new constants
 * added here must also be added to the prt_constants routine 
 * this is in the func_test case in the auto test suite 

Definition in file physconst.h.


Variable Documentation

const double ATOMIC_MASS_UNIT = 1.660538782e-24

atomic mass unit, gram

Definition at line 89 of file physconst.h.

Referenced by AbundancesPrt(), collision_strength_VF01(), CS_l_mixing_PS64(), dense_fabden(), iso_collide(), iso_satellite_update(), iso_state_lifetime(), mie_write_opc(), newspecies(), pah1_fun(), prt_constants(), rfield_opac_malloc(), ritoa(), and velset().

const double AU = 1.49597870691e13

astronomical unit, cm, nearly the length of the semimajor axis of the Earth's elliptical orbit around the sun

Definition at line 81 of file physconst.h.

Referenced by dense_fabden(), and prt_constants().

const double BOHR_RADIUS_CM = FINE_STRUCTURE/(PI4*RYD_INF)

Bohr radius in cm, 5.29177249e-9

Definition at line 205 of file physconst.h.

Referenced by collision_strength_VF01(), CS_PercivalRichards78(), prt_constants(), and ritoa().

const double BOLTZMANN = 1.3806504e-16

this is the Boltzmann factor, erg/K

Definition at line 98 of file physconst.h.

Referenced by AgeCheck(), collision_strength_VF01(), Fe4Lev12(), GrainDrift(), GrainElecRecomb1(), GrainScreen(), grn_abs(), HydroRecCool(), lines_molecules(), ParseConstant(), prt_constants(), PrtZone(), qheat_init(), rfield_opac_malloc(), uderiv(), and velset().

const double COLL_CONST = SAHA*BOLTZMANN/HPLANCK

this is the square of the value roughly equal to 8.629e-6 that appears in converting collision strengths to rates. The constant is h^2/((2PI*me)^3/2 * k^1/2).

Definition at line 212 of file physconst.h.

Referenced by CoolCarb(), CS_l_mixing_PS64(), Hydcs123(), hydro_vs_deexcit(), HyperfineCreate(), iso_collide(), lines_lv1_li_ne(), prt_constants(), and Punch1LineData().

const double EE = 2.718281828459045235360287

the number e

Definition at line 23 of file physconst.h.

Referenced by ConvInitSolution(), Hydcs123(), and prt_constants().

const double ELECTRIC_CONST = 1.e11/(PI4*pow2(SPEEDLIGHT))

electric constant, in F/m, 8.854e-12

Definition at line 136 of file physconst.h.

Referenced by prt_constants().

const double ELECTRON_MASS = 9.10938215e-28

electron mass, gram

Definition at line 92 of file physconst.h.

Referenced by CO_Init(), CS_l_mixing_PS64(), GrainElecRecomb1(), highen(), hv(), Hydcs123(), iso_collide(), iso_satellite_update(), iso_state_lifetime(), prt_constants(), and ritoa().

const double ELEM_CHARGE = 1.602176487e-19

elementary charge, in C in SI units, to use this must convert to cgs

Definition at line 110 of file physconst.h.

Referenced by GrainScreen(), one_elec(), and prt_constants().

const double ELEM_CHARGE_ESU = ELEM_CHARGE*SPEEDLIGHT/10.

elementary charge, in ESU, 4.8032e-10

Definition at line 133 of file physconst.h.

Referenced by collision_strength_VF01(), and prt_constants().

const double EN1EV = EN1RYD/EVRYD

ergs per eV, 1.602176e-012

Definition at line 175 of file physconst.h.

Referenced by ChargTranSumHeat(), CS_l_mixing_S62(), and prt_constants().

const double EN1RYD = HPLANCK*SPEEDLIGHT*RYD_INF

ergs per inf mass Ryd, 2.180e-11

Definition at line 162 of file physconst.h.

Referenced by AGN_Hemis(), atom_level3(), collision_strength_VF01(), ContSetIntensity(), FillExtraLymanLine(), GammaBn(), GammaBnPL(), GammaK(), GetProbDistr_HighLimit(), GetProbDistr_LowLimit(), GrainDrift(), highen(), hmole_reactions(), hv(), HydroCSInterp(), ion_collis(), iso_create(), lindst(), lines_continuum(), OscStr_f(), OscStr_f_log10(), outsum(), prt_constants(), PrtHeader(), PrtZone(), PunchSpecial(), qheat(), qheat_init(), RebinQHeatResults(), RT_line_static(), Spec_cont(), sumcon(), TryDoubleStep(), and uderiv().

const double ERG1CM = HPLANCK*SPEEDLIGHT

number of ergs per wavenumber, 1.9864e-16

Definition at line 150 of file physconst.h.

Referenced by atmdat_Chianti_readin(), atmdat_lamda_readin(), DoSatelliteLines(), FeIIPoint(), iso_satellite(), iso_satellite_update(), and prt_constants().

const double EULER = 0.577215664901532860606512090082

the Euler constant (aka Euler-Mascheroni constant or gamma)

Definition at line 26 of file physconst.h.

Referenced by expn(), and prt_constants().

const double EVDEGK = ELEM_CHARGE*1.e7/BOLTZMANN

Kelvins per eV, 1.1604e4

Definition at line 169 of file physconst.h.

Referenced by Hydcs123(), hydro_vs_ioniz(), ligbar(), ParseConstant(), prt_constants(), S62_Therm_ave_coll_str(), and tfidle().

const double EVRYD = HPLANCK*SPEEDLIGHT*RYD_INF/ELEM_CHARGE*1.e-7

eV per inf mass Ryd, 13.606

Definition at line 172 of file physconst.h.

Referenced by AnuUnit(), atmdat_dielrec_fe(), t_ADfA::coll_ion(), cross_section(), CS_l_mixing_S62(), CS_VS80(), elec_esc_length(), GrainCharge(), GrainDrive(), Hion_coll_ioniz_ratecoef(), Hydcs123(), hydro_vs_coll_str(), hydro_vs_deexcit(), hydro_vs_ioniz(), t_yield::init_yield(), InitBinAugerData(), ipShells(), iso_create(), one_elec(), Opacity_iso_photo_cs(), OpacityCreate1Element(), OpacityCreateReilMan(), pah1_fun(), prt_constants(), punchFITS_EnergyData(), ReadAugerData(), ReadTable(), S62_Therm_ave_coll_str(), Therm_ave_coll_str_int_VS80(), UpdatePot1(), y0b(), Yan_H2_CS(), and Yfunc().

const double FINE_STRUCTURE = pow2(ELEM_CHARGE_ESU)/SPEEDLIGHT/H_BAR

the fine-structure constant a= 2pi e^2/hc 7.297 352 533 x 10-3

Definition at line 199 of file physconst.h.

Referenced by iso_state_lifetime(), and prt_constants().

const double FINE_STRUCTURE2 = pow2(FINE_STRUCTURE)

the square of the fine-structure constant

Definition at line 202 of file physconst.h.

const double FR1RYD = SPEEDLIGHT*RYD_INF

frequency of one Ryd for infinite mass nuclei, 3.289842e15

Definition at line 178 of file physconst.h.

Referenced by ContSetIntensity(), GrainMakeDiffuse(), lgCompileAtmosphere(), lgValidModel(), ParseF_nu(), PlanckIntegral(), plankf(), prt_constants(), and rayleh().

const double FR1RYDHYD = SPEEDLIGHT*RYD_INF*HIONPOT

frequency of ionization potential of H (not inf mass), 3.288087e15 - never used

Definition at line 184 of file physconst.h.

Referenced by prt_constants().

const double FREQ_1EV = SPEEDLIGHT*RYD_INF/EVRYD

the frequency of one eV, 2.418e14

Definition at line 196 of file physconst.h.

Referenced by prt_constants().

const double GRAV_CONST = 6.67428e-8

Gravitational constant, cm^3/g/s^2

Definition at line 107 of file physconst.h.

Referenced by InitSimPostparse(), and prt_constants().

const double H_BAR = HPLANCK/(2.*PI)

h/2pi = 1.05457e-27

Definition at line 130 of file physconst.h.

Referenced by collision_strength_VF01(), iso_state_lifetime(), prt_constants(), and ritoa().

const double HBAReV = H_BAR/EN1EV

H_BAR in eV sec, 6.582e-16

Definition at line 187 of file physconst.h.

Referenced by prt_constants().

const double HION_LTE_POP = pow2(HPLANCK)/(PI2*BOLTZMANN*ELECTRON_MASS)

this is the factor that appears in front of Boltzmann factor to get LTE level populations for hydrogenic ions. It is given in the first parts of section 5 of part 2 of hazy, and is actually ( planck^2 / (2 pi m_e k ) )^3/2, but cannot evaluate powers here, must raise this to 3/2 when used, HION_LTE_POP = 5.556e-11 cm^2 K

Definition at line 143 of file physconst.h.

Referenced by iso_collide(), iso_satellite_update(), and prt_constants().

const double HIONPOT = 0.999466508219

ionization potential of real hydrogen atom, in inf mass ryd, based on Codata 98, uncertainty 24e-12, calculated by Peter van Hoof

Definition at line 117 of file physconst.h.

Referenced by ParseBackgrd(), ParseF_nu(), ParseRangeOption(), and prt_constants().

const double HNU3C2 = 2.*HPLANCK*SPEEDLIGHT*pow3(RYD_INF)

2 h FR1RYD^3 / c^2 for infinite mass nucleus, 0.5250

Definition at line 181 of file physconst.h.

Referenced by ContSetIntensity(), and prt_constants().

const double HPLANCK = 6.62606896e-27

Planck's constant

Definition at line 104 of file physconst.h.

Referenced by H_Einstein_A_lin(), H_Einstein_A_log10(), PlanckIntegral(), PressureRadiationLine(), and prt_constants().

const double LN_TEN = 2.302585092994045684017991

ln(10)

Definition at line 53 of file physconst.h.

Referenced by prt_constants().

const double LN_TWO = 0.6931471805599453094172321

ln(2)

Definition at line 50 of file physconst.h.

Referenced by ASINH(), GetProbDistr_HighLimit(), GetProbDistr_LowLimit(), prt_constants(), and qheat().

const double LOG10_E = 0.4342944819032518276511289

log(e)

Definition at line 56 of file physconst.h.

Referenced by bhG_mx(), and prt_constants().

const double MILNE_CONST = SPEEDLIGHT*sqrt(pow3(FINE_STRUCTURE2)*pow3(TE1RYD)/PI)

this is the square of the value roughly equal to 4.123e11 that appears in the integration of photoionization cross-sections to obtain recombination coefficients.

Definition at line 216 of file physconst.h.

Referenced by helike_radrecomb_from_cross_section(), hlike_radrecomb_from_cross_section(), prt_constants(), and RT_diffuse().

const double OPTDEP2EXTIN = 1.085736204758129569127822

factor that converts optical depth into extinction in mags, 2.5 log e

Definition at line 60 of file physconst.h.

Referenced by prt_constants().

const double PARSEC = AU*RADIAN*3600.

parsec in cm, 3.085678e18

Definition at line 127 of file physconst.h.

Referenced by dense_fabden(), ParseRadius(), ParseStop(), and prt_constants().

const double PI = 3.141592653589793238462643

pi

Definition at line 29 of file physconst.h.

Referenced by AngerJ(), bhG(), bhG_mx(), collision_strength_VF01(), CS_PercivalRichards78(), DoFSMixing(), FreeFreeGaunt(), GetProbDistr_HighLimit(), GrainElecRecomb1(), GrainScreen(), grn_abs(), HydroOscilStr(), Jint(), L_mix_integrand_VF01(), lgCompileAtmosphereCoStar(), mie_auxiliary(), mie_cs(), pah1_fun(), prt_constants(), RauchInitializeSub(), ritoa(), RT_LineWidth(), S62_Therm_ave_coll_str(), scqdri(), StarkCollTransProb_VF01(), Stognienko(), velset(), and vib_evap().

const double PI2 = 6.283185307179586476925287

2*pi

Definition at line 32 of file physconst.h.

Referenced by collision_strength_VF01(), ParseMagnet(), and prt_constants().

const double PI4 = 12.56637061435917295385057

4*pi

Definition at line 35 of file physconst.h.

Referenced by ContSetIntensity(), FeIIPoint(), FillExtraLymanLine(), GrainMakeDiffuse(), iso_collapsed_lifetimes_update(), iso_satellite_update(), Magnetic_evaluate(), mie_calc_ial(), ParseBackgrd(), ParseTable(), PlanckIntegral(), and prt_constants().

const double PI8 = 25.13274122871834590770115

8*pi

Definition at line 38 of file physconst.h.

Referenced by Magnetic_evaluate(), PressureRadiationLine(), and prt_constants().

const double PROTON_MASS = 1.672621637e-24

proton mass, gram

Definition at line 95 of file physconst.h.

Referenced by collision_strength_VF01(), hv(), Hydcs123(), and prt_constants().

const double RADIAN = 57.29577951308232087679815

180/pi

Definition at line 63 of file physconst.h.

Referenced by prt_constants().

const double RYD_INF = 1.0973731568527e5

infinite mass rydberg constant, in cm^-1

Definition at line 113 of file physconst.h.

Referenced by AnuUnit(), FillExtraLymanLine(), he_1trans(), helike_energy(), helike_quantum_defect(), hydro_vs_coll_str(), hydro_vs_deexcit(), HydroEinstA(), iso_create(), pah1_fun(), prt_constants(), and RT_diffuse().

const double RYDLAM = 1.e8/RYD_INF

Angstrom per infinite mass Ryd, 911.2671

Definition at line 159 of file physconst.h.

Referenced by AnuUnit(), AtlasCompile(), ContBandsCreate(), CoStarCompile(), iso_satellite(), lgCompileAtmosphere(), lgCompileAtmosphereCoStar(), lindst(), PntForLine(), prt_constants(), PunchNewContinuum(), PunchSpecial(), read_hm05(), and ReadTable().

const double RYDLAMHYD = RYDLAM/HIONPOT

wavelength (A) of ionization potential of Hydrogen, 911.7535 - never used

Definition at line 190 of file physconst.h.

Referenced by prt_constants().

const double SAHA = sqrt(pow3(HION_LTE_POP))

SAHA is ( h^2/2/pi/m/k )^3/2, is correct constant for free electron SAHA = 4.14132e-16 cm^3 K^(3/2)

Definition at line 147 of file physconst.h.

Referenced by H2_X_coll_rate_evaluate(), hmole_reactions(), and prt_constants().

const double SOLAR_LUMINOSITY = 3.846e33

solar luminosity erg s-1 >>refer phys const http://pdg.lbl.gov/2000/astrorpp.pdf

Definition at line 75 of file physconst.h.

Referenced by prt_constants(), and PrtHeader().

const double SOLAR_MASS = 1.9889e33

solar mass in gram >>refer phys const http://pdg.lbl.gov/2000/astrorpp.pdf

Definition at line 71 of file physconst.h.

Referenced by InitSimPostparse(), and prt_constants().

const double SPEEDLIGHT = 2.99792458e10

speed of light, cm/s

Definition at line 101 of file physconst.h.

Referenced by AnuUnit(), GrainDrift(), GrainMakeDiffuse(), highen(), InitDefaultsPreparse(), iso_state_lifetime(), lgCompileAtmosphere(), ParseSet(), PlanckIntegral(), prt_constants(), PrtZone(), rfield_opac_malloc(), and ritoa().

const double SQRT2 = 1.414213562373095048801689

sqrt(2)

Definition at line 41 of file physconst.h.

Referenced by prt_constants().

const double SQRTPI = 1.772453850905516027298167

sqrt(pi)

Definition at line 44 of file physconst.h.

Referenced by esc_CRDwing_1side(), esca0k2(), prt_constants(), RT_DestProb(), RTesc_lya(), and RTesc_lya_1side().

const double SQRTPIBY2 = 1.253314137315500251207883

sqrt(pi/2)

Definition at line 47 of file physconst.h.

Referenced by bhG(), bhG_mx(), and prt_constants().

const double STEFAN_BOLTZ = pow2(PI*pow2(BOLTZMANN))/(60.*pow3(H_BAR)*pow2(SPEEDLIGHT))

Stefan-Boltzmann constant, 5.6704e-5

Definition at line 193 of file physconst.h.

Referenced by lgValidModel(), ParseBackgrd(), ParseBlackbody(), ParseCMB(), and prt_constants().

const double T1CM = HPLANCK*SPEEDLIGHT/BOLTZMANN

degrees kelvin per unit wavenumber, 1.4388

Definition at line 153 of file physconst.h.

Referenced by atmdat_Chianti_readin(), atmdat_lamda_readin(), atmol_popsolve(), atom_levelN(), atom_pop5(), CoolCarb(), Fe4Lev12(), FeIIPoint(), H21_cm_pops(), H2_Level_low_matrix(), hmole_reactions(), iso_satellite(), prt_constants(), and tfidle().

const double TE1RYD = HPLANCK*SPEEDLIGHT*RYD_INF/BOLTZMANN

the temperature of 1 Rydberg te1ryd is h/k is temp of 1 Rydberg, 1.579e5

Definition at line 166 of file physconst.h.

Referenced by atmdat_dielrec_fe(), t_ADfA::coll_ion(), cont_gaunt_calc(), CS_l_mixing_VF01(), CS_ThermAve_PR78(), CS_VS80(), dgaunt(), ffun1(), FillExtraLymanLine(), FreeFreeGaunt(), GrainDrift(), GrainMakeDiffuse(), helike_radrecomb_from_cross_section(), Hion_coll_ioniz_ratecoef(), hlike_radrecomb_from_cross_section(), Hydcs123(), hydro_vs_deexcit(), iso_cool(), iso_create(), ParseAgn(), ParseBlackbody(), ParsePowerlawContinuum(), PlanckIntegral(), prt_constants(), PunchGaunts(), Recomb_Seaton59(), RT_recom_effic(), tfidle(), Therm_ave_coll_str_int_VS80(), and UpdatePot2().

const double TRANS_PROB_CONST = PI4*HPLANCK*FINE_STRUCTURE/ELECTRON_MASS

This is the constant used in converting oscillator strengths to As. The formula is Aul = TRANS_PROB_CONST * f(u,l) * wavenumber^2. TRANS_PROB_CONST is 0.667025

Definition at line 220 of file physconst.h.

Referenced by CS_l_mixing_S62(), DoFSMixing(), eina(), GetGF(), he_1trans(), HydroEinstA(), prt_constants(), and ritoa().

const double TWO_PHOT_CONST = 9.*pow3(FINE_STRUCTURE2)*FR1RYD/2048.

the two photon constant as defined by Breit & Teller, as in equation 4 of Spitzer & Greenstein 51, 2.18313

Definition at line 208 of file physconst.h.

Referenced by prt_constants().

const double WAVNRYD = 1./RYD_INF

number of Ryd per wavenumber, 9.11267e-6

Definition at line 156 of file physconst.h.

Referenced by ColStrGBar(), CS_l_mixing_S62(), FeIIPoint(), FillExtraLymanLine(), H2_ContPoint(), iso_create(), MakeCS(), mie_calc_ial(), ParsePunch(), prt_constants(), pun1Line(), and tbl_fun().


Generated on Mon Feb 16 12:11:19 2009 for cloudy by  doxygen 1.4.7