Simulations of Crystal Nucleation from Solution at Constant Chemical Potential
Pith reviewed 2026-05-25 00:13 UTC · model grok-4.3
The pith
A method adapts constant chemical potential molecular dynamics to study crystal nucleation from solution without depleting the solute concentration.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
By adopting the constant chemical potential molecular dynamics approach, the method allows determining the crystal nucleus size and nucleation rates at constant supersaturation, as shown in the example of homogeneous nucleation of sodium chloride from its supersaturated aqueous solution.
What carries the argument
Constant chemical potential molecular dynamics, which maintains fixed solute chemical potential by particle exchange with a reservoir during the nucleation process.
If this is right
- Crystal nucleus size can be extracted under unchanging supersaturation.
- Nucleation rates can be computed without progressive depletion of the solution.
- Thermodynamic conditions remain those of the target supersaturated state throughout the run.
- The technique applies directly to homogeneous nucleation of ionic salts from aqueous solution.
Where Pith is reading between the lines
- The method could make direct comparison with experiments performed at fixed concentration more straightforward.
- It opens the possibility of mapping how nucleation barriers vary with supersaturation in a single consistent framework.
- Similar constant-potential control might be tested on systems where solute depletion is even more pronounced, such as larger organic molecules.
Load-bearing premise
The constant chemical potential molecular dynamics approach can be directly adopted to nucleation studies while preserving accurate thermodynamic conditions and without introducing new artifacts specific to the nucleation process.
What would settle it
A simulation in which the solution concentration still changes appreciably during nucleus growth despite the constant-chemical-potential protocol, or nucleation rates that fail to match independent constant-supersaturation measurements.
read the original abstract
A widely spread method of crystal preparation is to precipitate it from a supersaturated solution. In such a process, control of solution concentration is of paramount importance. Nucleation process, polymorph selection, and crystal habits depend crucially on this thermodynamic parameter. When performing simulations in the canonical ensemble as the crystalline phase is deposited the solution is depleted of solutes. This unavoidable modification of the thermodynamic conditions leads to significant artifact. Here we adopt the idea of the constant chemical potential molecular dynamics approach of Perego et al. [J. Chem. Phys. 2015, 142, 144113] to the study of nucleation. Our method allows determining the crystal nucleus size and nucleation rates at constant supersaturation. As an example we study the homogeneous nucleation of sodium chloride from its supersaturated aqueous solution.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript adapts the constant chemical potential molecular dynamics (CμMD) scheme of Perego et al. (J. Chem. Phys. 2015) to nucleation from solution. The central claim is that the adapted method maintains constant supersaturation during crystal growth, thereby eliminating depletion artifacts present in canonical-ensemble simulations and enabling reliable computation of nucleus size and nucleation rates; NaCl homogeneous nucleation from aqueous solution is presented as the working example.
Significance. If the adaptation preserves accurate thermodynamic control without introducing new dynamical artifacts, the approach would be a useful methodological contribution for precipitation studies where supersaturation must remain fixed. The paper correctly identifies the solute-depletion problem in standard MD and attempts to transplant an existing constant-μ technique; however, the significance is tempered by the absence of explicit checks that the reservoir coupling leaves nucleation kinetics unaltered.
major comments (2)
- [Methods / Results] The manuscript provides no systematic invariance tests (e.g., variation of reservoir volume, exchange frequency, or coupling strength) showing that nucleus-size distributions or mean-first-passage times remain unchanged. Such tests are required to establish that the CμMD reservoir does not modify attachment/detachment rates at the crystal-solution interface, which directly affects the central claim of artifact-free nucleation at constant supersaturation.
- [Results] No comparison is reported between the adapted CμMD nucleation rates or critical nucleus sizes and either (i) independent constant-supersaturation calculations or (ii) the known limiting behavior at very large reservoir volumes. Without such benchmarks the claim that the method “allows determining the crystal nucleus size and nucleation rates at constant supersaturation” rests on the unverified assumption that the Perego et al. scheme transplants without new artifacts.
minor comments (1)
- [Methods] Notation for the chemical-potential control parameters (e.g., reservoir size, bias strength) should be defined explicitly in the text rather than only by reference to Perego et al.
Simulated Author's Rebuttal
We thank the referee for their careful reading and constructive comments. We address the major comments point by point below, agreeing where additional validation is warranted and outlining the revisions.
read point-by-point responses
-
Referee: [Methods / Results] The manuscript provides no systematic invariance tests (e.g., variation of reservoir volume, exchange frequency, or coupling strength) showing that nucleus-size distributions or mean-first-passage times remain unchanged. Such tests are required to establish that the CμMD reservoir does not modify attachment/detachment rates at the crystal-solution interface, which directly affects the central claim of artifact-free nucleation at constant supersaturation.
Authors: We agree that explicit invariance tests are necessary to confirm the reservoir coupling does not alter interface kinetics. In the revised manuscript we will add a dedicated subsection with simulations at multiple reservoir volumes, exchange frequencies, and coupling strengths. These will demonstrate that nucleus-size distributions and mean-first-passage times remain unchanged within statistical uncertainty, directly supporting that attachment/detachment rates are unaffected. revision: yes
-
Referee: [Results] No comparison is reported between the adapted CμMD nucleation rates or critical nucleus sizes and either (i) independent constant-supersaturation calculations or (ii) the known limiting behavior at very large reservoir volumes. Without such benchmarks the claim that the method “allows determining the crystal nucleus size and nucleation rates at constant supersaturation” rests on the unverified assumption that the Perego et al. scheme transplants without new artifacts.
Authors: We acknowledge the benefit of explicit benchmarks. The revised manuscript will include additional runs at substantially larger reservoir volumes to verify convergence toward the expected limiting behavior. We will also place our computed nucleus sizes and rates in context with the constant-supersaturation limit and any relevant literature values for NaCl, thereby providing the requested comparisons. revision: yes
Circularity Check
No circularity; method adopted from independent prior reference
full rationale
The paper's core contribution is the direct adoption and application of the constant-chemical-potential MD scheme introduced in Perego et al. (JCP 2015). No equations, parameters, or predictions in the present work are shown to reduce by construction to quantities defined or fitted within this manuscript. The derivation chain consists of citing an external method and then using it to compute nucleus sizes and rates for NaCl; this is self-contained against the cited benchmark and does not exhibit self-definition, fitted-input-as-prediction, or load-bearing self-citation chains. Score 0 is the appropriate default when the central claim rests on independently published methodology rather than tautological re-derivation.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The constant chemical potential MD framework of Perego et al. can be adapted to nucleation without loss of accuracy in thermodynamic driving force.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.