Experimental evidence for the long-time decay of the velocity-autocorrelation in liquid sodium

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.

Velocity autocorrelation function in supercooled liquids: Long-time tails and anomalous shear-wave propagation

Molecular dynamic simulations are performed to reveal the long-time behavior of the velocity autocorrelation function (VAF) by utilizing the finite-size effect in a Lennard-Jones binary mixture. Whereas in normal liquids the classical positive t^{-3/2} long-time tail is observed, we find in supercooled liquids a negative tail. It is strongly influenced by the transfer of the transverse current wave across the period boundary. The t^{-5/2} decay of the negative long-time tail is confirmed in the spectrum of VAF. Modeling the long-time transverse current within a generalized Maxwell model, we reproduce the negative long-time tail of the VAF, but with a slower algebraic t^{-2} decay.

Experimental evidence for a dynamical crossover in liquid aluminium

The temperature dependence of the dynamic structure factor at next-neighbour distances has been investigated for liquid aluminium. This correlation function is a sensitive parameter for changes in the local environment and its Fourier transform was measured in a coherent inelastic neutron scattering experiment. The zero frequency amplitude decreases in a nonlinear way and indicates a change in dynamics around 1.4 ∙ Tmelting. From that amplitude a generalized viscosity can be derived which is a measure of local stress correlations on next-neighbour distances. The derived generalized longitudinal viscosity shows a changing slope at the same temperature range. At this temperature the freezing out of degrees of freedom for structural relaxation upon cooling sets in which can be understood as a precursor towards the solid state. That crossover in dynamics of liquid aluminium shows the same signatures as previously observed in liquid rubidium and lead, indicating an universal character.

Transition from hydrodynamic to viscoelastic propagation of sound in molten RbBr

Inelastic neutron scattering was applied to measure the acoustic-type excitations in the molten alkali halide rubidium bromide. For molten RbBr neutron scattering is mainly sensitive to the number density fluctuation spectrum and is not influenced by charge fluctuations. Utilizing a dedicated Brillouin scattering spectrometer, we focused on the small-wave-vector range. From inelastic excitations in the spectra a dispersion relation was obtained, which shows a large positive dispersion effect. This frequency enhancement is related to a viscoelastic response of the liquid at high frequencies. Towards small wave vectors we identify the transition to hydrodynamic behavior. This observation is supported by a transition of the sound velocity from a viscoelastic enhanced value to the adiabatic speed of sound for the acoustic-type excitations. Furthermore, the spectrum transforms into a line shape compatible with a prediction from hydrodynamics.

Superdiffusion in quasi-two-dimensional Yukawa liquids

The emergence and vanishing of superdiffusion in quasi-two-dimensional Yukawa systems are investigated by molecular dynamics simulations. Using both the asymptotic behavior of the mean-squared displacement of the particles and the long-time tail of the velocity autocorrelation function as indicators of superdiffusion, we confirm the existence of a transition from normal diffusion to superdiffusion in systems changing from a three-dimensional to a two-dimensional character. A connection between superdiffusion and dimensionality is established by the behavior of the projected pair distribution function.

Stochastic dynamics of bionanosystems: Multiscale analysis and specialized ensembles

An approach for simulating bionanosystems such as viruses and ribosomes is presented. This calibration-free approach is based on an all-atom description for bionanosystems, a universal interatomic force field, and a multiscale perspective. The supramillion-atom nature of these bionanosystems prohibits the use of a direct molecular dynamics approach for phenomena such as viral structural transitions or self-assembly that develop over milliseconds or longer. A key element of these multiscale systems is the cross-talk between, and consequent strong coupling of processes over many scales in space and time. Thus, overall nanoscale features of these systems control the relative probability of atomistic fluctuations, while the latter mediate the average forces and diffusion coefficients that induce the dynamics of these nanoscale features. This feedback loop is overlooked in typical coarse-grained methods. We elucidate the role of interscale cross-talk and overcome bionanosystem simulation difficulties with (1) automated construction of order parameters (OPs) describing suprananometer scale structural features, (2) construction of OP-dependent ensembles describing the statistical properties of atomistic variables that ultimately contribute to the entropies driving the dynamics of the OPs, and (3) the derivation of a rigorous equation for the stochastic dynamics of the OPs. As the OPs capture hydrodynamic modes in the host medium, "long-time tails" in the correlation functions yielding the generalized diffusion coefficients do not emerge. Since the atomic-scale features of the system are treated statistically, several ensembles are constructed that reflect various experimental conditions. Attention is paid to the proper use of the Gibbs hypothesized equivalence of long-time and ensemble averages to accommodate the varying experimental conditions. The theory provides a basis for a practical, quantitative bionanosystem modeling approach that preserves the cross-talk between the atomic and nanoscale features. A method for integrating information from nanotechnical experimental data in the derivation of equations of stochastic OP dynamics is also introduced.

Short-time inertial response of viscoelastic fluids measured with Brownian motion and with active probes

We have directly observed short-time stress propagation in viscoelastic fluids using two optically trapped particles and a fast interferometric particle-tracking technique. We have done this both by recording correlations in the thermal motion of the particles and by measuring the response of one particle to the actively oscillated second particle. Both methods detect the vortexlike flow patterns associated with stress propagation in fluids. This inertial vortex flow propagates diffusively for simple liquids, while for viscoelastic solutions the pattern spreads superdiffusively, depending on the shear modulus of the medium.

Long-time behavior of the velocity autocorrelation function at low densities and near the critical point of simple fluids

Numerous theoretical and numerical works have been devoted to the study of the algebraic decrease at large times of the velocity autocorrelation function of particles in a fluid. The derivation of this behavior, the so-called long-time tail, generally based on linearized hydrodynamics, makes no reference to any specific characteristic of the particle interactions. However, in the literature doubts have been expressed about the possibility that by numerical simulations the long-time tail can be observed in the whole fluid phase domain of systems in which the particles interact by soft-core and attractive pair potentials. In this work, extensive and accurate molecular-dynamics simulations establish that the predicted long-time tail of the velocity autocorrelation function exists in a low-density fluid of particles interacting by a soft-repulsive potential and near the liquid-gas critical point of a Lennard-Jones system. These results contribute to the confirmation that the algebraic decay of the velocity autocorrelation function is universal in these fluid systems.

Excess entropy scaling for the diffusion coefficient in expanded liquid metals

Molecular-dynamics simulation is used to compute the pair correlation function and the velocity autocorrelation function of Cs and Rb along the liquid-vapor coexistence curve, from which the excess entropy S(ex) and the diffusion coefficient D are deduced. The numerical results of both physical properties are correlated and a scaling law between the excess entropy and the reduced diffusion coefficient D(*)(=D/D(0)) is investigated for different expressions of the reduction parameter D(0). The choice of thermodynamic states along the liquid--vapor coexistence curve gives us the possibility to extend the investigation of the relation between the reduced diffusion coefficient and the excess entropy over a wide area and to test the adequacy of the scaling law confidently.

Long-time behavior of the velocity autocorrelation function for moderately dense, soft-repulsive, and Lennard-Jones fluids

The long-time behavior of the velocity autocorrelation function (VAF), for hard disk and sphere systems, has been extensively explored. Its behavior for systems interacting via a soft repulsive or attractive potential is less well known. We explore the conditions under which the nonexponential, long-time tail in the velocity autocorrelation function of a tagged atom, in soft-repulsive sphere (Weeks-Chandler-Andersen) and Lennard-Jones atomic fluids, may be readily observed by the molecular dynamics method. The effect of changing the system size, the fluid density, the form of the interatomic force and the mass of the tagged atoms are investigated. We were able to observe this long-time tail only for systems of moderate density. At low density the effect, if it exists, is at longer times than we can currently simulate owing to limitations of system size and at higher densities these tails were not observed possibly due to other effects dominating the behavior of the VAF and masking this behavior. Under the physical conditions that are simulated here attractive forces have very little effect on the behavior of the VAF. However, as the mass of the tagged particles is increased the time at which the long-time tail commences is lengthened and its magnitude is significantly increased. This later effect suggests that by increasing the mass of the tagged particles one may be able to study more readily the behavior, nature and physical origin of long-time behavior of the VAF both by computational and by experimental techniques.