This paper elaborates the integral transformation technique and uses it for the case of the non‐relativistic kinetic and Coulomb potential energy operators, as well as for the relativistic mass‐velocity and Darwin terms. The techniques are tested for the ground electronic state of the helium atom and perturbative relativistic energies are reported for the ground electronic state of the H3+ molecular ion near its equilibrium structure.

A recently developed lower bound theory for Coulombic problems (E. Pollak, R. Martinazzo, J. Chem. Theory Comput. 2021Comput. , 17, 1535 is further developed and applied to the highly accurate calculation of the ground-state energy of two-(He, Li + , and H − ) and three-(Li) electron atoms. The method has been implemented with explicitly correlated many-particle basis sets of Gaussian type, on the basis of the highly accurate (Ritz) upper bounds they can provide with relatively small numbers of functions. The use of explicitly correlated Gaussians is developed further for computing the variances, and the necessary modifications are here discussed. The computed lower bounds are of submilli-Hartree (parts per million relative) precision and for Li represent the best lower bounds ever obtained.Although not yet as accurate as the corresponding (Ritz) upper bounds, the computed bounds are orders of magnitude tighter than those obtained with other lower bound methods, thereby demonstrating that the proposed method is viable for lower bound calculations in quantum chemistry applications. Among several aspects, the optimization of the wave function is shown to play a key role for both the optimal solution of the lower bound problem and the internal check of the theory.

Despite the fact that most quantum chemistry basis sets are designed for accurately modeling valence chemistry, these general-purpose basis sets continue to be widely used to model core-dependent properties. Core-specialized basis sets are designed with specific features to accurately represent the behavior of the core region. This design typically incorporates Gaussian primitives with higher exponents to capture core behavior effectively, as well as some decontraction of basis functions to provide flexibility in describing the core electronic wave function. The highest Gaussian exponent and the degree of contraction for both s- and p-basis functions effectively characterize these design aspects. In this study, we compare the design and performance of general-purpose basis sets against several literature-based basis sets specifically designed for three core-dependent properties: J coupling constants, hyperfine coupling constants, and magnetic shielding constants (used for calculating chemical shifts). Our findings consistently demonstrate a significant reduction in error when employing core-specialized basis sets, often at a marginal increase in computational cost compared to the popular 6-31G** basis set. Notably, for expedient calculations of J coupling, hyperfine coupling, and magnetic shielding constants, we recommend the use of the pcJ-1, EPR-II, and pcSseg-1 basis sets, respectively. For higher accuracy, the pcJ-2, EPR-III, and pcSseg-2 basis sets are recommended.