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 (Formula presented.) odd and (Formula presented.) 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 (Formula presented.) and (Formula presented.) 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 mu, SCOPUS: ar.j, info:eu-repo/semantics/published

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 (Formula presented.) and (Formula presented.) states in (Formula presented.) Li 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., SCOPUS: ar.j, info:eu-repo/semantics/published

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., SCOPUS: ar.j, info:eu-repo/semantics/published