We outline how principal component analysis can be applied to particle configuration data to detect a variety of phase transitions in off-lattice systems, both in and out of equilibrium. Specifically, we discuss its application to study (1) the nonequilibrium random organization (RandOrg) model that exhibits a phase transition from quiescent to steady-state behavior as a function of density, (2) orientationally and positionally driven equilibrium phase transitions for hard ellipses, and (3) a compositionally driven demixing transition in the non-additive binary Widom-Rowlinson mixture.
Tuning nonequilibrium phase transitions with inertia
In striking contrast to equilibrium systems, inertia can profoundly alter the structure of active systems. Here, we demonstrate that driven systems can exhibit effective equilibrium-like states with increasing particle inertia, despite rigorously violating the fluctuation–dissipation theorem. Increasing inertia progressively eliminates motility-induced phase separation and restores equilibrium crystallization for active Brownian spheres. This effect appears to be general for a wide class of active systems, including those driven by deterministic time-dependent external fields, whose nonequilibrium patterns ultimately disappear with increasing inertia. The path to this effective equilibrium limit can be complex, with finite inertia sometimes acting to accentuate nonequilibrium transitions. The restoration of near equilibrium statistics can be understood through the conversion of active momentum sources to passive-like stresses. Unlike truly equilibrium systems, the effective temperature is now density dependent, the only remnant of the nonequilibrium dynamics. This density-dependent temperature can in principle introduce departures from equilibrium expectations, particularly in response to strong gradients. Our results provide additional insight into the effective temperature ansatz while revealing a mechanism to tune nonequilibrium phase transitions.
more » « less- Award ID(s):
- 1734030
- NSF-PAR ID:
- 10440220
- Publisher / Repository:
- American Institute of Physics
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 158
- Issue:
- 7
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
-
Living systems are composed of a select number of biopolymers and minerals yet exhibit an immense diversity in materials properties. The wide-ranging characteristics, such as enhanced mechanical properties of skin and bone, or responsive optical properties derived from structural coloration, are a result of the multiscale, hierarchical structure of the materials. The fields of materials and polymer chemistry have leveraged equilibrium concepts in an effort to mimic the structure complex materials seen in nature. However, realizing the remarkable properties in natural systems requires moving beyond an equilibrium perspective. An alternative method to create materials with multiscale structures is to approach the issue from a kinetic perspective and utilize chemical processes to drive phase transitions. This Account features an active area of research in our group, reaction-induced phase transitions (RIPT), which uses chemical reactions such as polymerizations to induce structural changes in soft material systems. Depending on the type of phase transition (e.g., microphase versus macrophase separation), the resulting change in state will occur at different length scales (e.g., nm – μm), thus dictating the structure of the material. For example, the in situ formation of either a block copolymer or a homopolymer initially in a monomer mixture during a polymerization will drive nanoscale or macroscale transitions, respectively. Specifically, three different examples utilizing reaction-driven phase changes will be discussed: 1) in situ polymer grafting from block copolymers, 2) multiscale polymer nanocomposites, and 3) Lewis adduct-driven phase transitions. All three areas highlight how chemical changes via polymerizations or specific chemical binding result in phase transitions that lead to nanostructural and multiscale changes. Harnessing kinetic chemical processes to promote and control material structure, as opposed to organizing pre-synthesized molecules, polymers, or nanoparticles within a thermodynamic framework, is a growing area of interest. Trapping nonequilibrium states in polymer materials has been primarily focused from a polymer chain conformation viewpoint in which synthesized polymers are subjected to different thermal and processing conditions. The impact of reaction kinetics and polymerization rate on final polymer material structure is starting to be recognized as a new way to access different morphologies not available through thermodynamic means. Furthermore, kinetic control of polymer material structure is not specific to polymerizations and encompasses any chemical reaction that induce morphology transitions. Kinetically driven processes to dictate material structure directly impact a broad range of areas including separation membranes, biomolecular condensates, cell mobility, and the self-assembly of polymers and colloids. Advancing polymer material syntheses using kinetic principles such as RIPT opens new possibilities for dictating material structure and properties beyond what is currently available with traditional self-assembly techniques.more » « less
-
more » « less
An outstanding problem in statistical mechanics is the determination of whether prescribed functional forms of the pair correlation function g2(r) [or equivalently, structure factor S(k)] at some number density ρ can be achieved by many-body systems in d-dimensional Euclidean space. The Zhang–Torquato conjecture states that any realizable set of pair statistics, whether from a nonequilibrium or equilibrium system, can be achieved by equilibrium systems involving up to two-body interactions. To further test this conjecture, we study the realizability problem of the nonequilibrium iso-g2 process, i.e., the determination of density-dependent effective potentials that yield equilibrium states in which g2 remains invariant for a positive range of densities. Using a precise inverse algorithm that determines effective potentials that match hypothesized functional forms of g2(r) for all r and S(k) for all k, we show that the unit-step function g2, which is the zero-density limit of the hard-sphere potential, is remarkably realizable up to the packing fraction ϕ = 0.49 for d = 1. For d = 2 and 3, it is realizable up to the maximum “terminal” packing fraction ϕc = 1/2d, at which the systems are hyperuniform, implying that the explicitly known necessary conditions for realizability are sufficient up through ϕc. For ϕ near but below ϕc, the large-r behaviors of the effective potentials are given exactly by the functional forms exp[ − κ(ϕ)r] for d = 1, r−1/2 exp[ − κ(ϕ)r] for d = 2, and r−1 exp[ − κ(ϕ)r] (Yukawa form) for d = 3, where κ−1(ϕ) is a screening length, and for ϕ = ϕc, the potentials at large r are given by the pure Coulomb forms in the respective dimensions as predicted by Torquato and Stillinger [Phys. Rev. E 68, 041113 (2003)]. We also find that the effective potential for the pair statistics of the 3D “ghost” random sequential addition at the maximum packing fraction ϕc = 1/8 is much shorter ranged than that for the 3D unit-step function g2 at ϕc; thus, it does not constrain the realizability of the unit-step function g2. Our inverse methodology yields effective potentials for realizable targets, and, as expected, it does not reach convergence for a target that is known to be non-realizable, despite the fact that it satisfies all known explicit necessary conditions. Our findings demonstrate that exploring the iso-g2 process via our inverse methodology is an effective and robust means to tackle the realizability problem and is expected to facilitate the design of novel nanoparticle systems with density-dependent effective potentials, including exotic hyperuniform states of matter.
-
more » « less
Abstract This study investigates the emergence of hurricane‐like vortices in idealized simulations of rotating moist convection. A Boussinesq atmosphere with simplified thermodynamics for phase transitions is forced by prescribing the temperature and humidity at the upper and lower boundaries. The governing equations are solved numerically using a variable‐density incompressible Navier‐Stokes solver with adaptive mesh refinement to explore the behavior of moist convection under a broad range of conditions. In the absence of rotation, convection aggregates into active patches separated by large unsaturated regions. Rotation modulates this statistical equilibrium state so that the self‐aggregated convection organizes hurricane‐like vortices. The warm and saturated air converges to the center of the vortices, and the latent heat released through the upwelling, forms the warm core structure. These hurricane‐like vortices share characteristics similar to tropical cyclones in the earth's atmosphere. The hurricane‐like vortices occur under conditionally unstable conditions where the potential energy given at the boundaries is large enough, corresponding to a moderate rate of rotation. This regime shares many similar characteristics to the tropical atmosphere indicating that the formation of intense meso‐scale vortices is a general characteristic of rotating moist convection. The model used here does not include any interactions with radiation, wind‐evaporation feedback, or cloud microphysics, indicating that, while these processes may be relevant for tropical cyclogenesis in the Earth atmosphere, they are not its primary cause. Instead, our results confirm that the formation and maintenance of hurricane‐like vortices involve a combination of atmospheric dynamics under the presence of rotation and of phase transitions.
-
more » « less
Chemically driven fluids can demix to form condensed droplets that exhibit phase behaviors not observed at equilibrium. In particular, nonequilibrium interfacial properties can emerge when the chemical reactions are driven differentially between the interior and exterior of the phase-separated droplets. Here, we use a minimal model to study changes in the interfacial tension between coexisting phases away from equilibrium. Simulations of both droplet nucleation and interface roughness indicate that the nonequilibrium interfacial tension can either be increased or decreased relative to its equilibrium value, depending on whether the driven chemical reactions are accelerated or decelerated within the droplets. Finally, we show that these observations can be understood using a predictive theory based on an effective thermodynamic equilibrium.
