Hydrodynamic anomalies in supercritical fluid

Properties of aqueous 1,4-dioxane solution via molecular dynamics

Polyethers are promising compounds for the creation of electrochemical energy storage systems. The molecular dynamics method can facilitate the search of compounds that have the most potential. However, the application of this method requires verification of the force fields. We perform molecular dynamics calculations of the physical properties of the aqueous 1,4-dioxane solution (density, enthalpy of mixing, and viscosity) and compare them to the available experimental data. In addition, we confirm the idea that the solution structure depends on the dioxane molar fraction, proposed in the experiment of Takamuku et al. [J. Mol. Liq. 83(1-3), 163-177 (1999)]. The hydrogen bonds between dioxane and water are analyzed. The correlation between the excess viscosity and enthalpy of mixing is demonstrated.

Transport coefficients of model lubricants up to 400 MPa from molecular dynamics

In this paper, the predictive power of molecular dynamics methods is demonstrated for the cases of model paraffinic and aromatic lubricant liquids at pressures up to 400 MPa. The shear viscosity and self-diffusion coefficients are calculated for 2,2,4-trimethylpentane (C₈H₁₈) at 298 K and 1,1-diphenylethane (C₁₄H₁₄) at 333 K. Three force fields with different levels of accuracy are compared by the ability to predict the experimental data. The Stokes-Einstein correlation between viscosity and self-diffusion is demonstrated for both compounds.

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.

Thermodynamics, dynamics, and structure of supercritical water at extreme conditions

Molecular dynamics (MD) simulations to understand the thermodynamic, dynamic, and structural changes in supercritical water across the Frenkel line and the melting line have been performed.

Structural behavior of fluids from the vapor and liquid region to the supercritical phase

A metric (χ) is introduced to quantify the relative proportion of particles having a specified number of near neighbors that are characteristic of liquid-phase properties. It can be used as a simple alternative to other methods for the investigation of some aspects of percolation behavior. Values of χ are obtained from molecular-dynamics simulations spanning the heterogeneous vapor and liquid region and the supercritical phase of the Lennard-Jones fluid. The supercritical phase can be delineated into regions of different structural properties. At different isochoric subcritical conditions, the temperature versus χ behavior shows evidence of inflections, which are associated with the onset of transitions from the vapor and liquid region to the supercritical phase. The analysis suggests a phenomenological requirement for the critical point in terms of a near-equal proportion of near neighbors with gaslike and liquidlike characteristics.

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.

A simple closure procedure for the study of velocity autocorrelation functions in fluids as a "bridge" between different theoretical approaches

Velocity autocorrelation functions (VAFs) of the fluids are studied on short- and long-time scales within a unified approach. This approach is based on an effective summation of the infinite continued fraction at a reasonable assumption about convergence of relaxation times of the high order memory functions, which have a purely kinetic origin. The VAFs obtained within our method are compared with computer simulation data for the liquid Ne at different densities and the results, which follow from the Markovian approximation for the highest order kinetic kernels. It is shown that in all the thermodynamic points and at the chosen level of the hierarchy, our results agree much better with the molecular dynamic data than those of the Markovian approximation. The density dependence of the transition time, needed for the fluid to attain the hydrodynamic stage of evolution, is evaluated. The common and distinctive features of our method are discussed in their relations to the generalized collective mode theory, the mode coupling theory, and some other theoretical approaches.

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.

Self-consistent molecular dynamics calculation of diffusion in higher n-alkanes

Diffusion is one of the key subjects of molecular modeling and simulation studies. However, there is an unresolved lack of consistency between Einstein-Smoluchowski (E-S) and Green-Kubo (G-K) methods for diffusion coefficient calculations in systems of complex molecules. In this paper, we analyze this problem for the case of liquid n-triacontane. The non-conventional long-time tails of the velocity autocorrelation function (VACF) are found for this system. Temperature dependence of the VACF tail decay exponent is defined. The proper inclusion of the long-time tail contributions to the diffusion coefficient calculation results in the consistency between G-K and E-S methods. Having considered the major factors influencing the precision of the diffusion rate calculations in comparison with experimental data (system size effects and force field parameters), we point to hydrogen nuclear quantum effects as, presumably, the last obstacle to fully consistent n-alkane description.

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.

Self-assembly of the decagonal quasicrystalline order in simple three-dimensional systems

Using molecular dynamics simulations we show that a one-component system can be driven to a three-dimensional decagonal (10-fold) quasicrystalline state just by purely repulsive, isotropic and monotonic interaction pair potential with two characteristic length scales; no attraction is needed. We found that self-assembly of a decagonal quasicrystal from a fluid can be predicted by two dimensionless effective parameters describing the fluid structure. We demonstrate stability of the results under changes of the potential by obtaining the decagonal order for three particle systems with different interaction potentials, both purely repulsive and attractive, but with the same values of the effective parameters. Our results suggest that soft matter quasicrystals with decagonal symmetry can be experimentally observed for the same systems demonstrating the dodecagonal order for an appropriate tuning of the effective parameters.

Thermophysical properties of supercritical water and bond flexibility

Molecular dynamics results are reported for the thermodynamic properties of supercritical water using examples of both rigid (TIP4P/2005) and flexible (TIP4P/2005f) transferable interaction potentials. Data are reported for pressure, isochoric and isobaric heat capacities, the thermal expansion coefficient, isothermal and adiabatic compressibilities, Joule-Thomson coefficient, speed of sound, self-diffusion coefficient, viscosities, and thermal conductivity. Many of these properties have unusual behavior in the supercritical phase such as maximum and minimum values. The effectiveness of bond flexibility on predicting these properties is determined by comparing the results to experimental data. The influence of the intermolecular potential on these properties is both variable and state point dependent. In the vicinity of the critical density, the rigid and flexible potentials yield very different values for the compressibilities, heat capacities, and thermal expansion coefficient, whereas the self-diffusion coefficient, viscosities, and thermal conductivities are much less potential dependent. Although the introduction of bond flexibility is a computationally expedient way to improve the accuracy of an intermolecular potential, it can be counterproductive in some cases and it is not an adequate replacement for incorporating the effects of polarization.