Abstract:We developed and implemented a method-independent, fully numerical, finite difference approach to calculating nuclear magnetic resonance shieldings, using gauge-including atomic orbitals. The resulting capability can be used to explore non-standard methods, given only the energy as a function of finite-applied magnetic fields and nuclear spins. For example, standard second-order Møller-Plesset theory (MP2) has well-known efficacy for 1H and 13C shieldings and known limitations for other nuclei such as 15N and 17O. It is, therefore, interesting to seek methods that offer good accuracy for 15N and 17O shieldings without greatly increased compute costs, as well as exploring whether such methods can further improve 1H and 13C shieldings. Using a small molecule test set of 28 species, we assessed two alternatives: κ regularized MP2 (κ-MP2), which provides energy-dependent damping of large amplitudes, and MP2.X, which includes a variable fraction, X, of third-order correlation (MP3). The aug-cc-pVTZ basis was used, and coupled cluster with singles and doubles and perturbative triples [CCSD(T)] results were taken as reference values. Our κ-MP2 results reveal significant improvements over MP2 for 13C and 15N, with the optimal κ value being element-specific. κ-MP2 with κ = 2 offers a 30% rms error reduction over MP2. For 15N, κ-MP2 with κ = 1.1 provides a 90% error reduction vs MP2 and a 60% error reduction vs CCSD. On the other hand, MP2.X with a scaling factor of 0.6 outperformed CCSD for all heavy nuclei. These results can be understood as providing renormalization of doubles amplitudes to partially account for neglected triple and higher substitutions and offer promising opportunities for future applications.

Low-Scaling, Efficient and Memory Optimized Computation of Nuclear Magnetic Resonance Shieldings within the Random Phase Approximation Using Cholesky-Decomposed Densities and an Attenuated Coulomb Metric

An in-silico NMR laboratory for nuclear magnetic shieldings computed via finite fields: Exploring nucleus-specific renormalizations of MP2 and MP3

Sub-quadratic scaling real-space random-phase approximation correlation energy calculations for periodic systems with numerical atomic orbitals

Accurate NMR Shieldings with σ-Functionals

Stochastic Resolution-of-the-Identity Auxiliary-Field Quantum Monte Carlo: Scaling Reduction without Overhead

Application of the Adiabatic Connection Random Phase Approximation to Electron-Nucleus Hyperfine Coupling Constants

Static Subspace Approximation for Random Phase Approximation Correlation Energies: Implementation and Performance

Accurate Prediction of Nuclear Magnetic Resonance Shielding Constants: an Extension of the Focal-Point Analysis Method for Magnetic Parameter Calculations (FPA-M) with Improved Efficiency.

Efficient Implementation of the Random Phase Approximation with Domain-based Local Pair Natural Orbitals

Accurate Prediction of Nuclear Magnetic Resonance Shielding Constants: Towards the Accuracy of CCSD(T) Complete Basis Set Limit

Systematic studies on the computation of nuclear magnetic resonance shielding constants and chemical shifts: the density functional models.

LibRPA: A Software Package for Low-scaling First-principles Calculations of Random Phase Approximation Electron Correlation Energy Based on Numerical Atomic Orbitals

Revving up 13C NMR shielding predictions across chemical space: Benchmarks for atoms-in-molecules kernel machine learning with new data for 134 kilo molecules

Low-Scaling Algorithm for the Random Phase Approximation using Tensor Hypercontraction with k-point Sampling

Sparse tensor based nuclear gradients for periodic Hartree-Fock and low-scaling correlated wave function methods in the CP2K software package: A massively parallel and GPU accelerated implementation

Localized Resolution of Identity Approach to the Analytical Gradients of Random-Phase Approximation Ground-State Energy: Algorithm and Benchmarks

On the calculation of second-order magnetic properties using subsystem approaches in the relativistic framework

Massive-Parallel Implementation of the Resolution-of-Identity Coupled-Cluster Approaches in the Numeric Atom-Centered Orbital Framework for Molecular Systems

Calculation of Magnetic Shielding Constants with meta-GGA Functionals Employing the Multipole-Accelerated Resolution of the Identity: Implementation and Assessment of Accuracy and Efficiency

The importance of nuclear quantum effects for NMR crystallography

Efficient Random Phase Approximation for Diradicals