From b99cdc98e3901866600f1370c5c9d3ffa2b552c2 Mon Sep 17 00:00:00 2001 From: Nuwan Yapa Date: Fri, 21 Feb 2025 10:40:27 -0500 Subject: [PATCH] Fixed meson_E (dimensions and couplings) --- mesons.jl | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/mesons.jl b/mesons.jl index 35e6045..2eb9ca0 100644 --- a/mesons.jl +++ b/mesons.jl @@ -73,10 +73,10 @@ end function meson_E(s::system) int = 0.0 for (r, Φ0, W0, B0, A0, ρ_sp, ρ_vp, ρ_sn, ρ_vn) in zip(rs(s), s.Φ0, s.W0, s.B0, s.A0, s.ρ_sp, s.ρ_vp, s.ρ_sn, s.ρ_vn) - E_σ = (1/2) * Φ0 * (ρ_sp + ρ_sn) - ((κ_ss/ħc)/12 * (Φ0/ħc)^3 + (λ/24) * (Φ0/ħc)^4) - E_ω = -(1/2) * W0 * (ρ_vp + ρ_vn) + (ζ/24) * (W0/ħc)^4 + E_σ = (1/2) * (Φ0/ħc) * (ρ_sp + ρ_sn) - ((κ_ss/ħc)/12 * (Φ0/ħc)^3 + (λ/24) * (Φ0/ħc)^4) + E_ω = -(1/2) * (W0/ħc) * (ρ_vp + ρ_vn) + (ζ/24) * (W0/ħc)^4 E_ρ = -(1/4) * (2B0/ħc) * (ρ_vp - ρ_vn) - E_γ = -(1/2) * g2_γ * A0 * ρ_vp + E_γ = -(1/2) * (A0/ħc) * ρ_vp E_ωρ = Λv * (W0/ħc)^2 * (2B0/ħc)^2 int += (E_σ + E_ω + E_ρ + E_γ + E_ωρ) * r^2 end