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 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]) plt.show()