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.

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.

The Effect of Embedded Nanoparticles on the Phonon Spectrum of Ice: An Inelastic X-ray Scattering Study.

As a contribution to the ongoing effort toward high-frequency sound manipulation in composite materials, we use Inelastic X-ray Scattering to probe the phonon spectrum of ice, either in a pure form or with a sparse amount of nanoparticles embedded in it. The study aims at elucidating the ability of nanocolloids to condition the collective atomic vibrations of the surrounding environment. We observe that a nanoparticle concentration of about 1 % in volume is sufficient to visibly affect the phonon spectrum of the icy substrate, mainly canceling its optical modes and adding nanoparticle phonon excitations to it. We highlight this phenomenon thanks to the lineshape modeling based on a Bayesian inference, which enables us to capture the finest detail of the scattering signal. The results of this study can empower new routes toward the shaping of sound propagation in materials through the control of their structural heterogeneity.

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.

Universal Matsubara time decay of quantum autocorrelations for Boltzmann particles.

The general properties of time dependent autocorrelations in many-body quantum systems are here analyzed at thermodynamic equilibrium in the Boltzmann canonical ensemble at temperature T, by means of the exponential expansion theory (EET). It is shown that the Kubo-Martin-Schwinger (KMS) symmetry applied to the exponential expansion of the correlation leads to the existence of two different sets of decay modes (channels) here indicated as "Matsubara modes" and "system modes," respectively. The Matsubara modes are a series of pure decay channels with time constants representing a direct action of the thermostat upon the correlation, with a characteristic principal decay time τ_{1}=ℏ/(2πk_{B}T), where ℏ and k_{B} are the Planck and Boltzmann constants, and T is the temperature. Moreover, the KMS condition implies that the amplitudes pertaining to the even and odd contribution of the system modes to the quantum correlation are not independent. These two properties are quantum mechanical in nature and "universal," in the sense that they are present for any autocorrelation of a quantum system at equilibrium at a temperature T. The Matsubara modes' contribution to the time behavior of a quantum correlation is limited to times of the order of τ_{1}, which however can be comparable with some of the characteristic decay times of the system modes. In addition, since the parameters representing the overall time behavior of the quantum correlation can be given in terms of the parameters of its Kubo transform, the EET representation turns out to be useful in calculations exploiting the outputs of some widespread quantum simulation methods. A discussion of the properties of these relations is described in detail with numerical examples. The case of the velocity autocorrelation function of para hydrogen at low temperature is also reported as a final example for a real system.

Altering Terahertz Sound Propagation in a Liquid upon Nanoparticle Immersion.

One of the grand challenges of new generation Condensed Matter physicists is the development of novel devices enabling the control of sound propagation at terahertz frequency. Indeed, phonon excitations in this frequency window are the leading conveyor of heat transfer in insulators. Their manipulation is thus critical to implementing heat management based on the structural design. To explore the possibility of controlling the damping of sound waves, we used high spectral contrast Inelastic X-ray Scattering (IXS) to comparatively study terahertz acoustic damping in a dilute suspension of 50 nm nanospheres in glycerol and on pure glycerol. Bayesian inference-based modeling of measured spectra indicates that, at sufficiently large distances, the spectral contribution of collective modes in the glycerol suspension becomes barely detectable due to the enhanced damping, the weakening, and the slight softening of the dominant acoustic mode.

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}

The damping of terahertz acoustic modes in aqueous nanoparticle suspensions.

In this work, we investigate the possibility of controlling the acoustic damping in a liquid when nanoparticles are suspended in it. To shed light on this topic, we performed Inelastic X-Ray Scattering (IXS) measurements of the terahertz collective dynamics of aqueous suspensions of nanospheres of various materials, size, and relative concentration, either charged or neutral. A Bayesian analysis of measured spectra indicates that the damping of the two acoustic modes of water increases upon nanoparticle immersion. This effect seems particularly pronounced for the longitudinal acoustic mode, which, whenever visible at all, rapidly damps off when increasing the exchanged wavevector. Results also indicate that the observed effect strongly depends on the material the immersed nanoparticles are made of.

Searching for a second excitation in the inelastic neutron scattering spectrum of a liquid metal: a Bayesian analysis.

When probed at nanometer and picosecond scales, the properties of a liquid present striking analogies with the ones of the corresponding solid, one of the most surprising is the ability of supporting shear wave propagation, as a rigid medium. Although this evidence is being reported by a growing number of terahertz scattering measurements, it remains an open question whether it is universal or rather typical of some liquids only. Furthermore, given its elusive signatures in the scattering signal, the detection of this effect appears as a typical case where an unintentional "bias of confirmation" can mislead experimentalists. We thus decided to use a Bayesian inference approach to achieve a probabilistically grounded and evidence-based lineshape modeling of the inelastic neutron scattering spectra from liquid silver, whose simulated density autocorrelations bear evidence of a shear mode propagation over very short distances. The result of our analysis indicates that the observation of any additional, non-longitudinal, acoustic modes in this simple system goes beyond the accuracy of the used scattering method.

Onset of interfacial waves in the terahertz spectrum of a nanoparticle suspension.

We used inelastic x-ray scattering to gain insight into the complex terahertz dynamics of a diluted Au-nanoparticle suspension in glycerol. We observe that, albeit sparse, Au nanoparticles leave clear signatures on the dynamic response of the system, the main one being an additional mode propagating at the nanoparticle-glycerol interface. A Bayesian inferential analysis of the line shape reveals that such a mode, at variance with conventional acoustic modes, keeps a hydrodynamiclike behavior well beyond the continuous limit and down to subnanometer distances.

A Structurally Simple Vaccine Candidate Reduces Progression and Dissemination of Triple-Negative Breast Cancer.

The Tn antigen is a well-known tumor-associated carbohydrate determinant, often incorporated in glycopeptides to develop cancer vaccines. Herein, four copies of a conformationally constrained mimetic of the antigen TnThr (GalNAc-Thr) were conjugated to the adjuvant CRM197, a protein licensed for human use. The resulting vaccine candidate, mime[4]CRM elicited a robust immune response in a triple-negative breast cancer mouse model, correlated with high frequency of CD4+ T cells and low frequency of M2-type macrophages, which reduces tumor progression and lung metastasis growth. Mime[4]CRM-mediated activation of human dendritic cells is reported, and the proliferation of mime[4]CRM-specific T cells, in cancer tissue and peripheral blood of patients with breast cancer, is demonstrated. The locked conformation of the TnThr mimetic and a proper presentation on the surface of CRM197 may explain the binding of the conjugate to the anti-Tn antibody Tn218 and its efficacy to fight cancer cells in mice.

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.

The Terahertz Dynamics of an Aqueous Nanoparticle Suspension: An Inelastic X-ray Scattering Study.

We used the high-resolution Inelastic X-ray Scattering beamline of the Advanced Photon Source at Argonne National Laboratory to measure the terahertz spectrum of pure water and a dilute aqueous suspension of 15 nm diameter spherical Au nanoparticles (Au-NPs). We observe that, despite their sparse volume concentration of about 0.5%, the immersed NPs strongly influence the collective molecular dynamics of the hosting liquid. We investigate this effect through a Bayesian inference analysis of the spectral lineshape, which elucidates how terahertz transport properties of water change upon Au-NP immersion. In particular, we observe a nearly complete disappearance of the longitudinal acoustic mode and a mildly decreased ability to support shear wave propagation.

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.

Interpreting the Terahertz Spectrum of Complex Materials: The Unique Contribution of the Bayesian Analysis.

In the last few decades, experimental studies of the terahertz spectrum of density fluctuations have considerably improved our knowledge of the mesoscopic dynamics of disordered materials, which also have imposed new demands on the data modelling and interpretation. Indeed, lineshape analyses are no longer limited to the phenomenological observation of inelastic features, as in the pioneering stage of Neutron or X-ray spectroscopy, rather aiming at the extraction from their shape of physically relevant quantities, as sound velocity and damping, relaxation times, or other transport coefficients. In this effort, researchers need to face both inherent and practical obstacles, respectively stemming from the highly damped nature of terahertz modes and the limited energy resolution, accessible kinematic region and statistical accuracy of the typical experimental outcome. To properly address these challenges, a global reconsideration of the lineshape modelling and the enforcement of evidence-based probabilistic inference is becoming crucial. Particularly compelling is the possibility of implementing Bayesian inference methods, which we illustrated here through an in-depth discussion of some results recently obtained in the analysis of Neutron and X-ray scattering results.

Model-free description of polymer-coated gold nanoparticle dynamics in aqueous solutions obtained by Bayesian analysis of neutron spin echo data.

We present a neutron spin echo study of the nanosecond dynamics of polyethylene glycol (PEG) functionalized nanosized gold particles dissolved in D_{2}O at two temperatures and two different PEG molecular weights (400D and 2000D). The analysis of the neutron spin echo data was performed by applying a Bayesian approach to the description of time correlation function decays in terms of exponential terms, recently proved to be theoretically rigorous. This approach, which addresses in a direct way the fundamental issue of model choice in any dynamical analysis, provides here a guide to the most statistically supported way to follow the decay of the intermediate scattering functions I(Q,t) by basing on statistical grounds the choice of the number of terms required for the description of the nanosecond dynamics of the studied systems. Then, the presented analysis avoids from the start resorting to a preselected framework and can be considered as model free. By comparing the results of PEG-coated nanoparticles with those obtained in PEG2000 solutions, we were able to disentangle the translational diffusion of the nanoparticles from the internal dynamics of the polymer grafted to them, and to show that the polymer corona relaxation follows a pure exponential decay in agreement with the behavior predicted by coarse grained molecular dynamics simulations and theoretical models. This methodology has one further advantage: in the presence of a complex dynamical scenario, I(Q,t) is often described in terms of the Kohlrausch-Williams-Watts function that can implicitly represent a distribution of relaxation times. By choosing to describe the I(Q,t) as a sum of exponential functions and with the support of the Bayesian approach, we can explicitly determine when a finer-structure analysis of the dynamical complexity of the system exists according to the available data without the risk of overparametrization. The approach presented here is an effective tool that can be used in general to provide an unbiased interpretation of neutron spin echo data or whenever spectroscopy techniques yield time relaxation data curves.

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.

Switching off hydrogen-bond-driven excitation modes in liquid methanol.

Hydrogen bonding plays an essential role on intermolecular forces, and consequently on the thermodynamics of materials defined by this elusive bonding character. It determines the property of a vital liquid as water as well as many processes crucial for life. The longstanding controversy on the nature of the hydrogen bond (HB) can be settled by looking at the effect of a vanishing HB interaction on the microscopic properties of a given hydrogen-bonded fluid. This task suits the capabilities of computer simulations techniques, which allow to easily switch off HB interactions. We then use molecular dynamics to study the microscopic properties of methanol, a prototypical HB liquid. Fundamental aspects of the dynamics of methanol at room temperature were contextualised only very recently and its rich dynamics was found to have striking analogies with that of water. The lower temperature (200 K) considered in the present study led us to observe that the molecular centre-of-mass dynamics is dominated by four modes. Most importantly, the computational ability to switch on and off hydrogen bonds permitted us to identify which, among these modes, have a pure HB-origin. This clarifies the role of hydrogen bonds in liquid dynamics, disclosing new research opportunities and unexplored interpretation schemes.