from imports import *
from central_problem import CentralProblem
@sp.vectorize
def V(r, l, s):
return - 1. / r
problem = CentralProblem("Hydrogen atom")
problem.potential = V
problem.mass = 1
problem.omega_interval = (0, 2)
from imports import *
from central_problem import CentralProblem
@sp.vectorize
def V(r, l, j):
"""Woods-Saxon potential for the He5 nucleus"""
#V0 = -70 # MeV
Vso = -7.5 # MeV
r0 = 2. # fm
d = .65 # fm
f = 1/(1 + sp.e**((r - r0)/d))
spin_orbit = .5*(j*(j + 1) - l*(l + 1) - .75)
return f * (problem.V0 - 4*Vso*spin_orbit*(f - 1) / (d * r))
problem = CentralProblem("Helium-5")
problem.potential = V
problem.mass = 0.019272
problem.eV_factor = 1e6
problem.omega_interval = (0, 50)
problem.V0 = -70 # MeV
def plotV():
x = sp.linspace(0.1, 10, 100)
plt.plot(x, V(x, 0, 0.5), 'k', x, V(x, 1, 1.5), 'r', x, V(x, 1, .5), 'b')
# plt.axis([7, 10, -2, 2])
