Search

Hyperfine structure constants have many applications, but are often hard to calculate accurately due to large and canceling contributions from different terms of the hyperfine interaction operator, and also from different closed and spherically symmetric core subshells that break up due to electron correlation effects. In multiconfiguration calculations, the wave functions are expanded in terms of configuration state functions (CSFs) built from sets of one-electron orbitals. The orbital sets are typically enlarged within the layer-by-layer approach. The calculations are energy-driven, and orbitals in each new layer of correlation orbitals are spatially localized in regions where the weighted total energy decreases the most, overlapping and breaking up different closed core subshells in an irregular pattern. As a result, hyperfine structure constants, computed as expectation values of the hyperfine operators, often show irregular or oscillating convergence patterns. Large orbital sets, and associated large CSF expansions, are needed to obtain converged values of the hyperfine structure constants. We analyze the situation for the states of the {2s22p3,2s22p23p,2s22p24p} odd and {2s22p23s,2s2p4,2s22p24s,2s22p23d} even configurations in N I, and show that the convergence with respect to the increasing sets of orbitals is radically improved by introducing separately optimized orbital sets targeted for describing the spin- and orbital-polarization effects of the 1s and 2s core subshells that are merged with, and orthogonalized against, the ordinary energy-optimized orbitals. In the layer-by-layer approach, the spectroscopic orbitals are kept frozen from the initial calculation and are not allowed to relax in response to the introduced layers of correlation orbitals. To compensate for this lack of variational freedom, the orbitals are transformed to natural orbitals prior to the final calculation based on single and double substitutions from an increased multireference set. The use of natural orbitals has an important impact on the states of the 2s22p23s configuration, bringing the corresponding hyperfine interaction constants in closer agreement with experiment. Relying on recent progress in methodology, the multiconfiguration calculations are based on configuration state function generators, cutting down the time for spin-angular integration by factors of up to 50, compared to ordinary calculations.

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.

The latest published version of GRASP (General-purpose Relativistic Atomic Structure Package), i.e., GRASP2018, retains a few suboptimal subroutines/algorithms, which reflect the limited memory and file storage of computers available in the 1980s. Here we show how the efficiency of the relativistic self-consistent-field (SCF) procedure of the multiconfiguration-Dirac–Hartree–Fock (MCDHF) method and the relativistic configuration-interaction (RCI) calculations can be improved significantly. Compared with the original GRASP codes, the present modified version reduces the CPU times by factors of a few tens or more. The MPI performances for all the original and modified codes are carefully analyzed. Except for diagonalization, all computational processes show good MPI scaling.

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.

In multiconfiguration Dirac–Hartree–Fock (MCDHF) calculations, there is a strong coupling between the localization of the orbital set and the configuration state function (CSF) expansion used to determine it. Furthermore, it is well known that an orbital set resulting from calculations, including CSFs describing core–core correlation and other effects, which aims to lower the weighted energies of a number of targeted states as much as possible, may be inadequate for building CSFs that account for correlation effects that are energetically unimportant but decisive for computed properties, e.g., hyperfine structures or transition rates. This inadequacy can be traced in irregular or oscillating convergence patterns of the computed properties as functions of the increasing orbital set. In order to alleviate the above problems, we propose a procedure in which the orbital set is obtained by merging several separately optimized, and mutually non-orthogonal, orbital sets. This computational strategy preserves the advantages of capturing electron correlation on the total energy through the variational MCDHF method and allows to target efficiently the correlation effects on the considered property. The orbital sets that are merged are successively orthogonalized against each other to retain orthonormality. The merged orbital set is used to build CSFs that efficiently lower the energy and also adequately account for the correlation effects that are important for the property. We apply the procedure to compute the hyperfine structure constants for the 1s22s2S1/2 and 1s22p2P1/2,3/2o states in 7Li and show that it leads to considerably improved convergence patterns with respect to the increasing orbital set compared to standard calculations based on a single orbital set, energy-optimized in the variational procedure. The perspectives of the new procedure are discussed in a broader context in the summary.

We modified the Hfs92 code of the GRASP package in order to describe the magnetic octupole hyperfine interaction. To illustrate the utility of the modified code, we carried out state-of-the-art calculations of the electronic factors of the magnetic octupole hyperfine interaction constants for levels in the ground configuration of the Bi atom. The nuclear magnetic octupole moment of the 209Bi isotope was extracted by combining old measurements of the hyperfine structures of 6p34S3/2o [Hull, R.; Brink, G. Phys. Rev. A 1970, 1, 685] and 2P3/2o [Landman, D.A.; Lurio, A. Phys. Rev. A 1970, 1, 1330] using the atomic-beam magnetic-resonance technique with our theoretical electronic factors. The present extracted octupole moment was consistent with all the available values but the one obtained in the single-particle nuclear shell model approximation. This observation supports the previous finding that nuclear many-body effects, such as the core polarization, significantly contribute to the nuclear magnetic octupole moment in the case of 209Bi.

In the case of the one-electron Dirac equation with a point nucleus, the virial theorem (VT) states that the ratio of the kinetic energy to potential energy is exactly −1, a ratio that can be an independent test of the accuracy of a computed solution. This paper studies the virial theorem for subshells of equivalent electrons and their interactions in many-electron atoms. This shows that the linear scaling of the dilation is achieved through the balancing of the contributions to the potential of an electron from inner and outer regions that some Slater integrals impose conditions on a single subshell, but others impose conditions between subshells. The latter slows the rate of convergence of the self-consistent field process in which radial functions are updated one at a time. Several cases are considered. Results are also extended to the nonrelativistic case.