yapanuwan
/
NuclearRMF
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Account for the 2x factor in B0

This commit is contained in:
Nuwan Yapa 2025-02-03 16:52:15 -05:00
parent 20e6b5ba84
commit 88345aab2c
2 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
mesons.jl
View File

@ -57,8 +57,8 @@ function solveMesonWfs(ρ_sp, ρ_vp, ρ_sn, ρ_vn, r_max, divs, iterations=3; in
for _ in 1:iterations for _ in 1:iterations
@. src_Φ0 = g2_s * ((κ/ħc)/2 * (Φ0/ħc)^2 + (λ/6) * (Φ0/ħc)^3) - g2_s * (ρ_sp + ρ_sn) @. src_Φ0 = g2_s * ((κ/ħc)/2 * (Φ0/ħc)^2 + (λ/6) * (Φ0/ħc)^3) - g2_s * (ρ_sp + ρ_sn)
@. src_W0 = g2_v * ((ζ/6) * (W0/ħc)^3 + 2Λv * (B0/ħc)^2 * (W0/ħc)) - g2_v * (ρ_vp + ρ_vn) @. src_W0 = g2_v * ((ζ/6) * (W0/ħc)^3 + 2Λv * (2B0/ħc)^2 * (W0/ħc)) - g2_v * (ρ_vp + ρ_vn)
@. src_B0 = 2Λv * g2_ρ * (W0/ħc)^2 * (B0/ħc) - g2_ρ/4 * (ρ_vp - ρ_vn) @. src_B0 = (2Λv * g2_ρ * (W0/ħc)^2 * (2B0/ħc) - g2_ρ/2 * (ρ_vp - ρ_vn)) / 2
Φ0 .= solveKG(m_s, src_Φ0, r_max) Φ0 .= solveKG(m_s, src_Φ0, r_max)
W0 .= solveKG(m_ω, src_W0, r_max) W0 .= solveKG(m_ω, src_W0, r_max)
B0 .= solveKG(m_ρ, src_B0, r_max) B0 .= solveKG(m_ρ, src_B0, r_max)

2
nucleons.jl
View File

@ -16,7 +16,7 @@ const r_reg = 1E-6 # fm # regulator for the centrifugal term
Reference: P. Giuliani, K. Godbey, E. Bonilla, F. Viens, and J. Piekarewicz, Frontiers in Physics 10, (2023)" Reference: P. Giuliani, K. Godbey, E. Bonilla, F. Viens, and J. Piekarewicz, Frontiers in Physics 10, (2023)"
function dirac!(du, u, (κ, p, E, Φ0, W0, B0, A0), r) function dirac!(du, u, (κ, p, E, Φ0, W0, B0, A0), r)
M = p ? M_p : M_n M = p ? M_p : M_n
common1 = E - W0(r) - (p - 0.5) * B0(r) - p * A0(r) common1 = E - W0(r) - (p - 0.5) * 2B0(r) - p * A0(r)
common2 = M - Φ0(r) common2 = M - Φ0(r)
(g, f) = u (g, f) = u
du[1] = -(κ/(r + r_reg)) * g + (common1 + common2) * f / ħc du[1] = -(κ/(r + r_reg)) * g + (common1 + common2) * f / ħc