Brownian self-driven particles on the surface of a sphere

Effect of confinement and topology: 2-TIPS vs. MIPS

2-TIPS (two temperature induced phase separation) refers to the phase separation phenomenon observed in mixtures of active and passive particles which are modelled using scalar activity. The active particles are connected to a thermostat at high temperature while the passive particles are connected to the thermostat at low temperature and the relative temperature difference between "hot" and "cold" particles is taken as the measure of the activity χ of the non-equilibrium system. The study of such binary mixtures of hot and cold particles under various kinds of confinement is an important problem in many physical and biological processes. The nature and extent of phase separation are heavily influenced by the geometry of confinement, activity, and density of the non-equilibrium binary mixture. Investigating such 3D binary mixtures confined by parallel walls, we observe that the active and passive particles phase separate, but the extent of phase separation is reduced compared to bulk phase separation at high densities and enhanced at low densities. However, when the binary mixture of active and passive particles is confined inside a spherical cavity, the phase separation is radial for small radii of the confining sphere and the extent of phase separation is higher compared to their bulk counterparts. Confinement leads to interesting properties in the passive (cold) region like enhanced layering and high compression in the direction parallel to the confining wall. In 2D, both the bulk and confined systems of the binary mixture show a significant decrement in the extent of phase separation at higher densities. This observation is attributed to the trapping of active particles inside the passive cluster, which increases with density. Thus the 2D systems show structures more akin to dense-dilute phase co-existence, which is observed in motility induced phase separation in 2D active systems. The binary mixture constrained on the spherical surface also shows similar phase co-existence. Our analyses reveal that the coexistent densities observed in 2-TIPS on the spherical surface agree with the findings of previous studies on MIPS in active systems on a sphere.

Fractional and scaled Brownian motion on the sphere: The effects of long-time correlations on navigation strategies

We analyze fractional Brownian motion and scaled Brownian motion on the two-dimensional sphere S^{2}. We find that the intrinsic long-time correlations that characterize fractional Brownian motion collude with the specific dynamics (navigation strategies) carried out on the surface giving rise to rich transport properties. We focus our study on two classes of navigation strategies: one induced by a specific set of coordinates chosen for S^{2} (we have chosen the spherical ones in the present analysis), for which we find that contrary to what occurs in the absence of such long-time correlations, nonequilibrium stationary distributions are attained. These results resemble those reported in confined flat spaces in one and two dimensions [Guggenberger et al. New J. Phys. 21, 022002 (2019)1367-263010.1088/1367-2630/ab075f; Vojta et al. Phys. Rev. E 102, 032108 (2020)2470-004510.1103/PhysRevE.102.032108]; however, in the case analyzed here, there are no boundaries that affect the motion on the sphere. In contrast, when the navigation strategy chosen corresponds to a frame of reference moving with the particle (a Frenet-Serret reference system), then the equilibrium distribution on the sphere is recovered in the long-time limit. For both navigation strategies, the relaxation times toward the stationary distribution depend on the particular value of the Hurst parameter. We also show that on S^{2}, scaled Brownian motion, distinguished by a time-dependent diffusion coefficient with a power-scaling, is independent of the navigation strategy finding a good agreement between the analytical calculations obtained from the solution of a time-dependent diffusion equation on S^{2}, and the numerical results obtained from our numerical method to generate ensemble of trajectories.

Diffusion coefficients and MSD measurements on curved membranes and porous media

We study some geometric aspects that influence the transport properties of particles that diffuse on curved surfaces. We compare different approaches to surface diffusion based on the Laplace-Beltrami operator adapted to predict concentration along entire membranes, confined subdomains along surfaces, or within porous media. Our goal is to summarize, firstly, how diffusion in these systems results in different types of diffusion coefficients and mean square displacement measurements, and secondly, how these two factors are affected by the concavity of the surface, the shape of the possible barriers or obstacles that form the available domains, the sinuosity, tortuosity, and constrictions of the trajectories and even how the observation plane affects the measurements of the diffusion. In addition to presenting a critical and organized comparison between different notions of MSD, in this review, we test the correspondence between theoretical predictions and numerical simulations by performing finite element simulations and illustrate some situations where diffusion theory can be applied. We briefly reviewed computational schemes for understanding surface diffusion and finally, discussed how this work contributes to understanding the role of surface diffusion transport properties in porous media and their relationship to other transport processes.

Curvature-controlled geometrical lensing behavior in self-propelled colloidal particle systems

In many biological systems, the curvature of the surfaces cells live on influences their collective properties. Curvature should likewise influence the behavior of active colloidal particles. We show using molecular simulation of self-propelled active particles on surfaces of Gaussian curvature (both positive and negative) how curvature sign and magnitude can alter the system's collective behavior. Curvature acts as a geometrical lens and shifts the critical density of motility-induced phase separation (MIPS) to lower values for positive curvature and higher values for negative curvature, which we explain theoretically by the nature of parallel lines in spherical and hyperbolic space. Curvature also fluidizes dense MIPS clusters due to the emergence of defect patterns disrupting the crystalline order inside the clusters. Using our findings, we engineer three confining surfaces that strategically combine regions of different curvature to produce a host of novel dynamical behaviors, including cyclic MIPS on spherocylinders, directionally biased cyclic MIPS on spherocones, and position dependent cluster fluctuations on metaballs.

Activity-induced polar patterns of filaments gliding on a sphere

Active matter systems feature the ability to form collective patterns as observed in a plethora of living systems, from schools of fish to swimming bacteria. While many of these systems move in a wide, three-dimensional environment, several biological systems are confined by a curved topology. The role played by a non-Euclidean geometry on the self-organization of active systems is not yet fully understood, and few experimental systems are available to study it. Here, we introduce an experimental setup in which actin filaments glide on the inner surface of a spherical lipid vesicle, thus embedding them in a curved geometry. We show that filaments self-assemble into polar, elongated structures and that, when these match the size of the spherical geometry, both confinement and topological constraints become relevant for the emergent patterns, leading to the formation of polar vortices and jammed states. These results experimentally demonstrate that activity-induced complex patterns can be shaped by spherical confinement and topology.

A hard-sphere model of protein corona formation on spherical and cylindrical nanoparticles

A nanoparticle (NP) immersed in biological media rapidly forms a corona of adsorbed proteins, which later controls the eventual fate of the particle and the route through which adverse outcomes may occur. The composition and timescale for the formation of this corona are both highly dependent on both the NP and its environment. The deposition of proteins on the surface of the NP can be imitated by a process of random sequential adsorption, and, based on this model, we develop a rate-equation treatment for the formation of a corona represented by hard spheres on spherical and cylindrical NPs. We find that the geometry of the NP significantly alters the composition of the corona through a process independent of the rate constants assumed for adsorption and desorption of proteins, with the radius and shape of the NP both influencing the corona. We further investigate the roles of protein mobility on the surface of the NP and changes in the concentration of proteins.

Engineered swift equilibration for arbitrary geometries

Engineered swift equilibration (ESE) is a class of driving protocols that enforce an equilibrium distribution with respect to external control parameters at the beginning and end of rapid state transformations of open, classical nonequilibrium systems. ESE protocols have previously been derived and experimentally realized for Brownian particles in simple, one-dimensional, time-varying trapping potentials; one recent study considered ESE in two-dimensional Euclidean configuration space. Here we extend the ESE framework to generic, overdamped Brownian systems in arbitrary curved configuration space and illustrate our results with specific examples not amenable to previous techniques. Our approach may be used to impose the necessary dynamics to control the full temporal configurational distribution in a wide variety of experimentally realizable settings.

Effects of the curvature gradient on the distribution and diffusion of colloids confined to surfaces

The properties and behavior of colloids confined to move on curved surfaces offer a fertile ground for analysis since the geometric constraints induce specific features that are not available in flat spaces. Given their pertinence for biological and physicochemical processes, both with potential useful applications, the development of the concepts and methodology necessary for a deeper understanding of these unconventional systems is indeed an essential pursuit. The present study discusses a general and rigorous algorithm for the implementation of Brownian dynamics simulations that solves underlying difficulties and shortcomings inherent to conventional first-order schemes. Still based on the Ermak-McCammon recipe, our approach complements it with the higher-order geodesical projections of the elementary jumps generated on the associated tangent plane. This strategy, which warrants the locally isotropic propagation of non-interacting particles, is tested with a model system of colloidal particles interacting through a screened Coulomb potential while confined to move on ellipsoidal surfaces. This allows us to measure the effects prompted by the curvature gradient on the static and dynamic properties of this system. The varying curvature thus induces energetically favorable configurations in which the particles maximize their Euclidean distancing by crowding the regions with the largest Gaussian curvature, while withdrawing from those with the lowest. In turn, these inhomogeneous distributions provoke the anisotropic self-diffusion of the confined colloids, which is examined by exploiting the pertinent geodesic radial coordinates. The proficient methods under consideration thus allows dealing with the rich and remarkable new phenomena generated by any distinctive surface geometry.

Binary mixtures of active and passive particles on a sphere

We study the cooperation and segregation dynamics of binary mixtures of active and passive particles on a sphere. According to the competition between rotational diffusion and polar alignment, we find three distinct phases: a mixed phase and two different demixed phases. When rotational diffusion dominates the dynamics, the demixing is due to the aggregation of passive particles, where active and passive particles respectively occupy two hemispheres. When polar alignment is dominated, the demixing is caused by the aggregation of active particles, where active particles occupy the equator of the sphere and passive particles occupy the two poles of the sphere. In this case, there exist a circulating band cluster and two cambered surface clusters, which is a purely curvature-driven effect with no equivalent in the planar model. When rotational diffusion and polar alignment are comparable, particles are completely mixed. Our findings are relevant to the experimental pursuit of segregation dynamics of binary mixtures on curved surfaces.

Active matter on Riemannian manifolds

We formulate the dynamics of overdamped Brownian active particles (swimmers) moving on any Riemannian 2-manifold. To characterize such dynamics at short times, an analytical expression for the variance of swimmers diffusing on any Riemmanian 2-manifold is derived. To show the generality of the present work, we apply the latter dynamics to swimmers moving on the surface of a spheroid and a torus, and offer analytical expressions for both their long-time variances and steady angular marginal probability density functions. Finally, Brownian dynamics simulations are used to validate our theoretical findings.

Nanoparticles confined to a spherical surface in the presence of an external field: Interaction forces and induced microstructure

The structural response of a set of charged nanoparticles confined to move on a spherical surface under the influence of an external field is studied by Brownian Dynamics (BD) simulations and by an integral equation approach (IEA). Considering an identical nanoparticle as the source of the external field, we analyze the force exerted by the N confined particles on the external one, as well as the corresponding potential energy, focusing on their dependence on the distance of the external particle to the center of the sphere r₀. The connection of the force and potential to the equilibrium local distribution of the adsorbed particles, that is, the microstructure within the spherical monolayer induced by the external nanoparticle, which is also dependent on r₀, is elucidated by this analysis. It is found that the external particle needs to surmount a considerable potential barrier when moving toward the spherical surface, although much smaller than the one generated by a uniform surface distribution with an equivalent amount of charge. This is understood in terms of the correlation hole within the confined monolayer induced by the external particle. Another interesting conclusion is that the IEA provides an accurate, almost quantitative, description of the main features observed in the BD results, yet it is much less computationally demanding. The connection of these results with the overall chemical equilibrium of charged surfactant nanoparticles in the context of Pickering emulsions is also briefly discussed.

Active motion on curved surfaces

A theoretical analysis of active motion on curved surfaces is presented in terms of a generalization of the telegrapher equation. Such a generalized equation is explicitly derived as the polar approximation of the hierarchy of equations obtained from the corresponding Fokker-Planck equation of active particles diffusing on curved surfaces. The general solution to the generalized telegrapher equation is given for a pulse with vanishing current as initial data. Expressions for the probability density and the mean squared geodesic displacement are given in the limit of weak curvature. As an explicit example of the formulated theory, the case of active motion on the sphere is presented, where oscillations observed in the mean squared geodesic displacement are explained.

Active crystals on a sphere

Two-dimensional crystals on curved manifolds exhibit nontrivial defect structures. Here we consider "active crystals" on a sphere, which are composed of self-propelled colloidal particles. Our work is based on a phase-field-crystal-type model that involves a density and a polarization field on the sphere. Depending on the strength of the self-propulsion, three different types of crystals are found: a static crystal, a self-spinning "vortex-vortex" crystal containing two vortical poles of the local velocity, and a self-translating "source-sink" crystal with a source pole where crystallization occurs and a sink pole where the active crystal melts. These different crystalline states as well as their defects are studied theoretically here and can in principle be confirmed in experiments.