An official website of the United States government
We introduce a method for solving the “inverse” phase equilibria problem: How should the interactions among a collection of molecular species be designed in order to achieve a target phase diagram? Using techniques from convex optimization theory, we show how to solve this problem for phase diagrams containing a large number of components and many coexisting phases with prescribed compositions. We apply our approach to commonly used mean-field models of multicomponent fluids and then use molecular simulations to verify that the designed interactions result in the target phase diagrams. Our approach enables the rational design of “programmable” fluids, such as biopolymer and colloidal mixtures, with complex phase behavior.
more »
« less
Phase separation in fluids with many interacting components
Fluids in natural systems, like the cytoplasm of a cell, often contain thousands of molecular species that are organized into multiple coexisting phases that enable diverse and specific functions. How interactions between numerous molecular species encode for various emergent phases is not well understood. Here, we leverage approaches from random-matrix theory and statistical physics to describe the emergent phase behavior of fluid mixtures with many species whose interactions are drawn randomly from an underlying distribution. Through numerical simulation and stability analyses, we show that these mixtures exhibit staged phase-separation kinetics and are characterized by multiple coexisting phases at steady state with distinct compositions. Random-matrix theory predicts the number of coexisting phases, validated by simulations with diverse component numbers and interaction parameters. Surprisingly, this model predicts an upper bound on the number of phases, derived from dynamical considerations, that is much lower than the limit from the Gibbs phase rule, which is obtained from equilibrium thermodynamic constraints. We design ensembles that encode either linear or nonmonotonic scaling relationships between the number of components and coexisting phases, which we validate through simulation and theory. Finally, inspired by parallels in biological systems, we show that including nonequilibrium turnover of components through chemical reactions can tunably modulate the number of coexisting phases at steady state without changing overall fluid composition. Together, our study provides a model framework that describes the emergent dynamical and steady-state phase behavior of liquid-like mixtures with many interacting constituents.
more »
« less
- Award ID(s):
- 1764269
- PAR ID:
- 10328865
- Date Published:
- Journal Name:
- Proceedings of the National Academy of Sciences
- Volume:
- 118
- Issue:
- 45
- ISSN:
- 0027-8424
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
The separation of substances into different phases is ubiquitous in nature and important scientifically and technologically. This phenomenon may become drastically different if the species involved, whether molecules or supramolecular assemblies, interconvert. In the presence of an external force large enough to overcome energetic differences between the interconvertible species (forced interconversion), the two alternative species will be present in equal amounts, and the striking phenomenon of steady-state, restricted phase separation into mesoscales is observed. Such microphase separation is one of the simplest examples of dissipative structures in condensed matter. In this work, we investigate the formation of such mesoscale steady-state structures through Monte Carlo and molecular dynamics simulations of three physically distinct microscopic models of binary mixtures that exhibit both equilibrium (natural) interconversion and a nonequilibrium source of forced interconversion. We show that this source can be introduced through an internal imbalance of intermolecular forces or an external flux of energy that promotes molecular interconversion, possible manifestations of which could include the internal nonequilibrium environment of living cells or a flux of photons. The main trends and observations from the simulations are well captured by a nonequilibrium thermodynamic theory of phase transitions affected by interconversion. We show how a nonequilibrium bicontinuous microemulsion or a spatially modulated state may be generated depending on the interplay between diffusion, natural interconversion, and forced interconversion.more » « less
-
Abstract Placed slightly out of dynamical equilibrium, an isolated stellar system quickly returns toward a steady virialized state. We study this process of collisionless relaxation using the matrix method of linear response theory. We show that the full phase-space distribution of the final virialized state can be recovered directly from the disequilibrium initial conditions, without the need to compute the time evolution of the system. This shortcut allows us to determine the final virialized configuration with minimal computational effort. Complementing this result, we develop tools to model the system's full time evolution in the linear approximation. In particular, we show that moments of the velocity distribution can be efficiently computed using a generalized moment matrix. We apply our linear methods to study the relaxation of energy-truncated Hernquist spheres, mimicking the tidal stripping of a cuspy dark matter subhalo. Comparison of our linear predictions against controlled, isolatedN-body simulations shows agreement at percent level for the parts of the system where a linear response to the perturbation is expected. We find that relaxation generates a tangential velocity anisotropy in the intermediate regions, despite the initial disequilibrium state having isotropic kinematics. Our results also strengthen the case for relaxation depleting the amplitude of the density cusp, without affecting its asymptotic slope. Finally, we compare the linear theory against anN-body simulation of tidal stripping on a radial orbit, confirming that the theory still accurately predicts density and velocity dispersion profiles for most of the system.more » « less
-
The moving contact line between a fluid, liquid and solid is a ubiquitous phenomenon, and determining the maximum speed at which a liquid can wet/dewet a solid is a practically important problem. Using continuum models, previous studies have shown that the maximum speed of wetting/dewetting can be found by calculating steady solutions of the governing equations and locating the critical capillary number, $$Ca_{{crit}}$$ , above which no steady-state solution can be found. Below $$Ca_{{crit}}$$ , both stable and unstable steady-state solutions exist and if some appropriate measure of these solutions is plotted against $Ca$ , a fold bifurcation appears where the stable and unstable branches meet. Interestingly, the significance of this bifurcation structure to the transient dynamics has yet to be explored. This article develops a computational model and uses ideas from dynamical systems theory to show the profound importance of the unstable solutions on the transient behaviour. By perturbing the stable state by the eigenmodes calculated from a linear stability analysis it is shown that the unstable branch is an ‘edge’ state that is responsible for the eventual dynamical outcomes and that the system can become transient when $$Ca< Ca_{{crit}}$$ due to finite-amplitude perturbations. Furthermore, when $$Ca>Ca_{{crit}}$$ , we show that the trajectories in phase space closely follow the unstable branch.more » « less
-
Abstract Metabolic cross-feeding plays vital roles in promoting ecological diversity. While some microbes depend on exchanges of essential nutrients for growth, the forces driving the extensive cross-feeding needed to support the coexistence of free-living microbes are poorly understood. Here we characterize bacterial physiology under self-acidification and establish that extensive excretion of key metabolites following growth arrest provides a collaborative, inter-species mechanism of stress resistance. This collaboration occurs not only between species isolated from the same community, but also between unrelated species with complementary (glycolytic vs. gluconeogenic) modes of metabolism. Cultures of such communities progress through distinct phases of growth-dilution cycles, comprising of exponential growth, acidification-triggered growth arrest, collaborative deacidification, and growth recovery, with each phase involving different combinations of physiological states of individual species. Our findings challenge the steady-state view of ecosystems commonly portrayed in ecological models, offering an alternative dynamical view based on growth advantages of complementary species in different phases.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1073/pnas.2108551118
