From 3bd5e4f53bd9432d86cfd5bc2f671307f38b598d Mon Sep 17 00:00:00 2001
From: Nuwan Yapa <yapanuwan@gmail.com>
Date: Thu, 18 Dec 2025 16:30:38 -0500
Subject: [PATCH] Center-of-mass correction

---
 NuclearRMF.jl | 5 ++++-
 1 file changed, 4 insertions(+), 1 deletion(-)

diff --git a/NuclearRMF.jl b/NuclearRMF.jl
index ffd91b5..5d3f331 100644
--- a/NuclearRMF.jl
+++ b/NuclearRMF.jl
@@ -3,7 +3,10 @@ include("nucleons.jl")
 include("mesons.jl")
 
 "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"
 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)))