## MAPLE program for modelling hydrogen-like atoms in quantum mechanics with non-negative distribution function

**Authors:**Zorin A.V., Tretyakov N.P.**Issue:**Vol 26, No 4 (2018)**Pages:**343-356**Section:**Mathematical Modeling**URL:**http://journals.rudn.ru/miph/article/view/20224**DOI:**http://dx.doi.org/10.22363/2312-9735-2018-26-4-343-356

#### Abstract

The program is proposed for a realization of the symbolic algorithm based on the quantum mechanics with non-negative probability distribution function (QDF) and for calculations of energy levels for hydrogen-like atoms. The program is written up in the language MAPLE. In the framework of the algorithm an original Maple package for calculations of necessary functions, such as hydrogen wave functions, Sturmian functions and their Fourier-transforms, Clebsch-Gordan coefficients, etc. is proposed. Operators of observables are calculated on the basis of the QDF quantization rule. According to the Ritz method, eigenvalues of Ritz matrices represent spectral values of the quantity under investigation, i.e. energy. As an example, energy levels of hydrogen-like atoms are calculated and compared with experimental data retrieved from the NIST Atomic Spectra Database Levels Data. It turns out that this theory seems to be equivalent to the traditional quantum mechanics in regard to predictions of experimental values. However, the existence of a phase-space probabilistic quantum theory may be an important advance towards the explanation and interpretation of quantum mechanics.

1. Introduction In [1] several computational techniques for investigation of characteristics of atomic structures were proposed. This method goes far beyond the techniques described in [2]. A deficiency of the single-configurational approximation was noted in [3]. D. Layzer [4] suggested the use of multy-configurational approximation restricting to the case of configurations belonging to the same complex (i.e. with the same parities and the same principal quantum numbers). By contrast, Eissner and Nussbaumer [1] consider interactions of configurations belonging to different complexes. Nussbaumer [5] introduced some refinement of the technique for calculation of radial functions as compared with variational approaches (for instance, Thomas-Fermi method). In particular, the procedure includes an effective charge � (as compared with the real one �): � - �1 < � < � which had been introduced in [4]. Calculations of radial functions were used in calculations of oscillators strengths and transitional probabilities. It is worth noting that the convergence of calculational results to experimental data while more and more complexes are taken into account was strictly demonstrated in [5]. The electron correlation effects on the oscillator strengths for resonance transitions in alkaline earth atoms have been studied by many authors. The conventional approach Received 27th October, 2018. consists in improving Hartree-Fock-like wavefunctions by introducing configuration interaction (CI) while using the usual dipole operator for the oscillator strength. This approach was used by Burke et al. [6] for Be, Zare [7] and Weis [8] for Mg, and Friedrich and Trefftz [9] for Ca and Ba. An alternative approach was adopted in [10], which uses an effective operator in addition to the common dipole operator and in so doing keeps a rather simple structure of corresponding wave functions. This method had been previously proposed by Hameed [11, 12] and it works well with alkaline metals and provides slightly worse results when applied to alkaline earth elements. However, the modification of the method proposed in [10], works very well for alkaline earths as well. Relativistic corrections introduced in [13, 14] were included in the extension [15] of the package [5]. It is worth noting that the whole cycle of the above-mentioned investigations has been elaborated in detail by Russian authors as well. The books [16-18] present a good introduction and an excellent review of the problem. Among other things they suggest that radial one-particle functions obtained using semi-empirical potentials, are best suited for calculations of transitional probabilities. Quantum mechanics with non-negative distribution function (QDF) is an alternative physical theory with a developed theoretical formalism [19-28]. However, the weak point in this theory lays in the fact that until the present time no truly extensive investigations of properties of real quantum systems have been conducted. This is due to the much more complex structure of the QDF theory as compared with traditional quantum mechanics. Generations of operators of observables include calculations of complex multidimensional integrals. Furthermore, this leads to rather cumbersome expressions for operators, resulting in the absence of analytical approaches to calculations of spectra. For purposes of effective computer calculations, it is necessary to combine analytical and numerical methods, i.e. creation of a package of programs mating with each other. The evidence for the capacity of the QDF theory to obtain the same results as conventional quantum mechanics would mean an important advance in the interpretation of quantum mechanics as a statistical theory. The present work solves the analytical part of the above-mentioned problem. The MAPLE package is used to generate necessary functions, their Fourier-transforms, to calculate integrals, operators and matrices according to formulae of QDF. Besides, numerical potentialities of MAPLE are exploited as well in calculations of eigenvalues of Ritz matrices. This program has been being developed for several years by us and co-authors, starting with [27]. The computation time raises as the dimension of Ritz matrices does. On the other hand, for small dimensions of matrices, the accuracy of the computations is insufficient. The hardware (processor and memory, for example) available for the computation may restrict either the degree of accuracy or the computation time. As a result of comparisons and tests of various program variants, we came to the conclusion that the dimension of the matrices 91 × 91 is optimal. Thus, this article sums up the many years of work on the creation of a package of programs and poses challenges for future research. 2. Quantum Mechanics with Non-Negative Distribution Function Operators in quantum mechanics with non-negative QDF are defined up to an arbitrary set of integrable functions of configuration space and time �� (�, �), normalized by the condition: ∑︁ ∫︁ � |�� (�, �)|2 d� = 1. (1) By virtue of square integrability the functions �� allow the Fourier transformation ∫︁ � ′ �(�)�(�, �) = (2�n)-� �(� - �, � - �)�n (�-� )��(�, �, �)�(�′, �) d� d� d�d�′, (2) where n is an arbitrary action-dimensioned constant. In what follows it is convenient to introduce a new function of phase space and time, constructed on the basic set (1) in the following manner: �(�, �, �) = (2�n)- 2 �- n �� ∑︁ � (�, �)�˜* (�, �). (3) � � � � � Then the rule for construction of quantum-mechanical operators with non-negative QDF may be formulated as follows: a classical function �(�, �, �) corresponds to a linear operator �(�), which acts on a generic Fourier-transformable function �(�, �) in the following way: ∫︁ � ′ �(�)�(�, �) = (2�n)-� �(� - �, � - �)�n (�-� )��(�, �, �)�(�′, �) d�d�d�d�′. (4) It follows that the formulae for calculation of means are in fact equivalent: ∫︁ and where the expression ⟨�⟩ = ∫︁ ⟨�⟩� = �*(�, �)�(�)�(�, �) d� (5) �(�, �, �)��(�, �, �) d� d�, (6) ⃒∫︁ � ⃒2 �� (�, �, �) = (2�n)-� ∑︁ ⃒ �*(� - �, �)�(�, �)�- n (�,�)⃒ (7) ⃒ � ⃒ ⃒ ⃒ defines the relation between wave functions and phase space probability densities. In this work we use the so called sturmian functions as auxiliary functions 1. The sturmian functions satisfy the equation [29] (︂ d2 -d�2 + �(� + 1) )︂ �2 + ���� (�) - �0 ���(�) = 0, (8) with the Coulomb potential � (�) = -2/� and with boundary conditions ���(0) = 0, ���(∞) = 0. Here � = 1, 2, 3, . . .; � = 0, 1, . . . , � - 1; �0 < 0; ��� = ��; � = √-�0. The solution of the equation (8) may be expressed in terms of the associated Laguerre polynomials: �+� ���(�) = ����-�� (2��)�+1�2�+1(2��), where � = √-�0 and the normalization constant are: � 2 1 {︂(� - � - 1)!}︂ . �� = 2 [(� + �)!]3 In our package SourceFunction we use sturmian functions devided by � and normalized to 1. In the case of � = 1, � = 0, sturmian and usual Coulomb wave functions coincide exactly, but as � grows, the amplitude and number of harmonics of sturmian functions differ from the Coulomb ones. In the course of calculations some formulae are very useful, such as products of spherical harmonics √︃ ∑︁ (2�1 + 1)(2�2 + 1) �0 �� ��1 �1 (�, �)��2 �2 (�, �) = �,� 4�(2� + 1) ��10 �20��1 �1 �2 �2 ��� (�, �), (9) as well as the formula for the integral of the triple product of spherical functions 2� � ∫︁ ∫︁ d� 0 0 �3 �3 d���1 �1 (�, �)��2 �2 (�, �)� * (�, �) = √︃(2�1 + 1)(2�2 + 1) = 4�(2�3 + 1) � �30 �10 �20 � �3 �3 �1 �1 �2 �2 ��� (�, �) (10) and the expansion: 1 ∞ = 4� ∑︁ 1 �� - �+1 � �� ∑︁ ���(�� , �� )� * (��, ��), (11) |⃗� - ⃗�| � 2� + 1 �=0 + �=-� where �- = min(|⃗�|, |⃗�|), �+ = max(|⃗�|, |⃗�|). The main goal of the program is a numerical calculation of the spectrum of hydrogen and some alkali metals, using the Ritz method with 91 × 91 matrices. This corresponds to the use of the basis of functions with the values of the principal quantum number � = 1, 2, 3, 4, 5, 6. The obtained values are then compared with experimental data to find the optimal parameters that provide the best fit for the data. In addition, an algorithm was developed for the numerical calculation of transition probabilities for hydrogen-like atoms and a procedure for comparing these data with those obtained experimentally. According to the second theorem of Kato, if � (�) → 0 for |⃗�| → ∞, then the essential spectrum of the operator �ˆ coincides with the essential spectrum of the operator �ˆ0, namely, it is half-line �+ = [0, ∞). The essential spectrum of an operator consists of all points of the spectrum, with the exception of isolated eigenvalues of finite multiplicity. Thus, the spectrum consists of two disjoint sets: eigenvalues of finite multiplicity and the essential spectrum. The Hamilton operator of a hydrogen-like atom in a quantum mechanics with nonnegative distribution function is �ˆ = �ˆ0 + ��ˆ + � (⃗�). The function � (⃗�) satisfies the conditions of the second Kato theorem. Therefore, the discrete spectrum of finite multiplicity lies below the essential spectrum and belongs to the interval [�- 1, �). Consequently, for the Hamilton operator �(�), we apply the minimax Rayleigh-Ritz method, and the eigenvalues �� for �-dimensional Ritz matrices converge to the eigenvalues �∞ = �� � � of this operator. The Hamiltonian operator �(�) of a hydrogen-like atom in quantum mechanics with nonnegative quantum distribution function is defined by the relations �{�� }(�) = �ˆ + �{�� }�ˆ + �{�� }(⃗�) for a set of valid functions {�� ∈ �2(�3)}. The operator �(�) is self-adjoint and bounded below by a constant (� - 1). The constant � depends on the choice of auxiliary functions. In view of the nonorthogonality of the used system of coordinate functions, to find the spectrum, it is necessary to solve a generalized eigenvalue problem � ⃗� = ��⃗�, (12) where � is the Ritz matrix �(�), and � is the matrix of scalar products of coordinate functions. The elements of the Ritz matrix are calculated by the formula: � � ∫︁ �0 [︂ (︂ �⃗ 2 )︂ (︂ �eff �2 )︂]︂ � �� = �� (⃗�) �� 2� + �� - � �� 0 d⃗�. (13) The result is a matrix that depends on the parameters �0, �� , �� and �eff. These parameters were used to find the spectrum that best fits the experimental data. The calculations were carried out using the Ritz matrices of dimensions 55 and 91, which correspond to the first five or six principal quantum numbers �. Nevertheless, even for such limited dimensions, it was possible to obtain a spectrum of hydrogen and some alkali metals with high accuracy. To assess the quality of the model, we used the function of the quadratic discrepancy �Ω = ∑︁ (�� - �exp)2 / (�exp)2 . (14) � � �∈Ω As shown in our works [22,24,25,27,28], for a hydrogen atom, the first 30 spectral values may be obtained with high accuracy. Moreover, the only sensitive parameter was the value �0. For parameters �� in the case of a hydrogen atom, the optimal value regardless of the number of spectral values in the residual is the value �� = 0. As to the coefficients �� , only �1 = 1, and the remaining coefficients in the optimization turned to zero. In contrast to the hydrogen atom, in the case of alkali metals it is impossible to build a single effective potential capable of describing the energy levels for different quantum numbers �, �. This occurs due to the interaction of the external electron with internal electrons, and the nature of this interaction is fundamentally different for different values of the orbital number �. To verify the quality of the model and the correctness of the choice of auxiliary functions, transition probabilities may be calculated. In the described program, the transition probabilities are calculated by the Galerkin method with the Sturm functions of the hydrogen atom as coordinate functions, which allows us to reduce the calculations to algebraic operations with matrix elements that are calculated analytically. 3. The Program The program is written up in the language MAPLE according to the algorithm described above. Figure 1 presents the general scheme of the program. It consists of nine stages. The first stage represents the original package for calculations of necessary functions. The second stage is devoted to calculations of potential energy operators. The third stage is devoted to calculations of kinetic energy operators. The fourth stage is devoted to calculations of Ritz matrices. The fifth stage is devoted to generation of codes for Ritz matrices for later use in other programs. The sixth stage is devoted to calculation of spectra of the Hamilton operator for an electron in hydrogen-like atoms. The seventh stage is devoted to optimization of parameters of the model for hydrogen. The final stage is devoted to optimization of parameters of the model for alkaline metals. Finally, the ninth part is the calculation of transition probabilities for model verification. The structures of the parts of the software package are shown in more detail: for calculating auxiliary functions (Figure 2), for calculating matrices and energy levels (Figure 3), and for calculating transition probabilities (Figure 4). Figure 1. The global structure of the program, QDF01 Figure 2. The structure of the part of the program for calculating auxiliary functions Figure 3. The structure of the part of the program for calculating matrices and energy levels Figure 4. The structure of the part of the program for calculating transition probabilities The subpackage “SourceFunctions” consists of the following parts. Calculation of spherical harmonics according to well-known formulae. The Legendre polynomials are generated using the Rodrigues formula. Calculations of sturmian functions, further these functions are normalized. Then follow calculations of auxiliary QDF functions (as products of sturmian functions and spherical harmonics) and calculations of fourier transforms of auxiliary functions. Let us note that all integrals are calculated analitically. Then follow calculation of Clebsh-Gordan ��(3) coefficients according to common wellknown formulae. We do not intend to carry out extensive calculations and investigations in this field, since there exist excellent Racah and Bethe packages and their extensions [30-33]. Calculations of products of spherical harmonics according to (9) etc. Matrix elements of Ritz matrices (13) are computed and written to external files. According to the Ritz method, eigenvalues of Ritz matrices represent spectral values of the quantity under investigation, i.e. energy. As stated above, the algorithm consists in solving the generalized eigenvalues problem � ⃗� = �� ⃗� where � is a Ritz matrix and � is the matrix of pair-to-pair scalar products of coordinate functions. The program allows calculating of Ritz matrices of arbitrary dimension. However, there exist apparatus restrictions due to the degree of the dimension. Currently, calculations with dimensions 55 and 91 are available. The parameter �0 is adjustable. For example, Figure 5 shows the dependence of the residual on �0 for the lithium atom and � + � = 5. Figure 5. Graph of the quadratic discrepancy on the parameter �0 in the case of the group of atomic levels of the lithium atom with � + � = 5. The minimum is equal to 0.088 and corresponds to �0 = -0.01414. Other options are �1 = 3.33; �2 = 0.001 = �3 = �4 = 0.001; �eff = 1; �� = 1 Residuals for smaller numbers of levels have similar deep minimax. But for alkali metal atoms, unlike a hydrogen atom, it is impossible to construct a unique effective potential capable of describing energy levels for different quantum numbers �, � [27]. One of the ways for describing energy levels in atoms of alkali metals is the so-called quantum defect model [16]. In this model, in a hydrogen-like spectrum, either the effective principal quantum number or, which is equivalent, the effective nucleus charge �eff (which differs from the true charge because of the screening by the inner electrons) is used. Figure 6 shows the dependence of the residual on the effective charge �eff for the lithium atom and � + � = 5. It is also possible to vary the effective electron mass. Figure 7 shows the dependence of the residual on the effective mass �� for the lithium atom and �+ � = 5. Figure 6. Graph of the quadratic discrepancy on the parameter �eff in the case of the group of atomic levels of the lithium atom with � + � = 5. The minimum is equal to 0.085 and corresponds to �eff = -1.242. Other options are �1 = 3.33; �2 = 0.001 = �3 = �4 = 0.001; �0 = -0.1414; �� = 1 Figure 7. Graph of the quadratic discrepancy on the parameter �� in the case of the group of atomic levels of the lithium atom with � + � = 5. The minimum is equal to 0.085 and corresponds to �� = 0.99. Other options are �1 = 3.33; �2 = 0.001 = �3 = �4 = 0.001; �0 = -0.1414; �eff = 1.242 4. Conclusion Quantum mechanics with a non-negative distribution function ceases to be an exotic theory and becomes a real design scheme that has its advantages. This happened after formulating it as a theory of quantum measurements. The used fitting parameters have a clear physical meaning and allow one to take into account and describe real phenomena, for example, the screening of the nuclear field by internal electrons and the interaction of an external electron with them. The proposed program allows one to effectively carry out such calculations, and also permits generalization to other quantum systems.

### Alexander V Zorin

Peoples’ Friendship University of Russia (RUDN University)
**Author for correspondence.**

Email: zorin@mx.rudn.ru

6 Miklukho-Maklaya str., Moscow, 117198, Russian Federation

Associate Professor, Candidate of Sciences in Physics and Mathematics, Associate Professor of Peoples’ Friendship University of Russia (RUDN University)

### Nikolay P Tretyakov

The Russian Presidential Academy of National Economy and Public Administration; Russian State Social University
Email: trn11@rambler.ru

Prospect Vernadskogo, 82, Moscow 119571, Russian Federation; Wilhelm Pieck str., 4-1, Moscow, 119571, Russian Federation

Associate Professor, Candidate of Sciences in Physics and Mathematics, Associate Professor of Department of Mathematics and Information Technologies of Russian State Social University

- W. Eissner, H. Nussbaumer, Resonances in Cross Sections for Excitation of Forbidden Lines in O2+, Journal of Physics B: Atomic, Molecular and Optical Physics 2 (3) (1969) 1028-1043. doi: 10.1088/0022-3700/2/3/305.
- J. C. Slater, Atomic Shielding Constants, Physical Review 36 (1930) 57-64. doi: 10.1103/PhysRev.36.57.
- G. H. S. E. U. Condon, The Theory of Atomic Spectra, Cambridge University Press, Cambridge, 1970.
- D. Layzer, On a Screening Theory of Atomic Spectra, Annals of Physics 8 (1959) 271-296. doi: 10.1016/0003-4916(59)90023-5.
- H. Nussbaumer, Improved Bound Wave Functions for Complex Atoms, Journal of Physics B: Atomic, Molecular and Optical Physics 5 (10) (1972) 1837-1843. doi: 10.1088/0022-3700/5/10/012.
- P. G. Burke, A. Hibbert, W. D. Robb, Wavefunctions and Oscillator Strengths of the Beryllium Iso-Electronic Sequence, Journal of Physics B: Atomic, Molecular and Optical Physics 5 (1) (1972) 37-43. doi: 10.1088/0022-3700/5/1/013.
- R. N. Zare, Correlation Effects in Complex Spectra. II. Transition Probabilities for the Magnesium Isoelectronic Sequence, Journal of Chemical Physics 47 (1967) 3561-72. doi: 10.1063/1.1712423.
- A. W. Weiss, Theoretical multiplet strengths for Mg I, Al II, and Si III, Journal of Chemical Physics 47 (1967) 3573-3578. doi: 10.1063/1.1712424.
- H. Friedrich, E. Trefftz, Configuration Mixing and Oscillator Strengths for Some Two-Electron Spectra (Ca I, Ba I, and Others), Journal of Quantitative Spectroscopy and Radiative Transfer 9 (1969) 333-359. doi: 10.1016/0022-4073(69)90030-2.
- Y.-K. Kim, P. S. Bagus, Oscillator Strengths for the Resonance Transitions in Alkaline Earth Atoms, Journal of Physics B: Atomic, Molecular and Optical Physics 5 (10) (1972) L193-L195. doi: 10.1088/0022-3700/5/10/001.
- S. Hameed, A. Herzenberg, M. G. James, Core Polarization Corrections to Oscillator Strengths in the Alkali Atoms, Journal of Physics B: Atomic, Molecular and Optical Physics 1 (5) (1968) 822-830. doi: 10.1088/0022-3700/1/5/308.
- S. Hameed, Core Polarization Corrections to Oscillator Strengths and Singlet-Triplet Splittings in Alkaline Earth Atoms, Journal of Physics B: Atomic, Molecular and Optical Physics 9 (4) (1972) 746-760. doi: 10.1088/0022-3700/5/4/009.
- M. Jones, Relativistic Corrections to Atomic Energy Levels, Journal of Physics B: Atomic, Molecular and Optical Physics 3 (12) (1972) 1571-1592. doi: 10.1088/00223700/3/12/003.
- M. Jones, Mutual Spin-Orbit and Spin-Spin Interactions in Atomic Structure Calculations, Journal of Physics B: Atomic, Molecular and Optical Physics 4 (11) (1971) 1422-1439. doi: 10.1088/0022-3700/4/11/006.
- W. Eissner, M. Jones, H. Nussbaumer, Techniques for the Calculation of Atomic Structures and Radiative Data Including Relativistic Corrections, Computer Physics Communications 8 (1974) 270-306. doi: 10.1016/0010-4655(74)90019-8.
- S. E. Frisch, Optical Atomic Spectrs, Fizmatgiz, Moscow, 1963, in Russian.
- M. G. Veselov, L. N. Labzovski, Theory of Atom. Structure of Electron Scells, Nauka, Moscow, 1986, in Russian.
- I. I. Sobelman, Atomic Spectra and Radiative Transitions, 2nd Edition, Springer, Berlin, 1996.
- A. V. Zorin, A. L. Sevastianov, L. A. Sevastianov, Application of the Noncommutative Theory of Statistical Decisions to the Modeling of Quantum Communication Channels, IEEE, 2017, pp. 26-31. doi: 10.1109/ICUMT.2017.8255195.
- V. V. Kuryshkin, La mechanique quantique avec une function nonnegative de distribution dans l’espace des phases, Annales de l’I.H.P. Physique th´eorique 17 (1972) 81-95.
- V. V. Kuryshkin, Some Problems of Quantum Mechanics Possessing a Non-Negative Phase-Space Distribution Function, International Journal of Theoretical Physics 7 (1973) 451-466. doi: 10.1007/BF00713247.
- A. V. Gorbachev, L. A. Sevastianov, A. V. Zorin, Kuryshkin-Wodkiewicz Model of Quantum Measurements for Atoms and Ions with One Valence Electron, Bulletin of Peoples’ Friendship University of Russia. Series: Mathematics. Information Sciences. Physics 2 (2016) 44-52.
- A. V. Zorin, V. V. Kuryshkin, L. A. Sevastianov, Description of the Spectrum of a Hydrogen-Like Atom, Bulletin of Peoples’ Friendship University of Russia. Series: Physics 6 (1998) 62-66, in Russian.
- A. V. Zorin, L. A. Sevastianov, Hydrogen-Like Atom with Nonnegative Quantum Distribution Function, Physics of Atomic Nuclei 70 (2007) 792-799. doi: 10.1134/S1063778807040229.
- L. Sevastianov, A. Zorin, A. Gorbachev, Pseudo-Differential Operators in the Operational Model of a Quantum Measurement of Observables, Lecture Notes in Computer Science 7152 (2012) 174-181. doi: 10.1007/978-3-642-28212-6. 17.
- A. V. Zorin, The Operational Model of Quantum Measurement of Kuryshkin- Wodkiewicz, Bulletin of Peoples’ Friendship University of Russia. Series: Mathematics. Information Sciences. Physics 2 (2012) 42-54, in Russian.
- A. V. Zorin, L. A. Sevastianov, N. P. Tretyakov, Computer Modelling of Hydrogen-Like Atoms in Quantum Mechanics with Nonnegative Distribution Function, Programming and Computer Software 33 (2007) 94-104. doi: 10.1134/S0361768807020077.
- L. A. Sevastianov, A. V. Zorin, The Computer-Based Model of Quantum Measurement, Physics of Atomic Nuclei 80 (2017) 774-780. doi: 10.1134/S1063778817040238.
- M. Rotenberg, Theory and Application of Sturmian Functions, Advances in Atomic and Molecular Physics 6 (1970) 233-268. doi: 10.1016/S0065-2199(08)60206-7.
- K. Rykhlinskaya, S. Fritzsche, Use of Group Theory for the Analysis of Vibrational Spectra, Computer Physics Communications 162 (2004) 124-142. doi: 10.1016/j.cpc.2004.06.088.
- K. Rykhlinskaya, S. Fritzsche, Generation of Molecular Symmetry Orbitals for the Point and Double Groups, Computer Physics Communications 171 (2005) 119-132. doi: 10.1016/j.cpc.2005.03.112.
- K. Rykhlinskaya, S. Fritzsche, Generation of Clebsch-Gordan Coefficients for the Point and Double Groups, Computer Physics Communications 174 (2006) 903-913. doi: 10.1016/j.cpc.2006.01.001.
- O. Gaigalas, O. Scharf, S. Fritzsche, Maple Procedures for the Coupling of Angular Momenta. VIII. Spin-Angular Coefficients for Single-Shell Configurations, Computer Physics Communications 166 (2005) 141-169. doi: 10.1016/j.cpc.2004.11.003.