Extreme case of density scaling: The Weeks-Chandler-Andersen system at low temperatures

Scaling Properties of Liquid Dynamics Predicted from a Single Configuration: Pseudoisomorphs for Harmonic-Bonded Molecules

Isomorphs are curves in the thermodynamic phase diagram of invariant excess entropy, structure, and dynamics, while pseudoisomorphs are curves of invariant structure and dynamics, but not of the excess entropy. The latter curves have been shown to exist in molecular models with flexible bonds [Olsen, A. E. J. Chem. Phys. 2016, 145, 241103]. We here present three force-based methods to trace out pseudoisomorphs based on a single configuration and test them on the asymmetric dumbbell and 10-bead Lennard-Jones chain models with bonds modeled as harmonic springs. The three methods are based on requiring that particle forces, center-of-mass forces, and torques, respectively, are invariant in reduced units. For each of the two investigated models we identify a method that works well for tracing out pseudoisomorphs, but these methods are not the same. Overall, we find that the more internal degrees of freedom there are in the molecule studied, the less they appear to affect the gross dynamical behavior. Moreover, the "internal" degrees of freedom (including rotation) do not significantly affect the scaling behavior of the dynamical/transport coefficients provided some 'quenching' is performed.

A comparative study of the correlation between the structure and the dynamics for systems interacting via attractive and repulsive potentials

We present the study of the structure-dynamics correlation for systems interacting via attractive Lennard-Jones (LJ) and its repulsive counterpart, the Weeks-Chandler-Andersen (WCA) potentials. The structural order parameter (SOP) is related to the microscopic mean-field caging potential. At a particle level, the SOP shows a distribution. Although the two systems have similar pair structures, their average SOP differs. However, this difference alone is insufficient to explain the well known slowing down of the dynamics in the LJ system at low temperatures. The slowing down can be explained in terms of a stronger coupling between the SOP and the dynamics. To understand the origin of this system specific coupling, we study the difference in the microscopic structure between the hard and soft particles. We find that for the LJ system, the structural differences of the hard and soft particles are more significant and have a much stronger temperature dependence compared to the WCA system. Thus, the study suggests that attractive interaction creates more structurally different communities. This broader difference in the structural communities is probably responsible for stronger coupling between the structure and dynamics. Thus, the system specific structure-dynamics correlation, which also leads to a faster slowing down in the dynamics, appears to have a structural origin. A comparison of the predictive power of our SOP with the local energy and two body excess entropy in determining the dynamics shows that in the LJ system, the enthalpy plays a dominant role and in the WCA system, the entropy plays a dominant role, and our SOP can capture both these contributions.

The origin of the density scaling exponent for polyatomic molecules and the estimation of its value from the liquid structure

In this article, we unravel the problem of interpreting the density scaling exponent for the polyatomic molecules representing the real van der Waals liquids. Our studies show that the density scaling exponent is a weighted average of the exponents of the repulsive terms of all interatomic interactions that occur between molecules, where the potential energy of a given interaction represents its weight. It implies that potential energy is a key quantity required to calculate the density scaling exponent value for real molecules. Finally, we use the well-known method for potential energy estimation and show that the density scaling exponent could be successfully predicted from the liquid structure for fair representatives of the real systems.

Departures from perfect isomorph behavior in Lennard-Jones fluids and solids

Isomorphs are lines on a fluid or solid phase diagram along which the microstructure is invariant on affine density scaling of the molecular coordinates. Only inverse power (IP) and hard sphere potential systems are perfectly isomorphic. This work provides new theoretical tools and criteria to determine the extent of deviation from perfect isomorphicity for other pair potentials using the Lennard-Jones (LJ) system as a test case. A simple prescription for predicting isomorphs in the fluid range using the freezing line as a reference is shown to be quite accurate for the LJ system. The shear viscosity and self-diffusion coefficient scale well are calculated using this method, which enables comments on the physical significance of the correlations found previously in the literature to be made. The virial–potential energy fluctuation and the concept of an effective IPL system and exponent, n′, are investigated, particularly with reference to the LJ freezing and melting lines. It is shown that the exponent, n′, converges to the value 12 at a high temperature as ∼ T−1/2, where T is the temperature. Analytic expressions are derived for the density, temperature, and radius derivatives of the radial distribution function along an isomorph that can be used in molecular simulation. The variance of the radial distribution function and radial fluctuation function are shown to be isomorph invariant.

Hidden scale invariance in the Gay-Berne model. II. Smectic-B phase

This paper complements a previous study of the isotropic and nematic phases of the Gay-Berne liquid-crystal model [Mehri et al., Phys. Rev. E 105, 064703 (2022)2470-004510.1103/PhysRevE.105.064703] with a study of its smectic-B phase found at high density and low temperatures. We find also in this phase strong correlations between the virial and potential-energy thermal fluctuations, reflecting hidden scale invariance and implying the existence of isomorphs. The predicted approximate isomorph invariance of the physics is confirmed by simulations of the standard and orientational radial distribution functions, the mean-square displacement as a function of time, and the force, torque, velocity, angular velocity, and orientational time-autocorrelation functions. The regions of the Gay-Berne model that are relevant for liquid-crystal experiments can thus fully be simplified via the isomorph theory.

Configurational temperature in active matter. I. Lines of invariant physics in the phase diagram of the Ornstein-Uhlenbeck model

This paper shows that the configurational temperature of liquid-state theory, T_{conf}, defines an energy scale, which can be used for adjusting model parameters of active Ornstein-Uhlenbeck particle (AOUP) models in order to achieve approximately invariant structure and dynamics upon a density change. The required parameter changes are calculated from the variation of a single configuration's T_{conf} for a uniform scaling of all particle coordinates. The resulting equations are justified theoretically for models involving a potential-energy function with hidden scale invariance. The validity of the procedure is illustrated by computer simulations of the Kob-Andersen binary Lennard-Jones AOUP model, showing the existence of lines of approximate invariance of the reduced-unit radial distribution function and time-dependent mean-square displacement.

Isomorphs in sheared binary Lennard-Jones glass: Transient response

We have studied shear deformation of binary Lennard-Jones glasses to investigate the extent to which the transient part of the stress strain curves is invariant when the thermodynamic state point is varied along an isomorph. Shear deformations were carried out on glass samples of varying stability, determined by cooling rate, and at varying strain rates, at state points deep in the glass. Density changes up to and exceeding a factor of two were made. We investigated several different methods for generating isomorphs but none of the previously developed methods could generate sufficiently precise isomorphs given the large density changes and nonequilibrium situation. Instead, the temperatures for these higher densities were chosen to give state points isomorphic to the starting state point by requiring the steady-state flow stress for isomorphic state points to be invariant in reduced units. In contrast to the steady-state flow stress, we find that the peak stress on the stress strain curve is not invariant. The peak stress decreases by a few percent for each ten percent increase in density, although the differences decrease with increasing density. Analysis of strain profiles and nonaffine motion during the transient phase suggests that the root of the changes in peak stress is a varying tendency to form shear bands, with the largest tendency occurring at the lowest densities. We suggest that this reflects the effective steepness of the potential; a higher effective steepness gives a greater tendency to form shear bands.

Comparing four hard-sphere approximations for the low-temperature WCA melting line

By combining interface-pinning simulations with numerical integration of the Clausius–Clapeyron equation, we accurately determine the melting-line coexistence pressure and fluid/crystal densities of the Weeks–Chandler–Andersen system, covering four decades of temperature. The data are used for comparing the melting-line predictions of the Boltzmann, Andersen–Weeks–Chandler, Barker–Henderson, and Stillinger hard-sphere approximations. The Andersen–Weeks–Chandler and Barker–Henderson theories give the most accurate predictions, and they both work excellently in the zero-temperature limit for which analytical expressions are derived here.

Hidden scale invariance in the Gay-Berne model

This paper presents a numerical study of the Gay-Berne liquid crystal model with parameters corresponding to calamitic (rod-shaped) molecules. The focus is on the isotropic and nematic phases at temperatures above unity, where we find strong correlations between the virial and potential-energy thermal fluctuations, reflecting the hidden scale invariance symmetry. This implies the existence of isomorphs, which are curves in the thermodynamic phase diagram of approximately invariant physics. We study numerically one isomorph in the isotropic phase and one in the nematic phase. In both cases, good invariance of the dynamics is demonstrated via data for the mean-square displacement and the reduced-unit time-autocorrelation functions of the velocity, angular velocity, force, torque, and first- and second-order Legendre polynomial orientational order parameters. Deviations from isomorph invariance are observed at short times for the orientational time-autocorrelation functions, which reflects the fact that the moment of inertia is assumed to be constant and thus not isomorph-invariant in reduced units. Structural isomorph invariance is demonstrated from data for the radial distribution functions of the molecules and their orientations. For comparison, all quantities were also simulated along an isochore of similar temperature variation, in which case invariance is not observed. We conclude that the thermodynamic phase diagram of the calamitic Gay-Berne model is essentially one-dimensional in the studied regions as predicted by isomorph theory, a fact that potentially allows for simplifications of future theories and numerical studies.