/home66/gary/public_html/cloudy/c13_branch/source/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 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 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 226 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().

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

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

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

const double FR1RYD = SPEEDLIGHT*RYD_INF

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

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

Referenced by prt_constants().

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

Referenced by ContSetIntensity(), and prt_constants().

const double HPLANCK = 6.62606957e-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
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 230 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 204 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 207 of file physconst.h.

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

const double T1CM = HPLANCK*SPEEDLIGHT/BOLTZMANN

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

Referenced by prt_constants().

const double WAVNRYD = 1./RYD_INF
 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines

Generated on 15 Nov 2012 for cloudy by  doxygen 1.6.1