Constants.H

Go to the documentation of this file.

 
#ifndef _CONSTANTS_H_
#define _CONSTANTS_H_
 
#include <gsl/gsl_const_cgsm.h>
#include <gsl/gsl_const_mksa.h>
#include <gsl/gsl_const_num.h>
#include "../Definitions.H"
#include "Types.H"
#include "Units.H"
 
namespace criptic {
 
  namespace constants {
    
#ifdef CRIPTIC_UNITS_CGS
  
    // CGS units
    static constexpr Real c
    = GSL_CONST_CGSM_SPEED_OF_LIGHT;  
    static constexpr Real mp
    = GSL_CONST_CGSM_MASS_PROTON;     
    static constexpr Real me
    = GSL_CONST_CGSM_MASS_ELECTRON;   
    static constexpr Real sigmaT
    = GSL_CONST_CGSM_THOMSON_CROSS_SECTION;  
    static constexpr Real kB =
      GSL_CONST_CGSM_BOLTZMANN;       
    static constexpr Real h =
      GSL_CONST_CGSM_PLANCKS_CONSTANT_H;  
    static constexpr Real hBar =
      GSL_CONST_CGSM_PLANCKS_CONSTANT_HBAR;  
    static constexpr Real e =
      GSL_CONST_CGSM_ELECTRON_CHARGE * c; 
    static constexpr Real a0 =
      GSL_CONST_CGSM_BOHR_RADIUS; 
    static constexpr Real Ryd =
      GSL_CONST_CGSM_RYDBERG; 
    static constexpr Real sigmaSB =
      GSL_CONST_CGSM_STEFAN_BOLTZMANN_CONSTANT; 
    static constexpr Real r0 =
      e * e / (me * c * c); 
    static constexpr Real Msun =
      GSL_CONST_CGSM_SOLAR_MASS; 
    static constexpr Real Lsun =
      3.828e33; 
    static constexpr Real amu = 
      GSL_CONST_CGSM_UNIFIED_ATOMIC_MASS; 
#elif defined(CRIPTIC_UNITS_MKS)
 
    // MKS units
    static constexpr Real c
    = GSL_CONST_MKSA_SPEED_OF_LIGHT;  
    static constexpr Real mp
    = GSL_CONST_MKSA_MASS_PROTON;     
    static constexpr Real me
    = GSL_CONST_MKSA_MASS_ELECTRON;   
    static constexpr Real sigmaT
    = GSL_CONST_MKSA_THOMSON_CROSS_SECTION;  
    static constexpr Real kB =
      GSL_CONST_MKSA_BOLTZMANN;         
    static constexpr Real h =
      GSL_CONST_MKSA_PLANCKS_CONSTANT_H;  
    static constexpr Real hBar =
      GSL_CONST_MKSA_PLANCKS_CONSTANT_HBAR;  
    static constexpr Real e =
      GSL_CONST_MKSA_ELECTRON_CHARGE; 
    static constexpr Real a0 =
      GSL_CONST_MKSA_BOHR_RADIUS; 
    static constexpr Real Ryd =
      GSL_CONST_MKSA_RYDBERG; 
    static constexpr Real sigamSB =
      GSL_CONST_MKSA_STEFAN_BOLTZMANN_CONSTANT; 
    static constexpr Real eps0 =
      GSL_CONST_MKSA_VACUUM_PERMITTIVITY; 
    static constexpr Real mu0 =
      GSL_CONST_MKSA_VACUUM_PERMEABILITY;  
    static constexpr Real r0 =
      e * e / (4.0 * M_PI * eps0 * me * c * c); 
    static constexpr Real Msun =
      GSL_CONST_MKSA_SOLAR_MASS; 
    static constexpr Real Lsun =
      3.828e26; 
    static constexpr Real amu = 
      GSL_CONST_MKSA_UNIFIED_ATOMIC_MASS; 
#endif
 
    // Additional constants derived from the above
    static constexpr Real c2 = c*c;        
    static constexpr Real c3 = c*c*c;      
    static constexpr Real e2 = e*e;        
    static constexpr Real e3 = e*e*e;      
    static constexpr Real mp_c = mp * c;   
    static constexpr Real mp_c2 = mp * c2; 
    static constexpr Real me_c = me * c;   
    static constexpr Real me_c2 = me * c2; 
    static constexpr Real alpha = 
      GSL_CONST_NUM_FINE_STRUCTURE; 
    static constexpr Real aR = 4.0 * sigmaSB / c; 
    // Particle masses and related quantities -- from CODATA 2018
    static constexpr Real mPi0 = 134.9758 * units::MeV / c2; 
    static constexpr Real mPic = 139.57061 * units::MeV / c2; 
    static constexpr Real mMu = 105.6583755 * units::MeV / c2; 
    static constexpr Real Tpth = 2 * mPi0 * c2 +
      mPi0 * c2 * mPi0 * c2 /
      (2 * mp_c2); 
    static constexpr Real rMuPi = mMu / mPic; 
    // Masses of astrophysically-important light isotopes -- all taken from
    // https://physics.nist.gov/cgi-bin/Compositions/stand_alone.pl. Note
    // the namine convention for isotopes: mXn is the mass of the isotope of
    // element X with atomic number n. When one isotope is overwhelmingly
    // dominant, for example carbon 12, mX without the trailing number gives
    // the mass of the dominant isotope. The only except to this is for
    // hydrogen, where the mass of deuterium is mDeut -- this is to avoid
    // the confusion that would result from calling it mH2, which could be
    // either the mass of deuterium, or of the H2 molecule.
    static constexpr Real mH = 1.007825032239 * amu; 
    static constexpr Real mDeut = 2.0141017781212 * amu; 
    static constexpr Real mHe3 = 3.016029320125 * amu; 
    static constexpr Real mHe4 = 4.002603254136 * amu; 
    static constexpr Real mHe = mHe4; 
    static constexpr Real mLi6 = 6.015122887416 * amu; 
    static constexpr Real mLi7 = 7.016003436645 * amu; 
    static constexpr Real mBe = 9.01218306582 * amu; 
    static constexpr Real mB10 = 10.0129369541 * amu; 
    static constexpr Real mB11 = 11.0093053645 * amu; 
    static constexpr Real mC12 = 12.0 * amu; 
    static constexpr Real mC13 = 13.0033548350723 * amu; 
    static constexpr Real mC14 = 14.003241988440 * amu; 
    static constexpr Real mC = mC12; 
    static constexpr Real mN14 = 14.0030740044320 * amu; 
    static constexpr Real mN15 = 15.0001088988864 * amu; 
    static constexpr Real mN = mN14;
    static constexpr Real mO16 = 15.9949146195717 * amu; 
    static constexpr Real mO17 = 16.9991317565069 * amu; 
    static constexpr Real mO18 = 17.9991596128676 * amu; 
    static constexpr Real mO = mO16; 
    // Nuclear binding energies -- in units of mp c²
    static constexpr Real BE_carbon12 = 92.16175014 * units::MeV / constants::mp_c2;
    static constexpr Real BE_boron10 = 64.7507 * units::MeV / constants::mp_c2;   
    static constexpr Real BE_nitrogen14 = 104.658 * units::MeV / constants::mp_c2;
    static constexpr Real BE_boron11 = 76.204 * units::MeV / constants::mp_c2;
    static constexpr Real BE_oxygen16 = 127.619 * units::MeV / constants::mp_c2;
  }
}
 
 
#endif
// _CONSTANTS_H_