physconst.h

Go to the documentation of this file.

/* This file is part of Cloudy and is copyright (C)1978-2013 by Gary J. Ferland and
 * others.  For conditions of distribution and use see copyright notice in license.txt */

#ifndef PHYSCONST_H_
#define PHYSCONST_H_

/*#include "physconst.h"*/

/*********************************************************************
 * first come math constants                                         *
 *********************************************************************/

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;

/*********************************************************************
 * astronomical constants go here                                    *
 *********************************************************************/

const double SOLAR_MASS = 1.9884e33;

const double SOLAR_LUMINOSITY = 3.8427e33;

/* >>refer      phys    const   http://pdg.lbl.gov/2010/reviews/rpp2010-rev-astrophysical-constants.pdf */
const double AU = 1.49597870700e13;

/*********************************************************************
 * fundamental constants go next, eventually rest should be defined  *
 * in terms of these, these are Codata 2010 values.                  *
 *********************************************************************/

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;

/*********************************************************************
 * below here should be derived constants                            *
 *                                                                   *
 * NB - explicit values in comments are approximate                  *
 *      and are not maintained !                                     *
 *********************************************************************/

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;

#endif /* PHYSCONST_H_ */