Nuclei-selected NMR shielding calculations: a sublinear-scaling quantum-chemical method.

An ab initio method for the direct calculation of NMR shieldings for selected nuclei at the Hartree-Fock and density-functional theory level is presented. Our method shows a computational effort scaling only sublinearly with molecular size, as it is motivated by the physical consideration that the chemical shielding is dominated by its local environment. The key feature of our method is to avoid the conventionally performed calculation of all NMR shieldings but instead to solve directly for specific nuclear shieldings. This has important implications not only for the study of large molecules, but also for the simulation of solvent effects and molecular dynamics, since often just a few shieldings are of interest. Our theory relies on two major aspects both necessary to provide a sublinear scaling behavior: First, an alternative expression for the shielding tensor is derived, which involves the response density matrix with respect to the nuclear magnetic moment instead of the response to the external magnetic field. Second, as unphysical long-range contributions occur within the description of distributed gauge origin methods that do not influence the final expectation value, we present a screening procedure to truncate the B-field dependent basis set, which is crucial in order to ensure an early onset of the sublinear scaling. The screening is in line with the r(-2) distance decay of Biot-Savarts law for induced magnetic fields. Our present truncation relies on the introduced concept of "individual gauge shielding contributions" applied to a reformulated shielding tensor, the latter consisting of gauge-invariant terms. The presented method is generally applicable and shows typical speed-ups of about one order of magnitude; moreover, due to the reduced scaling behavior of O(1) as compared to O(N), the wins become larger with increasing system size. We illustrate the validity of our method for several test systems, including ring-current dominated systems and biomolecules with more than 1000 atoms.

[1]  Johannes Grotendorst,et al.  Modern methods and algorithms of quantum chemistry , 2000 .

[2]  R. Ditchfield,et al.  Self-consistent perturbation theory of diamagnetism , 1974 .

[3]  F. Jensen.  Introduction to Computational Chemistry , 1998 .

[4]  J. Autschbach,et al.  Calculation of optical rotation with time-periodic magnetic-field-dependent basis functions in approximate time-dependent density-functional theory. , 2005, The Journal of chemical physics.

[5]  R. S. Mulliken.  Criteria for the Construction of Good Self‐Consistent‐Field Molecular Orbital Wave Functions, and the Significance of LCAO‐MO Population Analysis , 1962 .

[6]  Aage E. Hansen,et al.  Localized orbital/local origin method for calculation and analysis of NMR shieldings. Applications to 13C shielding tensors , 1985 .

[7]  Jörg Kussmann,et al.  Linear-scaling method for calculating nuclear magnetic resonance chemical shifts using gauge-including atomic orbitals within Hartree-Fock and density-functional theory. , 2007, The Journal of chemical physics.

[8]  The helix-hinge-helix structural motif in human apolipoprotein A-I determined by NMR spectroscopy. , 1997, Biochemistry.

[9]  Christian Ochsenfeld,et al.  A reformulation of the coupled perturbed self-consistent field equations entirely within a local atomic orbital density matrix-based scheme , 1997 .

[10]  Christian Ochsenfeld,et al.  Efficient linear-scaling calculation of response properties: density matrix-based Laplace-transformed coupled-perturbed self-consistent field theory. , 2008, The Journal of chemical physics.

[11]  Christian Ochsenfeld,et al.  Linear scaling exchange gradients for Hartree–Fock and hybrid density functional theory , 2000 .

[12]  J. Kussmann,et al.  Linear‐Scaling Methods in Quantum Chemistry , 2007 .

[13]  L. Landau,et al.  Lehrbuch der theoretischen Physik , 2007 .

[14]  Benny G. Johnson,et al.  THE CONTINUOUS FAST MULTIPOLE METHOD , 1994 .

[15]  Christian Ochsenfeld,et al.  Linear and sublinear scaling formation of Hartree-Fock-type exchange matrices , 1998 .

[16]  Dennis R. Salahub,et al.  Calculations of NMR shielding constants by uncoupled density functional theory , 1993 .

[17]  W. Lamb.  Internal Diamagnetic Fields , 1941 .

[18]  Trygve Helgaker,et al.  Nuclear shielding constants by density functional theory with gauge including atomic orbitals , 2000 .

[19]  Paul Adrien Maurice Dirac,et al.  Note on the Interpretation of the Density Matrix in the Many-Electron Problem , 1931, Mathematical Proceedings of the Cambridge Philosophical Society.

[20]  Jörg Kussmann,et al.  A density matrix-based method for the linear-scaling calculation of dynamic second- and third-order properties at the Hartree-Fock and Kohn-Sham density functional theory levels. , 2007, The Journal of chemical physics.

[21]  J. Gauss.  Calculation of NMR chemical shifts at second-order many-body perturbation theory using gauge-including atomic orbitals , 1992 .

[22]  Arash A. Mostofi,et al.  A converse approach to the calculation of NMR shielding tensors , 2007, 0709.4429.

[23]  H. Shinbori,et al.  Calculation of NMR chemical shifts. VI. Gauge invariant and Hermitian condition , 1985 .

[24]  Jörg Kussmann,et al.  Ab initio NMR spectra for molecular systems with a thousand and more atoms: a linear-scaling method. , 2004, Angewandte Chemie.

[25]  J. V. Vleck,et al.  The theory of electric and magnetic susceptibilities , 1934, The Mathematical Gazette.

[26]  J. F. Hinton,et al.  Nuclear magnetic resonance spectroscopy applied to Li+ complexation by small ligands: An ab initio and experimental study of 7Li chemical shifts , 1980 .

[27]  Werner Kutzelnigg,et al.  Theory of Magnetic Susceptibilities and NMR Chemical Shifts in Terms of Localized Quantities , 1982 .

[28]  J. Gauss,et al.  Electron-correlated approaches for the calculation of nmr chemical shifts , 2003 .

[29]  Matt Challacombe,et al.  A simplified density matrix minimization for linear scaling self-consistent field theory , 1999 .

[30]  E. Purcell,et al.  Resonance Absorption by Nuclear Magnetic Moments in a Solid , 1946 .

[31]  A. Dalgarno.  Perturbation theory for atomic systems , 1959, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[32]  J. Pople,et al.  High-resolution nuclear magnetic resonance , 1959 .

[33]  Peter Pulay,et al.  Efficient implementation of the gauge-independent atomic orbital method for NMR chemical shift calculations , 1990 .

[34]  Ian M. Mills,et al.  Force Constants and Dipole-Moment Derivatives of Molecules from Perturbed Hartree-Fock Calculations. I , 1968 .

[35]  Nicholas C. Handy,et al.  The density functional calculation of nuclear shielding constants using London atomic orbitals , 1995 .

[36]  Henry F. Schaefer,et al.  Analytic second derivative techniques for self-consistent-field wave functions. A new approach to the solution of the coupled perturbed hartree-fock equations , 1983 .

[37]  Hans W. Horn,et al.  Direct computation of second-order SCF properties of large molecules on workstation computers with an application to large carbon clusters , 1992 .

[38]  Henry F. Schaefer,et al.  On the evaluation of analytic energy derivatives for correlated wave functions , 1984 .

[39]  J. Pople.  The theory of chemical shifts in nuclear magnetic resonance I. Induced current densities , 1957, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[40]  Norman Ramsey.  Magnetic Shielding of Nuclei in Molecules , 1950 .

[41]  W. Kutzelnigg,et al.  Theory of magnetic susceptibilities and NMR chemical shifts in terms of localized quantities. II. Application to some simple molecules , 1982 .

[42]  F. London,et al.  Théorie quantique des courants interatomiques dans les combinaisons aromatiques , 1937 .

[43]  Christian Ochsenfeld,et al.  An atomic orbital-based reformulation of energy gradients in second-order Møller-Plesset perturbation theory. , 2008, The Journal of chemical physics.

[44]  H. F. Hameka.  On the nuclear magnetic shielding in the hydrogen molecule , 1958 .