BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Holger Stark

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

ResearcherUnclaimed source record
cond-mat.softphysics.flu-dynBiological Physicscond-mat.stat-mechCell Behaviorcond-mat.mes-hallphysics.chem-phphysics.class-phphysics.opticsPopulations and Evolution

Catalog footprint

What is connected

39works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

Characterizing elastic turbulence in the three-dimensional von Karman swirling flow using the Oldroyd-B model

We present the full three-dimensional numerical investigation of the von Karman swirling flow between two parallel plates using the Oldroyd-B model and characterize the onset and development of elastic turbulence. We quantify the flow state with the secondary-flow strength, a measure of the average strength of the velocity fluctuations, and then define an order parameter as the time average of the secondary-flow strength. The order parameter displays a subcritical transition from the laminar to a bistable flow that switches between weakly chaotic flow and elastic turbulence. The transition to the bistable flow occurs at the critical Weissenberg number $\mathrm{Wi_c}=12$. Above $\mathrm{Wi_c}$, in the elastic turbulent state, we observe a strong increase in velocity fluctuations and flow resistance, which we define as the total work performed on the fluid. Upon starting simulations in the turbulent state and subsequently lowering $\mathrm{Wi}$ below its critical value, we observe hysteretic behavior in the order parameter and the flow resistance, which is a common feature of a subcritical transition. Hysteresis has also been found in experiments. Additionally, we find power-law scaling in the spatial and temporal power spectra of the velocity fluctuations, characteristic for elastic turbulence. The maximum values of the power-law exponents in our simulations are $α_t= 3.41$ for the temporal exponent and $α_s= 3.17$ for the spatial exponent, which are remarkably close to the values obtained in experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Gyrotactic cluster formation of bottom-heavy squirmers

Squirmers that are bottom-heavy experience a torque that aligns them along the vertical so that they swim upwards. In a suspension of many squirmers, they also interact hydrodynamically via flow fields that are initiated by their swimming motion and by gravity. Swimming under the combined action of flow field vorticity and gravitational torque is called gyrotaxis. Using the method of multi-particle collision dynamics, we perform hydrodynamic simulations of a many-squirmer system floating above the bottom surface. Due to gyrotaxis they exhibit pronounced cluster formation with increasing gravitational torque. The clusters are more volatile at low values but compactify to smaller clusters at larger torques. The mean distance between clusters is mainly controlled by the gravitational torque and not the global density. Furthermore, we observe that neutral squirmers form clusters more easily, whereas pullers require larger gravitational torques due to their additional force-dipole flow fields. We do not observe clustering for pusher squirmers. Adding a rotlet dipole to the squirmer flow field induces swirling clusters. At high gravitational strengths, the hydrodynamic interactions with the no-slip boundary create an additional vertical alignment for neutral squirmers, which also supports cluster formation.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Turbulence-induced clustering in compressible active fluids

We study a novel phase of active polar fluids, which is characterized by the continuous creation and destruction of dense clusters due to self-sustained turbulence. This state arises due to the interplay of the self-advection of the aligned swimmers and their defect topology. The typical cluster size is determined by the characteristic vortex size. Our results are obtained by investigating a continuum model of compressible polar active fluids, which incorporates typical experimental observations in bacterial suspensions by assuming a non-monotone dependence of speed on density.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Emergent collective dynamics of bottom-heavy squirmers under gravity

We present the results of hydrodynamic simulations using the method of multi-particle collision dynamics for a system of squirmer microswimmers moving under the influence of gravity at low Reynolds numbers. In addition, the squirmers are bottom-heavy so that they experience a torque which aligns them along the vertical. The squirmers interact hydrodynamically by the flow fields of a stokeslet and rotlet, which are initiated by the acting gravitational force and torque, respectively, and by their own flow fields. By varying the ratio of swimming to bulk sedimentation velocity and the torque, we determine state diagrams for the emergent collective dynamics of neutral squirmers as well as strong pushers and pullers. For low swimming velocity and torque we observe conventional sedimentation, while the sedimentation profile becomes inverted when their values are increased. For neutral squirmers we discover convective rolls of circulating squirmers between both sedimentation states, which sit at the bottom of the system and are fed by plumes made of collectively sinking squirmers. At larger torques porous clusters occur that spawn single squirmers. The two latter states can also occur transiently starting from a uniform squirmer distribution and then disappear in the long-time limit. For strong pushers and pullers only weak plume formation is observed.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Optimal control of colloidal trajectories in inertial microfluidics using the Saffman effect

In inertial microfluidics colloidal particles in a Poiseuille flow experience the Segré-Silberberg lift force, which drives them to specific positions in the channel cross section. Due to the Saffman effect an external force applied along the microchannel induces a cross-streamline migration to a new equilibrium position. We apply optimal control theory to design the time protocol of the axial control force in order to steer a single particle as precisely as possible from a channel inlet to an outlet at a chosen target position. We discuss the influence of particle radius and channel length and show that optimal steering is cheaper than using a constant control force. Using a single optimized control-force protocol, we demonstrate that even a pulse of particles spread along the channel axis can be steered to a target and that particles of different radii can be separarted most efficiently.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Particle pairs and trains in inertial microfluidics

Staggered and linear multi-particle trains constitute characteristic structures in inertial microfluidics. Using lattice-Boltzmann simulations, we investigate their properties and stability, when flowing through microfluidic channels. We confirm the stability of cross-streamline pairs by showing how they contract or expand to their equilibrium axial distance. In contrast, same-streamline pairs quickly expand to a characteristic separation but even at long times slowly drift apart. We reproduce the distribution of particle distances with its characteristic peak as measured in experiments. Staggered multi-particle trains initialized with an axial particle spacing larger than the equilibrium distance contract non-uniformly due to collective drag reduction. Linear particle trains, similar to pairs, rapidly expand towards a value about twice the equilibrium distance of staggered trains and then very slowly drift apart non-uniformly. Again, we reproduce the statistics of particle distances and the characteristic peak observed in experiments. Finally, we thoroughly analyze the damped displacement pulse traveling as a microfluidic phonon through a staggered train and show how a defect strongly damps its propagation.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Optimal steering of a smart active particle

We formulate the theory for steering an active particle with optimal travel time between two locations and apply it to the Mexican hat potential without brim. For small heights the particle can cross the potential barrier, while for large heights it has to move around it. Thermal fluctuations in the orientation strongly affect the path over the barrier. Then we consider a smart active particle and apply reinforcement learning. We show how the active particle learns in repeating episodes to move optimally. The optimal steering is stored in the optimized action-value function, which is able to rectify thermal fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Active Brownian motion of emulsion droplets: Coarsening dynamics at the interface and rotational diffusion

A micron-sized droplet of bromine water immersed in a surfactant-laden oil phase can swim (S. Thutupalli, R. Seemann, S. Herminghaus, New J. Phys. 13 073021 (2011)). The bromine reacts with the surfactant at the droplet interface and generates a surfactant mixture. It can spontaneously phase-separate due to solutocapillary Marangoni flow, which propels the droplet. We model the system by a diffusion-advection-reaction equation for the mixture order parameter at the interface including thermal noise and couple it to fluid flow. Going beyond previous work, we illustrate the coarsening dynamics of the surfactant mixture towards phase separation in the axisymmetric swimming state. Coarsening proceeds in two steps: an initially slow growth of domain size followed by a nearly ballistic regime. On larger time scales thermal fluctuations in the local surfactant composition initiates random changes in the swimming direction and the droplet performs a persistent random walk, as observed in experiments. Numerical solutions show that the rotational correlation time scales with the square of the inverse noise strength. We confirm this scaling by a perturbation theory for the fluctuations in the mixture order parameter and thereby identify the active emulsion droplet as an active Brownian particle.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Active Brownian particles and run-and-tumble particles separate inside a maze

A diverse range of natural and artificial self-propelled particles are known and are used nowadays. Among them, active Brownian particles (ABPs) and run-and-tumble particles (RTPs) are two important classes. We numerically study non-interacting ABPs and RTPs strongly confined to different maze geometries in two dimensions. We demonstrate that by means of geometrical confinement alone, ABPs are separable from RTPs. By investigating Matryoshka-like mazes with nested shells, we show that a circular maze has the best filtration efficiency. Results on the mean first-passage time reveal that ABPs escape faster from the center of the maze, while RTPs reach the center from the rim more easily. According to our simulations and a rate theory, which we developed, ABPs in steady state accumulate in the outermost region of the Matryoshka-like mazes, while RTPs occupy all locations within the maze with nearly equal probability. These results suggest a novel technique for separating different types of self-propelled particles by designing appropriate confining geometries without using chemical or biological agents.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Active Brownian particles moving in a random Lorentz gas

Biological microswimmers often inhabit a porous or crowded environment such as soil. In order to understand how such a complex environment influences their spreading, we numerically study non-interacting active Brownian particles (ABPs) in a two-dimensional random Lorentz gas. Close to the percolation transition in the Lorentz gas, they perform the same subdiffusive motion as ballistic and diffusive particles. However, due to their persistent motion they reach their long-time dynamics faster than passive particles and also show superdiffusive motion at intermediate times. While above the critical obstacle density $η_c$ the ABPs are trapped, their long-time diffusion below $η_c$ is strongly influenced by the propulsion speed $v_0$. With increasing $v_0$, ABPs are stuck at the obstacles for longer times. Thus, for large propulsion speed, the long-time diffusion constant decreases more strongly in a denser obstacle environment than for passive particles. This agrees with the behavior of an effective swimming velocity and persistence time, which we extract from the velocity autocorrelation function.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dynamics of a bacterial flagellum under reverse rotation

To initiate tumbling of an E. coli, one of the helical flagella reverses its sense of rotation. It then transforms from its normal form first to the transient semicoiled state and subsequently to the curly-I state. The dynamics of polymorphism is effectively modeled by describing flagellar elasticity through an extended Kirchhoff free energy. However, the complete landscape of the free energy remains undetermined because the ground state energies of the polymorphic forms are not known. We investigate how variations in these ground state energies affect the dynamics of a reversely rotated flagellum of a swimming bacterium. We find that the flagellum exhibits a number of distinct dynamical states and comprehensively summarize them in a state diagram. As a result, we conclude that tuning the landscape of the extended Kirchhoff free energy alone cannot generate the intermediate full-length semicoiled state. However, our model suggests an ad hoc method to realize the sequence of polymorphic states as observed for a real bacterium. Since the elastic properties of bacterial flagella are similar, our findings can easily be extended to other peritrichous bacteria.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Emergent behavior in active colloids

Active colloids are microscopic particles, which self-propel through viscous fluids by converting energy extracted from their environment into directed motion. We first explain how articial microswimmers move forward by generating near-surface flow fields via self-phoresis or the self-induced Marangoni effect. We then discuss generic features of the dynamics of single active colloids in bulk and in confinement, as well as in the presence of gravity, field gradients, and fluid flow. In the third part, we review the emergent collective behavior of active colloidal suspensions focussing on their structural and dynamic properties. After summarizing experimental observations, we give an overview on the progress in modeling collectively moving active colloids. While active Brownian particles are heavily used to study collective dynamics on large scales, more advanced methods are necessary to explore the importance of hydrodynamic and phoretic particle interactions. Finally, the relevant physical approaches to quantify the emergent collective behavior are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Kinetic Monte Carlo Simulations for Birefringence Relaxation of Photo-Switchable Molecules on a Surface

Recent experiments have demonstrated that in a dense monolayer of photo-switchable dye Methyl-Red molecules the relaxation of an initial birefringence follows a power-law decay, typical for glass-like dynamics. The slow relaxation can efficiently be controlled and accelerated by illuminating the monolayer with circularly polarized light, which induces trans-cis isomerization cycles. To elucidate the microscopic mechanism, we develop a two-dimensional molecular model in which the trans and cis isomers are represented by straight and bent needles, respectively. As in the experimental system, the needles are allowed to rotate and to form overlaps but they cannot translate. The out-of-equilibrium rotational dynamics of the needles is generated using kinetic Monte Carlo simulations. We demonstrate that, in a regime of high density and low temperature, the power-law relaxation can be traced to the formation of spatio-temporal correlations in the rotational dy- namics, i.e., dynamic heterogeneity. We also show that the nearly isotropic cis isomers can prevent dynamic heterogeneity from forming in the monolayer and that the relaxation then becomes exponential.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Phase Separation and Coexistence of Hydrodynamically Interacting Microswimmers

A striking feature of the collective behavior of spherical microswimmers is that for sufficiently strong self-propulsion they phase-separate into a dense cluster coexisting with a low-density dis- ordered surrounding. Extending our previous work, we use the squirmer as a model swimmer and the particle-based simulation method of multi-particle collision dynamics to explore the influence of hydrodynamics on their phase behavior in a quasi-two-dimensional geometry. The coarsening dynamics towards the phase-separated state is diffusive in an intermediate time regime followed by a final ballistic compactification of the dense cluster. We determine the binodal lines in a phase diagram of Péclet number versus density. Interestingly, the gas binodals are shifted to smaller densities for increasing mean density or dense-cluster size, which we explain using a recently introduced pressure balance [S. C. Takatori et al., Phys. Rev. Lett. 113, 028103 (2014)] extended by a hydrodynamic contribution. Furthermore, we find that for pushers and pullers the binodal line is shifted to larger Péclet numbers compared to neutral squirmers. Finally, when lowering the Péclet number, the dense phase transforms from a hexagonal "solid" to a disordered "fluid" state.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Taylor line swimming in microchannels and cubic lattices of obstacles

Microorganisms naturally move in microstructured fluids. Using the simulation method of multi-particle collision dynamics, we study an undulatory Taylor line swimming in a two-dimensional microchannel and in a cubic lattice of obstacles, which represent simple forms of a microstructured environment. In the microchannel the Taylor line swims at an acute angle along a channel wall with a clearly enhanced swimming speed due to hydrodynamic interactions with the bounding wall. While in a dilute obstacle lattice swimming speed is also enhanced, a dense obstacle lattice gives rise to geometric swimming. This new type of swimming is characterized by a drastically increased swimming speed. Since the Taylor line has to fit into the free space of the obstacle lattice, the swimming speed is close to the phase velocity of the bending wave traveling along the Taylor line. While adjusting its swimming motion within the lattice, the Taylor line chooses a specific swimming direction, which we classify by a lattice vector. When plotting the swimming velocity versus the magnitude of the lattice vector, all our data collapse on a single master curve. Finally, we also report more complex trajectories within the obstacle lattice.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Detention times of microswimmers close to surfaces: Influence of hydrodynamic interactions and noise

After colliding with a surface, microswimmers reside there during the detention time. They accumulate and may form complex structures such as biofilms. We introduce a general framework to calculate the distribution of detention times using the method of first-passage times and study how rotational noise and hydrodynamic interactions influence the escape from a surface. We compare generic swimmer models to the simple active Brownian particle. While the respective detention times of source dipoles are smaller, the ones of pullers are larger by up to several orders of magnitude, and pushers show both trends. We apply our results to the more realistic squirmer model, for which we use lubrication theory, and validate them by simulations with multi-particle collision dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Frequency-Dependent Selection at Rough Expanding Fronts

Microbial colonies are experimental model systems for studying the colonization of new territory by biological species through range expansion. We study a generalization of the two-species Eden model, which incorporates local frequency-dependent selection, in order to analyze how social interactions between two species influence surface roughness of growing microbial colonies. The model includes several classical scenarios from game theory. We then concentrate on an expanding public goods game, where either cooperators or defectors take over the front depending on the system parameters. We analyze in detail the critical behavior of the nonequilibrium phase transition between global cooperation and defection and thereby identify a new universality class of phase transitions dealing with absorbing states. At the transition, the number of boundaries separating sectors decays with a novel power law in time and their superdiffusive motion crosses over from Eden scaling to a nearly ballistic regime. In parallel, the width of the front initially obeys Eden roughening and, at later times, passes over to selective roughening.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Marangoni flow at droplet interfaces: Three-dimensional solution and applications

The Marangoni effect refers to fluid flow induced by a gradient in surface tension at a fluid-fluid interface. We determine the full three-dimensional Marangoni flow generated by a non-uniform surface tension profile at the interface of a self-propelled spherical emulsion droplet. For all flow fields inside, outside, and at the interface of the droplet, we give analytical formulas. We also calculate the droplet velocity vector $\mathbf{v}^D$, which describes the swimming kinematics of the droplet, and generalize the squirmer parameter $β$, which distinguishes between different swimmer types called neutral, pusher, or puller. In the second part of this paper, we present two illustrative examples, where the Marangoni effect is used in active emulsion droplets. First, we demonstrate how micelle adsorption can spontaneously break the isotropic symmetry of an initially surfactant-free emulsion droplet, which then performs directed motion. Second, we think about light-switchable surfactants and laser light to create a patch with a different surfactant type at the droplet interface. Depending on the setup such as the wavelength of the laser light and the surfactant type in the outer bulk fluid, one can either push droplets along unstable trajectories or pull them along straight or oscillatory trajectories regulated by specific parameters. We explore these cases for strongly absorbing and for transparent droplets.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Phase ordering of zig-zag and bow-shaped hard needles in two dimensions

We perform extensive Monte Carlo simulations of a two-dimensional bent hard-needle model in both its chiral zig-zag and its achiral bow-shape configurations and present their phase diagrams. We find evidence for a variety of stable phases: isotropic, quasi-nematic, smectic-C, anti-ferromorphic smectic-A, and modulated-nematic. This last phase consists of layers formed by supramolecular arches. They create a periodic modulation of the molecular polarity whose period is sensitively controlled by molecular geometry. We identify transition densities using correlation functions together with appropriately defined order parameters and compare them with predictions from Onsager theory. The contribution of the molecular excluded area to deviations from Onsager theory and simple liquid crystal phase morphology is discussed. We demonstrate the isotropic--quasi-nematic transition to be consistent with a Kosterlitz-Thouless disclination unbinding scenario.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Zipping and Entanglement in Flagellar Bundle of E. Coli: Role of Motile Cell Body

The course of a peritrichous bacterium such as E. coli crucially depends on the level of synchronization and self-organization of several rotating flagella. However, the rotation of each flagellum generates counter body movements which in turn affect the flagellar dynamics. Using a detailed numerical model of an E. coli, we demonstrate that flagellar entanglement, besides fluid flow relative to the moving body, dramatically changes the dynamics of flagella from that compared to anchored flagella. In particular, bundle formation occurs through a zipping motion in a remarkably rapid time, affected little by initial flagellar orientation. A simplified analytical model supports our observations. Finally, we illustrate how entanglement, hydrodynamic interactions, and body movement contribute to zipping and bundling.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Controlling inertial focussing using rotational motion

In inertial microfluidics lift forces cause a particle to migrate across streamlines to specific positions in the cross section of a microchannel. We control the rotational motion of a particle and demonstrate that this allows to manipulate the lift-force profile and thereby the particle's equilibrium positions. We perform two-dimensional simulation studies using the method of multi-particle collision dynamics. Particles with unconstrained rotational motion occupy stable equilibrium positions in both halfs of the channel while the center is unstable. When an external torque is applied to the particle, two equilibrium positions annihilate by passing a saddle-node bifurcation and only one stable fixpoint remains so that all particles move to one side of the channel. In contrast, non-rotating particles accumulate in the center and are pushed into one half of the channel when the angular velocity is fixed to a non-zero value.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dynamic clustering and chemotactic collapse of self-phoretic active particles

Recent experiments with self-phoretic particles at low concentrations show a pronounced dynamic clustering [I. Theurkauff \emph{et al.}, Phys.\ Rev.\ Lett.\ \textbf{108}, 268303 (2012)]. We model this situation by taking into account the translational and rotational diffusiophoretic motion, which the active particles perform in their self-generated chemical field. Our Brownian dynamics simulations show pronounced dynamic clustering only when these two phoretic contributions give rise to competing attractive and repulsive interactions, respectively. We identify two dynamic clustering states and characterize them by power-law-exponential distributions. In case of mere attraction a chemotactic collaps occurs directly from the gas-like into the collapsed state, which we also predict by mapping our Langevin dynamics on the Keller-Segel model for bacterial chemotaxis.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Feedback control of inertial microfluidics using axial control forces

Inertial microfluidics is a promising tool for many lab-on-a-chip applications. Particles in channel flows with Reynolds numbers above one undergo cross-streamline migration to a discrete set of equilibrium positions in square and rectangular channel cross sections. This effect has been used extensively for particle sorting and the analysis of particle properties. Using the lattice Boltzmann method, we determine equilibrium positions in square and rectangular cross sections and classify their types of stability for different Reynolds numbers, particle sizes, and channel aspect ratios. Our findings thereby help to design microfluidic channels for particle sorting. Furthermore, we demonstrate how an axial control force, which slows down the particles, shifts the stable equilibrium position towards the channel center. Ultimately, the particles then stay on the centerline for forces exceeding a threshold value. This effect is sensitive to particle size and channel Reynolds number and therefore suggests an efficient method for particle separation. In combination with a hysteretic feedback scheme, we can even increase particle throughput.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Self-induced polar order of active Brownian particles in a harmonic trap

Hydrodynamically interacting active particles in an external harmonic potential form a self-assembled fluid pump at large enough Péclet numbers. Here, we give a quantitative criterion for the formation of the pump and show that particle orientations align in the self-induced flow field in surprising analogy to ferromagnetic order where the active Péclet number plays the role of inverse temperature. The particle orientations follow a Boltzmann distribution $Φ(\mathbf{p}) \sim \exp(A p_z)$ where the ordering mean field $A$ scales with active Péclet number and polar order parameter. The mean flow field in which the particles' swimming directions align corresponds to a regularized stokeslet with strength proportional to swimming speed. Analytic mean-field results are compared with results from Brownian dynamics simulations with hydrodynamic interactions included and are found to capture the self-induced alignment very well.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Stability of liquid films covered by a carpet of self-propelled surfactant particles

We consider a carpet of self-propelled particles at the liquid-gas interface of a liquid film on a solid substrate. The particles excert an excess pressure on the interface and also move along the interface while the swimming direction performs rotational diffusion. We study the intricate influence of these self-propelled insoluble surfactants on the stability of the film surface and show that depending on the strength of in-surface rotational diffusion and the absolute value of the in-surface swimming velocity several characteristic instability modes can occur. In particular, rotational diffusion can either stabilize the film or induce instabilities of different character.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Hydrodynamics determines collective motion and phase behavior of active colloids

We study the collective motion of confined spherical microswimmers such as active colloids which we model by so-called squirmers. To simulate hydrodynamic flow fields including thermal noise, we use the method of multi-particle collision dynamics. We demonstrate that hydrodynamic near fields acting between squirmers as well as squirmers and bounding walls crucially determine their collective motion. In particular, with increasing density we observe a clear phase separation into a gas-like and cluster phase for neutral squirmers whereas strong pushers and pullers more gradually approach the hexagonal cluster state.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Rectification of self-propelled particles by symmetric barriers

The motion of self-propelled particles can be rectified by asymmetric or ratchet-like periodic patterns in space. Here we show that a non-zero average drift can already be induced in a periodic potential with symmetric barriers when the self-propulsion velocity is also symmetric and periodically modulated but phase-shifted against the potential. In the adiabatic limit of slow rotational diffusion we determine the mean drift analytically and discuss the influence of temperature. In the presence of asymmetric barriers modulating the self-propulsion can largely enhance the mean drift or even reverse it.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Sedimentation and polar order of active bottom-heavy particles

Self-propelled particles in an external gravitational field have been shown to display both an increased sedimentation length and polar order even without particle interactions. Here, we investigate self-propelled particles which additionally are bottom-heavy, that is they feel a torque aligning them to swim against the gravitational field. For bottom-heavy particles the gravitational field has the two opposite effects of i) sedimentation and ii) upward alignment of the particles' swimming direction. We perform a multipole expansion of the one-particle distribution with respect to orientation and derive expressions for sedimentation length and mean particle orientation which we check against Brownian Dynamics simulations. For large strength of gravity or small particle speeds and aligning torque, we observe sedimentation with increased sedimentation length compared with passive colloids but also active colloids without bottom-heaviness. Increasing, for example, swimming speed the sedimentation profile is inverted and the particles swim towards the top wall of the enclosing box. We find maximal orientational order at intermediate swimming speeds for both cases of particles with bottom-heaviness and those without. Ordering unsurprisingly is increased for the bottom-heavy particles, but this difference disappears at higher levels of activity and for very high activities ordering goes to zero in both cases.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Diffusive transport of light in three-dimensional disordered Voronoi structures

The origin of diffusive transport of light in dry foams is still under debate. In this paper, we consider the random walks of photons as they are reflected or transmitted by liquid films according to the rules of ray optics. The foams are approximately modeled by three-dimensional Voronoi tessellations with varying degree of disorder. We study two cases: a constant intensity reflectance and the reflectance of thin films. Especially in the second case, we find that in the experimentally important regime for the film thicknesses, the transport-mean-free path does not significantly depend on the topological and geometrical disorder of the Voronoi foams including the periodic Kelvin foam. This may indicate that the detailed structure of foams is not crucial for understanding the diffusive transport of light. Furthermore, our theoretical values for transport-mean-free path fall in the same range as the experimental values observed in dry foams. One can therefore argue that liquid films contribute substantially to the diffusive transport of light in {dry} foams.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Modelling bacterial flagellar growth

The growth of bacterial flagellar filaments is a self-assembly process where flagellin molecules are transported through the narrow core of the flagellum and are added at the distal end. To model this situation, we generalize a growth process based on the TASEP model by allowing particles to move both forward and backward on the lattice. The bias in the forward and backward jump rates determines the lattice tip speed, which we analyze and also compare to simulations. For positive bias, the system is in a non-equilibrium steady state and exhibits boundary-induced phase transitions. The tip speed is constant. In the no-bias case we find that the length of the lattice grows as $N(t)\propto\sqrt{t}$, whereas for negative drift $N(t)\propto\ln{t}$. The latter result agrees with experimental data of bacterial flagellar growth.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Motor-Driven Bacterial Flagella and Buckling Instabilities

Many types of bacteria swim by rotating a bundle of helical filaments also called flagella. Each filament is driven by a rotary motor and a very flexible hook transmits the motor torque to the filament. We model it by discretizing Kirchhoff's elastic-rod theory and develop a coarse-grained approach for driving the helical filament by a motor torque. A rotating flagellum generates a thrust force, which pushes the cell body forward and which increases with the motor torque. We fix the rotating flagellum in space and show that it buckles under the thrust force at a critical motor torque. Buckling becomes visible as a supercritical Hopf bifurcation in the thrust force. A second buckling transition occurs at an even higher motor torque. We attach the flagellum to a spherical cell body and also observe the first buckling transition during locomotion. By changing the size of the cell body, we vary the necessary thrust force and thereby obtain a characteristic relation between the critical thrust force and motor torque. We present a sophisticated analytical model for the buckling transition based on a helical rod which quantitatively reproduces the critical force-torque relation. Real values for motor torque, cell body size, and the geometry of the helical filament suggest that buckling should occur in single bacterial flagella. We also find that the orientation of pulling flagella along the driving torque is not stable and comment on the biological relevance for marine bacteria.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Nonlinear dynamics of a microswimmer in Poiseuille flow

We study the three-dimensional dynamics of a spherical microswimmer in cylindrical Poiseuille flow which can be mapped onto a Hamiltonian system. Swinging and tumbling trajectories are identified. In 2D they are equivalent to oscillating and circling solutions of a mathematical pendulum. Hydrodynamic interactions between the swimmer and confining channel walls lead to dissipative dynamics and result in stable trajectories, different for pullers and pushers. We demonstrate this behavior in the dipole approximation of the swimmer and with simulations using the method of multi-particle collision dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Periodic and Quasiperiodic Motion of an Elongated Microswimmer in Poiseuille Flow

We study the dynamics of a prolate spheroidal microswimmer in Poiseuille flow for different flow geometries. When moving between two parallel plates or in a cylindrical microchannel, the swimmer performs either periodic swinging or periodic tumbling motion. Although the trajectories of spherical and elongated swimmers are qualitatively similar, the swinging and tumbling frequency strongly depends on the aspect ratio of the swimmer. In channels with reduced symmetry the swimmers perform quasiperiodic motion which we demonstrate explicitely for swimming in a channel with elliptical cross section.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Swimming active droplet: A theoretical analysis

Recently, an active microswimmer was constructed where a micron-sized droplet of bromine water was placed into a surfactant-laden oil phase. Due to a bromination reaction of the surfactant at the interface, the surface tension locally increases and becomes non-uniform. This drives a Marangoni flow which propels the squirming droplet forward. We develop a diffusion-advection-reaction equation for the order parameter of the surfactant mixture at the droplet interface using a mixing free energy. Numerical solutions reveal a stable swimming regime above a critical Marangoni number M but also stopping and oscillating states when M is increased further. The swimming droplet is identified as a pusher whereas in the oscillating state it oscillates between being a puller and a pusher.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Archimedean-like colloidal tilings on substrates with decagonal and tetradecagonal symmetry

Two-dimensional colloidal suspensions subject to laser interference patterns with decagonal symmetry can form an Archimedean-like tiling phase where rows of squares and triangles order aperiodically along one direction [J. Mikhael et al., Nature 454, 501 (2008)]. In experiments as well as in Monte-Carlo and Brownian dynamics simulations, we identify a similar phase when the laser field possesses tetradecagonal symmetry. We characterize the structure of both Archimedean-like tilings in detail and point out how the tilings differ from each other. Furthermore, we also estimate specific particle densities where the Archimedean-like tiling phases occur. Finally, using Brownian dynamics simulations we demonstrate how phasonic distortions of the decagonal laser field influence the Archimedean-like tiling. In particular, the domain size of the tiling can be enlarged by phasonic drifts and constant gradients in the phasonic displacement. We demonstrate that the latter occurs when the interfering laser beams are not adjusted properly.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Force-extension curves of bacterial flagella

Bacterial flagella assume different helical shapes during the tumbling phase of a bacterium but also in response to varying environmental conditions. Force-extension measurements by Darnton and Berg explicitly demonstrate a transformation from the coiled to the normal helical state [N.C. Darnton and H.C. Berg, Biophys. J. {92}, 2230 (2007)]. We here develop an elastic model for the flagellum based on Kirchhoff's theory of an elastic rod that describes such a polymorphic transformation and use resistive force theory to couple the flagellum to the aqueous environment. We present Brownian dynamics simulations that quantitatively reproduce the force-extension curves and study how the ratio $Γ$ of torsional to bending rigidity and the extensional rate influence the response of the flagellum. An upper bound for $Γ$ is given. Using clamped flagella, we show in an adiabatic approximation that the mean extension, where a local coiled-to-normal transition occurs first, depends on the logarithm of the extensional rate.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Metachronal waves in a chain of rowers with hydrodynamic interactions

Filaments on the surface of a microorganism such as Paramecium or Ophalina beat highly synchronized and form so-called metachronal waves that travel along the surfaces. In order to study under what principal conditions these waves form, we introduce a chain of beads, called rowers, each periodically driven by an external force on a straight line segment. To implement hydrodynamic interactions between the beads, they are considered point-like. Two beads synchronize in antiphase or in phase depending on the positive or negative curvature of their driving-force potential. Concentrating on in-phase synchronizing rowers, we find that they display only transient synchronization in a bulk fluid. On the other hand, metachronal waves with wavelengths of 7-10 rower distances emerge, when we restrict the range of hydrodynamic interactions either artificially to nearest neighbors or by the presence of a bounding surface as in any relevant biological system.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Superdiffusion in a Honeycomb Billiard

We investigate particle transport in the honeycomb billiard that consists of connected channels placed on the edges of a honeycomb structure. The spreading of particles is superdiffusive due to the existence of ballistic trajectories which we term perfect paths. Simulations give a time exponent of 1.72 for the mean square displacement and a starlike, i.e., anisotropic particle distribution. We present an analytical treatment based on the formalism of continuous-time random walks and explain both the time exponent and the anisotropic distribution. In billiards with randomly distributed channels, conventional diffusion is always observed in the long-time limit, although for small disorder transient superdiffusional behavior exists. Our simulation results are again supported by an analytical analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Surmounting Barriers: The Benefit of Hydrodynamic Interactions

We experimentally and theoretically investigate the collective behavior of three colloidal particles that are driven by a constant force along a toroidal trap. Due to hydrodynamic interactions, a characteristic limit cycle is observed. When we additionally apply a periodic sawtooth potential, we find a novel caterpillar-like motional sequence that is dominated by hydrodynamic interactions and promotes the surmounting of potential barriers by the particles.

0 citations0 reviewsNo code linkedOpen access