cloudy  trunk
Variables
physconst.h File Reference

Go to the source code of this file.

Variables

const UNUSED double EE = 2.718281828459045235360287
 
const UNUSED double EULER = 0.577215664901532860606512090082
 
const UNUSED double PI = 3.141592653589793238462643
 
const UNUSED double PI2 = 6.283185307179586476925287
 
const UNUSED double PI4 = 12.56637061435917295385057
 
const UNUSED double PI8 = 25.13274122871834590770115
 
const UNUSED double SQRT2 = 1.414213562373095048801689
 
const UNUSED double SQRTPI = 1.772453850905516027298167
 
const UNUSED double SQRTPIBY2 = 1.253314137315500251207883
 
const UNUSED double LN_TWO = 0.6931471805599453094172321
 
const UNUSED double LN_TEN = 2.302585092994045684017991
 
const UNUSED double LOG10_E = 0.4342944819032518276511289
 
const UNUSED double OPTDEP2EXTIN = 1.085736204758129569127822
 
const UNUSED double RADIAN = 57.29577951308232087679815
 
const UNUSED double SOLAR_MASS = 1.9884e33
 
const UNUSED double SOLAR_LUMINOSITY = 3.8427e33
 
const UNUSED double AU = 1.49597870700e13
 
const UNUSED double ATOMIC_MASS_UNIT = 1.660538921e-24
 
const UNUSED double ELECTRON_MASS = 9.10938291e-28
 
const UNUSED double PROTON_MASS = 1.672621777e-24
 
const UNUSED double BOLTZMANN = 1.3806488e-16
 
const UNUSED double SPEEDLIGHT = 2.99792458e10
 
const UNUSED double HPLANCK = 6.62606957e-27
 
const UNUSED double AVOGADRO = 6.0221415e23
 
const UNUSED double GRAV_CONST = 6.67384e-8
 
const UNUSED double ELEM_CHARGE = 1.602176565e-19
 
const UNUSED double RYD_INF = 1.0973731568539e5
 
const UNUSED double HIONPOT = 0.999466508345
 
const UNUSED double AS1RAD = RADIAN*3600.
 
const UNUSED double SQAS1SR = pow2(AS1RAD)
 
const UNUSED double SQAS_SKY = PI4*SQAS1SR
 
const UNUSED double PARSEC = AU*AS1RAD
 
const UNUSED double MEGAPARSEC = 1.e6*PARSEC
 
const UNUSED double H_BAR = HPLANCK/(2.*PI)
 
const UNUSED double ELEM_CHARGE_ESU = ELEM_CHARGE*SPEEDLIGHT/10.
 
const UNUSED double ELECTRIC_CONST = 1.e11/(PI4*pow2(SPEEDLIGHT))
 
const UNUSED double HION_LTE_POP = pow2(HPLANCK)/(PI2*BOLTZMANN*ELECTRON_MASS)
 
const UNUSED double SAHA = sqrt(pow3(HION_LTE_POP))
 
const UNUSED double ERG1CM = HPLANCK*SPEEDLIGHT
 
const UNUSED double T1CM = HPLANCK*SPEEDLIGHT/BOLTZMANN
 
const UNUSED double KJMOL1CM = ERG1CM*AVOGADRO/1e10
 
const UNUSED double WAVNRYD = 1./RYD_INF
 
const UNUSED double RYDLAM = 1.e8/RYD_INF
 
const UNUSED double EN1RYD = HPLANCK*SPEEDLIGHT*RYD_INF
 
const UNUSED double TE1RYD = HPLANCK*SPEEDLIGHT*RYD_INF/BOLTZMANN
 
const UNUSED double EVDEGK = ELEM_CHARGE*1.e7/BOLTZMANN
 
const UNUSED double EVRYD = HPLANCK*SPEEDLIGHT*RYD_INF/ELEM_CHARGE*1.e-7
 
const UNUSED double EN1EV = EN1RYD/EVRYD
 
const UNUSED double FR1RYD = SPEEDLIGHT*RYD_INF
 
const UNUSED double HNU3C2 = 2.*HPLANCK*SPEEDLIGHT*pow3(RYD_INF)
 
const UNUSED double FR1RYDHYD = SPEEDLIGHT*RYD_INF*HIONPOT
 
const UNUSED double HBAReV = H_BAR/EN1EV
 
const UNUSED double RYDLAMHYD = RYDLAM/HIONPOT
 
const UNUSED double STEFAN_BOLTZ = pow2(PI*pow2(BOLTZMANN))/(60.*pow3(H_BAR)*pow2(SPEEDLIGHT))
 
const UNUSED double FREQ_1EV = SPEEDLIGHT*RYD_INF/EVRYD
 
const UNUSED double FINE_STRUCTURE = pow2(ELEM_CHARGE_ESU)/SPEEDLIGHT/H_BAR
 
const UNUSED double FINE_STRUCTURE2 = pow2(FINE_STRUCTURE)
 
const UNUSED double BOHR_RADIUS_CM = FINE_STRUCTURE/(PI4*RYD_INF)
 
const UNUSED double TWO_PHOT_CONST = 9.*pow3(FINE_STRUCTURE2)*FR1RYD/2048.
 
const UNUSED double COLL_CONST = SAHA*BOLTZMANN/HPLANCK
 
const UNUSED double MILNE_CONST = SPEEDLIGHT*sqrt(pow3(FINE_STRUCTURE2)*pow3(TE1RYD)/PI)
 
const UNUSED 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

◆ AS1RAD

const UNUSED 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().

◆ ATOMIC_MASS_UNIT

const UNUSED double ATOMIC_MASS_UNIT = 1.660538921e-24

◆ AU

const UNUSED 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().

◆ AVOGADRO

const UNUSED 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().

◆ BOHR_RADIUS_CM

const UNUSED double BOHR_RADIUS_CM = FINE_STRUCTURE/(PI4*RYD_INF)

◆ BOLTZMANN

const UNUSED double BOLTZMANN = 1.3806488e-16

◆ COLL_CONST

const UNUSED 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().

◆ EE

const UNUSED double EE = 2.718281828459045235360287

the number e

Definition at line 23 of file physconst.h.

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

◆ ELECTRIC_CONST

const UNUSED 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().

◆ ELECTRON_MASS

const UNUSED double ELECTRON_MASS = 9.10938291e-28

◆ ELEM_CHARGE

const UNUSED 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().

◆ ELEM_CHARGE_ESU

const UNUSED 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().

◆ EN1EV

const UNUSED double EN1EV = EN1RYD/EVRYD

◆ EN1RYD

const UNUSED double EN1RYD = HPLANCK*SPEEDLIGHT*RYD_INF

◆ ERG1CM

const UNUSED double ERG1CM = HPLANCK*SPEEDLIGHT

◆ EULER

const UNUSED 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().

◆ EVDEGK

const UNUSED double EVDEGK = ELEM_CHARGE*1.e7/BOLTZMANN

◆ EVRYD

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

◆ FINE_STRUCTURE

const UNUSED 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().

◆ FINE_STRUCTURE2

const UNUSED double FINE_STRUCTURE2 = pow2(FINE_STRUCTURE)

the square of the fine-structure constant

Definition at line 219 of file physconst.h.

◆ FR1RYD

const UNUSED double FR1RYD = SPEEDLIGHT*RYD_INF

◆ FR1RYDHYD

const UNUSED 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().

◆ FREQ_1EV

const UNUSED 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().

◆ GRAV_CONST

const UNUSED 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().

◆ H_BAR

const UNUSED 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().

◆ HBAReV

const UNUSED 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().

◆ HION_LTE_POP

const UNUSED 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().

◆ HIONPOT

const UNUSED 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().

◆ HNU3C2

const UNUSED 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().

◆ HPLANCK

const UNUSED double HPLANCK = 6.62606957e-27

◆ KJMOL1CM

const UNUSED 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().

◆ LN_TEN

const UNUSED double LN_TEN = 2.302585092994045684017991

ln(10)

Definition at line 53 of file physconst.h.

Referenced by prt_constants().

◆ LN_TWO

const UNUSED double LN_TWO = 0.6931471805599453094172321

◆ LOG10_E

const UNUSED double LOG10_E = 0.4342944819032518276511289

log(e)

Definition at line 56 of file physconst.h.

Referenced by bhG_mx(), and prt_constants().

◆ MEGAPARSEC

const UNUSED double MEGAPARSEC = 1.e6*PARSEC

megaparsec in cm, 3.085678e24

Definition at line 141 of file physconst.h.

Referenced by GetHubbleFactor(), and GravitationalPressure().

◆ MILNE_CONST

const UNUSED 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().

◆ OPTDEP2EXTIN

const UNUSED 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().

◆ PARSEC

const UNUSED 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().

◆ PI

const UNUSED double PI = 3.141592653589793238462643

◆ PI2

const UNUSED double PI2 = 6.283185307179586476925287

2*pi

Definition at line 32 of file physconst.h.

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

◆ PI4

const UNUSED double PI4 = 12.56637061435917295385057

◆ PI8

const UNUSED double PI8 = 25.13274122871834590770115

◆ PROTON_MASS

const UNUSED 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().

◆ RADIAN

const UNUSED double RADIAN = 57.29577951308232087679815

180/pi

Definition at line 63 of file physconst.h.

Referenced by ParseIlluminate(), and prt_constants().

◆ RYD_INF

const UNUSED double RYD_INF = 1.0973731568539e5

◆ RYDLAM

const UNUSED double RYDLAM = 1.e8/RYD_INF

◆ RYDLAMHYD

const UNUSED 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().

◆ SAHA

const UNUSED 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().

◆ SOLAR_LUMINOSITY

const UNUSED 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().

◆ SOLAR_MASS

const UNUSED double SOLAR_MASS = 1.9884e33

◆ SPEEDLIGHT

const UNUSED double SPEEDLIGHT = 2.99792458e10

◆ SQAS1SR

const UNUSED 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().

◆ SQAS_SKY

const UNUSED 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().

◆ SQRT2

const UNUSED double SQRT2 = 1.414213562373095048801689

sqrt(2)

Definition at line 41 of file physconst.h.

Referenced by ConvTempEdenIoniz(), and prt_constants().

◆ SQRTPI

const UNUSED double SQRTPI = 1.772453850905516027298167

◆ SQRTPIBY2

const UNUSED double SQRTPIBY2 = 1.253314137315500251207883

sqrt(pi/2)

Definition at line 47 of file physconst.h.

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

◆ STEFAN_BOLTZ

const UNUSED 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().

◆ T1CM

const UNUSED double T1CM = HPLANCK*SPEEDLIGHT/BOLTZMANN

◆ TE1RYD

const UNUSED double TE1RYD = HPLANCK*SPEEDLIGHT*RYD_INF/BOLTZMANN

◆ TRANS_PROB_CONST

const UNUSED 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().

◆ TWO_PHOT_CONST

const UNUSED 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().

◆ WAVNRYD

const UNUSED double WAVNRYD = 1./RYD_INF