Solid-liquid transition in polydisperse Lennard-Jones systems

NVU view on energy polydisperse Lennard-Jones systems

When energy polydispersity is introduced into the Lennard-Jones (LJ) system, there is little effect on structure and dynamics [T. S. Ingebrigtsen and J. C. Dyre, J. Phys. Chem. B 127, 2837 (2023)10.1021/acs.jpcb.3c00346]. For instance, at a given state point both the radial distribution function and the mean-square displacement as a function of time are virtually unaffected by even large energy polydispersity, which is in stark contrast to what happens when size polydispersity is introduced. We here argue-and validate by simulations of up to 30% polydispersity-that this almost invariance of structure and dynamics reflects an approximate invariance of the constant-potential-energy surface. Because NVU dynamics defined as geodesic motion at constant potential energy is equivalent to Newtonian dynamics in the thermodynamic limit, the approximate invariance of the constant-potential-energy surface implies virtually the same physics of energy polydisperse LJ systems as of the standard single-component version. In contrast, the constant-potential-energy surface is significantly affected by introducing size polydispersity.

Dynamic heterogeneity in polydisperse systems: A comparative study of the role of local structural order parameter and particle size

In polydisperse systems, describing the structure and any structural order parameter (SOP) is not trivial as it varies with the number of species we use to describe the system, M. Depending on the degree of polydispersity, there is an optimum value of M = M0 where we show that the mutual information of the system increases. However, surprisingly, the correlation between a recently proposed SOP and the dynamics is highest for M = 1. This effect increases with polydispersity. We find that the SOP at M = 1 is coupled with the particle size, σ, and this coupling increases with polydispersity and decreases with an increase in M. Careful analysis shows that at lower polydispersities, the SOP is a good predictor of the dynamics. However, at higher polydispersity, the dynamics is strongly dependent on σ. Since the coupling between the SOP and σ is higher for M = 1, it appears to be a better predictor of the dynamics. We also study the Vibrality, an order parameter independent of structural information. Compared to SOP, at high polydispersity, we find Vibrality to be a marginally better predictor of the dynamics. However, this high predictive power of Vibrality, which is not there at lower polydispersity, appears to be due to its stronger coupling with σ. Therefore, our study suggests that for systems with high polydispersity, the correlation of any order parameter and σ will affect the correlation between the order parameter and dynamics and need not project a generic predictive power of the order parameter.

Even Strong Energy Polydispersity Does Not Affect the Average Structure and Dynamics of Simple Liquids

Size-polydisperse liquids have become standard models for avoiding crystallization, thereby enabling studies of supercooled liquids and glasses formed, e.g., by colloidal systems. Purely energy-polydisperse liquids have been studied much less, but provide an interesting alternative. We here study numerically the difference in structure and dynamics obtained by introducing these two kinds of polydispersity into systems of particles interacting via the Lennard-Jones and EXP pair potentials. To a very good approximation, the average pair structure and dynamics are unchanged even for strong energy polydispersity, which is not the case for size-polydisperse systems. When the system at extreme energy polydispersity undergoes a continuous phase separation into lower and higher particle-energy regions whose structure and dynamics are different from the average, the average structure and dynamics are still virtually the same as for the monodisperse system. Our findings are consistent with the fact that the distribution of forces on the individual particles do not change when energy polydispersity is introduced, while they do change in the case of size polydispersity. A theoretical explanation remains to be found, however.

Effective structure of a system with continuous polydispersity

In a system of N particles, with continuous size polydispersity, there exists an N(N - 1) number of partial structure factors, making it analytically less tractable. A common practice is to treat the system as an effective one component system, which is known to exhibit an artificial softening of the structure. The aim of this study is to describe the system in terms of M pseudospecies such that we can avoid this artificial softening but, at the same time, have a value of M ≪ N. We use potential energy and pair excess entropy to estimate an optimum number of species, M₀. We then define the maximum width of polydispersity, Δσ₀, that can be treated as a monodisperse system. We show that M₀ depends on the degree and type of polydispersity and also on the nature of the interaction potential, whereas Δσ₀ weakly depends on the type of polydispersity but shows a stronger dependence on the type of interaction potential. Systems with a softer interaction potential have a higher tolerance with respect to polydispersity. Interestingly, M₀ is independent of system size, making this study more relevant for bigger systems. Our study reveals that even 1% polydispersity cannot be treated as an effective monodisperse system. Thus, while studying the role of polydispersity by using the structure of an effective one component system, care must be taken in decoupling the role of polydispersity from that of the artificial softening of the structure.

Hidden Scale Invariance in Polydisperse Mixtures of Exponential Repulsive Particles

Polydisperse systems of particles interacting by the purely repulsive exponential (EXP) pair potential are studied in regard to how structure and dynamics vary along isotherms, isochores, and isomorphs. The sizable size polydispersities of 23%, 29%, 35%, and 40%, as well as energy polydispersity 35%, were considered. For each system an isomorph was traced out covering about one decade in density. For all systems studied, the structure and dynamics vary significantly along the isotherms and isochores but are invariant to a good approximation along the isomorphs. We conclude that the single-component EXP system's hidden scale invariance (implying isomorph invariance of structure and dynamics) is maintained even when a sizable polydispersity is introduced into the system.

Effect of Liquid State Organization on Nanostructure and Strength of Model Multicomponent Solids

When a multicomponent liquid composed of particles with random interactions is slowly cooled below the freezing temperature, the fluid reorganizes in order to increase (decrease) the number of strong (weak) attractive interactions and solidifies into a structure composed of domains of strongly and of weakly interacting particles. Using Langevin dynamics simulations of a model system we find that the tensile strength, mode of fracture, and thermal stability of such solids differ from those of one-component solids and that these properties can be controlled by the method of preparation.

Effect of Energy Polydispersity on the Nature of Lennard-Jones Liquids

In the companion paper [ Ingebrigtsen , T. S. ; Tanaka , H. J. Phys. Chem. B 2015 , 119 , 11052 ] the effect of size polydispersity on the nature of Lennard-Jones (LJ) liquids, which represent most molecular liquids without hydrogen bonds, was studied. More specifically, it was shown that even highly size polydisperse LJ liquids are Roskilde-simple (RS) liquids. RS liquids are liquids with strong correlation between constant volume equilibrium fluctuations of virial and potential energy and are simpler than other types of liquids. Moreover, it was shown that size polydisperse LJ liquids have isomorphs to a good approximation. Isomorphs are curves in the phase diagram of RS liquids along which structure, dynamics, and some thermodynamic quantities are invariant in dimensionless (reduced) units. In this paper, we study the effect of energy polydispersity on the nature of LJ liquids. We show that energy polydisperse LJ liquids are RS liquids. However, a tendency of particle segregation, which increases with the degree of polydispersity, leads to a loss of strong virial-potential energy correlation but is mitigated by increasing temperature and/or density. Isomorphs are a good approximation also for energy polydisperse LJ liquids, although particle-resolved quantities display a somewhat poorer scaling compared to the mean quantities along the isomorph.

Effect of Size Polydispersity on the Nature of Lennard-Jones Liquids

Polydisperse fluids are encountered everywhere in biological and industrial processes. These fluids naturally show a rich phenomenology exhibiting fractionation and shifts in critical point and freezing temperatures. We study here the effect of size polydispersity on the basic nature of Lennard-Jones (LJ) liquids, which represent most molecular liquids without hydrogen bonds, via two- and three-dimensional molecular dynamics computer simulations. A single-component liquid constituting spherical particles and interacting via the LJ potential is known to exhibit strong correlations between virial and potential energy equilibrium fluctuations at constant volume. This correlation significantly simplifies the physical description of the liquid, and these liquids are now known as Roskilde-simple (RS) liquids. We show that this simple nature of the single-component LJ liquid is preserved even for very high polydispersities (above 40% polydispersity for the studied uniform distribution). We also investigate isomorphs of moderately polydisperse LJ liquids. Isomorphs are curves in the phase diagram of RS liquids along which structure, dynamics, and some thermodynamic quantities are invariant in dimensionless units. We find that isomorphs are a good approximation even for polydisperse LJ liquids. The theory of isomorphs thus extends readily to size polydisperse fluids and can be used to improve even further the understanding of these intriguing systems.