34#ifndef HEED_PHYSICAL_CONSTANTS_H
35#define HEED_PHYSICAL_CONSTANTS_H
46static constexpr double Avogadro = 6.0221367e+23 / mole;
52static constexpr double c_light = 2.99792458e+8 * m / s;
53static constexpr double c_squared = c_light * c_light;
60static constexpr double h_Planck = 6.6260755e-34 * joule * s;
61static constexpr double hbar_Planck = h_Planck / twopi;
62static constexpr double hbarc = hbar_Planck * c_light;
63static constexpr double hbarc_squared = hbarc * hbarc;
68static constexpr double electron_charge = -eplus;
69static constexpr double e_squared = eplus * eplus;
75static constexpr double electron_mass_c2 = 0.51099906 * MeV;
76static constexpr double proton_mass_c2 = 938.27231 * MeV;
77static constexpr double neutron_mass_c2 = 939.56563 * MeV;
78static constexpr double amu_c2 = 931.49432 * MeV;
79static constexpr double amu = amu_c2 / c_squared;
85static constexpr double mu0 = 4 * pi * 1.e-7 * henry / m;
86static constexpr double epsilon0 = 1. / (c_squared * mu0);
91static constexpr double elm_coupling = e_squared / (4 * pi * epsilon0);
92static constexpr double fine_structure_const = elm_coupling / hbarc;
93static constexpr double classic_electr_radius = elm_coupling / electron_mass_c2;
94static constexpr double electron_Compton_length = hbarc / electron_mass_c2;
95static constexpr double Bohr_radius =
96 electron_Compton_length / fine_structure_const;
98static constexpr double alpha_rcl2 =
99 fine_structure_const * classic_electr_radius * classic_electr_radius;
101static constexpr double twopi_mc2_rcl2 =
102 twopi * electron_mass_c2 * classic_electr_radius * classic_electr_radius;
106static constexpr double k_Boltzmann = 8.617385e-11 * MeV / kelvin;
111static constexpr double STP_Temperature = 273.15 * kelvin;
112static constexpr double STP_Pressure = 1. * atmosphere;
113static constexpr double kGasThreshold = 10. * mg / cm3;
118static constexpr double universe_mean_density = 1.e-25 * g / cm3;