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

Center-of-mass correction

This commit is contained in:
Nuwan Yapa 2025-12-18 16:30:38 -05:00
parent 12acf3297e
commit 3bd5e4f53b
1 changed files with 4 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
NuclearRMF.jl
View File

@ -3,7 +3,10 @@ include("nucleons.jl")
include("mesons.jl") include("mesons.jl")
"Total binding energy of the system" "Total binding energy of the system"
total_E(s::system) = -(nucleon_E(s) + meson_E(s)) function total_E(s::system)
E_cm = (3/4) * 41.0 * A(s)^(-1/3) # Center-of-mass correction [Ring and Schuck, Sec. 2.3]
return -(nucleon_E(s) + meson_E(s)) + E_cm
end
"Normalized Woods-Saxon form used for constructing an initial solution" "Normalized Woods-Saxon form used for constructing an initial solution"
Woods_Saxon(r::Float64; R::Float64=7.0, a::Float64=0.5) = -1 / (8π * a^3 * reli3(-exp(R / a)) * (1 + exp((r - R) / a))) Woods_Saxon(r::Float64; R::Float64=7.0, a::Float64=0.5) = -1 / (8π * a^3 * reli3(-exp(R / a)) * (1 + exp((r - R) / a)))