CODEWITHSUNDEEP
wolfram-mathematicaphysicsdifferential-equationspde
Kol
Kol

Reputation: 33

Negative Eigenvalues of the Woods-Saxon Nuclear Potential for Oxygen ND Eigensystem

I am trying to use ND Eigensystem to find the eigenvalues of the Woods-Saxon Potential Schroedinger Equation for Oxygen. I've defined the potential and operator, but I know the functions given are incorrect because when I compare them to the harmonic oscillator solutions they are incorrect.

h = 1/2; V[r_] := (-V0)/(1 + E^((r - R)/a))
\[ScriptCapitalL] = -h^2*u''[r] + V[r]*u[r];

{vals, funs} = NDEigensystem[\[ScriptCapitalL], u[r], {r, 0, 10}, 10];
Show[Plot[Evaluate[h*funs + vals], {r, 0, 10}], 
Plot[V[r], {r, 0, 10}], PlotRange -> {{0, 10}, {-50, 5}}, 
AxesOrigin -> {0, 0}]

Upvotes: 0

Views: 36

Answers (0)

Related Questions