Pith. sign in

REVIEW 3 major objections 5 minor 54 references

A Simple Iterative Approach for Constant Chemical Potential Simulations at Interfaces

T0 review · 3 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read An iterative run-measure-update loop that rescales ion counts by the ratio of target to measured bulk concentration reaches a stable 1 M bulk concentration in two iterations in all three interfacial systems tested, including a…

desk verdict A practical constant-concentration MD loop, overbranded as constant chemical potential; the μ claim is untested and the core value is the simple workflow plus MLIP integration. read the letter →

arxiv 2506.01050 v1 pith:LF5JVGRS submitted 2025-06-01 physics.chem-ph cond-mat.mtrl-sci

Ademola Soyemi , Khagendra Baral , Tibor Szilvasi This is my paper
classification physics.chem-phcond-mat.mtrl-sci
keywords iterativeconstantchemicalpotentialmoleculardynamicssolid-liquidinterfaceliquid-airmachinelearninginteratomicbulkconcentrationcontrolelectrolyte
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

Conventional molecular dynamics fixes the number of ions, so ions can migrate to an interface and deplete the bulk; the paper introduces iterative constant chemical potential molecular dynamics (iCμMD), which periodically measures the bulk ion concentration from a short MD run, then adds or removes ions to steer that concentration to a preset target. The claim is that a simple run-measure-update loop can reproduce a constant-chemical-potential condition in just two iterations, without a reservoir region or an on-the-fly particle-insertion scheme. The paper demonstrates this on NaCl(aq) against graphite and against air using classical force fields, and on Na2SO4(aq) between graphene sheets using a machine-learned interatomic potential trained on DFT data. If correct, the method gives a cheap way to hold a chosen bulk salt concentration in interfacial simulations and brings constant-potential-style boundary conditions to MLIP-driven MD.

What carries the argument

The load-bearing mechanism is Eq. (1), a linear scaling of the total ion count by the ratio of target to measured bulk concentration, $N_t^{i+1} = (C_b^t / C_b^i) N_t^i$. The bulk region is defined as a slab away from the interface (typically ±5 or ±7 Å around the cell center), and an iteration is considered converged when the slope of a 5 ns moving average of bulk concentration stays within $\pm 2 \times 10^{-3}$ M/ns for at least 15 ns. That converged measurement feeds the update rule, and the whole loop repeats until the bulk concentration is within ±0.03 M of the target.

What would settle it

Run iCμMD on a 1 M or higher electrolyte where the activity coefficient is known to deviate strongly from unity, and independently compute the chemical potential of the ion (for example, by thermodynamic integration or osmotic-pressure measurement) at the initial and converged ion counts; if the converged bulk concentration corresponds to a chemical potential measurably different from the target, the constant-μ claim fails. A simpler check is to compare the bulk concentration reached when the same target is approached from above versus from below: if the final states differ beyond the ±0.03 M window, the linear update rule is not converging to a unique chemical potential.

Watch

Extended reading notes

Core claim

The central claim is that a target bulk ion concentration can be reached by iterating a proportional update: after each converged MD run, the total number of target ions is scaled by the ratio of target to measured bulk concentration, $N_t^{i+1} = (C_b^t / C_b^i) N_t^i$ (Eq. 1). In the three test systems the first iteration depleted or enriched the bulk relative to the nominal 1 M (0.85 M for NaCl–graphite, 1.07 M for NaCl–air, 1.07 M for Na2SO4–graphene), one proportional update corrected the count, and the second iteration stabilized within ±0.03 M of 1 M for hundreds of nanoseconds. The paper takes this as evidence that iCμMD achieves a constant-μ condition in two iterations and that its only requirement is a measurable bulk region.

Load-bearing premise

The method's control target is a bulk molar concentration, and its correction is a straight proportion, so the whole procedure stands on the assumption that bulk concentration is a faithful proxy for chemical potential — that the solution behaves ideally with a fixed activity coefficient across the ion-count changes.

Editorial extensions

If this is right

  • In the three demonstrator systems, the method reaches a stable 1 M bulk concentration in two iterations, so a user can expect convergence after roughly 50–100 ns of MD rather than an on-the-fly insertion scheme.
  • Because the loop only requires measuring a bulk concentration and restarting the simulation with a new ion count, it drops into existing MD workflows that already use classical or MLIP potentials.
  • The paper states that the framework generalizes to solid-gas, liquid-liquid, and mixed-solute systems, provided the interfacial simulation can be converged with a measurable bulk region.
  • For MLIP-driven simulations, iCμMD supplies a route to constant-μ boundary conditions with DFT-level accuracy, which matters for double-layer and electrochemical studies.

Reading between the lines

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

  • The paper never computes the chemical potential itself, so a natural next test is to convert the converged 1 M bulk state into an excess chemical potential via free-energy perturbation and check whether the iCμMD operating point equals a thermodynamic target.
  • Because Eq. (1) is a direct proportion, the method should be easiest to validate at low concentration; at 1 M and above it may need an activity-coefficient correction or a concentration-dependent proportionality factor.
  • The same loop could be adapted to control multiple species simultaneously, for example pH or ionic strength, by updating each species count with its own measured bulk concentration, though the paper's demonstrations treat a single electrolyte.
  • The stopping criterion relies on a user-chosen slope tolerance, so porting iCμMD to a new system will require a short convergence check on the moving-average window before trusting the two-iteration claim.
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

3 major / 5 minor

Summary. The manuscript introduces iCμMD, an iterative 'run-measure-update' scheme in which the total number of solute species in a molecular dynamics simulation is adjusted based on the measured bulk concentration, with the aim of achieving a target bulk molar concentration. The method is demonstrated on three interfacial systems: NaCl(aq) on graphite, NaCl(aq) at the air interface, and Na2SO4(aq) on graphene, the last using a machine-learned interatomic potential. The authors report that two iterations suffice to reach a target bulk concentration of 1 M within ±0.03 M, and they claim that this procedure realizes a constant chemical potential condition.

Significance. The proposed algorithm is simple, requires no reservoir region, and is compatible with both classical force fields and MLIPs, which is a practical advantage. The authors provide reproducible training data, pseudocode, and open-source Zenodo deposition for the MLIP, and they demonstrate a genuine MLIP integration for a polyatomic divalent electrolyte. However, the central advertised contribution—'constant chemical potential simulations'—is not supported by the evidence presented. What is demonstrated is iterative control of bulk concentration; the equivalence between a target molarity and a fixed chemical potential is assumed rather than tested. If the authors reframe the contribution accordingly and add a thermal-validation check, the method would still be useful as a concentration-control tool, but in its current form the significance is substantially weaker than claimed.

major comments (3)
  1. [Section 2; Conclusions] The manuscript claims that achieving the target bulk concentration realizes a 'constant chemical potential condition' and a 'constant-mu condition', but no thermodynamic quantity relating concentration to chemical potential is ever computed. The paper does not calculate the chemical potential, an activity coefficient, or a Kirkwood-Buff integral, nor does it compare against CμMD or GCMC on the same system. For a 1 M electrolyte such as NaCl, the mean ionic activity coefficient is considerably below unity and concentration-dependent, so the equivalence between a target molarity and a fixed chemical potential requires an ideal-solution assumption that is neither stated nor tested. The abstract and conclusion therefore overstate the result; the established contribution is iterative bulk-concentration control.
  2. [Equation (1); Sections 4.1 and 4.2] The update rule N_t^{i+1} = (C_b^t / C_b^i) N_t^i presumes a strictly linear relation between the total number of target species and the bulk concentration. This ignores nonlinear partitioning between interface and bulk, such as saturation of interfacial adsorption sites or concentration-dependent interfacial affinity. All three demonstrations start from a nominal 1 M solution and make only small corrections (0.85→1 M, 1.07→1 M, and 1.07→1 M), so the claimed 'two-iteration' convergence is only shown for starting points already near the target. A test from a distant starting concentration (e.g., 0.1 M or 3 M) would be needed to support the general convergence claim.
  3. [Section 4; Concluding remarks] The validation of the method is circular: the controlled variable is the bulk concentration, and the convergence criterion is agreement with the target bulk concentration. Since the algorithm adjusts particle counts to force that agreement, the loop converges by construction. No independent observable is used to confirm that the resulting state corresponds to a fixed chemical potential, nor is a direct comparison with the external benchmark CμMD (which is discussed in the introduction) presented for the same systems. A comparison of interfacial concentration profiles or of a computed excess chemical potential would break this circularity.
minor comments (5)
  1. [Abstract] The abstract uses 'iCuMD' while the main text uses 'iCμMD'; please standardize the notation to one form throughout.
  2. [References] References 8 and 10 cite the identical paper (Finney et al., Chemical Science, 2021, 12, 11166–11180); this duplicate should be removed and the numbering adjusted.
  3. [Section 4.1] In the paragraph discussing the NaCl(aq)-air system, the text refers to 'Figure 3 (right)' when describing the Z-direction ion distribution; this should be 'Figure 4 (right)' for consistency with the figure numbering.
  4. [Section 3.2] The sentence 'the and with electronic occupations smeared...' contains a typo; it should read 'and the electronic occupations are smeared...'.
  5. [Conclusions] The phrase 'in the the NaClaq-graphite' contains a doubled article and should be corrected.

Circularity Check

2 steps flagged · score 6.0 of 10

iCμMD's constant-chemical-potential claim reduces to its own bulk-concentration control target.

  1. self definitional [Abstract; Sec. 2 (The iCμMD Approach); Sec. 3.1/3.3 stopping criteria; Conclusions]
    "adjusting the number of species in the solution to reach a target concentration (chemical potential) ... By iterating this ‘run-measure-update’ loop until the bulk concentration converges to the desired target concentration, iCμMD can achieve a constant‑μ condition in a few iterations ... we stop iterating when the bulk concentration is within ±0.03 M of the target concentration."

    The control objective and the validation metric are the same quantity: the target bulk concentration. The abstract explicitly equates the target concentration with chemical potential, and the conclusion calls convergence of the bulk concentration a 'constant-μ condition.' Because no chemical potential, activity coefficient, or independent thermodynamic observable is computed, the central advertised claim is a restatement of the loop's own stopping rule rather than an independent result.

  2. fitted input called prediction [Eq. 1 (Sec. 2.2) and Sec. 4.1 (NaClaq-graphite result)]
    "𝑁𝑡𝑖+1 = (𝐶𝑏𝑡/𝐶𝑏𝑖) 𝑁𝑡𝑖 ... Using equation 1, we then calculated the required number of ions to reach the target concentration and added 39 more ion pairs ... resulting in a bulk concentration of 0.99 M which remains stable (within ±0.03 M of the targeted 1 M bulk concentration)."

    The update rule rescales the total ion count by the ratio of target to measured bulk concentration, which implicitly assumes C_b scales linearly with N_t. Under that assumption the next measured bulk concentration is forced toward the target. Reporting that the second iteration lands within ±0.03 M is therefore a consequence of the update/stopping construction, not an independent prediction. No CμMD, GCMC, or chemical-potential comparison is run to break the self-reference.

full rationale

The paper's core run-measure-update loop is internally consistent as a concentration-control scheme: Eq. 1 rescales the ion count by the ratio of target to measured bulk concentration, and the stopping rule terminates when the bulk concentration is within ±0.03 M. The circularity is in the advertised physical interpretation: the abstract identifies chemical potential with the target concentration ('target concentration (chemical potential)'), and the conclusion calls the same bulk-concentration convergence a 'constant-μ condition.' No activity coefficient, excess chemical potential, Kirkwood-Buff integral, or independent CμMD/GCMC calculation is reported, so the central constant-chemical-potential claim is not independently established; it is a renaming of the controlled variable. The 'two iterations' demonstrations are also close to trivial: all start at a nominal 1 M target, and after Eq. 1 the next bulk concentration is forced toward 1 M under the linear C_b-N_t relation implicit in the update. The intro's assertion 'We show results can be obtained equivalent to CμMD simulation within just two iterations' is not backed by any CμMD simulation in the paper, so CμMD equivalence is an unsupported promise rather than an external check. No load-bearing self-citation chain is present. Overall score 6: partial circularity because the method's success metric is its own control target, but the algorithm itself retains some nontrivial content (e.g., the empirical observation that the linear update stabilizes at the target in two iterations for these systems).

Assumptions & free parameters 5 free parameters · 3 assumptions · 0 invented entities

The method depends on several user-chosen convergence thresholds and on the unstated assumption that bulk concentration is a valid proxy for chemical potential. No new physical entities are introduced.

free parameters (5)
  • bulk_region_half_width = 5 Å (classical), 7 Å (MLIP)
    User-defined region at the center of the solution where concentration is measured; the choice affects the measured concentration and the convergence behavior. Defined in Sections 3.1 and 3.3.
  • moving_average_window = 5 ns
    Used to smooth the bulk concentration time series before checking the slope; user-chosen in Sections 3.1 and 3.3.
  • slope_threshold = ±2e-3 M/ns
    Defines when the moving-average slope is considered zero; user-chosen in Sections 3.1 and 3.3.
  • minimum_duration_within_threshold = 15 ns
    Required time for the slope to remain within threshold before declaring convergence; user-chosen to avoid transient convergence, Sections 3.1 and 3.3.
  • concentration_tolerance = ±0.03 M
    Stopping tolerance for the iterative loop; based on Finney et al. uncertainty, Sections 3.1 and 3.3.
assumptions (3)
  • domain assumption Bulk concentration is a valid proxy for chemical potential
    The method regulates the bulk molar concentration, not the thermodynamic chemical potential. It is assumed that the target concentration defines the target chemical potential, which requires ideal-solution behavior and a fixed activity coefficient. This enters in Section 2 and is the basis for calling the method 'constant chemical potential.'
  • domain assumption Bulk concentration scales linearly with total particle number
    The update rule Eq. (1) assumes C_b ∝ N_t. This holds only if interfacial excess adsorption is proportional to bulk concentration and there is no saturation. No test of this scaling is provided, Section 2.2.
  • domain assumption The measured bulk concentration is converged and reproducible
    The method assumes a relatively short MD run (50 ns) yields a reliable estimate of the bulk steady-state concentration. No block averaging or statistical test is reported for the classical simulations, Sections 3.1 and 4.1.

how reviews work

0 comments
Cite this review

Pith. "Pith review of A Simple Iterative Approach for Constant Chemical Potential Simulations at Interfaces." pith.science (2026). https://pith.science/paper/LF5JVGRS

@misc{pith2026250601050,
  author       = {Pith},
  title        = {Pith review of: A Simple Iterative Approach for Constant Chemical Potential Simulations at Interfaces},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/LF5JVGRS}},
  note         = {Machine review of arXiv:2506.01050}
}
read the original abstract

Chemical potential of species in solution is essential for understanding various chemical processes at interfaces. Molecular dynamics (MD) simulations, constrained by fixed compositions, cannot satisfy a constant chemical potential condition as solute species can migrate to the interface and deplete the bulk due to solute-interface interactions. In this study, we introduce a simple and computationally efficient approach named iterative constant chemical potential molecular dynamics (iCuMD) simulation, which helps simulate targeted molar concentrations of species in solution. iCuMD overcomes the limitations of conventional MD by adjusting the number of species in the solution to reach a target concentration (chemical potential). We demonstrate our approach using solid-liquid and liquid-air interfacial systems as case studies. Specifically, we perform classical force field-based MD simulations of NaCl(aq)-air and NaCl(aq)-graphite interfaces and machine learning interatomic potential (MLIP)-based MD simulations of the Na2SO4(aq)-graphene interface. Our results show that the iCuMD approach efficiently achieves the desired bulk ion concentration within two iterations and can also be integrated with MLIP-driven simulations which enable constant potential simulations with DFT-level accuracy. We show that iCuMD offers a robust and simple computational framework for constant chemical potential simulations as its only requirement is to be able to converge interfacial simulations with a measurable bulk region.

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

54 extracted references · 51 canonical work pages

  1. [1]

    European journal of physics 2006, 27 (2), 353

    Job, G.; Herrmann, F., Chemical potential—a quantity in search of recognition. European journal of physics 2006, 27 (2), 353

    work page 2006

  2. [2]

    Mrs Bulletin 2019, 44 (7), 520-523

    Chen, L.-Q., Chemical potential and Gibbs free energy. Mrs Bulletin 2019, 44 (7), 520-523

    work page 2019

  3. [3]

    R.; Wolf, R

    Ray, J. R.; Wolf, R. J., Monte Carlo simulations at constant chemical potential and pressure. The Journal of chemical physics 1993, 98 (3), 2263-2267

    work page 1993

  4. [4]

    Journal of computational chemistry 2007, 28 (10), 1763-1773

    Eslami, H.; Müller‐Plathe, F., Molecular dynamics simulation in the grand canonical ensemble. Journal of computational chemistry 2007, 28 (10), 1763-1773

    work page 2007

  5. [5]

    The Journal of chemical physics 2003, 119 (2), 978-987

    Delgado-Buscalioni, R.; Coveney, P., USHER: an algorithm for particle insertion in dense fluids. The Journal of chemical physics 2003, 119 (2), 978-987

    work page 2003

  6. [6]

    The Journal of chemical physics 1995, 102 (2), 925-931

    Lo, C.; Palmer, B., Alternative Hamiltonian for molecular dynamics simulations in the grand canonical ensemble. The Journal of chemical physics 1995, 102 (2), 925-931

    work page 1995

  7. [7]

    D.; Lu pkowski, M.; van Swol, F., Molecular dynamics and Monte Carlo simulations in the grand canonical ensemble: Local versus global control

    Papadopoulou, A.; Becker, E. D.; Lu pkowski, M.; van Swol, F., Molecular dynamics and Monte Carlo simulations in the grand canonical ensemble: Local versus global control. The Journal of chemical physics 1993, 98 (6), 4897-4908

    work page 1993

  8. [8]

    R.; McPherson, I

    Finney, A. R.; McPherson, I. J.; Unwin, P. R.; Salvalaglio, M., Electrochemistry, ion adsorption and dynamics in the double layer: a study of NaCl (aq) on graphite. Chemical science 2021, 12 (33), 11166- 11180

    work page 2021

Show all 54 references
  1. [9]

    R.; Salvalaglio, M., Bridging the gap between mesoscopic and molecular models of solid/liquid interfaces out-of-equilibrium

    Finney, A. R.; Salvalaglio, M., Bridging the gap between mesoscopic and molecular models of solid/liquid interfaces out-of-equilibrium. Chemical Engineering Research and Design 2022, 180, 285-295

    2022

  2. [10]

    R.; McPherson, I

    Finney, A. R.; McPherson, I. J.; Unwin, P. R.; Salvalaglio, M., Electrochemistry, ion adsorption and dynamics in the double layer: a study of NaCl(aq) on grap hite. Chemical Science 2021, 12 (33), 11166- 11180

    2021

  3. [11]

    R.; Elliott, J

    Di Pasquale, N.; Finney, A. R.; Elliott, J. D.; Carbone, P.; Salvalaglio, M., Constant chemical potential–quantum mechanical–molecular dynamics simulations of the graphene–electrolyte double layer. The Journal of Chemical Physics 2023, 158 (13), 134714

    2023

  4. [12]

    Journal of the American Chemical Society 2024, 146 (18), 12423-12430

    Murke, S.; Chen, W.; Pezzotti, S.; Havenith, M., Tuning Acid –Base Chemistry at an Electrified Gold/Water Interface. Journal of the American Chemical Society 2024, 146 (18), 12423-12430

    2024

  5. [13]

    Z.; Y ue, S., Dynamics of Aqueous Electrolyte S olutions: Challenges for Simulations

    Panagiotopoulos, A. Z.; Y ue, S., Dynamics of Aqueous Electrolyte S olutions: Challenges for Simulations. The Journal of Physical Chemistry B 2023, 127 (2), 430-437

    2023

  6. [14]

    T.; Baer, M

    Duignan, T. T.; Baer, M. D.; Mundy, C. J., Ions interacting in solution: Moving from intrinsic to collective properties. Current Opinion in Colloid & Interface Science 2016, 23, 58-65

    2016

  7. [15]

    T.; Mundy, C

    Duignan, T. T.; Mundy, C. J.; Schenter, G. K.; Zhao, X. S., Method for Accurately Predicting Solvation Structure. Journal of Chemical Theory and Computation 2020, 16 (8), 5401-5409

    2020

  8. [16]

    Avula, N. V . S.; Klein, M. L.; Balasubramanian, S., Understanding the Anomalous Diffusion of Water in Aqueous Electrolytes Using Machine Learned Potentials. The Journal of Physical Chemistry Letters 2023, 14 (42), 9500-9507

    2023

  9. [17]

    D.; Sumić, B.; O’Neill, N.; Schran, C.; Grey, C

    Fong, K. D.; Sumić, B.; O’Neill, N.; Schran, C.; Grey, C. P.; Michaelides, A., The Interplay of Solvation and Polarization Effects on Ion Pairing in Nanoconfined Electrolytes. Nano Letters 2024, 24 (16), 5024-5030

    2024

  10. [18]

    X.; Fong, K.; Michaelides, A.; Schran, C., To Pair or not to Pair? Machine - Learned Explicitly-Correlated Electronic Structure for NaCl in Water

    O’Neill, N.; Shi, B. X.; Fong, K.; Michaelides, A.; Schran, C., To Pair or not to Pair? Machine - Learned Explicitly-Correlated Electronic Structure for NaCl in Water. The Journal of Physical Chemistry Letters 2024, 15 (23), 6081-6091

    2024

  11. [19]

    The Journal of Physical Chemistry Letters 2024, 15 (14), 3740-3747

    Maxson, T.; Szilvási, T., Transferable Water Potentials Using Equivariant Neural Networks. The Journal of Physical Chemistry Letters 2024, 15 (14), 3740-3747

    2024

  12. [20]

    Maxson, T.; Soyemi, A.; Chen, B. W. J.; Szilvási, T., Enhancing the Quality and Reliability of Machine Learning Interatomic Potentials through Better Reporting Practices. The Journal of Physical Chemistry C 2024, 128 (16), 6524-6537. 14

    2024

  13. [21]

    F.; Goldsmith, Z

    Zhang, C.; Calegari Andrade, M. F.; Goldsmith, Z. K.; Raman, A. S.; Li, Y .; Piaggi, P. M.; Wu, X.; Car, R.; Selloni, A., Molecular-scale insights into the electrical double layer at oxide-electrolyte interfaces. Nature Communications 2024, 15 (1), 10270

    2024

  14. [22]

    Rashmi, R.; Balogun, T. O.; Azom, G.; Agnew, H.; Kumar, R.; Paesani, F., Revealing the Water Structure at Neutral and Charged Graphene/Water Interfaces through Quantum Simulations of Sum Frequency Generation Spectra. ACS Nano 2025, 19 (4), 4876-4886

    2025

  15. [23]

    The Journal of Chemical Physics 2015, 142 (14), 144113

    Perego, C.; Salvalaglio, M.; Parrinello, M., Molecular dynamics simulations of solutions at constant chemical potential. The Journal of Chemical Physics 2015, 142 (14), 144113

    2015

  16. [24]

    G.; Martínez, J

    Martínez, L.; Andrade, R.; Birgin, E. G.; Martínez, J. M., PACKMOL: A package for building initial configurations for molecular dynamics simulations. Journal of Computational Chemistry 2009, 30 (13), 2157-2164

    2009

  17. [25]

    Berendsen, H. J. C.; van der Spoel, D.; van Drunen, R., GROMACS: A message -passing parallel molecular dynamics implementation. Computer Physics Communications 1995, 91 (1), 43-56

    1995

  18. [26]

    https://github.com/aaronrfinney/CmuMD-NaCl_at_graphite (accessed 03/28)

  19. [27]

    J.; Grigera, J

    Berendsen, H. J.; Grigera, J. -R.; Straatsma, T. P., The missing term in effective pair potentials. Journal of Physical Chemistry 1987, 91 (24), 6269-6271

    1987

  20. [28]

    Hess, B.; Bekker, H.; Berendsen, H. J. C.; Fraaije, J. G. E. M., LINCS: A linear constraint solve r for molecular simulations. Journal of Computational Chemistry 1997, 18 (12), 1463-1472

    1997

  21. [29]

    S.; Cheatham III, T

    Joung, I. S.; Cheatham III, T. E., Determination of alkali and halide monovalent ion parameters for use in explicitly solvated biomolecular simulations. The journal of physical chemistry B 2008, 112 (30), 9020-9041

    2008

  22. [30]

    The Journal of Physical Chemistry B 2013, 117 (29), 8802-8813

    Wu, Y .; Aluru, N., Graphitic carbon–water nonbonded interaction parameters. The Journal of Physical Chemistry B 2013, 117 (29), 8802-8813

    2013

  23. [31]

    D.; Dix, J.; Troisi, A.; Carbone, P., Effective polarization in pairwise potentials at the graphene–electrolyte interface

    Williams, C. D.; Dix, J.; Troisi, A.; Carbone, P., Effective polarization in pairwise potentials at the graphene–electrolyte interface. The Journal of Physical Chemistry Letters 2017, 8 (3), 703-708

    2017

  24. [32]

    Physical Review B 1993, 47 (1), 558-561

    Kresse, G.; Hafner, J., Ab initio molecular dynamics for liquid metals. Physical Review B 1993, 47 (1), 558-561

    1993

  25. [33]

    Physical review B 1996, 54 (16), 11169

    Kresse, G.; Furthmüller, J., Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Physical review B 1996, 54 (16), 11169

    1996

  26. [34]

    B.; Nørskov, J

    Hammer, B.; Hansen, L. B.; Nørskov, J. K., Improved adsorption energetics within density - functional theory using revised Perdew -Burke-Ernzerhof functionals. Physical Review B 1999, 59 (11), 7413-7421

    1999

  27. [35]

    The Journal of Chemical Physics 2014, 141 (6), 064501

    Forster-Tonigold, K.; Groß, A., Dispersion corrected RPBE studies of liquid water. The Journal of Chemical Physics 2014, 141 (6), 064501

    2014

  28. [36]

    Physical Chemistry Chemical Physi cs 2020, 22 (19), 10431-10437

    Sakong, S.; Groß, A., Water structures on a Pt(111) electrode from ab initio molecular dynamic simulations for a variety of electrochemical conditions. Physical Chemistry Chemical Physi cs 2020, 22 (19), 10431-10437

    2020

  29. [37]

    The Journal of Chemical Physics 2010, 132 (15), 154104

    Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H., A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT -D) for the 94 elements H -Pu. The Journal of Chemical Physics 2010, 132 (15), 154104

    2010

  30. [38]

    E.; Martinez-Seara, H.; Jungwirth, P., Common Cations Are Not Polarizable: Effects of Dispersion Correction on Hydration Structures from Ab Initio Molecular Dynamics

    Kostal, V .; Mason, P. E.; Martinez-Seara, H.; Jungwirth, P., Common Cations Are Not Polarizable: Effects of Dispersion Correction on Hydration Structures from Ab Initio Molecular Dynamics. The Journal of Physical Chemistry Letters 2023, 14 (19), 4403-4408

    2023

  31. [39]

    L.; Rowe, P.; Müller, E

    Schran, C.; Thiemann, F. L.; Rowe, P.; Müller, E. A.; Marsalek, O.; Michaelides, A., Machine learning potentials for complex aqueous systems made simple. Proceedings of the National Academy of Sciences 2021, 118 (38), e2110077118

    2021

  32. [40]

    M.; Kovács, D

    Batatia, I.; Benner, P.; Chiang, Y .; Elena, A. M.; Kovács, D. P.; Riebesell, J.; Advincula, X. R.; Asta, M.; Avaylon, M.; Baldwin, W. J., A foundation model for atomistic materials chemistry. arXiv preprint arXiv:2401.00096 2023

    2023 arXiv

  33. [41]

    E.; Christensen, R.; Dułak, M.; Friis, J.; Groves, M

    Hjorth Larsen, A.; Jørgen Mortensen, J.; Blomqvist, J.; Castelli, I. E.; Christensen, R.; Dułak, M.; Friis, J.; Groves, M. N.; Hammer, B.; Hargus, C.; Hermes, E. D.; Jennings, P. C.; Bjerre Jensen, P.; 15 Kermode, J.; Kitchin, J. R.; Leonhard Kolsbjerg, E.; Kubal, J.; Kaasbjer...

    2017

  34. [42]

    Applied Geochemistry 2017, 86, 105-120

    Mao, S.; Peng, Q.; Wang, M.; Hu, J.; Peng, C.; Zhang, J., The PVTx properties of aqueous electrolyte solutions containing Li+, Na+, K+, Mg2+, Ca 2+, Cl− and SO42− under conditions of CO2 capture and sequestration. Applied Geochemistry 2017, 86, 105-120

    2017

  35. [43]

    J.; Rondinelli, J

    Waters, M. J.; Rondinelli, J. M., Energy contour exploration with potentiostatic kinematics. Journal of Physics: Condensed Matter 2021, 33 (44), 445901

    2021

  36. [44]

    P.; Kornbluth, M.; Molinari, N.; Smidt, T

    Batzner, S.; Musaelian, A.; Sun, L.; Geiger, M.; Mailoa, J. P.; Kornbluth, M.; Molinari, N.; Smidt, T. E.; Kozinsky, B., E(3) -equivariant graph neural networks for data -efficient and accurate interatomic potentials. Nature Communications 2022, 13 (1), 2453

    2022

  37. [45]

    P.; Abreu, C

    Eastman, P.; Galvelis, R.; Peláez, R. P.; Abreu, C. R. A.; Farr, S. E.; Gallicchio, E.; Gorenko, A.; Henry, M. M.; Hu, F.; Huang, J.; Krämer, A.; Michel, J.; Mitchell, J. A.; Pande, V . S.; Rodrigues, J. P.; Rodriguez-Guerra, J.; Simmonett, A. C.; Singh, S.; Swails, J.; Turner...

    2024

  38. [46]

    A.; Parrinello, M., PLUMED: A portable plugin for free -energy calculations with molecular dynamics

    Bonomi, M.; Branduardi, D.; Bussi, G.; Camilloni, C.; Provasi, D.; Raiteri, P.; Donadio, D.; Marinelli, F.; Pietrucci, F.; Broglia, R. A.; Parrinello, M., PLUMED: A portable plugin for free -energy calculations with molecular dynamics. Computer Physics Communications 2009, 180...

    2009

  39. [47]

    A.; Bonomi, M.; Branduardi, D.; Camilloni, C.; Bussi, G., PLUMED 2: New feathers for an old bird

    Tribello, G. A.; Bonomi, M.; Branduardi, D.; Camilloni, C.; Bussi, G., PLUMED 2: New feathers for an old bird. Computer Physics Communications 2014, 185 (2), 604-613

    2014

  40. [48]

    A.; Banáš, P.; Barducci, A.; Bernetti, M.; Bolhuis, P

    Bonomi, M.; Bussi, G.; Camilloni, C.; Tribello, G. A.; Banáš, P.; Barducci, A.; Bernetti, M.; Bolhuis, P. G.; Bottaro, S.; Branduardi, D.; Capelli, R.; Carloni, P.; Ceriotti, M.; Cesari, A.; Chen, H.; Chen, W.; Colizzi, F.; De, S.; De La Pierre, M.; Donadio, D.; Drobot, V .; E...

    2019

  41. [49]

    The Journal of Chemical Physics 2020, 153 (10), 104105

    Schran, C.; Brezina, K.; Marsalek, O., Committee neural network potentials control generalization errors and enable active learning. The Journal of Chemical Physics 2020, 153 (10), 104105

    2020

  42. [50]

    S.; Yang, P.; Baskin, A.; Lawson, J

    Yang, S.; Zhao, X.; Lu, Y .-H.; Barnard, E. S.; Yang, P.; Baskin, A.; Lawson, J. W.; Prendergast, D.; Salmeron, M., Nature of the Electrical Double Layer on Suspended Graphene Electrodes. Journal of the American Chemical Society 2022, 144 (29), 13327-13333

    2022

  43. [51]

    Chemical Reviews 1993, 93 (3), 1157-1204

    Ohtaki, H.; Radnai, T., Structure and dynamics of hydrated ions. Chemical Reviews 1993, 93 (3), 1157-1204

    1993

  44. [52]

    R.; Pettitt, B

    Cannon, W. R.; Pettitt, B. M.; McCammon, J. A., Sulfate Anion in Water: Model Structural, Thermodynamic, and Dynamic Properties. The Journal of Physical Chemistry 1994, 98 (24), 6225-6230

    1994

  45. [53]

    P.; Vega, C.; Gallo, P., Impact of sulfate salts on water structure: insights from molecular dynamics

    Lamas, C. P.; Vega, C.; Gallo, P., Impact of sulfate salts on water structure: insights from molecular dynamics. Molecular Physics, e2406260

  46. [54]

    M.; Abascal, J

    Zeron, I. M.; Abascal, J. L. F.; Vega, C., A force field of Li+, Na+, K+, Mg2+, Ca2+, Cl−, and SO42− in aqueous solution based on the TIP4P/2005 water model and scaled charges for the ions. The Journal of Chemical Physics 2019, 151 (13), 134504

    2005

Pith tools

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