Scaling of the dynamics of flexible Lennard-Jones chains. II. Effects of harmonic bonds

The previous paper [Veldhorst et al., J. Chem. Phys. 141, 054904 (2014)] demonstrated that the isomorph theory explains the scaling properties of a liquid of flexible chains consisting of ten Lennard-Jones particles connected by rigid bonds. We here investigate the same model with harmonic bonds. The introduction of harmonic bonds almost completely destroys the correlations in the equilibrium fluctuations of the potential energy and the virial. According to the isomorph theory, if these correlations are strong a system has isomorphs, curves in the phase diagram along which structure, dynamics and the excess entropy are invariant. The Lennard-Jones chain liquid with harmonic bonds does have curves in the phase diagram along which the structure and dynamics are invariant. The excess entropy is not invariant on these curves, which we refer to as "pseudoisomorphs". In particular this means that Rosenfeld's excess-entropy scaling (the dynamics being a function of excess entropy only) does not apply for the Lennard-Jones chain with harmonic bonds.Comment: 11 pages, 11 figure

Stability of supercooled binary liquid mixtures

Recently the supercooled Wahnstrom binary Lennard-Jones mixture was partially crystallized into ${\rm MgZn_2}$ phase crystals in lengthy Molecular Dynamics simulations. We present Molecular Dynamics simulations of a modified Kob-Andersen binary Lennard-Jones mixture that also crystallizes in lengthy simulations, here however by forming pure fcc crystals of the majority component. The two findings motivate this paper that gives a general thermodynamic and kinetic treatment of the stability of supercooled binary mixtures, emphasizing the importance of negative mixing enthalpy whenever present. The theory is used to estimate the crystallization time in a Kob-Andersen mixture from the crystallization time in a series of relared systems. At T=0.40 we estimate this time to be 5$\times 10^{7}$ time units ($\approx 1. ms$). A new binary Lennard-Jones mixture is proposed that is not prone to crystallization and faster to simulate than the two standard binary Lennard-Jones mixtures; this is obtained by removing the like-particle attractions by switching to Weeks-Chandler-Andersen type potentials, while maintaining the unlike-particle attraction

Thermodynamics of condensed matter with strong pressure-energy correlations

We show that for any liquid or solid with strong correlation between its $NVT$ virial and potential-energy equilibrium fluctuations, the temperature is a product of a function of excess entropy per particle and a function of density, $T=f(s)h(\rho)$. This implies that 1) the system's isomorphs (curves in the phase diagram of invariant structure and dynamics) are described by $h(\rho)/T={\rm Const.}$, 2) the density-scaling exponent is a function of density only, 3) a Gr{\"u}neisen-type equation of state applies for the configurational degrees of freedom. For strongly correlating atomic systems one has $h(\rho)=\sum_nC_n\rho^{n/3}$ in which the only non-zero terms are those appearing in the pair potential expanded as $v(r)=\sum_n v_n r^{-n}$. Molecular dynamics simulations of Lennard-Jones type systems confirm the theory

Crystallization of the Wahnstr\"om Binary Lennard-Jones Liquid

We report observation of crystallization of the glass-forming binary Lennard-Jones liquid first used by Wahnstr\"om [G. Wahnstr\"om, Phys. Rev. A 44, 3752 (1991)]. Molecular dynamics simulations of the metastable liquid on a timescale of microseconds were performed. The liquid crystallized spontaneously. The crystal structure was identified as MgZn_2. Formation of transient crystallites is observed in the liquid. The crystallization is investigate at different temperatures and compositions. At high temperature the rate of crystallite formation is the limiting factor, while at low temperature the limiting factor is growth rate. The melting temperature of the crystal is estimated to be T_m=0.93 at rho=0.82. The maximum crystallization rate of the A_2B composition is T=0.60+/-0.02.Comment: 4 pages, 4 figures; corrected typo

The Geometry of Slow Structural Fluctuations in a Supercooled Binary Alloy

The liquid structure of a glass-forming binary alloy is studied using molecular dynamics simulations. The analysis combines common neighbour analysis with the geometrical approach of Frank and Kasper to establish that the supercooled liquid contains extended clusters characterised by the same short range order as the crystal. Fluctuations in these clusters exhibit strong correlations with fluctuations in the inherent structure energy. The steep increase in the heat capacity on cooling is, thus, directly coupled to the growing fluctuations of the Frank-Kasper clusters. The relaxation of particles in the clusters dominates the slow tail of the self-intermediate scattering function

Strong pressure-energy correlations in liquids as a configuration space property: Simulations of temperature down jumps and crystallization

Computer simulations recently revealed that several liquids exhibit strong correlations between virial and potential energy equilibrium fluctuations in the NVT ensemble [U. R. Pedersen {\it et al.}, Phys. Rev. Lett. {\bf 100}, 015701 (2008)]. In order to investigate whether these correlations are present also far from equilibrium constant-volume aging following a temperature down jump from equilibrium was simulated for two strongly correlating liquids, an asymmetric dumbbell model and Lewis-Wahnstr{\"o}m OTP, as well as for SPC water that is not strongly correlating. For the two strongly correlating liquids virial and potential energy follow each other closely during the aging towards equilibrium. For SPC water, on the other hand, virial and potential energy vary with little correlation as the system ages towards equilibrium. Further proof that strong pressure-energy correlations express a configuration space property comes from monitoring pressure and energy during the crystallization (reported here for the first time) of supercooled Lewis-Wahnstr{\"o}m OTP at constant temperature

NVU dynamics. I. Geodesic motion on the constant-potential-energy hypersurface

An algorithm is derived for computer simulation of geodesics on the constant potential-energy hypersurface of a system of N classical particles. First, a basic time-reversible geodesic algorithm is derived by discretizing the geodesic stationarity condition and implementing the constant potential energy constraint via standard Lagrangian multipliers. The basic NVU algorithm is tested by single-precision computer simulations of the Lennard-Jones liquid. Excellent numerical stability is obtained if the force cutoff is smoothed and the two initial configurations have identical potential energy within machine precision. Nevertheless, just as for NVE algorithms, stabilizers are needed for very long runs in order to compensate for the accumulation of numerical errors that eventually lead to "entropic drift" of the potential energy towards higher values. A modification of the basic NVU algorithm is introduced that ensures potential-energy and step-length conservation; center-of-mass drift is also eliminated. Analytical arguments confirmed by simulations demonstrate that the modified NVU algorithm is absolutely stable. Finally, simulations show that the NVU algorithm and the standard leap-frog NVE algorithm have identical radial distribution functions for the Lennard-Jones liquid

Exponential distributions of collective flow-event properties in viscous liquid dynamics

We study the statistics of flow events in the inherent dynamics in supercooled two- and three-dimensional binary Lennard-Jones liquids. Distributions of changes of the collective quantities energy, pressure and shear stress become exponential at low temperatures, as does that of the event "size" $S\equiv\sum {d_i}^2$. We show how the $S$-distribution controls the others, while itself following from exponential tails in the distributions of (1) single particle displacements $d$, involving a Lindemann-like length $d_L$ and (2) the number of active particles (with $d>d_L$).Comment: Accepter version (PRL

Estimating the density scaling exponent of viscous liquids from specific heat and bulk modulus data

It was recently shown by computer simulations that a large class of liquids exhibits strong correlations in their thermal fluctuations of virial and potential energy [Pedersen et al., Phys. Rev. Lett. 100, 015701 (2008)]. Among organic liquids the class of strongly correlating liquids includes van der Waals liquids, but excludes ionic and hydrogen-bonding liquids. The present note focuses on the density scaling of strongly correlating liquids, i.e., the fact their relaxation time tau at different densities rho and temperatures T collapses to a master curve according to the expression tau propto F(rho^gamma/T) [Schroder et al., arXiv:0803.2199]. We here show how to calculate the exponent gamma from bulk modulus and specific heat data, either measured as functions of frequency in the metastable liquid or extrapolated from the glass and liquid phases to a common temperature (close to the glass transition temperature). Thus an exponent defined from the response to highly nonlinear parameter changes may be determined from linear response measurements