Pith. sign in

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 →

arxiv 2411.17844 v1 pith:FE6MKQE5 submitted 2024-11-26 physics.chem-ph physics.comp-ph

classification physics.chem-phphysics.comp-ph
keywords QM/MMsubsystemDFTorbital-freedensityfunctionaltheorypolarizableembeddingchemicalaccuracyaqueoussolvationnonadditivefunctionalsPaulirepulsion
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper claims that the slow convergence of QM/MM calculations on solvated systems can be cured by giving the classical MM region a real electron density and treating the interaction between QM and MM regions with density functionals, not fitted point-charge potentials. The method assigns each MM atom a smooth Gaussian density built from the force field's charges and induced dipoles, then uses orbital-free density functionals for exchange, correlation, and Pauli repulsion across the QM-MM boundary. The authors show that, once the MM density widths are fitted to reproduce subsystem-DFT interaction energies of one water molecule in bulk water, the approach reaches chemical accuracy in water clusters, bulk water, solvated glucose, a palladium aqua ion, and a wet monolayer of molybdenum disulfide. Crucially, the interaction energy converges to the full quantum reference with only the first solvation shell (or fewer) in the QM region. If true, this would make QM/MM simulations of aqueous systems accurate and substantially cheaper than current practice.

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.

Watch

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

Editorial extensions of the paper, not claims the author makes directly.

  • 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.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

5 major / 4 minor

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)
  1. [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.
  2. [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.
  3. [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.
  4. [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.
  5. [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)
  1. [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.
  2. [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.
  3. [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.
  4. [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

1 steps flagged · score 6.0 of 10

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.

  1. 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 3 free parameters · 5 assumptions · 0 invented entities

The method's accuracy rests on three empirical pillars: the polarizable force field (MB-PBE/MB-Pol) fitted to reference data, the Gaussian density reconstruction with fitted widths, and the dipole self-energy correction with fitted constants. These are not derived from first principles, despite the density-functionalized framing. No new physical entities are introduced.

free parameters (3)
  • Gaussian width for MM oxygen (sigma_O) = Not given (Table ?? in supplementary)
    Controls the spatial extent of the reconstructed MM valence density for oxygen; fitted to reproduce QM/QM interaction energies of a water molecule in bulk water.
  • Gaussian width for MM hydrogen (sigma_H) = Not given (Table ?? in supplementary)
    Same fit as for oxygen, applied to hydrogen sites in the reconstructed MM density.
  • Dipole self-energy constants (k_SE) = Not given
    Element-dependent proportionality constants added to the induced dipole self-energy to correct for the non-point nature of the dipoles; the fitting procedure is not detailed in the main text.
assumptions (5)
  • domain assumption The total electron density is the sum of QM and MM subsystem densities (Eq. 1).
    This is the standard subsystem DFT ansatz; it assumes the QM and MM densities are additive and that inter-subsystem correlation is captured entirely by the nonadditive functionals.
  • standard math The energy is decomposable into additive subsystem energies plus a nonadditive functional (Eq. 2).
    This is the standard sDFT decomposition; it relies on the existence of accurate approximations for T_s^nad and E_xc^nad.
  • domain assumption The MB-PBE and MB-Pol force fields faithfully represent the underlying DFT (PBE) or CCSD(T) interaction energies for water.
    The method's consistency argument depends on the MM force field being the classical image of the QM functional. This is a strong assumption for non-water systems.
  • 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.
    The widths are empirical parameters fitted to bulk water interaction energies; no first-principles derivation is provided.
  • 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.
    The authors state they keep only the classical electrostatic terms for polarizing the MM subsystem; this is a stated approximation whose impact is not quantified.

how reviews work

0 comments
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 reproduced from arXiv: 2411.17844 by the authors.

Figure 1
Figure 1. Workflow of the QM/MM method with emphasis on the software implementation. [PITH_FULL_IMAGE:figures/full_fig_p010_1.png] view at source ↗
Figure 2
Figure 2. Water dimer energy curve (structure shown) for O–O distances ranging from 2.3 [PITH_FULL_IMAGE:figures/full_fig_p012_2.png] view at source ↗
Figure 3
Figure 3. Polarization density, defined as ρ(r) − ρiso(r), where ρiso is the sum of the elec￾tron densities of the isolated water monomers, for the water dimer at the equilibrium O–O distance. In each panel, we present isosurfaces (top) and contour plots (bottom) generated with a cutoff of ±0.0007 e · a −3 0 . The key aspect of this system lies in its asymmetry, as one monomer is a hydrogen bond donor, and the other is an acc… view at source ↗
Figures from the paper (3 more)
Figure 4
Figure 4. Figure 4: Panels (a), (b) and (c): correlation plots of the interaction energy (in kcal/mol) of [PITH_FULL_IMAGE:figures/full_fig_p019_4.png]
Figure 5
Figure 5. Figure 5: Polarization density (defined as the difference of the embedded and isolated [PITH_FULL_IMAGE:figures/full_fig_p020_5.png]
Figure 6
Figure 6. Figure 6: Convergence of the glucose (left panel) and the [Pd(H [PITH_FULL_IMAGE:figures/full_fig_p023_6.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

104 extracted references · 79 canonical work pages

  1. [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

  2. [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

  3. [3]

    M.; Thiel, W

    Senn, H. M.; Thiel, W. QM/MM studies of enzymes. Current Opinion in Chemical Biology 2007, 11, 182--187

  4. [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

  5. [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

  6. [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

  7. [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

  8. [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

Show all 104 references
  1. [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

  2. [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

  3. [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

  4. [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

  5. [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

  6. [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

  7. [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

  8. [16]

    Electrostatic embedding of machine learning potentials

    Zinovjev, K. Electrostatic embedding of machine learning potentials. J. Chem. Theory Comput. 2023, 19, 1888--1897

  9. [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,

  10. [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

  11. [19]

    Introduction to QM/MM simulations

    Groenhof, G. Introduction to QM/MM simulations. Biomolecular simulations: methods and protocols 2013, 43--66

  12. [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

  13. [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

  14. [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

  15. [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

  16. [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

  17. [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

  18. [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

  19. [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

  20. [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...

  21. [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

  22. [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

  23. [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

  24. [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

  25. [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

  26. [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

  27. [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

  28. [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

  29. [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

  30. [38]

    Karelina, M.; Kulik, H. J. Systematic quantum mechanical region determination in QM/MM simulation. J. Chem. Theory Comput. 2017, 13, 563--576

  31. [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

  32. [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

  33. [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

  34. [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

  35. [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

  36. [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

  37. [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

  38. [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

  39. [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

  40. [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

  41. [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

  42. [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

  43. [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

  44. [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

  45. [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

  46. [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

  47. [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

  48. [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

  49. [57]

    R.; Neugebauer, J

    Jacob, C. R.; Neugebauer, J. Subsystem density-functional theory. WIREs: Comp. Mol. Sci. 2014, 4, 325--362

  50. [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

  51. [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

  52. [60]

    R.; Neugebauer, J

    Jacob, C. R.; Neugebauer, J. Subsystem density-functional theory (update). Wiley Interdisciplinary Reviews: Computational Molecular Science 2024, 14, e1700

  53. [61]

    See supplementary materials document at [url-to-be-inserted]

    Supplementary Materials. See supplementary materials document at [url-to-be-inserted]

  54. [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

  55. [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

  56. [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

  57. [65]

    Nonlocal subsystem density functional theory

    Mi, W.; Pavanello, M. Nonlocal subsystem density functional theory. J. Phys. Chem. Lett. 2019, 11, 272--279

  58. [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

  59. [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

  60. [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

  61. [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

  62. [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

  63. [71]

    Focke, K.; Jacob, C. R. Coupled-cluster density-based many-body expansion. The Journal of Physical Chemistry A 2023, 127, 9139--9148

  64. [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

  65. [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

  66. [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

  67. [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

  68. [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

  69. [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

  70. [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

  71. [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)

  72. [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)

  73. [81]

    Giannozzi, P. et al. Advanced capabilities for materials modelling with QUANTUM ESPRESSO. Journal of Physics: Condensed Matter 2017, 29, 465901

  74. [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

  75. [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

  76. [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

  77. [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

  78. [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

  79. [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

  80. [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

  81. [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...

  82. [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

  83. [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

  84. [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

  85. [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

  86. [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

  87. [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

  88. [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

  89. [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

  90. [98]

    Nature 1905, 72, 318--319

    Lord Rayleigh , The problem of the random walk. Nature 1905, 72, 318--319

  91. [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

  92. [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

  93. [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

  94. [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

  95. [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

  96. [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...

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.