Quantum capping potentials with point charges: a simple QM/MM approach for the calculation of large-molecule NMR shielding tensors

A simple quantum-mechanics/molecular-mechanics (QM/MM) approach for calculating NMR shielding tensors (sigma) is presented. The method involves capping the QM region with quantum capping potentials (QCPs) and representing the MM region with point charges. Test calculations on simple systems without MM charges show that calculated sigma values improve relative to the full QM results with increasing distance between the capped bond and chromophore. Calculations on the histidine amino acid and cytosine monophosphate (CMP) nucleic acid show that the use of QCPs with point charges result in mean errors in the isotropic component of sigma that are less than 1.6 ppm. The results also reveal that, contrary to previous work, the explicit effect of point charges on sigma through coupling with gauge factors, as in the gauge including atomic orbital approach, is minimal for the CMP molecule. The present QM/MM approach for calculating sigma is easy to apply and requires no code modification.