pith. sign in

arxiv: 2606.04242 · v1 · pith:SYDKOQ5Znew · submitted 2026-06-02 · ⚛️ physics.chem-ph · quant-ph

Spin dynamics and ortho-para conversion in H₂O at the gas-ice phase transition in external magnetic fields

Chrysovalantis S. Kannis , Ralf Engels , Nicolas Faatz , Simon J. P\"utz , Markus B\"uscher This is my paper

Pith reviewed 2026-06-28 07:42 UTC · model grok-4.3

classification ⚛️ physics.chem-ph quant-ph
keywords spin dynamicsortho-para conversionwater icemagnetic fieldsnuclear spin polarizationgas-ice transitiondensity operator
0
0 comments X

The pith

External magnetic fields can convert initial para-H2O to over 90% ortho population right after the gas-to-ice transition.

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

The paper extends a four-spin model of nearest-neighbor water molecules to include coupling with external magnetic fields and tracks the evolution of ortho and para populations in the density-operator formalism during the first tens of milliseconds after freezing. Static homogeneous fields are shown to suppress dipolar depolarization and raise the ortho fraction to roughly 50% from a pure para starting state. Suitably chosen sinusoidal fields in relative motion with the molecules push the ortho population above 90%. These conversions matter for any scheme that tries to preserve or control nuclear-spin polarization while water molecules deposit as ice.

Core claim

By adding magnetic-field interaction terms to the four-spin nearest-neighbor Hamiltonian of two water molecules, the time evolution under static fields limits ortho population growth to about 50%, whereas time-dependent sinusoidal fields can drive it beyond 90% from an initial all-para condition, all within the density-operator description of the first tens of milliseconds after the gas-to-solid transition.

What carries the argument

The four-spin nearest-neighbor model extended by magnetic-field coupling terms and evolved in the density-operator formalism.

If this is right

  • Static homogeneous fields reduce the rate of dipolar-induced depolarization after deposition.
  • Time-dependent sinusoidal fields allow ortho populations to exceed those reachable by static fields.
  • The ortho/para ratio can be set during the deposition step itself rather than afterward.
  • Nuclear-spin polarization can be manipulated on the timescale of tens of milliseconds during ice formation.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same field-control approach might be tested on other molecular ices whose spin dynamics are also dominated by nearest-neighbor dipolar couplings.
  • If the conversion works in the lab, it could be used to prepare ice samples with defined ortho fractions for low-temperature NMR or neutron-scattering studies.
  • The model assumes homogeneous or uniformly sinusoidal fields; spatially inhomogeneous fields during real deposition remain an open extension.

Load-bearing premise

The four-spin nearest-neighbor description continues to capture the spin dynamics accurately immediately after the gas-to-solid phase transition once external magnetic fields are added.

What would settle it

An experiment that deposits water ice while applying the predicted sinusoidal field configuration and then measures an ortho population below 90% would falsify the central claim.

Figures

Figures reproduced from arXiv: 2606.04242 by Chrysovalantis S. Kannis, Markus B\"uscher, Nicolas Faatz, Ralf Engels, Simon J. P\"utz.

Figure 1
Figure 1. Figure 1: Distances between the four protons of two nearest-neighbor water molecules in [PITH_FULL_IMAGE:figures/full_fig_p007_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Time evolution of the populations ⟨Poo⟩, ⟨Ppp⟩, and ⟨Pop⟩ along with the average total spin projection components ⟨Mx⟩, ⟨My⟩, and ⟨Mz⟩ for magnetic field values of 0 mT (a), 1 mT (b), and 2 mT (c). All molecules are initially prepared in the para state. different, as explained in Ref. [7]. This is expected as the Hamiltonian in the minimal model involves the dipolar interactions between the neighbor￾ing mo… view at source ↗
Figure 3
Figure 3. Figure 3: Time evolution of the populations ⟨Poo⟩, ⟨Ppp⟩, and ⟨Pop⟩ along with the average total spin projection components ⟨Mx⟩, ⟨My⟩, and ⟨Mz⟩ for magnetic field values of 0 mT (a), 0.5 mT (b), and 1 mT (c). The initial state corresponds to full spin polarization along the z axis. field strengths. Figures 4a and 4b show these quantities assuming the para state is initially fully populated, for magnetic fields up t… view at source ↗
Figure 4
Figure 4. Figure 4: Populations ⟨Poo⟩, ⟨Ppp⟩, and ⟨Pop⟩ (a) and average total spin projection com￾ponents ⟨Mx⟩, ⟨My⟩, and ⟨Mz⟩ (b) at time 6T2 as functions of the magnetic field. All molecules are initially prepared in the para state. interactions are of comparable strength, namely for fields below about 2 mT, more pronounced effects appear. In particular, the ortho population reaches a maximum of 52.6% at 1.48 mT, exceeding … view at source ↗
Figure 5
Figure 5. Figure 5: Populations ⟨Poo⟩, ⟨Ppp⟩, and ⟨Pop⟩ (a) and average total spin projection com￾ponents ⟨Mx⟩, ⟨My⟩, and ⟨Mz⟩ (b) at time 6T2 as functions of the magnetic field. The initial state corresponds to full spin polarization along the z axis. depolarization for an initially nuclear-spin-polarized beam. At intermediate fields (< 2 mT), where the two interactions are comparable, the ortho pop￾ulation can increase to ∼… view at source ↗
Figure 6
Figure 6. Figure 6: Schematic illustration of an H2O pellet moving with constant velocity v through a longitudinal sinusoidal magnetic field. The wavelength λ of the field is assumed to be much larger than the pellet radius rp (λ ≫ rp). where B0 is the field amplitude and λ its spatial period. In the following, we restrict ourselves to the regime r ≪ λ, such that the longitudinal component Bz dominates and the radial contribu… view at source ↗
Figure 7
Figure 7. Figure 7: Magnetic field configurations in the pellet rest frame corresponding to frequencies [PITH_FULL_IMAGE:figures/full_fig_p019_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Populations ⟨Poo⟩, ⟨Ppp⟩, and ⟨Pop⟩ at time 6T2 as functions of the magnetic field amplitude B0. All molecules are initially prepared in the para state and are subjected to a sinusoidal magnetic field with f = 100 Hz applied at t = 0 (see Fig. 7a). Figures 8–10 present the spin-state populations at time 6T2 for the mag￾netic field pulses depicted in [PITH_FULL_IMAGE:figures/full_fig_p020_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Populations ⟨Poo⟩, ⟨Ppp⟩, and ⟨Pop⟩ at time 6T2 as functions of the magnetic field amplitude B0. All molecules are initially prepared in the para state and are subjected to a sinusoidal magnetic field with f = 200 Hz at t = 15 ms (see Fig. 7b). 〈Poo〉 〈Ppp〉 〈Pop〉 sum 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 0.0 0.2 0.4 0.6 0.8 1.0 B0 (mT) 〈 P 〉 (a) 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 0.0 0.2 0.4 0.6 0.8 1.0 B0 (mT) 〈 P … view at source ↗
Figure 10
Figure 10. Figure 10: Populations ⟨Poo⟩, ⟨Ppp⟩, and ⟨Pop⟩ at time 6T2 as functions of the magnetic field amplitude B0. All molecules are initially prepared in the para state and are subjected to a sinusoidal magnetic field with f = 300 Hz at t = 0 (see Fig. 7c). 83.9% (at 3.6 mT), 86.8% (at 3.1 mT), and 91.9% (at 2.1 mT), respectively. To illustrate the last case in more detail, [PITH_FULL_IMAGE:figures/full_fig_p021_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Time evolution of the populations ⟨Poo⟩, ⟨Ppp⟩, and ⟨Pop⟩ (a), along with the average total spin projection components ⟨Mx⟩, ⟨My⟩, and ⟨Mz⟩ (b) for the magnetic field shown in Fig. 7c with B0 = 2.1 mT. All molecules are initially prepared in the para state. efficiencies, as suggested by the nearly 92% conversion obtained in the last of the examined cases. Moreover, as discussed in Ref. [14], the sinusoida… view at source ↗
read the original abstract

The spin dynamics of water ice in the presence of external magnetic fields are investigated. The employed model builds upon the approach introduced by Buntkowsky et al. [Z. Phys. Chem. 222, 1049 (2008)], which considers two nearest-neighbor water molecules and yields a four-spin system, as the abundant oxygen isotope has zero nuclear spin. The model is extended to include coupling to external magnetic fields, allowing us to analyze the interplay between magnetic dipole-dipole interactions and magnetic field coupling. Two types of configurations are examined: (i) static, homogeneous fields, corresponding to a time-independent interaction, and (ii) spatially varying sinusoidal fields in relative motion with the molecules, leading to a time-dependent interaction. All computations are performed within the density operator formalism. The ortho/para populations and the total spin projections are evaluated during the first tens of milliseconds following the gas-to-solid phase transition. For static homogeneous fields, we show that increasing field strength suppresses dipolar-induced depolarization. Assuming that all molecules are initially in the para state, we show that static homogeneous fields can drive the ortho population up to approximately $50\%$, whereas suitably chosen sinusoidal-field configurations can increase it beyond $90\%$. These results are relevant for schemes aiming to preserve or manipulate nuclear-spin polarization during deposition.

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 / 0 minor

Summary. The manuscript extends the four-spin nearest-neighbor model of Buntkowsky et al. to incorporate external magnetic fields (static homogeneous and time-dependent sinusoidal) and employs the density-operator formalism to evolve ortho/para populations and spin projections over the first tens of milliseconds after the gas-to-solid transition. Assuming an initial all-para state, it reports that static fields suppress dipolar depolarization and drive the ortho population to approximately 50%, while suitably chosen sinusoidal configurations can push it above 90%. The results are framed as relevant to schemes for preserving or manipulating nuclear-spin polarization during deposition.

Significance. If the four-spin model remains quantitatively accurate once external-field terms are added and when the system is at the gas-ice interface, the concrete population predictions would supply a theoretical route to magnetic control of ortho-para conversion on the 10-100 ms timescale, with possible implications for hyperpolarization and spin-polarized ice experiments. The work supplies explicit numerical outputs from unitary evolution under an extended Hamiltonian, which could be directly confronted with measurement.

major comments (2)
  1. [Abstract] Abstract: the headline claims that static fields reach ~50% ortho and sinusoidal fields exceed 90% ortho rest on unitary evolution of the four-spin density operator; no explicit test is supplied that enlarging the spin cluster to include additional proton neighbors or adding relaxation channels (whose timescales overlap the reported window) leaves these values unchanged.
  2. [Abstract] Computations (as described in the abstract): the reported population numbers are presented without error bars, convergence checks with respect to integration step or Hilbert-space truncation, or quantitative comparison against measured ortho-para conversion rates in ice under comparable conditions.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments on the abstract and the computational presentation. We address each point below and will revise the manuscript accordingly to improve clarity and acknowledge model limitations.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the headline claims that static fields reach ~50% ortho and sinusoidal fields exceed 90% ortho rest on unitary evolution of the four-spin density operator; no explicit test is supplied that enlarging the spin cluster to include additional proton neighbors or adding relaxation channels (whose timescales overlap the reported window) leaves these values unchanged.

    Authors: The four-spin nearest-neighbor model is the established framework introduced by Buntkowsky et al. and adopted here to enable direct comparison with prior work. Enlarging the cluster to include additional protons would expand the Hilbert space from 16 to 64 dimensions or higher, rendering the density-operator propagation over tens of milliseconds computationally infeasible with current resources. Relaxation terms are omitted because the study targets coherent dipolar evolution in the early post-transition window; their inclusion would require separate phenomenological rates not constrained by the present data. We will add a dedicated paragraph in the Discussion section explicitly stating these limitations and outlining how larger-cluster or open-system extensions could be pursued in future work. revision: partial

  2. Referee: [Abstract] Computations (as described in the abstract): the reported population numbers are presented without error bars, convergence checks with respect to integration step or Hilbert-space truncation, or quantitative comparison against measured ortho-para conversion rates in ice under comparable conditions.

    Authors: The reported populations result from exact unitary evolution of the closed four-spin density operator; consequently no statistical error bars are applicable. Numerical convergence was verified by halving the integration time step (results agree to better than 0.1 %) and by confirming that the 16-dimensional space is fully spanned with no truncation. Direct quantitative comparison to experimental ortho-para rates is limited because existing measurements are predominantly at zero or low field and do not replicate the gas-ice deposition protocol under controlled external fields; the present numbers are therefore forward predictions. We will revise the abstract to note the convergence checks and insert a short comparison to literature zero-field conversion timescales in the Results section. revision: yes

Circularity Check

0 steps flagged

No circularity; results follow from time evolution of externally cited Hamiltonian extended by standard Zeeman/dipolar terms

full rationale

The derivation begins from the Buntkowsky et al. four-spin nearest-neighbor Hamiltonian (external citation, different authors), augments it with magnetic-field coupling terms via the density-operator formalism, and computes ortho/para populations by unitary time evolution over 10–100 ms. No parameters are fitted to the reported populations, no self-citation chain supports the central claim, and no ansatz or uniqueness theorem is imported from the present authors' prior work. The initial para-state assumption is explicit and the outputs are direct consequences of the extended Schrödinger dynamics rather than redefinitions of the inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the validity of the two-molecule four-spin truncation and on the assumption that the phase-transition dynamics can be captured by coherent evolution under the extended dipolar-plus-Zeeman Hamiltonian for the first tens of milliseconds.

axioms (2)
  • domain assumption Oxygen-16 has zero nuclear spin, reducing the system to four proton spins.
    Stated in the abstract as the basis for the four-spin model.
  • domain assumption The nearest-neighbor two-molecule truncation introduced by Buntkowsky et al. remains sufficient when external fields are added.
    The paper builds directly upon that model without additional justification in the abstract.

pith-pipeline@v0.9.1-grok · 5790 in / 1349 out tokens · 13881 ms · 2026-06-28T07:42:50.204519+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

24 extracted references · 15 canonical work pages

  1. [1]

    Leger, J

    A. Leger, J. Klein, S. de Cheveigne, C. Guinet, D. Defourneau, M. Belin, The3.1µmabsorptioninmolecularcloudsisprobablyduetoamorphous H2O ice., A & A 79 (1-2) (1979) 256–259. URLhttps://ui.adsabs.harvard.edu/abs/1979A&A....79..256L

    1979

  2. [2]

    H. A. Weaver, M. J. Mumma, H. P. Larson, Infrared investigation of water in comet P/Halley, A & A 187 (1-2) (1987) 411–418. URLhttps://ui.adsabs.harvard.edu/abs/1987A&A...187..411W 26

    1987

  3. [3]

    M. J. Mumma, H. A. Weaver, H. P. Larson, The ortho-para ratio of water vapor in comet P/ Halley, A & A 187 (1987) 419. URLhttps://ui.adsabs.harvard.edu/abs/1987A&A...187..419M

    1987

  4. [4]

    Cernicharo, R

    J. Cernicharo, R. Bachiller, E. Gonzalez-Alfonso, Water emission at 183 GHz from HH7-11 and other low-mass star-forming regions., A & A 305 (1996) L5. URLhttps://ui.adsabs.harvard.edu/abs/1996A&A...305L...5C

    1996

  5. [5]

    J. C. Raich, R. H. Good, Jr., Ortho-Para Transition in Molecular Hy- drogen., A & A 139 (1964) 1004.doi:10.1086/147835

    work page doi:10.1086/147835 1964

  6. [6]

    Dodelson, Relativistic treatment of ortho-para-H 2 transitions, J

    S. Dodelson, Relativistic treatment of ortho-para-H 2 transitions, J. Phys. B: At. Mol. Phys. 19 (18) (1986) 2871.doi:10.1088/0022-3700/ 19/18/017

    work page doi:10.1088/0022-3700/ 1986

  7. [7]

    Buntkowsky, H.-H

    G. Buntkowsky, H.-H. Limbach, B. Walaszek, A. Adamczyk, Y. Xu, H. Breitzke, A. Schweitzer, T. Gutmann, M. Wächtler, N. Amadeu, D. Tietze, B. Chaudret, Mechanisms of dipolar ortho/para-H2O conver- sion in ice, Z. Phys. Chem. 222 (7) (2008) 1049–1063.doi:10.1524/ zpch.2008.5359

    arXiv 2008

  8. [8]

    Chlan, K

    A. Miani, J. Tennyson, Can ortho–para transitions for water be ob- served?, J. Chem. Phys. 120 (6) (2004) 2732–2739.doi:10.1063/1. 1633261

    work page doi:10.1063/1 2004

  9. [9]

    C. S. Kannis, T. P. Rakitzis, Macroscopic production of highly nuclear- spin-polarized molecules from IR-excitation and photodissociation of molecular beams, Chem. Phys. Lett. 784 (2021) 139092.doi:10.1016/ j.cplett.2021.139092

    arXiv 2021

  10. [10]

    C. S. Kannis, T. P. Rakitzis, Production of spin-polarized molecular beams via microwave or infrared rotational excitation, Chem. Phys. 607 (2026) 113200.doi:10.1016/j.chemphys.2026.113200

    work page doi:10.1016/j.chemphys.2026.113200 2026

  11. [11]

    C.M.Edwards, D.Zhou, N.S.Sullivan, Unusuallow-temperatureeffects on the NMR line shapes in solid hydrogen, Phys. Rev. B 34 (1986) 6540– 6542.doi:10.1103/PhysRevB.34.6540. 27

    work page doi:10.1103/physrevb.34.6540 1986

  12. [12]

    Ciullo, R

    G. Ciullo, R. Engels, M. Büscher, A. Vasilyev (Eds.), Nuclear Fusion with Polarized Fuel, Springer Proceedings in Physics 187, Springer In- ternational Publishing, 2016.doi:10.1007/978-3-319-39471-8

    work page doi:10.1007/978-3-319-39471-8 2016

  13. [13]

    Baylor, A

    L. Baylor, A. Deur, N. Eidietis, W. Heidbrink, G. Jackson, J. Liu, M. Lowry, G. Miller, D. Pace, A. Sandorfi, S. Smith, S. Tafti, K. Wei, X. Wei, X. Zheng, Polarized fusion and potential in situ tests of fuel polarization survival in a tokamak plasma, Nucl. Fusion 63 (7) (2023) 076009.doi:10.1088/1741-4326/acc3ae

    work page doi:10.1088/1741-4326/acc3ae 2023

  14. [14]

    C.S.Kannis, R.Engels, T.El-Kordy, N.Faatz, S.J.Pütz, V.Verhoeven, T. P. Rakitzis, M. Büscher, Spin manipulation and nuclear polarization enhancement in particle beams with static magnetic fields, Phys. Rev. A 112 (2025) 012801.doi:10.1103/4nr6-xt7m

    work page doi:10.1103/4nr6-xt7m 2025

  15. [15]

    Coudert, Analysis of the rotational levels of water, J

    L. Coudert, Analysis of the rotational levels of water, J. Mol. Spectrosc. 154 (2) (1992) 427–442.doi:10.1016/0022-2852(92)90220-I

    work page doi:10.1016/0022-2852(92)90220-i 1992

  16. [16]

    Buntkowsky, personal communication (2025)

    G. Buntkowsky, personal communication (2025)

    2025

  17. [17]

    Buntkowsky, B

    G. Buntkowsky, B. Walaszek, A. Adamczyk, Y. Xu, H.-H. Limbach, B. Chaudret, Mechanism of nuclear spin initiated para-H2 to ortho- H2 conversion, Phys. Chem. Chem. Phys. 8 (2006) 1929–1935.doi: 10.1039/B601594H

    work page doi:10.1039/b601594h 2006

  18. [18]

    Buntkowsky, I

    G. Buntkowsky, I. Sack, H. H. Limbach, B. Kling, J. Fuhrhop, Struc- ture elucidation of amide bonds with dipolar chemical shift NMR spectroscopy, J. Phys. Chem. B 101 (51) (1997) 11265–11272.doi: 10.1021/jp971904i

    work page doi:10.1021/jp971904i 1997

  19. [19]

    Tiesinga, P

    E. Tiesinga, P. J. Mohr, D. B. Newell, B. N. Taylor, CODATA recom- mended values of the fundamental physical constants: 2018, Rev. Mod. Phys. 93 (2021) 025010.doi:10.1103/RevModPhys.93.025010

    work page doi:10.1103/revmodphys.93.025010 2018

  20. [20]

    Modig, B

    K. Modig, B. Halle, Proton magnetic shielding tensor in liquid wa- ter, J. Am. Chem. Soc. 124 (40) (2002) 12031–12041.doi:10.1021/ ja026981s

    2002

  21. [21]

    Sofikitis, T

    D. Sofikitis, T. P. Rakitzis, Mesoscopic production of hyperpolarized 15N2OandH 2Ovia optical excitation, Phys. Rev. A 92 (2015) 032507. doi:10.1103/PhysRevA.92.032507. 28

    work page doi:10.1103/physreva.92.032507 2015

  22. [22]

    Reistad, B

    D. Reistad, B. Gålnander, T. Lofnes, Y.-N. Rao, Experiences of oper- ating CELSIUS with a hydrogen pellet target, Nucl. Instrum. Methods Phys. Res. Sect. A 532 (1) (2004) 118–122.doi:10.1016/j.nima.2004. 06.037

    work page doi:10.1016/j.nima.2004 2004

  23. [23]

    Abragam, The Principles of Nuclear Magnetism, 1st Edition, Oxford at Clarendon Press, Oxford, England, 1961

    A. Abragam, The Principles of Nuclear Magnetism, 1st Edition, Oxford at Clarendon Press, Oxford, England, 1961

    1961

  24. [24]

    Buntkowsky, H.-H

    G. Buntkowsky, H.-H. Limbach, H-solid state NMR studies of tunneling phenomena and isotope effects in transition metal dihydrides, J. Low Temp. Phys. 143 (3) (2006) 55–114.doi:10.1007/s10909-006-9211-y. 29

    work page doi:10.1007/s10909-006-9211-y 2006