From 0684059359238c8a328ddbd2467c1b66de056fa3 Mon Sep 17 00:00:00 2001
From: Nuwan Yapa <yapanuwan@gmail.com>
Date: Fri, 22 Nov 2024 15:30:30 -0500
Subject: [PATCH] Constants taken from Hartree.f code

---
 mesons.jl | 23 +++++++++++++----------
 1 file changed, 13 insertions(+), 10 deletions(-)

diff --git a/mesons.jl b/mesons.jl
index 7842f98..9c73d44 100644
--- a/mesons.jl
+++ b/mesons.jl
@@ -2,17 +2,20 @@ using DifferentialEquations
 
 const ħc = 197.327 # ħc in MeVfm
 
-# Values taken from B. G. Todd and J. Piekarewicz, Relativistic Mean-Field Study of Neutron-Rich Nuclei, Phys. Rev. C 67, 044317 (2003)
-const m_ρ = 763
+# Values defined in C. J. Horowitz and J. Piekarewicz, Phys. Rev. Lett. 86, 5647 (2001)
+# Values taken from Hartree.f (FSUGarnet) and
+const m_ρ = 763.0
 const m_ω = 782.5
-const m_s = 508.194
-const g2_s = 104.3871
-const g2_v = 165.5854
-const g2_ρ = 79.6000
-const κ = 3.8599
-const λ = -0.01591
-const ζ = 0.00
-const Λv = 0 # value unknown
+const m_s = 496.939473213388
+const m_γ = 0.000001000 # not defined in paper
+const g2_s = 110.349189097820
+const g2_v = 187.694676506801
+const g2_ρ = 192.927428365698
+const g2_g = 0.091701236 # not defined in paper
+const κ = 3.260178893447
+const λ = -0.003551486718 # LambdaSS
+const ζ = 0.023499504053 # LambdaVV
+const Λv = 0.043376933644 # LambdaVR
 
 const r_reg = 1E-6 # regulator for the centrifugal term in fm