Pith. sign in

REVIEW 5 major objections 5 minor 32 references

Electrically driven first-order phase transition of a 2D ionic crystal at the electrode/electrolyte interface

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

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

Pith's one-line read Applying a voltage drives the adsorbed LiCl layer at an aluminum electrode through two distinct ordering stages, with the final step a genuine first-order phase transition.

desk verdict Plausible two-stage ordering picture, but the 'first-order' classification rests on thin finite-size evidence; worth reviewing with a demand for error bars and a structural order parameter. read the letter →

arxiv 2507.19087 v1 pith:Y5HZE4DX submitted 2025-07-25 physics.chem-ph cond-mat.mtrl-sci

Federica Angiolari , Alessandro Coretti , Mathieu Salanne , Sara Bonella This is my paper · ORCID
classification physics.chem-phcond-mat.mtrl-sci
keywords moltensaltelectrolyteelectrode/electrolyteinterfacetwo-dimensionalioniccrystalvoltage-drivenphasetransitionfirst-orderconstant-potentialmoleculardynamicsinterfacialcapacitancefinite-sizeeffects
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

The paper claims that a voltage applied across a molten LiCl/aluminum interface does not order the adsorbed ion layer in one step. Instead, as the electrode potential rises, the layer first turns from a liquid-like state into a polycrystalline mosaic, and only then locks into a single square-lattice crystal. The first of these transitions is continuous (or weakly first-order); the second is first-order, with a latent energy of about 3.4 kJ/mol and a free-energy barrier that grows from about 0.25 to 2 kcal/mol when the simulation cell is enlarged. The final step also shows a sharp, size-dependent peak in the differential capacitance. This matters because it predicts observable signatures—voltage shift, capacitance spike, and hysteresis—for voltage-controlled ordering in molten-salt electrochemistry.

What carries the argument

The argument is carried by treating the electrode surface charge density $\sigma$ as the order parameter. At fixed applied voltage $\Delta\Psi$, the simulations build the probability distribution $P(\sigma|\Delta\Psi)$, either directly or by histogram reweighting (WHAM), and the free energy $F(\sigma)=-k_B T\ln P(\sigma|\Delta\Psi)$; two coexisting minima with a barrier between them signal a first-order transition, and the barrier grows because creating an interface between coexisting phases costs more in larger cells. The differential capacitance $C_{\mathrm{diff}} = \frac{S}{k_B T}\langle(\delta\sigma)^2\rangle$ is the conjugate response function, and its peak sharpening and growth with system size provide the thermodynamic signature. Structural identification uses two-dimensional structure-factor peak intensities and orientational correlation functions $g_4(r)$ and $g_6(r)$, which distinguish square-crystal, hexatic/liquid-like, and mixed polycrystalline states.

What would settle it

Run the large system (or a larger one) at the transition voltage for much longer than the reported runs, or use an enhanced-sampling bias on a different reaction coordinate: if the system switches between polycrystal and monocrystal repeatedly with a barrier that does not keep growing with cell area, or if the capacitance peak stops sharpening, the first-order classification is wrong.

Watch

Extended reading notes

Core claim

On the paper's own terms, the central discovery is a two-stage crystallization path for the adsorbed chloride sublattice at the positive aluminum electrode. At low potentials the Cl− layer is liquid-like with local sixfold correlations; over a wide intermediate potential window it becomes a mosaic of small tetratic and hexatic domains (a polycrystal); above a threshold voltage the mosaic abruptly converts to a single square-lattice monocrystal. The paper's explicit classification is that the liquid-like to polycrystal transition is continuous (or weakly first-order), while the polycrystal-to-monocrystal transition is first-order. The first-order assignment is supported by finite-size analysis: the free-energy barrier between the two phases grows from about 0.25 kcal/mol in the small cell to about 2 kcal/mol in the medium cell, the large cell remains trapped in the phase it starts in, and the differential capacitance develops a sharper, higher peak at the transition as the system grows. The latent energy of the final transition, 3.4 kJ/mol, is about six times smaller than bulk LiCl melting, consistent with a two-dimensional crystal.

Load-bearing premise

The first-order classification rests on the assumption that the large cell's inability to switch phases and the growing barrier are signs of a genuine nucleation barrier, rather than slow relaxation or an order parameter that misses the true reaction coordinate.

Editorial extensions

If this is right

  • At intermediate potentials the adsorbed layer would appear as a polycrystal—a mosaic of tetratic and hexatic domains—so experiments should see preordered patches before full crystallization, not a single disorder-to-monocrystal jump.
  • The measured onset voltage would depend on electrode size: small cells order near 0.4–0.8 V while the medium cell orders near 0.75 V, so predictions for real electrodes must account for finite-size shifts.
  • The differential capacitance would develop a voltage peak at the monocrystallization transition, and the peak would sharpen and grow with system size, making capacitance measurements a direct way to locate the transition.
  • The latent heat of the 2D ordering is low (3.4 kJ/mol, roughly six times below bulk LiCl melting), which limits supercapacitor energy storage but the accompanying charge jump could drive sensors or actuators.

Reading between the lines

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

  • If the first-order classification carries to macroscopic electrodes, the barrier growth implies pronounced hysteresis in cyclic voltammetry: forward and reverse crystallizations would occur at noticeably different voltages, a testable experimental signature.
  • The intermediate polycrystal with coexisting tetratic and hexatic domains is richer than standard two-dimensional melting pictures; a natural extension is to test whether other polarizable ionic adsorbates, such as iodide or bromide on metals, show the same mosaic stage.
  • Because the surface charge density is the thermodynamic conjugate of the voltage, a first-order transition implies a discontinuous jump in electrode charge, so the interface may act as a voltage-triggered switch; this is an inference beyond the paper's explicit claims.
  • The asymmetry between chloride ordering and lithium non-ordering suggests the phenomenon depends on ion polarizability and surface interaction; replacing LiCl with a more symmetric salt might lead to ordering at both electrodes, a possibility the authors themselves flag as speculative.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

5 major / 5 minor

Summary. The paper reports constant-potential molecular dynamics simulations of molten LiCl on an Al(100) electrode at three system sizes, studying the voltage-driven ordering of the adsorbed layer. Through 2D structure factors and orientational correlation functions, the authors identify a two-stage transition: from a liquid-like layer to a polycrystalline mosaic, then to a monocrystalline layer. Using histogram reweighting on the surface charge density sigma, they compute free energy surfaces at the transition potentials, estimate a latent energy for the second transition, and analyze the differential capacitance. The central claim is that the polycrystal-to-monocrystal transition is first-order, supported by an increasing free energy barrier with system size (0.25 to 2 kcal/mol) and the failure of the largest cell to switch phases in brute-force MD.

Significance. If established, the two-stage mechanism and first-order final crystallization would be an important contribution to the understanding of voltage-driven ordering at electrochemical interfaces, with implications for molten salt technologies. The simulation approach—constant-potential MD with polarizable force fields and histogram reweighting—is appropriate, and the structural analysis via g4/g6 orientational correlations provides a clear qualitative picture. However, the evidence for the first-order classification is currently underdetermined: the free energy barriers lack error bars, the reaction-coordinate validity of sigma is not verified, and the large-cell observation is a null result. The paper's contributions would be significantly strengthened by quantitative error analysis and a validation of the order parameter; as written, the conclusions exceed the support.

major comments (5)
  1. [II.C, Fig. 5] The free energy barriers reported for the small and medium systems are single WHAM estimates with no statistical uncertainty. The small-system barrier of approximately 0.25 kcal/mol is more than an order of magnitude below kBT approximately 2.4 kcal/mol at 1200 K, so the two minima are not resolved and the statement that the two states are 'equally probable' at ΔΨ = 0.80 V is not supported. Please report error bars (e.g., from g_wham bootstrap or block averaging) and quantify the resolution limit of the histogram method.
  2. [II.C, last paragraph] The inference that the large-cell metastability 'points to a further increase of the free energy barrier, consistent with the first-order character' is based on a null observation: the system remained in one phase over the simulated time at a single potential (1.375 V). This behavior is equally compatible with slow relaxation, near-critical slowing down at a continuous transition, or simply insufficient sampling; it does not by itself establish a growing thermodynamic barrier. An enhanced-sampling calculation (e.g., umbrella sampling along a structural order parameter) or a measurement of the transition rate for the large system would be needed to support this claim. In addition, the capacitance value for the large system (88.6 µF/cm²) is estimated by assuming the transition occurs at the same voltage as the medium system, an assumption that is unjustified given the documented shift of the transition voltage with system size.
  3. [II.B] The free energy is computed from P(σ|ΔΨ), with σ the surface charge density, but the polycrystal-to-monocrystal transition is defined by orientational order (g4). It is not demonstrated that σ distinguishes the polycrystal and monocrystal states or that the two minima in F(σ) correspond to these structural phases. Please validate σ as a reaction coordinate, for example by computing the joint distribution P(σ, ψ4) or by plotting the average bond-orientational order parameter along the WHAM path. Without this, the barrier and latent heat cannot be attributed to the structural transition.
  4. [II.B, Eq. (4)] The latent energy estimate ΔE_T2 = 0.038 J/m² (L = 3.4 kJ/mol) is a single number with no uncertainty, and the formula includes only the change in electrostatic energy of the capacitor, not the full enthalpy difference between the two phases (e.g., changes in ion-ion interaction energy). The comparison with the bulk melting enthalpy (19.9 kJ/mol) is suggestive but does not prove first-order character, since a continuous transition can also involve a small energy release. Please provide error estimates and discuss which energy contributions are included in this latent energy.
  5. [II.C] The reported barrier increase (approximately 0.25 kcal/mol to approximately 2 kcal/mol) over a factor-of-2 increase in lateral length is not compared to any finite-size scaling law. For a first-order transition with a line-tension barrier, one would expect the barrier to scale with the interface length (or with a nontrivial power of system size for nucleation); a factor-of-8 increase over a factor-of-2 length increase is not matched to such a law. A quantitative scaling analysis, or at least a discussion of the expected scaling, is needed to support the first-order interpretation.
minor comments (5)
  1. [Abstract and Conclusion] The word 'proves' is used in the abstract, Significance Statement, and Conclusion for the first-order character of the transition; this is stronger than the evidence presented. Suggest replacing with 'indicates' or 'is consistent with' throughout.
  2. [II.A] The inability to distinguish exponential vs power-law decay of g6 means the 'liquid-like' assignment is provisional; please state this more explicitly in the main text and figure captions, rather than only in a parenthetical remark.
  3. [IV.A] The total simulation times (200-600 ps of dynamics) are short relative to typical nucleation timescales. Please provide a table or list of trajectory lengths for each system and potential so the reader can assess the metastability claim, especially for the large system.
  4. [IV.B and Fig. 4] The histogram reweighting/WHAM calculation is described only briefly; the number of windows, the overlap between histograms, and convergence criteria are not reported. Please include these details in the Methods or SI.
  5. [Introduction, ref. [8]] The term 'mean solidification length' is used without definition; please define it or rephrase the sentence for clarity.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the first-order classification is inferred from computed free-energy surfaces, barrier trends, and metastability observations, not from fitted parameters or self-citations.

full rationale

The paper's central claim—that the polycrystal-to-monocrystal transition is first-order—is derived from simulation data rather than assumed or fitted. The free energy F(σ) = −kBT ln P(σ|ΔΨ) is obtained via histogram reweighting (WHAM) from constant-potential MD trajectories, and the barrier heights (0.25 kcal/mol in the small cell, 2 kcal/mol in the medium cell) are computed, not imposed. The latent heat is estimated from the observed charge-density jump across the transition, ΔE_T2 = [σ(ΔΨ+_T2) − σ(ΔΨ−_T2)] ΔΨ_T2, using values read from the simulated probability distributions. The capacitance is computed from the fluctuation formula C_diff = S/(kBT) ⟨(δσ)^2⟩, which is an identity in the constant-potential ensemble, and its peak sharpening is a measured trend across system sizes. The force field is an independent external input obtained from DFT-based fitting (Refs. 10, 24), and no parameter is tuned to reproduce the transition order, the barrier height, or the latent heat. Self-citations appear only as methodological references (MetalWalls, MaZe, the WHAM protocol, charge-fluctuation formalism) and are not load-bearing for the physical conclusion. The authors' own stated limitations—the small-system barrier is below kBT, the large system failed to switch phases, and the order parameter σ may not capture the full reaction coordinate—affect the confidence and strength of the first-order inference, but they do not make the derivation circular, because the classification is not built into the inputs. The analysis is a standard, self-contained computational thermodynamics study; no step reduces by construction to its own inputs.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

The central claims rest on the classical constant-potential MD model and its force field, not on any new fitted parameter introduced here. No free parameters are introduced in this paper; the order-parameter cutoff is a standard structural choice. The listed axioms are the model assumptions the conclusions inherit.

assumptions (4)
  • domain assumption The polarizable classical force field of Refs. [10,24], fitted to DFT, accurately models LiCl adsorption on Al(100) at 1200 K.
    All conclusions are conditional on this force field; the system is a model system, with the electrode atoms fixed while real Al would melt at 1200 K, as stated in the Introduction.
  • domain assumption The constant-potential fluctuating-charge ensemble (Siepmann-Sprik) and Mass-Zero constrained dynamics sample equilibrium distributions of electrode charge and melt polarization.
    Used in Methods IV.A; the free-energy analysis assumes this ensemble's statistics.
  • domain assumption Applied potential DeltaPsi serves as a valid biasing variable, and WHAM/histogram reweighting converges with 200-600 ps trajectories.
    Methods IV.B; the free-energy barriers and transition voltages are computed this way, without reported statistical error bars.
  • domain assumption The structural order parameters (2D structure factor peaks and psi4/psi6 with a 5 A cutoff) are sufficient to distinguish the liquid-like, polycrystalline, and monocrystalline states.
    Used in Section II.A; the hexatic vs liquid distinction is admitted to be unresolved due to small size.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Electrically driven first-order phase transition of a 2D ionic crystal at the electrode/electrolyte interface." pith.science (2026). https://pith.science/paper/Y5HZE4DX

@misc{pith2026250719087,
  author       = {Pith},
  title        = {Pith review of: Electrically driven first-order phase transition of a 2D ionic crystal at the electrode/electrolyte interface},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/Y5HZE4DX}},
  note         = {Machine review of arXiv:2507.19087}
}
read the original abstract

Liquid electrolytes adsorbed at the surface of metallic electrodes display a multitude of structures that can largely differ from the parent bulk system, both in terms of composition and local organization. In particular, the existence of disorder-order or order-order transitions has been increasingly reported in experimental and simulation studies, and the electrode potential identified as the corresponding driving force. The microscopic mechanisms and the stages of the process are, however, poorly understood, and the free energy variation during the transition remain insufficiently characterized. To fill this gap, we investigate the crystallization of the adsorbed layer in a prototypical molten salt-metal interface. We demonstrate that the transition from the disordered to ordered structures proceeds in two stages. Pre-ordering effects are observed across a wide range of potentials, resulting in the formation of a poly-crystalline structure on the electrode surface, before an abrupt ordering transition finally occurs. The pre-ordering displays signature of a continuous transition. On the other hand, finite-size effects analysis proves the first-order character of the transition towards the mono-crystalline state. Upon increasing the system size, a shift in the onset applied voltage is observed, accompanied by a dramatic increase in the free energy barrier. The latter reflects in the interfacial capacitance, which displays a peak that sharpens with increasing system size.

Figures

Figures reproduced from arXiv: 2507.19087 by the authors.

Figure 1
Figure 1. Instantaneous in-plane structure factors for the adsorbed layer on the positive electrode [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Top: Variation with time of the intensities of three characteristic 2D structure factor [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. Orientational correlation functions of the adsorbed chloride ions, [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figures from the paper (2 more)
Figure 4
Figure 4. Figure 4: (A) Probability P(σ | ∆Ψ) of the charge density on the electrodes at a given applied potential. (B) Differential capacitance as a function of applied voltage. 10 [PITH_FULL_IMAGE:figures/full_fig_p010_4.png]
Figure 5
Figure 5. Figure 5: Free energy profile for the medium (A) and small (B) size systems computed at the [PITH_FULL_IMAGE:figures/full_fig_p012_5.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

32 extracted references · 32 canonical work pages

  1. [1]

    Magnussen

    O.M. Magnussen. Ordered anion adlayers on metal electrode surfaces. Chem. Rev., 102(3):679– 726, 2002

    work page 2002

  2. [2]

    Pan and W

    G.-B. Pan and W. Freyland. 2D phase transition of PF 6 adlayers at the electrified ionic liquid/Au(1 1 1) interface. Chem. Phys. Lett. , 427(1–3):96–100, 2006

    work page 2006

  3. [3]

    Su, Y.-C

    Y.-Z. Su, Y.-C. Fu, J.-W. Yan, Z.-B. Chen, and B.-W. Mao. Double layer of Au(100)/ionic liquid interface and its stability in imidazolium-based ionic liquids. Angew. Chem., Int. Ed. , 48:5148–5151, 2009

    work page 2009

  4. [4]

    R. Wen, B. Rahn, and O. M. Magnussen. Potential-dependent adlayer structure and dynamics at the ionic liquid/au(111) interface: A molecular-scale in situ video-stm study.Angew. Chem., Int. Ed. , 54:6062–6026, 2015

    work page 2015

  5. [5]

    G. Shi, L. Chen, Y. Yang, D. Li, Z. Qian, S. Liang, L. Yan, L. H. Li, M. Wu, and H. Fang. Two-dimensional Na–Cl crystals of unconventional stoichiometries on graphene surface from dilute solution at ambient conditions. Nat. Chem. , 10:776–779, 2018

    work page 2018

  6. [6]

    R. Yi, J. Jiang, Y. Yang, Y. Zhang, S. Gao, Y. Zhao, J. Hu, X. Su, X. Cia, B. Peng, F. Dai, P. Li, Z. Guan, H. Yang, F. Zhu, J. Cao, Z. Wang, H. Fang, L. Zhang, and L. Chen. Two- dimensional anion-rich NaCl2 crystal under ambient conditions. Nat. Commun., 16:464, 2025

    work page 2025

  7. [7]

    J. Wang, H. Li, M. Tavakol, A. Serva, B. Nener, G. Parish, M. Salanne, G. G. Warr, K. Vo¨ ıtchovsky, and R. Atkin. Ions adsorbed at amorphous solid/solution interfaces form wigner crystal-like structures. ACS Nano, 18:1181–1194, 2024. 16

    work page 2024

  8. [8]

    Comtet, A

    J. Comtet, A. Nigu` es, V. Kaiser, B. Coasne, L. Bocquet, and A. Siria. Nanoscale capillary freezing of ionic liquids confined between metallic interfaces and the role of electronic screening. Nat. Mater. , 16:634–639, 2017

    work page 2017

Show all 32 references
  1. [9]

    Wandlowski, J

    T. Wandlowski, J. X. Wang, and B. M. Ocko. an equilibrium lattice-gas model. J. Electroanal. Chem., 500:418–434, 2001

    2001

  2. [10]

    Pounds, S

    M. Pounds, S. Tazi, M. Salanne, and PA. Madden. Ion adsorption at a metallic electrode: an ab initio based simulation study. J Phys Condens Matter. , 2009

    2009

  3. [11]

    Kirchner, T

    K. Kirchner, T. Kirchner, V. Ivanistsev, and M. V. Fedorov. Electrical double layer in ionic liquids: structural transitions from multilayer to monolayer structure at the interface. Elec- trochim. Acta, 110:762—771, 2013

    2013

  4. [12]

    Limmer, Mathieu Salanne, Ren´ e van Roij, Paul A

    C´ eline Merlet, David T. Limmer, Mathieu Salanne, Ren´ e van Roij, Paul A. Madden, David Chandler, and Benjamin Rotenberg. The electric double layer has a life of its own. The Journal of Physical Chemistry C , 118(32):18291–18298, 2014

    2014

  5. [13]

    Ilja Siepmann and Michiel Sprik

    J. Ilja Siepmann and Michiel Sprik. Influence of surface topology and electrostatic potential on water/electrode systems. The Journal of Chemical Physics , 102(1):511–524, January 1995

    1995

  6. [14]

    Reed, Oliver J

    Stewart K. Reed, Oliver J. Lanning, and Paul A. Madden. Electrochemical interface be- tween an ionic liquid and a model metallic electrode. The Journal of Chemical Physics , 126(8):084704, February 2007

    2007

  7. [15]

    Coretti, L

    A. Coretti, L. Scalfi, C. Bacon, B. Rotenberg, R. Vuilleumier, G. Ciccotti, and M. Salanne. Mass-zero constrained molecular dynamics for electrode charges in simulations of electrochem- ical systems. J. Chem. Phys. , 152:194701, 2020

    2020

  8. [16]

    Coretti, C

    A. Coretti, C. Bacon, R. Berthin, A. Serva, L. Scalfi, I. Chubak, K. Goloviznina, M. Hae- fele, A. Marin-Lafl` eche, B. Rotenberg, S. Bonella, and M. Salanne. MetalWalls: simulating electrochemical interfaces between polarizable electrolytes and metallic electrodes. J. Chem. P...

    2022

  9. [17]

    Reed, Alessandra Serva, Roxanne Berthin, Camille Bacon, Sara Bonella, Benjamin Rotenberg, Paul A

    Abel Marin-Lafl` eche, Matthieu Haefele, Laura Scalfi, Alessandro Coretti, Thomas Dufils, Guillaume Jeanmairet, Stewart K. Reed, Alessandra Serva, Roxanne Berthin, Camille Bacon, Sara Bonella, Benjamin Rotenberg, Paul A. Madden, and Mathieu Salanne. Metalwalls: A classical mol...

    2020

  10. [18]

    J. M. Kosterlitz and D. J. Thouless. Ordering, metastability and phase transitions in two- dimensional systems. J. Phys. C: Solid State Phys. , 6:1181, 1973

    1973

  11. [19]

    B. I. Halperin and D. R. Nelson. Theory of two-dimensional melting. Phys. Rev. Lett., 41(121), 1978

    1978

  12. [20]

    A. P. Young. Melting and the vector coulomb gas in two dimensions. Phys. Rev. B , 19:1855, 1979

    1979

  13. [21]

    B. K. Clark, M. Casula, and D. M. Ceperley. Hexatic and mesoscopic phases in a 2d quantum coulomb system. Phys. Rev. Lett. , 103:055701, 2009

    2009

  14. [22]

    Y.-W. Li, Y. Yao, and M. Pica Ciamarra. Two-dimensional melting of two- and three- component mixtures. Phys. Rev. Lett. , 130:258202, 2023

    2023

  15. [23]

    L¨ offler, L

    R. L¨ offler, L. Siedentop, and P. Keim. Tetratic phase in 2D crystals of squares. Soft Matter , 21:2026, 2025

    2026

  16. [24]

    S. Tazi, M. Salanne, C. Simon, P. Turq, M. Pounds, and P.A. Madden. Potential-induced ordering transition of the adsorbed layer at the ionic liquid/electrified metal interface. The Journal of Physical Chemistry B , 114(25):8453–8459, 2010

    2010

  17. [25]

    D. T. Limmer, C. Merlet, M. Salanne, D. Chandler, P. A. Madden, R. van Roij, and B. Roten- berg. Charge fluctuations in nanoscale capacitors. Phys. Rev. Lett. , 111:106102, 2013

    2013

  18. [26]

    David R. Lide. CRC Handbook of Chemistry and Physics . CRC Press, Taylor & Francis, Boca Raton, FL, 86th edition, 2005. p. 6-110

    2005

  19. [27]

    Scalfi, D

    L. Scalfi, D. T. Limmer, A. Coretti, S. Bonella, P. A. Madden, M. Salanne, and B. Rotenberg. Charge fluctuations from molecular simulations in the constant-potential ensemble. Phys. Chem. Chem. Phys. , 22:10480–10489, 2020

    2020

  20. [28]

    K. Binder. Finite size effects of phase transitions. Ferroelectrics, 73:43–67, 1987

    1987

  21. [29]

    M. Sprik. Computer simulation of the dynamics of induced polarization fluctuations in water. The Journal of Physical Chemistry , 95(6):2283–2291, 1991

    1991

  22. [30]

    Ryckaert, A

    J-P. Ryckaert, A. Bellemans, and G. Ciccotti and. The rotation-translation coupling in di- atomic molecules. Molecular Physics, 44(4):979–996, 1981

    1981

  23. [31]

    Coretti, S

    A. Coretti, S. Bonella, and G. Ciccotti. Communication: Constrained molecular dynamics for polarizable models. The Journal of Chemical Physics , 149:191102, 11 2018. 18

    2018

  24. [32]

    Hub, Bert L

    Jochen S. Hub, Bert L. de Groot, and David van der Spoel. g wham—a free weighted his- togram analysis implementation including robust error and autocorrelation estimates. Journal of Chemical Theory and Computation , 6(12):3713–3720, 2010. 19 Supporting Information for Electric...

    2010

Pith tools

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