andam 0.1.2

A comprehensive Rust library for cosmological calculations, visualizations and simulations
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130

//! Physical and cosmological constants
//!
//! All constants are in SI units unless otherwise specified

use std::f64::consts::PI;

/// Speed of light in vacuum [m/s]
pub const C: f64 = 2.99792458e8;

/// Gravitational constant [m³/(kg·s²)]
pub const G: f64 = 6.67430e-11;

/// Planck constant [J·s]
pub const H_PLANCK: f64 = 6.62607015e-34;

/// Reduced Planck constant (ℏ) [J·s]
pub const HBAR: f64 = 1.054571817e-34;

/// Boltzmann constant [J/K]
pub const K_B: f64 = 1.380649e-23;

/// Stefan-Boltzmann constant [W/(m²·K⁴)]
pub const SIGMA_SB: f64 = 5.670374419e-8;

/// Electron mass \[kg\]
pub const M_E: f64 = 9.1093837015e-31;

/// Proton mass \[kg\]
pub const M_P: f64 = 1.67262192369e-27;

/// Neutron mass \[kg\]
pub const M_N: f64 = 1.67492749804e-27;

/// Atomic mass unit \[kg\]
pub const AMU: f64 = 1.66053906660e-27;

/// Elementary charge \[C\]
pub const E_CHARGE: f64 = 1.602176634e-19;

/// Thomson scattering cross-section \[m²\]
pub const SIGMA_T: f64 = 6.6524587321e-29;

/// Fine structure constant (dimensionless)
pub const ALPHA: f64 = 7.2973525693e-3;

// Astronomical constants

/// Astronomical unit \[m\]
pub const AU: f64 = 1.495978707e11;

/// Parsec \[m\]
pub const PARSEC: f64 = 3.0856775814913673e16;

/// Solar mass \[kg\]
pub const M_SUN: f64 = 1.98847e30;

/// Solar luminosity \[W\]
pub const L_SUN: f64 = 3.828e26;

/// Solar radius \[m\]
pub const R_SUN: f64 = 6.957e8;

// Cosmological constants

/// Hubble constant (Planck 2018) [km/s/Mpc]
pub const H0_PLANCK: f64 = 67.66;

/// Critical density coefficient [kg/(m³·(km/s/Mpc)²)]
pub const RHO_CRIT_COEFF: f64 = 3.0 / (8.0 * PI * G);

/// CMB temperature today \[K\]
pub const T_CMB: f64 = 2.7255;

/// CMB photon number density today [m⁻³]
pub const N_GAMMA: f64 = 4.105e8;

/// Neutrino temperature today \[K\]
pub const T_NU: f64 = 1.95;

// Energy conversions

/// Electron volt to Joules
pub const EV_TO_J: f64 = 1.602176634e-19;

/// Joules to electron volt
pub const J_TO_EV: f64 = 6.241509074e18;

// Time conversions

/// Year in seconds
pub const YEAR: f64 = 3.15576e7;

/// Gigayear in seconds
pub const GYR: f64 = 3.15576e16;

/// Hubble time \[s\]
pub fn hubble_time(h0_km_s_mpc: f64) -> f64 {
    PARSEC * 1e6 / (h0_km_s_mpc * 1e3)
}

/// Critical density today [kg/m³]
pub fn critical_density(h0_km_s_mpc: f64) -> f64 {
    let h_si = h0_km_s_mpc * 1e3 / (PARSEC * 1e6); // Convert to SI
    3.0 * h_si * h_si / (8.0 * PI * G)
}

#[cfg(test)]
mod tests {
    use super::*;
    use approx::assert_relative_eq;

    #[test]
    fn test_speed_of_light() {
        assert_eq!(C, 2.99792458e8);
    }

    #[test]
    fn test_hubble_time() {
        let t_h = hubble_time(70.0);
        let expected = 4.4e17; // approximately 14 Gyr in seconds
        assert_relative_eq!(t_h, expected, epsilon = 1e16);
    }

    #[test]
    fn test_critical_density() {
        let rho_c = critical_density(70.0);
        assert!(rho_c > 0.0);
        assert!(rho_c < 1e-25); // Should be very small
    }
}