Automated fragmentation quantum mechanical calculation of C-13 and H-1 chemical shifts in molecular crystals

In this work, the automated fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) approach was applied to calculate the C-13 and H-1 nuclear magnetic resonance (NMR) chemical shifts in molecular crystals. Two benchmark sets of molecular crystals were selected to calculate the NMR chemical shifts. Systematic investigation was conducted to examine the convergence of AF-QM/MM calculations and the impact of various density functionals with different basis sets on the NMR chemical shift prediction. The result demonstrates that the calculated NMR chemical shifts are close to convergence when the distance threshold for the QM region is larger than 3.5 angstrom. For C-13 chemical shift calculations, the mPW1PW91 functional is the best density functional among the functionals chosen in this study (namely, B3LYP, B3PW91, M06-2X, M06-L, mPW1PW91, OB98, and OPBE), while the OB98 functional is more suitable for the H-1 NMR chemical shift prediction of molecular crystals. Moreover, with the B3LYP functional, at least a triple-zeta basis set should be utilized to accurately reproduce the experimental C-13 and H-1 chemical shifts. The employment of diffuse basis functions will further improve the accuracy for C-13 chemical shift calculations, but not for the H-1 chemical shift prediction. We further proposed a fragmentation scheme of dividing the central molecule into smaller fragments. By comparing with the results of the fragmentation scheme using the entire central molecule as the core region, the AF-QM/MM calculations with the fragmented central molecule can not only achieve accurate results but also reduce the computational cost. Therefore, the AF-QM/MM approach is capable of predicting the C-13 and H-1 NMR chemical shifts for molecular crystals accurately and effectively, and could be utilized for dealing with more complex periodic systems such as macromolecular polymers and biomacromolecules. The AF-QM/MM program for molecular crystals is available at https://github.com/shiman1995/NMR.

Automated summary

The AF-QM/MM approach was applied to accurately calculate the C-13 and H-1 NMR chemical shifts in molecular crystals, showing that the results are close to convergence under certain conditions. Different density functionals and basis sets have varying impacts on the predictions, with the choice affecting the accuracy of the results. By selecting suitable combinations, the accuracy of the predictions can be improved.