Phase behavior of Lennard-Jones particles in two dimensions

Comparison of integral equation theories of the liquid state

The Ornstein-Zernike equation is a powerful tool in liquid state theory for predicting structural and thermodynamic properties of fluids. Combined with a suitable closure, it has been shown to reproduce, e.g., the static structure factor, pressure, and compressibility of liquids to a great degree of accuracy. However, out of the multitude of closures that exist for the Ornstein-Zernike equation, it is hard to predict a priori which closure will give the most accurate predictions for the system at hand. To alleviate this problem, we compare the predictive power of many closures on a curated set of representative benchmark systems, including those with hard-sphere, inverse power-law, Gaussian core, and Lennard-Jones particles, in three and two dimensions. For example, we find that the well-known and highly used Percus-Yevick closure gives significantly worse predictions than lesser-known closures of equal complexity in all cases studied. We anticipate that the trends observed in our results will aid in making more informed decisions regarding closure choices. To facilitate the adoption of more modern closure theories, we also have packaged, documented, and distributed the code necessary to numerically solve the equations for a given closure and pair interaction potential.

A critical edge number revealed for phase stabilities of two-dimensional ball-stick polygons

Phase behaviours of two-dimensional (2D) systems constitute a fundamental topic in condensed matter and statistical physics. Although hard polygons and interactive point-like particles are well studied, the phase behaviours of more realistic molecular systems considering intermolecular interaction and molecular shape remain elusive. Here we investigate by molecular dynamics simulation phase stabilities of 2D ball-stick polygons, serving as simplified models for molecular systems. Below the melting temperature Tm, we identify a critical edge numbern c = 4 , at which a distorted square lattice emerges; when n < n c , the triangular system stabilizes at a spin-ice-like glassy state; when n > n c , the polygons stabilize at crystalline states. Moreover, in the crystalline state, Tm is higher for polygons with more edges at higher pressures but exhibits a crossover for hexagon and octagon at low pressures. A theoretical framework taking into account the competition between entropy and enthalpy is proposed to provide a comprehensive understanding of our results, which is anticipated to facilitate the design of 2D materials.

Phase Coexistence and Edge Currents in the Chiral Lennard-Jones Fluid

We study a model chiral fluid in two dimensions composed of Brownian disks interacting via a Lennard-Jones potential and a nonconservative transverse force, mimicking colloids spinning at a given rate. The system exhibits a phase separation between a chiral liquid and a dilute gas phase that can be characterized using a thermodynamic framework. We compute the equations of state and show that the surface tension controls interface corrections to the coexisting pressure predicted from the equal-area construction. Transverse forces increase surface tension and generate edge currents at the liquid-gas interface. The analysis of these currents shows that the rotational viscosity introduced in chiral hydrodynamics is consistent with microscopic bulk mechanical measurements. Chirality can also break the solid phase, giving rise to a dense fluid made of rotating hexatic patches. Our Letter paves the way for the development of the statistical mechanics of chiral particles assemblies.

Liquid-hexatic transition for soft disks

We study the liquid-hexatic transition of soft disks with massively parallel simulations and determine the equation of state as a function of system size. For systems with interactions decaying as the inverse mth power of the separation, the liquid-hexatic phase transition is continuous for m=12 and m=8, while it is of first order for m=24. The critical power m for the transition between continuous and first-order behavior is larger than previously reported. The continuous transition for m=12 implies that the two-dimensional Lennard-Jones model has a continuous liquid-hexatic transition at high temperatures. We also study the Weeks-Chandler-Andersen model and find a continuous transition at high temperatures that is consistent with the soft-disk case for m=12. Pressure data as well as our implementation are available from an open-source repository.

Collective behavior of passive and active circle swimming particle mixtures

We present a numerical study on a binary mixture of passive and circle swimming, self-propelling particles which interact via the Lennard-Jones (LJ) potential in two dimensions. Using Brownian Dynamics (BD) simulations, we present state diagrams using the control parameters such as attraction strength, angular velocity, self-propulsion velocity and composition. In a symmetric mixture, the system undergoes a transition from a mixed gel to a rotating passive cluster state and finally to a homogeneous fluid state as translational activity increases. The formation of the rotating cluster of passive particles surrounded by active and passive monomers is attributed to the combined effect of composition, activity and strength of attraction of the active particles. Different phases are characterized using radial distribution functions, bond order parameters, cluster fraction and probability distribution of local volume fractions. The present study addresses comprehensively the intricate role of activity, angular velocity, inter-particle interaction and compositional variation on the phase behavior. The predictions presented in the study can be experimentally realized in synthetic colloidal swimmers and motile bacterial suspensions.

Collective dynamics and shattering of disturbed two-dimensional Lennard-Jones crystals

Elucidating collective dynamics in crystalline systems is a common scientific question in multiple fields. In this work, by combination of high-precision numerical approach and analytical normal mode analysis, we systematically investigate the dynamical response of two-dimensional Lennard-Jones crystal as a purely classical mechanical system under random disturbance of varying strength, and reveal rich microscopic dynamics. Specifically, we observe highly symmetric velocity field composed of sharply divided coherent and disordered regions, and identify the order-disorder dynamical transition of the velocity field. Under stronger disturbance, we reveal the vacancy-driven shattering of the crystal. This featured disruption mode is fundamentally different from the dislocation-unbinding scenario in two-dimensional melting. We also examine the healing dynamics associated with vacancies of varying size. The results in this work advance our understanding about the formation of collective dynamics and crystal disruption, and may have implications in elucidating relevant non-equilibrium behaviors in a host of crystalline systems. Microscopic dynamics and underlying topological defects in the disruption of the Lennard-Jones lattice under random disturbance.

Precise determination of pair interactions from pair statistics of many-body systems in and out of equilibrium

The determination of the pair potential v(r) that accurately yields an equilibrium state at positive temperature T with a prescribed pair correlation function g_{2}(r) or corresponding structure factor S(k) in d-dimensional Euclidean space R^{d} is an outstanding inverse statistical mechanics problem with far-reaching implications. Recently, Zhang and Torquato [Phys. Rev. E 101, 032124 (2020)2470-004510.1103/PhysRevE.101.032124] conjectured that any realizable g_{2}(r) or S(k) corresponding to a translationally invariant nonequilibrium system can be attained by a classical equilibrium ensemble involving only (up to) effective pair interactions. Testing this conjecture for nonequilibrium systems as well as for nontrivial equilibrium states requires improved inverse methodologies. We have devised an optimization algorithm to precisely determine effective pair potentials that correspond to pair statistics of general translationally invariant disordered many-body equilibrium or nonequilibrium systems at positive temperatures. This methodology utilizes a parameterized family of pointwise basis functions for the potential function whose initial form is informed by small-, intermediate- and large-distance behaviors dictated by statistical-mechanical theory. Subsequently, a nonlinear optimization technique is utilized to minimize an objective function that incorporates both the target pair correlation function g_{2}(r) and structure factor S(k) so that the small intermediate- and large-distance correlations are very accurately captured. To illustrate the versatility and power of our methodology, we accurately determine the effective pair interactions of the following four diverse target systems: (1) Lennard-Jones system in the vicinity of its critical point, (2) liquid under the Dzugutov potential, (3) nonequilibrium random sequential addition packing, and (4) a nonequilibrium hyperuniform "cloaked" uniformly randomized lattice. We found that the optimized pair potentials generate corresponding pair statistics that accurately match their corresponding targets with total L_{2}-norm errors that are an order of magnitude smaller than that of previous methods. The results of our investigation lend further support to the Zhang-Torquato conjecture. Furthermore, our algorithm will enable one to probe systems with identical pair statistics but different higher-body statistics, which will shed light on the well-known degeneracy problem of statistical mechanics.

The role of attraction in the phase diagrams and melting scenarios of generalized 2D Lennard-Jones systems

Monolayer and two-dimensional (2D) systems exhibit rich phase behavior, compared with 3D systems, in particular, due to the hexatic phase playing a central role in melting scenarios. The attraction range is known to affect critical gas–liquid behavior (liquid–liquid in protein and colloidal systems), but the effect of attraction on melting in 2D systems remains unstudied systematically. Here, we have revealed how the attraction range affects the phase diagrams and melting scenarios in a 2D system. Using molecular dynamics simulations, we have considered the generalized Lennard-Jones system with a fixed repulsion branch and different power indices of attraction from long-range dipolar to short-range sticky-sphere-like. A drop in the attraction range has been found to reduce the temperature of the gas–liquid critical point, bringing it closer to the gas–liquid–solid triple point. At high temperatures, attraction does not affect the melting scenario that proceeds through the cascade of solid–hexatic (Berezinskii–Kosterlitz–Thouless) and hexatic–liquid (first-order) phase transitions. In the case of dipolar attraction, we have observed two triple points inherent in a 2D system: hexatic–liquid–gas and crystal–hexatic–gas, the temperature of the crystal–hexatic–gas triple point is below the hexatic–liquid–gas triple point. This observation may have far-reaching consequences for future studies, since phase diagrams determine possible routes of self-assembly in molecular, protein, and colloidal systems, whereas the attraction range can be adjusted with complex solvents and external electric or magnetic fields. The results obtained may be widely used in condensed matter, chemical physics, materials science, and soft matter.

Viscoelastic response of fluid trapped between two dissimilar van der Waals surfaces

Employing grand canonical Monte-Carlo and molecular dynamics simulations, the viscoelastic response of trapped fluid under molecularly thin confinement by walls having different wall-fluid interaction strengths, is investigated. With increase in slit asymmetry, given by the ratio of interaction strengths of the wall having strong wall-fluid interaction to that of the wall with weak wall-fluid interaction, a crossover in effective density of the fluid film, from rarer (R) to denser (D) than the bulk density is observed. Upon increasing asymmetry further, the dense fluid (F) layers undergo bond-orientational (S) ordering. The variation of viscoelastic relaxation time with scaled asymmetry shows a universal behavior, independent of slit width, with two distinct regimes. Below a critical value of asymmetry, the viscoelastic relaxation time is a slowly varying function of asymmetry, comparable with the structural relaxation time. Beyond the critical asymmetry, on the other hand, viscoelastic response time shows a sharp increase upon increasing asymmetry, deviating markedly from the structural relaxation time. Interestingly the critical asymmetry value is found to correlate with R to D crossover. The microscopic origin of the two-regime universal behavior of viscoelastic response time is found to stem from the fact that below critical asymmetry, the overall viscoelastic behaviour of the slit is dominated by that of the fast relaxing layer close to the weakly attracting surface, while above the critical asymmetry, the relaxation behaviour is guided by the dense fluid layer adjacent to the strongly attracting wall. In vicinity of fluid to ordering transition, the loss and storage moduli merge for low frequencies as in gel-like mechanical behaviour. The storage modulus takes over the loss modulus in the phase co-existence region even before the long ranged order sets in. Our findings bear important implications for fluid transport in hetero-structured geometry in nanotechnology.

Unified analysis of topological defects in 2D systems of active and passive disks

We provide a comprehensive quantitative analysis of localized and extended topological defects in the steady state of 2D passive and active repulsive Brownian disk systems. We show that, both in and out-of-equilibrium, the passage from the solid to the hexatic is driven by the unbinding of dislocations, in quantitative agreement with the KTHNY singularity. Instead, extended clusters of defects largely dominate below the solid-hexatic critical line. The latter percolate in the liquid phase very close to the hexatic-liquid transition, both for continuous and discontinuous transitions, in the homogeneous liquid regime. At critical percolation the clusters of defects are fractal with statistical and geometric properties that are independent of the activity and compatible with the universality class of uncorrelated critical percolation. We also characterize the spatial organization of point-like defects and we show that the disclinations are not free, but rather always very near more complex defect structures. At high activity, the bulk of the dense phase generated by Motility-Induced Phase Separation is characterized by a density of point-like defects, and statistics and morphology of defect clusters, set by the amount of activity and not the packing fraction. Hexatic domains within the dense phase are separated by grain-boundaries along which a finite network of topological defects resides, interrupted by gas bubbles in cavitation. This structure is dynamic in the sense that the defect network allows for an unzipping mechanism that leaves free space for gas bubbles to appear, close, and even be released into the dilute phase.

Signatures of continuous hexatic-liquid transition in two-dimensional melting

Recent studies have shown that the melting of two-dimensional crystals can be either continuous or discontinuous, relying on multiple parameters such as particle stiffness, density, and particle size dispersity. However, what determines the continuity or discontinuity of the two-dimensional melting remains elusive. Here we study the two-dimensional melting of binary mixtures of soft-core particles. The two particle species are different in either particle size or particle stiffness. Starting with the mono-component systems which exhibit discontinuous hexatic-liquid transition, we gradually increase the particle size or stiffness dispersity and find that the hexatic-liquid coexistent region shrinks and eventually vanishes above a critical dispersity. Therefore, the growth of disorder caused by the particle size or stiffness dispersity leads to the discontinuous-continuous transition of the two-dimensional melting. We further find that as long as the melting is continuous the defect concentrations on the boundary between hexatic and liquid phases remain almost constant, accompanied by an almost constant correlation length. These characteristic defect concentrations and correlation length are universal and independent of particle interactions, temperature, and type of particle dispersity, which act as signatures of the continuous two-dimensional melting.