We combine percolation theory and Monte Carlo simulation to study in two dimensions the connectivity of an equilibrium lattice model of interacting Janus disks which self-assemble into an orientationally ordered stripe phase at low temperature. As the patch size is increased or the temperature is lowered, clusters of patch-connected disks grow, and a percolating cluster emerges at a threshold. In the stripe phase, the critical clusters extend longer in the direction parallel to the stripes than in the perpendicular direction, and percolation is thus anisotropic. It is found that the critical behavior of percolation in the Janus system is consistent with that of standard isotropic percolation, when an appropriate spatial rescaling is made. The rescaling procedure can be applied to understand other anisotropic systems, such as the percolation of aligned rigid rods and of the $q$-state Potts model with anisotropic interactions.

We study the energy cost of flocking in the active Ising model (AIM) and show that, besides the energy cost for self-propelled motion, an additional energy dissipation is required to power the alignment of spins. We find that this additional alignment dissipation reaches its maximum at the flocking transition point in the form of a cusp with a discontinuous first derivative with respect to the control parameter. To understand this singular behavior, we analytically solve the two- and three-site AIM models and obtain the exact dependence of the alignment dissipation on the flocking order parameter and control parameter, which explains the cusped dissipation maximum at the flocking transition. Our results reveal a trade-off between the energy cost of the system and its performance measured by the flocking speed and sensitivity to external perturbations. This trade-off relationship provides a new perspective for understanding the dynamics of natural flocks and designing optimal artificial flocking systems.

We study numerically the Toner-Tu field theory where the density field is maintained constant, a limit case of “Malthusian” flocks for which the asymptotic scaling of correlation functions in the ordered phase is known exactly. While we confirm these scaling laws, we also show that such constant-density flocks are metastable to the nucleation of a specific defect configuration, and are replaced by a globally disordered phase consisting of asters surrounded by shock lines that constantly evolves and remodels itself. We demonstrate that the main source of disorder lies along shock lines, rendering this active foam fundamentally different from the corresponding equilibrium system. We thus show that in the context of active matter also, a result obtained at all orders of perturbation theory can be superseded by nonperturbative effects, calling for a different approach.

The sequential exchange of filament composition to increase filament curvature was proposed as a mechanism for how some biological polymers deform and cut membranes. The relationship between the filament composition and its mechanical effect is lacking. We develop a kinetic model for the assembly of composite filaments that includes protein–membrane adhesion, filament mechanics and membrane mechanics. We identify the physical conditions for such a membrane remodeling and show this mechanism of sequential polymer assembly lowers the energetic barrier for membrane deformation.

We use a continuum model to examine the effect of activity on a phase-separating mixture of an extensile active nematic and a passive fluid. We highlight the distinct role of (i) previously considered interfacial active stresses and (ii) bulk active stresses that couple to liquid crystalline degrees of freedom. Interfacial active stresses can arrest phase separation, as previously demonstrated. Bulk extensile active stresses can additionally strongly suppress phase separation by sustained self-stirring of the fluid, substantially reducing the size of the coexistence region in the temperature-concentration plane relative to that of the passive system. The phase-separated state is a dynamical emulsion of continuously splitting and merging droplets, as suggested by recent experiments. Using scaling analysis and simulations, we identify various regimes for the dependence of droplet size on activity. These results can provide a criterion for identifying the mechanisms responsible for arresting phase separation in experiments.

We study the critical features of the order parameter’s fluctuations near the threshold of mixed-order phase transitions in randomly interdependent spatial networks. Remarkably, we find that although the structure of the order parameter is not scale invariant, its fluctuations are fractal up to a well-defined correlation length ${\xi }^{\prime }$ that diverges when approaching the mixed-order transition threshold. We characterize the self-similar nature of these critical fluctuations through their effective fractal dimension $d_f^{\prime }=3d/4$, and correlation length exponent ${\nu }^{\prime }=2/d$, where $d$ is the dimension of the system. By analyzing percolation and magnetization, we demonstrate that $d_f^{\prime }$ and ${\nu }^{\prime }$ are the same for both, i.e., independent of the symmetry of the process for any $d$ of the underlying networks.

Active nematics can be modeled using phenomenological continuum theories that account for the dynamics of the nematic director and fluid velocity through partial differential equations (PDEs). While these models provide a statistical description of the experiments, the relevant terms in the PDEs and their parameters are usually identified indirectly. We adapt a recently developed method to automatically identify optimal continuum models for active nematics directly from spatiotemporal data, via sparse regression of the coarse-grained fields onto generic low order PDEs. After extensive benchmarking, we apply the method to experiments with microtubule-based active nematics, finding a surprisingly minimal description of the system. Our approach can be generalized to gain insights into active gels, microswimmers, and diverse other experimental active matter systems.

Hydrodynamic interactions (HIs) play a critical role in the self-organization of colloidal suspensions and biological solutions. However, their roles have remained elusive particularly for charged soft matter systems. Here we consider the role of HIs in the self-assembly of oppositely charged colloidal particles, which is a promising candidate for electrical tunable soft materials. We employ the fluid particle dynamics method to consider many-body HIs and the coupling between the colloid, ion, and fluid motions. We find that, under a constant electric field, oppositely charged colloidal particles form clusters and percolate into a gel network, unlike bundlelike aggregates aligned in the field direction observed by Brownian dynamics simulations neglecting HIs. We reveal that the cluster-forming tendency originates from the incompressibility-induced “inverse squeezing flow” effect that dramatically slows down the disaggregation of attached colloids. Our findings indicate that the HI selects a unique kinetic pathway to the nonequilibrium colloidal self-assembly.

We demonstrate nonequilibrium scaling laws for the aging and equilibration dynamics in glass formers that emerge from combining a relaxation equation for the static structure with the equilibrium scaling laws of glassy dynamics. Different scaling regimes are predicted for the evolution of the structural relaxation time $\tau$ with age (waiting time $t_w$), depending on the depth of the quench from the liquid into the glass: “simple” aging ($\tau \sim t_w$) applies for quenches close to the critical point of mode-coupling theory (MCT) and implies “subaging” ($\tau \approx t_w^{\delta }$ with $\delta <1$) as a broad equilibration crossover for quenches to nearly arrested equilibrium states; “hyperaging” (or superaging, $\tau \sim t_w^{{\delta }^{\prime }}$ with ${\delta }^{\prime }>1$) emerges for quenches deep into the glass. The latter is cut off by non-mean-field fluctuations that we account for within a recent extension of MCT, the stochastic $\beta$-relaxation theory (SBR). We exemplify the scaling laws with a schematic model that quantitatively fits simulation data.

Cytoskeletal networks form complex intracellular structures. Here we investigate a minimal model for filament-motor mixtures in which motors act as depolymerases and thereby regulate filament length. Combining agent-based simulations and hydrodynamic equations, we show that resource-limited length regulation drives the formation of filament clusters despite the absence of mechanical interactions between filaments. Even though the orientation of individual remains fixed, collective filament orientation emerges in the clusters, aligned orthogonal to their interfaces.

Diffusion of proteins on length scales of their size is crucial for understanding the machinery of living cells. X-ray photon correlation spectroscopy (XPCS) is currently the only way to access long-time collective diffusion on these length scales, but radiation damage so far limits the use in biological systems. We apply a new approach to use XPCS to measure cage relaxation in crowded $\alpha$-crystallin solutions. This allows us to correct for radiation effects, obtain missing information on long time diffusion, and support the fundamental analogy between protein and colloid dynamical arrest.

The Onsager reciprocity relations were formulated in the context of irreversible thermodynamics, but they are based on assumptions that have a wider applicability. Here, we present simulations testing the Onsager relations between surface-coupled diffusive and bulk fluxes in a system prepared in a nonequilibrium steady state. The system consists of a mixture of two identical species maintained at different temperatures inside a channel. In order to tune the friction of the two species with the walls independently, while keeping the particle-wall interaction potentials the same, we allow the kinematics of particle-wall collisions to be different: “bounce-back” ($B$) or “specular” ($S$). In the $BB$ case, diffusio-capillary transport can only take place if the two species have different temperatures. We find that the Onsager reciprocity relations are obeyed in the linear regime, even in the $BB$ case where all fluxes are the result of perturbing the system from a nonequilibrium steady state in a way that does not satisfy time-reversal symmetry. Our Letter provides a direct, numerical illustration of the validity of the Onsager relations outside their original range of application, and suggests their relevance for transport in driven or active systems.

Polycatenanes, macrochains of topologically interlocked rings with unique physical properties have recently gained considerable interest in supramolecular chemistry, biology, and soft matter. Most of the work has been, so far, focused on linear chains and on their variety of conformational properties compared to standard polymers. Here we go beyond the linear case and show that, by circularizing such macrochains, one can exploit the topology of the local interlockings to store twist in the system, significantly altering its metric and local properties. Moreover, by properly defining the twist (Tw) and writhe (Wr) of these macrorings we show the validity of a relation equivalent to the Călugăreanu-White-Fuller theorem $Tw+Wr=\mathrm{const}$, originally proved for ribbonlike structures such as double stranded DNA. Our results suggest that circular polycatenanes with storable and tunable twist can form a new category of highly designable multiscale structures with potential applications in supramolecular chemistry and material science.

Upon loading, amorphous solids can exhibit brittle yielding, with the abrupt formation of macroscopic shear bands leading to fracture, or ductile yielding, with a multitude of plastic events leading to homogeneous flow. It has been recently proposed, and subsequently questioned, that the two regimes are separated by a sharp critical point, as a function of some control parameter characterizing the intrinsic disorder strength and the degree of stability of the solid. In order to resolve this issue, we have performed extensive numerical simulations of athermally driven elastoplastic models with long-range and anisotropic realistic interaction kernels in two and three dimensions. Our results provide clear evidence for a finite-disorder critical point separating brittle and ductile yielding, and we provide an estimate of the critical exponents in 2D and 3D.

We investigate the dynamics of a single chiral active particle subject to an external torque due to the presence of a gravitational field. Our computer simulations reveal an arbitrarily strong increase of the long-time diffusivity of the gravitactic agent when the external torque approaches the intrinsic angular drift. We provide analytic expressions for the mean-square displacement in terms of eigenfunctions and eigenvalues of the noisy-driven-pendulum problem. The pronounced maximum in the diffusivity is then rationalized by the vanishing of the lowest eigenvalues of the Fokker-Planck equation for the angular motion as the rotational diffusion decreases and the underlying classical bifurcation is approached. A simple harmonic-oscillator picture for the barrier-dominated motion provides a quantitative description for the onset of the resonance while its range of validity is determined by the crossover to a critical-fluctuation-dominated regime.

The out of equilibrium nature of active systems can be exploited for the design of information-based engines. We design two types of an active Szilard engine that use a Maxwell demon to extract work from an active bath composed of noninteracting active Brownian particles. The two engines exploit either the quasistatic active pressure of active Brownian particles or the long correlation time of their velocities. For both engines the active bath allows us to overcome the Landauer principle and to extract larger work compared to conventional Szilard engines operating in quasithermal equilibrium. For both of our engines, we identify the optimal regimes at which the work extracted and the efficiency are maximized. Finally, we discuss them in the context of synthetic and biological active systems.

Inhomogeneity of ion correlation widely exists in many physicochemical, soft matter, and biological systems. Here, we apply the modified Gaussian renormalized fluctuation theory to study the classic example of the vapor-liquid interface of ionic fluids. The ion correlation is decomposed into a short-range contribution associated with the local electrostatic environment and a long-range contribution accounting for the spatially varying ionic strength and dielectric permittivity. For symmetric salt, both the coexistence curve and the interfacial tension predicted by our theory are in quantitative agreement with simulation data reported in the literature. Furthermore, we provide the first theoretical prediction of interfacial structure for asymmetric salt, highlighting the importance of capturing local charge separation.