cloudy  trunk
 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Variables
physconst.h File Reference

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.9884e33
 
const double SOLAR_LUMINOSITY = 3.8427e33
 
const double AU = 1.49597870700e13
 
const double ATOMIC_MASS_UNIT = 1.660538921e-24
 
const double ELECTRON_MASS = 9.10938291e-28
 
const double PROTON_MASS = 1.672621777e-24
 
const double BOLTZMANN = 1.3806488e-16
 
const double SPEEDLIGHT = 2.99792458e10
 
const double HPLANCK = 6.62606957e-27
 
const double AVOGADRO = 6.0221415e23
 
const double GRAV_CONST = 6.67384e-8
 
const double ELEM_CHARGE = 1.602176565e-19
 
const double RYD_INF = 1.0973731568539e5
 
const double HIONPOT = 0.999466508345
 
const double AS1RAD = RADIAN*3600.
 
const double SQAS1SR = pow2(AS1RAD)
 
const double SQAS_SKY = PI4*SQAS1SR
 
const double PARSEC = AU*AS1RAD
 
const double MEGAPARSEC = 1.e6*PARSEC
 
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 KJMOL1CM = ERG1CM*AVOGADRO/1e10
 
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 AS1RAD = RADIAN*3600.

number of arcsec in 1 radian, 206264.806

Definition at line 129 of file physconst.h.

Referenced by conorm(), optimize_func(), prt_constants(), and PrtFinal().

const double ATOMIC_MASS_UNIT = 1.660538921e-24
const double AU = 1.49597870700e13

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

Definition at line 80 of file physconst.h.

Referenced by prt_constants().

const double AVOGADRO = 6.0221415e23

Avogadro constant – CoData 2002, http://physics.nist.gov/cgi-bin/cuu/Value?na|search_for=avogadro

Definition at line 106 of file physconst.h.

Referenced by t_mole_local::chem_heat(), and mole_partition_function().

const double BOHR_RADIUS_CM = FINE_STRUCTURE/(PI4*RYD_INF)
const double BOLTZMANN = 1.3806488e-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 229 of file physconst.h.

Referenced by ChiantiCollRate(), CoolCarb(), CS_l_mixing_PS64(), dBase_solve(), EdenChange(), FeIICollRatesBoltzmann(), Hydcs123(), hydro_vs_deexcit(), HyperfineCreate(), iso_collide(), lines_lv1_li_ne(), prt_constants(), Save1LineData(), StoutCollRate(), 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 150 of file physconst.h.

Referenced by mie_read_opc(), and prt_constants().

const double ELECTRON_MASS = 9.10938291e-28
const double ELEM_CHARGE = 1.602176565e-19

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

Definition at line 112 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 147 of file physconst.h.

Referenced by collision_strength_VF01(), ContSetIntensity(), EdenChange(), and prt_constants().

const double EN1EV = EN1RYD/EVRYD
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
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 216 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 219 of file physconst.h.

const double FR1RYD = SPEEDLIGHT*RYD_INF
const double FR1RYDHYD = SPEEDLIGHT*RYD_INF*HIONPOT

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

Definition at line 201 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 213 of file physconst.h.

Referenced by prt_constants().

const double GRAV_CONST = 6.67384e-8

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

Definition at line 109 of file physconst.h.

Referenced by GravitationalPressure(), highen(), InitSimPostparse(), PresTotCurrent(), and prt_constants().

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

h/2pi = 1.05457e-27

Definition at line 144 of file physconst.h.

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

const double HBAReV = H_BAR/EN1EV

H_BAR in eV sec, 6.582e-16

Definition at line 204 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 157 of file physconst.h.

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

const double HIONPOT = 0.999466508345

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

Definition at line 119 of file physconst.h.

Referenced by InitDefaultsPreparse(), ParseBackgrd(), ParseF_nu(), ParseRangeOption(), ParseTable(), 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 198 of file physconst.h.

Referenced by ContSetIntensity(), and prt_constants().

const double HPLANCK = 6.62606957e-27
const double KJMOL1CM = ERG1CM*AVOGADRO/1e10

kJ/mol per unit wavenumber

Definition at line 170 of file physconst.h.

Referenced by t_mole_global::make_species().

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
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 MEGAPARSEC = 1.e6*PARSEC

megaparsec in cm, 3.085678e24

Definition at line 141 of file physconst.h.

Referenced by GetHubbleFactor(), and GravitationalPressure().

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 233 of file physconst.h.

Referenced by iso_radrecomb_from_cross_section(), prt_constants(), RT_diffuse(), and RT_iso_integrate_RRC().

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*AS1RAD

parsec in cm, 3.085678e18

Definition at line 138 of file physconst.h.

Referenced by GravitationalPressure(), ParseDistance(), ParseGravity(), 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(), ParseMagnet(), prt_constants(), and ZoneStart().

const double PI4 = 12.56637061435917295385057
const double PI8 = 25.13274122871834590770115
const double PROTON_MASS = 1.672621777e-24

proton mass, gram

Definition at line 94 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 ParseIlluminate(), and prt_constants().

const double RYD_INF = 1.0973731568539e5
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 207 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 161 of file physconst.h.

Referenced by diatomics::H2_LevelPops(), mole_create_react(), mole_h_reactions(), and prt_constants().

const double SOLAR_LUMINOSITY = 3.8427e33

solar luminosity erg s-1 >>refer phys const http://pdg.lbl.gov/2010/reviews/rpp2010-rev-astrophysical-constants.pdf

Definition at line 75 of file physconst.h.

Referenced by prt_constants(), and PrtHeader().

const double SOLAR_MASS = 1.9884e33
const double SPEEDLIGHT = 2.99792458e10
const double SQAS1SR = pow2(AS1RAD)

number of square arcsec in 1 steradian, 4.254517e10

Definition at line 132 of file physconst.h.

Referenced by prt_constants().

const double SQAS_SKY = PI4*SQAS1SR

number of square arcsec in the whole sky, 5.3463838e11

Definition at line 135 of file physconst.h.

Referenced by conorm(), Flux::p_get(), Flux::p_set(), and prt_constants().

const double SQRT2 = 1.414213562373095048801689

sqrt(2)

Definition at line 41 of file physconst.h.

Referenced by ConvTempEdenIoniz(), and prt_constants().

const double SQRTPI = 1.772453850905516027298167
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 210 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 237 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 225 of file physconst.h.

Referenced by prt_constants().

const double WAVNRYD = 1./RYD_INF