cloudy trunk
Loading...
Searching...
No Matches
physconst.h File Reference

Go to the source code of this file.

Variables

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

◆ AS1RAD

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

◆ ATOMIC_MASS_UNIT

◆ AU

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

◆ AVOGADRO

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

◆ BOHR_RADIUS_CM

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

◆ BOLTZMANN

◆ COLL_CONST

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

◆ EE

UNUSED const double EE = 2.718281828459045235360287

the number e

Definition at line 23 of file physconst.h.

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

◆ ELECTRIC_CONST

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

◆ ELECTRON_MASS

◆ ELEM_CHARGE

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

◆ ELEM_CHARGE_ESU

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

◆ EN1EV

UNUSED const double EN1EV = EN1RYD/EVRYD

◆ EN1RYD

◆ ERG1CM

◆ EULER

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

◆ EVDEGK

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

◆ EVRYD

◆ FINE_STRUCTURE

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

◆ FINE_STRUCTURE2

UNUSED const double FINE_STRUCTURE2 = pow2(FINE_STRUCTURE)

the square of the fine-structure constant

Definition at line 219 of file physconst.h.

◆ FR1RYD

UNUSED const double FR1RYD = SPEEDLIGHT*RYD_INF

◆ FR1RYDHYD

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

◆ FREQ_1EV

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

◆ GRAV_CONST

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

◆ H_BAR

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

◆ HBAReV

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

◆ HION_LTE_POP

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

◆ HIONPOT

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

◆ HNU3C2

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

◆ HPLANCK

UNUSED const double HPLANCK = 6.62606957e-27

◆ KJMOL1CM

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

◆ LN_TEN

UNUSED const double LN_TEN = 2.302585092994045684017991

ln(10)

Definition at line 53 of file physconst.h.

Referenced by prt_constants().

◆ LN_TWO

UNUSED const double LN_TWO = 0.6931471805599453094172321

◆ LOG10_E

UNUSED const double LOG10_E = 0.4342944819032518276511289

log(e)

Definition at line 56 of file physconst.h.

Referenced by bhG_mx(), and prt_constants().

◆ MEGAPARSEC

UNUSED const 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

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

◆ OPTDEP2EXTIN

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

◆ PARSEC

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

◆ PI

◆ PI2

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

◆ PI4

◆ PI8

UNUSED const double PI8 = 25.13274122871834590770115

◆ PROTON_MASS

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

◆ RADIAN

UNUSED const double RADIAN = 57.29577951308232087679815

180/pi

Definition at line 63 of file physconst.h.

Referenced by ParseIlluminate(), and prt_constants().

◆ RYD_INF

◆ RYDLAM

◆ RYDLAMHYD

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

◆ SAHA

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

◆ SOLAR_LUMINOSITY

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

◆ SOLAR_MASS

◆ SPEEDLIGHT

◆ SQAS1SR

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

◆ SQAS_SKY

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

◆ SQRT2

UNUSED const double SQRT2 = 1.414213562373095048801689

sqrt(2)

Definition at line 41 of file physconst.h.

Referenced by ConvTempEdenIoniz(), and prt_constants().

◆ SQRTPI

UNUSED const double SQRTPI = 1.772453850905516027298167

◆ SQRTPIBY2

UNUSED const 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

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

◆ T1CM

◆ TE1RYD

◆ TRANS_PROB_CONST

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

◆ TWO_PHOT_CONST

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

◆ WAVNRYD