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

Use Hartree.h convention for A0 (still wrong)

This commit is contained in:
Nuwan Yapa 2025-01-29 13:15:31 -05:00
parent 4801e9070d
commit 15d23ead77
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
mesons.jl
View File

@ -11,7 +11,7 @@ const m_γ = 0.000001000 # MeV/c2 # not defined in paper
const g2_s = 110.349189097820 # units? const g2_s = 110.349189097820 # units?
const g2_v = 187.694676506801 # units? const g2_v = 187.694676506801 # units?
const g2_ρ = 192.927428365698 # units? const g2_ρ = 192.927428365698 # units?
const g2_g = 0.091701236 # units? # not defined in paper const g2_γ = 0.091701236 # dimensionless # equal to 4πα
const κ = 3.260178893447 # units? const κ = 3.260178893447 # units?
const λ = -0.003551486718 # units? # LambdaSS const λ = -0.003551486718 # units? # LambdaSS
const ζ = 0.023499504053 # units? # LambdaVV const ζ = 0.023499504053 # units? # LambdaVV
@ -44,7 +44,7 @@ end
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 solveMesonWfs(ρ_sp, ρ_vp, ρ_sn, ρ_vn, r_max, divs, iterations=3) function solveMesonWfs(ρ_sp, ρ_vp, ρ_sn, ρ_vn, r_max, divs, iterations=3)
# A0 doesn't need iterations # A0 doesn't need iterations
src_A0 = ρ_vp ./ ħc src_A0 = -sqrt(g2_γ) .* ρ_vp
A0 = solveKG(0, src_A0, r_max) A0 = solveKG(0, src_A0, r_max)
(Φ0, W0, B0) = (zeros(1 + divs) for _ in 1:3) (Φ0, W0, B0) = (zeros(1 + divs) for _ in 1:3)