Computation of Atomic Properties with the Hy-CI Method

## Summary:

Impressive advances have been made throughout the years in the study of atomic structure, at both the experimental and theoretical levels. For atomic hydrogen and other equivalent two-body systems, exact analytical solutions to the nonrelativistic Schroedinger equation are known. It is now possible to calculate essentially exact nonrelativistic energies for helium (He) and other three-body (two-electron) systems as well. Even for properties other than the nonrelativistic energy, the precision of the calculation has been referred to as *essentially exact for all practical purposes*, i.e., the precision goes well beyond what can be achieved experimentally. These high-precision results for two-electron systems have been produced using wave functions which include interelectronic coordinates, a trademark of the classic Hylleraas (Hy) calculations in the 1920s. The Hy-CI method is a hybrid method which includes interelectronic coordinates in the wave function to mimic the high precision of Hy methods, but also includes configurational terms that are the trademark of the conventional Configuration Interaction (CI) methods employed in calculating energies for many-electron atomic (and molecular) systems.

## Additional Technical Details:

**Why is the Hy-CI Method Important?**

In any attempt to get very precise energies, large basis sets have to be employed, which means that linear dependence in the basis set is never very far away. To proceed to tens of thousands of terms in a wave function, extended precision arithmetic is needed to obviate the linear dependence problem, which in turn leads to higher CPU costs. The use of tens of thousands of terms in a wave function leads to memory problems arising from storage of the matrix elements prior to the matrix diagonalization step. And this is true already for the two-electron case of helium (He) and its isoelectronic ions (He-like systems). Where three electron atomic systems (lithium (Li) and other members of its isoelectronic series) have been treated essentially as accurately as He-like systems, demand on computer resources has increased by 6000 fold. Because of these computational difficulties, until recently in the four electron case (beryllium (Be) and other members of its isoelectronic series) there were no calculations of the ground or excited states with an error of less than 0.1 microhartrees (0.0000001 hartree). The challenge for computational scientists is to extend the phenomenal He accomplishments to three, four, and more electron atomic systems. This is where the Hy-CI method becomes important, because the use of configurations wherever possible leads to less difficult integrals than in a purely Hy method, and if one restricts the wave function to at most a single interelectronic coordinate to the first power, then the most difficult integrals are already dealt with at the four electron level and the calculation retains the precision of Hy techniques, but is greatly simplified.

## Applications:

The Hy-CI method has been used to compute not only energy levels, but also other atomic properties such as ionization potentials, electron affinities, electric polarizabilities, and transition probabilities of two, three, and four electron atoms and other members of their isoelectronic sequences. In the early parallel work, Hy-CI method calculations with up to 4648 expansion terms were carried out for the ground ^{1}S state of neutral helium and 4 of its isoelectronic ions, H-, Li+, Be++, and B+++. This has resulted in high-precision results for the nonrelativistic energies that are believed to be accurate to 20 decimal digits. This work employs a very novel wave function, namely, one consisting of at most a single r12 raised to the first power combined with a conventional non-orthogonal configuration interaction (CI) basis. This technique can be extended to multielectron systems (more than 3 or 4 electrons). The combination of computational simplicity of this form of the wave function, compared to other wave functions of comparable accuracy, as well as the use of parallel processing and extended precision arithmetic, make it possible to achieve levels of accuracy comparable to what has been achieved for He (a 2 electron atom), for atoms with more than 2 electrons. Indeed, the technique has been used to obtain results for the nonrelativistic energies for four excited states of the lithium atom which represented at the time the highest level of accuracy ever reached (10^{-10} hartree) in atomic quantum computations with more than two electrons. Recent work has included demonstrating what precision can be obtained for atomic beryllium, the key to multielectron (more than 4 electrons) systems, since the integrals that arise for more than 4 electrons are of the same type as the ones that arise in the 4 electron systems. One of our publications discusses the most difficult integral arising in Hy-CI calculations, the three-electron triangle integral. We found that convergence acceleration techniques could be used to speed up the convergence of the integral and hence remove the real bottleneck to highly accurate Hy-CI calculations. Another of our publications tackles making the four electron integrals, the real bottleneck at the 4 and more electron level, as efficient as possible. The utility of these techniques has been demonstrated in calculations of the nonrelativistic energies of the four electron atomic ground state with nuclear charge ranging from 4 through 113 (110 states) with record levels of accuracy. This work has demonstrated that doing really accurate calculations on atoms with N_{e} (number of electrons) greater than or equal to 5 is a real possibility.

## Publications:

- James S. Sims and Stanley A. Hagstrom, Mathematical and Computational Science Issues in high precision Hylleraas-Configuration Interaction variational calculations: III. Four-electron integrals,
*J. Phys. B: At. Mol. Opt. Phys*. 48, 175003(2015). - James S. Sims and Stanley A. Hagstrom, Hylleraas-configuration interaction nonrelativistic energies for the
^{1}S ground states of the beryllium isoelectronic sequence,*Journal of Chemical Physics*, DOI 10.1063/1.4881639, June 11, 2014. - James S. Sims and Stanley A. Hagstrom, Hylleraas-configuration interaction study of the
^{1}S ground state of neutral beryllium,*Physical Review A*, 83, 2011. ID: 032518.**Note:**DOI: 10.1103/PhysRevA.83.032518 - James S. Sims and Stanley A. Hagstrom, Hy-CI Study of the 2
^{2}S Ground State of Neutral Lithium and the First Five Excited^{2}S States, Physical Review A, 80, 2009. ID: 052507. ( DOI: 10.1103/PhysRevA.80.052507 ) - James S. Sims and Stanley A. Hagstrom, Math and computational science issues in high-precision Hy-CI calculations II. Kinetic Energy and electron-nucleus int eraction integrals, Journal of Physics B: Atomic, Molecular, and Optical Physics, 40, 2007, pp. 1575-1587.
- James S. Sims and Stanley A. Hagstrom, James S. Sims and Stanley A. Hagstrom, High Precision Variational Born-Oppenheimer Energies of the Ground State of the Hydrogen Molecule , Journal of Chemical Physics, 124 (9) , 3/7/2006, p. 7. (09410)
- James S. Sims and Stanley A. Hagstrom, Math and computational science issues in high-precision Hy-CI calculations I. Three-electron integrals, Journal of Physics B: Atomic, Molecular, and Optical Physics, 37 (7) , 2004, pp. 1519-1540.
- James S. Sims and Stanley A. Hagstrom, Analytic Value of the Atomic Three-electron Integral with Slater Wave Functions, Physical Review A, 68, 2003, p. 016501. ( Comment: Phys. Rev. A 44,5492(1991) )
- James S. Sims and Stanley A. Hagstrom, High Precision Hy-CI Variational Calculations for the Ground State of Neutral Helium and Heliumlike Ions, International Journal of Quantum Chemistry, 90 (6) , 2002, pp. 1600-1609. Note: DOI 10.1002/qua.10344
- James S. Sims, William L. George, Steven G. Satterfield, Howard K. Hung, John G. Hagedorn, Peter M. Ketcham, Terence J. Griffin, Stanley A. Hagstrom, Julien C. Franiatte, Garnett W. Bryant, W. Jaskolski, Nicos S. Martys, Charles E. Bouldin, Vernon Simmons, Olivier P. Nicolas, James A. Warren, Barbara A. am Ende, John E. Koontz, B. James Filla, Vital G. Pourprix, Stefanie R. Copley, Robert B. Bohn, Adele P. Peskin, Yolanda M. Parker and Judith E. Devaney, Accelerating Scientific Discovery Through Computation and Visualization II, NIST Journal of Research, 107 (3) , May-June, 2002, pp. 223-245.