Simultaneous Learning of Static and Dynamic Charges
Pith reviewed 2026-05-21 16:45 UTC · model grok-4.3
The pith
Despite their physical link, learning static and dynamic charges independently is the more practical approach for modeling condensed-phase and cluster systems.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In comparisons on water systems, coupled learning with a learned environment-dependent screening factor restores accuracy for dynamical charges but provides no significant accuracy gain over independent models while increasing computational expense, indicating that independent modeling of static and dynamic charges is more practical despite their formal connection.
What carries the argument
The comparison of independent versus coupled learning of static charges and Born effective charges, using either a single global coupling constant or a local environment-dependent screening factor to account for dielectric screening.
If this is right
- Independent models achieve comparable accuracy for both bulk water and water clusters without added complexity.
- Coupled models require correction for dielectric screening to be effective, yet the common assumption of homogeneous isotropic screening fails in heterogeneous clusters.
- Environment-dependent screening improves coupled predictions for dynamical charges but still provides negligible benefit over separate models.
- Separate models for static and dynamic charges are the recommended practical choice for applications in both condensed-phase and isolated cluster systems.
Where Pith is reading between the lines
- The finding that independent modeling is preferable may extend to other polar or heterogeneous materials where local screening varies strongly.
- Future tests on non-aqueous systems could confirm whether the negligible gain from coupling is general or specific to water-like hydrogen-bonded networks.
- Developers of machine-learned potentials may simplify their architectures by maintaining separate predictors for different charge types rather than enforcing physical couplings.
Load-bearing premise
Water clusters and bulk water are sufficiently representative of the heterogeneous condensed-phase systems where these charge models would be applied.
What would settle it
Repeating the accuracy and cost comparison on a different heterogeneous system such as a solvated biomolecule or an ionic liquid interface and observing whether coupled learning then yields clear accuracy gains that outweigh the added cost.
Figures
read the original abstract
Long-range interactions and electric response are essential for accurate modeling of condensed-phase systems, but capturing them efficiently remains a challenge for atomistic machine learning. Traditionally, these two phenomena can be represented by static charges, that participate in Coulomb interactions between atoms, and dynamic charges such as atomic polar tensors - aka Born effective charges - describing the response to an external electric field. We critically compare different approaches to learn both types of charges, taking bulk water and water clusters as paradigmatic examples: (1) Learning them independently; (2) Coupling static and dynamic charges based on their physical relationship with a single global coupling constant to account for dielectric screening; (3) Coupled learning with a local, environment-dependent screening factor. In the coupled case, correcting for dielectric screening is essential, yet the common assumption of homogeneous, isotropic screening breaks down in heterogeneous systems such as water clusters. A learned, environment-dependent screening restores high accuracy for the dynamical charges. However, the accuracy gain over independent dynamic predictions is negligible, while the computational cost increases compared to using separate models for static and dynamical charges. This suggests that, despite the formal connection between the two charge types, modeling them independently is the more practical choice for both condensed-phase and isolated cluster systems.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper compares three strategies for learning static atomic charges and dynamic charges (Born effective charges) in atomistic machine learning: independent training of separate models, coupled training with a single global coupling constant to account for dielectric screening, and coupled training with a learned local environment-dependent screening factor. Using bulk water (homogeneous) and water clusters (heterogeneous) as paradigmatic systems, the authors show that homogeneous screening breaks down in clusters, local screening recovers accuracy for dynamic charges, but the accuracy gain over independent learning remains negligible while increasing computational cost. They conclude that independent modeling is the more practical choice for both condensed-phase and isolated cluster systems despite the formal physical connection between charge types.
Significance. If the empirical comparisons hold, the work offers practical guidance for efficient modeling of long-range electrostatics and electric response in machine learning potentials for chemistry. Highlighting the failure of homogeneous screening assumptions in heterogeneous environments and showing limited benefit from coupling can streamline model development and reduce training overhead. The focus on water as a testbed provides clear illustrations of the tradeoffs, and the empirical comparison of training strategies adds to the literature on charge learning in condensed-phase systems.
major comments (2)
- [Conclusions] Conclusions and discussion of practicality: the central recommendation that independent modeling is more practical for condensed-phase and cluster systems because coupled approaches yield 'negligible' accuracy gains rests on the representativeness of bulk water and water clusters; in other heterogeneous condensed-phase environments (e.g., liquid-solid interfaces or solvated biomolecules) the environment dependence could be stronger and the error reduction from local screening larger, altering the cost-accuracy tradeoff. This assumption underpins the generalization and requires either additional tests or a more qualified statement.
- [Results] Results section on water clusters: the statement that local screening 'restores high accuracy' yet the gain over independent dynamic predictions 'is negligible' is load-bearing for the practicality claim; without explicit numerical values (e.g., MAE or RMSE for Born charges, with error bars or statistical tests) and direct cost comparisons in the relevant table or figure, the magnitude of the tradeoff cannot be independently assessed.
minor comments (2)
- [Methods] Methods: specify whether the same ML architecture and hyperparameters are used across the three training strategies to ensure the comparison isolates the effect of coupling versus independence.
- [Theory] Notation: define the local screening factor with an explicit equation or formula when first introduced to improve clarity for readers unfamiliar with the coupling approach.
Simulated Author's Rebuttal
We thank the referee for the constructive and detailed report. We address each major comment below and have revised the manuscript accordingly to strengthen the presentation and qualify our conclusions.
read point-by-point responses
-
Referee: [Conclusions] Conclusions and discussion of practicality: the central recommendation that independent modeling is more practical for condensed-phase and cluster systems because coupled approaches yield 'negligible' accuracy gains rests on the representativeness of bulk water and water clusters; in other heterogeneous condensed-phase environments (e.g., liquid-solid interfaces or solvated biomolecules) the environment dependence could be stronger and the error reduction from local screening larger, altering the cost-accuracy tradeoff. This assumption underpins the generalization and requires either additional tests or a more qualified statement.
Authors: We agree that the generalization of our practicality recommendation is limited by the choice of paradigmatic systems. Bulk water and water clusters were selected to contrast homogeneous and heterogeneous environments, and the results clearly show the breakdown of global screening in the latter. We have revised the conclusions section to include a more qualified statement acknowledging that stronger environment dependence in other systems (such as liquid-solid interfaces or solvated biomolecules) could alter the observed accuracy-cost tradeoff. Additional tests in those systems are beyond the scope of the present work. revision: partial
-
Referee: [Results] Results section on water clusters: the statement that local screening 'restores high accuracy' yet the gain over independent dynamic predictions 'is negligible' is load-bearing for the practicality claim; without explicit numerical values (e.g., MAE or RMSE for Born charges, with error bars or statistical tests) and direct cost comparisons in the relevant table or figure, the magnitude of the tradeoff cannot be independently assessed.
Authors: We thank the referee for this observation. In the revised manuscript we have added explicit MAE and RMSE values (with error bars from five independent training runs using different random seeds) for the Born effective charges on water clusters in the results section and the associated table. We have also included a direct comparison of training and inference wall times for the independent versus coupled models, reported in the supplementary information with a reference in the main text. These additions make the magnitude of the accuracy gain and computational overhead directly quantifiable. revision: yes
Circularity Check
No circularity: empirical comparison of independent vs coupled charge models
full rationale
The paper conducts an empirical benchmark of three training strategies (independent, global-coupling, local-screening) on bulk water and water-cluster data. All reported accuracy gains, cost differences, and the final practicality conclusion are obtained by direct numerical evaluation of trained models against reference charges and forces; no central result is obtained by substituting a fitted parameter back into itself or by reducing a prediction to a quantity defined solely by the training procedure. Self-citations, if present, are not load-bearing for the comparison itself.
Axiom & Free-Parameter Ledger
free parameters (2)
- global coupling constant
- local environment-dependent screening factor
axioms (2)
- domain assumption Static and dynamic charges are physically related through dielectric screening.
- domain assumption Homogeneous isotropic screening holds in bulk but breaks down in heterogeneous clusters.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We critically compare different approaches to learn both types of charges... (1) Learning them independently; (2) Coupling static and dynamic charges based on their physical relationship with a single global coupling constant... (3) Coupled learning with a local, environment-dependent screening factor.
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
the common assumption of homogeneous, isotropic screening breaks down in heterogeneous systems such as water clusters
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Forward citations
Cited by 1 Pith paper
-
Mimyria: Machine learned vibrational spectroscopy for aqueous systems made simple
Mimyria learns atom-resolved atomic polar tensors and a new polarizability gradient tensor to produce accurate IR and Raman spectra from MD trajectories using small training sets from ab initio data.
Reference graph
Works this paper leans on
-
[1]
describes the charge flux: the change in the dipole moment due to a collective alter- ation of the charge distribution as the i-th atom is dis- placed by an infinitesimal amount. Because Z i encodes these charge redistribution effects, BECs are also called “dynamic charges”. B. Long-Ranged Machine Learned Potentials Via Learned Static Charges As briefly expla...
-
[2]
Learning Long-Range Representations with Equivariant Messages
Here, the electrostatic potential is processed by a learnable function fθ , yielding U lr = ∑ i fθ (Vi({qj}j, ξ i) . (6) This modification, referred to as “Learning Long-Range Representations with Equivariant Messages” (LOREM) has been shown to significantly reduce errors across mul- tiple prediction tasks. We note that in this work we re- strict the archit...
-
[3]
( 5) and retains interpretability of the pseudo charges
When the learnable function is set to be the identity, fθ = id, the model re- duces to eq. ( 5) and retains interpretability of the pseudo charges. We refer to this latter version as the physical long-range architecture. The simplest extension of long-ranged MLIPs to also predict BECs is to treat the BECs as independent out- puts without connection to the...
-
[4]
( 3) is the basis for the second family of MLIPs that we study, the coupled models
and ( 8) to predict Z i via eq. ( 3) is the basis for the second family of MLIPs that we study, the coupled models. These mod- els, shown in fig. 1B, use the predicted qi(ξi) to construct the BECs. We explore two variants of these coupled mod- els. (i) In the coupled, global approach, the MLIP is first trained only on energies and forces. Afterwards, a sin-...
-
[5]
and ( 7) can −1 0 1 Z pred α,α [e] A Uncoupl. Coupl. γ Coupl. γi −1.5 −1.0 −0.5 0.0 0.5 1.0 1.5 Z DFT α,α [e] −1 0 1 Z pred α,α [e] B FIG. 2: Scatter plot of diagonal components of the BECs Zα,α comparing DFT values with model predictions for bulk water (A) and water clusters (B). result in ill-defined values, due to the periodic bound- aries. Zhong et al....
-
[7]
may perform even better, we also trained models with the LOREM formulation. As can be seen in fig. 4, this modification improves performance on all targets with a particularly strong improvement in bulk potential energy errors. The parity plot (similar to fig. 2) for the LOREM models are presented in the supporting information (Figure S2), showing qualitativ...
-
[8]
The water cluster structures are composed from three sources: (1) the water cluster subset of Hobza’s benchmark en- ergy and geometry data base (BEGDB)[ 54], (2) the W A- TER27 component of the GMTKN24 and GMTKN30 benchmark suites[ 55] and (3) dimer and trimer structures from Ref. 56. These sets contain very distinct geometri- cal structures from 1 to 20 ...
work page 2000
-
[9]
Provided P is known for a reference configuration (e.g
with respect to ri is thus given by: Fi,α = − ∂ ∂r i,α ∑ j U j model + ∑ β Z i α,β Eext β , (12) neglecting second order derivatives of Z i with respect to ri. Provided P is known for a reference configuration (e.g. via an ab-initio calculation for the starting struc- ture), the enthalpy H is well-defined through Pref during a simulation run. This p...
-
[10]
0 ˚ A and spherical harmonics up to degree 6, with eight spherical channels
All models use a cutoff radius of rc = 5. 0 ˚ A and spherical harmonics up to degree 6, with eight spherical channels. The radial dependence is repre- sented using 32 radial basis functions, and the chemical embedding uses eight channels. We employ a cosine cut- off to ensure smooth behavior at the cutoff distance, and perform the radial expansion using basi...
-
[11]
for non-periodic structures and eq. ( 9) for periodic structures) we deter- mine a single γ value which minimizes the mean average error on BEC predictions on the entire validation set, in- cluding clusters and bulk structures simultaneously. This value is then used for all subsequent evaluations and sim- ulations of the model. C. Molecular Dynamics and I...
work page 2075
-
[12]
Y. Pan, X. Wang, W. Zhang, L. Tang, Z. Mu, C. Liu, B. Tian, M. Fei, Y. Sun, H. Su, et al., Nature Communi- cations 13, 3063 (2022), ISSN 2041-1723
work page 2022
- [13]
- [14]
-
[15]
O. T. Unke, S. Chmiela, H. E. Sauceda, M. Gastegger, I. Poltavsky, K. T. Sch¨ utt, A. Tkatchenko, and K.-R. M¨ uller, Chem. Rev.121, 10142 (2021), ISSN 0009-2665
work page 2021
-
[16]
K. Sch¨ utt, P.-J. Kindermans, H. E. Sauceda Felix, S. Chmiela, A. Tkatchenko, and K.-R. M¨ uller, Advances in neural information processing systems 30, 0 (2017)
work page 2017
-
[17]
Tensor field networks: Rotation- and translation-equivariant neural networks for 3D point clouds
N. Thomas, T. Smidt, S. Kearnes, L. Yang, L. Li, K. Kohlhoff, and P. Riley, arXiv preprint arXiv:1802.08219 (2018)
work page internal anchor Pith review Pith/arXiv arXiv 2018
-
[18]
S. Batzner, A. Musaelian, L. Sun, M. Geiger, J. P. Mailoa, M. Kornbluth, N. Molinari, T. E. Smidt, and B. Kozinsky, Nature communications 13, 2453 (2022)
work page 2022
-
[19]
I. Batatia, D. P. Kovacs, G. Simm, C. Ortner, and G. Cs´ anyi, Advances in neural information processing systems 35, 11423 (2022)
work page 2022
-
[20]
S. Falletta, A. Cepellotti, A. Johansson, C. W. Tan, A. Musaelian, C. J. Owen, and B. Kozinsky, Uni- fied Differentiable Learning of Electric Response (2024), 2403.17207
- [21]
-
[22]
M. Veit, D. M. Wilkins, Y. Yang, R. A. DiStasio, Jr., and M. Ceriotti, The Journal of Chemical Physics 153, 024113 (2020), ISSN 0021-9606
work page 2020
-
[23]
R. Resta and D. Vanderbilt, in Physics of Ferroelectrics (Springer Berlin Heidelberg, Berlin, Heidelberg, 2007), vol. 105, pp. 31–68, ISBN 978-3-540-34590-9
work page 2007
- [24]
- [25]
-
[26]
J. Li, L. Knijff, Z.-Y. Zhang, L. Andersson, and C. Zhang, Journal of Chemical Theory and Computation 21, 1382 (2025), ISSN 1549-9618
work page 2025
-
[27]
Schienbein, Journal of Chemical Theory and Compu- tation 19, 705 (2023), ISSN 1549-9618
P. Schienbein, Journal of Chemical Theory and Compu- tation 19, 705 (2023), ISSN 1549-9618
work page 2023
-
[28]
K. Joll, P. Schienbein, K. M. Rosso, and J. Blumberger, Nature Communications 15, 8192 (2024), ISSN 2041- 1723
work page 2024
-
[29]
N. Bergmann, N. Bonnet, N. Marzari, K. Reuter, and N. G. H¨ ormann, Phys. Rev. Lett. 135, 146201 (2025)
work page 2025
-
[30]
A. P. Bart´ ok, M. C. Payne, R. Kondor, and G. Cs´ anyi, Physical Review Letters 104, 136403 (2010), ISSN 0031- 9007, 1079-7114
work page 2010
-
[31]
Z. Deng, C. Chen, X.-G. Li, and S. P. Ong, npj Compu- tational Materials 5, 1 (2019), ISSN 2057-3960
work page 2019
-
[32]
N. Artrith, T. Morawietz, and J. Behler, Physical Review B 83, 153101 (2011)
work page 2011
-
[33]
T. Morawietz, V. Sharma, and J. Behler, The Journal of Chemical Physics 136, 064103 (2012), ISSN 0021-9606
work page 2012
-
[34]
K. Yao, J. E. Herr, D. W. Toth, R. Mckintyre, and J. Parkhill, Chemical Science 9, 2261 (2018), ISSN 2041- 6539
work page 2018
-
[35]
T. W. Ko, J. A. Finkler, S. Goedecker, and J. Behler, Nature Communications 12, 398 (2021), ISSN 2041-1723
work page 2021
- [36]
-
[37]
Cheng, npj Computational Materials 11, 1 (2025), ISSN 2057-3960
B. Cheng, npj Computational Materials 11, 1 (2025), ISSN 2057-3960
work page 2025
- [38]
-
[39]
K. C. Gross, P. G. Seybold, and C. M. Hadad, Inter- national Journal of Quantum Chemistry 90, 445 (2002), ISSN 1097-461X
work page 2002
-
[40]
A. E. Reed, R. B. Weinstock, and F. Weinhold, The Jour- nal of Chemical Physics 83, 735 (1985), ISSN 0021-9606
work page 1985
-
[41]
F. L. Hirshfeld, Theoretica chimica acta 44, 129 (1977), ISSN 1432-2234
work page 1977
- [43]
-
[44]
A. Milani and C. Castiglioni, Journal of Molecular Struc- ture: THEOCHEM 955, 158 (2010), ISSN 0166-1280
work page 2010
-
[45]
R. S. Mulliken, The Journal of Chemical Physics 23, 1833 (1955), ISSN 0021-9606
work page 1955
- [46]
-
[47]
L.-P. Wang, T. J. Martinez, and V. S. Pande, J. Phys. Chem. Lett. 5, 1885 (2014)
work page 2014
-
[48]
K. K. Huguenin-Dumittan, P. Loche, N. Haoran, and M. Ceriotti, J. Phys. Chem. Lett. pp. 9612–9618 (2023)
work page 2023
-
[50]
E. Rumiantsev, M. F. Langer, T.-E. Sodjargal, M. Ceri- otti, and P. Loche, Learning Long-Range Representations with Equivariant Messages (2025), 2507.19382
- [51]
-
[52]
B. Schmiedmayer and G. Kresse, The Journal of Chemi- cal Physics 161, 084703 (2024), ISSN 0021-9606
work page 2024
-
[53]
I. V. Leontyev and A. A. Stuchebrukhov, The Journal of Chemical Physics 130, 085102 (2009), ISSN 0021-9606
work page 2009
-
[54]
I. V. Leontyev and A. A. Stuchebrukhov, Journal of Chemical Theory and Computation 6, 3153 (2010), ISSN 1549-9618
work page 2010
-
[55]
A. A. Kornyshev, Journal of electroanalytical chemistry and interfacial electrochemistry 204, 79 (1986)
work page 1986
-
[56]
H. A. Stern and S. E. Feller, The Journal of Chemical Physics 118, 3401 (2003), ISSN 0021-9606
work page 2003
-
[57]
V. Ballenegger and J.-P. Hansen, The Journal of Chem- ical Physics 122, 114711 (2005), ISSN 0021-9606
work page 2005
-
[58]
P. St¨ ark, H. Stooß, P. Loche, D. J. Bonthuis, R. R. Netz, and A. Schlaich, Static Dielectric Permittivity Profiles and Coarse-graining Approaches for Water in Graphene Slit Pores (2025), 2512.07548
work page internal anchor Pith review Pith/arXiv arXiv 2025
-
[59]
M. F. Langer, J. T. Frank, and F. Knoop, The Journal of Chemical Physics 159, 174105 (2023), ISSN 0021-9606
work page 2023
-
[62]
C. P´ erez, M. T. Muckle, D. P. Zaleski, N. A. Seifert, B. Temelso, G. C. Shields, Z. Kisiel, and B. H. Pate, Science 336, 897 (2012)
work page 2012
-
[63]
Y. Wang and J. M. Bowman, J. Phys. Chem. Lett. 4, 1104 (2013), ISSN 1948-7185, 1948-7185
work page 2013
- [64]
-
[65]
B. Temelso, K. A. Archer, and G. C. Shields, The Journal of Physical Chemistry A 115, 12034 (2011), ISSN 1089- 5639
work page 2011
-
[66]
V. S. Bryantsev, M. S. Diallo, A. C. T. van Duin, and W. A. I. Goddard, Journal of Chemical Theory and Com- putation 5, 1016 (2009), ISSN 1549-9618
work page 2009
-
[67]
T. T. Nguyen, E. Sz´ ekely, G. Imbalzano, J. Behler, G. Cs´ anyi, M. Ceriotti, A. W. G¨ otz, and F. Paesani, The Journal of Chemical Physics 148, 241725 (2018), ISSN 0021-9606
work page 2018
-
[68]
A. Mazitov, F. Bigi, M. Kellner, P. Pegolo, D. Tisi, G. Fraux, S. Pozdnyakov, P. Loche, and M. Ceriotti, PET-MAD, a universal interatomic potential for ad- vanced materials modeling (2025), 2503.14118
-
[69]
A. Mazitov, S. Chorna, G. Fraux, M. Bercx, G. Pizzi, S. De, and M. Ceriotti, Scientific Data 12, 1857 (2025), ISSN 2052-4463
work page 2025
-
[70]
T. D. K¨ uhne, M. Iannuzzi, M. Del Ben, V. V. Ry- bkin, P. Seewald, F. Stein, T. Laino, R. Z. Khaliullin, O. Sch¨ utt, F. Schiffmann, et al., The Journal of Chemi- cal Physics 152, 194103 (2020), ISSN 0021-9606
work page 2020
- [71]
- [72]
-
[73]
M. Stengel, N. A. Spaldin, and D. Vanderbilt, Nature Physics 5, 304 (2009), ISSN 1745-2481
work page 2009
-
[74]
S. A. Barr and A. Z. Panagiotopoulos, Physical Review E 86, 016703 (2012)
work page 2012
-
[75]
G. G. Lorentz, Bernstein Polynomials (American Mathe- matical Society, New York, NY, 1986), ISBN 978-0-8284- 0323-8
work page 1986
-
[76]
A. Hjorth Larsen, J. Jørgen Mortensen, J. Blomqvist, I. E. Castelli, R. Christensen, M. Du/suppress lak, J. Friis, M. N. Groves, B. Hammer, C. Hargus, et al., Journal of Physics: Condensed Matter 29, 273002 (2017), ISSN 0953-8984
work page 2017
-
[77]
G. Bussi, D. Donadio, and M. Parrinello, The Journal of Chemical Physics 126, 014101 (2007), ISSN 0021-9606. Supplemental Information: Simultaneous Learning of Static and Dynamic Charg es Philipp St¨ ark,1, 2 Henrik Stooß, 3 Marcel F. Langer, 4 Egor Rumiantsev, 4 Alexander Schlaich, 3, ∗ Michele Ceriotti, 4, † and Philip Loche 4, ‡ 1Stuttgart Center for S...
- [78]
-
[79]
M. F. Langer, J. T. Frank, and F. Knoop, The Journal of Chemica l Physics 159, 174105 (2023), ISSN 0021-9606. 5 10−8 10−6 10−4 10−2 100 Eext [V/ ˚A] −1.0 −0.5 0.0 0.5 T rZ i T rZ Oxygen T rZ Hydrogen FIG. S5: Stability of finite differences estimates for the BECs a s a function of applied fields Eext in a finite differences scheme. The boxplots show median, fir...
work page 2023
-
[80]
M. F. Langer, Physical Review B 108, 0 (2023)
work page 2023
-
[81]
S. Carlson, F. N. Br¨ unig, P. Loche, D. J. Bonthuis, and R. R. Netz, The Journal of Physical Chemistry A 124, 5599 (2020), ISSN 1089-5639
work page 2020
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.