Search

Humanity has long sought inspiration from nature to innovate materials and devices. As science advances, nature-inspired materials are becoming part of our lives. Animate materials, characterized by their activity, adaptability, and autonomy, emulate properties of living systems. While only biological materials fully embody these principles, artificial versions are advancing rapidly, promising transformative impacts across various sectors. This roadmap presents authoritative perspectives on animate materials across different disciplines and scales, highlighting their interdisciplinary nature and potential applications in diverse fields including nanotechnology, robotics and the built environment. It underscores the need for concerted efforts to address shared challenges such as complexity management, scalability, evolvability, interdisciplinary collaboration, and ethical and environmental considerations. The framework defined by classifying materials based on their level of animacy can guide this emerging field encouraging cooperation and responsible development. By unravelling the mysteries of living matter and leveraging its principles, we can design materials and systems that will transform our world in a more sustainable manner.

Sand is a highly dissipative system, where the local spatial arrangements and densities depend strongly on the applied forces, resulting in fluid-like or solid-like behaviour. This makes sand swimming challenging and intriguing, raising questions about the nature of the motion and how to optimize the design of artificial swimmers able to swim in sand. Recent experiments suggest that lateral undulatory motion enables efficient locomotion, with a non-monotonic dependence of the swimming speed on the undulatory frequency and the height of the sediment bed. Here, we propose a quasi-2D granular model, where the effect of the bed height is modeled by a coarse-grained frictional force with the substrate. We show that the optimal frequency coincides with the second vibrational mode of the swimmer and explain the underlying mechanism through a characterization of the rheology of the medium. Potential implications in the design of artificial swimmers are discussed.

Collective cellular movement plays a crucial role in many processes fundamental to health, including development, reproduction, infection, wound healing, and cancer. The emergent dynamics that arise in these systems are typically thought to depend on how cells interact with one another and the mechanisms used to drive motility, both of which exhibit remarkable diversity across different biological systems. Here, we report experimental evidence of a universal feature in the patterns of flow that spontaneously emerges in groups of collectively moving cells. Specifically, we demonstrate that the flows generated by collectively moving dog kidney cells, human breast cancer cells, and by two different strains of pathogenic bacteria, all exhibit conformal invariance. Remarkably, not only do our results show that all of these very different systems display robust conformal invariance, but we also discovered that the precise form of the invariance in all four systems is described by the Schramm-Loewner Evolution (SLE), and belongs to the percolation universality class. A continuum model of active matter can recapitulate both the observed conformal invariance and SLE form found in experiments. The presence of universal conformal invariance reveals that the macroscopic features of living biological matter exhibit universal translational, rotational, and scale symmetries that are independent of the microscopic properties of its constituents. Our results show that the patterns of flows generated by diverse cellular systems are highly conserved and that biological systems can unexpectedly be used to experimentally test predictions from the theories for conformally invariant structures

Persistent motion of microswimmers near boundaries is known to result in surface accumulation. Recently it was shown experimentally that surface accumulation of microswimmers is impacted primarily by steric forces and short-ranged hydrodynamics. A way to control surface accumulation is by reducing the contact surface area between swimmers and surface by modifying its topography, typically through the application of microscale structure. In this work, we introduce a three-dimensional(3D) model of a microswimmer navigating a volume bounded by a top and bottom surface. We describe the swimmer-surface interaction with an effective short-ranged hydrodynamic alignment force, and study numerically the effect of surface textures, modelled by convex obstacles, on the surface accumulation of chiral and non-chiral microswimmers. We find that, depending on the angular velocity of the swimmer, and the alignment force, convex obstacles can either hinder or enhance surface accumulation. We discuss potential applications to sorting of microswimmers by their angular velocity.

With the rapid development of studies involving droplet microfluidics, drug delivery, cell detection, and microparticle synthesis, among others, many scientists have invested significant efforts to model the flow of these fluid-filled bodies. Motivated by the intricate coupling between hydrodynamics and the interactions of fluid-filled bodies, several methods have been developed. The objective of this review is to present a compact foundation of the methods used in the literature in the context of lattice Boltzmann methods. For hydrodynamics, we focus on the lattice-Boltzmann method due to its specific ability to treat time and spatial-dependent boundary conditions and to incorporate new physical models in a computationally efficient way. We split the existing methods into two groups with regard to the interfacial boundary: fluid-structure and fluid-fluid methods. The fluid-structure methods are characterised by the coupling between fluid dynamics and mechanics of the flowing body, often used in applications involving membranes and similar flexible solid boundaries. We further divide fluid-structure-based methods into two subcategories, those which treat the fluid-structure boundary as a continuum medium and those that treat it as a discrete collection of individual springs and particles. Next, we discuss the fluid-fluid methods, particularly useful for the simulations of fluid-fluid interfaces. We focus on models for immiscible droplets and their interaction in a suspending fluid and describe benchmark tests to validate the models for fluid-filled bodies.

Inverse design aims at the development of elementary building blocks that organize spontaneously into target shapes. In self-assembly, the blocks diffuse to their target position. Alternatively, recent experiments point to a more robust process in which the shape is formed from the self-folding of a planar template. To control the folding of templates with competing folded structures, we propose the inclusion of bond specificity. We consider a template that can fold into an octahedron or a boat shell and find the minimal design capable of targeting either shell or switching between the two through an external stimulus, adding a new dimension to the design of shape-changing materials.

Kirigami-inspired designs can enable self-folding three-dimensional materials from flat, two-dimensional sheets. Hierarchical designs of connected levels increase the diversity of possible target structures, yet they can lead to longer folding times in the presence of fluctuations. Here, we study the effect of rotational coupling between levels on the self-folding of two-level kirigami designs driven by thermal noise in a fluid. Naturally present due to hydrodynamic resistance, we find that optimization of this coupling as control parameter can significantly improve a structure's self-folding rate and yield.

The formation of groups of interacting individuals improves performance and fitness in many decentralised systems, from micro-organisms to social insects, from robotic swarms to artificial intelligence algorithms. Often, group formation and high-level coordination in these systems emerge from individuals with limited information-processing capabilities implementing low-level rules of communication to signal to each other. Here, we show that, even in a community of clueless individuals incapable of processing information and communicating, a dynamic environment can coordinate group formation by transiently storing memory of the earlier passage of individuals. Our results identify a new mechanism of indirect coordination via shared memory that is primarily promoted and reinforced by dynamic environmental factors, thus overshadowing the need for any form of explicit signalling between individuals. We expect this pathway to group formation to be relevant for understanding and controlling self-organisation and collective decision making in both living and artificial active matter in real-life environments.

Several natural and industrial processes involve the extraction or retention of a solute by a fluid invading a network of channels. Examples include aquifer contamination, chemical filtration, and coffee extraction. We propose a continuum equation to model these processes, parametrized by the P\'eclet number and the rate of mass transfer between the solid and the fluid. We study the time dependence of the extracted mass for different values of the parameter space. The continuum description is validated by combining extraction experiments with coffee and computational fluid dynamics. An analytical solution is derived for the limit of slow mass transfer, which is corroborated by numerical simulations.

The flow of deformable particles, such as droplets, dragged by a fluid, through a network of narrow pores inside rocks or other porous media is key in a range of applications, from enhanced oil recovery and water filtration to lab on a chip sorting of cells. The collective dynamics and its impact on the flow are poorly understood. Here, using droplets as a prototype, we show that collective transport can occur for conditions under which a single particle would get trapped at a pore channel. When a series of droplets gets trapped, the fluids flow is affected significantly, leading to an increase of the pressure difference across the pore channels, which in turn squeezes the particles through the channels. We analyze the conditions for a single droplet to flow through one pore and derive the corresponding Bond number. We also obtain a rule for the collective flow of droplets in porous media.

The study of the interactions of living adherent cells with mechanically stable (visco)elastic materials enables understanding and exploiting physiological phenomena mediated by cell-extracellular communication. However, insight on the interaction of cells and surrounding objects with different stability patterns upon cell contact might unveil cell responses that may be engineered for innovative applications. Here, it is hypothesized that the efficiency of cell attachment, spreading and movement across a free-packed granular bed of microparticles depend on microparticle diameter, raising the possibility of a necessary minimum traction force for the reinforcement of cell-particle bonds, and long-term cell adhesion. The results suggest that microparticles with 14-20 {\mu}m are prone to cell-mediated mobility, holding the potential of inducing early cell detachment, while objects with diameters from 38-85 {\mu}m enable long-lasting cell adhesion and proliferation. An in-silico hybrid particle-based model that addresses time-dependent biological mechanisms of cell adhesion is proposed, providing inspiration for engineering platforms to address healthcare-related challenges.

The presence of obstacles is intuitively expected to hinder the diffusive transport of micro-swimmers. However, for chiral micro-swimmers, a low density of obstacles near a surface can enhance their diffusive behavior, due to the rectification of the chiral motion by the obstacles. Here, we study numerically the role that disorder plays in determining the transport dynamics of chiral micro-swimmers on surfaces with obstacles. We consider different densities of regularly spaced obstacles and distinct types of disorder: noise in the dynamics of the micro-swimmer, quenched noise in the positions of the obstacles as well as obstacle size polydispersity. We show that, depending on the type and strength of the disorder, the presence of obstacles can either enhance or hinder transport, and discuss implications for the control of active transport in disordered media.

Self-organisation is the spontaneous emergence of spatio-temporal structures and patterns from the interaction of smaller individual units. Examples are found across many scales in very different systems and scientific disciplines, from physics, materials science and robotics to biology, geophysics and astronomy. Recent research has highlighted how self-organisation can be both mediated and controlled by confinement. Confinement occurs through interactions with boundaries, and can function as either a catalyst or inhibitor of self-organisation. It can then become a means to actively steer the emergence or suppression of collective phenomena in space and time. Here, to provide a common framework for future research, we examine the role of confinement in self-organisation and identify overarching scientific challenges across disciplines that need to be addressed to harness its full scientific and technological potential. This framework will not only accelerate the generation of a common deeper understanding of self-organisation but also trigger the development of innovative strategies to steer it through confinement, with impact, e.g., on the design of smarter materials, tissue engineering for biomedicine and crowd management.

Efficient nutrient mixing is crucial for the survival of bacterial colonies and other living systems. This raises the question of whether the optimization of mixing through the emergence of active turbulent motion in bacterial swarms played a role in the evolution of bacterial shapes. Here, to address this question, we solve the hydrodynamic equation for active nematics coupled with an advection-diffusion equation for the nutrients. The latter models a conserved activity field and mimics the conservation of nutrients in bacterial swarms. At the interface between active and passive nematic phases, in addition to diffusion, the activity is transported by interfacial flows and in turn modifies them through active stresses. We find that the interfacial dispersion of the conserved activity is subdiffusive due to the emergence of a barrier of negative defects at the active-passive interface, which hinders the propagation of the motile positive defects. Furthermore, we observe a non-monotonic dependence of the generalized diffusion coefficient on the aligning parameter, which is related to the shape of the particles. Our simulations suggest that there is an optimal shape that maximizes the dispersion of conserved activity at the active-passive nematic interface.

We study the competitive irreversible adsorption of a binary mixture of monomers and square-shaped particles of linear size $R$ on the square lattice. With the random sequential adsorption model, we investigate how the jamming coverage and percolation properties depend on the size ratio $R$ and relative flux $F$. We find that the onset of percolation of monomers is always lower for the binary mixture than in the case with only monomers ($R=1$). Moreover, for values $F$ below a critical value, the higher is the flux or size of the largest species, the lower is the value of the percolation threshold for monomers., Comment: Submitted to a special issue of Journal of Physics A: Mathematical and Theoretical to mark the 70th Birthday of Robert M. Ziff

On a daily commute, the shortest route from home to work rarely overlaps completely the shortest way back. We analyze this asymmetry for several cities and show that it exists even without traffic, due to a non-negligible fraction of one-way streets. For different pairs of origin-destination ($\rm OD$), we compute the log-ratio $r=\ln(\ell_{\rm D}/\ell_{\rm O})$, where $\ell_{\rm O}$ and $\ell_{\rm D}$ are the lengths of the shortest routes from $\rm O$ to $\rm D$ and from $\rm D$ to $\rm O$, respectively. While its average is zero, the amplitude of the fluctuations decays as a power law of the $\rm OD$ shortest path length, $r\sim \ell_{\rm O}^{-\beta}$. Similarly, the fraction of one-way streets in a shortest route also decays as $\ell_{\rm O}^{-\alpha}$. Based on semi-analytic arguments, we show that $\beta=(1+\alpha)/2$. Thus, the value of the exponent $\beta$ is related to correlations in the structure of the underlying street network., Comment: 6 pages, 4 figures

In vivo and in vitro cells rely on the support of an underlying biocompatible substrate, such as the extracellular matrix or a culture substrate, to spread and proliferate. The mechanical and chemical properties of such structures play a central role in the dynamical and statistical properties of the tissue. At the cell scale, these substrates are highly disordered. Here, we investigate how spatial heterogeneities of the cell-substrate interaction influence the motility of the cells in a model confluent tissue. We use the Self-Propelled Voronoi model and describe the disorder as a spatially dependent preferred geometry of the individual cells. We found that when the characteristic length scale of the preferred geometry is smaller than the cell size, the tissue is less rigid than its homogeneous counterpart, with a consequent increase in cell motility. This result is in sharp contrast to what has been reported for tissues with heterogeneity in the mechanical properties of the individual cells, where the disorder favors rigidity. Using the fraction of rigid cells, we observe a collapse of the motility data for different model parameters and provide evidence that the rigidity transition in the model tissue is accompanied by the emergence of a spanning cluster of rigid cells.

We consider a simple theoretical model to investigate the impact of inheritances on the wealth distribution. Wealth is described as a finite resource, which remains constant over different generations and is divided equally among offspring. All other sources of wealth are neglected. We consider different societies characterized by a different offspring probability distribution. We find that, if the population remains constant, the society reaches a stationary wealth distribution. We show that inequality emerges every time the number of children per family is not always the same. For realistic offspring distributions from developed countries, the model predicts a Gini coefficient of $G\approx 0.3$. If we divide the society into wealth classes and set the probability of getting married to depend on the distance between classes, the stationary wealth distribution crosses over from an exponential to a power-law regime as the number of wealth classes and the level of class distinction increase., Comment: 7 pages, 7 figures

Soft particles suspended in fluids flowing through microchannels are often encountered in biological flows such as cells in vessels. They can deform under flow or when subject to mechanical stresses as they interact with themselves. Deformability adds distinct characteristics to the nature of the flow of these particles. We simulate deformable particles suspended in a fluid at a high volume fraction flowing through a 3D wetting channel. We find a discontinuous shear thinning transition, which depends on the particle deformability. The capillary number is the main parameter that controls the transition. We also notice significant differences in the velocity profiles when compared to their 2D counterpart. To obtain these results, we improved and extended to 3D a multicomponent Lattice Boltzmann method which prevents the coalescence between the droplets., Comment: 8 pages, 7 figures

The flow through a porous medium strongly depends on the boundary conditions, very often assumed to be static. Here, we consider changes in the medium due to swelling and erosion and extend existing Lattice-Boltzmann models to include both. We study two boundary conditions: a constant pressure drop and a constant flow rate. For a constant flow rate, the steady state depends solely on the erosion dynamics while for a constant pressure drop it depends also on the timescale of swelling. We analyze the competition between swelling and erosion and identify a transition between regimes where either swelling or erosion dominate., Comment: 20 pages, 7 figures

Activity in nematics drives interfacial flows that lead to preferential alignment that is tangential or planar for extensile systems (pushers) and perpendicular or homeotropic for contractile ones (pullers). This alignment is known as active anchoring and has been reported for a number of systems and described using active nematic hydrodynamic theories. The latter are based on the one-elastic constant approximation, i.e., they assume elastic isotropy of the underlying passive nematic. Real nematics, however, have different elastic constants, which lead to interfacial anchoring. In this paper, we consider elastic anisotropy in multiphase and multicomponent hydrodynamic models of active nematics and investigate the competition between the interfacial alignment driven by the elastic anisotropy of the passive nematic and the active anchoring. We start by considering systems with translational invariance to analyze the alignment at flat interfaces and, then, consider two-dimensional systems and active nematic droplets. We investigate the competition of the two types of anchoring over a wide range of the other parameters that characterize the system. The results of the simulations reveal that the active anchoring dominates except at very low activities, when the interfaces are static. In addition, we found that the elastic anisotropy does not affect the dynamics but changes the active length that becomes anisotropic.

The implementation of boundary conditions is among the most challenging parts of modeling fluid flow through channels and complex media. Here, we show that the existing methods to deal with liquid-wall interactions using multicomponent Lattice Boltzmann are accurate when the wall is aligned with the main axes of the lattice but fail otherwise. To solve this problem, we extend a strategy previously developed for multiphase models. As an example, we study the coalescence of two droplets on a curved surface in two dimensions. The strategy proposed here is of special relevance for binary flows in complex geometries., Comment: 12 pages, 8 figures

The structure and dynamics of tissue cultures depend strongly on the physical and chemical properties of the underlying substrate. Inspired by previous advances in the context of inorganic materials, the use of patterned culture surfaces has been proposed as an effective way to induce space-dependent properties in cell tissues. However, cells move and diffuse and the transduction of external stimuli to biological signals is not instantaneous. Here, we show that the fidelity of patterns depends on the relation between the diffusion ($\tau_D$) and adaptation ($\tau$) times. Numerical results for the self-propelled Voronoi model reveal that the fidelity decreases with $\tau/\tau_D$, a result that is reproduced by a continuum reaction-diffusion model. We derive a minimum length scale for the patterns that depends on $\tau/\tau_D$ and can be much larger than the cell size.

Motivated by multiphase flow in reservoirs, we propose and study a two-species sandpile model in two dimensions. A pile of particles becomes unstable and topples if, at least one of the following two conditions is fulfilled: 1) the number of particles of one species in the pile exceeds a given threshold or 2) the total number of particles in the pile exceeds a second threshold. The latter mechanism leads to the invasion of one species through regions dominated by the other species. We studied numerically the statistics of the avalanches and identified two different regimes. For large avalanches the statistics is consistent with ordinary Bak-Tang-Weisenfeld model. Whereas, for small avalanches, we find a regime with different exponents. In particular, the fractal dimension of the external perimeter of avalanches is $D_f=1.47\pm 0.02$ and the exponent of their size distribution exponent is $\tau_s=0.95\pm 0.03$, which are significantly different from $D_f=1.25\pm 0.01$ and $\tau_s=1.26\pm 0.04$, observed for large avalanches.

Motivated by results for the propagation of active-passive interfaces of bacterial Serratia marcescens swarms [Nat. Comm., 9, 5373 (2018)] we use a hydrodynamic multiphase model to investigate the propagation of interfaces of active nematics on substrates. We characterize the active nematic phase of the model and discuss its description of the statistical dynamics of the swarms. We show that the velocicity correlation functions and the energy spectrum of the active turbulent phase of the model scale with the active length, for a range of activities. In addition, the energy spectrum exhibits two power-law regimes with exponents close to those reported for other models and bacterial swarms. Although the exponent of the rising branch of the spectrum (small wavevector) appears to be independent of the activity, the exponent of the decay changes with activity systematically, albeit slowly. We characterize also the propagation of circular and flat active-passive interfaces. We find that the closing time of the circular passive domain decays quadratically with the activity and that the structure factor of the flat interface is similar to that reported for swarms, with an activity dependent exponent. Finally, the effect of the substrate friction was investigated. We found an activity dependent threshold, above which the turbulent active nematic forms isolated islands that shrink until the system becomes isotropic and below which the active nematic expands, with a well defined propagating interface. The interface may be stopped by a fricion gradient.

We use a hybrid method of lattice Boltzmann and finite differences to simulate flat and curved interfaces between the nematic and isotropic phases of a liquid crystal described by the Landau-de Gennes theory. For the flat interface, we measure the interfacial velocity at different temperatures around the coexistence. We show that the interface is completely static at the coexistence temperature and that the profile width is in line with the theoretical predictions. The interface is stable in a range of temperatures around coexistence and disappears when one of the two phases becomes mechanically unstable. We stabilize circular nematic domains by a shift in temperature, related to the Laplace pressure, and estimate the spurious velocities of these lattice Boltzmann simulations., Comment: Proceedings of the Portuguese Conference on Condensed Matter Physics

We use numerical simulations to investigate the hydrodynamic behavior of the interface between nematic (N) and isotropic (I) phases of a confined active liquid crystal. At low activities, a stable interface with constant shape and velocity is observed separating the two phases. For nematics in homeotropic channels, the velocity of the interface at the NI transition increases from zero (i) linearly with the activity for contractile systems and (ii) quadratically for extensile ones. Interestingly, the nematic phase expands for contractile systems while it contracts for extensile ones, as a result of the active forces at the interface. Since both activity and temperature affect the stability of the nematic, for active nematics in the stable regime the temperature can be tuned to observe static interfaces, providing an operational definition for the coexistence of active nematic and isotropic phases. At higher activities, beyond the stable regime, an interfacial instability is observed for extensile nematics. In this regime defects are nucleated at the interface and move away from it. The dynamics of these defects is regular and persists asymptotically for a finite range of activities. We used an improved hybrid model of finite differences and lattice Boltzmann method with multi-relaxation-time collision operator, the accuracy of which allowed us to characterize the dynamics of the distinct interfacial regimes., Comment: 11 pages, 13 figures

Structural defects are ubiquitous in condensed matter, and not always a nuisance. For example, they underlie phenomena such as Anderson localization and hyperuniformity, and they are now being exploited to engineer novel materials. Here, we show experimentally that the density of structural defects in a 2D binary colloidal crystal can be engineered with a random potential. We generate the random potential using an optical speckle pattern, whose induced forces act strongly on one species of particles (strong particles) and weakly on the other (weak particles). Thus, the strong particles are more attracted to the randomly distributed local minima of the optical potential, leaving a trail of defects in the crystalline structure of the colloidal crystal. While, as expected, the crystalline ordering initially decreases with increasing fraction of strong particles, the crystalline order is surprisingly recovered for sufficiently large fractions. We confirm our experimental results with particle-based simulations, which permit us to elucidate how this non-monotonic behavior results from the competition between the particle-potential and particle-particle interactions.

How motile bacteria move near a surface is a problem of fundamental biophysical interest and is key to the emergence of several phenomena of biological, ecological and medical relevance, including biofilm formation. Solid boundaries can strongly influence a cell's propulsion mechanism, thus leading many flagellated bacteria to describe long circular trajectories stably entrapped by the surface. Experimental studies on near-surface bacterial motility have, however, neglected the fact that real environments have typical microstructures varying on the scale of the cells' motion. Here, we show that micro-obstacles influence the propagation of peritrichously flagellated bacteria on a flat surface in a non-monotonic way. Instead of hindering it, an optimal, relatively low obstacle density can significantly enhance cells' propagation on surfaces due to individual forward-scattering events. This finding provides insight on the emerging dynamics of chiral active matter in complex environments and inspires possible routes to control microbial ecology in natural habitats.

The random sequential adsorption (RSA) model is a classical model in Statistical Physics for adsorption on two-dimensional surfaces. Objects are deposited sequentially at random and adsorb irreversibly on the landing site, provided that they do not overlap any previously adsorbed object. The kinetics of adsorption ceases when no more objects can be adsorbed (jamming state). Here, we investigate the role of post-relaxation on the jamming state and percolation properties of RSA of dimers on a two-dimensional lattice. We consider that, if the deposited dimer partially overlaps with a previously adsorbed one, a sequence of dimer displacements may occur to accommodate the new dimer. The introduction of this simple relaxation dynamics leads to a more dense jamming state than the one obtained with RSA without relaxation. We also consider the anisotropic case, where one dimer orientation is favored over the other, finding a non-monotonic dependence of the jamming coverage on the strength of anisotropy. We find that the density of adsorbed dimers at which percolation occurs is reduced with relaxation, but the value depends on the strength of anisotropy.

An extension of the Random Sequential Adsorption (RSA) model has been proposed recently, motivated by the coverage of oil droplets by DNA-functionalized colloidal particles. Particles arrive to a flat substrate with a uniform flux F but they can only adsorb on patches. Patches diffuse on the substrate with a diffusion coefficient D if they are free and they remain immobile when attached to an adsorbed particle. The adsorption is considered irreversible and particles cannot adsorb on top of each other. Thus, the system reaches a jammed state, consisting of a monolayer where no more particles can adsorb. We performed Monte Carlo simulations to study the adsorption kinetics and jammed-state morphology on a one-dimensional lattice. We show that, while the time-dependence of the coverage depends on F and D, the jammed-state coverage depends solely on the ratio F/D. This result is grasped by a simple mean-field calculation. We also report two different regimes for the functional dependence of the jammed-state coverage on the size of the particles, for low and high density of patches.

We study the percolation properties for a system of functionalized colloids on patterned substrates via Monte Carlo simulations. The colloidal particles are modeled as hard disks with three equally-distributed attractive patches on their perimeter. We describe the patterns on the substrate as circular potential wells of radius $R_p$ arranged in a regular square or hexagonal lattice. We find a nonmonotonic behavior of the percolation threshold (packing fraction) as a function of $R_p$. For attractive wells, the percolation threshold is higher than the one for clean (non-patterned) substrates if the circular wells are non-overlapping and can only be lower if the wells overlap. For repulsive wells we find the opposite behavior. In addition, at high packing fractions the formation of both structural and bond defects suppress percolation. As a result, the percolation diagram is reentrant with the non-percolated state occurring at very low and intermediate densities.

We consider mixtures of two species of spherical colloidal particles that differ in their hydrodynamic radii, but are otherwise identical, in the presence of an external field. Since the particle-particle and particle-field interactions are the same for both species, they are completely mixed in the thermodynamic limit in the presence of any static field. Here, we combine Brownian Dynamics and Dynamic Density Functional theory of fluids to show that for sufficiently large differences in the hydrodynamic radius of the particles (and corresponding differences in their electrophoretic mobilities) dynamical demixing is observed. These demixed states are transient but, under certain conditions, packing effects compromise the relaxation towards the thermodynamic states and the lifetime of the demixed phases increases significantly.

The self-assembly of functionalized (patchy) particles with directional interactions into target structures is still a challenge, despite the significant experimental advances on their synthesis. The self-assembly pathways are typically characterized by high energy barriers that hinder the access to stable (equilibrium) structures. A possible strategy to tackle this challenge is to perform annealing cycles. By periodically switching on and off the inter-particle bonds, one expects to smooth-out the kinetic pathways and favor the assembly of targeted structures. Preliminary results have shown that the efficiency of annealing cycles depends strongly on their frequency. Here, we study numerically how this frequency-dependence scales with the strength of the directional interactions (size of the patch $\sigma$). We use analytical arguments to show that the scaling results from the statistics of a random walk in configurational space.

Active systems are inherently out of equilibrium, as they collect energy from their surroundings and transform it into directed motion. A recent theoretical study suggests that binary mixtures of active particles with distinct effective diffusion coefficients exhibit dynamical demixing above a threshold ratio of the diffusion coefficients. Here, we show that this threshold may be reduced drastically in the presence of external fields. We investigate the demixing as a function of the ratio of the diffusion coefficients and discuss the implications of the results for active systems.

Connectivity - or the lack thereof - is crucial for the function of many man-made systems, from financial and economic networks over epidemic spreading in social networks to technical infrastructure. Often, connections are deliberately established or removed to induce, maintain, or destroy global connectivity. Thus, there has been a great interest in understanding how to control percolation, the transition to large-scale connectivity. Previous work, however, studied control strategies assuming unlimited resources. Here, we depart from this unrealistic assumption and consider the effect of limited resources on the effectiveness of control. We show that, even for scarce resources, percolation can be controlled with an efficient intervention strategy. We derive this strategy and study its implications, revealing a discontinuous transition as an unintended side-effect of optimal control., Comment: 5 pages, 4 figures, additional supplemental material (19 pages)

By analyzing a unique dataset of more than 270,000 scientists, we discovered substantial gender differences in scientific collaborations. While men are more likely to collaborate with other men, women are more egalitarian. This is consistently observed over all fields and regardless of the number of collaborators a scientist has. The only exception is observed in the field of engineering, where this gender bias disappears with increasing number of collaborators. We also found that the distribution of the number of collaborators follows a truncated power law with a cut-off that is gender dependent and related to the gender differences in the number of published papers. Considering interdisciplinary research, our analysis shows that men and women behave similarly across fields, except in the case of natural sciences, where women with many collaborators are more likely to have collaborators from other fields.

We report a study of reversible adsorption of DNA-coated colloids on complementary functionalized oil droplets. We show that it is possible to control the surface coverage of oil droplets by colloidal particles, by exploiting the fact that during slow adsorption, compositional arrest takes place well before structural arrest occurs. As a consequence, we can prepare colloid-coated oil droplets with a `frozen' degree of loading, but with fully ergodic colloidal dynamics on the droplets. We illustrate the equilibrium nature of the adsorbed colloidal phase by exploring the quasi two-dimensional (2d) phase behaviour of the adsorbed colloids under the influence of depletion interactions. We present simulations of a simple model that illustrates the nature of the compositional arrest and the structural ergodicity., Comment: 9 pages, 4 figures

Advances in mathematical epidemiology have led to a better understanding of the risks posed by epidemic spreading and informed strategies to contain disease spread. However, a challenge that has been overlooked is that, as a disease becomes more prevalent, it can limit the availability of the capital needed to effectively treat those who have fallen ill. Here we use a simple mathematical model to gain insight into the dynamics of an epidemic when the recovery of sick individuals depends on the availability of healing resources that are generated by the healthy population. We find that epidemics spiral out of control into "explosive" spread if the cost of recovery is above a critical cost. This can occur even when the disease would die out without the resource constraint. The onset of explosive epidemics is very sudden, exhibiting a discontinuous transition under very general assumptions. We find analytical expressions for the critical cost and the size of the explosive jump in infection levels in terms of the parameters that characterize the spreading process. Our model and results apply beyond epidemics to contagion dynamics that self-induce constraints on recovery, thereby amplifying the spreading process., Comment: 24 pages, 6 figures

The saturation length of aeolian sand transport ($L_s$), characterizing the distance needed by wind-blown sand to adapt to changes in the wind shear, is essential for accurate modeling of the morphodynamics of Earth's sandy landscapes and for explaining the formation and shape of sand dunes. In the last decade, it has become a widely-accepted hypothesis that $L_s$ is proportional to the characteristic distance needed by transported particles to reach the wind speed (the ``drag length''). Here we challenge this hypothesis. From extensive numerical Discrete Element Method simulations, we find that, for medium and strong winds, $L_s\propto V_s^2/g$, where $V_s$ is the saturated value of the average speed of sand particles traveling above the surface and $g$ the gravitational constant. We show that this proportionality is consistent with a recent analytical model, in which the drag length is just one of four similarly important length scales relevant for sand transport saturation., Comment: 5 pages, 2 figures

A multi-scale scheme for the invasion percolation of rock fracture networks with heterogeneous fracture aperture fields is proposed. Inside fractures, fluid transport is calculated on the finest scale and found to be localized in channels as a consequence of the aperture field. The channel network is characterized and reduced to a vectorized artificial channel network (ACN). Different realizations of ACNs are used to systematically calculate efficient apertures for fluid transport inside differently sized fractures as well as fracture intersection and entry properties. Typical situations in fracture networks are parameterized by fracture inclination, flow path length along the fracture and intersection lengths in the entrance and outlet zones of fractures. Using these scaling relations obtained from the finer scales, we simulate the invasion process of immiscible fluids into saturated discrete fracture networks, which were studied in previous works., Comment: 21 pages

We investigate the metallic breakdown of a substrate on which highly conducting particles are adsorbed and desorbed with a probability that depends on the local electric field. We find that, by tuning the relative strength $q$ of this dependence, the breakdown can change from continuous to explosive. Precisely, in the limit in which the adsorption probability is the same for any finite voltage drop, we can map our model exactly onto the $q$-state Potts model and thus the transition to a jump occurs at $q=4$. In another limit, where the adsorption probability becomes independent of the local field strength, the traditional bond percolation model is recovered. Our model is thus an example of a possible experimental realization exhibiting a truly discontinuous percolation transition.

Resilience of most critical infrastructures against failure of elements that appear insignificant is usually taken for granted. The World Airline Network (WAN) is an infrastructure that reduces the geographical gap between societies, both small and large, and brings forth economic gains. With the extensive use of a publicly maintained data set that contains information about airports and alternative connections between these airports, we empirically reveal that the WAN is a redundant and resilient network for long distance air travel, but otherwise breaks down completely due to removal of short and apparently insignificant connections. These short range connections with moderate number of passengers and alternate flights are the connections that keep remote parts of the world accessible. It is surprising, insofar as there exists a highly resilient and strongly connected core consisting of a small fraction of airports (around 2.3%) together with an extremely fragile star-like periphery. Yet, in spite of their relevance, more than 90% of the world airports are still interconnected upon removal of this core. With standard and unconventional removal measures we compare both empirical and topological perceptions for the fragmentation of the world. We identify how the WAN is organized into different classes of clusters based on the physical proximity of airports and analyze the consequence of this fragmentation., Comment: 8 pages, 4 figures, Suupporting Info

Power grids, road maps, and river streams are examples of infrastructural networks which are highly vulnerable to external perturbations. An abrupt local change of load (voltage, traffic density, or water level) might propagate in a cascading way and affect a significant fraction of the network. Almost discontinuous perturbations can be modeled by shock waves which can eventually interfere constructively and endanger the normal functionality of the infrastructure. We study their dynamics by solving the Burgers equation under random perturbations on several real and artificial directed graphs. Even for graphs with a narrow distribution of node properties (e.g., degree or betweenness), a steady state is reached exhibiting a heterogeneous load distribution, having a difference of one order of magnitude between the highest and average loads. Unexpectedly we find for the European power grid and for finite Watts-Strogatz networks a broad pronounced bimodal distribution for the loads. To identify the most vulnerable nodes, we introduce the concept of node-basin size, a purely topological property which we show to be strongly correlated to the average load of a node.