Microscopic self-dynamics in liquid Ne-D_{2} mixtures: Quantum features and itinerant oscillators reexamined.

In this paper, we report the results of a centroid molecular dynamics (CMD) study of the canonical velocity autocorrelation functions (VACFs) in liquid Ne-D_{2} mixtures at a temperature of T=30K and in the full D_{2}-concentration range (0%≤x_{D_{2}}≤100%). This binary system was selected because of its moderate, although sizable, quantum effects which, as far as its equilibrium properties are concerned, are fully described by the path integral Monte Carlo (PIMC) simulations that have been also implemented. A comprehensive test of the VACF spectral moments carried out using three physical quantities (namely, mean kinetic energy, Einstein frequency, and mean-squared force) obtained from PIMC was performed revealing the potentialities, as well as the limitations, of the CMD approach to the single-particle dynamics in these low-T liquid mixtures. Additional physical information was extracted from the canonical VACFs by fitting their spectra via two distinct methods: the Levesque-Verlet model (LV, very flexible but highly heuristic) and the itinerant oscillator model (IO, based on the physical ground of a single particle rattling inside a short-lived diffusing pseudocage). Both provided good fits of the CMD outputs, with LV being always more adequate than IO in the case of the Ne VACFs, while, as for the D_{2} VACFs, the LV superiority is evident only at high x_{D_{2}} values. However, a peculiar and systematic effect was found after analyzing the IO-fitted parameters: the estimated pseudocage masses turned out to be at least one order of magnitude lower than the corresponding values inferred from the PIMC simulations. This outcome concerns both the Ne and the D_{2} rattling molecule and, as we also discovered, had already been observed (but promptly forgotten) in purely classical simulations of liquid Ar. The possible physical origins of this finding have been finally discussed in some detail, also in connection with the result of the more recent exponential expansion theory (EET), which manages to shed more light on the concept of single particles rattling inside short-lived pseudocages, ultimately demonstrating its untenability.

Large-cage occupation and quantum dynamics of hydrogen molecules in sII clathrate hydrates.

Hydrogen clathrate hydrates are ice-like crystalline substances in which hydrogen molecules are trapped inside polyhedral cages formed by the water molecules. Small cages can host only a single H2 molecule, while each large cage can be occupied by up to four H2 molecules. Here, we present a neutron scattering study on the structure of the sII hydrogen clathrate hydrate and on the low-temperature dynamics of the hydrogen molecules trapped in its large cages, as a function of the gas content in the samples. We observe spectral features at low energy transfer (between 1 and 3 meV), and we show that they can be successfully assigned to the rattling motion of a single hydrogen molecule occupying a large water cage. These inelastic bands remarkably lose their intensity with increasing the hydrogen filling, consistently with the fact that the probability of single occupation (as opposed to multiple occupation) increases as the hydrogen content in the sample gets lower. The spectral intensity of the H2 rattling bands is studied as a function of the momentum transfer for partially emptied samples and compared with three distinct quantum models for a single H2 molecule in a large cage: (i) the exact solution of the Schrödinger equation for a well-assessed semiempirical force field, (ii) a particle trapped in a rigid sphere, and (iii) an isotropic three-dimensional harmonic oscillator. The first model provides good agreement between calculations and experimental data, while the last two only reproduce their qualitative trend. Finally, the radial wavefunctions of the three aforementioned models, as well as their potential surfaces, are presented and discussed.

Role of the single-particle dynamics in the transverse current autocorrelation function of a liquid metal.

A recent simulation study of the transverse current autocorrelation of the Lennard-Jones fluid [Guarini et al., Phys. Rev. E 107, 014139 (2023)] revealed that this function can be perfectly described within the exponential expansion theory [Barocchi et al., Phys. Rev. E 85, 022102 (2012)]. However, above a certain wavevector Q, not only transverse collective excitations were found to propagate in the fluid, but a second oscillatory component of unclear origin (therefore called X) must be considered to fully account for the time dependence of the correlation function. Here, we present an extended investigation of the transverse current autocorrelation of liquid Au as obtained by ab initio molecular dynamics in the very wide range of wavevectors 5.7 ≤ Q ≤ 32.8 nm-1 in order to also follow the behavior of the X component, if present, at large Q values. A joint analysis of the transverse current spectrum and its self-portion indicates that the second oscillatory component arises from the longitudinal dynamics, as suggested by its close resemblance with the previously determined component accounting for the longitudinal part of the density of states. We conclude that such a mode, albeit featuring a merely transverse property, fingerprints the effect of longitudinal collective excitations on single-particle dynamics, rather than arising from a possible coupling between transverse and longitudinal acoustic waves.

Onset of collective excitations in the transverse dynamics of simple fluids.

A thorough analysis of the transverse current autocorrelation function obtained by molecular dynamics simulations of a dense Lennard-Jones fluid reveals that even such a simple system is characterized by a varied dynamical behavior with changing length scale. By using the exponential expansion theory, we provide a full account of the time correlation at wavevectors Q between the upper boundary of the hydrodynamic region and Q_{p}/2, with Q_{p} being the position of the main peak of the static structure factor. In the Q range studied, we identify and accurately locate the wavevector at which shear wave propagation starts to take place, and show clearly how this phenomenon may be represented by a damped harmonic oscillator changing, in a continuous way, from an overdamped to an underdamped condition. The decomposition into exponential modes allows one to convincingly establish not only the crossover related to the onset of transverse waves but, surprisingly, also the existence of a second pair of modes equivalent to another oscillator that undergoes, at higher Q values, a similarly smooth over to underdamped transition.

Microscopic collective dynamics in liquid neon-deuterium mixtures: Inelastic neutron scattering and quantum simulations.

In this paper a combined neutron scattering and quantum simulation study of the collective dynamics in liquid Ne-D_{2} mixtures, at a temperature of T=30K and in the wave-vector transfer range 4 nm^{-1}

Proton Dynamics in Palladium-Silver: An Inelastic Neutron Scattering Investigation.

Proton dynamics in Pd77Ag23 membranes is investigated by means of various neutron spectroscopic techniques, namely Quasi Elastic Neutron Scattering, Incoherent Inelastic Neutron Scattering, Neutron Transmission, and Deep Inelastic Neutron Scattering. Measurements carried out at the ISIS spallation neutron source using OSIRIS, MARI and VESUVIO spectrometers were performed at pressures of 1, 2, and 4 bar, and temperatures in the 330-673 K range. The energy interval spanned by the different instruments provides information on the proton dynamics in a time scale ranging from about 102 to 10-4 ps. The main finding is that the macroscopic diffusion process is determined by microscopic jump diffusion. In addition, the vibrational density of states of the H atoms in the metal lattice has been determined for a number of H concentrations and temperatures. These measurements follow a series of neutron diffraction experiments performed on the same sample and thus provide a complementary information for a thorough description of structural and dynamical properties of H-loaded Pd-Ag membranes.

Time dependence of quantum correlation functions.

In the past few years, the exponential expansion analysis of time autocorrelation functions has provided profound insight into the leading microscopic processes driving the atomic-scale dynamics and has made it possible to highlight the presence and the role of various relaxation channels through which the fundamental correlation functions decay with time. Here we apply this method to the determination of the full time dependence of a correlation function c(t) in a quantum system at nonzero temperature, by making explicit its relationship with its Kubo transform c_{K}(t), which in some cases can be approximately computed with the presently available quantum simulation techniques. We obtain an exact expression for c(t) in terms of the exponential modes that describe the time behavior of c_{K}(t). The relative importance of the various modes in determining the overall shape of c(t) can then be studied in detail. This work extends to the full time domain the results of a previous paper [Guarini et al., Phys. Rev. Lett. 123, 135301 (2019)PRLTAO0031-900710.1103/PhysRevLett.123.135301], in which we employed the same method to calculate the zero time value of the velocity autocorrelation function, to obtain a microscopic description of the quantum mean kinetic energy in a fluid. In particular, we show that the decay constants and the frequencies of the dominant microscopic modes of c(t) are the same as those of c_{K}(t), but the dynamics of the quantum system also contains an additional term decaying on a time scale determined solely by temperature of the system.

Dynamical Origin of the Total and Zero-Point Kinetic Energy in a Quantum Fluid.

By applying an exponential mode analysis to ring polymer molecular dynamics simulations of dense fluid parahydrogen, we find that the dynamical processes establishing the time behavior of the Kubo velocity autocorrelation function display the same nature as those already observed in high-density classical fluids. This result permits us to demonstrate that the exponential mode decomposition is a unique tool to identify which dynamical processes lead to one of the most notable properties of quantum fluids: the large value of the mean kinetic energy per particle and the importance of the zero-temperature quantum effects in determining it.

Hydrogen self-dynamics in diluted liquid mixtures with neon: An inelastic neutron scattering study.

We have measured the dynamic structure factor of liquid neon-hydrogen mixtures (T=30.1 K) at two different H_{2} concentration levels (namely, 3.4% and 10%) making use of inelastic neutron scattering. This system has been selected since the presence of heavy Ne atoms strongly influences the self-dynamics of the H_{2} centers of mass via the formation of short-lived cages, which act both on the vibrational and the diffusive parts of the single-particle motion. After operating a standard data reduction and the subtraction of the Ne signal, experimental neutron spectra were analyzed through a generalization of the Young and Koppel model, and the H_{2} center-of-mass self-dynamic structure factor was finally extracted for the two liquid samples. Important physical quantities (namely, single-particle mean kinetic energy and self-diffusion coefficient) were estimated from the experimental data and then compared with quantum dynamical calculations, which also provided simulations of the velocity autocorrelation functions for Ne atoms and H_{2} centers of mass. The latter estimates, in the framework of the well-known Gaussian approximation, were used to simulate the H_{2} center-of-mass self-dynamic structure factor in the same kinematic range and thermodynamic conditions of the neutron scattering one. The comparison between measured and calculated spectra turned out to be qualitatively good, but some discrepancies, especially in the low-frequency part, seem to reinforce the idea of the existence of relevant non-Gaussian effects as in the case of pure hydrogen and H_{2}-D_{2} mixtures.

Density dependence of the dynamical processes governing the velocity autocorrelation function of a quantum fluid.

We present an exponential mode analysis of the dynamical processes determining the time behavior of the Kubo velocity autocorrelation function (KVAF) of fluid para-H_{2}, as obtained by ring polymer molecular dynamics simulations at various fluid densities. The mechanisms contributing to the decay of the KVAF are thoroughly characterized at a slightly supercritical temperature, in a density interval ranging from the critical point to the fluid-solid transition. We show that the quantum nature of the system does not influence the specific phenomena and decay channels through which a loss in velocity correlation takes place, since these are the same as found in classical fluids. Similarly, a dynamical crossover is observed with increasing density, signaling the onset of a transverse-like dynamics like in classical systems. We also investigate the effect of density on the processes contributing to the most relevant property of a quantum fluid, namely, the large values of the total and zero-point kinetic energy arising through the Heisenberg uncertainty principle.

VESPA: The vibrational spectrometer for the European Spallation Source.

VESPA, Vibrational Excitation Spectrometer with Pyrolytic-graphite Analysers, aims to probe molecular excitations via inelastic neutron scattering. It is a thermal high resolution inverted geometry time-of-flight instrument designed to maximise the use of the long pulse of the European Spallation Source. The wavelength frame multiplication technique was applied to provide simultaneously a broad dynamic range (about 0-500 meV) while a system of optical blind choppers allows to trade flux for energy resolution. Thanks to its high flux, VESPA will allow the investigation of dynamical and in situ experiments in physical chemistry. Here we describe the design parameters and the corresponding McStas simulations.

Impact of the Condensed-Phase Environment on the Translation-Rotation Eigenstates and Spectra of a Hydrogen Molecule in Clathrate Hydrates.

We systematically investigate the manifestations of the condensed-phase environment of the structure II clathrate hydrate in the translation-rotation (TR) dynamics and the inelastic neutron scattering (INS) spectra of an H2 molecule confined in the small dodecahedral cage of the hydrate. The aim is to elucidate the extent to which these properties are affected by the clathrate water molecules beyond the confining cage and the proton disorder of the water framework. For this purpose, quantum calculations of the TR eigenstates and INS spectra are performed for H2 inside spherical clathrate domains of gradually increasing radius and the number of water molecules ranging from 20 for the isolated small cage to more than 1800. For each domain size, several hundred distinct hydrogen-bonding topologies are constructed in order to simulate the effects of the proton disorder. Our study reveals that the clathrate-induced splittings of the j = 1 rotational level and the translational fundamental of the guest H2 are influenced by the condensed-phase environment to a dramatically different degree, the former very strongly and the latter only weakly.

Hydrogen self-dynamics in liquid H(2)-D(2) mixtures studied through inelastic neutron scattering.

We have measured the dynamic structure factor of liquid para-hydrogen mixed with normal deuterium (T=20 K) at two different concentration levels using incoherent inelastic neutron scattering. This choice has been made since the presence of D(2} modifies the self-dynamics of H(2) in a highly nontrivial way, acting both on its pseudophononic and its diffusive parts in a tunable way. After an accurate data reduction, recorded neutron spectra were studied through the modified Young and Koppel model and the H(2) center-of-mass self-dynamics structure factor was finally extracted for the two mixtures. Some physical quantities (i.e., self-diffusion coefficient and mean kinetic energy) were determined and compared with accurate quantum calculations, which, in addition, also provided estimates of the velocity autocorrelation function for the H(2) centers of mass. These estimates, in conjunction with the Gaussian approximation, were used to simulate the H(2) center-of-mass self-dynamics structure factor in the same range as the experimental one. The agreement between measured and calculated spectra was globally good, but some discrepancies proved the unquestionable breakdown of the Gaussian approximation in these semiquantum systems at a level comparable to that already observed in pure liquid para-hydrogen.

The HD molecule in small and medium cages of clathrate hydrates: quantum dynamics studied by neutron scattering measurements and computation.

We report inelastic neutron scattering (INS) measurements on molecular hydrogen deuteride (HD) trapped in binary cubic (sII) and hexagonal (sH) clathrate hydrates, performed at low temperature using two different neutron spectrometers in order to probe both energy and momentum transfer. The INS spectra of binary clathrate samples exhibit a rich structure containing sharp bands arising from both the rotational transitions and the rattling modes of the guest molecule. For the clathrates with sII structure, there is a very good agreement with the rigorous fully quantum simulations which account for the subtle effects of the anisotropy, angular and radial, of the host cage on the HD microscopic dynamics. The sH clathrate sample presents a much greater challenge, due to the uncertainties regarding the crystal structure, which is known only for similar crystals with different promoter, but nor for HD (or H2) plus methyl tert-butyl ether (MTBE-d12).

Experimental inelastic neutron scattering spectrum of hydrogen hexagonal clathrate-hydrate compared with rigorous quantum simulations.

We have performed high-resolution inelastic neutron scattering (INS) measurements on binary hydrogen clathrate hydrates exhibiting the hexagonal structure (sH). Two samples, differing only in the ortho/para fraction of hydrogen, were prepared using heavy water and methyl tert-butyl ether as the promoter in its perdeuterated form. The INS spectrum of the translation-rotation (TR) excitations of the guest H2 molecule was obtained by subtracting the very weak signal due to the D2O lattice modes. By means of a subtraction procedure, it has been possible to obtain separately the spectra of caged p-H2 and o-H2. sH clathrates are comprised of three distinct types of cages, two of which, differing in shape and size, are each occupied by one H2 molecule only. Both contribute to the measured INS spectrum which is, therefore, rather complex and challenging to assign unambiguously. To assist with the interpretation, the INS spectra are calculated accurately utilizing the quantum methodology which incorporates the coupled five-dimensional TR energy levels and wave functions of the H2 molecule confined in each type of nanocage. The computed INS spectra are highly realistic and reflect the complexity of the coupled TR dynamics of the guest H2 in the anisotropic confining environment. The simulated INS spectra of p-H2 and o-H2 in the small and medium cages are compared with the experimental data, and are indispensable for their interpretation.

Neutron scattering measurements and computation of the quantum dynamics of hydrogen molecules trapped in the small and large cages of clathrate hydrates.

We report inelastic neutron scattering (INS) measurements on molecular hydrogen trapped in simple (D2O) and binary (D2O plus perdeuterated tetrahydrofuran) clathrate hydrates, performed at a low temperature using two different neutron spectrometers to probe both energy and momentum transfer. The INS spectra of binary clathrate samples exhibit a rich structure containing sharp bands arising from both the rotational transitions and the rattling modes of the guest H2 molecule. They agree well with the rigorous fully quantum simulations, which account for the subtle effects of the anisotropy, angular and radial, of the host cage on the H2 microscopic dynamics and the resulting spectra. The simple clathrate samples present a much greater challenge, due to the multiple H2 occupancy of the large cages, which makes the quantum calculations an extremely difficult task. In addition, we discuss in detail various physical aspects of the experimental and simulated INS spectra, such as their temperature dependence, the effects of the cage geometry, and the different features associated with the ortho-hydrogen and para-hydrogen species.