diff --git a/mesons.jl b/mesons.jl
index 47c3918..75bf59a 100644
--- a/mesons.jl
+++ b/mesons.jl
@@ -11,7 +11,7 @@ const m_γ = 0.000001000 # MeV/c2 # not defined in paper
 const g2_s = 110.349189097820 # units?
 const g2_v = 187.694676506801 # 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 λ = -0.003551486718 # units? # LambdaSS
 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)"
 function solveMesonWfs(ρ_sp, ρ_vp, ρ_sn, ρ_vn, r_max, divs, iterations=3)
     # A0 doesn't need iterations
-    src_A0 = ρ_vp ./ ħc
+    src_A0 = -sqrt(g2_γ) .* ρ_vp
     A0 = solveKG(0, src_A0, r_max)
 
     (Φ0, W0, B0) = (zeros(1 + divs) for _ in 1:3)