Active-parameter polydispersity in the 2d ABP Yukawa model.

In experiments and simulations of passive as well as active matter the most commonly studied kind of parameter polydispersity is that of varying particles size. This paper investigates by simulations the effects of introducing polydispersity in other parameters for two-dimensional active Brownian particles with Yukawa pair interactions. Polydispersity is studied separately in the translational and rotational diffusion coefficients, as well as in the swim velocityv0. Uniform and binary parameter distributions are considered in the homogeneous and the motility-induced phase-separation (MIPS) phases. We find only minute changes in structure and dynamics upon the introduction of parameter polydispersity, even for situations involving 50% polydispersity. The reason for this is not clear. An exception is the case ofv0polydispersity for which the average radial distribution function with changing polydispersity shows significant variations in the MIPS phase. Even in this case, however, the dynamics is only modestly affected. As a possible application of our findings, we suggest that a temporary introduction of polydispersity into a single-component active-matter model characterized by a very long equilibration time, i.e. a glass-forming active system, may be used to equilibrate the system efficiently by particle swaps.

Active-matter isomorphs in the size-polydisperse Ornstein-Uhlenbeck Lennard-Jones model.

This paper studies size-polydisperse Lennard-Jones systems described by active Ornstein-Uhlenbeck particle (AOUP) dynamics. The focus is on the existence of isomorphs (curves of invariant structure and dynamics) in the model's three-dimensional phase diagram. Isomorphs are traced out from a single steady-state configuration by means of the configurational-temperature method. Good isomorph invariance of the reduced-unit radial distribution function and the mean-square displacement as a function of time is demonstrated for three uniform-distribution polydispersities,12%, 23%, and 29%. Comparing to active-matter isomorphs generated by the analytical direct-isomorph-check method, the latter have poorer invariance of the structure, but better invariance of the dynamics. We conclude that both methods can be used to quickly get an overview of the phase diagram of polydisperse AOUP models involving a potential-energy function obeying the hidden-scale-invariance property required for isomorph theory to apply.

Configurational temperature in active matter. II. Quantifying the deviation from thermal equilibrium.

This paper proposes using the configurational temperature T_{conf} for quantifying how far an active-matter system is from thermal equilibrium. We measure this "distance" by the ratio of the systemic temperature T_{s} to T_{conf}, where T_{s} is the canonical-ensemble temperature for which the average potential energy is equal to that of the active-matter system. T_{conf} is "local" in the sense that it is the average of a function, which depends only on how the potential energy varies in the vicinity of a given configuration. In contrast, T_{s} is a global quantity. The quantity T_{s}/T_{conf} is straightforward to evaluate in a computer simulation; equilibrium simulations in conjunction with a single steady-state active-matter configuration are enough to determine T_{s}/T_{conf}. We validate the suggestion that T_{s}/T_{conf} quantifies the deviation from thermal equilibrium by data for the radial distribution function of the 3D Kob-Andersen and 2D Yukawa active-matter models with active Ornstein-Uhlenbeck and active Brownian Particle dynamics. Moreover, we show that T_{s}/T_{conf}, structure, and dynamics of the homogeneous phase are all approximately invariant along the motility-induced phase separation boundary in the phase diagram of the 2D Yukawa model. The measure T_{s}/T_{conf} is not limited to active matter and can be used for quantifying how far any system involving a potential-energy function, e.g., a driven Hamiltonian system, is from thermal equilibrium.

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.

Isomorph Invariance of Higher-Order Structural Measures in Four Lennard-Jones Systems.

In the condensed liquid phase, both single- and multicomponent Lennard-Jones (LJ) systems obey the "hidden-scale-invariance" symmetry to a good approximation. Defining an isomorph as a line of constant excess entropy in the thermodynamic phase diagram, the consequent approximate isomorph invariance of structure and dynamics in appropriate units is well documented. However, although all measures of the structure are predicted to be isomorph invariant, with few exceptions only the radial distribution function (RDF) has been investigated. This paper studies the variation along isomorphs of the nearest-neighbor geometry quantified by the occurrence of Voronoi structures, Frank-Kasper bonds, icosahedral local order, and bond-orientational order. Data are presented for the standard LJ system and for three binary LJ mixtures (Kob-Andersen, Wahnström, NiY2). We find that, while the nearest-neighbor geometry generally varies significantly throughout the phase diagram, good invariance is observed along the isomorphs. We conclude that higher-order structural correlations are no less isomorph invariant than is the RDF.

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

Combining the recent Piskulich-Thompson approach [Z. A. Piskulich and W. H. Thompson, J. Chem. Phys. 152, 011102 (2020)JCPSA60021-960610.1063/1.5135932] 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 ranges where the van't Hoff assumption that energetic and entropic forces are temperature independent is valid to a good approximation. 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 van't Hoff assumption in the whole range of temperatures can be determined in the whole liquid region of the phase diagram from a single simulation.

Erratum: Structural Origin of Enhanced Dynamics at the Surface of a Glassy Alloy [Phys. Rev. Lett. 119, 245501 (2017)].

This corrects the article DOI: 10.1103/PhysRevLett.119.245501.

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

We analyze 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 strain rates 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 strain rate: At the lower end of the unstable region, force chains merge to form giant bands that span the system in the compressional direction and propagate in the 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.

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 nonuniform 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.

Nonaffine displacements and the nonlinear response of a strained amorphous solid.

We demonstrate that irreversible structural reorganization is not necessary for the observation of yield behavior in an amorphous solid. While the majority of solids strained to their yield point do indeed undergo an irreversible reorganization, we find that a significant fraction of solids exhibits 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 nonaffine motions driven by the applied strain. The empirical similarity of the dependence of the ratio of stress over strain on the nonaffine mean-square displacement to that for the shear modulus obtained from quiescent liquid at nonzero temperature supports the proposition that rigidity depends on the size of the sampled configurational space only and is insensitive to how this space is sampled.

Role of density modulation in the spatially resolved dynamics of strongly confined liquids.

Confinement by walls usually produces a strong modulation in the density of dense liquids near the walls. Using molecular dynamics simulations, we examine the effects of the density modulation on the spatially resolved dynamics of a liquid confined between two parallel walls, using a resolution of a fraction of the interparticle distance in the liquid. The local dynamics is quantified by the relaxation time associated with the temporal autocorrelation function of the local density. We find that this local relaxation time varies in phase with the density modulation. The amplitude of the spatial modulation of the relaxation time can be quite large, depending on the characteristics of the wall and thermodynamic parameters of the liquid. To disentangle the effects of confinement and density modulation on the spatially resolved dynamics, we compare the dynamics of a confined liquid with that of an unconfined one in which a similar density modulation is induced by an external potential. We find several differences indicating that density modulation alone cannot account for all the features seen in the spatially resolved dynamics of confined liquids. We also examine how the dynamics near a wall depends on the separation between the two walls and show that the features seen in our simulations persist in the limit of large wall separation.

Rigidity in Condensed Matter and Its Origin in Configurational Constraint.

Motivated by the formal argument that a nonzero shear modulus is the result of averaging over a constrained configuration 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.

Suspensions of polymer-grafted nanoparticles with added polymers-Structure and effective pair-interactions.

We present the results of combined experimental and theoretical (molecular dynamics simulations and integral equation theory) studies of the structure and effective interactions of suspensions of polymer grafted nanoparticles (PGNPs) in the presence of linear polymers. Due to the absence of systematic experimental and theoretical studies of PGNPs, it is widely believed that the structure and effective interactions in such binary mixtures would be very similar to those of an analogous soft colloidal material-star polymers. In our study, polystyrene-grafted gold nanoparticles with functionality f = 70 were mixed with linear polystyrene (PS) of two different molecular weights for obtaining two PGNP:PS size ratios, ξ = 0.14 and 2.76 (where, ξ = Mg/Mm, Mg and Mm being the molecular weights of grafting and matrix polymers, respectively). The experimental structure factor of PGNPs could be modeled with an effective potential (Model-X), which has been found to be widely applicable for star polymers. Similarly, the structure factor of the blends with ξ = 0.14 could be modeled reasonably well, while the structure of blends with ξ = 2.76 could not be captured, especially for high density of added polymers. A model (Model-Y) for effective interactions between PGNPs in a melt of matrix polymers also failed to provide good agreement with the experimental data for samples with ξ = 2.76 and high density of added polymers. We tentatively attribute this anomaly in modeling the structure factor of blends with ξ = 2.76 to the questionable assumption of Model-X in describing the added polymers as star polymers with functionality 2, which gets manifested in both polymer-polymer and polymer-PGNP interactions especially at higher fractions of added polymers. The failure of Model-Y may be due to the neglect of possible many-body interactions among PGNPs mediated by matrix polymers when the fraction of added polymers is high. These observations point to the need for a new framework to understand not only the structural behavior of PGNPs but also possibly their dynamics and thermo-mechanical properties as well.

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.

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.

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 nonmonotonic 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 always find 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.