deep_causality_physics 0.6.0

Standard library of physics formulas and engineering primitives for DeepCausality.
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

/*
 * SPDX-License-Identifier: MIT
 * Copyright (c) 2023 - 2026. The DeepCausality Authors and Contributors. All Rights Reserved.
 */

// =============================================================================
// Weak Isospin Structure (Matter Field Helper)
// =============================================================================

use crate::{PhysicsError, SIN2_THETA_W};

/// Weak isospin representation for a fermion doublet.
#[derive(Debug, Clone, Copy)]
pub struct WeakIsospin {
    pub isospin: f64,
    pub i3: f64,
    pub hypercharge: f64,
}

impl WeakIsospin {
    pub fn new(isospin: f64, i3: f64, charge: f64) -> Result<Self, PhysicsError> {
        if i3.abs() > isospin + 1e-10 {
            return Err(PhysicsError::DimensionMismatch(format!(
                "I₃ = {} must satisfy |I₃| ≤ I = {}",
                i3, isospin
            )));
        }
        let hypercharge = 2.0 * (charge - i3);
        Ok(Self {
            isospin,
            i3,
            hypercharge,
        })
    }

    pub fn lepton_doublet() -> Self {
        Self {
            isospin: 0.5,
            i3: -0.5,
            hypercharge: -1.0,
        }
    }

    pub fn neutrino() -> Self {
        Self {
            isospin: 0.5,
            i3: 0.5,
            hypercharge: -1.0,
        }
    }

    pub fn up_quark() -> Self {
        Self {
            isospin: 0.5,
            i3: 0.5,
            hypercharge: 1.0 / 3.0,
        }
    }

    pub fn down_quark() -> Self {
        Self {
            isospin: 0.5,
            i3: -0.5,
            hypercharge: 1.0 / 3.0,
        }
    }

    pub fn right_handed(charge: f64) -> Self {
        Self {
            isospin: 0.0,
            i3: 0.0,
            hypercharge: 2.0 * charge,
        }
    }

    pub fn electric_charge(&self) -> f64 {
        self.i3 + self.hypercharge / 2.0
    }

    pub fn vector_coupling(&self) -> f64 {
        let q = self.electric_charge();
        self.i3 - 2.0 * q * SIN2_THETA_W
    }

    pub fn axial_coupling(&self) -> f64 {
        self.i3
    }
    pub fn left_coupling(&self) -> f64 {
        self.vector_coupling() + self.axial_coupling()
    }
    pub fn right_coupling(&self) -> f64 {
        self.vector_coupling() - self.axial_coupling()
    }
}

impl Default for WeakIsospin {
    fn default() -> Self {
        Self::lepton_doublet()
    }
}