Nuclear Science
Volume 2, Issue 1, March 2017, Pages: 1-4

Energy Levels Calculations of 24Al and 25Al Isotopes

Saeed Mohammadi, Banafsheh Nemati Giv, Nafiseh Shayan Shakib

Department of Physics, Payame Noor University, Tehran, Iran

Email address:

(S. Mohammadi)

To cite this article:

Saeed Mohammadi, Banafsheh Nemati Giv, Nafiseh Shayan Shakib. Energy Levels Calculations of 24Al and 25Al Isotopes. Nuclear Science. Vol. 2, No. 1, 2017, pp. 1-4. doi: 10.11648/j.ns.20170201.11

Received: December 10, 2016; Accepted: December 28, 2016; Published: January 20, 2017

Abstract: In this article, the energy levels of Aluminum isotopes 24Al and25Al are calculated using OXBASH shell model code. The calculations were carried out in the SD model space with the W and CW Hamiltonian [1, 2] using the code OXBASH for windows PC [3]. This code which is based on one of the most applicable nuclear models, the shell model, deals with evaluating energy levels in nuclei. Using this code to calculate the energy levels of an isotope, produces several files that contains a set of data. The ground state energy of a nucleus and also the probable energies of excitation levels can be calculated by OXBASH. Programs will be reliable only when results meet experimental procedures. A compilation of SD-shell energy levels calculated with the USD Hamiltonian was published around 1988. A comparison had been made between our results and the available experimental data to test theoretical shell model description of nuclear structure in Aluminum isotopes. The calculated energy spectrum is in good agreement with the available experimental data.

Keywords: Aluminum Isotopes, OXBASH Code, Shell Model Structure, USD Interaction

1. Introduction

Obtaining the nuclear structure and energy levels of nuclei is one of the criteria to improve investigations of nuclei properties. Nuclear models have the property to help us to better understanding of nuclear structure which contains main physical properties of nuclei, and shell-model is one of the most prominent and successful nuclear models [4-11]. This model can be compared with the electron shell model for atoms. As atomic behavior and properties can be described with valance electrons which exist out of a closed shell, similarly, valance nucleons (protons or neutrons) in a nucleus which are placed out of close shells (with magic numbers 2,8,20,28,50,82 and 126) play important roles in determining nuclear properties. Nuclei with magic numbers are very stable and have completely different properties comparing with their neighbors.

Existence of spatial levels is determined by Pauli Exclusion Principle. By knowing nuclear potentials for all nucleons in a nucleus, one can calculate energy levels.

For light nuclei, there are several "standard" effective interactions for the p and SD shells, respectively [12,13]. Analysis of neutron-rich SD nuclei has been of high interest in recent years as they present new aspects of nuclear structure [14]. Traditional shell-model studies have recently received a renewed interest through large scale shell-model computing in no-core calculations for light nuclei, the 1s0d shell, the 1p0f shell and the 3s2d1g7/2shell with the inclusion of the 0h11/2 intruder state as well. It is now therefore fully possible to work to large-scale shell-model data and study the excitation levels for large systems. In these systems, inter core is assumed and space is determined by considering shell gaps. Figure 1 shows the PF and SD model spaces according to shell model theory.

Figure 1. Sd and pf shell model spaces.

2. Theory and Calculation

In order to calculate the nuclear structure properties of both ground and excited states based on the nuclear shell model one needs to have wave functions of those states. These wave functions are obtained by using the shell-model code OXBASH [12]. OXBASH (Oxford- Buenos Aires Shell Model Code) is a powerful computer code to calculate the energy levels of light and medium nuclei. By using it, we can measure the energy levels of the nucleus and compare it with experimental data. One of the most important challenges in nuclear physics is to calculate and measure energy levels of medium and heavy nuclei by using shell model theory. As a result, it is important to check computer codes in this area. There are many codes written for this purpose such as: NutShell, REDSTICK, ANTOINE, OXBASH, etc. In this paper we are using OXBASH code, a code which has been developed in this field for many years.

Experimental studies on the stability of some nuclei show that nuclei with protons or neutrons numbers (2, 8, 20, 50, 82, etc.) are more numerous and stable. In other words, the existence of these magic numbers suggests shell model structure in nuclear physics similar to atomic physics. In this paper, the energy levels of Na and Mg isotopes have been measured using the code OXBASH. The program includes a set of computational code which is based on the ability to measure the energy levels by forming ground state matrices with dimensions up to 2000000 and JT matrix with dimension up to 100,000. The version of this code is 2005-8 which can be installed and used on any operating system without using any other additional software.

Space model SD

Considering the number of valence nucleons for sodium and magnesium isotopes, SD model space is the suitable model for these calculations. This model consists of،و  valence orbitals. The code OXBASH for Windows has been used to calculate the nuclear structure for Aluminum nucleus, by employing the SD (independent charges) and SDPN (depending charges) model space with three effective interactions [3]. The first interaction for the lower part of the SD-shell is Chung-Wildenthal particle interaction (CW), secondly, the Universal SD-shell Hamiltonian (USD interaction). In the third interaction the New Universal SD-shell Hamiltonian (USDAPN) is used. Richter et a1 [15] used this shell model successfully in the p-shell, and fp-shell [16], [17] and Wildenthal [2] and Brown et al[18] in the SD-shell to describe the systematic observed in the spectra and transition strengths. In the present work, we focus our attention on the description of energy levels of SD shell of Aluminum isotopes 24Al and 25Al which have configurations 0d5/2, 1s1/2and0d3/2.

The basis of Hamiltonian in shell model calculations is normalized Green matrix. For SD model space in OXBASH code, 13 different interaction potentials to perform computations on different nuclei is included [3].

The first normalized G matrix was introduced in 1960 which its results were in good agreement with experimental data for A = 18 to A = 38 nuclei [17,18], however, with increasing nucleon numbers, solving this matrix was very difficult [19].In 1976 for the bottom layer of SD model space, the Hamiltonian CWH [3,17] was introduced and in 1984 the potential W was introduced for the entire layer SD [21].

OXBASH code is a computer program that is described with a set of model spaces and interactions to apply in shell model calculations with high dimensions.

In order to use this code one should specify the model space and interaction. In other words after choosing appropriate model space which is chosen considering valance nucleons, this code constructs a set of possible ground states and then makes JT matrix based on linear combination of ground states which give suitable T and J values. Finally by choosing the desirable interaction potential it constructs the Hamiltonian of the problem and carries out the calculations and as a default gives 10 lowest energies.

The applied model space illustrates the orbitals which are considered in calculations using main shells in shell-model theory.

In the second line of approach the two body matrix elements are treated as parameters, and their values are obtained from best fit to experimental data [18].

Brown and coworkers [19] have carried extensive studies of energy level and spectroscopic properties of SD-shell nuclei in terms of a unified Hamiltonian applied in full SD-shell model space. The universal Hamiltonian was obtained from a least square fit of 380 energy data with experimental errors of 0.2MeV or less from 66 nuclei. The USD Hamiltonian is defined by 63 SD-shell two body matrix element and their single particle energies. In more recent work Brown and co-workers have modified USD type Hamiltonian to USDA and USDB based on updated set of binding energy and energy levels of O, F, Ne, Na, Mg and P isotopes [22].

3. Results

The isotopes of Aluminum 24Al and 25Al provide a unique work frame for examining the foundation of SD shell model calculations. The nucleons of the core of these isotopes are 8 protons and 8 neutrons which are in the (1s1/2,0p3/2,0p1/2)J=0,T=0 ground state configuration of 16O and the remaining nucleons are distributed over all possible combinations of the 0d5/2,1s1/2 and 0d3/2 orbits according to Pauli Exclusion Principle. The package of program called "SHELL" was used to generate the One Body Density Matrix Element (OBDME), and the package of program called "LPE" is used to calculate the wave functions and energy levels.

We present here some results concerning ground and excitation energies of the Al isotopes for which recent data has been reported in the literature. Table 1 and figure 2 shows data for 24AL isotope and table 2 and figure 3 shows data for 25Alisotope. E (OXBASH) are energies calculated in this work and E(exp) are experimental data.

Figure 2. Comparison Calculated data with Experimental data for 24Al.

Figure 3. Comparison Calculated data with Experimental data for 25Al.

Table 1. Data for 24Al. All energies are in MeV.

E(exp) E(OXBASH) J/AL24
1.110 1.155 1→J=2
1.130 1.155 2→J=2
2.380 2.378 3→J=2
4.340 4.369 4→J=2
3.060 3.031 5→J=3
4.770 4.736 6→J=3
4.258 4.259 7→J=4
4.770 4.728 8→J=4
1.110 1.101 9→J=5
6.810 6.822 10→J=5
6.810 6.808 11→J=6
8.250 8.205 12→J=7
8.250 8.183 13→J=8
3.350 3.331 14→J=3

4. Conclusions

As seen from all figures, very good agreement is obtained for most of energy levels of Aluminium isotopes, and the ordering of levels is correctly reproduced. Unfortunately, not enough experimental data were available for 24AL; regarding the USD data we can judge that almost all calculations reasonably reproduce the observed level structure.

Table 2. Data for 25Al. All energies are in MeV.

E(exp) E(OXBASH) J/AL25
4.514 4.546 1→J=1/2
4.583 4.546 2→J=1/2
8.026 8.025 3→J=1/2
9.065 9.078 4→J=1/2
3.744 3.715 5→J=3/2
6.288 6.330 6→J=3/2
6.327 6.330 7→J=3/2
7.300 7.344 8→J=3/2
1.306 1.352 9→J=5/2
4.196 4.129 10→J=5/2
5.232 5.200 11→J=5/2
5.232 5.247 12→J=5/2
6.909 6.969 13→J=5/2
2.236 2.237 14→J=7/2
5.232 5.244 15→J=7/2
6.327 6.393 16→J=7/2
6.881 6.883 17→J=7/2
7.022 7.015 18→J=7/2
7.112 7.139 19→J=7/2
3.251 3.254 20→J=9/2
4.026 4.034 21→J=9/2
5.582 5.563 22→J=9/2
5.809 5.809 23→J=9/2
6.645 6.637 24→J=9/2
7.240 7.261 25→J=9/2

In general the best and most complete results are found with the largest model space while calculations in an infinite space are not possible and the computation time increases exponentially with model space size so some truncation is required. Also the interaction used must be appropriate for the model space. The empirical interactions are (usually) better determined for smaller model spaces. The model space in OXBASH is defined by the active valance nucleon orbits and our calculated results are reasonably consistent with experimental data, although the structure of odd-even nuclei is much more complicated than their odd-odd neighbours.


  1. B. H. Wildenthal and W. Chung, in "Mesons in Nuclei", edited by M. Rho and D. Wilkinson (North-Holland, Amsterdam, 1979), Vol. II., p. 722.
  2. B. H. Wildenthal, Prog. Part. Nucl. Phys. 11 (1984).
  3. OXBASH for Windows, B. A. Brown, A. Etchegoyen, N. S. Godwin, W. D. M. Rae, W. A. Richter, W. E. Ormand, E. K. Warburton, J. S. Winfield, L. Zhao and C. H. Zimmerman, MSU_NSCL report number 1289.
  4. Saeed Mohammadi, Sara Sirjani. Energy Levels Calculations of 29S, 30S and 31S Isotopes. American Journal of Modern Physics.SpecialIssue: Many Particle Simulations. Vol. 4, No. 3-1, 2015, pp. 5-9. doi: 10.11648/j.ajmp.s.2015040301.12.
  5. S. Mohammadi, S. Arbab, E. Tavakoli. Energy Levels Calculations of 22-23Na and 24-26Mg Isotopes Using Shell Model Code OXBASH. American Journal of Modern Physics.Special Issue: Many Particle Simulations. Vol. 4, No. 3-1, 2015, pp. 27-31. doi: 10.11648/j.ajmp.s.2015040301.16.
  6. Saeed Mohammadi, Monna Chobbdar. Energy Levels Calculations of 32Cl and 33Cl Isotopes.American Journal of Modern Physics.Special Issue: Many Particle Simulations. Vol. 4, No. 3-1, 2015, pp. 10-14. doi: 10.11648/j.ajmp.s.2015040301.13.
  7. S. Mohammadi, A. Heydarzade, V. Ragheb. Energy Levels Calculations of 28, 30Na and 26, 28Al Isotopes Using Shell Model Code OXBASH, American Journal of Modern Physics. Special Issue: Many Particle Simulations. Vol. 4, No. 3-1, 2015, pp. 32-35. doi: 10.11648/j.ajmp.s.2015040301.17.
  8. S. Mohammadi, F. Bakhshabadi. Calculation of the Energy Levels of Phosphorus Isotopes (A=31 to 35) by Using OXBASH Code. American Journal of Modern Physics.Special Issue: Many Particle Simulations. Vol. 4, No. 3-1, 2015, pp. 15-22. doi: 10.11648/j.ajmp.s.2015040301.14.
  9. Saeed Mohammadi, Maryam Mounesi. Energy Levels Calculations of 20Ne and 21Ne Isotopes. American Journal of Modern Physics. Special Issue: Many Particle Simulations. Vol. 4, No. 3-1, 2015, pp. 36-39. doi: 10.11648/j.ajmp.s.2015040301.18.
  10. Saeed Mohammadi, Hassan Rostam Nezhad. Energy Levels Calculations of 26Al and 29Al Isotopes.American Journal of Modern Physics. Special Issue: Many Particle Simulations. Vol. 4, No. 3-1, 2015, pp. 1-4. doi: 10.11648/j.ajmp.s.2015040301.11.
  11. S. Mohammadi, S. Hesami Rostami, A. Rashed Mohasel and M. Ghamary. Energy Levels Calculations of 36-37-38Ar Isotopes Using Shell Model Code OXBASH. American Journal of Modern Physics. Special Issue: Many Particle Simulations. Vol. 4, No. 3-1, 2015, pp. 23-26.doi: 10.11648/j.ajmp.s.2015040301.15.
  12. Nuclear Physics Graduate Course: Nuclear Structure Notes, Oxford, (1980).
  13. Effective interactions and the nuclear shell-model, D. J. Dean,T. Engeland M. Hjorth-Jensen, M. P. Kartamyshev, E. Osnes, Progress in Particle and Nuclear Physics 53 419–500(2004).
  14. N. A. Smirnova, Shell structure evolution and effective in-medium NN interaction, Chemin du Solarium, BP 120, 33175 Gradignan, France.
  15. W. A. Richter, M. G. van der Merwe, R. E. Julies, B. A. Brown, Nucl. Phys. A 523, 325, (1991).
  16. R. E. Julies, PhD. thesis, University of Stelenbosch, (1990).
  17. M. G. van der Merwe, PhD. thesis, University of Stelenbosch (1992).
  18. B. A. Brown and B. H. Wildenthal, Ann. Rev. of Nucl. Part. Sci. 38, 29 (1988).
  19. B. A. Brown, W. A. Richter, New "USD" Hamiltonians for the sd shell, Physical Review C, 74 (2006).
  20. B. H. Wildenthal, Prog. Part. Nucl. Phys. 11, 5 (1984), and references therein.
  21. M. Homma, B. A. Brown, T. Mizusaki and T. Otsuka, Nucl. Phys. A704, 134c (2002).
  22. B. A. Brown and W. A. Richter, Phys. Rev. C 74, 034315, (2006).

Article Tools
Follow on us
Science Publishing Group
NEW YORK, NY 10018
Tel: (001)347-688-8931