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

Fixed Green's function issue

This commit is contained in:
Nuwan Yapa 2025-01-29 18:24:09 -05:00
parent 15d23ead77
commit 4f61e8a51f
1 changed files with 9 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
mesons.jl
View File

@ -20,19 +20,23 @@ const Λv = 0.043376933644 # units? # LambdaVR
const r_reg = 1E-9 # fm # regulator for Green's functions
"Green's function for Klein-Gordon equation in natural units"
greensFunctionKG(m, r, rp) = -(1 / m) * rp / (r + r_reg) * sinh(m * min(r, rp)) * exp(-m * max(r, rp))
greensFunctionKG(m, r, rp) = -1 / (m * (r + r_reg) * (rp + r_reg)) * sinh(m * min(r, rp)) * exp(-m * max(r, rp))
"Green's function for Poisson's equation in natural units"
greensFunctionP(r, rp) = -rp^2 / (max(r, rp) + r_reg)
greensFunctionP(r, rp) = -1 / (max(r, rp) + r_reg)
"Solve the Klein-Gordon equation (or Poisson's equation when m=0) for a given a source function (as an array) where
m is the mass of the meson in MeV/c2,
r_max is the r-cutoff in fm."
function solveKG(m, source, r_max)
Δr = r_max / (length(source) - 1)
rs = range(0, r_max; length=length(source))
int_measure = Δr .* rs .^ 2
greensFunction = m == 0 ? greensFunctionP : (r, rp) -> greensFunctionKG(m / ħc, r, rp)
mesh = range(0, r_max; length=length(source))
greenMat = greensFunction.(mesh, transpose(mesh))
return greenMat * source
greenMat = greensFunction.(rs, transpose(rs))
return greenMat * (int_measure .* source)
end
"Iteratively solve meson equations and return the wave functions u(r)=[Φ0(r), W0(r), B0(r), A0(r)] where