Temperature-dependent dynamic structure factors for liquid water inferred from inelastic neutron scattering measurements

Abstract

Temperature-dependent dynamic structure factors S(Q, ω) for liquid water have been calculated using a composite model, which is based on the decoupling approximation of the mean square displacement of the water molecules into diffusion and solid-like vibrational parts. The solid-like vibrational part Svib(Q, ω) is calculated with the phonon expansion method established in the framework of the incoherent Gaussian approximation. The diffusion part Sdiff(Q, ω) relies on the Egelstaff-Schofield translational diffusion model corrected for jump diffusions and rotational diffusions with the Singwi-Sjölander random model and Sears expansion, respectively. Systematics of the model parameters as a function of temperature were deduced from quasi-elastic neutron scattering data analysis reported in the literature and from molecular dynamics (MD) simulations relying on the TIP4P/2005f model. The resulting S(Q, ω) values are confronted by means of Monte Carlo simulations to inelastic neutron scattering data measured with IN4, IN5, and IN6 time-of-flight spectrometers of the Institut Laue-Langevin (ILL) (Grenoble, France). A modest range of temperatures (283-494 K) has been investigated with neutron wavelengths corresponding to incident neutron energies ranging from 0.57 to 67.6 meV. The neutron-weighted multiphonon spectra deduced from the ILL data indicate a slight overestimation by the MD simulations of the frequency shift and broadening of the librational band. The descriptive power of the composite model was suited for improving the comparison to experiments via Bayesian updating of prior model parameters inferred from MD simulations. The reported posterior temperature-dependent densities of state of hydrogen in H2O would represent valuable insights for studying the collective coupling interactions in the water molecule between the inter- and intramolecular degrees of freedom.