Skip to content

Isotope effects in water: differences of structure, dynamics, spectrum, and proton transport between heavy and light water from ReaxFF reactive force field simulations

Non-primary PDF

This corpus path is an AIP proof PDF for the JPCL letter. For stable pagination, figures, and final copy-editing, prefer [[20180000-0002-5255-7340-j-phys-chem-isotope-effects]] when that slug’s pdf_path points at the version-of-record file (Zhang_HeavyWater_JPC_Letters_2018.pdf). See docs/corpus/NON_PRIMARY_ARTICLE_PAPER_SLUGS.md (proof duplicate pattern).

Summary

Heavy water (D\(_2\)O) and light water (H\(_2\)O) differ in hydrogen bonding, diffusion, spectroscopic signatures, and proton (or deuteron) transport mechanisms—effects that are partly quantum-mechanical in origin. This letter reports ReaxFF reactive molecular dynamics with an isotope-aware parametrization that encodes H versus D differences in selected force-field terms so that classical simulations can recover key isotope effects without full path-integral molecular dynamics. The study compares bulk H\(_2\)O and D\(_2\)O structure (radial distribution functions), dynamics (diffusion), vibrational spectra, and Grotthuss-style proton/deuteron transport, and it highlights solvation differences for H\(_3\)O\(^+\) versus D\(_3\)O\(^+\) in the bulk solution picture presented in the abstract.

The practical motivation is selective labeling: many experiments compare H and D in the same environment, so a simulation framework that yields zero isotope splitting by construction cannot be aligned with those measurements. The authors therefore augment the aqueous ReaxFF description with explicit mass- and interaction-level distinctions intended to reproduce macroscopic ordering (for example diffusion constants) while remaining computationally lighter than PIMD.

Methods

Force field. The work builds on the CHON-2017_weak-class ReaxFF aqueous framework with H/D-specific adjustments documented in the letter and Supporting Information.

Simulations. Bulk H\(_2\)O and D\(_2\)O cells are equilibrated and sampled under protocols stated in the article (ensemble, timestep, cutoffs, and duration); spectral analysis and transport diagnostics follow standard post-processing for ReaxFF trajectories.

Comparators. Results are discussed relative to experimental ordering of diffusion constants and spectroscopic trends, as cited in the JPCL text.

Reproducibility note. Use the VOR PDF on the sibling page for figure numbers and SI pointers; proof PDFs can differ in layout.

Implementation detail. Reproduce the letter’s H vs D comparisons using identical statistical ensembles and identical numerical settings aside from isotope labels, so differences trace to the force field and mass choices rather than to different thermostats or initial configurations.

MD protocol (same letter as VOR)

  • Engine / code: LAMMPS molecular dynamics with ReaxFF (CHON-2017_weak-derived H/D parameters).
  • System size & composition: Bulk cubic supercells of light and heavy water; atom counts on 20180000-0002-5255-7340-j-phys-chem-isotope-effects / SI.
  • Boundaries / periodicity: 3D PBC.
  • Ensemble: NVT for paired H\(_2\)O vs D\(_2\)O runs.
  • Timestep / duration: Femtosecond timestep; equilibration and production lengths in SI.
  • Thermostat: Nose–Hoover or Berendsen as published.
  • Barostat: N/A — constant-volume NVT without NPT barostat unless SI states otherwise.
  • Temperature: K setpoints identical between isotopes aside from labels.
  • Pressure: N/A — no applied pressure targets for these cells.

Findings

ReaxFF MD captures major light vs heavy water differences emphasized in the abstract, including diffusion ordering and aspects of vibrational and transport behavior. Structural distinctions between D\(_3\)O\(^+\) and H\(_3\)O\(^+\) appear in the simulated bulk environment. The authors describe Grotthuss hopping for proton transport as captured appropriately within their parametrization, framing the model for isotope-labeled studies in complex aqueous environments.

Limitation stated in spirit of the letter. The approach improves classical isotope sensitivity but does not replace full nuclear quantization; users should still cross-check PIMD or experiment when ZPE and tunneling dominate (for example very low temperatures or ultrafast spectroscopies).

Navigation. Always reconcile parameter tables and figure numbering with the VOR sibling page [[20180000-0002-5255-7340-j-phys-chem-isotope-effects]] before citing specific Supporting Information entries in manuscripts or theses.

Limitations

Letter-length exposition compresses numerical settings; proof pagination may not match the final issue. Nuclear quantum effects beyond the isotope-aware classical correction are not treated at the PIMD level.

Relevance to group

Zhang / van Duin aqueous ReaxFF line extended to isotope effects; this slug is duplicate proof registration for manifest provenance—cite the VOR sibling for scholarly reference lists.

Citations and evidence anchors

Reader notes (navigation)

  • Version of record page: [[20180000-0002-5255-7340-j-phys-chem-isotope-effects]] Proof duplicate. Use [[20180000-0002-5255-7340-j-phys-chem-isotope-effects]] for VOR pagination and Supporting Information pointers for H/D ReaxFF settings.

Reproducibility and corpus locators

This note documents where to find primary evidence in-repo; it does not add new scientific claims beyond the cited publication.

Normalized layer. When present, normalized/papers/{slug}.json mirrors manifest hashes, bibliography fields, and extraction pointers; if pdf_path or PDF bytes change, follow AGENTS.md and docs/PHASE3_RUNBOOK.md to re-profile rather than editing PDFs in place.

Authority chain. For numerical settings (cutoffs, timesteps, ensembles, kinetics), use the peer-reviewed PDF (and publisher Supporting Information) as the authoritative source; this wiki summarizes for navigation and retrieval.