Structure of the Lennard-Jones liquid estimated from a single simulation

Combining the recent Piskulich-Thompson approach [Z. A. Piskulich and W. H. Thompson, {it J. Chem. Phys.} {bf 152}, 011102 (2020)] with isomorph theory, from a single simulation, the structure of a single-component Lennard-Jones (LJ) system is obtained at an arbitrary state point in almost the whole liquid region of the temperature-density phase diagram. The LJ system exhibits two temperature range where the vant Hoffs assumption that energetic and entropic forces are temperature independent is valid. A method to evaluate the structure at an arbitrary state point along an isochore from the knowledge of structures at two temperatures on the isochore is also discussed. We argue that, in general, the structure of any hidden scale-invariant system obeying the vant Hoffs assumption in the whole range of temperatures can be determined in the whole liquid region of the phase diagram from only a single simulation.

Unsteady flow, clusters and bands in a model shear-thickening fluid

We analyse the flow curves of a two-dimensional assembly of granular particles which are interacting via frictional contact forces. For packing fractions slightly below jamming, the fluid undergoes a large scale instability, implying a range of stress and strainrates where no stationary flow can exist. Whereas small systems were shown previously to exhibit hysteretic jumps between the low and high stress branches, large systems exhibit continuous shear thickening arising from averaging unsteady, spatially heterogeneous flows. The observed large scale patterns as well as their dynamics are found to depend on strainrate: At the lower end of the unstable region, force chains merge to form giant bands that span the system in compressional direction and propagate in dilational direction. At the upper end, we observe large scale clusters which extend along the dilational direction and propagate along the compressional direction. Both patterns, bands and clusters, come in with infinite correlation length similar to the sudden onset of system-spanning plugs in impact experiments.

The Structural Origin of Enhanced Dynamics at the Surface of a Glassy Alloy

The enhancement of mobility at the surface of an amorphous alloy is studied using a combination of molecular dynamic simulations and normal mode analysis of the non-uniform distribution of Debye-Waller factors. The increased mobility at the surface is found to be associated with the appearance of Arrhenius temperature dependence. We show that the transverse Debye-Waller factor exhibits a peak at the surface. Over the accessible temperature range, we find that the bulk and surface diffusion coefficients obey the same empirical relationship with the respective Debye-Waller factors. Extrapolating this relationship to lower T, we argue that the observed decrease in the constraint at the surface is sufficient to account for the experimentally observed surface enhancement of mobility.

Non-Affine Displacements and the Non-Linear Response of a Strained Amorphous Solid

We demonstrate that irreversible structural reorganization is not necessary for the observation of yield behaviour in an amorphous solid. While the majority of solids strained to their yield point do indeed undergo an irreversible reorganization, we find a significant fraction of solids exhibit yield via a reversible strain. We also demonstrate that large instantaneous strains in excess of the yield stress can result in complete stress relaxation, a result of the large non-affine motions driven by the applied strain. The empirical similarity of the dependence of the ratio of stress over strain on the non-affine mean squared displacement with that for the shear modulus obtained from quiescent liquid at non-zero temperature supports the proposition that rigidity depends on the size of the sampled configurational space only, and is insensitive as to how this space is sampled.

The Geometric Mean Squared Displacement and the Stokes-Einstein Scaling in a Supercooled Liquid

It is proposed that the rate of relaxation in a liquid is better described by the geometric mean of the van Hove distribution function, rather than the standard arithmetic mean used to obtain the mean squared displacement. The difference between the two means is shown to increase significantly with an increase in the non-Gaussian character of the displacement distribution. Preliminary results indicate that the geometric diffusion constant results in a substantial reduction of the deviation from Stokes-Einstein scaling.

Rigidity in Condensed Matter and its Origin in Configurational Constraint

Motivated by the formal argument that a non-zero shear modulus is the result of averaging over a constrained configurations space, we demonstrate that the shear modulus calculated over a range of temperatures and averaging times can be expressed (relative to its infinite frequency value) as a single function of the mean squared displacement. This result is shown to hold for both a glass-liquid and a crystal-liquid system.

Spatial modulation of the composition of a binary liquid near a repulsive wall

When a binary liquid is confined by a strongly repulsive wall, the local density is depleted near the wall and an interface similar to that between the liquid and its vapor is formed. This analogy suggests that the composition of the binary liquid near this interface should exhibit spatial modulation similar to that near a liquid-vapor interface even if the interactions of the wall with the two components of the liquid are the same. The Guggenheim adsorption relation quantifies the concentrations of two components of a binary mixture near a liquid-vapor interface and qualitatively states that the majority (minority) component enriches the interface for negative (positive) mixing energy if the surface tensions of the two components are not very different. From molecular dynamics simulations of binary mixtures with different compositions and interactions, we find that the Guggenheim relation is qualitatively satisfied at wall-induced interfaces for systems with negative mixing energy at all state points considered. For systems with positive mixing energy, this relation is found to be qualitatively valid at low densities, while it is violated at state points with high density where correlations in the liquid are strong. This observation is validated by a calculation of the density profiles of the two components of the mixture using density functional theory with the Ramakrishnan-Yussouff free-energy functional. Possible reasons for the violation of the Guggenheim relation are discussed.

Liquid-liquid critical point in supercooled silicon

A novel liquid-liquid phase transition has been proposed and investigated in a wide variety of pure substances recently, including water, silica and silicon. From computer simulations using the Stillinger-Weber classical empirical potential, Sastry and Angell [1] demonstrated a first order liquid-liquid transition in supercooled silicon, subsequently supported by experimental and simulation studies. Here, we report evidence for a liquid-liquid critical end point at negative pressures, from computer simulations using the SW potential. Compressibilities exhibit a growing maximum upon lowering temperature below 1500 K and isotherms exhibit density discontinuities below 1120 K, at negative pressure. Below 1120 K, isotherms obtained from constant volume-temperature simulations exhibit non-monotonic, van der Waals-like behavior signaling a first order transition. We identify Tc ~ 1120 +/- 12 K, Pc -0.60 +/- 0.15 GPa as the critical temperature and pressure for the liquid-liquid critical point. The structure of the liquid changes dramatically upon decreasing the temperature and pressure. Diffusivities vary over 4 orders of magnitude, and exhibit anomalous pressure dependence near the critical point. A strong relationship between local geometry quantified by the coordination number, and diffusivity, is seen, suggesting that atomic mobility in both low and high density liquids can usefully be analyzed in terms of defects in the tetrahedral network structure. We have constructed the phase diagram of supercooled silicon. We identify the lines of compressibility, density extrema (maxima and minima) and the spinodal which reveal the interconnection between thermodynamic anomalies and the phase behaviour of the system as suggested in previous works [2-9]

Computer Simulation Study of the Phase Behavior and Structural Relaxation in a Gel-Former Modeled by Three Body Interactions

We report a computer simulation study of a model gel-former obtained by modifying the three-body interactions of the Stillinger-Weber potential for silicon. This modification reduces the average coordination number and consequently shifts the liquid-gas phase coexistence curve to low densities, thus facilitating the formation of gels without phase separation. At low temperatures and densities, the structure of the system is characterized by the presence of long linear chains interconnected by a small number of three coordinated junctions at random locations. At small wave-vectors the static structure factor shows a non-monotonic dependence on temperature, a behavior which is due to the competition between the percolation transition of the particles and the stiffening of the formed chains. We compare in detail the relaxation dynamics of the system as obtained from molecular dynamics with the one obtained from Monte Carlo dynamics. We find that the bond correlation function displays stretched exponential behavior at moderately low temperatures and densities, but exponential relaxation at low temperatures. The bond lifetime shows an Arrhenius behavior, independent of the microscopic dynamics. For the molecular dynamics at low temperatures, the mean squared displacement and the (coherent and incoherent) intermediate scattering function display at intermediate times a dynamics with ballistic character and we show that this leads to compressed exponential relaxation. For the Monte Carlo dynamics we find always an exponential or stretched exponential relaxation. Thus we conclude that the compressed exponential relaxation observed in experiments is due to the out-of-equilibrium dynamics.

Structural Relaxation of a Gel Modeled by Three Body Interactions

We report a molecular dynamics simulation study of a model gel whose interaction potential is obtained by modifying the three body Stillinger-Weber model potential for silicon. The modification reduces the average coordination number, and suppresses the liquid-gas phase coexistence curve. The low density, low temperature equilibrium gel that can thus form exhibits interesting dynamical behavior, including compressed exponential relaxation of density correlations. We show that motion responsible for such relaxation has ballistic character, and arises from the motion of chain segments in the gel without the restructuring of the gel network.