pith. sign in

arxiv: 2605.01045 · v1 · submitted 2026-05-01 · ❄️ cond-mat.mtrl-sci

Understanding the lifetime of water with dynamic network analysis: the case of CsOH.H2O

Graeme J. Ackland , Ciprian G. Pruteanu , John S. Loveday , Keishiro Yamashita This is my paper

Pith reviewed 2026-05-09 18:50 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci
keywords caesium hydroxide monohydrateproton exchangehydrogen bond interconversionRaman spectroscopyionic conductionorder-disorder transitionwater lifetimedynamic network analysis
0
0 comments X

The pith

Water and hydroxide ions in CsOH·H2O continually interconvert by exchanging protons rather than rotating or diffusing as whole molecules.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper uses dynamic network analysis on simulations of caesium hydroxide monohydrate to track how the atoms move in its layered structure. Oxygen atoms sit in a honeycomb lattice where each participates in three hydrogen bonds, but the covalent and hydrogen bonds switch identities through ongoing proton transfers. This means individual water molecules and hydroxide ions exchange their identities on short timescales. The order-disorder transition therefore happens by this chemical reaction, and the same process creates a combined stretch-plus-exchange Raman signal that the authors model. The result also accounts for rapid vacancy diffusion that produces fast ionic conduction.

Core claim

In CsOH·H2O the oxygen atoms form a quasi-2D honeycomb lattice held by three hydrogen bonds per oxygen. Covalent and hydrogen bonds continually interconvert through proton exchange, so water molecules and hydroxide ions switch identities. The order-disorder transition therefore occurs by this chemical reaction instead of molecular rotation or diffusion. Proton exchange also creates a new Raman activity involving both stretch and exchange processes, which produces a broad peak covering H2O and OH stretches and a low-frequency peak that appears at higher temperature.

What carries the argument

Dynamic network analysis of continual covalent-hydrogen bond interconversion via proton exchange in the quasi-2D oxygen honeycomb lattice.

Load-bearing premise

The molecular dynamics simulations and the network analysis built on them accurately reflect the real atomic motions without artifacts from the force field or system size.

What would settle it

Raman spectra at low temperature that show two distinct narrow peaks for separate H2O and OH stretches, instead of one broad combined peak, would falsify the continual proton-exchange claim.

Figures

Figures reproduced from arXiv: 2605.01045 by Ciprian G. Pruteanu, Graeme J. Ackland, John S. Loveday, Keishiro Yamashita.

Figure 1
Figure 1. Figure 1: FIG. 1. Illustrating the vicinal molecular dynamics cell. Top view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Unstable I4 view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Top left: The arrangement of OH and H view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Partial RDFs at temperatures 200K, 300K, 400K, view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Anharmonic potential model for H view at source ↗
read the original abstract

We describe the atomic-level motions in caesium hydroxide monohydrate (CsOH$\cdot$H$_2$O), which is a chemical compound containing layers of water and hydroxide ions. At this composition, each oxygen is involved in three hydrogen bonds which, in the hexagonal structure, form a quasi-2D honeycomb lattice. While oxygen and caesium atoms form a typical crystal lattice, the dynamics of the hydrogen atoms are more complex. Here we show that the covalent and hydrogen bonds are continually interconverting, meaning that the water and hydroxyl are interconverting by proton exchange. The order-disorder transition of the water and hydroxyl proceeds by chemical reaction rather than rotation or diffusion of the molecules. A hydrogen can rotate out of the layer, leaving a vacant site in the 2D layer. Such a hydrogen vacancy can diffuse rapidly by single molecule rotation, leading to fast-ionic conduction. The proton exchange leads to a novel type of Raman activity combining stretch and exchange processes, for which we develop a theoretical model. This would manifest in a broad single peak associated with both H$_2$O and OH stretches and a low frequency peak appearing at elevated temperature.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript applies dynamic network analysis to CsOH·H2O, a layered hydroxide monohydrate with a quasi-2D honeycomb hydrogen-bond lattice. It claims that covalent O–H bonds and hydrogen bonds continually interconvert, so that water molecules and hydroxide ions exchange protons; the order-disorder transition therefore occurs by chemical reaction rather than molecular rotation or diffusion. Vacancy diffusion by single-molecule rotation is said to produce fast ionic conduction. A theoretical model is developed for a novel Raman process that combines stretch and exchange motions, predicted to yield a broad single peak encompassing both H2O and OH stretches plus a low-frequency feature that appears at elevated temperature.

Significance. If the central claims are substantiated, the work would supply a concrete mechanistic picture of proton exchange in a model 2D hydrogen-bonded system and a spectroscopic signature that could be tested experimentally. The combination of network-based lifetime analysis with an explicit Raman model is a potentially reusable approach for other fast-ion conductors or disordered hydrogen-bond networks.

major comments (2)
  1. [Abstract and Methods] The central claim that the order-disorder transition proceeds by proton exchange (rather than rotation or diffusion) rests on the dynamic network analysis correctly identifying genuine covalent-bond breaking/forming events. The abstract supplies no description of the underlying trajectories (ab initio MD, reactive force field, or classical non-reactive MD), the precise network definition (distance/angle cutoffs, lifetime thresholds), or any validation that observed interconversions are not analysis artifacts arising from librations or H-bond rearrangements. This information is load-bearing for the chemical-reaction interpretation.
  2. [Theoretical model section] The theoretical Raman model is presented as combining stretch and exchange processes and yielding a broad single peak plus a low-frequency feature. No equations, parameter values, or derivation steps are referenced in the abstract, so it is impossible to assess whether the predicted lineshape is a direct consequence of the network analysis or contains additional assumptions that could alter the temperature dependence.
minor comments (2)
  1. [Abstract] The abstract states qualitative conclusions without any numerical values, error bars, or comparison to experimental Raman spectra or conductivity data; adding at least one quantitative result (e.g., exchange lifetime or diffusion coefficient) would strengthen the presentation.
  2. [Methods] Notation for the hydrogen-bond network (e.g., how nodes and edges are defined) should be introduced with a short diagram or explicit equations in the main text rather than left implicit.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments and the opportunity to clarify key aspects of our work. We address the major comments point by point below.

read point-by-point responses

  1. Referee: [Abstract and Methods] The central claim that the order-disorder transition proceeds by proton exchange (rather than rotation or diffusion) rests on the dynamic network analysis correctly identifying genuine covalent-bond breaking/forming events. The abstract supplies no description of the underlying trajectories (ab initio MD, reactive force field, or classical non-reactive MD), the precise network definition (distance/angle cutoffs, lifetime thresholds), or any validation that observed interconversions are not analysis artifacts arising from librations or H-bond rearrangements. This information is load-bearing for the chemical-reaction interpretation.

    Authors: We agree that the abstract would benefit from a concise description of the computational approach to make the chemical-reaction interpretation immediately accessible. The trajectories are ab initio molecular dynamics simulations. The dynamic network analysis uses distance and angle cutoffs to identify covalent versus hydrogen bonds together with lifetime thresholds to detect genuine interconversions. We have validated that these events correspond to proton transfers rather than librational artifacts through direct inspection of bond-length time series and trajectory visualizations. We will revise the abstract to include a brief statement on the simulation method, network definition, and validation procedure, and we will expand the Methods section with the exact parameter values and supporting checks. revision: yes

  2. Referee: [Theoretical model section] The theoretical Raman model is presented as combining stretch and exchange processes and yielding a broad single peak plus a low-frequency feature. No equations, parameter values, or derivation steps are referenced in the abstract, so it is impossible to assess whether the predicted lineshape is a direct consequence of the network analysis or contains additional assumptions that could alter the temperature dependence.

    Authors: The Raman model is derived in the dedicated Theoretical model section, where the stretch frequencies are combined with the proton-exchange rates obtained directly from the network analysis. The broad single peak results from exchange-induced averaging of the H2O and OH stretches, while the low-frequency feature arises from the temperature dependence of the exchange process. We will add a short sentence to the abstract that references the key equations and notes that the lineshape follows from the network-derived exchange rates, allowing readers to evaluate the temperature dependence without consulting the full text. revision: yes

Circularity Check

0 steps flagged

No significant circularity; claims derive from independent network analysis of trajectories

full rationale

The paper's central claims rest on dynamic network analysis applied to (presumably ab initio or reactive) molecular dynamics trajectories of CsOH·H2O. The interconversion of covalent and hydrogen bonds, the proton-exchange mechanism for the order-disorder transition, and the derived Raman model are presented as direct outputs of that analysis rather than as quantities fitted to the same data or defined in terms of themselves. No self-citations are invoked as load-bearing uniqueness theorems, no parameters are fitted on a subset and then 'predicted' on the remainder, and no ansatz is smuggled via prior work. The derivation chain therefore remains self-contained against external benchmarks (simulation trajectories and network definitions) and does not reduce to its inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Only the abstract is available, so the ledger is necessarily incomplete; no explicit free parameters, invented entities, or non-standard axioms are stated beyond the implicit use of standard condensed-matter simulation assumptions.

axioms (1)
  • domain assumption Standard assumptions underlying molecular dynamics or network analysis of hydrogen-bonded systems hold for this material.
    The description of bond interconversion and vacancy diffusion implies reliance on computational modeling whose validity is not detailed.

pith-pipeline@v0.9.0 · 5525 in / 1359 out tokens · 54150 ms · 2026-05-09T18:50:50.169426+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

28 extracted references · 28 canonical work pages

  1. [1]

    In this structure each O has three neighbouring O with a hydrogen between

    Tetragonal structure - 3D network We set up a unit cell with theI4 1/amdsymmetry[6] (fig.2a) and relaxed it (Table.III). In this structure each O has three neighbouring O with a hydrogen between. The formal symmetry places each hydrogen midway be- tween two oxygen atoms, but in calculations this is highly unstable with respect to molecule formation. There...

    work page

  2. [2]

    wrong bonds

    Hexagonal structure - 2D network An alternative set of structures are based on close- packed layers of Cs and oxygen occupying the interstitial sites. It is helpful to consider first the high symmetry structureP6/mmm. If one ignores the hydrogens this is strukturbericht C32 (AlB 2-type). The unit cell has just one formula unit: the Cs is arranged in a tri...

    work page

  3. [3]

    We aver- aged over a credible range of barrier parameters forV(x) FIG

    The second excited state is centred on the other potential well, and has ap- proximately the same energy for both H and D. We aver- aged over a credible range of barrier parameters forV(x) FIG. 1. Illustrating the vicinal molecular dynamics cell. Top figure shows all 60/120 independent O/H atoms as a single layer. Lower figure shows side view of the layer...

    work page

  4. [4]

    Shapley, Of stars and men; the human response to an expanding universe., Boston (1964)

    H. Shapley, Of stars and men; the human response to an expanding universe., Boston (1964)

    work page 1964

  5. [5]

    Given those 1047 molecules on earth, if the hydrogen mixes the chance of a given oxygen having the same two hydrogens is 10−94

    Typically about 10 25 molecules in the cup, 10 22 cups in the oceans, so the standard answer is about 10 3. Given those 1047 molecules on earth, if the hydrogen mixes the chance of a given oxygen having the same two hydrogens is 10−94. So that means of order 10 −69 molecules in both cups. Which, OK, isn’t exactly zero

    work page

  6. [6]

    Hermann, High-pressure phase transitions in rubidium and caesium hydroxides, Phys

    A. Hermann, High-pressure phase transitions in rubidium and caesium hydroxides, Phys. Chem. Chem. Phys.18, 16527 (2016)

    work page 2016

  7. [7]

    Jacobs, B

    H. Jacobs, B. Harbrecht, P. M¨ uller, and W. Bronger, Struktur und eigenschaften von caesiumhydroxid- monohydrat–einer verbindung, die in ihrer hochtemperatur-form schichtenf¨ ormige [h3o2-]- polyanionen enth¨ alt [1, 2], Zeitschrift f¨ ur anorganische und allgemeine Chemie491, 154 (1982)

    work page 1982

  8. [8]

    R. Marx, H. Dachs, and R. Ibberson, Time-of-flight neu- tron diffraction studies on the different phases of csoh, h2o, The Journal of chemical physics93, 5972 (1990)

    work page 1990

  9. [9]

    ˇCern` y, V

    R. ˇCern` y, V. Favre-Nicolin, and B. Bertheville, A tetrag- onal polymorph of caesium hydroxide monohydrate, csoh·h2o, from x-ray powder data, Crystal Structure Communications58, i31 (2002)

    work page 2002

  10. [10]

    H. P. Rodenburg, F. Stainer, K. M. Draijer, H. Ni, J. Spy- chala, N. Artrith, H. M. R. Wilkening, and P. Ngene, Superprotonic conductivity in hexagonal and tetragonal cesium hydroxide hydrate, Advanced Functional Materi- als35, 2412219 (2025)

    work page 2025

  11. [11]

    H. Lutz, J. Henning, H. Jacobs, and B. Harbrecht, Hy- drogen Bonding and Phase Transitions of RbOH·H2O and CsOH·H2O Studied by IR and Raman Spectroscopy, Berichte der Bunsengesellschaft f¨ ur Phys. Chemie92, 1557 (1988)

    work page 1988

  12. [12]

    Stahn, R

    M. Stahn, R. E. Lechner, H. Dachs, and H. E. Jacobs, Dy- namics of the hydrogen bond in CsH3O2 from quasielas- tic neutron scattering, J. Phys. C Solid State Phys.16, 5073 (1983)

    work page 1983

  13. [13]

    S. J. Clark, M. D. Segall, C. J. Pickard, P. J. Hasnip, M. J. Probert, K. Refson, and M. C. Payne, First princi- ples methods using castep, Zeitschrift f¨ ur kristallographie , 567 (2005)

    work page 2005

  14. [14]

    A. P. Bart´ ok and J. R. Yates, Regularized scan func- tional, The Journal of chemical physics150(2019)

    work page 2019

  15. [15]

    J. P. Perdew, K. Burke, and M. Ernzerhof, Generalized gradient approximation made simple, Physical review let- ters77, 3865 (1996)

    work page 1996

  16. [16]

    CASTEP eliminates the divergent dipole-dipole energy using an Ewald sum

    work page

  17. [17]

    C. H. Loach and G. J. Ackland, Stacking characteristics of close packed materials, Physical Review Letters119, 205701 (2017)

    work page 2017

  18. [18]

    Whalley, J

    E. Whalley, J. Heath, and D. Davidson, Ice ix: an antifer- roelectric phase related to ice iii, The Journal of Chemical Physics48, 2362 (1968)

    work page 1968

  19. [19]

    W. L. Jorgensen, J. Chandrasekhar, J. D. Madura, R. W. 8 Impey, and M. L. Klein, Comparison of simple poten- tial functions for simulating liquid water, The Journal of chemical physics79, 926 (1983)

    work page 1983

  20. [20]

    Kuhs and M

    W. Kuhs and M. Lehmann, The structure of ice-ih, Water Sci. Rev2(1986)

    work page 1986

  21. [21]

    W. F. Kuhs and M. Lehmann, The geometry and ori- entation of the water molecule in ice ih, Le Journal de Physique Colloques48, C1 (1987)

    work page 1987

  22. [22]

    Haus and T

    J. Haus and T. Tanaka, Model for the ice VII-ice VIII transition, Physical Review B16, 2148 (1977)

    work page 1977

  23. [23]

    G. J. Ackland, A. Hermann, K. Komatsu, K. Yamashita, and J. Loveday, Distinction between ice phases vii, viii, and x, Physical Review B112, 094116 (2025)

    work page 2025

  24. [24]

    Fletcher,Practical methods of optimization(John Wi- ley & Sons, 2013)

    R. Fletcher,Practical methods of optimization(John Wi- ley & Sons, 2013)

    work page 2013

  25. [25]

    M. F. Sykes and J. W. Essam, Exact critical percolation probabilities for site and bond problems in two dimen- sions, Journal of Mathematical Physics5, 1117 (1964)

    work page 1964

  26. [26]

    Tajima, T

    Y. Tajima, T. Matsuo, and H. Suga, Phase transition in koh-doped hexagonal ice, Nature299, 810–812 (1982)

    work page 1982

  27. [27]

    C. G. Salzmann, P. G. Radaelli, A. Hallbrucker, E. Mayer, and J. L. Finney, The preparation and struc- tures of hydrogen ordered phases of ice, Science311, 1758 (2006)

    work page 2006

  28. [28]

    A. Gao, G. Li, B. Peng, Y. Xie, and H. F. Schaefer, The symmetric exchange reaction OH+ H2O→H2O+ OH: convergent quantum mechanical predictions, The Journal of Physical Chemistry A120, 10223 (2016)

    work page 2016