Track papers across labs, topics, and trends.
Search papers, labs, and topics across Lattice.
This paper introduces a scalable framework for computing nuclear magnetic resonance (NMR) shielding tensors using finite magnetic field differentiation of gauge-including atomic-orbital (GIAO) based self-consistent field (SCF) calculations. By leveraging hybrid MPI and OpenMP parallelization alongside a resolution-of-identity (RI) approximation, the method achieves significant efficiency improvements over traditional analytical implementations. The results show that this approach yields highly precise DFT-level NMR shielding constants, paving the way for future advancements in quantum chemistry methodologies.
Achieving DFT-level precision in NMR shielding calculations with a scalable, efficient framework could revolutionize how we interpret NMR spectra across diverse chemical systems.
Absolute nuclear magnetic shielding constants and relative chemical shifts are second-derivative response properties that underpin the interpretation of Nuclear Magnetic Resonance (NMR) spectra and provide critical insights into the structural and electronic environments of diverse chemical systems. However, their accurate computation is often constrained by the need for method-specific, analytical response implementation, and therefore is particularly challenging for non-variational, correlated wavefunction methods where analytical second derivatives are frequently unavailable. Here we present a scalable framework that computes NMR shielding tensors through finite magnetic field differentiation of complex-valued, gauge-including atomic-orbital (GIAO) based self-consistent field (SCF) calculations. By combining hybrid MPI and OpenMP parallelization and resolution-of-identity (RI) approximation-based Coulomb (J) and Exchange (K) implementation and mixed numerical/analytical derivatives computational scheme, we achieve very good efficiency for hybrid density functional theory (DFT) NMR shielding, relative to existing state-of-the-art analytic implementations across realistic chemical systems of various sizes. Forward differentiation with optimal step sizes derived from rigorous error analysis retains favorable numerical errors much smaller than experimental uncertainty or intrinsic DFT errors. The results demonstrate that RI-accelerated finite magnetic field calculations can obtain DFT-level NMR shielding constants very precisely, providing a scalable foundation for extensions to more advanced quantum chemistry methods in future work.