Active matter represents a broad class of systems that evolve far from equilibrium due to the local injection of energy. Like their passive analogs, transformations between distinct metastable states in active matter proceed through rare fluctuations; however, their detailed balance violating dynamics renders these events difficult to study. Here, we present a simulation method for evaluating the rate and mechanism of rare events in generic nonequilibrium systems and apply it to study the conformational changes of a passive solute in an active fluid. The method employs a variational optimization of a control force that renders the rare event a typical one, supplying an exact estimate of its rate as a ratio of path partition functions. Using this method we find that increasing activity in the active bath can enhance the rate of conformational switching of the passive solute in a manner consistent with recent bounds from stochastic thermodynamics.

Very thin elastic sheets, even at zero temperature, exhibit nonlinear elastic response by virtue of their dominant bending modes. Their behavior is even richer at finite temperature. Here, we use molecular dynamics to study the vibrations of a thermally fluctuating two-dimensional elastic sheet with one end clamped at its zero-temperature length. We uncover a tilted phase in which the sheet fluctuates about a mean configuration inclined with respect to the horizontal, thus breaking reflection symmetry. We determine the phase behavior as a function of the aspect ratio of the sheet and the temperature. We show that tilt may be viewed as a type of transverse buckling instability induced by clamping coupled to thermal fluctuations and develop an analytic model that predicts the tilted and untilted regions of the phase diagram. Qualitative agreement is found with the molecular dynamics simulations. Unusual response driven by control of purely geometric quantities like the aspect ratio, as opposed to external fields, offers a very rich playground for two-dimensional mechanical metamaterials.

We study the phase behavior of a quasi-two-dimensional cholesteric liquid crystal shell. We characterize the topological phases arising close to the isotropic-cholesteric transition and show that they differ in a fundamental way from those observed on a flat geometry. For spherical shells, we discover two types of quasi-two-dimensional topological phases: finite quasicrystals and amorphous structures, both made up of mixtures of polygonal tessellations of half-skyrmions. These structures generically emerge instead of regular double twist lattices because of geometric frustration, which disallows a regular hexagonal tiling of curved space. For toroidal shells, the variations in the local curvature of the surface stabilizes heterogeneous phases where cholesteric patterns coexist with hexagonal lattices of half-skyrmions. Quasicrystals and amorphous and heterogeneous structures could be sought experimentally by self-assembling cholesteric shells on the surface of emulsion droplets.

We investigate the hopping dynamics of a colloidal particle across a potential barrier and within a viscoelastic, i.e., non-Markovian, bath and report two clearly separated timescales in the corresponding waiting time distributions. While the longer timescale exponentially depends on the barrier height, the shorter one is similar to the relaxation time of the fluid. This short timescale is a signature of the storage and release of elastic energy inside the bath that strongly increases the hopping rate. Our results are in excellent agreement with numerical simulations of a simple Maxwell model.

We investigate the origin of yield stress aging in semidense, saline, and turbid suspensions in which structural evolution is rapidly arrested by the formation of thermally irreversible roll-resisting interparticle contacts. By performing optical tweezer three-point bending tests on particle rods, we show that these contacts yield by overcoming a rolling threshold, the critical bending moment of which grows logarithmically with time. We demonstrate that this time-dependent contact-scale rolling threshold controls the suspension yield stress and its aging kinetics. We identify a simple constitutive relation between the contact-scale flexural rigidity and rolling threshold, which transfers to macroscopic scales. This leads us to establishing a constitutive relation between macroscopic shear modulus and yield stress that is generic for an array of colloidal systems.

Epithelial monolayers are subjected to various mechanical forces, such as stretching, shearing, and compression. Thus, its mechanical response to external loadings is essential for its biological functions. However, the mechanism of the fracture failure of the epithelial monolayer remains poorly understood. Here, by introducing a new type of topological transition, i.e., detach transition or $T4$ transition, we develop a modified cellular vertex model to investigate the rupture of the cell monolayer. Interestingly, we find a brittle-to-ductile transition in epithelial monolayers, which is controlled by the mechanical properties of single cells and cell-cell contacts. We reveal that the external loadings can activate cell rearrangement in ductile cell monolayers. The plastic deformation results from the nucleation and propagation of “pentagon-heptagon defects” in analogy with the topological defects commonly seen in 2D materials. By using a simplified four-cell model, we further demonstrate that the brittle-to-ductile transition is induced by the competition between cell rearrangement and cell detachment. Our work provides a new theoretical framework to study the rupture of living tissues and may have important implications for many other biological processes, such as wound healing and tissue morphogenesis.

We use molecular simulations to examine how the dynamics of a coarse-grained polymer melt are altered by additives of variable size and interaction strength with the polymer matrix. The effect of diluent size $\sigma$ on polymer dynamics changes significantly when its size is comparable to the polymer segment size. For each $\sigma$, we show that the localization model (LM) quantitatively describes the dependence of the segmental relaxation time $\tau$ on temperature $T$ in terms of dynamic free volume, quantified by the Debye-Waller factor $⟨u^2⟩$. Within this model, we show that the additive size alone controls the functional form of the $T$ dependence. The LM parameters reach asymptotic values when the diluent size exceeds the monomer size, converging to a limit applicable to macroscopic interfaces.

The prediction of flow profiles of slowly sheared granular materials is a major geophysical and industrial challenge. Understanding the role of gravity is particularly important for future planetary exploration in varying gravitational environments. Using the principle of minimization of energy dissipation, and combining experiments and variational analysis, we disentangle the contributions of the gravitational acceleration, confining pressure, and layer thickness on shear strain localization induced by moving fault boundaries at the bottom of a granular layer. The flow profile is independent of the gravity for geometries with a free top surface. However, under a confining pressure or if the sheared layer withstands the weight of the upper layers, increasing gravity promotes the transition from closed shear zones buried in the bulk to open ones that intersect the top surface. We show that the center position and width of the shear zone and the axial angular velocity at the top surface follow universal scaling laws when properly scaled by the gravity, applied pressure, and layer thickness. Our finding that the flow profiles lie on a universal master curve opens the possibility to predict the quasistatic shear flow of granular materials in varying gravitational environments.

Dense active matter is gaining widespread interest due to its remarkable similarity with conventional glass-forming materials. However, active matter is inherently out of equilibrium and even simple models such as active Brownian particles (ABPs) and active Ornstein-Uhlenbeck particles (AOUPs) behave markedly differently from their passive counterparts. Controversially, this difference has been shown to manifest itself via either a speedup, slowdown, or nonmonotonic change of the glassy relaxation dynamics. Here we rationalize these seemingly contrasting views on the departure from equilibrium by identifying the ratio of the short-time length scale to the cage length, i.e., the length scale of local particle caging, as a vital and unifying control parameter for active glassy matter. In particular, we explore the glassy dynamics of both thermal and athermal ABPs and AOUPs upon increasing the persistence time. We find that for all studied systems there is an optimum of the dynamics; this optimum occurs when the cage length coincides with the corresponding short-time length scale of the system, which is either the persistence length for athermal systems or a combination of the persistence length and a diffusive length scale for thermal systems. This new insight, for which we also provide a simple physical argument, allows us to reconcile and explain the manifestly disparate departures from equilibrium reported in many previous studies of dense active materials.

The success of deep learning has revealed the application potential of neural networks across the sciences and opened up fundamental theoretical problems. In particular, the fact that learning algorithms based on simple variants of gradient methods are able to find near-optimal minima of highly nonconvex loss functions is an unexpected feature of neural networks. Moreover, such algorithms are able to fit the data even in the presence of noise, and yet they have excellent predictive capabilities. Several empirical results have shown a reproducible correlation between the so-called flatness of the minima achieved by the algorithms and the generalization performance. At the same time, statistical physics results have shown that in nonconvex networks a multitude of narrow minima may coexist with a much smaller number of wide flat minima, which generalize well. Here, we show that wide flat minima arise as complex extensive structures, from the coalescence of minima around “high-margin” (i.e., locally robust) configurations. Despite being exponentially rare compared to zero-margin ones, high-margin minima tend to concentrate in particular regions. These minima are in turn surrounded by other solutions of smaller and smaller margin, leading to dense regions of solutions over long distances. Our analysis also provides an alternative analytical method for estimating when flat minima appear and when algorithms begin to find solutions, as the number of model parameters varies.

The determination of the coil-globule transition of a polymer is generally based on the reconstruction of scaling laws, implying the need for samples from a rather wide range of different polymer lengths $N$. The spectral point of view developed in this work allows for a very parsimonious description of all the aspects of the finite-size coil-globule transition on the basis of the first two Rouse (cosine) modes only, shedding new light on polymer theory and reintroducing well-established spectral methods that have been surprisingly neglected in this field so far. Capturing the relevant configuration path features, the proposed approach enables one to determine the state of a polymer without the need of any information about the polymer length or interaction strength. Importantly, we propose an experimental implementation of our analysis that can be easily performed with modern fluorescent imaging techniques and would allow differentiation of coil or globule conformations by simply recording the positions of at least three discernible loci on the polymer.

Using computed x-ray tomography we determine the three dimensional (3D) structure of binary hard sphere mixtures as a function of composition and size ratio of the particles $q$. Using a recently introduced four-point correlation function we reveal that this 3D structure has on intermediate and large length scales a surprisingly regular order, the symmetry of which depends on $q$. The related structural correlation length has a minimum at the composition at which the packing fraction is highest. At this composition also the number of different local particle arrangements has a maximum, indicating that efficient packing of particles is associated with a structure that is locally maximally disordered.

We present the hydrodynamic theory of coherent collective motion (“flocking”) at a solid-liquid interface, and many of its predictions for experiment. We find that such systems are stable, and have long-range orientational order, over a wide range of parameters. When stable, these systems exhibit “giant number fluctuations,” which grow as the $3/4$th power of the mean number. Stable systems also exhibit anomalous rapid diffusion of tagged particles suspended in the passive fluid along any directions in a plane parallel to the solid-liquid interface, whereas the diffusivity along the direction perpendicular to the plane is not anomalous. In the remaining parameter space, the system becomes unstable.

Suspensions of active agents with nematic interactions exhibit complex spatiotemporal dynamics such as mesoscale turbulence. Since the Reynolds number of microscopic flows is very small on the scale of individual agents, inertial effects are typically excluded in continuum theories of active nematic turbulence. Whether active stresses can collectively excite inertial flows is currently unclear. To address this question, we investigate a two-dimensional continuum theory for active nematic turbulence. In particular, we compare mesoscale turbulence with and without the effects of advective inertia. We find that inertial effects can influence the flow already close to the onset of the turbulent state and, moreover, give rise to large-scale fluid motion for strong active driving. A detailed analysis of the kinetic energy budget reveals an energy transfer to large scales mediated by inertial advection. While this transfer is small in comparison to energy injection and dissipation, its effects accumulate over time. The inclusion of friction, which is typically present in experiments, can compensate for this effect. The findings suggest that the inclusion of inertia and friction may be necessary for dynamically consistent theories of active nematic turbulence.

Despite there being an infinite variety of types of flow, most rheological studies focus on a single type such as simple shear. Using discrete element simulations, we explore bulk granular systems in a wide range of flow types at large strains and characterize invariants of the stress tensor for different inertial numbers and interparticle friction coefficients. We identify a strong dependence on the type of flow, which grows with increasing inertial number or friction. Standard models of yielding, repurposed to describe the dependence of the stress on flow type in steady-state flow and at finite rates, are compared with data.