Fingerprints of hydrogen bonding in the terahertz dynamics of ethanol and water: An inelastic x-ray scattering study.

We used inelastic x-ray scattering methods to measure the terahertz spectrum of density fluctuations of ethanol in both liquid and solid phases. The results of a Bayesian inference-based lineshape analysis with a multiple excitation model and the comparison with a previous similar analysis on water indicate that the different structures induced by hydrogen bonds in ethanol and water have a profound influence on the respective dynamic responses, the latter being characterized by longer living and better resolved high-frequency acoustic excitations. In addition, we compare these findings with those obtained with an alternative approach based on the exponential expansion theory and ensuring sum rules fulfillment, demonstrating that the model's choice directly impacts the number of spectral modes detected.

Ice phonon spectra and Bayes inference: A gateway to a new understanding of terahertz sound propagation in water.

Understanding how molecules engage in collective motions in a liquid where a network of bonds exists has both fundamental and applied relevance. On the one hand, it can elucidate the "ordering" role of long-range correlations and inspire new avenues to control such order to implement sound manipulation. Water represents an ideal investigation case to unfold these general aspects, and, across the decades, it has been the focus of thorough scrutiny. Despite this investigative effort, the spectrum of terahertz density fluctuations of water largely remains a puzzle for condensed matter physicists. To unravel it, we compare previous scattering measurements of water spectra with new ones on ice. Owing to the unique asset of Bayesian inference, we draw a more detailed portrayal of the phonon response of ice. The comparison with the one of liquid water challenges the current understanding of density fluctuations in water, or more in general, of any networked liquid.

Density and time scaling effects on the velocity autocorrelation function of quantum and classical dense fluid para-hydrogen.

We report the results of a ring polymer molecular dynamics study of the Kubo velocity autocorrelation function of a quantum fluid as para-hydrogen aimed at the comparison with its classical counterpart. Quite different density conditions were considered for both the classical and quantum cases, in order to compare the two systems before and after the dynamical crossover typically undergone by the velocity autocorrelation function (VAF) of fluids at densities around the triple point, where the shape of the function changes from a monotonic to an oscillatory behavior with a negative minimum. A detailed study of the phase diagram of classical para-hydrogen was necessary for a reasonable choice of the classical states to be taken into consideration, in the spirit of the classical principle of corresponding states. The shape of the quantum and classical VAF was thoroughly analyzed, exhibiting at all studied densities clear differences that might be taken as evidence of quantum effects. We show that these differences are substantially reduced by applying a state-dependent time scaling with respect to a reference time identified with the inverse of the collision rate. An even better coincidence in shape is found by comparing the two systems at slightly non-corresponding reduced densities, suggesting that the quantum system behaves almost like the classical one, but at systematically less dense reduced states of the latter. We also find an unexpected and quite interesting density trend of the collision rate of both classical and quantum para-hydrogen, which accounts for the effectiveness of the scaling throughout the explored density range. The mean kinetic energy and the diffusion coefficients are also discussed in some detail.

Density of states from mode expansion of the self-dynamic structure factor of a liquid metal.

We show that by exploiting multi-Lorentzian fits of the self-dynamic structure factor at various wave vectors it is possible to carefully perform the Q→0 extrapolation required to determine the spectrum Z(ω) of the velocity autocorrelation function of a liquid. The smooth Q dependence of the fit parameters makes their extrapolation to Q=0 a simple procedure from which Z(ω) becomes computable, with the great advantage of solving the problems related to resolution broadening of either experimental or simulated self-spectra. Determination of a single-particle property like the spectrum of the velocity autocorrelation function turns out to be crucial to understanding the whole dynamics of the liquid. In fact, we demonstrate a clear link between the collective mode frequencies and the shape of the frequency distribution Z(ω). In the specific case considered in this work, i.e., liquid Au, analysis of Z(ω) revealed the presence, along with propagating sound waves, of lower frequency modes that were not observed before by means of dynamic structure factor measurements. By exploiting ab initio simulations for this liquid metal we could also calculate the transverse current-current correlation spectra and clearly identify the transverse nature of the above mentioned less energetic modes. Evidence of propagating transverse excitations has actually been reported in various works in the recent literature. However, in some cases, like the present one, these modes are difficult to detect in density fluctuation spectra. We show here that the analysis of the single-particle dynamics is able to unveil their presence in a very effective way. The properties here shown to characterize Z(ω), and the information in it contained therefore allow us to identify it with the density of states (DoS) of the liquid. We demonstrate that only nonhydrodynamic modes contribute to the DoS, thus establishing its purely microscopic origin. Finally, as a by-product of this work, we provide our estimate of the self-diffusion coefficient of liquid gold just above melting.

Density of states and dynamical crossover in a dense fluid revealed by exponential mode analysis of the velocity autocorrelation function.

Extending a preceding study of the velocity autocorrelation function (VAF) in a simulated Lennard-Jones fluid [Phys. Rev. E 92, 042166 (2015)PLEEE81539-375510.1103/PhysRevE.92.042166] to cover higher-density and lower-temperature states, we show that the recently demonstrated multiexponential expansion method allows for a full account and understanding of the basic dynamical processes encompassed by a fundamental quantity as the VAF. In particular, besides obtaining evidence of a persisting long-time tail, we assign specific and unambiguous physical meanings to groups of exponential modes related to the longitudinal and transverse collective dynamics, respectively. We have made this possible by consistently introducing the interpretation of the VAF frequency spectrum as a global density of states in fluids, generalizing a solid-state concept, and by giving to specific spectral components, obtained through the VAF exponential expansion, the corresponding meaning of partial densities of states relative to specific dynamical processes. The clear identification of a high-frequency oscillation of the VAF with the near-top excitation frequency in the dispersion curve of acoustic waves is a neat example of the power of the method. As for the transverse mode contribution, its analysis turns out to be particularly important, because the multiexponential expansion reveals a transition marking the onset of propagating excitations when the density is increased beyond a threshold value. While this finding agrees with the recent literature debating the issue of dynamical crossover boundaries, such as the one identified with the Frenkel line, we can add detailed information on the modes involved in this specific process in the domains of both time and frequency. This will help obtain a still missing full account of transverse dynamics, in both its nonpropagating and propagating aspects which are linked through dynamical transitions depending on both the thermodynamic states and the excitation wave vectors.

Bayesian approach to the analysis of neutron Brillouin scattering data on liquid metals.

When the dynamics of liquids and disordered systems at mesoscopic level is investigated by means of inelastic scattering (e.g., neutron or x ray), spectra are often characterized by a poor definition of the excitation lines and spectroscopic features in general and one important issue is to establish how many of these lines need to be included in the modeling function and to estimate their parameters. Furthermore, when strongly damped excitations are present, commonly used and widespread fitting algorithms are particularly affected by the choice of initial values of the parameters. An inadequate choice may lead to an inefficient exploration of the parameter space, resulting in the algorithm getting stuck in a local minimum. In this paper, we present a Bayesian approach to the analysis of neutron Brillouin scattering data in which the number of excitation lines is treated as unknown and estimated along with the other model parameters. We propose a joint estimation procedure based on a reversible-jump Markov chain Monte Carlo algorithm, which efficiently explores the parameter space, producing a probabilistic measure to quantify the uncertainty on the number of excitation lines as well as reliable parameter estimates. The method proposed could turn out of great importance in extracting physical information from experimental data, especially when the detection of spectral features is complicated not only because of the properties of the sample, but also because of the limited instrumental resolution and count statistics. The approach is tested on generated data set and then applied to real experimental spectra of neutron Brillouin scattering from a liquid metal, previously analyzed in a more traditional way.

Time dependence of the velocity autocorrelation function of a fluid: An eigenmode analysis of dynamical processes.

The velocity autocorrelation function (VAF), a key quantity in the atomic-scale dynamics of fluids, has been the first paradigmatic example of a long-time tail phenomenon, and much work has been devoted to detecting such long-lasting correlations and understanding their nature. There is, however, much more to the VAF than simply the evidence of this long-time dynamics. A unified description of the VAF from very short to long times, and of the way it changes with varying density, is still missing. Here we show that an approach based on very general principles makes such a study possible and opens the way to a detailed quantitative characterization of the dynamical processes involved at all time scales. From the analysis of molecular dynamics simulations for a slightly supercritical Lennard-Jones fluid at various densities, we are able to evidence the presence of distinct fast and slow decay channels for the velocity correlation on the time scale set by the collision rate. The density evolution of these decay processes is also highlighted. The method presented here is very general, and its application to the VAF can be considered as an important example.

Comment on 'An effective fitting scheme for the dynamic structure of pure liquids'.

The occurrence of a propagation gap in the acoustic excitations of a liquid is excluded by Wax and Bryk (2013 J. Phys.: Condens. Matter 25 325104). The requirement of a finite second frequency moment for the dynamic structure factor is used to come to this conclusion. We show here that this requirement does not conflict with the existence of overdamped, non-propagating modes which give rise to spectra that do not contain inelastic components. Such a behaviour has indeed been detected in the analysis of the collective dynamics of several liquids, carried out by using well-established sum-rule-compliant S(q, ω) models.

Orientational and translational correlations of liquid methane over the nanometer-picosecond scales by molecular dynamics simulation and inelastic neutron scattering.

Five models for the site-site intermolecular pair interactions of methane are compared in some detail and used to investigate both structural and dynamical properties of the dense liquid deuteromethane by means of molecular dynamics (MD) simulations. The orientational distribution probabilities of molecular pairs are carefully analyzed for each anisotropic potential model. We propose a revision of existing classification methods used to group the innumerable relative orientations of methane-methane pairs into six basic geometries. With this new approach, our results for the probability of the six basic categories as a function of the intermolecular distance are different from the ones present in the literature, where the role of the angular spread on the anisotropic interaction energy is not taken in full consideration and certain configurations with no significant change in the pair-potential are assigned to different categories. The analysis of the static orientational correlations in liquid methane and the prevalence of certain configurations in different ranges guide the subsequent discussion of the MD model-dependent results for the dynamic structure factor. Comparison with our inelastic neutron scattering results for liquid CD(4) at the nanometer and picosecond space and time scales allows us to confirm the full adequacy of the Tsuzuki, Uchimaru and Tanabe model of 1998 with respect to more recent potentials.

Structure of liquid n-hexane.

The structure of deuterated liquid n-hexane has been investigated at room temperature by neutron diffraction and molecular dynamics simulations. By carrying out a careful analysis of the measurements, experimental data were obtained in very good agreement with the simulated data. This allowed a thorough analysis of the simulation results aiming at the evaluation of the partial, intra-, and intermolecular components of the n-hexane structure. We finally compare the intramolecular differential cross sections calculated from the most probable n-hexane molecular configurations with the measured and simulated data.

Inelastic neutron scattering and molecular dynamics determination of the interaction potential in liquid CD4.

Anisotropic interactions of liquid CD4 are studied in detail by comparison of inelastic neutron Brillouin scattering data with molecular dynamics simulations using up to four different models of the methane site-site potential. We demonstrate that the experimental dynamic structure factor S(Q,omega) acts as a highly discriminating quantity for possible interaction schemes. In particular, the Q evolution of the spectra enables a selective probing of the short- and medium-range features of the anisotropic potentials. We show that the preferential configuration of methane dimers at liquid densities can thus be discerned by analyzing the orientation-dependent model potential curves, in light of the experimental and simulation results.

Neutron Brillouin scattering study of collective dynamics in a dense He-Ne gaseous mixture.

The dynamic structure factor S(k,omega) of a 77% He and 23% Ne gaseous mixture at T = 39.3 K and total number density n = 15.8 nm(-3) has been measured by inelastic neutron scattering at small angles. In the range of wave vectors studied, 0.7

Structure of gaseous Kr in the low-q region by neutron scattering and interaction potentials.

Accurate experimental information on the long-range triplet interactions in noble fluids, as well as on the two-body potential, can be obtained from the low-density behavior of the static structure factor S(q) in the small-q region. The results here reported of a recent low-q neutron diffraction investigation in Kr, devoted to undercritical densities in the range 2.4