yapanuwan
/
BergEC-jl
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Gaussian potential implemented

This commit is contained in:
Nuwan Yapa 2024-02-21 19:32:11 -05:00
parent da91031577
commit 6b7b759979
1 changed files with 5 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
ho_basis.jl
View File

@ -1,7 +1,11 @@
using NuclearToolkit
using SpecialFunctions
# Gaussian potentials in HO space
gl(R, l, n) = throw("unimplemented")
N_nl(n, l) = sqrt(2*factorial(n) / gamma(l + n + 3/2)^3)
prefactor(n, l, k) = (-1)^k * binomial(n + l + 1/2, n - k) / factorial(k)
Talmi(l, R, k1, k2) = (1/2) / (1 + 1/R^2)^(3/2 + l + k1 + k2) * gamma(3/2 + l + k1 + k2)
V_Gaussian(R, l, n1, n2) = N_nl(n1, l) * N_nl(n2, l) * sum([prefactor(n1, l, k1) * prefactor(n2, l, k2) * Talmi(l, R, k1, k2) for (k1, k2) in Iterators.product(0:n1, 0:n2)])
function get_sp_basis(E_max)
Es = Int[]