Abstract: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.

Tensorized orbitals for computational chemistry

Computationally Efficient Molecular Integrals of Solid Harmonic Gaussian Orbitals Using Quantum Entanglement of Angular Momentum

Virtual Orbital Many-Body Expansions: A Possible Route towards the Full Configuration Interaction Limit

Direct determination of optimal real-space orbitals for correlated electronic structure of molecules

Optimizing the Hartree-Fock Orbitals by the DMRG(2)(New Development of Numerical Simulations in Low-Dimensional Quantum Systems: from Density Matrix Renormalization Group to Tensor Network Formulations)

Even Faster Exact Exchange for Solids via Tensor Hypercontraction

Qubit Count Reduction by Orthogonally-Constrained Orbital Optimization for Variational Quantum Excited States Solvers

The influence of basis sets and ansatze building to quantum computing in chemistry

A quantum eigenvalue solver based on tensor networks

Efficient simulation of quantum chemistry problems in an enlarged basis set

Bi-orthonormal orbital optimization with a cheap core-electron free three-body correlation factor for Quantum Monte Carlo and Transcorrelation

Accurate and Interpretable Representation of Correlated Electronic Structure via Tensor Product Selected CI

Direct Unconstrained Optimization of Molecular Orbital Coefficients in Density Functional Theory

Direct energy functional minimization under orthogonality constraints

Natural Orbitals and Targeted Non-Orthogonal Orbital Sets for Atomic Hyperfine Structure Multiconfiguration Calculations

Quantum Solvers for Plane-Wave Hamiltonians: Abridging Virtual Spaces Through the Optimization of Pairwise Correlations

Robust linear-scaling optimization of compact localized orbitals in density functional theory

Full-Dimensional Schrödinger Wavefunction Calculations using Tensors and Quantum Computers: the Cartesian component-separated approach

The Coupled Cluster Singles, Doubles, and a Hybrid Treatment of Connected Triples Based on the Split Virtual Orbitals.

Self‐Consistent Molecular‐Orbital Methods. IX. An Extended Gaussian‐Type Basis for Molecular‐Orbital Studies of Organic Molecules