Source author record

Ignacio Pagonabarraga

Ignacio Pagonabarraga appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

39works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

39 published item(s)

preprint2022arXiv

A lattice Boltzmann model for self-diffusiophoretic particles near and at liquid-liquid interfaces

We introduce a novel mesoscopic computational model based on a multiphase-multicomponent lattice Boltzmann method for the simulation of self-phoretic particles in the presence of liquid-liquid interfaces. Our model features fully resolved solvent hydrodynamics and, thanks to its versatility, it can handle important aspects of the multiphysics of the problem, including particle wettability and differential solubility of the product in the two liquid phases. The method is extensively validated in simple numerical experiments, whose outcome is theoretically predictable, and then applied to the study of the behaviour of active particles next to and trapped at interfaces. We show that their motion can be variously steered by tuning relevant control parameters, such as the phoretic mobilities, the contact angle and the product solubility.

preprint2022arXiv

From Motility-Induced Phase-Separation to Glassiness in Dense Active Matter

Dense active systems are widespread in nature, examples range from bacterial colonies to biological tissues. Dense clusters of active particles can be obtained by increasing the packing fraction of the system or taking advantage of a peculiar phenomenon named motility-induced phase separation (MIPS). In this work, we explore the phase diagram of a two-dimensional model of active glass and show that disordered active materials develop a rich collective behaviour encompassing both MIPS and glassiness. We find that, although the glassy state is almost indistinguishable from that of equilibrium glasses, the mechanisms leading to its fluidization do not have any equilibrium counterpart. Our results can be rationalized in terms of a crossover between a low-activity regime, where glassy dynamics is controlled by an effective temperature, and a high-activity regime, which drives the system towards MIPS.

preprint2022arXiv

Impact of dipole-dipole interactions on motility-induced phase separation

We present a hydrodynamic theory for systems of dipolar active Brownian particles which, in the regime of weak dipolar coupling, predicts the onset of motility-induced phase separation (MIPS), consistent with Brownian dynamics (BD) simulations. The hydrodynamic equations are derived by explicitly coarse-graining the microscopic Langevin dynamics, thus allowing for a quantitative comparison of parameters entering the coarse-grained model and particle-resolved simulations. Performing BD simulations at fixed density, we find that dipolar interactions tend to hinder MIPS, as first reported in [Liao et al., Soft Matter, 2020, 16, 2208]. Here we demonstrate that the theoretical approach indeed captures the suppression of MIPS. Moreover, the analysis of the numerically obtained, angle-dependent correlation functions sheds light into the underlying microscopic mechanisms leading to the destabilization of the homogeneous phase.

preprint2022arXiv

Onsager reciprocal relations and chemo-mechanical coupling for chemically-active colloids

Similar to cells, bacteria, and other microorganisms, synthetic chemically-active colloids can harness the energy from their environment through a surface chemical reaction and use its energy to self-propel in fluidic environments. In this paper, we study the chemo-mechanical coupling that leads to the self-propulsion of chemically active colloids. The coupling between chemical reactions and momentum transport is a consequence of the Onsager reciprocal relations. They state that the velocity and the surface reaction rate are related to the mechanical and chemical affinities through a symmetric matrix. A consequence of the Onsager reciprocal relations is that, if a chemical reaction drives the motion of the colloid, then an external force generates a reaction rate. Here, we investigate the Onsager reciprocal relations for a spherical active colloid that catalyzes a reversible surface chemical reaction between two species. We solve the relevant transport equations using a perturbation expansion and numerical simulations to demonstrate the validity of the reciprocal relations around the equilibrium. Our results are consistent with previous studies and highlight the key role of solute advection in preserving the symmetry of the Onsager matrix. Finally, we show that the Onsager reciprocal relations break down around a nonequilibrium steady state, which has implications for the thermal fluctuations of the active colloids used in experiments.

preprint2021arXiv

Small Obstacle in a Large Polar Flock

We show that arbitrarily large polar flocks are susceptible to the presence of a single small obstacle. In a wide region of parameter space, the obstacle triggers counter-propagating dense bands leading to reversals of the flow. In very large systems, these bands interact yielding a never-ending chaotic dynamics that constitutes a new disordered phase of the system. While most of these results were obtained using simulations of aligning self-propelled particles, we find similar phenomena at the continuous level, not when considering the basic Toner-Tu hydrodynamic theory, but in simulations of truncations of the relevant Boltzmann equation.

preprint2020arXiv

Block copolymer-nanorod co-assembly in thin films: effects of rod-rod interaction and confinement

Simulations and experiments of nanorods (NRs) show that co-assembly with block copolymer (BCP) melts leads to the formation of a superstructure of side-to-side NRs perpendicular to the lamellar axis. A mesoscopic model is validated against scanning electron microscopy (SEM) images of CdSe NRs mixed with polystyrene-block-poly(methyl methacrylate). It is then used to study the co-assembly of anisotropic nanoparticles (NPs) with a length in the same order of magnitude as the lamellar spacing. The phase diagram of BCP/NP is explored as well as the time evolution of the NR. NRs that are slightly larger than the lamellar spacing are found to rotate and organise side-to-side with a tilted orientation with respect to the interface. Strongly interacting NPs are found to dominate the co-assembly while weakly interacting nanoparticles are less prone to form aggregates and tend to form well-ordered configurations.

preprint2020arXiv

Dynamics and Clogging of colloidal monolayers magnetically driven through an heterogeneous landscape

We combine experiments and numerical simulations to investigate the emergence of clogging in a system of interacting paramagnetic colloidal particles driven against a disordered landscape of larger obstacles. We consider a single aperture in a landscape of immobile silica particles which are irreversibly attached to the substrate. We use an external rotating magnetic field to generate a traveling wave potential which drives the magnetic particles against these obstacles at a constant and frequency tunable speed. Experimentally we find that the particles display an intermittent dynamics with power law distributions at high frequencies. We reproduce these results by using numerical simulations and shows that clogging in our system arises at large frequency, when the particles desynchronize with the moving landscape. Further, we use the model to explore the hidden role of flexibility in the obstacle displacements and the effect of hydrodynamic interactions between the particles. We also consider numerically the situation of a straight wall and investigate the range of parameters where clogging emerges in such case. Our work provides a soft matter test-bed system to investigate the effect of clogging in driven microscale matter.

preprint2020arXiv

Extraction of the propulsive speed of catalytic nano- and micro-motors under different motion dynamics

Motion of active particles, such as catalytic micro- and nano-motors, is usually characterized via either dynamic light scattering or optical microscopy. In both cases, speed of particles is obtained from the calculus of the mean square displacement (MSD) and typically, the theoretical formula of the MSD is derived from the motion equations of an active Brownian particle. One of the most commonly reported parameters is the speed of the particle, usually attributed to its propulsion, and is widely used to compare the motion efficiency of catalytic motors. However, it is common to find different methods to compute this parameter, which are not equivalent approximations and do not possess the same physical meaning. Here, we review the standard methods of speed analysis and focus on the errors that arise when analyzing the MSD of self-propelled particles. We analyze the errors from the computation of the instantaneous speed, as well as the propulsive speed and diffusion coefficient through fittings to parabolic equations, and we propose a revised formula for the motion analysis of catalytic particles moving with constant speed that can improve the accuracy and the amount of information obtained from the MSD. Moreover, we emphasize the importance of spotting the presence of different motion dynamics, such as particles with active angular speed or that move under the presence of drift, and how the breaking of ergodicity can completely change the analysis by considering particles with an exponentially decaying speed. In all cases, real data from enzymatically propelled micro-motors and simulations are used to back up the theories. Finally, we propose several analytical approaches and analyze limiting cases that will help to deal with these scenarios while still obtaining accurate results.

preprint2020arXiv

Fluctuation-Dissipation Relations in the absence of Detailed Balance: formalism and applications to Active Matter

We present a comprehensive study about the relationship between the way Detailed Balance is broken in non-equilibrium systems and the resulting violations of the Fluctuation-Dissipation Theorem. Starting from stochastic dynamics with both odd and even variables under Time-Reversal, we exploit the relation between entropy production and the breakdown of Detailed Balance to establish general constraints on the non-equilibrium steady-states (NESS), which relate the non-equilibrium character of the dynamics with symmetry properties of the NESS distribution. This provides a direct route to derive extended Fluctuation-Dissipation Relations, expressing the linear response function in terms of NESS correlations. Such framework provides a unified way to understand the departure from equilibrium of active systems and its linear response. We then consider two paradigmatic models of interacting self-propelled particles, namely Active Brownian Particles (ABP) and Active Ornstein-Uhlenbeck Particles (AOUP). We analyze the non-equilibrium character of these systems (also within a Markov and a Chapman-Enskog approximation) and derive extended Fluctuation-Dissipation Relations for them, clarifying which features of these active model systems are genuinely non-equilibrium.

preprint2020arXiv

Propulsion and energetics of a minimal magnetic microswimmer

In this manuscript we describe the realization of a minimal hybrid microswimmer, composed of a ferromagnetic nanorod and a paramagnetic microsphere. The unbounded pair is propelled in water upon application of a swinging magnetic field that induces a periodic relative movement of the two composing elements, where the nanorod rotates and slides on the surface of the paramagnetic sphere. When taken together, the processes of rotation and sliding describe a finite area in the parameter space, which increases with the frequency of the applied field. We develop a theoretical approach and combine it with numerical simulations, which allow us to understand the dynamics of the propeller and explain the experimental observations. Furthermore, we demonstrate a reversal of the microswimmer velocity by varying the length of the nanorod, as predicted by the model. Finally, we determine theoretically and in experiments the Lighthill's energetic efficiency of this minimal magnetic microswimmer.

preprint2020arXiv

Spontaneous polarization and locomotion of an active particle with surface-mobile enzymes

We examine a mechanism of locomotion of active particles whose surface is uniformly coated with mobile enzymes. The enzymes catalyze a reaction that drives phoretic flows but their homogeneous distribution forbids locomotion by symmetry. We find that the ability of the enzymes to migrate over the surface combined with self-phoresis can lead to a spontaneous symmetry breaking instability whereby the homogeneous distribution of enzymes polarizes and the particle propels. The instability is driven by the advection of enzymes by the phoretic flows and occurs above a critical Péclet number. The transition to polarized motile states occurs via a supercritical or subcritical pitchfork bifurcations, the latter of which enables hysteresis and coexistence of uniform and polarized states.

preprint2019arXiv

Asymmetric and Long Range Interactions in Shaken Granular Media

We use a computational model to investigate the emergence of interaction forces between pairs of intruders in a horizontally vibrated granular fluid. The time evolution of a pair of particles shows a maximum of the likelihood to find the pair at contact in the direction of shaking. This relative interaction is further studied by fixing the intruders in the simulation box where we identify effective mechanical forces, and torques between particles and quantify an emergent long range attractive force as a function of the shaking relative angle, amplitude, and the packing density of grains. We determine the local density and kinetic energy profiles of granular particles along the axis of the dimer to find no gradients in the density fields and additive gradients in the kinetic energies.

preprint2019arXiv

Linear response theory and Green-Kubo relations for active matter

We address the question of how interacting active systems in a non-equilibrium steady-state respond to an external perturbation. We establish an extended fluctuation-dissipation theorem for Active Brownian Particles (ABP) which highlights the role played by the local violation of detailed balance due to activity. By making use of a Markovian approximation we derive closed Green-Kubo expressions for the diffusivity and mobility of ABP and quantify the deviations from the Stokes-Einstein relation. We compute the linear response function to an external force using unperturbed simulations of ABP and compare the results with the analytical predictions of the transport coefficients. Our results show the importance of the interplay between activity and interactions in the departure from equilibrium linear response.

preprint2016arXiv

Colloidal Microworms Propelling via a Cooperative Hydrodynamic Conveyor Belt

We study propulsion arising from microscopic colloidal rotors dynamically assembled and driven in a viscous fluid upon application of an elliptically polarized rotating magnetic field. Close to a confining plate, the motion of this self-assembled microscopic worm results from the cooperative flow generated by the spinning particles which act as a hydrodynamic "conveyor belt." Chains of rotors propel faster than individual ones, until reaching a saturation speed at distances where induced-flow additivity vanishes. By combining experiments and theoretical arguments, we elucidate the mechanism of motion and fully characterize the propulsion speed in terms of the field parameters.

preprint2016arXiv

Emergent structures and dynamics of cell colonies by contact inhibition of locomotion

Cells in tissues can organize into a broad spectrum of structures according to their function. Drastic changes of organization, such as epithelial-mesenchymal transitions or the formation of spheroidal aggregates, are often associated either to tissue morphogenesis or to cancer progression. Here, we study the organization of cell colonies by means of simulations of self-propelled particles with generic cell-like interactions. The interplay between cell softness, cell-cell adhesion, and contact inhibition of locomotion (CIL) yields structures and collective dynamics observed in several existing tissue phenotypes. These include regular distributions of cells, dynamic cell clusters, gel-like networks, collectively migrating monolayers, and 3D aggregates. We give analytical predictions for transitions between noncohesive, cohesive, and 3D cell arrangements. We explicitly show how CIL yields an effective repulsion that promotes cell dispersal, thereby hindering the formation of cohesive tissues. Yet, in continuous monolayers, CIL leads to collective cell motion, ensures tensile intercellular stresses, and opposes cell extrusion. Thus, our work highlights the prominent role of CIL in determining the emergent structures and dynamics of cell colonies.

preprint2016arXiv

Entropically-induced asymmetric passage times of charged tracers across corrugated channels

We analyze the diffusion of charged and neutral tracers suspended in an electrolyte embedded in a channel of varying cross-section. Making use of systematic approximations, the diffusion equation governing the motion of tracers is mapped into an effective $1D$ equation describing the dynamics along the longitudinal axis of the channel where its varying-section is encoded as an effective entropic potential. This simplified approach allows us to characterize tracer diffusion under generic confinement by measuring their mean first passage time (MFPT). In particular, we show that the interplay between geometrical confinement and electrostatic interactions strongly affect the MFTP of tracers across corrugated channels hence leading to alternative means to control tracers translocation across charged pores. Finally, our results show that the MFPTs of a charged tracer in opposite directions along an asymmetric channel may differ. We expect our results to be relevant for biological as well synthetic devices whose dynamics is controlled by the detection of diluted tracers.

preprint2016arXiv

Morphology of clusters of attractive dry and wet self-propelled spherical particle suspensions

In order to asses the effect of hydrodynamics in the assembly of active attractive spheres, we simulate a semi-dilute suspension of attractive self-propelled spherical particles in a quasi two dimensional geometry comparing the case with and without hydrodynamics interactions. To start with, independently on the presence of hydrodynamics, we observe that depending on the ratio between attraction and propulsion, particles either coarsen or aggregate forming finite-size clusters. Focusing on the clustering regime, we characterize two different clusters parameters, i.e. their morphology and orientational order, and compare the case when active particles behave either as pushers or pullers (always in the regime where inter-particles attractions competes with self-propulsion). Studying cluster phases for squirmers with respect to those obtained for active Brownian disks (indicated as ABP), we have shown that hydrodynamics alone can sustain a cluster phase of active swimmers (pullers), while ABP form cluster phases due to the competition between attraction and self propulsion. The structural properties of the cluster phases of squirmers and ABP are similar, although squirmers show sensitivity to active stresses. Active Brownian disks resemble weakly pusher squirmer suspensions in terms of cluster size distribution, structure of the radius of gyration on cluster size and degree of cluster polarity.

preprint2016arXiv

Programmable Active Janus Droplets Driven by Water/Alcohol Phase Separation

We report the existence of self-propelled Janus droplets driven by phase separation, which are able to deliver cargo in a programmable manner. The self-propelling droplets are initially formed by a water/ethanol mixture in a squalane/monoolein solution, and evolve in up to three stages depending on ethanol concentration. In the first stage, the droplet propulsion is generated by Marangoni flow originating from the solubilization of ethanol in the oily phase. During this process the droplets absorb surfactant molecules; in combination with the continuous loss of ethanol this leads to a phase separation of the water/ethanol/monoolein mixture and the formation of Janus droplets, i.e. a water-rich droplet connected to an ethanol-rich droplet that is able to deliver cargo. We characterize the different evolution stages of self-propulsion by the flow field around the droplet that evolves from a weak pusher, over a neutral swimmer, to a dimer of neutral swimmers. Finally, we utilize this active system to deliver DNA as a cargo. Tuning the delay time before phase separation, by varying the chemical composition of the droplets, several different cargo delivery processes can be programmed.

preprint2016arXiv

Spontaneous aggregation and global polar ordering in squirmer suspensions

We have developed numerical simulations of three dimensional suspensions of active particles to characterize the capabilities of the hydrodynamic stresses induced by active swimmers to promote global order and emergent structures in active suspensions. We have considered squirmer suspensions embedded in a fluid modeled under a Lattice Boltzmann scheme. We have found that active stresses play a central role to decorrelate the collective motion of squirmers and that contractile squirmers develop significant aggregates.

preprint2015arXiv

Clustering and pattern formation in chemorepulsive active colloids

We demonstrate that migration away from self-produced chemicals (chemorepulsion) generates a generic route to clustering and pattern formation among self-propelled colloids. The clustering instability can be caused either by anisotropic chemical production, or by a delayed orientational response to changes of the chemical environment. In each case, chemorepulsion creates clusters of a self-limiting area which grows linearly with self-propulsion speed. This agrees with recent observations of dynamic clusters in Janus colloids (albeit not yet known to be chemorepulsive). More generally, our results could inform design principles for the self-assembly of chemorepulsive synthetic swimmers and/or bacteria into nonequilibrium patterns.

preprint2015arXiv

Geometrically-tuned channel permeability

We characterize the motion of charged as well as neutral tracers, in an electrolyte embedded in a varying section channel. We exploit a set of systematic approximations that allows us to simplify the problem, yet capturing the essential of the interplay between the geometrical confinement provided by the corrugated channel walls and the electrolyte properties. Our simplified approach allows us to characterize the transport properties of corrugated channels when a net flux of tracers is obtained by keeping the extrema of the channel at different chemical potentials. For highly diluted tracer suspensions, we have characterized tracers currents and we have estimated the net electric current which occurs when both positively and negatively charged tracers are considered.

preprint2015arXiv

Tracer diffusion of hard-sphere binary mixtures under nano-confinement

The physics of diffusion phenomena in nano and micro channels has attracted a lot of attention in recent years, due to its close connection with many technological, medical and industrial applications. In the present paper we employ a kinetic approach to investigate how the confinement in nanostructured geometries affects the diffusive properties of fluid mixtures and leads to the appearance of properties different from those of bulk systems. In particular, we derive an expression for the friction tensor in the case of a bulk fluid mixture confined to a narrow slit having undulated walls. The boundary roughness leads to a new mechanism for transverse diffusion, and can even lead to an effective diffusion along the channel larger than the one corresponding to a planar channel of equivalent section. Finally we discuss a reduction of the previous equation to a one dimensional effective diffusion equation in which an entropic term encapsulates the geometrical information on the channel shape.

preprint2014arXiv

Accounting for adsorption and desorption in lattice Boltzmann simulations

We report a Lattice-Boltzmann scheme that accounts for adsorption and desorption in the calculation of mesoscale dynamical properties of tracers in media of arbitrary complexity. Lattice Boltzmann simulations made it possible to solve numerically the coupled Navier-Stokes equations of fluid dynamics and Nernst-Planck equations of electrokinetics in complex, heterogeneous media. Associated to the moment propagation scheme, it became possible to extract the effective diffusion and dispersion coefficients of tracers, or solutes, of any charge, e.g. in porous media. Nevertheless, the dynamical properties of tracers depend on the tracer-surface affinity, which is not purely electrostatic, but also includes a species-specific contribution. In order to capture this important feature, we introduce specific adsorption and desorption processes in a Lattice-Boltzmann scheme through a modified moment propagation algorithm, in which tracers may adsorb and desorb from surfaces through kinetic reaction rates. The method is validated on exact results for pure diffusion and diffusion-advection in Poiseuille flows in a simple geometry. We finally illustrate the importance of taking such processes into account on the time-dependent diffusion coefficient in a more complex porous medium.

preprint2014arXiv

Entropic electrokinetics: recirculation, particle separation and negative mobility

We show that when particles are suspended in an electrolyte confined between corrugated charged surfaces, electrokinetic flows lead to a new set of phenomena such as particle separation, mixing for low-Reynolds micro- and nano-metric devices and negative mobility. Our analysis shows that such phenomena arise, for incompressible fluids, due to the interplay between the electrostatic double layer and the corrugated geometrical confinement and that they are magnified when the width of the channel is comparable to the Debye length. Our characterization allows us to understand the physical origin of such phenomena therefore shading light on their possible relevance in a wide variety of situations, ranging from nano- and micro-fluidic devices to biological systems.

preprint2014arXiv

On the role of composition entropies in the statistical mechanics of polydisperse systems

Polydisperse systems are commonly encountered when dealing with soft matter in general or any non simple fluid. Yet their treatment within the framework of statistical thermodynamics is a delicate task as the latter has been essentially devised for simple --- non fully polydisperse --- systems. In this paper, we address the issue of defining a non ambiguous combinatorial entropy for these systems. We do so by focusing on the general property of extensivity of the thermodynamic potentials and discussing a specific mixing experiment. This leads us to introduce the new concept of composition entropy for single phase systems that we do not assimilate to a mixing entropy. We then show that they do not contribute to the thermodynamics of the system at fixed composition and prescribe to substract $\ln N!$ from the free energy characterizing a system however polydisperse it can be. We then re-derive general expressions for the mixing entropy between any two polydisperse systems and interpret them in term of distances between probability distributions and show that one of these metrics relates naturally to a recent extension of Landauer's principle. We then propose limiting expressions for the mixing entropy in the case of mixing with equal proportion in the original compositions and finally address the challenging problem of chemical reactions.

preprint2014arXiv

Working under confinement

We analyze the performance of a Brownian ratchet in the presence of geometrical constraints. A two-state model that describes the kinetics of molecular motors is used to characterize the energetic cost when the motor proceeds under confinement, in the presence of an external force. We show that the presence of geometrical constraints has a strong effect on the performance of the motor. In particular, we show that it is possible to enhance the ratchet performance by a proper tuning of the parameters characterizing the environment. These results open the possibility of engineering entropically-optimized transport devices.

preprint2013arXiv

Confined Brownian ratchets

We analyze the dynamics of Brownian ratchets in a confined environment. The motion of the particles is described by a Fick-Jakobs kinetic equation in which the presence of boundaries is modeled by means of an entropic potential. The cases of a flashing ratchet, a two-state model and a ratchet under the influence of a temperature gradient are analyzed in detail. We show the emergence of a strong cooperativity between the inherent rectification of the ratchet mechanism and the entropic bias of the fluctuations caused by spatial confinement. Net particle transport may take place in situations where none of those mechanisms leads to rectification when acting individually. The combined rectification mechanisms may lead to bidirectional transport and to new routes to segregation phenomena. Confined Brownian ratchets (CBR) could be used to control transport in mesostructures and to engineer new and more efficient devices for transport at the nanoscale.

preprint2013arXiv

Theory of wetting-induced fluid entrainment by advancing contact lines on dry surfaces

We report on the onset of fluid entrainment when a contact line is forced to advance over a dry solid of arbitrary wettability. We show that entrainment occurs at a critical advancing speed beyond which the balance between capillary, viscous and contact line forces sustaining the shape of the interface is no longer satisfied. Wetting couples to the hydrodynamics by setting both the morphology of the interface at small scales and the viscous friction of the front. We find that the critical deformation that the interface can sustain is controlled by the friction at the contact line and the viscosity contrast between the displacing and displaced fluids, leading to a rich variety of wetting-entrainment regimes. We discuss the potential use of our theory to measure contact-line forces using atomic force microscopy, and to study entrainment under microfluidic conditions exploiting colloid-polymer fluids of ultra-low surface tension.

preprint2012arXiv

Cooperative rectification in confined Brownian ratchets

We analyze the rectified motion of a Brownian particle in a confined environment. We show the emergence of strong cooperativity between the inherent rectification of the ratchet mechanism and the entropic bias of the fluctuations caused by spatial confinement. Net particle transport may develop even in situations where separately the ratchet and the geometric restrictions do not give rise to particle motion. The combined rectification effects can lead to bidirectional transport depending on particle size, resulting in a new route for segregation. The reported mechanism can be used to control transport in mesostructures and nanodevices in which particles move in a reduced space.

preprint2012arXiv

Running faster together: huge speed up of thermal ratchets due to hydrodynamic coupling

We present simulations that reveal a surprisingly large effect of hydrodynamic coupling on the speed of thermal ratchet motors. The model that we use considers particles performing thermal ratchet motion in a hydrodynamic solvent. Using particle-based, mesoscopic simulations that maintain local momentum conservation, we analyze quantitatively how the coupling to the surrounding fluid affects ratchet motion. We find that coupling can increase the mean velocity of the moving particles by almost two orders of magnitude, precisely because ratchet motion has both a diffusive and a deterministic component. The resulting coupling also leads to the formation of aggregates at longer times. The correlated motion that we describe increases the efficiency of motor-delivered cargo transport and we speculate that the mechanism that we have uncovered may play a key role in speeding up molecular motor-driven intracellular transport.

preprint2011arXiv

Colloids electrophoresis for strong and weak ion diffusivity

We study the electrophoretic flow of suspensions of charged colloids with a mesoscopic method that allows to model generic experimental conditions. We show that for highly charged colloids their electrophoretic mobility increases significantly and displays a mobility maximum on increasing the colloidal charge for all salt concentrations. The electrophoretic mobility of highly charged colloids is also enhanced significantly when ion advection is dominant, leading to a strong heterogeneity in the local electrophoretic response especially at low salt concentration, when ion diffuse layers overlap.

preprint2011arXiv

Determination of the zeta potential for highly charged colloidal suspensions

We compute the electrostatic potential at the surface, or zeta potential $ζ$, of a charged particle embedded in a colloidal suspension using a hybrid mesoscopic model. We show that for weakly perturbing electric fields, the value of $ζ$ obtained at steady state during electrophoresis is statistically indistinguishable from $ζ$ in thermodynamic equilibrium. We quantify the effect of counterions concentration on $ζ$. We also evaluate the relevance of the lattice resolution for the calculation of $ζ$ and discuss how to identify the effective electrostatic radius.

preprint2011arXiv

Growth saturation of unstable thin films on transverse-striped hydrophilic-hydrophobic micropatterns

Using three-dimensional numerical simulations, we demonstrate the growth saturation of an unstable thin liquid film on micropatterned hydrophilic-hydrophobic substrates. We consider different transverse-striped micropatterns, characterized by the total fraction of hydrophilic coverage and the width of the hydrophilic stripes. We compare the growth of the film on the micropatterns to the steady states observed on homogeneous substrates, which correspond to a saturated sawtooth and growing finger configurations for hydrophilic and hydrophobic substrates, respectively. The proposed micropatterns trigger an alternating fingering-spreading dynamics of the film, which leads to a complete suppression of the contact line growth above a critical fraction of hydrophilic stripes. Furthermore, we find that increasing the width of the hydrophilic stripes slows down the advancing front, giving smaller critical fractions the wider the hydrophilic stripes are. Using analytical approximations, we quantitatively predict the growth rate of the contact line as a function of the covering fraction, and predict the threshold fraction for saturation as a function of the stripe width.

preprint2011arXiv

Instrinsic oscillations of treadmilling microtubules in a motor bath

We analyse the dynamics of overlapping antiparallel treadmilling microtubules in the presence of crosslinking processive motor proteins that counterbalance an external force. We show that coupling the force-dependent velocity of motors and the kinetics of motor exchange with a bath in the presence of treadmilling leads generically to oscillatory behavior. In addition we show that coupling the polymerization kinetics to the external force through the kinetics of the crosslinking motors can stabilize the oscillatory instability into finite-amplitude nonlinear oscillations and may lead to other scenarios, including bistability.

preprint2011arXiv

Lattice Boltzmann - Langevin simulations of binary mixtures

We report a hybrid numerical method for the solution of the model H fluctuating hydrodynamic equations for binary mixtures. The momentum conservation equations with Landau-Lifshitz stresses are solved using the fluctuating lattice Boltzmann equation while the order parameter conservation equation with Langevin fluxes are solved using the stochastic method of lines. Two methods, based on finite difference and finite volume, are proposed for spatial discretisation of the order parameter equation. Special care is taken to ensure that the fluctuation-dissipation theorem is maintained at the lattice level in both cases. The methods are benchmarked by comparing static and dynamic correlations and excellent agreement is found between analytical and numerical results. The Galilean invariance of the model is tested and found to be satisfactory. Thermally induced capillary fluctuations of the interface are captured accurately, indicating that the model can be used to study nonlinear fluctuations.

preprint2011arXiv

Phase segregation and transport in a two species multi-lane system

We present a two channel driven lattice gas model with oppositely directed species moving on two parallel lanes with lane switching processes. We study correlated lane switching mechanism for particles so that switching may occur with finite probability only when oppositely directed species meet on the same channel. The system is analyzed for closed ring with conserved total particle number. For asymmetric particle exchange between the lanes, the system exhibits unique polarization phenomenon with segregation of oppositely directed species between the two lanes. The polarization phenomenon can be understood as a consequence of existence of an absorbing steady state. For symmetric exchange rate of particles between the lanes, the system remains unpolarized, with equal particle density on both the lanes in the thermodynamic limit of large system size. We study the system using a combination of a Mean Field(MF) analysis and Monte Carlo simulations. The nature of phase segregation that we see for this system, is distinct from driven particle systems which are in contact with particle reservoir. The features observed for this minimal model will have ramifications for biofilament based intracellular transport, wherein cellular cargo, e.g; organelles and vesicles are transported by oppositely directed particles on multiple filament tracks.

preprint2010arXiv

Dynamics of gravity driven three-dimensional thin films on hydrophilic-hydrophobic patterned substrates

We investigate numerically the dynamics of unstable gravity driven three-dimensional thin liquid films on hydrophilic-hydrophobic patterned substrates of longitudinal stripes and checkerboard arrangements. The thin film can be guided preferentially on hydrophilic longitudinal stripes, while fingers develop on adjacent hydrophobic stripes if their width is large enough. On checkerboard patterns, the film fingering occurs on hydrophobic domains, while lateral spreading is favoured on hydrophilic domains, providing a mechanism to tune the growth rate of the film. By means of kinematical arguments, we quantitatively predict the growth rate of the contact line on checkerboard arrangements, providing a first step towards potential techniques that control thin film growth in experimental setups.

preprint2010arXiv

Lattice-gas model for active vesicle transport by molecular motors with opposite polarities

We introduce a multi-species lattice gas model for motor protein driven collective cargo transport on cellular filaments. We use this model to describe and analyze the collective motion of interacting vesicle cargoes being carried by oppositely directed molecular motors, moving on a single biofilament. Building on a totally asymmetric exclusion process (TASEP) to characterize the motion of the interacting cargoes, we allow for mass exchange with the environment, input and output at filament boundaries and focus on the role of interconversion rates and how they affect the directionality of the net cargo transport. We quantify the effect of the various different competing processes in terms of non-equilibrium phase diagrams. The interplay of interconversion rates, which allow for flux reversal and evaporation/deposition processes introduce qualitatively new features in the phase diagrams. We observe regimes of three-phase coexistence, the possibility of phase re-entrance and a significant flexibility in how the different phase boundaries shift in response to changes in control parameters. The moving steady state solutions of this model allows for different possibilities for the spatial distribution of cargo vesicles, ranging from homogeneous distribution of vesicles to polarized distributions, characterized by inhomogeneities or {\it shocks}. Current reversals due to internal regulation emerge naturally within the framework of this model. We believe this minimal model will clarify the understanding of many features of collective vesicle transport, apart from serving as the basis for building more exact quantitative models for vesicle transport relevant to various {\it in-vivo} situations.

preprint2003arXiv

Kinetics of Fragmenting Freely Evolving Granular Gases

We investigate the effect of fragmentation on the homogeneous free cooling of inelastic hard spheres, using Boltzmann kinetic theory and Direct Monte Carlo simulations. We analyze in detail a model where dissipative collisions may subsequently lead to a break-up of the grains. With a given probability, two off-springs are then created from one of the two colliding partners, with conservation of mass, momentum and kinetic energy. We observe a scaling regime characterized by a single collisional average, that quantifies the deviations from Gaussian behaviour for the joint size and velocity distribution function. We also discuss the possibility of a catastrophe whereby the number of particles diverges in a finite time. This phenomenon appears correlated to a ``shattering'' transition marked by a delta singularity at vanishingly small grains for the rescaled size distribution.