Scaling and Localization in Multipole-Conserving Diffusion.

We study diffusion in systems of classical particles whose dynamics conserves the total center of mass. This conservation law leads to several interesting consequences. In finite systems, it allows for equilibrium distributions that are exponentially localized near system boundaries. It also yields an unusual approach to equilibrium, which in d dimensions exhibits scaling with dynamical exponent z=4+d. Similar phenomena occur for dynamics that conserves higher moments of the density, which we systematically classify using a family of nonlinear diffusion equations. In the quantum setting, analogous fermionic systems are shown to form real-space Fermi surfaces, while bosonic versions display a real-space analog of Bose-Einstein condensation.

Defect turbulence in a dense suspension of polar, active swimmers.

We study the effects of inertia in dense suspensions of polar swimmers. The hydrodynamic velocity field and the polar order parameter field describe the dynamics of the suspension. We show that a dimensionless parameter R (ratio of the swimmer self-advection speed to the active stress invasion speed [Phys. Rev. X 11, 031063 (2021)2160-330810.1103/PhysRevX.11.031063]) controls the stability of an ordered swimmer suspension. For R smaller than a threshold R_{1}, perturbations grow at a rate proportional to their wave number q. Beyond R_{1} we show that the growth rate is O(q^{2}) until a second threshold R=R_{2} is reached. The suspension is stable for R>R_{2}. We perform direct numerical simulations to characterize the steady-state properties and observe defect turbulence for R

Metastability of Discrete-Symmetry Flocks.

We study the stability of the ordered phase of flocking models with a scalar order parameter. Using both the active Ising model and a hydrodynamic description, we show that droplets of particles moving in the direction opposite to that of the ordered phase nucleate and grow. We characterize analytically this self-similar growth and demonstrate that droplets spread ballistically in all directions. Our results imply that, in the thermodynamic limit, discrete-symmetry flocks-and, by extension, continuous-symmetry flocks with rotational anisotropy-are metastable in all dimensions.

Target searches of interacting Brownian particles in dilute systems.

We study the target searches of interacting Brownian particles in a finite domain, focusing on the effect of interparticle interactions on the search time. We derive the integral equation for the mean first-passage time and acquire its solution as a series expansion in the orders of the Mayer function. We analytically obtain the leading order correction to the search time for dilute systems, which are most relevant to target search problems and prove a universal relation given by the particle density and the second virial coefficient. Finally, we validate our theoretical prediction by Langevin dynamics simulations for the various types of the interaction potential.

Model-Free Measurement of Local Entropy Production and Extractable Work in Active Matter.

Time-reversal symmetry breaking and entropy production are universal features of nonequilibrium phenomena. Despite its importance in the physics of active and living systems, the entropy production of systems with many degrees of freedom has remained of little practical significance because the high dimensionality of their state space makes it difficult to measure. Here we introduce a local measure of entropy production and a numerical protocol to estimate it. We establish a connection between the entropy production and extractability of work in a given region of the system and show how this quantity depends crucially on the degrees of freedom being tracked. We validate our approach in theory, simulation, and experiments by considering systems of active Brownian particles undergoing motility-induced phase separation, as well as active Brownian particles and E.coli in a rectifying device in which the time-reversal asymmetry of the particle dynamics couples to spatial asymmetry to reveal its effects on a macroscopic scale.

Optimal searcher distribution for parallel random target searches.

We consider a problem of finding a target located in a finite d-dimensional domain, using N independent random walkers, when partial information about the target location is given as a probability distribution. When N is large, the first-passage time sensitively depends on the initial searcher distribution, which invokes the question of the optimal searcher distribution that minimizes the first-passage time. Here, we analytically derive the equation for the optimal distribution and explore its limiting expressions. If the target volume can be ignored, the optimal distribution is proportional to the target distribution to the power of one third. If we consider a target of a finite volume and the probability of the initial overlapping of searchers with the target cannot be ignored in the large N limit, the optimal distribution has a weak dependence on the target distribution, with its variation being proportional to the logarithm of the target distribution. Using Langevin dynamics simulations, we numerically demonstrate our predictions in one and two dimensions.

Disordered boundaries destroy bulk phase separation in scalar active matter.

We show that disordered boundaries destroy bulk phase separation in scalar active systems in dimension d

Disorder-Induced Long-Ranged Correlations in Scalar Active Matter.

We study the impact of quenched random potentials and torques on scalar active matter. Microscopic simulations reveal that motility-induced phase separation is replaced in two dimensions by an asymptotically homogeneous phase with anomalous long-ranged correlations and nonvanishing steady-state currents. Using a combination of phenomenological models and a field-theoretical treatment, we show the existence of a lower-critical dimension d_{c}=4, below which phase separation is only observed for systems smaller than an Imry-Ma length scale. We identify a weak-disorder regime in which the structure factor scales as S(q)∼1/q^{2}, which accounts for our numerics. In d=2, we predict that, at larger scales, the behavior should cross over to a strong-disorder regime. In d>2, these two regimes exist separately, depending on the strength of the potential.

Rodlike Counterions near Charged Cylinders: Counterion Condensation and Intercylinder Interaction.

We study a system composed of like-charged cylinders and dumbbell-like counterions, with the focus laid on the role of the internal structure of counterions, using Monte Carlo simulations. The dumbbell ions are found to exhibit novel counterion condensation behavior governed by their length. Effective electrostatic interactions mediated between charged parallel cylinders also turn out significantly different from the case of pointlike ions, as a result of the complex interplay between the spatially separated charge distribution in the dumbbell counterions, their orientation, and the curvature of the charged cylinder. We show that at a weak-to-moderate electrostatic coupling strength, where effective like-charge interactions are usually found to be repulsive, the intercylinder interaction can become attractive and display a distinctive sensitivity to the cylinder curvature and dumbbell size, proving the significant effect of ion structure.

Interactions between a fluctuating polymer barrier and transport factors together with enzyme action are sufficient for selective and rapid transport through the nuclear pore complex.

The nuclear pore complex, the only pathway for transport between the nucleus and cytoplasm, functions as a highly selective gate that blocks nonspecific macromolecules while allowing the rapid transport of tagged [transport factor (TF) bound] cargo up to an order of magnitude larger. The mechanism of this gate's operation is not yet fully understood and progress has been primarily hindered by the inherent complexity and multiscale nature of the problem. One needs to consider the hundreds of disordered proteins (phenylalanine glycine nucleoporins or FG nups) lining the pore, as well as their overall architecture and dynamics at the microsecond scale, while also accounting for transport at the millisecond scale across the entire pore. Here we formulate an approach that addresses transport properties over a large range of length and time scales. We do this by incorporating microscopic biophysical details, such as charge and specific TF-FG nup interactions, to compute the free energy landscape encountered by the cargo. We connect this to macroscopic transport by treating cargo translocation as a stochastic barrier crossing process and computing the current and the translocation time. We then identify distinct transport regimes (fast permeable, slow permeable, and impermeable) determined by the cargo size, TF affinity for FG nups, and the activity of the enzymes that cleave TFs from cargo. Our results, therefore provide an integrated picture of transport through the NPC, while highlighting how FG nup interactions with TFs and enzyme activity cooperate to produce selectivity and efficiency.

Parallel random target searches in a confined space.

We study a random target searching performed by N independent searchers in a d-dimensional domain of a large but finite volume. Considering the two initial distributions of searchers where searchers are either uniformly or point distributed, we estimate the mean time for the first of the searchers to reach the target and refer to it as searching time. The searching time for the uniformly distributed searchers exhibits a universal power-law dependence on N, irrespective of dimensionality and the target-to-domain size ratio. For point-distributed searching, the searching time has a logarithmic dependence on N in the large N limit, while in the small N limit, it shows qualitatively different behaviors depending upon r_{0}, the initial distance of the searchers from a target. We obtain a diagram by comparing the searching times of the two initial distributions in the parameter space (r_{0},N) and therein present the asymptotic lines separating three characteristic regions to explain numerical simulation results.

Chiral Separation by Flows: The Role of Flow Symmetry and Dimensionality.

Separation of enantiomers by flows is a promising chiral resolution method using cost-effective microfluidics. Notwithstanding a number of experimental and numerical studies, a fundamental understanding still remains elusive, and an important question as to whether it is possible to specify common physical properties of flows that induce separation has not been addressed. Here, we study the separation of rigid chiral objects of an arbitrary shape induced by a linear flow field at low Reynolds numbers. Based on a symmetry property under parity inversion, we show that the rate-of-strain field is essential to drift the objects in opposite directions according to chirality. From eigenmode analysis, we also derive an analytic expression for the separation conditions which shows that the flow field should be quasi-two-dimensional for the precise and efficient resolutions of microscopic enantiomers. We demonstrate this prediction by Langevin dynamics simulations with hydrodynamic interactions fully implemented. Finally, we discuss the practical feasibility of the linear flow analysis, considering separations by a vortex flow or an extensional flow under a confining potential.