Computational atomic physics continues to play a crucial role in both increasing the understanding of fundamental physics (e.g., quantum electrodynamics and correlation) and producing atomic data for interpreting observations from large-scale research facilities ranging from fusion reactors to high-power laser systems, space-based telescopes and isotope separators. A number of different computational methods, each with their own strengths and weaknesses, is available to meet these tasks. Here, we review the relativistic multiconfiguration method as it applies to the General Relativistic Atomic Structure Package [grasp2018, C. Froese Fischer, G. Gaigalas, P. Jönsson, J. Bieroń, Comput. Phys. Commun. (2018). DOI: 10.1016/j.cpc.2018.10.032]. To illustrate the capacity of the package, examples of calculations of relevance for nuclear physics and astrophysics are presented.

Grasp is a software package in Fortran 95, adapted to run in parallel under MPI, for research in atomic physics. The basic premise is that, given a wave function, any observed atomic property can be computed. Thus, the first step is always to determine a wave function. Different properties challenge the accuracy of the wave function in different ways. This software is distributed under the MIT Licence.

Atomic data, such as wavelengths, spectroscopic labels, broadening parameters and transition rates, are necessary for many applications, especially in plasma diagnostics, and for interpreting the spectra of distant astrophysical objects. The experiment with its limited resources is unlikely to ever be able to provide a complete dataset on any atomic system. Instead, the bulk of the data must be calculated. Based on fundamental principles and well-justified approximations, theoretical atomic physics derives and implements algorithms and computational procedures that yield the desired data. We review progress and recent developments in fully-relativistic multiconfiguration Dirac–Hartree–Fock methods and show how large-scale calculations can give transition energies of spectroscopic accuracy, i.e., with an accuracy comparable to the one obtained from observations, as well as transition rates with estimated uncertainties of a few percent for a broad range of ions. Finally, we discuss further developments and challenges.

The theoretical foundations of relativistic electronic structure theory within quantum electrodynamics (QED) and the computational basis of the atomic structure code GRASP are briefly surveyed. A class of four-component basis set is introduced, which we denote the CKG-spinor set, that enforces the charge-conjugation symmetry of the Dirac equation. This formalism has been implemented using the Gaussian function technology that is routinely used in computational quantum chemistry, including in our relativistic molecular structure code, BERTHA. We demonstrate that, unlike the kinetically matched two-component basis sets that are widely employed in relativistic quantum chemistry, the CKG-spinor basis is able to reproduce the well-known eigenvalue spectrum of point-nuclear hydrogenic systems to high accuracy for all atomic symmetry types. Calculations are reported of third- and higher-order vacuum polarization effects in hydrogenic systems using the CKG-spinor set. These results reveal that Gaussian basis set expansions are able to calculate accurately these QED effects without recourse to the apparatus of regularization and in agreement with existing methods. An approach to the evaluation of the electron self-energy is outlined that extends our earlier work using partial-wave expansions in QED. Combined with the treatment of vacuum polarization effects described in this article, these basis set methods suggest the development of a comprehensive ab initio approach to the calculation of radiative and QED effects in future versions of the GRASP code.

Large configuration interaction (CI) calculations can be performed if part of the interaction is treated perturbatively. To evaluate the combined CI and perturbative method, we compute excitation energies for the $3l3l^{\prime }$ , $3l4l^{\prime }$ and $3s5l$ states in Mg-like iron. Starting from a CI calculation including valence and core–valence correlation effects, it is found that the perturbative inclusion of core–core electron correlation halves the mean relative differences between calculated and observed excitation energies. The effect of the core–core electron correlation is largest for the more excited states. The final relative differences between calculated and observed excitation energies is 0.023%, which is small enough for the calculated energies to be of direct use in line identifications in astrophysical and laboratory spectra.