REVIEW 5 major objections 4 minor 104 references
Density-Functionalized QM/MM Delivers Chemical Accuracy For Solvated Systems
T0 review · 5 major / 4 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read Giving the force-field region a reconstructed electron density, and treating the QM-MM interaction with orbital-free density functionals, makes QM/MM reach chemical accuracy in solvated systems using only the first solvation shell as the…
desk verdict Worth a serious referee, but the chemical-accuracy claim currently runs ahead of the evidence: the fitted Gaussian widths are missing from the preprint and nearly all validations share the same water/O–H manifold. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the reconstructed MM electron density of Eq. (6): each permanent charge becomes a normalized Gaussian and each induced dipole becomes the gradient of a Gaussian, with element-specific widths fitted so that QM/MM interaction energies reproduce subsystem-DFT results for a single water molecule in bulk water. On top of that density, the nonadditive energy functional of Eq. (5), containing the nonadditive kinetic ($T_s^{\mathrm{nad}}$) and exchange-correlation ($E_{xc}^{\mathrm{nad}}$) functionals, supplies the embedding potential that polarizes the QM subsystem self-consistently with the MM dipoles through Eqs. (7)-(8). In the current implementation only the classical electrostatic terms of that embedding potential act on the MM dipoles, while the nonadditive terms act on the QM density; the fitted Gaussian widths and the consistency of the force field with the QM functional carry the accuracy of the interface.
What would settle it
Apply the published Gaussian widths to a system whose environment is not bulk water, such as a chloride or lithium ion in water or a water/methanol mixture, and compare the QM/MM interaction energy against a full subsystem-DFT reference as the QM region grows from the first shell to several shells. If the error stays outside roughly 1 kcal/mol with only the first solvation shell in QM, or if the per-boundary error grows significantly beyond 0.2-0.3 kcal/mol, the density ansatz and its transferability claim are falsified.
Extended reading notes
Core claim
The central claim is that QM/MM can be reformulated as a fully density-functional theory of two interacting subsystems: the QM subsystem is described by orbital-based DFT, while the MM subsystem is assigned a valence electron density reconstructed as a sum of Gaussians centered on each classical charge and dipole site. The QM-MM interaction is then evaluated with subsystem-DFT nonadditive functionals, including the nonadditive kinetic energy and exchange-correlation, so that electrostatics, charge penetration, exchange, correlation, and Pauli repulsion are all described from first principles rather than by ad hoc potentials. The authors report that the method reproduces subsystem-DFT reference interaction energies for water dimers within about 0.2 kcal/mol, matches bulk-water dipole moments with 0.29 D root mean square error, and achieves the +/-2 kcal/mol convergence target with roughly 13 to 14 QM water molecules for both a neutral solute and a doubly charged aqua ion; with first-shell QM regions the errors drop below about 1 kcal/mol. They conclude that mutual polarization at a density-based QM-MM interface, combined with a force field whose electrostatics are consistent with the QM functional, is what makes the fast convergence possible.
Load-bearing premise
The load-bearing premise is that a smooth Gaussian blob with a fitted width, placed on each classical charge and dipole, really captures the electron density of the environment where it matters, especially the Pauli repulsion and charge penetration at the QM-MM interface; if those shapes are wrong in a new environment, the claimed accuracy and fast convergence collapse.
Editorial extensions
If this is right
- If the central claim holds, QM/MM studies of aqueous solutes no longer need large QM regions: the first solvation shell suffices for chemical accuracy, cutting the cost of condensed-phase simulations.
- Because the MM density is generated from the same force-field charges and induced dipoles already used in simulation, the method can be dropped into existing polarizable force-field workflows without changing the MM model.
- The reported ~0.2-0.3 kcal/mol error per QM-MM boundary implies a simple rule of thumb for QM-region construction: estimate the total boundary error from the number of QM-MM interfaces, and enlarge the QM region only until that estimate falls below the target.
- Mutual polarization captured through the density interface fixes the known failure of electrostatic-embedding QM/MM in bulk water, where neglecting MM polarization raises the interaction-energy RMSE to about 7 kcal/mol.
Reading between the lines
- A testable extension the authors do not report: run the same fitted Gaussian widths on an ionic solute with a tightly bound first shell, such as Li+ or Cl- in water, and check whether the per-boundary error stays near 0.2-0.3 kcal/mol; if it does, the interface ansatz transfers beyond the oxygen and hydrogen sites it was fitted on.
- Since the current MM embedding potential drops the nonadditive kinetic and exchange-correlation terms when polarizing MM dipoles, a natural next step is to include those terms and re-fit the dipole self-energy; this could reduce the remaining dipole RMSE and improve the polarization density near the QM-MM interface.
- The linear relation between partition RMSE and the square root of the number of partition members suggests an uncertainty estimator for arbitrary QM/MM partitions: count the boundaries, multiply by the per-boundary error, and use the result to choose the QM region before running the expensive calculation.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript introduces a QM/MM method in which the MM subsystem is assigned a smooth Gaussian-based electron density (Eq. 6) and the QM-MM interaction is evaluated with orbital-free density functionals for the nonadditive kinetic, exchange-correlation, Coulomb, and Pauli-repulsion terms, using the sDFT framework implemented in eDFTpy. The MM force fields are MB-PBE and MB-Pol, and the QM-MM coupling includes mutual polarization. The method is tested on the water dimer, water hexamer, bulk water, first-shell water pentamers, solvated glucose, a Pd aqua ion, and a wet MoS2 monolayer, all benchmarked against sDFT (PBE-based) reference calculations. The central claims are that the approach reaches chemical accuracy for solvated systems and that including only the first solvation shell in the QM region is sufficient.
Significance. If the central claims are correct, this is a noteworthy advance for QM/MM: the paper presents a clear SCF workflow, a plausible route to include Pauli repulsion and charge penetration at the QM-MM interface, open-source software components, and a broader pilot validation set than many QM/MM method papers. The comparison with mutually polarizable QM/MM (AMOEBA) benchmarks in Ref. 41 is useful and the reported dimer and bulk-water errors are competitive. However, the paper's own equations and fitting procedure make the transferability claim load-bearing: the Gaussian widths and dipole self-energy constants are fitted to bulk-water sDFT interaction energies, and the fitted parameter values are not reported. Most validation systems are water or aqueous environments, so the evidence for transferability to non-aqueous or charged interfaces is currently weak, and the MoS2 result itself shows a factor-of-two deviation from the sDFT reference.
major comments (5)
- [Parameters Defining the MM Density] The fitted parameters are not reported. The text states that four parameters (sigma for O and H, for MB-PBE and MB-Pol, plus the k_SE self-energy constants) were fitted to reproduce sDFT interaction energies for 640 bulk-water snapshots, but the table containing the final values is empty ("Table ??"). Without these values, the results are not reproducible, and the reader cannot assess how strongly the subsequent water dimer, hexamer, bulk-water, and first-shell results are in-sample tests of the fit rather than independent validations.
- [Results and Discussion, Water Dimer and Bulk Water] The validation set is largely on the same manifold as the fitting data. The Gaussian widths are fitted to bulk-water sDFT interaction energies, and the dimer, hexamer, pentamer, bulk-water, glucose/water, Pd-aqua-ion/water, and wet-MoS2 systems all involve only O and H MM sites in hydrogen-bonded aqueous environments. A genuinely out-of-sample test, such as a non-aqueous solvent or an interface with different chemistry, is needed to support the claim that the fitted Gaussian densities are physically transferable rather than system-specific corrections to the PBE/revAPBEk functional errors.
- [Wet Surfaces] The MoS2 result contradicts the uniform-accuracy claim. The paper reports a QM/MM interaction energy of -0.11 kcal/mol per square angstrom versus a QM/QM reference of -0.23 kcal/mol per square angstrom, a factor-of-two discrepancy. This is not chemical accuracy, and it shows that the method's performance is not uniform even within the water-family benchmarks; the manuscript should either explain this discrepancy quantitatively or qualify the accuracy claim.
- [Details of the implementation and Eqs. (7)-(8)] The MM embedding potential used to polarize the MM dipoles retains only the electrostatic terms and neglects the nonadditive kinetic and exchange-correlation contributions, as explicitly acknowledged in the text. This means the MM density is not variationally relaxed with respect to the very nonadditive terms that the method advertises. The manuscript frames this as future work, but the current wording in the abstract and conclusions, which states that the method accounts for exchange, correlation, and Pauli repulsion, should be qualified to make clear that these terms enter the QM-MM energy but not the MM dipole response. A quantitative test of the neglected terms, e.g., by comparing MM polarization densities against sDFT for a charged or non-hydrogen-bonding interface, would strengthen the central claim.
- [Conclusions and Abstract] The phrase "chemical accuracy" is used relative to PBE-based sDFT reference calculations, not relative to experiment or high-level wavefunction theory. Since the MB-Pol force field is fitted to CCSD(T) data, some comparisons are implicitly high-level, but the central benchmarks are sDFT. The manuscript should state this qualification explicitly, because "chemical accuracy" generally implies agreement with experimental or converged quantum-chemical results, and the current evidence does not establish accuracy relative to those references.
minor comments (4)
- [Details of the implementation] In the paragraph describing the QM electrostatic field, the text says the field is derived from "the MM embedding potential in Eq. (7)", but the MM embedding potential is defined in Eq. (8); Eq. (7) defines the QM embedding potential. This cross-reference should be corrected.
- [Throughout] Several tables and figures are missing or referenced as empty or undefined (e.g., Table ??, Figures ?? and ??, supplementary Figure ??). The manuscript is not self-contained as submitted; these placeholders must be filled before the paper can be evaluated or reproduced.
- [Results and Discussion, Water Hexamer] The hexamer results are reported as RMSEs and as an inferred error per QM/MM boundary, but the actual RMSE values in Figures ?? and ?? are not visible in the text. Reporting the numerical values in the text would make the comparison with the dimer and bulk results more direct.
- [Parameters Defining the MM Density] The text introducing the dipole self-energy correction states that the correction is k_SE times the squared difference between the induced dipole and the dipole in the isolated MM subsystem, but it does not specify how k_SE is fitted or whether it is determined separately for MB-PBE and MB-Pol. This should be clarified, especially because this correction directly affects the MM response at the QM-MM interface.
Circularity Check
Bulk-water and first-shell validations reuse the exact 640-structure training set used to fit the Gaussian widths; the chemical-accuracy claim rests partly on in-sample error, though out-of-sample solute/surface tests give independent content.
-
fitted input called prediction
[Section 'Parameters Defining the MM Density' and Results 'Bulk and First Solvation Shell Water Environment' (Eq. 6, Fig. 4)]
"The parameters defining the MM density were fitted so as to reproduce the QM-QM interaction energies for a single water molecule in bulk water. Ten snapshots of 64-water molecule cubic systems were taken from Ref. 94. These provided 640 water-bulk interaction energies."
The four σi widths (Eq. 6) and kSE self-energy constants are optimized on the 640 bulk-water snapshots to reproduce sDFT single-water interaction energies; the paper's 'forward mapping' says the widths are chosen to give the best match between QM/MM and QM/QM interaction energies. The bulk validation in Fig. 4 reports RMSEs for exactly the same quantity on the same 10 snapshots/640 structures (partitions 1/63, 1/4, 5/59), so the RMSEs 1.32, 1.05, 1.11 kcal/mol are in-sample training errors, not predictions. The first-shell pentamers are also cut from these same snapshots with the same O/H MM sites, so the 'first solvation shell suffices' conclusion is partly fit-determined.
full rationale
The most load-bearing 'prediction'—chemical accuracy for water in bulk and for first-shell water environments—is evaluated on the very 640 bulk-water interaction energies used to fit the element-specific Gaussian widths and dipole self-energy constants in Eq. (6). Those RMSEs are therefore training-set errors, not independent validation. The dimer, prism hexamer, glucose, [Pd(H2O)4]2+, and wet-MoS2 tests are genuinely out-of-sample and give the method independent content, but they all use the same water O/H MM sites, so they cannot establish transferability beyond the fitted water manifold. The sDFT benchmark itself is supported in part by prior Pavanello-group citations, yet the QM/QM reference calculations are performed here with standard functionals, so that is normal citation rather than a circular derivation. The missing fitted parameter values (Table ??) and the factor-of-two discrepancy at wet MoS2 are reproducibility/correctness concerns, not additional circularity. Overall, the central water bulk/first-shell claim partially reduces to the fitting objective, while the method retains substantial independent content elsewhere; hence score 6.
Assumptions & free parameters
free parameters (3)
- Gaussian width for MM oxygen (sigma_O) =
Not given (Table ?? in supplementary)
- Gaussian width for MM hydrogen (sigma_H) =
Not given (Table ?? in supplementary)
- Dipole self-energy constants (k_SE) =
Not given
assumptions (5)
- domain assumption The total electron density is the sum of QM and MM subsystem densities (Eq. 1).
- standard math The energy is decomposable into additive subsystem energies plus a nonadditive functional (Eq. 2).
- domain assumption The MB-PBE and MB-Pol force fields faithfully represent the underlying DFT (PBE) or CCSD(T) interaction energies for water.
- ad hoc to paper The reconstructed Gaussian density (Eq. 6) with fitted widths captures the essential short-range physics (Pauli repulsion, charge penetration) at the QM-MM interface.
- ad hoc to paper Neglecting the nonadditive kinetic and exchange-correlation terms in the MM embedding potential (Eq. 7) does not qualitatively corrupt the MM dipole response.
Cite this review
Pith. "Pith review of Density-Functionalized QM/MM Delivers Chemical Accuracy For Solvated Systems." pith.science (2026). https://pith.science/paper/FE6MKQE5
@misc{pith2026241117844,
author = {Pith},
title = {Pith review of: Density-Functionalized QM/MM Delivers Chemical Accuracy For Solvated Systems},
year = {2026},
howpublished = {\url{https://pith.science/paper/FE6MKQE5}},
note = {Machine review of arXiv:2411.17844}
}
abstract
We present a reformulation of QM/MM as a fully quantum mechanical theory of interacting subsystems, all treated at the level of density functional theory (DFT). For the MM subsystem, which lacks orbitals, we assign an ad hoc electron density and apply orbital-free DFT functionals to describe its quantum properties. The interaction between the QM and MM subsystems is also treated using orbital-free density functionals, accounting for Coulomb interactions, exchange, correlation, and Pauli repulsion. Consistency across QM and MM subsystems is ensured by employing data-driven, many-body MM force fields that faithfully represent DFT functionals. Applications to water-solvated systems demonstrate that this approach achieves unprecedented, very rapid convergence to chemical accuracy as the size of the QM subsystem increases. We validate the method with several pilot studies, including water bulk, water clusters (prism hexamer and pentamers), solvated glucose, a palladium aqua ion, and a wet monolayer of MoS$_2$.
Figures
Figures from the paper (3 more)
Reference graph
Works this paper leans on
-
[1]
Development of Multiscale Models for Complex Chemical Systems: From H+H _2 to Biomolecules
Karplus, M. Development of Multiscale Models for Complex Chemical Systems: From H+H _2 to Biomolecules. Nobel Lecture, 2013; https://www.nobelprize.org/uploads/2018/06/karplus-lecture.pdf
2013
-
[2]
Implications of torsional potential of retinal isomers for visual excitation
Honig, B.; Karplus, M. Implications of torsional potential of retinal isomers for visual excitation. Nature 1971, 229, 558--560
1971
-
[3]
M.; Thiel, W
Senn, H. M.; Thiel, W. QM/MM studies of enzymes. Current Opinion in Chemical Biology 2007, 11, 182--187
2007
-
[4]
Kulik, H. J. Large-scale QM/MM free energy simulations of enzyme catalysis reveal the influence of charge transfer. Phys. Chem. Chem. Phys. 2018, 20, 20650--20660
2018
-
[5]
Drug Design in the Exascale Era: A Perspective from Massively Parallel QM/MM Simulations
Raghavan, B.; Paulikat, M.; Ahmad, K.; Callea, L.; Rizzi, A.; Ippoliti, E.; Mandelli, D.; Bonati, L.; De Vivo, M.; Carloni, P. Drug Design in the Exascale Era: A Perspective from Massively Parallel QM/MM Simulations. Journal of Chemical Information and Modeling 2023, 63, 3647--3658, PMID: 37319347
2023
-
[6]
W.; Wang, B.; Wu, X.-P.; Zheng, S.-W.; Gagliardi, L.; Truhlar, D
Lin, H.; Zhang, Y.; Pezeshki, S.; Duster, A. W.; Wang, B.; Wu, X.-P.; Zheng, S.-W.; Gagliardi, L.; Truhlar, D. G. QMMM 2023: A program for combined quantum mechanical and molecular mechanical modeling and simulations. Computer Physics Communications 2024, 295, 108987
2023
-
[7]
QMMMW: A wrapper for QM/MM simulations with Quantum ESPRESSO and LAMMPS
Ma, C.; Martin-Samos, L.; Fabris, S.; Laio, A.; Piccinin, S. QMMMW: A wrapper for QM/MM simulations with Quantum ESPRESSO and LAMMPS. Computer Physics Communications 2015, 195, 191--198
2015
-
[8]
A.; Paesani, F
Lambros, E.; Lipparini, F.; Cisneros, G. A.; Paesani, F. A Many-Body, Fully Polarizable Approach to QM/MM Simulations. Journal of Chemical Theory and Computation 2020, 16, 7462--7472
2020
Show all 104 references
-
[9]
A.; Darden, T
Gresh, N.; Cisneros, G. A.; Darden, T. A.; Piquemal, J.-P. Anisotropic, Polarizable Molecular Mechanics Studies of Inter- and Intramolecular Interactions and Ligand-Macromolecule Complexes. A Bottom-Up Strategy. Journal of Chemical Theory and Computation 2007, 3, 1960--1986
2007
-
[10]
Interfacing ab Initio Quantum Mechanical Method with Classical Drude Osillator Polarizable Model for Molecular Dynamics Simulation of Chemical Reactions
Lu, Z.; Zhang, Y. Interfacing ab Initio Quantum Mechanical Method with Classical Drude Osillator Polarizable Model for Molecular Dynamics Simulation of Chemical Reactions. Journal of Chemical Theory and Computation 2008, 4, 1237--1248
2008
-
[11]
Importance of MM Polarization in QM/MM Studies of Enzymatic Reactions: Assessment of the QM/MM Drude Oscillator Model
Ganguly, A.; Boulanger, E.; Thiel, W. Importance of MM Polarization in QM/MM Studies of Enzymatic Reactions: Assessment of the QM/MM Drude Oscillator Model. Journal of Chemical Theory and Computation 2017, 13, 2954--2961
2017
-
[12]
Hybrid QM/MM Molecular Dynamics with AMOEBA Polarizable Embedding
Loco, D.; Lagard è re, L.; Caprasecca, S.; Lipparini, F.; Mennucci, B.; Piquemal, J.-P. Hybrid QM/MM Molecular Dynamics with AMOEBA Polarizable Embedding. Journal of Chemical Theory and Computation 2017, 13, 4025--4033
2017
-
[13]
Polarizable embedding QM/MM: the future gold standard for complex (bio)systems? Phys
Bondanza, M.; Nottoli, M.; Cupellini, L.; Lipparini, F.; Mennucci, B. Polarizable embedding QM/MM: the future gold standard for complex (bio)systems? Phys. Chem. Chem. Phys. 2020, 22, 14433--14448
2020
-
[14]
Fast approximate but accurate QM/MM interactions for polarizable embedding
Reinholdt, P.; Kongsted, J.; Lipparini, F. Fast approximate but accurate QM/MM interactions for polarizable embedding. J. Chem. Theory Comput. 2021, 18, 344--356
2021
-
[15]
Polarizable QM/MM Approach with Fluctuating Charges and Fluctuating Dipoles: The QM/FQF Model
Giovannini, T.; Puglisi, A.; Ambrosetti, M.; Cappelli, C. Polarizable QM/MM Approach with Fluctuating Charges and Fluctuating Dipoles: The QM/FQF Model. Journal of Chemical Theory and Computation 2019, 15, 2233--2245, PMID: 30875213
2019
-
[16]
Electrostatic embedding of machine learning potentials
Zinovjev, K. Electrostatic embedding of machine learning potentials. J. Chem. Theory Comput. 2023, 19, 1888--1897
2023
-
[17]
J.; Zeng, J.; Lerew, L.; McCarthy, E.; Tao, Y.; Ekesan, S.; York, D
Giese, T. J.; Zeng, J.; Lerew, L.; McCarthy, E.; Tao, Y.; Ekesan, S.; York, D. M. Software Infrastructure for Next-Generation QM/MM- MLP Force Fields. J. Phys. Chem. B 2024,
2024
-
[18]
E.; Mohr, S.; Huhs, G.; Deutsch, T.; Masella, M.; Genovese, L
Ratcliff, L. E.; Mohr, S.; Huhs, G.; Deutsch, T.; Masella, M.; Genovese, L. Challenges in large scale quantum mechanical calculations. Wiley Interdisciplinary Reviews: Computational Molecular Science 2017, 7, e1290
2017
-
[19]
Introduction to QM/MM simulations
Groenhof, G. Introduction to QM/MM simulations. Biomolecular simulations: methods and protocols 2013, 43--66
2013
-
[20]
Theoretical studies of enzymic reactions: dielectric, electrostatic and steric stabilization of the carbonium ion in the reaction of lysozyme
Warshel, A.; Levitt, M. Theoretical studies of enzymic reactions: dielectric, electrostatic and steric stabilization of the carbonium ion in the reaction of lysozyme. Journal of molecular biology 1976, 103, 227--249
1976
-
[21]
General linear scaling implementation of polarizable embedding schemes
Lipparini, F. General linear scaling implementation of polarizable embedding schemes. Journal of Chemical Theory and Computation 2019, 15, 4312--4317
2019
-
[22]
The OpenMMPol library for polarizable QM/MM calculations of properties and dynamics
Bondanza, M.; Nottoli, T.; Nottoli, M.; Cupellini, L.; Lipparini, F.; Mennucci, B. The OpenMMPol library for polarizable QM/MM calculations of properties and dynamics. J. Chem. Phys. 2024, 160
2024
-
[23]
Olsen, J. M. H.; Bolnykh, V.; Meloni, S.; Ippoliti, E.; Bircher, M. P.; Carloni, P.; Rothlisberger, U. MiMiC: a novel framework for multiscale modeling in computational chemistry. Journal of chemical theory and computation 2019, 15, 3810--3823
2019
-
[24]
Cruzeiro, V. W. D.; Manathunga, M.; Merz Jr, K. M.; G \"o tz, A. W. Open-source multi-GPU-accelerated QM/MM simulations with AMBER and QUICK. Journal of Chemical Information and Modeling 2021, 61, 2109--2115
2021
-
[25]
P.; McDaniel, J
Pederson, J. P.; McDaniel, J. G. PyDFT-QMMM: A modular, extensible software framework for DFT-based QM/MM molecular dynamics. J. Chem. Phys. 2024, 161
2024
-
[26]
TINKTEP: A fully self-consistent, mutually polarizable QM/MM approach based on the AMOEBA force field
Dziedzic, J.; Mao, Y.; Shao, Y.; Ponder, J.; Head-Gordon, T.; Head-Gordon, M.; Skylaris, C.-K. TINKTEP: A fully self-consistent, mutually polarizable QM/MM approach based on the AMOEBA force field. J. Chem. Phys. 2016, 145
2016
-
[27]
S.; Purton, J
Lu, Y.; Sen, K.; Yong, C.; Gunn, D. S.; Purton, J. A.; Guan, J.; Desmoutier, A.; Nasir, J. A.; Zhang, X.; Zhu, L., et al. Multiscale QM/MM modelling of catalytic systems with ChemShell. Phys. Chem. Chem. Phys. 2023, 25, 21816--21835
2023
-
[28]
O.; Levi, G.; Mortensen, J
Dohn, A.; Jonsson, E. O.; Levi, G.; Mortensen, J. J.; Lopez-Acevedo, O.; Thygesen, K. S.; Jacobsen, K. W.; Ulstrup, J.; Henriksen, N. E.; M ller, K., et al. Grid-based projector augmented wave (GPAW) implementation of quantum mechanics/molecular mechanics (QM/MM) electrostatic...
2017
-
[29]
Nam, K.; Gao, J.; York, D. M. An Efficient Linear-Scaling Ewald Method for Long-Range Electrostatic Interactions in Combined QM/MM Calculations. J. Chem. Theory Comput. 2005, 1, 2--13, PMID: 26641110
2005
-
[30]
C.; Richard, R
Holden, Z. C.; Richard, R. M.; Herbert, J. M. Periodic boundary conditions for QM/MM calculations: Ewald summation for extended Gaussian basis sets . J. Chem. Phys. 2013, 139, 244108
2013
-
[31]
Analytic Gradients for the Electrostatic Embedding QM/MM Model in Periodic Boundary Conditions Using Particle-Mesh Ewald Sums and Electrostatic Potential Fitted Charge Operators
Bonfrate, S.; Ferr \'e , N.; Huix-Rotllant, M. Analytic Gradients for the Electrostatic Embedding QM/MM Model in Periodic Boundary Conditions Using Particle-Mesh Ewald Sums and Electrostatic Potential Fitted Charge Operators. J. Chem. Theory Comput. 2024, 20, 4338--4349
2024
-
[32]
P.; McDaniel, J
Pederson, J. P.; McDaniel, J. G. DFT-based QM/MM with particle-mesh Ewald for direct, long-range electrostatic embedding. J. Chem. Phys. 2022, 156
2022
-
[33]
An efficient linear-scaling electrostatic coupling for treating periodic boundary conditions in QM/MM simulations
Laino, T.; Mohamed, F.; Laio, A.; Parrinello, M. An efficient linear-scaling electrostatic coupling for treating periodic boundary conditions in QM/MM simulations. J. Chem. Theory Comput. 2006, 2, 1370--1378
2006
-
[34]
An efficient real space multigrid QM/MM electrostatic coupling
Laino, T.; Mohamed, F.; Laio, A.; Parrinello, M. An efficient real space multigrid QM/MM electrostatic coupling. J. Chem. Theory Comput. 2005, 1, 1176--1184
2005
-
[35]
Brandt, F.; Jacob, C. R. Systematic QM region construction in QM/MM calculations based on uncertainty quantification. J. Chem. Theory Comput. 2022, 18, 2584--2596
2022
-
[36]
Pedraza-Gonz \'a lez, L.; Mar \' n, M. D. C.; Jorge, A. N.; Ruck, T. D.; Yang, X.; Valentini, A.; Olivucci, M.; De Vico, L. Web-ARM: a web-based interface for the automatic construction of QM/MM models of rhodopsins. J. Chem. Inf. Model. 2020, 60, 1481--1493
2020
-
[37]
J.; Zhang, J.; Klinman, J
Kulik, H. J.; Zhang, J.; Klinman, J. P.; Mart \' nez, T. J. How large should the QM region be in QM/MM calculations? The case of catechol O-methyltransferase. J. Phys. Chem. B 2016, 120, 11381--11394
2016
-
[38]
Karelina, M.; Kulik, H. J. Systematic quantum mechanical region determination in QM/MM simulation. J. Chem. Theory Comput. 2017, 13, 563--576
2017
-
[39]
J.; Mandado, M
P \'e rez-Barcia, \'A .; C \'a rdenas, G.; Nogueira, J. J.; Mandado, M. Effect of the QM size, basis set, and polarization on QM/MM interaction energy decomposition analysis. J. Chem. Inf. Model. 2023, 63, 882--897
2023
-
[40]
Mutually polarizable QM/MM model with in situ optimized localized basis functions
Dziedzic, J.; Head-Gordon, T.; Head-Gordon, M.; Skylaris, C.-K. Mutually polarizable QM/MM model with in situ optimized localized basis functions . The Journal of Chemical Physics 2019, 150, 074103
2019
-
[41]
Unraveling Water Solvation Effects with Quantum Mechanics/Molecular Mechanics Semiclassical Vibrational Spectroscopy: The Case of Thymidine
Moscato, D.; Mandelli, G.; Bondanza, M.; Lipparini, F.; Conte, R.; Mennucci, B.; Ceotto, M. Unraveling Water Solvation Effects with Quantum Mechanics/Molecular Mechanics Semiclassical Vibrational Spectroscopy: The Case of Thymidine. J. Am. Chem. Soc. 2024, 146, 8179--8188
2024
-
[42]
Multiple facets of modeling electronic absorption spectra of systems in solution
G \'o mez, S.; Giovannini, T.; Cappelli, C. Multiple facets of modeling electronic absorption spectra of systems in solution. ACS Physical Chemistry Au 2022, 3, 1--16
2022
-
[43]
P.; Herbert, J
Coons, M. P.; Herbert, J. M. Quantum chemistry in arbitrary dielectric environments: Theory and implementation of nonequilibrium Poisson boundary conditions and application to compute vertical ionization energies at the air/water interface. Journal Chem. Phys. 2018, 148
2018
-
[44]
A.; Paetow, L.; Tolle, J.; Shao, X.; Ramos, P.; Neugebauer, J.; Pavanello, M
Martinez B, J. A.; Paetow, L.; Tolle, J.; Shao, X.; Ramos, P.; Neugebauer, J.; Pavanello, M. Which Physical Phenomena Determine the Ionization Potential of Liquid Water? J. Phys. Chem. B 2023, 127, 5470--5480
2023
-
[45]
Olsen, J. M. H.; Steinmann, C.; Ruud, K.; Kongsted, J. Polarizable density embedding: A new QM/QM/MM-based computational strategy. J. Phys. Chem. A 2015, 119, 5344--5355
2015
-
[46]
B.; Olsen, J
Kvedaraviciute, S.; Carrasco-Busturia, D.; M ller, K. B.; Olsen, J. M. H. Polarizable Embedding without Artificial Boundary Polarization. J. Chem. Theory Comput. 2023, 19, 5122--5141
2023
-
[47]
I.; Slipchenko, L
Viquez Rojas, C. I.; Slipchenko, L. V. Exchange repulsion in quantum mechanical/effective fragment potential excitation energies: Beyond polarizable embedding. J. Chem. Theory Comput. 2020, 16, 6408--6417
2020
-
[48]
A.; Piquemal, J.-P.; Cisneros, G
Gokcan, H.; Kratz, E.; Darden, T. A.; Piquemal, J.-P.; Cisneros, G. A. QM/MM simulations with the Gaussian electrostatic model: A density-based polarizable potential. J. Phys. Chem. Lett. 2018, 9, 3062--3067
2018
-
[49]
R.; Giese, T
Kuechler, E. R.; Giese, T. J.; York, D. M. Charge-dependent many-body exchange and dispersion interactions in combined QM/MM simulations. J. Chem. Phys. 2015, 143
2015
-
[50]
S.; Hattig, C.; Hofener, S
Tre , R. S.; Hattig, C.; Hofener, S. Employing pseudopotentials to tackle excited-state electron spill-out in frozen density embedding calculations. J. Chem. Theory Comput. 2022, 18, 1737--1747
2022
-
[51]
Energy benchmarks for water clusters and ice structures from an embedded many-body expansion
Gillan, M.; Alf \`e , D.; Bygrave, P.; Taylor, C.; Manby, F. Energy benchmarks for water clusters and ice structures from an embedded many-body expansion. J. Chem. Phys. 2013, 139
2013
-
[52]
Ben-Nun, M.; Mart \'i nez, T. J. Direct evaluation of the Pauli repulsion energy usingclassical'wavefunctions in hybrid quantum/classical potential energy surfaces. Chem. Phys. Lett. 1998, 290, 289--295
1998
-
[53]
A General Route to Include Pauli Repulsion and Quantum Dispersion Effects in QM/MM Approaches
Giovannini, T.; Lafiosca, P.; Cappelli, C. A General Route to Include Pauli Repulsion and Quantum Dispersion Effects in QM/MM Approaches. Journal of Chemical Theory and Computation 2017, 13, 4854--4870, PMID: 28898079
2017
-
[54]
Bedrov, D.; Piquemal, J.-P.; Borodin, O.; MacKerell, A. D. J.; Roux, B.; Schröder, C. Molecular Dynamics Simulations of Ionic Liquids and Electrolytes Using Polarizable Force Fields. Chemical Reviews 2019, 119, 7940--7995, PMID: 31141351
2019
-
[55]
Hr s ak, D.; Olsen, J. M. H.; Kongsted, J. Optimization and transferability of non-electrostatic repulsion in the polarizable density embedding model. Journal of computational chemistry 2017, 38, 2108--2117
2017
-
[56]
Subsystem real-time time dependent density functional theory
Krishtal, A.; Ceresoli, D.; Pavanello, M. Subsystem real-time time dependent density functional theory. J. Chem. Phys. 2015, 142
2015
-
[57]
R.; Neugebauer, J
Jacob, C. R.; Neugebauer, J. Subsystem density-functional theory. WIREs: Comp. Mol. Sci. 2014, 4, 325--362
2014
-
[58]
A.; Shedge, S.; Zhou, X
Wesolowski, T. A.; Shedge, S.; Zhou, X. Frozen-density embedding strategy for multilevel simulations of electronic structure. Chem. Rev. 2015, 115, 5891--5928
2015
-
[59]
Orbital-free density functional theory: An attractive electronic structure method for large-scale first-principles simulations
Mi, W.; Luo, K.; Trickey, S.; Pavanello, M. Orbital-free density functional theory: An attractive electronic structure method for large-scale first-principles simulations. Chem. Rev. 2023, 123, 12039--12104
2023
-
[60]
R.; Neugebauer, J
Jacob, C. R.; Neugebauer, J. Subsystem density-functional theory (update). Wiley Interdisciplinary Reviews: Computational Molecular Science 2024, 14, e1700
2024
-
[61]
See supplementary materials document at [url-to-be-inserted]
Supplementary Materials. See supplementary materials document at [url-to-be-inserted]
-
[62]
Avoiding fractional electrons in subsystem DFT based ab-initio molecular dynamics yields accurate models for liquid water and solvated OH radical
Genova, A.; Ceresoli, D.; Pavanello, M. Avoiding fractional electrons in subsystem DFT based ab-initio molecular dynamics yields accurate models for liquid water and solvated OH radical. J. Chem. Phys. 2016, 144
2016
-
[63]
Ab initio structure and dynamics of CO2 at supercritical conditions
Mi, W.; Ramos, P.; Maranhao, J.; Pavanello, M. Ab initio structure and dynamics of CO2 at supercritical conditions. J. Phys. Chem. Lett. 2019, 10, 7554--7559
2019
-
[64]
GGA-level subsystem DFT achieves Sub-kcal/mol accuracy intermolecular interactions by mimicking nonlocal functionals
Shao, X.; Mi, W.; Pavanello, M. GGA-level subsystem DFT achieves Sub-kcal/mol accuracy intermolecular interactions by mimicking nonlocal functionals. J. Chem. Theory Comput. 2021, 17, 3455--3461
2021
-
[65]
Nonlocal subsystem density functional theory
Mi, W.; Pavanello, M. Nonlocal subsystem density functional theory. J. Phys. Chem. Lett. 2019, 11, 272--279
2019
-
[66]
u ns, D.; Klahr, K.; M \
Schl \"u ns, D.; Klahr, K.; M \"u ck-Lichtenfeld, C.; Visscher, L.; Neugebauer, J. Subsystem-DFT potential-energy curves for weakly interacting systems. Phys. Chem. Chem. Phys. 2015, 17, 14323--14341
2015
-
[67]
Density embedding method for nanoscale molecule--metal interfaces
Shao, X.; Mi, W.; Pavanello, M. Density embedding method for nanoscale molecule--metal interfaces. J. Phys. Chem. Lett. 2022, 13, 7147--7154
2022
-
[68]
C.; Khan Musa, M
Shao, X.; Lopez, A. C.; Khan Musa, M. R.; Nouri, M. R.; Pavanello, M. Adaptive subsystem density functional theory. J. Chem. Theory Comput. 2022, 18, 6646--6655
2022
-
[69]
Unraveling the Hydration Shell Structure and Dynamics of Group 10 Aqua Ions
Chen, X.; Cifuentes-Lopez, A.; Shao, X.; Lin, L.; Prokopchuk, D.; Pavanello, M. Unraveling the Hydration Shell Structure and Dynamics of Group 10 Aqua Ions. J. Phys. Chem. Lett. 2024, 15, 5517--5528
2024
-
[70]
Schmitt-Monreal, D.; Jacob, C. R. Density-based many-body expansion as an efficient and accurate quantum-chemical fragmentation method: Application to water clusters. Journal of Chemical Theory and Computation 2021, 17, 4144--4156
2021
-
[71]
Focke, K.; Jacob, C. R. Coupled-cluster density-based many-body expansion. The Journal of Physical Chemistry A 2023, 127, 9139--9148
2023
-
[72]
Kevorkyants, R.; Dulak, M.; Wesolowski, T. A. Interaction energies in hydrogen-bonded systems: A testing ground for subsystem formulation of density-functional theory. J. Chem. Phys. 2006, 124
2006
-
[73]
Cisneros, G. A. Application of Gaussian Electrostatic Model (GEM) Distributed Multipoles in the AMOEBA Force Field. J. Chem. Theory Comput. 2012, 8, 5072--5080
2012
-
[74]
First Principles
Babin, V.; Leforestier, C.; Paesani, F. Development of a “First Principles” Water Potential with Flexible Monomers: Dimer Potential Energy Surface, VRT Spectrum, and Second Virial Coefficient. Journal of Chemical Theory and Computation 2013, 9, 5395--5403, PMID: 26592277
2013
-
[75]
First Principles
Babin, V.; Medders, G. R.; Paesani, F. Development of a “First Principles” Water Potential with Flexible Monomers. II: Trimer Potential Energy Surface, Third Virial Coefficient, and Small Clusters. Journal of Chemical Theory and Computation 2014, 10, 1599--1607, PMID: 26580372
2014
-
[76]
T.; Götz, A
Riera, M.; Lambros, E.; Nguyen, T. T.; Götz, A. W.; Paesani, F. Low-order many-body interactions determine the local structure of liquid water. Chem. Sci. 2019, 10, 8211--8218
2019
-
[77]
F.; Bennett, J
Garrity, K. F.; Bennett, J. W.; Rabe, K. M.; Vanderbilt, D. Pseudopotentials for high-throughput DFT calculations. Computational Materials Science 2014, 81, 446--452
2014
-
[78]
O (N N) scaling method to evaluate the ion--electron potential of crystalline solids
Shao, X.; Mi, W.; Xu, Q.; Wang, Y.; Ma, Y. O (N N) scaling method to evaluate the ion--electron potential of crystalline solids. J. Chem. Phys. 2016, 145
2016
-
[79]
eDFTpy: An Object-Oriented Platform for Density Embedding Simulations
Shao, X.; Mi, W.; Pavanello, M. eDFTpy: An Object-Oriented Platform for Density Embedding Simulations. available at http://edftpy.rutgers.edu, 2024; http://edftpy.rutgers.edu, (Last accessed on: 04/25/2024)
2024
-
[80]
QEpy: Quantum ESPRESSO in Python
Shao, X.; Andreussi, O.; Ceresoli, D.; Truscott, M.; Baczewski, A.; Campbell, Q.; Pavanello, M. QEpy: Quantum ESPRESSO in Python . Available at https://gitlab.com/shaoxc/qepy, 2024; https://gitlab.com/shaoxc/qepy, (Last accessed on: 04/25/2024)
2024
-
[81]
Giannozzi, P. et al. Advanced capabilities for materials modelling with QUANTUM ESPRESSO. Journal of Physics: Condensed Matter 2017, 29, 465901
2017
-
[82]
F.; Zhu, X.; Agnew, H.; Smith, D
Riera, M.; Knight, C.; Bull-Vulpe, E. F.; Zhu, X.; Agnew, H.; Smith, D. G. A.; Simmonett, A. C.; Paesani, F. MBX: A many-body energy and force calculator for data-driven many-body simulations . The Journal of Chemical Physics 2023, 159, 054802
2023
-
[83]
N rby, M.; Magnus Haugaard Olsen, J.; Kongsted, J.; Aagard Jensen, H
S. N rby, M.; Magnus Haugaard Olsen, J.; Kongsted, J.; Aagard Jensen, H. J. Multipole moments for embedding potentials: Exploring different atomic allocation algorithms. J. Comp. Chem. 2016, 37, 1887--1896
2016
-
[84]
Ferr \'e , N.; \'A ngy \'a n, J. G. Approximate electrostatic interaction operator for QM/MM calculations. Chem. Phys. Lett. 2002, 356, 331--339
2002
-
[85]
A.; Pavanello, M
Genova, A.; Ceresoli, D.; Krishtal, A.; Andreussi, O.; DiStasio Jr, R. A.; Pavanello, M. eQE: An open-source density functional embedding theory code for the condensed phase. Int. J. Quantum Chem. 2017, 117, e25401
2017
-
[86]
eQE 2.0: Subsystem DFT Beyond GGA Functionals
Mi, W.; Shao, X.; Genova, A.; Ceresoli, D.; Pavanello, M. eQE 2.0: Subsystem DFT Beyond GGA Functionals. Comp. Phys. Comm. 2021, 269, 108122
2021
-
[87]
Exploiting the locality of periodic subsystem density-functional theory: efficient sampling of the Brillouin zone
Genova, A.; Pavanello, M. Exploiting the locality of periodic subsystem density-functional theory: efficient sampling of the Brillouin zone. J. Phys.: Condens. Matter 2015, 27, 495501
2015
-
[88]
P.; Burke, K.; Ernzerhof, M
Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple [Phys. Rev. Lett. 77, 3865 (1996)]. Phys. Rev. Lett. 1997, 78, 1396--1396
1996
-
[89]
Laricchia, S.; Fabiano, E.; Constantin, L.; Della Sala, F. Generalized gradient approximations of the noninteracting kinetic energy from the semiclassical atom theory: Rationalization of the accuracy of the frozen density embedding theory for nonbonded interactions. J. Chem. T...
2011
-
[90]
Supramolecular binding thermodynamics by dispersion-corrected density functional theory
Grimme, S. Supramolecular binding thermodynamics by dispersion-corrected density functional theory. Chemistry--A European Journal 2012, 18, 9955--9964
2012
-
[91]
first principles
Babin, V.; Leforestier, C.; Paesani, F. Development of a “first principles” water potential with flexible monomers: Dimer potential energy surface, VRT spectrum, and second virial coefficient. J. Chem. Theory Comput. 2013, 9, 5395--5403
2013
-
[92]
Theory of Electric Polarization (Second Edition), second edition ed.; Theory of Electric Polarization; Elsevier: Amsterdam, 1973; p iii
BÖTTCHER, C., Ed. Theory of Electric Polarization (Second Edition), second edition ed.; Theory of Electric Polarization; Elsevier: Amsterdam, 1973; p iii
1973
-
[93]
Cooperation and Environment Characterize the Low-Lying Optical Spectrum of Liquid Water
Kumar P., S.; Genova, A.; Pavanello, M. Cooperation and Environment Characterize the Low-Lying Optical Spectrum of Liquid Water. J. Phys. Chem. Lett. 2017, 8, 5077--5083
2017
-
[94]
S.; Leininger, M
Tschumper, G. S.; Leininger, M. L.; Hoffman, B. C.; Valeev, E. F.; Schaefer, I., Henry F.; Quack, M. Anchoring the water dimer potential energy surface with explicitly correlated computations and focal point analyses . The Journal of Chemical Physics 2002, 116, 690--701
2002
-
[95]
o m, J.; Karlstr \
Holt, A.; Bostr \"o m, J.; Karlstr \"o m, G.; Lindh, R. A NEMO potential that includes the dipole--quadrupole and quadrupole--quadrupole polarizability. J. Comp. Chem. 2010, 31, 1583--1591
2010
-
[96]
Inclusion of the quadrupole moment when describing polarization
Holt, A.; Karlstr \"o m, G. Inclusion of the quadrupole moment when describing polarization. The effect of the dipole-quadrupole polarizability. J. Comp. Chem. 2008, 29, 2033--2038
2008
-
[97]
Structural Studies of the Water Hexamer
Hincapié, G.; Acelas, N.; Castaño, M.; David, J.; Restrepo, A. Structural Studies of the Water Hexamer. J. Phys. Chem. A 2010, 114, 7809--7814, PMID: 20590081
2010
-
[98]
Nature 1905, 72, 318--319
Lord Rayleigh , The problem of the random walk. Nature 1905, 72, 318--319
1905
-
[99]
Revealing electronic open quantum systems with subsystem TDDFT
Krishtal, A.; Pavanello, M. Revealing electronic open quantum systems with subsystem TDDFT. J. Chem. Phys. 2016, 144
2016
-
[100]
Local electric dipole moments for periodic systems via density functional theory embedding
Luber, S. Local electric dipole moments for periodic systems via density functional theory embedding. J. Chem. Phys. 2014, 141
2014
-
[101]
Tang, K.; Toennies, J. P. An improved simple model for the van der Waals potential based on universal damping functions for the dispersion coefficients. J. Chem. Phys. 1984, 80, 3726--3741
1984
-
[102]
Tkatchenko, A.; Ambrosetti, A.; DiStasio, R. A. Interatomic methods for the dispersion energy derived from the adiabatic connection fluctuation-dissipation theorem. J. Chem. Phys. 2013, 138
2013
-
[103]
M.; Potestio, R.; Donadio, D.; Bulo, R
Boereboom, J. M.; Potestio, R.; Donadio, D.; Bulo, R. E. Toward Hamiltonian adaptive qm/mm: accurate solvent structures using many-body potentials. J. Chem. Theory Comput. 2016, 12, 3441--3448
2016
-
[104]
yֈk 53 VK@q*c5S ǦW^5 k Kv ;N U
Chen, J.; Harper, J. B.; Ho, J. Improving the accuracy of quantum mechanics/molecular mechanics (QM/MM) models with polarized fragment charges. J. Chem. Theory Comput. 2022, 18, 5607--5617 mcitethebibliography document Bulk-2systems.pdf00006640000000000000000001054466147214213...
2022
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.