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.

Remembering the Work of Phillip L. Geissler: A Coda to His Scientific Trajectory.

Phillip L. Geissler made important contributions to the statistical mechanics of biological polymers, heterogeneous materials, and chemical dynamics in aqueous environments. He devised analytical and computational methods that revealed the underlying organization of complex systems at the frontiers of biology, chemistry, and materials science. In this retrospective we celebrate his work at these frontiers.

Computing Free Energies with Fluctuation Relations on Quantum Computers.

As a central thermodynamic property, free energy enables the calculation of virtually any equilibrium property of a physical system, allowing for the construction of phase diagrams and predictions about transport, chemical reactions, and biological processes. Thus, methods for efficiently computing free energies, which in general is a difficult problem, are of great interest to broad areas of physics and the natural sciences. The majority of techniques for computing free energies target classical systems, leaving the computation of free energies in quantum systems less explored. Recently developed fluctuation relations enable the computation of free energy differences in quantum systems from an ensemble of dynamic simulations. While performing such simulations is exponentially hard on classical computers, quantum computers can efficiently simulate the dynamics of quantum systems. Here, we present an algorithm utilizing a fluctuation relation known as the Jarzynski equality to approximate free energy differences of quantum systems on a quantum computer. We discuss under which conditions our approximation becomes exact, and under which conditions it serves as a strict upper bound. Furthermore, we successfully demonstrate a proof of concept of our algorithm using the transverse field Ising model on a real quantum processor. As quantum hardware continues to improve, we anticipate that our algorithm will enable computation of free energy differences for a wide range of quantum systems useful across the natural sciences.

Phase Diagram of Active Brownian Spheres: Crystallization and the Metastability of Motility-Induced Phase Separation.

Motility-induced phase separation (MIPS), the phenomenon in which purely repulsive active particles undergo a liquid-gas phase separation, is among the simplest and most widely studied examples of a nonequilibrium phase transition. Here, we show that states of MIPS coexistence are in fact only metastable for three-dimensional active Brownian particles over a very broad range of conditions, decaying at long times through an ordering transition we call active crystallization. At an activity just above the MIPS critical point, the liquid-gas binodal is superseded by the crystal-fluid coexistence curve, with solid, liquid, and gas all coexisting at the triple point where the two curves intersect. Nucleating an active crystal from a disordered fluid, however, requires a rare fluctuation that exhibits the nearly close-packed density of the solid phase. The corresponding barrier to crystallization is surmountable on a feasible timescale only at high activity, and only at fluid densities near maximal packing. The glassiness expected for such dense liquids at equilibrium is strongly mitigated by active forces, so that the lifetime of liquid-gas coexistence declines steadily with increasing activity, manifesting in simulations as a facile spontaneous crystallization at extremely high activity.

Entropy production fluctuations encode collective behavior in active matter.

We derive a general lower bound on distributions of entropy production in interacting active matter systems. The bound is tight in the limit that interparticle correlations are small and short-ranged, which we explore in four canonical active matter models. In all models studied, the bound is weak where collective fluctuations result in long-ranged correlations, which subsequently links the locations of phase transitions to enhanced entropy production fluctuations. We develop a theory for the onset of enhanced fluctuations and relate it to specific phase transitions in active Brownian particles. We also derive optimal control forces that realize the dynamics necessary to tune dissipation and manipulate the system between phases. In so doing, we uncover a general relationship between entropy production and pattern formation in active matter, as well as ways of controlling it.

Time reversal symmetry breaking and odd viscosity in active fluids: Green-Kubo and NEMD results.

Active fluids, which are driven at the microscale by non-conservative forces, are known to exhibit novel transport phenomena due to the breaking of time reversal symmetry. Recently, Epstein and Mandadapu [arXiv:1907.10041 (2019)] obtained Green-Kubo relations for the full set of viscous coefficients governing isotropic chiral active fluids, including the so-called odd viscosity, invoking Onsager's regression hypothesis for the decay of fluctuations in active non-equilibrium steady states. In this Communication, we test these Green-Kubo relations using molecular dynamics simulations of a canonical model system consisting of actively torqued dumbbells. We find the resulting odd and shear viscosity values from the Green-Kubo relations to be in good agreement with values measured independently through non-equilibrium molecular dynamics flow simulations. This provides a test of the Green-Kubo relations and lends support to the application of the Onsager regression hypothesis in relation to viscous behaviors of active matter systems.

Statistical mechanics of transport processes in active fluids. II. Equations of hydrodynamics for active Brownian particles.

We perform a coarse-graining analysis of the paradigmatic active matter model, active Brownian particles, yielding a continuum description in terms of balance laws for mass, linear and angular momentum, and energy. The derivation of the balance of linear momentum reveals that the active force manifests itself directly as a continuum-level body force proportional to an order parameter-like director field, which therefore requires its own evolution equation to complete the continuum description of the system. We derive this equation, demonstrating in the process that bulk currents may be sustained in homogeneous systems only in the presence of interparticle aligning interactions. Furthermore, we perform a second coarse-graining of the balance of linear momentum and derive the expression for active or swim pressure in the case of mechanical equilibrium.

Rare behavior of growth processes via umbrella sampling of trajectories.

We compute probability distributions of trajectory observables for reversible and irreversible growth processes. These results reveal a correspondence between reversible and irreversible processes, at particular points in parameter space, in terms of their typical and atypical trajectories. Thus key features of growth processes can be insensitive to the precise form of the rate constants used to generate them, recalling the insensitivity to microscopic details of certain equilibrium behavior. We obtained these results using a sampling method, inspired by the "s-ensemble" large-deviation formalism, that amounts to umbrella sampling in trajectory space. The method is a simple variant of existing approaches, and applies to ensembles of trajectories controlled by the total number of events. It can be used to determine large-deviation rate functions for trajectory observables in or out of equilibrium.

Phase separation and large deviations of lattice active matter.

Off-lattice active Brownian particles form clusters and undergo phase separation even in the absence of attractions or velocity-alignment mechanisms. Arguments that explain this phenomenon appeal only to the ability of particles to move persistently in a direction that fluctuates, but existing lattice models of hard particles that account for this behavior do not exhibit phase separation. Here we present a lattice model of active matter that exhibits motility-induced phase separation in the absence of velocity alignment. Using direct and rare-event sampling of dynamical trajectories, we show that clustering and phase separation are accompanied by pronounced fluctuations of static and dynamic order parameters. This model provides a complement to off-lattice models for the study of motility-induced phase separation.

Statistical mechanics of transport processes in active fluids: Equations of hydrodynamics.

The equations of hydrodynamics including mass, linear momentum, angular momentum, and energy are derived by coarse-graining the microscopic equations of motion for systems consisting of rotary dumbbells driven by internal torques. In deriving the balance of linear momentum, we find that the symmetry of the stress tensor is broken due to the presence of non-zero torques on individual particles. The broken symmetry of the stress tensor induces internal spin in the fluid and leads us to consider the balance of internal angular momentum in addition to the usual moment of momentum. In the absence of spin, the moment of momentum is the same as the total angular momentum. In deriving the form of the balance of total angular momentum, we find the microscopic expressions for the couple stress tensor that drives the spin field. We show that the couple stress contains contributions from both intermolecular interactions and the active forces. The presence of spin leads to the idea of balance of moment of inertia due to the constant exchange of particles in a small neighborhood around a macroscopic point. We derive the associated balance of moment of inertia at the macroscale and identify the moment of inertia flux that induces its transport. Finally, we obtain the balances of total and internal energy of the active fluid and identify the sources of heat and heat fluxes in the system.

Entropy Production and Fluctuation Theorems for Active Matter.

Active biological systems reside far from equilibrium, dissipating heat even in their steady state, thus requiring an extension of conventional equilibrium thermodynamics and statistical mechanics. In this Letter, we have extended the emerging framework of stochastic thermodynamics to active matter. In particular, for the active Ornstein-Uhlenbeck model, we have provided consistent definitions of thermodynamic quantities such as work, energy, heat, entropy, and entropy production at the level of single, stochastic trajectories and derived related fluctuation relations. We have developed a generalization of the Clausius inequality, which is valid even in the presence of the non-Hamiltonian dynamics underlying active matter systems. We have illustrated our results with explicit numerical studies.

Similarity of ensembles of trajectories of reversible and irreversible growth processes.

Models of bacterial growth tend to be "irreversible," allowing for the number of bacteria in a colony to increase but not to decrease. By contrast, models of molecular self-assembly are usually "reversible," allowing for the addition and removal of particles to a structure. Such processes differ in a fundamental way because only reversible processes possess an equilibrium. Here we show at the mean-field level that dynamic trajectories of reversible and irreversible growth processes are similar in that both feel the influence of attractors, at which growth proceeds without limit but the intensive properties of the system are invariant. Attractors of both processes undergo nonequilibrium phase transitions as model parameters are varied, suggesting a unified way of describing typical properties of reversible and irreversible growth. We also establish a connection at the mean-field level between an irreversible model of growth (the magnetic Eden model) and the equilibrium Ising model, supporting the findings made by other authors using numerical simulations.

Microscopic origin and macroscopic implications of lane formation in mixtures of oppositely driven particles.

Colloidal particles of two types, driven in opposite directions, can segregate into lanes [Vissers et al., Soft Matter 7, 2352 (2011)1744-683X10.1039/c0sm01343a]. This phenomenon can be reproduced by two-dimensional Brownian dynamics simulations of model particles [Dzubiella et al., Phys. Rev. E 65, 021402 (2002)1063-651X10.1103/PhysRevE.65.021402]. Here we use computer simulation to assess the generality of lane formation with respect to variation of particle type and dynamical protocol. We find that laning results from rectification of diffusion on the scale of a particle diameter: oppositely driven particles must, in the time taken to encounter each other in the direction of the drive, diffuse in the perpendicular direction by about one particle diameter. This geometric constraint implies that the diffusion constant of a particle, in the presence of those of the opposite type, grows approximately linearly with the Péclet number, a prediction confirmed by our numerics over a range of model parameters. Such environment-dependent diffusion is statistically similar to an effective interparticle attraction; consistent with this observation, we find that oppositely driven nonattractive colloids display features characteristic of the simplest model system possessing both interparticle attractions and persistent motion, the driven Ising lattice gas [Katz, Leibowitz, and Spohn, J. Stat. Phys. 34, 497 (1984)JSTPBS0022-471510.1007/BF01018556]. These features include long-ranged correlations in the disordered regime, a critical regime characterized by a change in slope of the particle current with the Péclet number, and fluctuations that grow with system size. By analogy, we suggest that lane formation in the driven colloid system is a phase transition in the macroscopic limit, but that macroscopic phase separation would not occur in finite time upon starting from disordered initial conditions.

Glutamate receptor desensitization is mediated by changes in quaternary structure of the ligand binding domain.

Glutamate receptor ion channels are membrane proteins that mediate excitatory synaptic transmission in the central nervous system of vertebrates. Insight into molecular mechanisms underlying glutamate receptor gating is limited by lack of structural information for receptors trapped in different conformational states. Here, we report the use of single-particle cryoelectron tomography to determine the structures, at ∼21 Å resolution, of full-length GluK2 kainate receptors trapped in antagonist-bound resting and agonist-bound desensitized states. The resting state, stabilized by the competitive antagonist LY466195, closely resembles the crystal structure of the AMPA receptor GluA2, with well-resolved proximal and distal subunits exhibiting cross-over between the twofold symmetric amino terminal domain and a twofold symmetric ligand binding domain (LBD) dimer of dimers assembly. In the desensitized state, the LBD undergoes a major rearrangement, resulting in a separation of the four subunits by ∼25 Å. However, the amino terminal domain, transmembrane, and cytoplasmic regions of the receptor have similar conformations in the resting and desensitized states. The LBD rearrangement was not anticipated in prior models based on crystal structures for soluble LBD dimer assemblies, and we speculate that subunit separation allows a better match to the fourfold symmetric ion channel domain. From fits of the amino terminal domain and LBD domains into the density map of the desensitized state we have derived a structural model for differences in quaternary conformation between the resting and desensitized states.

Promoters recognized by forkhead proteins exist for individual 21U-RNAs.

C. elegans 21U-RNAs are equivalent to the piRNAs discovered in other metazoans and have important roles in gametogenesis and transposon control. The biogenesis and molecular function of 21U-RNAs and piRNAs are poorly understood. Here, we demonstrate that transcription of each 21U-RNA is regulated separately through a conserved upstream DNA motif. We use genomic analysis to show that this motif is associated with low nucleosome occupancy, a characteristic of many promoters that drive expression of protein-coding genes, and that RNA polymerase II is localized to this nucleosome-depleted region. We establish that the most conserved 8-mer sequence in the upstream region of 21U-RNAs, CTGTTTCA, is absolutely required for their individual expression. Furthermore, we demonstrate that the 8-mer is specifically recognized by Forkhead family (FKH) transcription factors and that 21U-RNA expression is diminished in several FKH mutants. Our results suggest that thousands of small noncoding transcription units are regulated by FKH proteins.

Free energy of multiple overlapping chains.

How accurate is pair additivity in describing interactions between soft polymer-based nanoparticles? Using numerical simulations we compute the free energy cost required to overlap multiple chains in the same region of space, and provide a quantitative measure of the effectiveness of pair additivity as a function of chain number and length. Our data suggest that pair additivity can indeed become quite inadequate as the chain density in the overlapping region increases. We also show that even a scaling theory based on polymer confinement can only partially account for the complexity of the problem. In fact, we unveil and characterize an isotropic to star-polymer crossover taking place for large number of chains, and propose a revised scaling theory that better captures the physics of the problem.

Free energy of alternating two-component polymer brushes on cylindrical templates.

We use computer simulations to investigate the stability of a two-component polymer brush de-mixing on a curved template into phases of different morphological properties. It has been previously shown via molecular dynamics simulations that immiscible chains having different length and anchored to a cylindrical template will phase separate into stripes of different widths oriented perpendicularly to the cylindrical axis. We calculate free energy differences for a variety of stripe widths, and extract simple relationships between the sizes of the two polymers, N(1) and N(2), and the free energy dependence on the stripe width. We explain these relationships using simple physical arguments based upon previous theoretical work on the free energy of polymer brushes.