Single-particle energy states for a neutron and a proton are obtained by solving the time-independent Schrödinger equation for the mean-field Woods–Saxon potential along with the spin-orbit term. The wavefunctions are expanded as a linear combination of simple sine-wave basis states, which are eigenfunctions of the infinite spherical-well potential. The requisite algorithm based on matrix diagonalization is implemented in Free Open Source Software (FOSS) Scilab. Initial values for the simulation were taken from model parameters given in the book on Nuclear Structure by Bohr and Mottelson, which were then optimized to obtain the best convergence with the available experimental energy values of various nuclei with magic proton and neutron numbers. The level scheme, as well as the energy values for doubly magic nuclei 82208Pb and 2040Ca, which are obtained using our simulation, is presented in this paper. Finally, energy level diagrams for neutrons and protons with respect to mass number A were arrived at, based on those obtained for various magic nuclei. The evaluation method, which is based on the sine-wave basis, is akin to Fourier analysis. When done with the aid of FOSS Scilab, this technique becomes easily accessible to students at the under-graduate (UG) level and may be studied through small projects.

1.
F.
Marsiglio
, “
The harmonic oscillator in quantum mechanics: A third way
,”
Am. J. Phys.
77
,
253
258
(
2009
).
2.
O. S. K. S.
Sastri
,
Aditi
Sharma
,
Jyoti
Bhardwaj
,
Swapna
Gora
,
Vandana
Sharda
, and
B. P.
Jithin
, “
Numerical solution of square well potential with matrix method using worksheets
,”
Phys. Educ.
36
,
1
14
(
2019
).
3.
T.
Dauphinee
and
F.
Marsiglio
, “
Asymmetric wave functions from tiny perturbations
,”
Am. J. Phys.
83
,
861
866
(
2015
).
4.
R. L.
Pavelich
and
F.
Marsiglio
, “
The Kronig-Penney model extended to arbitrary potentials via numerical matrix mechanics
,”
Am. J. Phys.
83
,
773
780
(
2015
).
5.
Felipe
Le Vot
,
Juan J.
Meléndez
, and
Santos B.
Yuste
, “
Numerical matrix method for quantum periodic potentials
,”
Am. J. Phys.
84
,
426
433
(
2016
).
6.
B. A.
Jugdutt
and
F.
Marsiglio
, “
Solving for three-dimensional central potentials using numerical matrix methods
,”
Am. J. Phys.
81
,
343
350
(
2013
).
7.
Eugene
Meyer
, “
Shell model of nuclear structure
,”
Am. J. Phys.
36
,
250
257
(
1968
).
8.
Safqat
Hussain
,
Saeed
Anwar
, and
Muhammad Iqbal
Majoka
, “
Effect of peer group activity-based learning on students' academic achievement in physics at secondary level
,”
Int. J. Acad. Res.
3
,
940
944
(
2011
), see https://www.researchgate.net/profile/Shafqat_Hussain2/publication/234077357_EFFECT_OF_PEER_GROUP_ACTIVITY-BASED_LEARNING_ONSTUDENTS'_ACADEMIC_ACHIEVEMENT_IN_PHYSICSAT_SECONDARY_LEVEL/links/02bfe50ee2b9e9bebb000000.pdf.
9.
Muhammad
Khan
,
Niaz
Muhammad
,
Maqsood
Ahmed
,
Fiaza
Saeed
, and
Sher Aman
Khan
, “
Impact of activity-based teaching on students' academic achievements in physics at secondary level
,”
Acad. Res. Int.
3
,
146
(
2012
), see http://www.savap.org.pk/journals/ARInt./Vol.3(1)/2012(3.1-19).pdf.
10.
B. P.
Jithin
,
V. V. V.
Satyanarayana
,
S.
Gora
,
O. S. K. S.
Sastri
, and
B. P.
Ajith
, “
Measurement model of an alpha spectrometer for advanced undergraduate laboratories
,”
Phys. Educ.
35
,
1
14
(
2019
), see http://www.physedu.in/pub/2019/PE18-08-518.
11.
Swapna
Gora
,
B. P.
Jithin
,
V. V. V.
Satyanarayana
,
O. S. K. S.
Sastri
, and
B. P.
Ajith
, “
Alpha spectrum of 212 bi source prepared using electrolysis of non-enriched ThNO3 salt
,”
Phys. Educ.
35
,
1
16
(
2019
), see http://www.physedu.in/pub/2019/PE18-07-511.
12.
Kenneth S.
Krane
,
Introductory Nuclear Physics
(
John Wiley & Sons
,
New York
,
1988
).
13.
G. N.
Ghoshal
,
Atomic and Nuclear Physics
(
S. Chand and Co
,
New Delhi
,
1997
).
14.
Samuel S. M.
Wong
,
Introductory Nuclear Physics
(
Prentice Hall
,
New Jersey
,
1990
).
15.
David
Hestenes
, “
Toward a modeling theory of physics instruction
,”
Am. J. Phys.
55
,
440
454
(
1987
).
16.
Kris L. G.
Heyde
,
The Nuclear Shell Model
(
Springer
,
Berlin, Heidelberg
,
1994
).
17.
N.
Gauthier
and
S.
Sherrit
, “
Nuclear potential in the Yukawa model
,”
Am. J. Phys.
59
,
1144
1146
(
1991
).
18.
Ajoy K.
Ghatak
and
S.
Lokanathan
,
Quantum Mechanics: Theory and Applications
(
Macmillan
,
New York
,
2004
).
19.
O.
Bayrak
and
E.
Aciksoz
, “
Corrected analytical solution of the generalized Woods-Saxon potential for arbitrary states
,”
Phys. Scr.
90
,
015302
(
2014
).
20.
Aage
Bohr
and
Ben R.
Mottelson
,
Nuclear Structure
(
World Scientific
,
Singapore
,
1998
).
21.
Robert
Eisberg
and
Robert
Resnick
,
Quantum Physics of Atoms, Molecules, Solids, Nuclei, and Particles
, 2nd ed. (
John Wiley & Sons-VCH
,
New York
,
1985
).
22.
See <http://www.chemistry.wustl.edu/~dgs/chem436/Lectures/Chem-436-Lecture_notes.pdf> for solving Woods-Saxon potential using Runge-Kutta method.
23.
N.
Schwierz
,
I.
Wiedenhover
, and
A.
Volya
, “
Parameterization of the Woods-Saxon potential for shell-model calculations
,” preprint arXiv:0709.3525.
24.
M. R.
Pahlavani
and
S. A.
Alavi
, “
Solutions of Woods-Saxon potential with spin-orbit and centrifugal terms through Nikiforov-Uvarov method
,”
Commun. Theor. Phys.
58
,
739
743
(
2012
).
26.
S. B.
Patel
,
Nuclear Physics: An Introduction
(
New Age International
,
Chennai
,
1991
).
27.
H.
Grawe
, “
Shell model from a practitioner's point of view
,”
The Euroschool Lectures on Physics with Exotic Beams
(
Springer
,
Berlin, Heidelberg
,
2004
), pp.
33
75
.
AAPT members receive access to the American Journal of Physics and The Physics Teacher as a member benefit. To learn more about this member benefit and becoming an AAPT member, visit the Joining AAPT page.