**Journal of Modern Physics**

Vol.07 No.03(2016), Article ID:63482,14 pages

10.4236/jmp.2016.73029

Problems with Translational Invariance of Three-Particle Systems

Algirdas Deveikis

Department of Applied Informatics, Vytautas Magnus University, Kaunas, Lithuania

Copyright © 2016 by author and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

Received 6 January 2016; accepted 14 February 2016; published 17 February 2016

ABSTRACT

A new method for calculation of non-relativistic energy spectrum of Coulomb three-body systems with two identical particles has been developed. The novelty of the method is the introduction of an expansion of the wave function on harmonic oscillator (HO) functions with different sizes in the Jacobi coordinates instead of only one unique size parameter in the traditional approach. The method presented obeys the principles of antisymmetry and translational invariance. The theoretical formulation has been illustrated by evaluation of ground state energies of a number of Coulomb three-body systems with two identical particles for zero HO excitation energy. The analytical solution of this problem in case of only one size parameter has been derived. The obtained results show significant advantage of the base with different sizes over the traditional approach for investigation of the bound state problem of quantum systems.

**Keywords:**

Harmonic Oscillator Basis, Binding Energy, Variational Method, Quantum Systems with Coulomb Potential

1. Introduction

Obviously, the theoretical description of quantum systems should obey the principles of antisymmetry and translational invariance. Taking into account the principle of translational invariance is especially important for a few body systems. However, it is well known that this approach has to cope with the formidable calculation difficulties. A description of the many-particle wave functions in terms of the intrinsic (Jacobi) coordinates turns out to be unfeasible because of the antisymmetrization requirements of the Pauli principle. As a consequence, the best-known methods for an atoms and molecules description, such as the Shell Model or the Hartree-Fock Self-Consistent Field Method [1] , are employing only the non-translationally invariant wave-functions that are constructed by antisymmetrization procedure in single particle variables. Since electron and nuclei masses ratio is approximately equal to 10^{−}^{4}, this certainly restricts the precision of such kind of calculations. The problem is even more severe in atomic nuclei and hadrons calculations where neglecting of translational invariance will produce the unacceptable error of order 1/A, where A is baryon charge of corresponding nucleus. On the other hand, the methods that consistently line the principles of antisymmetrization and translational invariance are not yet able to perform calculations for systems with larger particle number. The well known representatives of these kinds of methods are based on the Faddeev [2] , Faddeev-Jakubovsky [3] , and functional-differential equations [4] or on the expansions in a large basis of harmonic oscillator functions (HO) [5] -[11] . So, the further progress in theoretical description of many-particle quantum systems should be supported by revision of traditional and development of new significantly more effective calculation methods. It should be pointed that there exist a wide variety of the methods that start with antisymmetrization in single particle variables and only then project out the excited states of the centrum of mass motion, thus ensuring the translational invariance of the calculations (see, e.g., [12] and references therein). However, this approach may cope with enormous technical problems since the projection procedure should be performed in very large single particle bases even for systems with small number of particles. This problem can be solved only by using the HO basis which is known of its exceptional property to ensure translational invariance in finite subspaces (i.e. complete Nℏω spaces) formed with transition to intrinsic coordinates. The distinct advantage of intrinsic coordinates is that they may be introduced with individual size parameters (oscillator lengths). In general, the number of different size parameters could amount up to the number of intrinsic coordinates. In this respect, the flaw of single particle bases is that since the identity of the particles constituting the system there must be used only one unique oscillator length.

The purpose of the present paper is to consider ab initio approach based on an expansion of the wave function on HO functions with different sizes in the Jacobi coordinates. The method under preposition revisits and extends the traditional HO basis expansion methods that are known of the employment the same oscillator length for all Jacobi coordinates. Introduction of the larger number of variational size parameters may produce intrinsic coordinates that are more adapted to better grasp the physical contents of the quantum systems. The proposed approach is implemented in a new HO basis expansion method for handling the non-relativistic energy spectrum calculation problem of Coulomb three-body systems with two identical particles. The three-particle systems are very abundant in the nature and of great interest in many fields of physics: molecular, atomic, nuclear and hadronic physics among others. The internal motions of particles in such systems are particularly often subjected to Coulomb interaction. Three-body Coulomb systems with two identical particles include the two electron atoms/ions, a variety of diatomic molecular ions, e.g. hydrogen and its isotops, as well as the exotic mesonic systems or positronium ion. It should be noted, that there exist a lot of different highly accurate approximations of wave functions and energy levels for three-body Coulomb systems states that are constructed by exponential expansion based variational approach (see, e.g., [13] and references therein). However, in spite of very accurate results, the wave functions in these types of calculations are constructed mostly with the intuitions. This circumstance severe restricts the perspectives of this approach in application to systems with larger number of particles. A serious advantage of the approach proposed in this paper is that it can be straightforward extended to the systems containing more particles. It should be stressed the universality of this approach, i.e. the calculation of the ground state and excited states may be performed by the same computational procedure. This is not the case for most of other methods. Although it is possible to deal with all types of forces in this formalism, we consider in this paper only Coulomb interactions. It is well known that the treatment of Coulomb potential in the harmonic oscillator basis is a hard computational problem due to the slow convergence of calculation results. So, the investigation of this case may provide sensitive test of the efficiency of the method and give valuable insight of its applicability in the solution of more complicated problems, e.g., systems with a larger number of particles, non-central potentials containing a strong attraction at short distances, and so on. The application of the theoretical formulation has been illustrated by calculation of ground state energies of a number of Coulomb three-body systems with two identical particles in zeroth order approach (i.e. zero HO excitation energy). The results obtained in basis with different sizes are compared with the ones calculated in traditional basis with the same oscillator length for each Jacobi coordinate.

The paper is organized in the following way. In Section 2 we present the basic equations of our formalism which lead to the compact expression for the intrinsic Hamiltonian matrix elements of the problem. In Section 3 we calculate the ground state energies of Coulomb three-body systems with two identical particles for zero HO excitation energy. The energies are compared to the values obtained in the traditional approach, when available. Section 4 is devoted for conclusions.

2. Coulomb Three-Body Systems with Two Identical Particles

The Hamiltonian of three-particle system, with masses and, and charges equal and interacting with Coulomb potential in laboratory reference frame can be written as sum of single-particle and two-particle operators

(2.1)

where α is fine structure constant, and the i-th particle radius vector. Let us assume that two (of the three) particles are identical. In order to get attractive interaction, one of the particles must be positive and the two others negative, or vice-versa. The two identical particles obviously have the equal masses. Therefore, it may be introduced the convenient Jacobi coordinates that not depend on the masses:

(2.2)

The reduced masses that correspond to the introduced Jacobi coordinates are defined by

(2.3)

The denominators of the expressions of the Coulomb interactions between the first and two identical particles in the Jacobi coordinates acquire the form:

(2.4)

So, the intrinsic Hamiltonian of the considered three-particle system in the Jacobi coordinates will be

(2.5)

Since the identity of two particles the mean values of two terms in the intrinsic Hamiltonian corresponding to the Coulomb interactions between the first and two identical particles are equal. Therefore the intrinsic Hamiltonian may be replaced by the equivalent one with only one doubled term of interaction:

(2.6)

In this paper, we introduce different sizes in the Jacobi coordinates instead of only one unique size parameter in the traditional approach. Therefore, in our method the standard procedure of transformation to dimensionless Jacobi coordinates takes the form:

(2.7)

Here are the dimensionless harmonic oscillator length parameters, one for each Jacobi coordinate (indicated with index) that implies a unique size for HO functions

(2.8)

where are the parameters with dimension of angular frequency and is the index of the reduced masses (2.3). The matrix of the intrinsic Hamiltonian will be calculated with the basis functions constructed of single-particle HO functions depending on the dimensionless Jacobi coordinates. In these coordinates the mean values of the kinetic energy terms of intrinsic Hamiltonian coincide with the mean values of the corresponding kinetic energy terms of single-particle HO Hamiltonian

(2.9)

The eigenfunctions of this Hamiltonian are orthogonal and normalized radial harmonic-oscillator wave functions that in dimensionless Jacobi coordinates have the form

(2.10)

here the quantum numbers for dimensionless Jacobi coordinates are displayed: is the number of harmonic oscillator quanta; is the principal oscillator quantum number; is the orbital angular momenta; is the Laguerre polynomial, and is the gamma function [14] . The intrinsic Hamiltonian in the dimensionless Jacobi coordinates takes the form

(2.11)

Let us rewrite the intrinsic Hamiltonian (2.11) in terms of dimensionless energy quantities. The greatest simplification may be achieved by introducing energy units. In this case, two new variational parameters are

(2.12)

These new variational parameters will replace the following ratios in the dimensionless expression of the intrinsic Hamiltonian

(2.13)

and

(2.14)

here.

Further simplification of the expression of dimensionless intrinsic Hamiltonian may be achieved taking use of the expression

. (2.15)

Finally, the dimensionless intrinsic Hamiltonian takes the form

(2.16)

Here, no assumptions are made that some parts of the Hamiltonian are negligible in comparison to others. It should be noted, that since the product is the constant, the mean values of the intrinsic Hamiltonian (2.11) and dimensionless intrinsic Hamiltonian (2.16) will have variational minimum at the same values of the variational parameters and.

The matrix of the obtained intrinsic Hamiltonian is calculated on the bound angular moment harmonic oscillator functions that can be expressed in terms of the first particle and two identical particle functions utilizing the angular momentum Clebsch-Gordan coefficients

(2.17)

here, and are total angular moment a quantum numbers of the states of the first particle, two identical particles and the three-particle system correspondingly; m is the magnetic quantum number of corresponding angular momentum; is spin variables of the i-th particle. Vector coupling of angular moment a is denoted by. The function depending on the first Jacobi coordinate is composed of angular momentum and spin functions

(2.18)

where is a spherical harmonic, is a spin-1/2 function in the spin space, is the spin quantum number of the first particle, and are the magnetic quantum numbers of orbital momentum and spin for the dimensionless Jacobi coordinate (associated with the first particle). The same self- explanatory structure have the wave function depending on the second Jacobi coordinate:

(2.19)

here and we use the simplified notation for the number of HO quanta and orbital angular moment a quantum numbers, i.e.. The spin function of two identical particles have the form

(2.20)

Let us denote the introduced vector coupled basis functions of three-particle system by ket vector

(2.21)

Let us list the conditions that should be imposed on the quantum numbers in these functions. At first, the wave functions (2.21) should be antisymmetric under permutation of two identical fermions. This imposes the well-known condition: should be an even number. At second, the total spin of the two identical particles with spin of 1/2 could take on only two values:. At third, the total angular momentum, is conserved, so its quantum numbers may take on only the integer values within the range:. Finally, the parity conservation imposes the condition:, here is the total parity of the three- particle wave function (the same for bra and ket functions in the Hamiltonian matrix element). As is usually accepted, the is related to the maximum number of total excitation oscillator quanta by formula.

The expectation values of kinetic energy operators of Hamiltonian (2.16) after application of the Wigner- Eckart theorem, may be expressed by reduced matrix element

(2.22)

represented as double-bar matrix element. The Clebsh-Gordan coefficient in this expression is Kronecker symbol for total angular momentum and its magnetic quantum numbers of the basis. The reduced matrix element of operator acting only on coordinates of first particle may be expressed in terms of 6j coefficient and the reduced matrix element between the single-particle functions depending only on coordinates of the first particle [15]

(2.23)

Since the considered Hamiltonian is a scalar in angular-momentum space, the tensorial rank of the kinetic energy operator is zero. After application of Equation (2.23), the expression of the reduced matrix element of the first kinetic energy term takes the form

(2.24)

Due to zero value of one of the parameters of the 6j coefficient, it may be significantly simplified [14] , therefore, the matrix element of the first kinetic energy term may be written as

(2.25)

Applying the same formula (2.23) for orbital momentum we may further simplify the reduced matrix element

(2.26)

After returning from the reduced to the ordinary matrix element

(2.27)

the matrix element of the first kinetic energy term assumes the form

(2.28)

The matrix element of the first kinetic energy operator in the HO basis depending on the coordinate:

(2.29)

The second kinetic energy term of Hamiltonian (2.16) may be evaluated similar way using the formula that expresses the reduced matrix element of operator acting only on coordinates of two identical particles of the three-particle coupled state in terms of 6j coefficient and the reduced matrix element between the single-particle functions depending only on coordinate [15] . Then for the matrix element of the second kinetic energy term, between the coupled three-particle harmonic oscillator functions, we find that

(2.30)

The single-particle Coulomb term as well as the second kinetic energy term of the Hamiltonian (2.16) depends on the same Jacobi coordinate, so its angular part will be calculated analogously. Then for the matrix element of the single-particle Coulomb term, between the coupled three-particle harmonic oscillator functions

(2.31)

The matrix element of the two-particle Coulomb term assumes the form

(2.32)

We calculate this matrix element by means of Talmi-Moshinsky brackets. In order to apply the Talmi- Moshinsky the moments of the coupled three-particle harmonic oscillator functions should be recoupled in the following form [15]

(2.33)

So, the reduced matrix element of the two-particle Coulomb term actually becomes the expansion in terms of the transformation matrices and reduced matrix elements with new momentum coupling scheme

(2.34)

The reduced matrix element of the two-particle Coulomb term on the functions with new coupling scheme may be further simplified

(2.35)

The expectation value of the two-particle Coulomb operator from Equation (2.16) can be obtained by expressing them in a single-particle form. This may be accomplished by means of new Jacobi coordinates. The orthogonal transformation to that Jacobi coordinates is

(2.36)

The same transformation to Jacobi coordinates should be accomplished for the wave functions. By means of Talmi-Moshinsky brackets the two particle harmonic oscillator functions may be expanded in terms of vector coupled products of the single Jacobi variable functions [16]

(2.37)

here and are oscillator quanta and angular momentum quantum numbers, associated with the center- of-mass motion, and are corresponding quantum numbers of relative motion. The Jacobi coordinate with the zero index is proportional to the center-of-mass coordinate, and Jacobi coordinate with the index equal to one is proportional to their relative coordinate. The various quantum numbers appearing in the Talmi- Moshinsky brackets are constrained by triangular inequalities, by parity conservation and by conservation of the number of quanta. The two-particle Coulomb operator will depend only on the one variable in the new basis if the parameter d of the Talmi-Moshinsky brackets takes the form

(2.38)

Then, in the denominator of the two-particle Coulomb operator may be expressed

(2.39)

here we introduce for brevity the factor C which is defined

(2.40)

Thus for the reduced matrix element of the two-particle term of the Hamiltonian (2.16) one obtains

(2.41)

The reduced matrix element in Equation (2.41) may be simplified similarly as the kinetic energy terms

(2.42)

here the radial harmonic-oscillator wave functions have the same form (2.10) and the integral may be obtained most easily by direct numerical integration. Then for the reduced matrix element of the two-particle Coulomb term, between the functions with coupled orbital angular momenta

(2.43)

Finally, the matrix element of the two-particle Coulomb term between the initial coupled three-particle harmonic oscillator functions takes the form

(2.44)

here, the summation parameters are under restrictions:

;;;

m_{j}―may take any allowed value for given j;; and are even numbers; (takes only two values); and;. In the expansion (2.44), we always consider all the basis states (2.37) with a number of quanta less or equal to a given number. This is the well-known prescription for HO bases to treat correctly the Pauli principle.

3. Ground State Energies Evaluation

A new method for calculation of non-relativistic energy spectrum of Coulomb three-body systems with two identical particles has been applied for calculation of ground state energies in zeroth order approach (with zero HO excitation energy). In this case, the Hamiltonian matrix has only one matrix element, therefore the application of variation principle for ground state energies calculations is significantly simplified. The traditional approach based on HO basis considers that the harmonic oscillator wave function have the same size (or the same scale) in all Jacobi coordinates. The proposed approach transforms to the traditional one with the simple substitution. In this case, the zeroth order matrix element of the dimensionless intrinsic Hamiltonian (2.16) takes the form

(3.1)

The minimization of this function may give an analytical expression for variable parameter, which is chosen to minimize the zeroth order ground-state eigenvalue

. (3.2)

So, for one variable parameter the considered variation eigenvalue problem could be solved analytically, that significantly simplify the analysis of bound state problem. The dependence of the zeroth order ground-state eigenvalue (with one variable parameter) of the dimensionless intrinsic Hamiltonian on the mass ratio of the particles and charge of the non-identical particle is illustrated in Figure 1.

Figure 1 shows that the ground-state eigenvalues have pronounced minima at some values of mass ratio. As charge of the non-identical particle increases, these ground-state eigenvalues lay lover. However, at large values of mass ratio of three-body Coulomb systems with two identical particles their binding energies calculated within the framework of the traditional approach become positive. Figure 1 displays calculation results only for electrons as identical particles, however the dependences of binding energies for other sort of identical particles with the same mass ratio will differ only with scaling factor. Therefore, the result that the binding energy of considered systems calculated within the framework of the traditional approach becomes positive at large values of mass ratio is quite general one. This is in strict contradiction to the well-established result that all Coulomb three-body systems consisting of two identical particles have at least one stable bound state against dissociation [17] . Some characteristic of the obtained dependencies of the zeroth order ground-state eigenvalues calculated within the framework of the traditional approach for the specified charges (with) are presented in Table 1. The first column tabulates the values of mass ratio up to which there exists at least one bound state for the considered three-body systems. In the following two rows the minima of the zeroth order ground-state eigenvalues are given as a function of the mass ratio. The last column contain the limit of the minima of the zeroth order ground-state eigenvalues when. Thus, for we have the “atomic” limit consisting of two light identical particles and a heavy one, on the other hand, for there is the “molecular” limit consisting of two identical heavy particles and a light one. So, the mass ratio represents the degree of adiabaticity of the considered system. The dependence of the zeroth order ground-state eigenvalue of the Hamiltonian (3.1) on the mass relation of the particles, charge of the non-identical particle, and the parameter chosen to minimize the ground-state eigenvalue is presented in Table 1.

Figure 1. The dependence of the zeroth order ground-state eigenvalue (in eV) of the Hamiltonian (3.1) on the mass ratio m_{2}/m_{1} of the particles and charge Z_{1} of the non-identical particle when identical particles are electrons. The parameter γ_{min} is chosen to minimize the ground-state eigenvalue. Every curve is labelled by their Z_{1} value and scaled by corresponding factors: 1 (Z_{1} = 1), 5.3 (Z_{1} = 2), 10.8 (Z_{1} = 3), and 16 (Z_{1} = 4).

Table 1. The characteristics of dependence of the zeroth order ground-state eigenvalue (with one variable parameter γ) on mass ratio m_{2}/m_{1} and charge Z_{1} of the non-identical particle when the charge of any of the identical particles is Z_{2} = −1. (m_{2}/m_{1})_{max} is the upper mass ratio limit for which there exists at least the one bound state of the Coulomb three-body systems with two identical particles (electrons) in zeroth order approach; E_{min} is the minimum of the zeroth order ground-state eigenvalue at specified value of mass ratio m_{2}/m_{1}; the last column displays the limit of the zeroth order ground-state eigenvalue for infinitely small mass m_{2} or infinitely large mass m_{1}.

It should be noted that the results presented in Table 1 remain unchanged if the signs of charges Z_{1} and Z_{2} are reversed.

In the presented method the minimization must be done in a two dimensional space (γ_{1} and γ_{2}). To search the minimum of the matrix element of Hamiltonian (2.16) for zero excitation HO energy, we apply our modified Golden section search method [18] alternatively for parameters γ_{1} and γ_{2}. This method is well-known of its stability and is fast enough in this case of one minimised function (computation time takes about few seconds). For the bound states, the parameters γ_{1} and γ_{2} are calculated by minimization of the corresponding zeroth ground- state eigenvalue of the Hamiltonian (2.17). Aiming to demonstrate the possibilities of our method, we have calculated the zero order binding energies for a large set of three-body Coulomb systems with two identical particles. These results are presented in Table 2.

The tabulated in Table 2 calculation results for considered Coulomb three-body systems with two identical particles consist of the following elementary particles: e-electron, e^{+}-positron, μ-and K^{±}-mesons, p-proton, d- deuteron, t-tritium, α-the nucleus of the helium atom, and Li^{+}-the ion of the lithium atom (^{7}Li). We have used the masses of these particles presented in the data tables [19] . In Table 2 all three-body systems are specified by their mass ratio m_{2}/m_{1}. The two columns indicated by presents the zeroth order ground-state eigenvalues calculated within the framework of the traditional approach and the parameter γ chosen to minimize the corresponding ground-state eigenvalues. The second three columns indicated by displays the zeroth order ground-state eigenvalues and values of variable parameters, γ_{1} and γ_{2}, which give the lowest energy of corresponding binding energies. Modern highly sophisticated and accurate variational computations of the three-body bound states allow one to determine up to 40 correct decimal digits in the total energy (see, e.g., [20] and references therein). These very accurate energy values we name numerically “exact” and present in the last column of the Table 2. The analysis of bound state problem of Coulomb systems with two identical particles with the presented method shows its distinct advantages over the traditional approach. At first, the zeroth order ground-state binding energies calculated with different sizes in the Jacobi coordinates are significantly lover than that calculated with only one size parameter for all Jacobi coordinates. At second, the proposed method correctly predicts at least one stable bound state for all Coulomb three-body systems consisting of two identical particles. This is not the case for traditionWal approach which produces unphysical positive eigenvalues in the “molecular” limit, i.e. when mass ratio m_{2}/m_{1} is large. The examples of such systems are epp and μpp which binding energies are not presented in Table 2. The proposed method correctly gives the bounding of ground-state for these systems.

Table 2. The non-relativistic binding (total) energies (in eV) for the ground-states of Coulomb systems with two identical particles. Energies are for values of parameters γ, γ_{1} and γ_{2}, which give the lowest energy for zero excitation HO energy. The “exact” numerical binding energies of ground-states are also indicated.

4. Conclusions

In this paper the new HO basis expansion method for calculation of non-relativistic energy spectrum of three- body systems with two identical particles has been developed. The proposed formalism consistently outlines the principles of antisymmetrization and translational invariance. Compared with the traditional approach, the method is based on an expansion of the wave function on HO functions with different sizes in the Jacobi coordinates. The advantage of this treatment steams from the introduction of more appropriate degrees of freedom that better reflects the basic physics of the system. The efficiency of the method is based on the very fast computation procedure for Talmi-Moshinsky brackets calculations [16] . The computation of many-particle operators matrix elements in HO basis can be significantly simplified using the Talmi transformations. One of the novelties developed in this paper is the calculation of the Hamiltonian matrix elements on the HO functions with two different oscillator lengths by means of Talmi-Moshinsky brackets. This formalism is particularly well suited for very asymmetric systems (for example one light and two heavy particles). The complications of proposed approach as compared to the traditional one are due to the necessity to perform a double minimization on both different sizes in the Jacobi coordinates instead of a single minimization on an unique size parameter. Nevertheless, since this treatment grasps better the physical contents of the system, it may be expected that for a given number of HO quanta, the precision achieved can be increased as compared to the traditional approach.

The method was applied for estimation of ground state energies of a number of Coulomb three-body systems with two identical particles for zero HO excitation energy. The binding of ground state is obtained for all Coulomb three-body systems with two identical particles. It is consistent with well established results [17] . This is not the case in the traditional approach with only one oscillator length where in zeroth order approach the negative energy of ground state for very asymmetric systems cannot be obtained. For considered three-body systems, this result is confirmed by the analysis of the derived compact solution for ground state energy with only one size parameter. The advantages of the developed approach with different sizes in the Jacobi coordinates let us expect that the proposed simple and efficient zeroth order calculation procedure may be successfully generalized to any many-particle quantum system for testing of the existence of bound ground states.

In this paper we focused on the three-body systems where the interaction potential was limited to the Coulomb interaction only. The method is quite general, however, and can be applied to problems if other interaction potentials are used. The proposed approach may be generalized to problems of more than three particles. The addition of another particle introduces one more Jacobi coordinate with corresponding size parameter (for example three size parameters in case of four-body systems). The generalization of this approach may provide consistent and effective treatment of principle of translational invariance in many-particle quantum systems. This could greatly improve the zeroth order calculation results and considerable increase the rate of convergence of HO basis expansion.

Acknowledgements

The author would like to thank Prof. Gintautas P. Kamuntavičius for stimulating this work and his careful reading of this paper.

Cite this paper

AlgirdasDeveikis, (2016) Problems with Translational Invariance of Three-Particle Systems. *Journal of Modern Physics*,**07**,290-303. doi: 10.4236/jmp.2016.73029

References

- 1. Friedrich, H.S. (2005) Theoretical Atomic Physics. 3rd Edition, Spinger, Berlin.
- 2. Faddeev, L.D. (1961) Soviet Physics JETP, 12, 1014-1019.
- 3. Yakubovsky, O.A. (1967) Soviet Journal of Nuclear Physics, 5, 937-951.
- 4. Kamuntavicius, G.P. (1986) Few-Body Systems, 1, 91-109.

http://dx.doi.org/10.1007/BF01277077 - 5. Kamuntavicius, G.P. (1989) Soviet Journal Particles and Nuclei, 20, 109-140.
- 6. Deveikis, A. and Kamuntavicius, G.P. (1995) Lithuanian Journal of Physics, 35, 14-19.
- 7. Deveikis, A. and Kamuntavicius, G.P. (1996) Lithuanian Journal of Physics, 36, 83-95.
- 8. Navrátil, P. and Barrett, B.R. (1996) Physical Review C, 54, 2986-2995.

http://dx.doi.org/10.1103/PhysRevC.54.2986 - 9. Deveikis, A. and Kamuntavicius, G.P. (1997) Lithuanian Journal of Physics, 37, 371-383.
- 10. Navrátil, P. and Barrett, B.R. (1998) Physical Review C, 57, 3119-3128.

http://dx.doi.org/10.1103/PhysRevC.57.3119 - 11. Navrátil, P., Kamuntavicius, G.P. and Barrett, B.R. (2000) Physical Review C, 61, Article ID: 044001.

http://dx.doi.org/10.1103/PhysRevC.61.044001 - 12. Kamuntavicius, G.P. (2014) Journal of Mathematical Physics, 55, Article ID: 042103.

http://dx.doi.org/10.1063/1.4870617 - 13. Ancarani, L.U., Rodriguez, K.V. and Gasaneo, G. (2010) EPJ Web of Conferences, 3, Article ID: 02009.
- 14. Abramowitz, M. and Stegun, I.A. (1972) Handbook of Mathematical Functions. 9th Edition, National Bureau of Standards, Washington DC.
- 15. Varshalovich, D.A., Moskalev, A.N. and Khersonskii, V.K. (1988) Quantum Theory of Angular Momentum. World Scientific, Singapore.

http://dx.doi.org/10.1142/0270 - 16. Kamuntavicius, G.P., Kalinauskas, R.K., Barrett, B.R., Mickevicius, S. and Germanas, D. (2001) Nuclear Physics A, 695, 191-201.

http://dx.doi.org/10.1016/S0375-9474(01)01101-0 - 17. Hill, R.N. (1977) Journal of Mathematical Physics, 18, 2316-2330. http://dx.doi.org/10.1063/1.523241
- 18. Press, W.H., Teukolsky, S.A., Vetterling, W.T. and Flannery, B.P. (2007) Numerical Recipes: The Art of Scientific Computing. 3rd Edition, Cambridge University Press, New York.
- 19. Mohr, P.J., Taylor, B.N. and Newell, D.B. (2010) CODATA Recommended Values of the Fundamental Physical Constants.

http://arxiv.org/abs/1203.5425 - 20. Liverts, E.Z. and Barnea, N. (2013) Computer Physics Communications, 184, 2596-2603.

http://dx.doi.org/10.1016/j.cpc.2013.06.013