arXiv:2610.41523
A Hybrid Message-Passing and Electric-Field Representation for Electrostatic QM/MM Embedding with Machine Learning Interatomic Potentials
Raphael T. Husistein, Shobhit S. Chaturvedi, Thomas Weymuth, Beatriz Piniello, Sélène Forget, Kresten Lindorff-Larsen, Gemma C. Solomon, Markus Reiher
Abstract
Within a hybrid quantum mechanics/molecular mechanics approach a small region of a large atomistic system is described quantum mechanically, while the rest of the system is represented by a classical force field. Despite this, the computational cost of the quantum calculation severely limits the accessible simulation and sampling time. Replacing the quantum calculation with a machine learning interatomic potential (MLIP) trained on quantum reference data can greatly increase the accessible simulation time. However, the accuracy and efficiency strongly depend on the specific way in which the MLIP incorporates information about the classical (electrostatic) environment. Here, we demonstrate that existing strategies based purely on features derived from an external electric field fail to distinguish certain environments, ultimately limiting their accuracy. As a solution, we propose an extension to the equivariant TensorNet architecture that mitigates this shortcoming by incorporating classical atoms close to the quantum region into the message-passing scheme, combined with features derived from the electric field for describing the environment farther away. We show that this combination achieves higher accuracy than either approach alone.
AI Summary
The paper addresses how machine-learning interatomic potentials can represent the classical electrostatic environment when replacing the quantum calculation in QM/MM simulations. It extends the equivariant TensorNet architecture by combining explicit message passing through nearby MM atoms with the electrostatic potential, electric field, and field gradient evaluated at QM atoms using all MM charges. Only QM-node features contribute directly to the predicted energy. Symmetric point-charge constructions, evaluated using CCSD(T)/aug-cc-pVQZ calculations for argon, demonstrate that truncated local field descriptors can coincide despite different interaction energies. Explicitly resolving nearby charges can distinguish these environments when they lie inside the graph cutoff. Ablations on capped arginine–glutamate in water, using 8,000 training/validation structures and 931 test structures, show energy/force MAEs of 1.94 kJ/mol and 1.15 kJ/mol/Å for the combined representation, versus 4.69/1.92 for field features alone and 4.32/2.25 for nearby MM atoms alone. Results are averaged over three seeds. Subtracting fixed-charge Coulomb interactions before training and restoring them at inference removes statistically detectable distance-dependent projected-force bias in the reported bins. Benchmarks on solvated benzene and uracil give test energy MAEs of 0.55 and 3.54 kJ/mol, respectively, improving on the reported non-delta-learning comparators. The representation still truncates far-environment field information and explicitly resolves only nearby MM atoms. Direct comparison with FieldMACE is withheld because matching data splits are unavailable; uracil errors change sharply between chronological and random splits, attributed to near-duplicate trajectory frames.
Error found: The second symmetry probe incorrectly attributes vanishing potential to symmetric placement alone: it also requires the unstated condition q+/d1=−q−/d2, although the electric field does vanish by symmetry.
Written by gpt-6.1-sol from the paper's full text. Automated: it can be wrong.
- Comments
- 23 pages, 3 figures