These are preliminary reports that have not been peer-reviewed. They should not be regarded as conclusive, guide clinical practice/health-related behavior, or be reported in news media as established information. For more information, please see our FAQs.
2 files

A Full Additive QM/MM Scheme for the Computation of Molecular Crystals with Extension to Many-Body Expansions

submitted on 12.11.2018, 10:24 and posted on 13.11.2018, 15:25 by Thorsten Teuteberg, Marco Eckhoff, Ricardo A. Mata
An additive quantum mechanics/molecular mechanics (QM/MM) model for the theoretical investigation of molecular crystals (AC-QM/MM) is presented. At the one-body level, a single molecule is chosen as the QM region. The MM region around it consists of a finite cluster of explicit MM atoms, represented by point charges and Lennard-Jones potentials, with additional background charges to mimic periodic electrostatics. Cluster charges are QM-derived and calculated self-consistently to ensure a polarizable embedding. We have also considered the extension to many-body QM corrections, calculating the interactions of a central molecule to neighbouring units in the crystal. Full gradient expressions have been derived, also including symmetry information. The scheme allows for the calculation of molecular properties as well as unconstrained optimisations of the molecular geometry and cell parameters with respect to the lattice energy. Benchmarking the approach with the X23 reference set confirms the convergence pattern of the many-body extension, although comparison to plane wave DFT reveals a systematic overestimation of cohesive energies by 6-16 kJ·mol−1 . While the scheme primarily aims to provide an inexpensive and flexible way to model a molecule in a crystal environment, it can also be used to reach highly accurate cohesive energies by the straightforward application of wave function correlated approaches. Calculations with local coupled cluster with singles, doubles, and perturbative triples, albeit limited to numerical gradients, show an impressive agreement with experimental estimates for small molecular crystals.


Email Address of Submitting Author


Georg-August-Universität Göttingen, Institut für Physikalische Chemie



ORCID For Submitting Author


Declaration of Conflict of Interest


Version Notes

The following article has been submitted to The Journal of Chemical Physics.