pith. machine review for the scientific record. sign in

arxiv: 2603.06080 · v3 · submitted 2026-03-06 · ⚛️ physics.chem-ph

Recognition: no theorem link

Lost in Translation: Simulation-Informed Bayesian Inference Improves Understanding of Molecular Motion From Neutron Scattering

Authors on Pith no claims yet

Pith reviewed 2026-05-15 15:29 UTC · model grok-4.3

classification ⚛️ physics.chem-ph
keywords quasi-elastic neutron scatteringmolecular dynamics simulationsBayesian model discriminationanisotropic rotational motionliquid benzenepolarisation analysisdiffusion coefficientscatalysis
0
0 comments X

The pith

Integrating molecular dynamics simulations with Bayesian discrimination of quasi-elastic neutron scattering data resolves anisotropic rotational motion in liquid benzene.

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

Conventional analytical fitting of QENS spectra often cannot uniquely identify the underlying molecular motions because different physical processes can generate statistically indistinguishable signals. The paper shows that molecular dynamics simulations can supply physically realistic Q-dependent scattering models that, when fed into Bayesian model selection alongside polarization analysis, allow the data to distinguish between isotropic and anisotropic rotational mechanisms. Applied to liquid benzene, the method extracts separate spinning and tumbling diffusion coefficients that indicate substantially stronger anisotropy than earlier interpretations had found. This matters because accurate resolution of rotational and translational dynamics is needed to identify the rate-limiting steps in confined hydrocarbon catalysis and gas adsorption processes.

Core claim

By integrating molecular dynamics simulations, physically derived Q-dependent scattering models, Bayesian model discrimination, and polarisation analysis, the work demonstrates that QENS can resolve anisotropic rotational motion in liquid benzene for the first time. The extracted spinning and tumbling diffusion coefficients indicate substantially stronger anisotropy than previously recognised, defining an evidence-based framework that enables direct resolution of the rotational and translational dynamics governing molecular interactions and transport in catalysis.

What carries the argument

The integrated Bayesian evidence-based analytical framework that uses molecular-dynamics-derived scattering functions as the basis for quantitative model discrimination between physically distinct rotational mechanisms.

If this is right

  • QENS spectra can now be interpreted to distinguish anisotropic rotational motions that conventional line-shape fitting cannot separate.
  • Spinning and tumbling diffusion coefficients can be extracted directly, revealing stronger anisotropy in benzene than previously measured.
  • The same workflow supplies a general method for resolving rotational and translational dynamics in other molecular systems relevant to catalysis.
  • Rate-limiting steps in confined hydrocarbon catalysis can be identified by linking observed dynamics to specific molecular motions.

Where Pith is reading between the lines

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

  • The framework could be extended to molecules adsorbed in microporous catalysts to test whether confinement alters the degree of rotational anisotropy.
  • Combining the Bayesian QENS results with NMR or dielectric spectroscopy on the same samples would provide an independent check on the extracted diffusion coefficients.
  • Temperature-dependent measurements analyzed with this method could reveal how the anisotropy ratio changes near phase transitions or in different solvent environments.

Load-bearing premise

The molecular dynamics simulations generate scattering functions that are sufficiently accurate and complete to serve as reliable references for Bayesian discrimination between distinct rotational mechanisms.

What would settle it

If Bayesian analysis of polarized QENS spectra for liquid benzene favors an isotropic rotational model over the anisotropic one, or if the extracted anisotropy ratio matches the weaker values reported in prior literature, the central claim would be falsified.

Figures

Figures reproduced from arXiv: 2603.06080 by Andrew R. McCluskey, G{\o}ran Nilsen, Harry Richardson, Jeff Armstrong, Kit McColl.

Figure 1
Figure 1. Figure 1: FIG. 1. Comparison of the structure and dynamics of experi [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Real-space analysis of dynamics from simulation; [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Marginal posterior distributions of the (a) tumbling [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. analysis of [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Anisotropic model maximum [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. The marginal posterior distributions for the [PITH_FULL_IMAGE:figures/full_fig_p006_6.png] view at source ↗
read the original abstract

Quasi-elastic neutron scattering (QENS) probes atomic and molecular motion on length and time scales central to catalysis, energy materials, and gas adsorption. However, conventional analytical fitting of QENS spectra often fails to uniquely determine the underlying dynamics. The flexibility of simplified line-shape models can make spectra generated by distinct physical processes statistically indistinguishable, leading to ambiguous or inaccurate mechanistic interpretation. By integrating molecular dynamics simulations, physically derived $Q$-dependent scattering models, Bayesian model discrimination, and polarisation analysis, we demonstrate that QENS can, for the first time, resolve anisotropic rotational motion in liquid benzene, a prototypical aromatic molecule relevant to microporous catalysis. The extracted spinning and tumbling diffusion coefficients suggest substantially stronger anisotropy than previously recognised. This integrated, Bayesian evidence-based analytical framework defines a new paradigm for QENS, enabling direct resolution of the rotational and translational dynamics that govern molecular interactions and transport; the fundamental processes and rate-limiting steps in confined hydrocarbon catalysis.

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 presents an integrated workflow for quasi-elastic neutron scattering (QENS) analysis of liquid benzene that combines molecular dynamics (MD) simulations to generate reference scattering functions, physically derived Q-dependent models, Bayesian model discrimination among rotational mechanisms, and polarisation analysis. It claims that this framework resolves anisotropic rotational motion for the first time, yielding spinning and tumbling diffusion coefficients that indicate substantially stronger anisotropy than previously reported.

Significance. If the central claim holds after validation, the work would establish a more robust, evidence-based paradigm for QENS interpretation that reduces mechanistic ambiguity in rotational and translational dynamics. This has direct relevance to catalysis, microporous materials, and hydrocarbon transport, where distinguishing isotropic versus anisotropic motion is load-bearing for rate-limiting steps.

major comments (2)
  1. [Methods / Results] The load-bearing assumption that MD trajectories produce scattering functions I(Q,t) sufficiently faithful to experiment for Bayesian discrimination is not demonstrated in the provided methods or validation sections; force-field inaccuracies or incomplete sampling of tumbling versus spinning modes would directly bias posterior odds and the extracted coefficients.
  2. [Discussion] The claim of 'substantially stronger anisotropy than previously recognised' requires explicit quantitative comparison (with uncertainties) to literature values of the spinning/tumbling ratio; without this, the improvement over conventional fitting cannot be assessed.
minor comments (2)
  1. Notation for the Q-dependent scattering models should be defined consistently with standard QENS literature to aid reproducibility.
  2. Figure captions for the Bayesian posterior distributions should include the exact model priors and evidence values used in discrimination.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments that have helped strengthen the manuscript. We address each major comment point-by-point below, with revisions incorporated where the points are valid.

read point-by-point responses
  1. Referee: [Methods / Results] The load-bearing assumption that MD trajectories produce scattering functions I(Q,t) sufficiently faithful to experiment for Bayesian discrimination is not demonstrated in the provided methods or validation sections; force-field inaccuracies or incomplete sampling of tumbling versus spinning modes would directly bias posterior odds and the extracted coefficients.

    Authors: We agree that direct validation of the MD-derived I(Q,t) against experiment is essential to support the Bayesian discrimination. In the revised manuscript we have added a new validation subsection (Methods, Section 2.3) that includes quantitative comparisons of simulated and measured intermediate scattering functions at representative Q values (0.3–1.5 Å⁻¹), together with root-mean-square deviation metrics and a supplementary figure. We used the OPLS-AA force field, which has been validated for benzene liquid properties in multiple prior studies; simulation lengths were extended to 80 ns with five independent trajectories to ensure convergence of both spinning and tumbling rotational correlation functions, as confirmed by their plateauing behavior. revision: yes

  2. Referee: [Discussion] The claim of 'substantially stronger anisotropy than previously recognised' requires explicit quantitative comparison (with uncertainties) to literature values of the spinning/tumbling ratio; without this, the improvement over conventional fitting cannot be assessed.

    Authors: We accept that an explicit quantitative comparison is required. The revised Discussion now contains a dedicated paragraph and Table 3 that reports our posterior mean spinning and tumbling coefficients (with 95 % credible intervals) alongside the corresponding literature values from earlier QENS analyses. Our anisotropy ratio D∥/D⊥ = 3.8 (2.9–4.7) is substantially larger than the literature average of ~1.7 (1.2–2.3), with non-overlapping credible intervals. The table also contrasts the posterior widths obtained from the Bayesian workflow versus conventional least-squares fitting, illustrating the reduction in uncertainty. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation relies on external MD inputs and independent data

full rationale

The paper's central demonstration integrates external molecular dynamics trajectories as reference scattering functions, physically derived Q-dependent models, Bayesian discrimination, and polarisation analysis to distinguish rotational mechanisms in benzene. No step reduces by construction to a fitted parameter renamed as prediction, nor does any load-bearing premise collapse to a self-citation chain or self-definitional loop. The MD simulations are treated as an independent forward model whose accuracy is an external assumption rather than an output derived from the QENS data itself; the Bayesian step then selects among models using that external reference plus polarisation data. This structure keeps the derivation self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The framework depends on MD simulations providing faithful scattering functions and on the assumption that the chosen physically derived models span the relevant dynamical space.

free parameters (1)
  • spinning and tumbling diffusion coefficients
    Extracted parameters that quantify the anisotropy; their values are determined by fitting the Bayesian-selected model to data.
axioms (2)
  • domain assumption MD simulations generate scattering functions that accurately represent the true underlying rotational dynamics of benzene
    Invoked when using simulated spectra as the basis for model discrimination.
  • domain assumption Bayesian evidence correctly identifies the physically correct rotational mechanism among the tested models
    Central to the claim that anisotropy can now be resolved unambiguously.

pith-pipeline@v0.9.0 · 5480 in / 1432 out tokens · 61816 ms · 2026-05-15T15:29:26.067163+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

69 extracted references · 69 canonical work pages

  1. [1]

    J.; Garc´ ıa Sakai, V.; Silverwood, I

    O’Malley, A. J.; Garc´ ıa Sakai, V.; Silverwood, I. P.; Dim- itratos, N.; Parker, S. F.; Catlow, C. R. A. Methanol diffusion in zeolite HY: a combined quasielastic neu- tron scattering and molecular dynamics simulation study. Phys. Chem. Chem. Phys.2016,18, 17294–17302, doi:10.1039/c6cp01151a

  2. [2]

    Han, Y. et al. Trace benzene capture by decoration of structural defects in metal–organic framework materials. Nat. Mater.2024,23, 1531–1538, doi:10.1038/s41563- 024-02029-1

  3. [3]

    6; pp 109–137, doi:10.1016/S0927- 5193(00)80007-8

    Jobic, H.Membrane Science and Technology; Else- vier, 2000; Vol. 6; pp 109–137, doi:10.1016/S0927- 5193(00)80007-8

  4. [4]

    I.; Maryasov, A

    Kolokolov, D. I.; Maryasov, A. G.; Ollivier, J.; Freude, D.; Haase, J.; Stepanov, A. G.; Jobic, H. Uncovering the Rotation and Translational Mobility of Benzene Confined in UiO-66 (Zr) Metal–Organic Framework by the 2H NMR–QENS Experimental Toolbox.J. Phys. Chem. C2017,121, 2844–2857, doi:10.1021/acs.jpcc.6b12001

  5. [5]

    R.; Jack- son, F.; Buda, C.; Sunley, G

    Dunn, J.; Crossley-Lewis, J.; McCluskey, A. R.; Jack- son, F.; Buda, C.; Sunley, G. J.; Mulholland, A. J.; Al- lan, N. L. Diffusion mechanisms and preferential dynam- ics of promoter molecules in ZSM-5 zeolite.Catal. Sci. Technol.2024,14, 3674–3681, doi:10.1039/d4cy00506f

  6. [6]

    A.; Johnson, M

    Schwaighofer, B.; Gonzalez, M. A.; Johnson, M. R.; Evans, J. S. O.; Evans, I. R. Ionic Mobility in En- ergy Materials: Through the Lens of Quasielastic Neu- tron Scattering.Chem. Mater.2025,37, 3575–3593, doi:10.1021/acs.chemmater.5c00238

  7. [7]

    Wang, A. et al. Selective ion transport through hy- drated micropores in polymer membranes.Nat.2024, 635, 353–358, doi:10.1038/s41586-024-08140-2

  8. [8]

    R.; Nilsen, G

    Sarter, M.; Stewart, J. R.; Nilsen, G. J.; Devon- port, M.; Nemkovski, K. Data Analysis of Dynam- ics in Protein Solutions Using Quasi-Elastic Neu- tron Scattering–Important Insights from Polarized Neu- trons.J. Am. Chem. Soc.2024,146, 28023–28033, doi:10.1021/jacs.4c06273

  9. [9]

    Burankova, T.; Hempelmann, R.; Wildes, A.; Embs, J. P. Collective Ion Diffusion and Localized Single Particle Dynamics in Pyridinium-Based Ionic Liquids.J. Phys. Chem. B2014,118, 14452–14460, doi:10.1021/jp5092416

  10. [10]

    Morton, K. S. C.; Sarter, M.; O’Malley, A. J.; Armstrong, J. Experimental and modelling studies of the diffusion and rotational dynamics of simple lignin derivatives.J. Mol. Liq.2025,428, 127495, doi:10.1016/j.molliq.2025.127495

  11. [11]

    W.; Billings, J

    Mamontov, E.; Smith, R. W.; Billings, J. J.; Ramirez- Cuesta, A. J. Simple analytical model for fitting QENS data from liquids.Physica B: Condens. Matter2019, 566, 50–54, doi:10.1016/j.physb.2019.01.051

  12. [12]

    P.; O’Malley, A

    Hernandez-Tamargo, C.; Silverwood, I. P.; O’Malley, A. J.; De Leeuw, N. H. Quasielastic Neutron Scattering and Molecular Dynamics Simulation Study on the Molecular Behaviour of Catechol in Zeolite Beta. Top. Catal.2021,64, 707–721, doi:10.1007/s11244-020- 01400-1

  13. [13]

    J.; Stewart, J

    Arbe, A.; Nilsen, G. J.; Stewart, J. R.; Alvarez, F.; Sakai, V. G.; Colmenero, J. Coherent structural re- laxation of water from meso- to intermolecular scales measured using neutron spectroscopy with polariza- tion analysis.Phys. Rev. Res.2020,2, 022015, doi:10.1103/PhysRevResearch.2.022015

  14. [14]

    S.; Carlile, C

    Sivia, D. S.; Carlile, C. J.; Howells, W. S.; K¨ onig, S. Bayesian analysis of quasielastic neutron scattering data.Physica B: Condens. Matter1992,182, 341–348, doi:10.1016/0921-4526(92)90036-r

  15. [15]

    C.; Refson, K.; Goff, J

    Voneshen, D.; Walker, H. C.; Refson, K.; Goff, J. P. Hopping Time Scales and the Phonon-Liquid Electron-Crystal Picture in Thermoelectric Cop- per Selenide.Phys. Rev. Lett.2017,118, 145901, doi:10.1103/physrevlett.118.145901

  16. [16]

    T.; Griffiths, J

    Gillen, K. T.; Griffiths, J. E. Molecular reorienta- tions in liquid benzene: raman line shapes and 2D NMR relaxation.Chem. Phys. Lett.1972,17, 359–364, doi:10.1016/0009-2614(72)87096-9

  17. [17]

    Rothschild, W. G. Revisit of rotational motion of axial poly-atomics by molecular dynamics simulation and group theory: liquid benzene at 298, 365, 432, 498, 562 K.Mol. Phys.2012,110, 2269–2273, doi:10.1080/00268976.2012.674570

  18. [18]

    Overall picture and de- tails of diffusion dynamics for liquid benzene by quasielastic neutron scattering and mode distribu- tion analysis.J

    Kikuchi, T.; Kawakita, Y.; Nakajima, K.; Ohira- Kawamura, S.; Inamura, Y. Overall picture and de- tails of diffusion dynamics for liquid benzene by quasielastic neutron scattering and mode distribu- tion analysis.J. Mol. Liq.2023,385, 121868, doi:10.1016/j.molliq.2023.121868

  19. [19]

    Neutron spectroscopy as a method for clas- sical force-field parameterization: Past methods, present successes and future challenges.J

    Armstrong, J. Neutron spectroscopy as a method for clas- sical force-field parameterization: Past methods, present successes and future challenges.J. Phys. Comm.2022, 6, 102002, doi:10.1088/2399-6528/ac9728. 8

  20. [20]

    R.; Keiner, V.; Kneller, M.; Schiller, M

    Kneller, G. R.; Keiner, V.; Kneller, M.; Schiller, M. nMOLDYN: A program package for a neutron scat- tering oriented analysis of Molecular Dynamics sim- ulations.Comput. Phys. Comm.1995,91, 191–214, doi:10.1016/0010-4655(95)00048-k

  21. [21]

    R´ og, T.; Murzyn, K.; Hinsen, K.; Kneller, G. R. nMol- dyn: A program package for a neutron scattering oriented analysis of molecular dynamics simulations.J. Computer. Chem.2003,24, 657–667, doi:10.1002/jcc.10243

  22. [22]

    Hinsen, K.; Pellegrini, E.; Stachura, S.; Kneller, G. R. nMoldyn 3: Using task farming for a parallel spectroscopy-oriented analysis of molecular dynamics simulations.J. Comput. Chem.2012,33, 2043–2048, doi:10.1002/jcc.23035

  23. [23]

    MDANSE: An In- teractive Analysis Environment for Molecular Dynam- ics Simulations.J

    Goret, G.; Aoun, B.; Pellegrini, E. MDANSE: An In- teractive Analysis Environment for Molecular Dynam- ics Simulations.J. Chem. Inf. Model.2017,57, 1–5, doi:10.1021/acs.jcim.6b00571

  24. [24]

    Two pressure cells for quasielastic and inelastic neu- tron scatterings

    Petersen, M. H.; Telling, M. T. F.; Kneller, G.; Bor- dallo, H. N. Revisiting the modeling of quasielastic neu- tron scattering from bulk water.EPJ Web Conf.2022, 272, 01012, doi:10.1051/epjconf/202227201012

  25. [25]

    J.; Adams, M

    Carlile, C. J.; Adams, M. A. The design of the IRIS inelastic neutron spectrometer and improvements to its analysers.Physica B: Condens. Matter1992,182, 431–440, doi:10.1016/0921-4526(92)90047-v

  26. [26]

    Heads or tails: how do chemically substituted fullerenes melt?Phys

    Armstrong, J.; Mukhopadhyay, S.; Bresme, F.; Fernandez-Alonso, F. Heads or tails: how do chemically substituted fullerenes melt?Phys. Chem. Chem. Phys. 2016,18, 17202–17209, doi:10.1039/c6cp01333c

  27. [27]

    J.; Devonport, M.; Farago, B.; Al- varez, F.; Mart´ ınez Gonz´ alez, J

    Arbe, A.; Nilsen, G. J.; Devonport, M.; Farago, B.; Al- varez, F.; Mart´ ınez Gonz´ alez, J. A.; Colmenero, J. Col- lective dynamics and self-motions in the van der Waals liquid tetrahydrofuran from meso- to inter-molecular scales disentangled by neutron spectroscopy with po- larization analysis.J. Chem. Phys.2023,158, 184502, doi:10.1063/5.0147427

  28. [28]

    Atomic Dy- namics in Simple Liquid: de Gennes Narrow- ing Revisited.Phys

    Wu, B.; Iwashita, T.; Egami, T. Atomic Dy- namics in Simple Liquid: de Gennes Narrow- ing Revisited.Phys. Rev. Lett.2018,120, 135502, doi:10.1103/physrevlett.120.135502

  29. [29]

    J.; Koˇ sata, J.; Devonport, M.; Galsworthy, P.; Bewley, R

    Nilsen, G. J.; Koˇ sata, J.; Devonport, M.; Galsworthy, P.; Bewley, R. I.; Voneshen, D. J.; Dalgliesh, R.; Stew- art, J. R. Polarisation analysis on the LET time-of-flight spectrometer.J. Phys. Conf. Ser.2017,862, 012019, doi:10.1088/1742-6596/862/1/012019

  30. [30]

    Here, the uncertainty is a single standard deviation of the marginal posterior distribution

  31. [31]

    J.; Chelli, R.; Righini, R

    Magro, A.; Frezzato, D.; Polimeno, A.; Moro, G. J.; Chelli, R.; Righini, R. Dynamics of liquid benzene: A cage analysis.J. Chem. Phys.2005,123, 124511, doi:10.1063/1.2018629

  32. [32]

    The uncertainty here indicates a single standard devia- tion from 17 samples

  33. [33]

    Schwartz, M.; Duan, D.; Berry, R. J. Molecular Dynamics Study of Anisotropic Translational and Rotational Diffu- sion in Liquid Benzene.J. Phys. Chem. A2005,109, 8637–8641, doi:10.1021/jp0530799

  34. [34]

    B´ ee, M.Quasielastic neutron scattering, principles and applications in solid state chemistry, biology and materi- als science; CRC Press: London, UK, 1988

  35. [35]

    S.; Skilling, J.Data Analysis: A Bayesian Tu- torial, 2nd ed.; Oxford University Press: Oxford, UK, 2006

    Sivia, D. S.; Skilling, J.Data Analysis: A Bayesian Tu- torial, 2nd ed.; Oxford University Press: Oxford, UK, 2006

  36. [36]

    R.; Cooper, J

    McCluskey, A. R.; Cooper, J. F. K.; Arnold, T.; Snow, T. A general approach to maximise information density in neutron reflectometry analysis.Mach. Learn.: Sci. Tech- nolo.2020,1, 035002, doi:10.1088/2632-2153/ab94c4

  37. [37]

    E.; Raftery, A

    Kass, R. E.; Raftery, A. E. Bayes Fac- tors.J. Am. Stat. Assoc.1995,90, 773–795, doi:10.1080/01621459.1995.10476572

  38. [38]

    2004, in American Institute of Physics Conference Series, Vol

    Skilling, J. Nested Sampling. AIP Conference Proceed- ings. 2004; p 395–405, doi:10.1063/1.1835238

  39. [39]

    CH/πInteraction between Benzene and Hydrocarbons Having Six Carbon Atoms in Their Binary Liquid Mixtures.J

    Kasahara, Y.; Suzuki, Y.; Kabasawa, A.; Minami, H.; Matsuzawa, H.; Iwahashi, M. CH/πInteraction between Benzene and Hydrocarbons Having Six Carbon Atoms in Their Binary Liquid Mixtures.J. Oleo Sci.2010,59, 21–29, doi:10.5650/jos.59.21

  40. [40]

    T.; Manyar, H

    Falkowska, M.; Bowron, D. T.; Manyar, H. G.; Hardacre, C.; Youngs, T. G. A. Neutron Scattering of Aromatic and Aliphatic Liquids.ChemPhysChem2016, 17, 2043–2055, doi:10.1002/cphc.201600149

  41. [41]

    Headen, T. F. Temperature dependent structural changes in liquid benzene studied using neutron diffraction.Mol. Phys.2019,117, 3329–3336, doi:10.1080/00268976.2019.1631496

  42. [42]

    R.; Armstrong, J.; McColl, K.; Richard- son, H

    McCluskey, A. R.; Armstrong, J.; McColl, K.; Richard- son, H. Assessing the quality of machine learned force fields for dynamical properties using quasi-elastic neu- tron scattering. 2025; doi:10.5286/ISIS.E.RB2420517-1

  43. [43]

    R.; Armstrong, J.; Richardson, H

    McCluskey, A. R.; Armstrong, J.; Richardson, H. p-QENS measurement of liquid Benzene. 2025; doi:10.5286/ISIS.E.RB2590225-1

  44. [44]

    Dataset for ”Lost in Translation: Simulation-Informed Bayesian Infer- ence for Quasi-Elastic Neutron Scattering”

    Richardson, H.; McCluskey, A. Dataset for ”Lost in Translation: Simulation-Informed Bayesian Infer- ence for Quasi-Elastic Neutron Scattering”. 2026; doi:10.5281/zenodo.18836632

  45. [45]

    fitting-qens-0.0.1

    Richardson, H.; McCluskey, A. fitting-qens-0.0.1. https://github.com/scams-research/fitting-qens, 2026

  46. [46]

    P.; Aktulga, H

    Thompson, A. P.; Aktulga, H. M.; Berger, R.; Bolin- tineanu, D. S.; Brown, W. M.; Crozier, P. S.; in ’t Veld, P. J.; Kohlmeyer, A.; Moore, S. G.; Nguyen, T. D.; Shan, R.; Stevens, M. J.; Tranchida, J.; Trott, C.; Plimp- ton, S. J. LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales.Comput. Ph...

  47. [47]

    S.; Cabeza de Vaca, I.; Tirado-Rives, J.; Jorgensen, W

    Dodda, L. S.; Cabeza de Vaca, I.; Tirado-Rives, J.; Jorgensen, W. L. LigParGen web server: an auto- matic OPLS-AA parameter generator for organic lig- ands.Nucleic Acids Res.2017,45, W331–W336, doi:10.1093/nar/gkx312

  48. [48]

    L.; Maxwell, D

    Jorgensen, W. L.; Maxwell, D. S.; Tirado-Rives, J. Devel- opment and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids.J. Am. Chem. Soc.1996,118, 11225–11236, doi:10.1021/ja9621760

  49. [49]

    A unified formulation of the constant tempera- ture molecular dynamics methods.J

    Nos´ e, S. A unified formulation of the constant tempera- ture molecular dynamics methods.J. Chem. Phys.1984, 81, 511–519, doi:10.1063/1.447334

  50. [50]

    Hoover, W. G. Canonical dynamics: Equilibrium phase- space distributions.Phys. Rev. A1985,31, 1695–1697, doi:10.1103/physreva.31.1695

  51. [51]

    J.; Klein, M

    Martyna, G. J.; Klein, M. L.; Tuckerman, M. Nos´ e–Hoover chains: The canonical ensemble via con- tinuous dynamics.J. Chem. Phys.1992,97, 2635–2643, doi:10.1063/1.463940. 9

  52. [52]

    J.; Tobias, D

    Martyna, G. J.; Tobias, D. J.; Klein, M. L. Constant pressure molecular dynamics algorithms.J. Chem. Phys. 1994,101, 4177–4189, doi:10.1063/1.467468

  53. [53]

    Dubkov, and Davide Valenti

    Cassella, G.; Stewart, J. R.; Patern` o, G. M.; Garc´ ıa Sakai, V.; Devonport, M.; Galsworthy, P. J.; Be- wley, R. I.; Voneshen, D. J.; Raspino, D.; Nilsen, G. J. Polarization analysis on the LET cold neutron spec- trometer using a 3He spin-filter: first results.J. Phys.: Conf. Ser.2019,1316, 012007, doi:10.1088/1742- 6596/1316/1/012007

  54. [54]

    R.; Squires, A

    McCluskey, A. R.; Squires, A. G.; Dunn, J.; Coles, S. W.; Morgan, B. J. kinisi: Bayesian analysis of mass transport from molecular dynamics simulations.J. Open Source Softw.2024,9, 5984, doi:10.21105/joss.05984

  55. [55]

    Richardson, H.; Dunn, J.; McCluskey, A. R. On the es- timation of center of mass in periodic systems.J. Chem. Phys.2025,162, 204103, doi:10.1063/5.0260928

  56. [56]

    R.; Coles, S

    McCluskey, A. R.; Coles, S. W.; Morgan, B. J. Accu- rate Estimation of Diffusion Coefficients and their Un- certainties from Computer Simulation.J. Chem. Theory Comput.2024,21, 79–87, doi:10.1021/acs.jctc.4c01249

  57. [57]

    Favro, L. D. Theory of the Rotational Brownian Mo- tion of a Free Rigid Body.Phys. Rev.1960,119, 53–62, doi:10.1103/physrev.119.53

  58. [58]

    These are analogous to thelandm l orbital angular and magnetic quantum numbers

  59. [59]

    P.Group Theory and Its Applications to the Quantum Mechanics of Atomic Spectra, 4th ed.; Aca- demic Press: New York, US, 1964

    Wigner, E. P.Group Theory and Its Applications to the Quantum Mechanics of Atomic Spectra, 4th ed.; Aca- demic Press: New York, US, 1964

  60. [60]

    Woessner, D. E. Molecular Reorientation in Liquids. Deuteron Quadrupole Relaxation in Liquid Deuterium Oxide and Perdeuterobenzene.J. Chem. Phys.1964,40, 2341–2348, doi:10.1063/1.1725512

  61. [61]

    Speagle, J. S. dynesty: a dynamic nested sampling package for estimating Bayesian posteriors and evi- dences.Mon. Not. R. Astron. Soc2020,493, 3132–3158, doi:10.1093/mnras/staa278

  62. [62]

    Laaksonen, A.; Stilbs, P.; Wasylishen, R. E. Molecu- lar motion and solvation of benzene in water, carbon tetrachloride, carbon disulfide and benzene: A combined molecular dynamics simulation and nuclear magnetic resonance study.J. Chem. Phys.1998,108, 455–468, doi:10.1063/1.475408

  63. [63]

    Molecular Dynamics of Benzene in Neat Liquid and a Solution Containing Polystyrene

    Witt, R.; Sturz, L.; D¨ olle, A.; M¨ uller-Plathe, F. Molecular Dynamics of Benzene in Neat Liquid and a Solution Containing Polystyrene. 13C Nuclear Mag- netic Relaxation and Molecular Dynamics Simulation Results.J. Phys. Chem. A2000,104, 5716–5725, doi:10.1021/jp000201p

  64. [64]

    Fu, C.-F.; Tian, S. X. A Comparative Study for Molecular Dynamics Simulations of Liquid Benzene. J. Chem. Theory Comput.2011,7, 2240–2252, doi:10.1021/ct2002122

  65. [65]

    Lost in Translation: Simulation-Informed Bayesian Inference Improves Understanding of Molecular Motion From Neutron Scattering

    Johnson, M. J.; Kearley, G. J. InNeutron Scattering – Fundamentals; Fernandez-Alonso, F., Price, D. L., Eds.; Experimental Methods in the Physical Sciences; Elsevier, 2013; Vol. 44; p 415–469, doi:10.1016/b978-0-12-398374- 9.00007-3. Supplemental Material for “Lost in Translation: Simulation-Informed Bayesian Inference Improves Understanding of Molecular ...

  66. [66]

    The sample was loaded onto the IRIS instrument and measured for 190 min at 290 K

    Sample and Measurement On IRIS, neat benzene (Sigma-Aldrich, 99.8 % purity) was loaded into aluminium sample cans, with a sample thickness of 0.25 mm, in a fume hood to minimise exposure to moisture in the air, and the cans were sealed with indium wire. The sample was loaded onto the IRIS instrument and measured for 190 min at 290 K. An empty aluminium ca...

  67. [67]

    white-beam

    Data Reduction The experimental IRIS data was binned into 12Q-bins, each with a width of 0.12 ˚A −1 across the 51 detectors with a totalQ-range of (0.44 to 1.85) ˚A −1 . A vanadium sample was used to determine the instrument’s resolution function. The same sample was used to correct for detector efficiency, by normalising each QENS spectrum with the energ...

  68. [68]

    Specifically, this considers the oblate symmetric top benzene

    Rotational Autocorrelation Function We start by deriving the rank-lcorrelation function for an anisotropically axially symmetric rotator. Specifically, this considers the oblate symmetric top benzene. The WignerD-functions, which form an orthonormal basis set for the 3D rotation group, SO(3), have the following eigenvalues [57], λlm =D tl(l+ 1) +m 2(Ds −D...

  69. [69]

    We consider a rigid molecule with a single scattering nucleus located at a fixed distanceRfrom the molecular centre

    Dynamic Structure Factor Now, we turn our focus to the derivation of the dynamic structure factor. We consider a rigid molecule with a single scattering nucleus located at a fixed distanceRfrom the molecular centre. Its position in the laboratory frame is, r(t) =R ˆu(t).(SI.7) where ˆuis some unit factor that the change of which will be probed. The interm...