cloudy  trunk
 All Data Structures Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Variables
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
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
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(), FeIICollRatesBoltzmann(), Hydcs123(), hydro_vs_deexcit(), HyperfineCreate(), iso_collide(), lines_lv1_li_ne(), prt_constants(), Punch1LineData(), and tfidle().

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 mie_read_opc(), and prt_constants().

const double ELECTRON_MASS = 9.10938215e-28
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(), GrainCollHeating(), H2_Cooling(), H2_Create(), highen(), hmole_step(), and prt_constants().

const double EN1RYD = HPLANCK*SPEEDLIGHT*RYD_INF
const double ERG1CM = HPLANCK*SPEEDLIGHT
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 hmole_step(), 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
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 conorm(), 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(), PresTotCurrent(), 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(), ParseTable(), prt_constants(), and zero().

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
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(), and radius_increment().

const double PARSEC = AU*RADIAN*3600.

parsec in cm, 3.085678e18

Definition at line 127 of file physconst.h.

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

const double PI = 3.141592653589793238462643
const double PI2 = 6.283185307179586476925287

2*pi

Definition at line 32 of file physconst.h.

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

const double PI4 = 12.56637061435917295385057
const double PI8 = 25.13274122871834590770115

8*pi

Definition at line 38 of file physconst.h.

Referenced by CoolEvaluate(), Magnetic_evaluate(), MeanInc(), 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
const double RYDLAM = 1.e8/RYD_INF
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_LevelPops(), 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(), PresTotCurrent(), prt_constants(), PrtComment(), and PrtFinal().

const double SPEEDLIGHT = 2.99792458e10
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
const double TE1RYD = HPLANCK*SPEEDLIGHT*RYD_INF/BOLTZMANN
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(), FeIICreate(), GetGF(), he_1trans(), HydroEinstA(), lines_setup(), 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

Generated for cloudy by doxygen 1.8.4