Next Article in Journal
Electric Octupole-Dependent Contributions to Optical Binding Energy
Previous Article in Journal
Scalar QED Model for Polarizable Particles in Thermal Equilibrium or in Hyperbolic Motion in Vacuum
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Three-Dimension Calculation for the Scattering Problem for Non-Spherical Potential

by
Pavel M. Krassovitskiy
1,* and
Fedor M. Pen’kov
1,2
1
Institute of Nuclear Physics, Almaty 050032, Kazakhstan
2
Department of Physics and Technology, Kazakh National University, Almaty 050040, Kazakhstan
*
Author to whom correspondence should be addressed.
Submission received: 5 December 2023 / Revised: 17 January 2024 / Accepted: 24 January 2024 / Published: 6 March 2024
(This article belongs to the Section Atomic Physics)

Abstract

:
The interaction of the 238U with a neutron is studied. Correct accounting for the non-spherical shape of the uranium nucleus is in focus. The optical potential is used as a model. It is shown that the spherically symmetric and non-spherical potentials give different scattering patterns, in particular different resonance features of the cross-section. The possibility of using the method as an extension of the particle–rotor model of the nucleus is illustrated.

1. Introduction

Progress in the use of numerical methods results in the need for the fullest possible inclusion of all possible aspects of nuclear interactions. The possibility of accounting for the non-spherical form of nuclear interaction, along with other necessary components, has recently been included in software packages for calculating nuclear reactions [1,2,3,4]. In a number of studies, the calculation of the interaction is based on the use of the time-dependent Hartree–Fock approximation [3] or the calculation of classical trajectories [4]. If the equation for stationary states with realistic deformed potential has been solved, then the coupled channel procedure [5] is strongly important for the description of the scattering at the selected states.
Therefore, on one hand, the non-sphericity effect is incorporated into the calculations, and on the other hand, the calculation of the interaction is based on approximate, sometimes not even quantum mechanical, schemes. New effects arising from scattering by a non-spherical potential within quantum mechanics are considered in the present work.
The mechanism of the “shape of the nucleus” in the few-body problem is well developed [6,7]. However, this area of science, which has received significant results in the description of the nucleon interaction, the simplest elements and some simple cluster nuclei, is still far from being widely used due to obvious difficulties.
The particle–rotor model of the nucleus (PRMN) [5] and its more developed models [8,9,10,11,12] are a trend that is not directly related to the description of scattering but coincides in many respects with this paper. There are a number of nuclei that can be thought of as a rotating axially symmetric quantum rotor model plus a single nucleon. Some of the nuclear states are described within formalism that is practically equivalent to that presented, for example, in Refs. [10,13]. Here, the states of the 11Be and 13C nuclei are studied as 10Be and 12C plus a neutron, respectively. However, the historical evolution of this trend, namely the use of the coupled channels method, led to ambiguous (not quite straightforward) features in the interpretation of the results. In particular, this is a description of levels with a certain total spin, angular momentum, etc. These quantities are not constants of motion in the case of non-spherical fields. In this work, such an interpretation is not applicable. The solution is used as a sum over the constant of motion, i.e., the projection of the angular momentum onto the axis of rotation. The question of the correspondence of the results obtained is open.
Nucleon–nucleus scattering is considered below to be the interaction of a point particle with an axially symmetric field. The interaction of a neutron with the 238U is chosen as a neutral example. It is one of the most widely used reactions in nuclear power. The interaction is described by the optical potential, which is one of the most popular simple ones. The potential parameters are taken from Ref. [14]. The aim of this work was not to describe the properties of the nucleus within the PRMN. However, due to the similarity of formalism, the presented results can also be interpreted as a representation of a new model for describing some properties of the 239U nucleus in the future, as well as other nuclei with similar structures.

2. Materials and Methods

2.1. Equation for a Potential with Axial Symmetry

The equation for the neutron interaction with a nucleus in the center-of-mass system, H 0 + V E Ψ r = 0 , with a zero impact parameter was written for the scattered wave Χ r , Ψ = e i k r + Χ in the following form:
H 0 + V E Χ r = V e i k r ,
where H0 is the free Hamiltonian, V is the interaction potential, E is the interaction energy (in the system center-of-mass), k is the wave vector corresponding to the incident plane wave. The coordinate system is rotated so that the potential does not depend on the azimuthal angle φ of the spherical coordinate system. We assumed that the wave vector lied in the ZOX plane (plane where φ = 0 ), making an angle θ with the OZ axis (polar axis, θ = 0 ). The above equation can be written in the form of an expansion in terms of a set of orthogonal wave functions, which are the eigenfunctions of the l 3 ^ operator:
Ψ = m ψ m   r , θ e i m φ 2 π = 1 2 π e i k r + m χ m e i m φ ,
where l 3 ^ is the operator of the third component of the angular momentum l 3 ^ e i m φ = m e i m φ , m is its eigenvalue. Substituting this expansion into the original equation, multiplying both sides by e i m φ 2 π and integrating over the angle φ , we obtain a set of equations:
H 0 ,   r θ + m 2 r 2 sin 2 θ   + V r , θ E χ m = V r , θ F r , θ , θ , m ,
where F r , θ , θ , m = 1 2 π 0 2 π e i k r cos θ ~ θ , θ , φ   e i m φ d φ , H 0 ,   r θ is the Hamilton operator corresponding to free motion in the two-dimensional space of the half-plane with fixed φ = c o n s t and θ ~ is the angle between the wave vector and the radius vector r .
The function F is expressed in terms of the Bessel function of the first kind. In the selected coordinate system, the following expression was used:
F r , θ , θ , m = e i k r cos θ cos θ + π 2 m J m k r sin θ sin θ .

2.2. Numerical Model

A numerical scheme for solving two-dimensional differential equations in the scattering problem was used. The solution was based on the finite element method, matrix sweep. The result was the scattering amplitude f θ , φ at distances where Ψ r is the asymptotic equivalent to the following expression:
e i k r cos θ ~ + f θ , φ e i k r r .
The two-dimension scattering program was modified. Previously, it was used to calculate the resonant diffraction of molecules, as well as preliminary calculations of scattering by a non-spherical nucleus; basic details can be found in Refs. [15,16,17].
The optical theorem [18] is valid for each angle θ ; in this case, it has the following form:
σ θ = f θ , φ 2 d Ω ~ I m   f θ = θ , φ = 0 ,
where σ θ is the scattering cross-section (with Ω being the solid angle). Its confirmation was a criterion for the accuracy of the calculations.
Assuming that the angle θ can take any values and taking into account the symmetry equation for θ and π θ , the final cross-section of scattering on the axially symmetric potential can be written by averaging over this angle:
σ = 0 π 2 σ θ sin θ d θ C j σ θ j sin θ j Δ θ ,
where C j is the factor for numerical integration.
The described algorithm was applied to scattering of a neutron at a 238U nucleus. The potential of the nucleus–neutron interaction could be considered axially symmetric in a fairly good approximation, and the higher deformation parameters, except β 2 , could be neglected.
For simplicity, we assumed that the shape of the nucleus was a spheroid with the ratio between the semi-major and semi-minor axes ( a b ) / R 0 = 1.06 β 2 , where R 0 is the radius of the nucleus [19]. The sign of the deformation parameter was considered to be positive for the uranium isotope with a mass of 238 [20]. Two models, namely oblate and prolate spheroids, were used in this study.
The optical potential, which was called global by the authors, was taken as the initial one. It showed good agreement between the calculations and the experimental data [14]. The Wood–Saxon potential with an energy-dependent depth is used in this study, the width of the potential well is considered constant. For the present study, the width parameter is considered depending on the polar angle, and the circle is converted to an ellipse on the plane. The parameters of this ellipse were calculated so that the volume of the spheroid obtained by rotating the ellipse was equal to the volume of the ball. The radius of the ball was determined according to the standard formula R 0 = r 0 A 1 3 , with the parameter r 0 derived from Ref. [14] and A being the atomic mass number.

3. Results

The following control and research calculation with the optical potential is performed:
V + i W + U ,
where V is the real potential, and W and U are the imaginary surface and volume absorption potentials taken from Ref. [14], respectively. The results shown in Figure 1 provide a satisfactory agreement between the cross-sections calculated and given in Ref. [14]. Some difference at low energies is apparently explained by the fact that in Ref. [14], the Hauser–Feshbach theory is included in the calculation of the elastic cross-section. Figure 1 shows different contributions from different parts of the potential. The above calculations demonstrate that the main contribution is made by the real potential, as well as by the imaginary volume and surface absorption ones. Further, the spin-orbital ls-interaction (Vls) is not taken into account in the calculations; its contribution does not significantly change the results.
The scattering of a neutron on uranium as on an axially symmetric nucleus is a little more difficult when in accordance with the above algorithm. Individual solutions are found for different angles θ . The examples of such calculations are illustrated in Figure 2.
Figure 3 shows the results of calculating the averaged cross-sections for neutron scattering on 238U.
Two models are presented: spherically symmetric potential and axially symmetric potential. The latter is used in two versions, for positive and negative internal quadrupole moments. The cross-section for scattering at the optical potential without the ls-interaction calculated in the first model is chosen as the control cross-section, which is in good agreement with the experimental one.
All scatter plots at the real potential quite differ from the control one. However, the cross-sections for scattering at axially symmetric potentials do not have pronounced maxima and minima, in contrast to the cross-section for the spherically symmetric potential. In particular, the scattering cross-section shifted by 2.5 bn on the prolate axially symmetric potential in the region above 0.3 MeV satisfactorily fits the control cross-section.
If the imaginary potential is included in the calculation (Figure 4), the cross-section decreases, and the resonances do not contribute, as for the centrally symmetric potential.
When real potentials, as well as full optical potentials, are used, there is a significant difference in the graphs for the models of the prolate (Figure 4, upper) and oblate (Figure 4, lower) spheroids. These calculations provide quite a simple mechanism for determining the deformation sign for some nuclei. It is possible to draw the following conclusion for the 238U and 235U nuclei based only on the available data. The calculations for 238U presented above also correspond quite well to 235U, and there is an insignificant difference in the potential depth (correction ≈ 5%) and width (less than a percent) due to the mass difference. Consequently, scattering by the 235 and 238 isotopes of uranium is described by similar potentials. Proceeding from the same cross-sections, the same sign of β 2 follows. Otherwise, the cross-sections would be significantly different.
The resonances are well seen in the scattering cross-section at the real potential. They come from the interaction of nuclear and centrifugal barriers. The latter has an analogue in the case of axial symmetry:
m 2 r 2 sin 2 θ ,  
This barrier of a more complex shape than in the case of the spherically symmetric potential gives rise to quasi-stationary states. Some of them are interpreted as states of the 239U nucleus within the PRMN. In particular, the bright double resonance at 0.3 MeV can correspond to certain states in the region of 5 MeV, if the methodology used in studies such as that in Ref. [10] is followed.
Figure 5 shows the results of calculating the cross-sections for different m channels at different values of θ . The sum of all of them, m = 0 , ± 1 ± 2 σ m   , results in plots similar to those shown in Figure 2. For example, Figure 5, left, displays the terms for Figure 2, most upper left, and Figure 5, right, displays the terms for Figure 2, most upper right.
It can be seen that individual resonances appear in individual m-channels. Thus, the double resonance at 0.3 MeV is actually the result of the addition of resonances at m = 1 and m = 2 , and the resonance contributes at angles θ ~ π / 4 and completely disappears at θ ~ π / 2 . Resonances for different third projections of the angular momentum appear differently in the averaged cross-section. For example, this is the resonance at E = 1.4 MeV and m = 2 ; its contribution to the total cross-section is large precisely at θ ~ π / 2 . Or vice versa, the contribution of the resonance at E = 4.3 MeV and m = 4 disappears through addition and subsequent averaging. The used matrix sweep method is sensitive enough to resonances [21]; therefore, the proposed method can be effective as an extension of the PRMN. The question of interpreting the results in generally accepted terms, as well as the question of the applicability of such concepts in the problem with the absence of spherical symmetry, remains open.

4. Discussion

A method for resolving nuclear physics problems based on the improved algorithm is developed. This method is available for the solution for scattering in an axially symmetric field as a set of two-dimensional independent fields.
Cross-sections for elastic scattering of a neutron by 238U are found. The cross-sections for centrally symmetric and axially symmetric potentials are significantly different. At the same time, the curves of the cross-sections in the case of axially symmetric real potentials are smoother than the centrally symmetric ones; there are no significant drops, which exist for the centrally symmetric real potential. This gives a hope for describing the elastic scattering process with a decrease in the contribution of imaginary potentials.
Cross-section resonances or quasi-stationary states can be used within an approach similar to the particle-rotor model of the nucleus. An important difference would be in the classification. A new description could be based on the third component of angular moment and optimal angle of position.

Author Contributions

Conceptualization F.M.P.; methodology, P.M.K.; software, P.M.K.; validation F.M.P. and P.M.K.; formal analysis, F.M.P. and P.M.K.; investigation, P.M.K.; resources, P.M.K.; data curation, P.M.K.; writing—original draft preparation, P.M.K.; writing—review and editing, F.M.P. and P.M.K.; visualization, P.M.K.; supervision, F.M.P.; project administration, P.M.K.; funding acquisition, P.M.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the grant No. BR23891530 (Development of complex scientific research in the field of nuclear and radiation physics on the basis of Kazakhstan accelerator complexes) of the Ministry of Energy of the Republic of Kazakhstan.

Data Availability Statement

The data presented in this study are available on request from the authors.

Acknowledgments

The authors are grateful to the participants of the extended seminar of the Low-Energy Nuclear Physics Sector of the JINR (Joint Institute for Nuclear Research, Dubna, Russia) Bogoliubov Laboratory of Theoretical Physics for their fruitful discussion of this work. A significant part of the calculations was performed on the JINR HYBRILIT cluster.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Maruhn, J.A.; Reinhard, P.-G.; Stevenson, P.D.; Umar, A.S. The TDHF code Sky3D. Comput. Phys. Commun. 2014, 185, 2195–2216. [Google Scholar] [CrossRef]
  2. Sekizawa, K. TDHF theory and its Extensions for the multinucleon transfer reaction: A mini review. Front. Phys. 2019, 7, 20. [Google Scholar] [CrossRef]
  3. Umar, A.S.; Oberacker, V.E. Time dependent Hartree-Fock fusion calculations for spherical, deformed systems. Phys. Rev. C 2006, 74, 024606. [Google Scholar] [CrossRef]
  4. Saiko, V.V.; Karpov, A.V. Analysis of multinucleon transfer reactions with spherical and statically deformed nuclei using a Langevin-type approach. Phys. Rev. C 2019, 99, 014613. [Google Scholar] [CrossRef]
  5. Bohr, A.N.; Mottelson, B.R. Nuclear Structure; World Scientific Publishing Company: Singapore, 1998; Volume II. [Google Scholar] [CrossRef]
  6. Elster, C.; Thomas, J.H.; Glöckle, W. Two-body T-matrices without angular momentum decomposition: Energy and momentum dependencies. Few-Body Syst. 1998, 24, 55–79. [Google Scholar] [CrossRef]
  7. Fachruddin, I.; Elster, C.; Glöckle, W. Nucleon-nucleon scattering in a three dimensional approach. Phys. Rev. C 2000, 62, 044002. [Google Scholar] [CrossRef]
  8. Bugrov, V.P.; Kadmensky, S.G. Proton decay of deformed nuclei. Sov. J. Nucl. Phys. 1989, 4, 967–972. [Google Scholar]
  9. Kadmensky, S.G.; Bugrov, V.P. Proton decay and shapes of neutron-deficient nuclei. Phys. At. Nucl. 1996, 59, 399–402. [Google Scholar]
  10. Esbensen, H.; Brown, B.A.; Sagawa, H. Positive parity states in 11Be. Phys. Rev. C 1995, 51, 1274–1279. [Google Scholar] [CrossRef] [PubMed]
  11. Davids, G.N.; Esbensen, H. Decay rate of triaxially deformed proton emitters. Phys. Rev. C 2004, 69, 034314. [Google Scholar] [CrossRef]
  12. Shneidman, T.M.; Adamian, G.G.; Antonenko, N.V.; Jolos, R.V.; Zhou, S.G. Cluster approach to the structure of 240Pu. Phys. Rev. C 2015, 92, 034302. [Google Scholar] [CrossRef]
  13. Ershov, S.N.; Vaagen, J.S.; Zhukov, M. Cluster model with core excitations. The 11Be example. Phys. At. Nucl. 2014, 77, 989–998. [Google Scholar] [CrossRef]
  14. Han, Y.; Xu, Y.; Liang, H.; Guo, H.; Shen, Q. Global phenomenological optical model potential for nucleon-actinide reactions at energies up to 300 MeV. Haiying Liang Phys. Rev. C 2010, 81, 024616. [Google Scholar] [CrossRef]
  15. Krassovitskiy, P.M.; Pen’kov, F.M. Quantum scattering by nonspherical objects. Bull. Russ. Acad. Sci. Phys. 2017, 81, 730–734. [Google Scholar] [CrossRef]
  16. Krassovitskiy, P.M.; Pen’kov, F.M. Asymptotic behavior of solutions in finite difference schemes. Bull. Russ. Acad. Sci. Phys. 2018, 82, 669–673. [Google Scholar] [CrossRef]
  17. Krassovitskiy, P.M.; Pen’kov, F.M. Determination of the scattering amplitude in a limited area of calculation. Bull. Russ. Acad. Sci. Phys. 2019, 83, 420–423. [Google Scholar] [CrossRef]
  18. Landau, L.D.; Lifshits, E.M. Quantum Mechanics. Nonrelativistic Theory; Pergamon Press Ltd.: Oxford, UK, 1991; Available online: https://archive.org/details/l-d-landau-e.-m.-lifshitz-quantum-mechanics-non-vol-3/ (accessed on 20 January 2024).
  19. Nemets, O.F.; Gofman, Y.V. Nuclear Physics Handbook; Naukova Dumka: Kiev, USSR, 1975. (In Russian) [Google Scholar]
  20. Yushkov, A.V. Surface of β(Z,N) Deformation for Nuclei with Z from 2 to 102. Phys. Elem. Part. At. Nucl. (PEPAN) 1993, 24, 348–408. Available online: http://www1.jinr.ru/Archive/Pepan/1993-v24/v-24-2/2.htm (accessed on 20 January 2024). (In Russian).
  21. Krassovitskiy, P.M.; Pen’kov, F.M. Contribution of resonance tunneling of molecule to physical observables. J. Phys. B At. Mol. Opt. Phys. 2014, 47, 225210. [Google Scholar] [CrossRef]
Figure 1. The calculation results of the cross-section of scattering on the spherically symmetric potential as a function of the energy of interaction. 1—the data from Ref. [14]; 2— V + i W + U ; 3— V , 4— V + V l s . See text for details.
Figure 1. The calculation results of the cross-section of scattering on the spherically symmetric potential as a function of the energy of interaction. 1—the data from Ref. [14]; 2— V + i W + U ; 3— V , 4— V + V l s . See text for details.
Physics 06 00024 g001
Figure 2. The calculation results of the cross-section of scattering on the non-spherical potential as a function of the energy of interaction at various polar angles θ as indicated. See text for details.
Figure 2. The calculation results of the cross-section of scattering on the non-spherical potential as a function of the energy of interaction at various polar angles θ as indicated. See text for details.
Physics 06 00024 g002
Figure 3. The averaged cross-sections for neutron scattering on 238U for two models, oblate (1) and prolate (2) spheroids. Scattering by the spherically symmetric potential is also shown, albeit only the real part (3) and full optical potential (4).
Figure 3. The averaged cross-sections for neutron scattering on 238U for two models, oblate (1) and prolate (2) spheroids. Scattering by the spherically symmetric potential is also shown, albeit only the real part (3) and full optical potential (4).
Physics 06 00024 g003
Figure 4. Comparison of the calculation of the cross-sections for scattering at the axially symmetric potential (real part—(1); full optical potential—(2)) and the spherically symmetric potential (3) for the prolate (upper) oblate (lower) spheroids.
Figure 4. Comparison of the calculation of the cross-sections for scattering at the axially symmetric potential (real part—(1); full optical potential—(2)) and the spherically symmetric potential (3) for the prolate (upper) oblate (lower) spheroids.
Physics 06 00024 g004
Figure 5. Calculation of the contribution to the cross-section for various m components (left) for θ = π / 4 and (right) for θ = π / 2 , where the numbers denote m values. See text for details.
Figure 5. Calculation of the contribution to the cross-section for various m components (left) for θ = π / 4 and (right) for θ = π / 2 , where the numbers denote m values. See text for details.
Physics 06 00024 g005
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Krassovitskiy, P.M.; Pen’kov, F.M. Three-Dimension Calculation for the Scattering Problem for Non-Spherical Potential. Physics 2024, 6, 368-375. https://0-doi-org.brum.beds.ac.uk/10.3390/physics6010024

AMA Style

Krassovitskiy PM, Pen’kov FM. Three-Dimension Calculation for the Scattering Problem for Non-Spherical Potential. Physics. 2024; 6(1):368-375. https://0-doi-org.brum.beds.ac.uk/10.3390/physics6010024

Chicago/Turabian Style

Krassovitskiy, Pavel M., and Fedor M. Pen’kov. 2024. "Three-Dimension Calculation for the Scattering Problem for Non-Spherical Potential" Physics 6, no. 1: 368-375. https://0-doi-org.brum.beds.ac.uk/10.3390/physics6010024

Article Metrics

Back to TopTop