Source author record

Hartmut Löwen

Hartmut Löwen 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

119works
19topics
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

119 published item(s)

preprint2023arXiv

Motility-induced shear thickening in dense colloidal suspensions

Phase transitions and collective dynamics of active colloidal suspensions are fascinating topics in soft matter physics, particularly for out-of-equilibrium systems, which can lead to rich rheological behaviours in the presence of steady shear flow. In this article, the role of self-propulsion in the rheological response of a dense colloidal suspension is investigated by using particle-resolved simulations. First, the interplay between activity and shear in the solid to the liquid transition of the suspension is analysed. While both self-propulsion and shear destroy order and melt the system by themselves above their critical values, self-propulsion lowers the stress barrier that needs to be overcome during the transition. Once the suspension reaches a non-equilibrium steady state the rheological response is analysed. While passive suspensions show a solid-like behaviour, turning on particle motility fluidises the system and, at low self-propulsion, the suspension behaves as a shear-thinning fluid. Increasing the self-propulsion of the colloids induces a transition from a shear-thinning to a shear-thickening behaviour, which we attribute to clustering in the suspensions induced by motility. This interesting phenomenon of motility-induced shear thickening (MIST) can be used to tailor the rheological response of colloidal suspensions.

preprint2022arXiv

A one-dimensional three-state run-and-tumble model with a `cell cycle'

We study a one-dimensional three-state run-and-tumble model motivated by the bacterium Caulobacter crescentus which displays a cell cycle between two non-proliferation mobile phase and a proliferating sedentary phase. Our model implements kinetic transitions between the two mobile and one sedentary states described in terms of their number densities, where mobility is allowed with different running speeds in forward and backward direction. We start by analyzing the stationary states of the system and compute the mean and squared-displacements for the distribution of all cells, as well as for the number density of settled cells. The latter displays a surprising super-ballistic scaling ~t^3 at early times. Including repulsive and attractive interactions between the mobile cell populations and the settled cells, we explore the stability of the system and then employ numerical methods to study structure formation in the fully nonlinear system. We find traveling waves of bacteria, whose occurrence is quantified in a nonequilibrium state diagram.

preprint2022arXiv

Active Brownian motion with memory delay induced by a viscoelastic medium

By now active Brownian motion is a well-established model to describe the motion of mesoscopic self-propelled particles in a Newtonian fluid. On the basis of the generalized Langevin equation, we present an analytic framework for active Brownian motion with memory delay assuming time-dependent friction kernels for both translational and orientational degrees of freedom to account for the time-delayed response of a viscoelastic medium. Analytical results are obtained for the orientational correlation function, mean displacement and mean-square displacement which we evaluate in particular for a Maxwell fluid characterized by a kernel which decays exponentially in time. Further, we identify a memory induced delay between the effective self-propulsion force and the particle orientation which we quantify in terms of a special dynamical correlation function. In principle our predictions can be verified for an active colloidal particle in various viscoelastic environments such as a polymer solution.

preprint2022arXiv

Active particles driven by competing spatially dependent self-propulsion and external force

We investigate how the competing presence of a nonuniform motility landscape and an external confining field affects the properties of active particles. We employ the active Ornstein-Uhlenbeck particle (AOUP) model with a periodic swim velocity profile to derive analytical approximations for the steady-state probability distribution of position and velocity, encompassing both the Unified Colored Noise Approximation and the theory of potential-free active particles with spatially dependent swim velocity recently developed. We test the theory by confining an active particle in a harmonic trap, which gives rise to interesting properties, such as a transition from a unimodal to a bimodal (and, eventually multimodal) spatial density, induced by decreasing the spatial period of the self propulsion. Correspondingly, the velocity distribution shows pronounced deviations from the Gaussian shape, even displaying a bimodal profile in the high-motility regions. Our results can be confirmed by real-space experiments on active colloidal Janus particles in external fields.

preprint2022arXiv

Active particles using reinforcement learning to navigate in complex motility landscapes

As the length scales of the smallest technology continue to advance beyond the micron scale it becomes increasingly important to equip robotic components with the means for intelligent and autonomous decision making with limited information. With the help of a tabular Q-learning algorithm, we design a model for training a microswimmer, to navigate quickly through an environment given by various different scalar motility fields, while receiving a limited amount of local information. We compare the performances of the microswimmer, defined via time of first passage to a target, with performances of suitable reference cases. We show that the strategy obtained with our reinforcement learning model indeed represents an efficient navigation strategy, that outperforms the reference cases. By confronting the swimmer with a variety of unfamiliar environments after the finalised training, we show that the obtained strategy generalises to different classes of random fields.

preprint2022arXiv

Brownian particles driven by spatially periodic noise

We discuss the dynamics of a Brownian particle under the influence of a spatially periodic noise strength in one dimension using analytical theory and computer simulations. In the absence of a deterministic force, the Langevin equation can be integrated formally exactly. We determine the short- and long-time behaviour of the mean displacement (MD) and mean-squared displacement (MSD). In particular we find a very slow dynamics for the mean displacement, scaling as $t^{-1/2}$ with time $t$. Placed under an additional external periodic force near the critical tilt value we compute the stationary current obtained from the corresponding Fokker-Planck equation and identify an essential singularity if the minimum of the noise strength is zero. Finally, in order to further elucidate the effect of the random periodic driving on the diffusion process, we introduce a phase factor in the spatial noise with respect to the external periodic force and identify the value of the phase shift for which the random force exerts its strongest effect on the long-time drift velocity and diffusion coefficient.

preprint2022arXiv

Collisions enhance self-diffusion in odd-diffusive systems

It is generally believed that collisions of particles reduce the self-diffusion coefficient. Here we show that in odd-diffusive systems, which are characterized by diffusion tensors with antisymmetric elements, collisions surprisingly can enhance the self-diffusion. In these systems, due to an inherent curving effect, the motion of particles is facilitated, instead of hindered by collisions leading to a mutual rolling effect. Using a geometric model, we analytically predict the enhancement of the self-diffusion coefficient with increasing density. This counterintuitive behaviour is demonstrated in the archetypal odd-diffusive system of Brownian particles under Lorentz force. We validate our findings by many body Brownian dynamics simulations in dilute systems.

preprint2022arXiv

Elastic deformations of loaded core-shell systems

Macroscopic elastic core-shell systems can be generated as toy models to be deformed and haptically studied by hand. On the mesoscale, colloidal core-shell particles and microgels are fabricated and investigated by different types of microscopy. We analyse, using linear elasticity theory, the response of spherical core-shell systems under the influence of a line density of force that is oriented radially and acts along the equator of the outer surface. Interestingly, deformational coupling of the shell to the core can determine the resulting overall appearance in response to the forces. We address various combinations of radii, stiffness, and Poisson ratio of core and shell and illustrate the resulting deformations. Macroscopically, the situation could be realized by wrapping a cord around the equator of a macroscopic model system and pulling it tight. On the mesoscale, colloidal microgel particles symmetrically confined to the interface between two immiscible fluids are pulled radially outward by surface tension.

preprint2022arXiv

Escape Dynamics in an Anisotropically Driven Brownian Magneto-System

Thermally activated escape of a Brownian particle over a potential barrier is well understood within Kramers theory. When subjected to an external magnetic field, the Lorentz force slows down the escape dynamics via a rescaling of the diffusion coefficient without affecting the exponential dependence on the barrier height. Here, we study the escape dynamics of a charged Brownian particle from a two-dimensional truncated harmonic potential under the influence of Lorentz force due to an external magnetic field. The particle is driven anisotropically by subjecting it to noises with different strengths along different spatial directions. We show that the escape time can largely be tuned by the anisotropic driving. While the escape process becomes anisotropic due to the two different noises, the spatial symmetry is restored in the limit of large magnetic fields. This is attributed to the Lorentz force induced coupling between the spatial degrees of freedom which makes the difference between two noises irrelevant at high magnetic fields. The theoretical predictions are verified by Brownian dynamics simulations. In principle, our predictions can be tested by experiments with a Brownian gyrator in the presence of a magnetic field.

preprint2022arXiv

Geometric criteria for the absence of effective many-body interactions in nonadditive hard particle mixtures

We consider a mixture of small and big classical particles in arbitrary spatial dimensions interacting via hard-body potentials with non-additive excluded-volume interactions. In particular, we focus on variants of the Asakura--Oosawa (AO) model where the interaction between the small particles is neglected but the big-small and big-big interactions are present and can be condensed into an effective depletion interaction among the big particles alone. The original AO model involves hard spherical particles in three spatial dimensions with interaction diameters $σ_\text{pp}=0$, $σ_\text{cc}>0$ and $σ_\text{pc}>σ_\text{cc}/2$ respectively, where $σ_{ij}$ with $\{i,j\}=\{\text{p},\text{c}\}$ (indicating the physical interpretation of the small and big particles as polymers (p) and colloids (c), respectively) is the minimum possible center-to-center distance between particle $i$ and particle $j$ allowed by the excluded-volume constraints. It is common knowledge that there are only pairwise effective depletion interactions between the big particles if the geometric condition $σ_\text{pc}/σ_\text{cc} < 1/\sqrt{3}$ is fulfilled. In this case, triplet and higher-order many body interactions are vanishing and the equilibrium statistics of the binary mixture can exactly be mapped onto that of an effective one-component system with the effective depletion pair-potential. Here we prove this geometric criterion rigorously and generalize it to polydisperse mixtures and to anisotropic particle shapes in any dimension, providing geometric criteria sufficient to guarantee the absence of triplet and higher-order many body interactions. For an external hard wall confining the full mixture, we also give criteria which guarantee that the system can be mapped onto one with effective external one-body interactions.

preprint2022arXiv

Hydrodynamics can Determine the Optimal Route for Microswimmer Navigation

Contrasting the well explored problem on how to steer a macroscopic agent like an airplane or a moon lander to optimally reach a target, "optimal microswimming", i.e. the quest for the optimal navigation strategy for microswimmers, remains unsolved. Here, we systematically explore this problem and show that the characteristic flow field of microswimmers crucially influences the required navigation strategy to reach a target fastest. The resulting optimal trajectories can have remarkable and non-intuitive shapes, which qualitatively differ from those of dry active particles or motile macroagents. Our results provide generic insights into the role of hydrodynamics and fluctuations on optimal navigation at the microscale and suggest that microorganisms might have survival advantages when strategically controlling their distance to remote walls. In particular, when fluctuations are present, choosing the optimal strategy, which appropriately respects hydrdynamics, can halve the time to reach the target compared to cases microswimmers head straight towards it.

preprint2022arXiv

Mutation induced infection waves in diseases like COVID-19

After more than 6 million deaths worldwide, the ongoing vaccination to conquer the COVID-19 disease is now competing with the emergence of increasingly contagious mutations, repeatedly supplanting earlier strains. Following the near-absence of historical examples of the long-time evolution of infectious diseases under similar circumstances, models are crucial to exemplify possible scenarios. Accordingly, in the present work we systematically generalize the popular susceptible-infected-recovered model to account for mutations leading to repeatedly occurring new strains, which we coarse grain based on tools from statistical mechanics to derive a model predicting the most likely outcomes. The model predicts that mutations can induce a super-exponential growth of infection numbers at early times, which self-amplify to giant infection waves which are caused by a positive feedback loop between infection numbers and mutations and lead to a simultaneous infection of the majority of the population. At later stage -- if vaccination progresses too slowly -- mutations can interrupt an ongoing decrease of infection numbers and can cause infection revivals which occur as single waves or even as whole wave trains featuring alternative periods of decreasing and increasing infection numbers. This panorama of possible mutation-induced scenarios should be tested in more detailed models to explore their concrete significance for specific infectious diseases. Further, our results might be useful for discussions regarding the importance of a release of vaccine-patents to reduce the risk of mutation-induced infection revivals but also to coordinate the release of measures following a downwards trend of infection numbers.

preprint2022arXiv

Steering self-organisation through confinement

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

preprint2022arXiv

Strongly enhanced dynamics of a charged Rouse dimer by an external magnetic field

While the dynamics of dimers and polymer chains in a viscous solvent is well understood within the celebrated Rouse model, the effect of an external magnetic field on the dynamics of a charged chain is much less understood. Here we generalize the Rouse model for a charged dimer to include the effect of an external magnetic field. Our analytically solvable model allows a fundamental insight into the magneto-generated dynamics of the dimer in the overdamped limit as induced by the Lorentz-force. Surprisingly, for a dimer of oppositely charged particles, we find an enormous enhancement of the dynamics of the dimer center which exhibits even a transient superballistic behavior. This is highly unusual in an overdamped system for there is neither inertia nor any internal or external driving. We attribute this to a significant translation and rotation coupling due to the Lorentz force. We also find that magnetic field reduces the mobility of a dimer along its orientation and its effective rotational diffusion coefficient. In principle, our predictions can be tested by experiments with colloidal particles and complex plasmas.

preprint2022arXiv

Structural correlations in highly asymmetric binary charged colloidal mixtures

We explore structural correlations of strongly asymmetric mixtures of binary charged colloids within the primitive model of electrolytes considering large charge and size ratios of 10 and higher. Using computer simulations with explicit microions, we obtain the partial pair correlation functions between the like-charged colloidal macroions. Interestingly the big-small correlation peak amplitude is smaller than that of the big-big and small-small macroion correlation peaks, which is unfamiliar for additive repulsive interactions. Extracting optimal effective microion-averaged pair interactions between the macroions, we find that on top of non-additive Yukawa-like repulsions an additional shifted Gaussian attractive potential between the small macroions is needed to accurately reproduce their correct pair correlations. For small Coulomb couplings, the behavior is reproduced in a coarse-grained theory with microion-averaged effective interactions between the macroions. However, the accuracy of the theory deteriorates with increasing Coulomb coupling. We emphasize the relevance of entropic interactions exerted by the microions on the macroions. Our results are experimentally verifiable in binary mixtures of micron-sized colloids and like-charge nanoparticles

preprint2022arXiv

The Parental Active Model: a unifying stochastic description of self-propulsion

We propose a new overarching model for self-propelled particles that flexibly generates a full family of "descendants". The general dynamics introduced in this paper, which we denote as "parental" active model (PAM), unifies two special cases commonly used to describe active matter, namely active Brownian particles (ABPs) and active Ornstein-Uhlenbeck particles (AOUPs). We thereby document the existence of a deep and close stochastic relationship between them, resulting in the subtle balance between fluctuations in the magnitude and direction of the self-propulsion velocity. Besides illustrating the relation between these two common models, the PAM can generate additional offspring, interpolating between ABP and AOUP dynamics, that could provide more suitable models for a large class of living and inanimate active matter systems, possessing characteristic distributions of their self-propulsion velocity. Our general model is evaluated in the presence of a harmonic external confinement. For this reference example, we present a two-state phase diagram which sheds light on the transition in the shape of the positional density distribution, from a unimodal Gaussian for AOUPs to a Mexican-hat-like profile for ABPs.

preprint2021arXiv

Dynamical density functional theory for "dry" and "wet" active matter

In the last 50 years, equilibrium density functional theory (DFT) has been proven to be a powerful, versatile and predictive approach for the statics and structure of classical particles. This theory can be extended to the nonequilibrium dynamics of completely overdamped Brownian colloidal particles towards so-called dynamical density functional theory (DDFT). The success of DDFT makes it a promising candidate for a first-principle description of active matter. In this lecture, we shall first recapitulate classical DDFT for passive colloidal particles typically described by Smoluchowski equation. After a basic derivation of DDFT from the Smoluchowski equation, we discuss orientational degrees of freedom and the effect of hydrodynamic interactions for passive particles. This brings us into an ideal position to generalize DDFT towards active matter. In particular we distinguish between "dry active matter" which is composed of self-propelled particles that contain no hydrodynamic flow effects of a surrounding solvent and "wet active matter" where the hydrodynamic flow fields generated by the microswimmers are taken into account. For the latter, DDFT is a tool which unifies thermal fluctuations, direct particle interactions, external driving fields and hydrodynamic effects arising from internal self-propulsion discriminating between "pushers" and "pullers". A number of recent applications is discussed including transient clustering of self-propelled rods and the spontaneous formation of a hydrodynamic pump in confined microswimmers.

preprint2021arXiv

Dynamics of active particles with space-dependent swim velocity

We study the dynamical properties of an active particle subject to a swimming speed explicitly depending on the particle position. The oscillating spatial profile of the swim velocity considered in this paper takes inspiration from experimental studies based on Janus particles whose speed can be modulated by an external source of light. We suggest and apply an appropriate model of an active Ornstein Uhlenbeck particle to the present case. This allows us to predict the stationary properties, by finding the exact solution of the steady-state probability distribution of particle position and velocity. From this, we obtain the spatial density profile and show that its form is consistent with the one found in the framework of other popular models. The reduced velocity distribution highlights the emergence of non-Gaussianity in our generalized AOUP model which becomes more evident as the spatial dependence of the velocity profile becomes more pronounced. Then, we focus on the time-dependent properties of the system. Velocity autocorrelation functions are studied in the steady-state combining numerical and analytical methods derived under suitable approximations. We observe a non-monotonic decay in the temporal shape of the velocity autocorrelation function which depends on the ratio between the persistence length and the spatial period of the swim velocity. Finally, we numerically and analytically study the mean square displacement and the long-time diffusion coefficient. The ballistic regime, observed in the small-time region, is deeply affected by the properties of the swim velocity landscape which induces also a crossover to a sub-ballistic but superdiffusive regime for intermediate times. Finally, the long-time diffusion coefficient decreases as the amplitude of the swim velocity oscillations increases because the diffusion is determined by those regions where the particles are slow.

preprint2021arXiv

Interface-induced hysteretic volume phase transition of microgels: simulation and experiment

Thermo-responsive microgel particles can exhibit a drastic volume shrinkage upon increasing the solvent temperature. Recently we found that the spreading of poly(N-isopropylacrylamide)(PNiPAm) microgels at a liquid interface under the influence of surface tension hinders the temperature-induced volume phase transition. In addition, we observed a hysteresis behavior upon temperature cycling, i.e. a different evolution in microgel size and shape depending on whether the microgel was initially adsorbed to the interface in expanded or collapsed state. Here, we model the volume phase transition of such microgels at an air/water interface by monomer-resolved Brownian dynamics simulations and compare the observed behavior with experiments. We reproduce the experimentally observed hysteresis in the microgel dimensions upon temperature variation. Our simulations did not observe any hysteresis for microgels dispersed in the bulk liquid, suggesting that it results from the distinct interfacial morphology of the microgel adsorbed at the liquid interface. An initially collapsed microgel brought to the interface and subjected to subsequent swelling and collapsing (resp. cooling and heating) will end up in a larger size than it had in the original collapsed state. Further temperature cycling, however, only shows a much reduced hysteresis, in agreement with our experimental observations. We attribute the hysteretic behavior to a kinetically trapped initial collapsed configuration, which relaxes upon expanding in the swollen state. We find a similar behavior for linear PNiPAm chains adsorbed to an interface. Our combined experimental - simulation investigation provides new insights into the volume phase transition of PNiPAm materials adsorbed to liquid interfaces.

preprint2021arXiv

Order-preserving dynamics in one dimension -- single-file diffusion and caging from the perspective of dynamical density functional theory

Dynamical density functional theory (DDFT) is a powerful variational framework to study the nonequilibrium properties of colloids by only considering a time-dependent one-body number density. Despite the large number of recent successes, properly modeling the long-time dynamics in interacting systems within DDFT remains a notoriously difficult problem, since structural information, accounting for temporary or permanent particle cages, gets lost. Here we address such a caging scenario by reducing it to a clean one-dimensional problem, where the particles are naturally ordered (arranged on a line) by perfect cages created by their two next neighbors. In particular, we construct a DDFT approximation based on an equilibrium system with an asymmetric pair potential, such that the corresponding one-body densities still carry the footprint of particle order. Applied to a system of confined hard rods, this order-preserving dynamics (OPD) yields exact results at the system boundaries, in addition to the imprinted correct long-time behavior of density profiles representing individual particles. In an open system, our approach correctly reproduces the reduced long-time diffusion coefficient and subdiffusion, characteristic for a single-file setup. These observations cannot be made using current forms of DDFT without particle order.

preprint2021arXiv

Particle-resolved topological defects of smectic colloidal liquid crystals in extreme confinement

Confined samples of liquid crystals are characterized by a variety of topological defects and can be exposed to external constraints such as extreme confinements with nontrivial topology. Here we explore the intrinsic structure of smectic colloidal layers dictated by the interplay between entropy and an imposed external topology. Considering an annular confinement as a basic example, a plethora of competing states is found with nontrivial defect structures ranging from laminar states to multiple smectic domains and arrays of edge dislocations which we refer to as Shubnikov states in formal analogy to the characteristic of type-II superconductors. Our particle-resolved results, gained by a combination of real-space microscopy of thermal colloidal rods and fundamental-measure-based density functional theory of hard anisotropic bodies, agree on a quantitative level.

preprint2020arXiv

Climbing two hills is faster than one: collective barrier-crossing by colloids driven through a microchannel

Ohm's law is one of the most central transport rules stating that the total resistance of sequential single resistances is additive. While this rule is most commonly applied to electronic circuits, it also applies to other transport phenomena such as the flow of colloids or nanoparticles through channels containing multiple obstacles, as long as these obstacles are sufficiently far apart. Here we explore the breakdown of Ohm's law for fluids of repulsive colloids driven over two energetic barriers in a microchannel, using real-space microscopy experiments, particle-resolved simulations, and dynamical density functional theory. If the barrier separation is comparable to the particle correlation length, the resistance is highly non-additive, such that the resistance added by the second barrier can be significantly higher or lower than that of the first. Surprisingly, in some cases the second barrier can even add a {\it negative} resistance, such that two identical barriers are easier to cross than a single one. We explain this counterintuitive observation in terms of the structuring of particles trapped between the barriers.

preprint2020arXiv

Dynamics of a microswimmer-microplatelet composite

Guiding active microswimmers by external fields to requested target locations is a promising strategy to realize complex transport on the microscale. To this end, one possibility consists of attaching the microswimmers to orientable passive components. Accordingly, we analyze theoretically, using a minimal model, the dynamics of a microswimmer when rigidly attached tangentially to a (significantly larger) microplatelet, here represented by a thin circular disk. On this way, we determine the flow field in the whole space induced by a Stokeslet that is located above the center of a spatially fixed rigid disk of no-slip surface conditions. Finally, we determine and analyze possible trajectories of the overall composite. To this end, the platelet is additionally endowed with a permanent magnetic moment, which allows to steer the motion of the whole composite by a homogeneous external magnetic field. As previous experimental studies suggest, related setups may be helpful to guide sperm cells to requested targets or for the purpose of coordinated drug delivery.

preprint2020arXiv

Persistent anti-correlations in Brownian dynamics simulations of dense colloidal suspensions revealed by noise suppression

Transport properties of a hard-sphere colloidal fluid are investigated by Brownian dynamics simulations. We implement a novel algorithm for the time-dependent velocity-autocorrelation function (VACF) essentially eliminating the noise of the bare random motion. The measured VACF reveals persistent anti-correlations manifested by a negative algebraic power-law tail $t^{-5/2}$ at all densities. At small packing fractions, the simulations fully agree with the analytic low-density prediction, yet the amplitude of the tail becomes dramatically suppressed as the packing fraction is increased. The mode-coupling theory of the glass transition provides a qualitative explanation for the strong variation in terms of the static compressibility as well as the slowing down of the structural relaxation.

preprint2020arXiv

Rototaxis: localization of active motion under rotation

The ability to navigate in complex, inhomogeneous environments is fundamental to survival at all length scales, giving rise to the rapid development of various subfields in bio-locomotion such as the well established concept of chemotaxis. In this work, we extend this existing notion of taxis to rotating environments and introduce the idea of roto-taxis to bio-locomotion. In particular, we explore both overdamped and inertial dynamics of a model synthetic self-propelled particle in the presence of constant global rotation, focusing on the particle's ability to localize near a rotation center as a survival strategy. We find that in the overdamped regime, the swimmer is in general able to generate a self restoring active torque that enables it to remain on stable epicyclical-like trajectories. On the other hand, for underdamped motion with inertial effects, the intricate competition between self-propulsion and inertial forces, in conjunction with the rototactic torque leads to complex dynamical behavior with non-trivial phase space of initial conditions which we reveal by numerical simulations. Our results are relevant for a wide range of setups, from vibrated granular matter on turntables to microorganisms or animals swimming near swirls or vortices.

preprint2020arXiv

Strategic Spatiotemporal Vaccine Distribution Increases the Survival Rate in an Infectious Disease like Covid-19

Covid-19 has caused hundred of thousands of deaths and an economic damage amounting to trillions of dollars, creating a desire for the rapid development of vaccine. Once available, vaccine is gradually produced, evoking the question on how to distribute it best. While official vaccination guidelines largely focus on the question to whom vaccines should be provided first (e.g. to risk groups), here we propose a strategy for their distribution in time and space, which sequentially prioritizes regions with a high local infection growth rate. To demonstrate this strategy, we develop a simple statistical model describing the time-evolution of infection patterns and their response to vaccination, for infectious diseases like Covid-19. For inhomogeneous infection patterns, locally well-mixed populations and basic reproduction numbers $R_0\sim 1.5-4$ the proposed strategy at least halves the number of deaths in our simulations compared to the standard practice of distributing vaccines proportionally to the population density. For $R_0\sim 1$ we still find a significant increase of the survival rate. The proposed vaccine distribution strategy can be further tested in detailed modelling works and could excite discussions on the importance of the spatiotemporal distribution of vaccines for official guidelines.

preprint2019arXiv

Active Brownian motion with orientation-dependent motility: theory and experiments

Combining experiments on active colloids, whose propulsion velocity can be controlled via a feedback loop, and theory of active Brownian motion, we explore the dynamics of an overdamped active particle with a motility that depends explicitly on the particle orientation. In this case, the active particle moves faster when oriented along one direction and slower when oriented along another, leading to an anisotropic translational dynamics which is coupled to the particle's rotational diffusion. We propose a basic model of active Brownian motion for orientation-dependent motility. Based on this model, we obtain analytic results for the mean trajectories, averaged over the Brownian noise for various initial configurations, and for the mean-square displacements including their anisotropic non-Gaussian behavior. The theoretical results are found to be in good agreement with the experimental data. Our findings establish a methodology to engineer complex anisotropic motilities of active Brownian particles, with potential impact in the study of the swimming behavior of microorganisms subjected to anisotropic driving fields.

preprint2019arXiv

Clustering-induced velocity-reversals of active colloids mixed with passive particles

Recent experiments have shown that colloidal suspensions can spontaneously self-assemble into dense clusters of various internal structures, sizes and dynamical properties when doped with active Janus particles. Characteristically, these clusters move ballistically during their formation, but dynamically revert their velocity and temporarily move opposite to the self-propulsion direction of the Janus particles they contain. Here we explore a simple effective model of colloidal mixtures which allows reproducing most aspects seen in experiments, including the morphology and the velocity-reversal of the clusters. We attribute the latter to the nonreciprocal phoretic attractions of the passive particles to the active colloids' caps, taking place even at close contact and pushing the active particles backwards. When the phoretic interactions are repulsive, in turn, they cause dynamical aggregation of passive colloids in the chemical density minima produced by the active particles, as recently seen in experiments; in other parameter regimes they induce travelling fronts of active particles pursued by passive ones coexisting with an active gas.

preprint2019arXiv

Inertial effects of self-propelled particles: from active Brownian to active Langevin motion

Active particles which are self-propelled by converting energy into mechanical motion represent an expanding research realm in physics and chemistry. For micron-sized particles moving in a liquid ("microswimmers"), most of the basic features have been described by using the model of overdamped active Brownian motion. However, for macroscopic particles or microparticles moving in a gas, inertial effects become relevant such that the dynamics is underdamped. Therefore, recently, active particles with inertia have been described by extending the active Brownian motion model to active Langevin dynamics which include inertia. In this perspective article recent developments of active particles with inertia ("microflyers") are summarized both for single particle properties and for collective effects of many particles. There include: inertial delay effects between particle velocity and self-propulsion direction, tuning of the long-time self-diffusion by the moment of inertia, effects of fictitious forces in non-inertial frames, and the influence of inertia on motility-induced phase separation. Possible future developments and perspectives are also proposed and discussed.

preprint2016arXiv

Active Particles in Complex and Crowded Environments

Differently from passive Brownian particles, active particles, also known as self-propelled Brownian particles or microswimmers and nanoswimmers, are capable of taking up energy from their environment and converting it into directed motion. Because of this constant flow of energy, their behavior can only be explained and understood within the framework of nonequilibrium physics. In the biological realm, many cells perform directed motion, for example, as a way to browse for nutrients or to avoid toxins. Inspired by these motile microorganisms, researchers have been developing artificial particles that feature similar swimming behaviors based on different mechanisms; these manmade micro- and nanomachines hold a great potential as autonomous agents for healthcare, sustainability, and security applications. With a focus on the basic physical features of the interactions of self-propelled Brownian particles with a crowded and complex environment, this comprehensive review will put the reader at the very forefront of the field, providing a guided tour through its basic principles, the development of artificial self-propelling micro- and nanoparticles, and their application to the study of nonequilibrium phenomena, as well as the open challenges that the field is currently facing.

preprint2016arXiv

Analysis of a Living Fluid Continuum Model

Generalized Navier-Stokes equations which were proposed recently to describe active turbulence in living fluids are analyzed rigorously. Results on wellposedness and stability in the $L^2(\mathbb{R}^n)$-setting are derived. Due to the presence of a Swift-Hohenberg term global wellposedness in a strong setting for arbitrary initial data in $L^2_σ(\mathbb{R}^n)$ is available. Based on the existence of global strong solutions, results on linear and nonlinear (in-) stability for the disordered steady state and the manifold of ordered polar steady states are derived, depending on the involved parameters.

preprint2016arXiv

Close packing of rods on spherical surfaces

We study the optimal packing of short, hard spherocylinders confined to lie tangential to a spherical surface, using simulated annealing and molecular dynamics simulations. For clusters of up to twelve particles, we map out the changes in the geometry of the closest-packed configuration as a function of the aspect ratio $L/D$, where $L$ is the cylinder length and $D$ the diameter of the rods. We find a rich variety of cluster structures. For larger clusters, we find that the best-packed configurations up to around 100 particles are highly dependent on the exact number of particles and aspect ratio. For even larger clusters, we find largely disordered clusters for very short rods ($L/D = 0.25$), while slightly longer rods ($L/D = 0.5$ or $1$) prefer a global baseball-like geometry of smectic-like domains, similar to the behavior of large-scale nematic shells. Our results provide predictions for experimentally realizable systems of colloidal rods trapped at the interface of emulsion droplets.

preprint2016arXiv

Colloidal suspensions of C-particles: Entanglement, percolation and microrheology

We explore structural and dynamical behavior of concentrated colloidal suspensions made up by C-shape particles using Brownian dynamics computer simulations and theory. In particular, we focus on the entanglement process between nearby particles for almost closed C-shapes with a small opening angle. Depending on the opening angle and the particle concentration, there is a percolation transition for the cluster of entangled particles which shows the classical scaling characteristics. In a broad density range below the percolation threshold, we find a stretched exponential function for the dynamical decorrelation of the entanglement process. Finally, we study a set-up typical in microrheology by dragging a single tagged particle with constant speed through the suspension. We measure the cluster connected to and dragged with this tagged particle. In agreement with a phenomenological theory, the size of the dragged cluster depends on the dragging direction and increases markedly with the dragging speed.

preprint2016arXiv

Dynamic Elastic Moduli in Magnetic Gels: Normal Modes and Linear Response

In the perspective of developing smart hybrid materials with customized features, ferrogels and magnetorheological elastomers allow a synergy of elasticity and magnetism. The interplay between elastic and magnetic properties gives rise to a unique reversible control of the material behavior by applying an external magnetic field. Albeit few works have been performed on the time-dependent properties so far, understanding the dynamic behavior is the key to model many practical situations, e.g. applications as vibration absorbers. Here we present a way to calculate the frequency-dependent elastic moduli based on the decomposition of the linear response to an external stress in normal modes. We use a minimal three-dimensional dipole-spring model to theoretically describe the magnetic and elastic interactions on the mesoscopic level. Specifically, the magnetic particles carry permanent magnetic dipole moments and are spatially arranged in a prescribed way, before they are linked by elastic springs. An external magnetic field aligns the magnetic moments. On the one hand, we study regular lattice-like particle arrangements to compare with previous results in the literature. On the other hand, we calculate the dynamic elastic moduli for irregular, more realistic particle distributions. Our approach measures the tunability of the linear dynamic response as a function of the particle arrangement, the system orientation with respect to the external magnetic field, as well as the magnitude of the magnetic interaction between the particles. The strength of the present approach is that it explicitly connects the relaxational modes of the system with the rheological properties as well as with the internal rearrangement of the particles in the sample, providing new insight into the dynamics of these remarkable materials.

preprint2016arXiv

Dynamical Density Functional Theory For Microswimmers

Dynamical density functional theory (DDFT) has been successfully derived and applied to describe on the one hand passive colloidal suspensions, including hydrodynamic interactions between individual particles. On the other hand, active "dry" crowds of self-propelled particles have been characterized using DDFT. Here we go one essential step further and combine these two approaches. We establish a DDFT for active microswimmer suspensions. For this purpose, simple minimal model microswimmers are introduced. These microswimmers self-propel by setting the surrounding fluid into motion. They hydrodynamically interact with each other through their actively self-induced fluid flows and via the common "passive" hydrodynamic interactions. An effective soft steric repulsion is also taken into account. We derive the DDFT starting from common statistical approaches. Our DDFT is then tested and applied by characterizing a suspension of microswimmers the motion of which is restricted to a plane within a three-dimensional bulk fluid. Moreover, the swimmers are confined by a radially symmetric trapping potential. In certain parameter ranges, we find rotational symmetry breaking in combination with the formation of a "hydrodynamic pumping state", which has previously been observed in the literature as a result of particle-based simulations. An additional instability of this pumping state is revealed.

preprint2016arXiv

Emerging activity in bilayered dispersions with wake-mediated interactions

In a bilayered system of particles with wake-mediated interactions, the action-reaction symmetry for the effective forces between particles of different layers is broken. Under quite general conditions we show that, if the interaction nonreciprocity exceeds a certain threshold, this creates an active dispersion of self-propelled clusters of Brownian particles. The emerging activity promotes unusual melting scenarios and an enormous diffusivity in the dense fluid. Our results are obtained by computer simulation and analytical theory, and can be verified in experiments with colloidal dispersions and complex plasmas.

preprint2016arXiv

Fission and fusion scenarios for magnetic microswimmer clusters

Fission and fusion processes of particles clusters occur in many areas of physics and chemistry from subnuclear to astronomic length scales. Here we study fission and fusion of magnetic microswimmer clusters as governed by their hydrodynamic and dipolar interactions. Rich scenarios are found which depend crucially on whether the swimmer is a pusher or a puller. In particular a linear magnetic chain of pullers is stable while a pusher chain shows a cascade of fission (or disassembly) processes as the self-propulsion velocity is increased. Contrarily, magnetic ring clusters show fission for any type of swimmer. Moreover, we find a plethora of possible fusion (or assembly) scenarios if a single swimmer collides with a ringlike cluster and two rings spontaneously collide. Our predictions are obtained by computer simulations and verifiable in experiments on active colloidal Janus particles and magnetotactic bacteria.

preprint2016arXiv

Following fluctuating signs: anomalous active superdiffusion of swimmers in anisotropic media

Active (i.e., self-propelled or swimming) particles moving through an isotropic fluid exhibit conventional diffusive behavior. We report anomalous diffusion of an active particle moving in an anisotropic, nematic background. Whilst the translational motion parallel to the nematic director shows ballistic behavior, the long-time transverse motion is super-diffusive, with an anomalous scaling $\propto t \ln t$ of the mean squared displacement with time $t$. This behavior is predicted by an analytical theory that we present here, and is corroborated by numerical simulation of active particle diffusion in a simple lattice model for a nematic liquid crystal. It is universal for any collection of self-propelled elements (e.g., bacteria or active rods) moving in a nematic background, provided only that the swimmers are sufficiently dilute that their interactions with each other can be neglected, and that they do not perform "hairpin" turns.

preprint2016arXiv

Hard rectangles near curved hard walls: tuning the sign of the Tolman length

Combining analytic calculations, computer simulations, and classical density functional theory we determine the interfacial tension of orientable two-dimensional hard rectangles near a curved hard wall. Both a circular cavity holding the particles and a hard circular obstacle surrounded by particles are considered. We focus on moderate bulk densities (corresponding to area fractions up to 50 percent) where the bulk phase is isotropic and vary the aspect ratio of the rectangles and the curvature of the wall. The Tolman length, which gives the leading curvature correction of the interfacial tension, is found to change sign at a finite density, which can be tuned via the aspect ratio of the rectangles.

preprint2016arXiv

Nucleation pathway and kinetics of phase-separating active Brownian particles

Suspensions of purely repulsive but self-propelled Brownian particles might undergo phase separation, a phenomenon that strongly resembles the phase separation of passive particles with attractions. Here we employ computer simulations to study the nucleation kinetics and the microscopic pathway active Brownian disks take in two dimensions when quenched from the homogeneous suspension to propulsion speeds beyond the binodal. We find the same qualitative behavior for the nucleation rate as a function of density as for a passive suspension undergoing liquid-vapor separation, suggesting that the scenario of an effective free energy also extends to the kinetics of phase separation. We study the transition in more detail through a committor analysis and find that transition states are best described by a combination of cluster size and the radial polarization of particles in the cluster.

preprint2016arXiv

Phototaxis of synthetic microswimmers in optical landscapes

Many microorganisms, with phytoplankton and zooplankton as prominent examples, display phototactic behaviour, that is, the ability to perform directed motion within a light gradient. Here we experimentally demonstrate that sensing of light gradients can also be achieved in a system of synthetic photo-activated microparticles being exposed to an inhomogeneous laser field. We observe a strong orientational response of the particles because of diffusiophoretic torques, which in combination with an intensity-dependent particle motility eventually leads to phototaxis. Since the aligning torques saturate at high gradients, a strongly rectified particle motion is found even in periodic asymmetric intensity landscapes. Our results are in excellent agreement with numerical simulations of a minimal model and should similarly apply to other particle propulsion mechanisms. Because light fields can be easily adjusted in space and time, this also allows to extend our approach to dynamical environments.

preprint2016arXiv

Superelastic stress-strain behavior in ferrogels of different types of magneto-elastic coupling

Colloidal magnetic particles embedded in an elastic polymer matrix constitute a smart material called ferrogel. It responds to an applied external magnetic field by changes in elastic properties, which can be exploited for various applications like dampers, vibration absorbers, or actuators. Under appropriate conditions, the stress-strain behavior of a ferrogel can display a fascinating feature: superelasticity, the capability to reversibly deform by a huge amount while barely altering the applied load. In a previous work, using numerical simulations, we investigated this behavior assuming that the magnetic moments carried by the embedded particles can freely reorient to minimize their magnetic interaction energy. Here, we extend the analysis to ferrogels where restoring torques by the surrounding matrix hinder rotations towards a magnetically favored configuration. For example, the particles can be chemically cross-linked into the polymer matrix and the magnetic moments can be fixed to the particle axes. We demonstrate that these systems still feature a superelastic regime. As before, the nonlinear stress-strain behavior can be reversibly tailored during operation by external magnetic fields. Yet, the different coupling of the magnetic moments causes different types of response to external stimuli. For instance, an external magnetic field applied parallel to the stretching axis hardly affects the superelastic regime but stiffens the system beyond it. Other smart materials featuring superelasticity, e.g. metallic shape-memory alloys, have already found widespread applications. Our soft polymer systems offer many additional advantages like a typically higher deformability and enhanced biocompatibility combined with high tunability.

preprint2015arXiv

Active dipole clusters: from helical motion to fission

The structure of a finite particle cluster is typically determined by total energy minimization. Here we consider the case where a cluster of soft sphere dipoles becomes active, i.e. when the individual particles exhibit an additional self-propulsion along their dipole moments. We numerically solve the overdamped equations of motion for soft-sphere dipoles in a solvent. Starting from an initial metastable dipolar cluster, the self-propulsion generates a complex cluster dynamics. The final cluster state has in general a structure widely different to the initial one, the details depend on the model parameters and on the protocol of how the self-propulsion is turned on. The center-of-mass of the cluster moves on a helical path, the details of which are governed by the initial cluster magnetization. An instantaneous switch to a high self-propulsion leads to fission of the cluster. However, fission does not occur if the self-propulsion is increased slowly to high strengths. Our predictions can be verified through experiments with self-phoretic colloidal Janus-particles and for macroscopic self-propelled dipoles in a highly viscous solvent.

preprint2015arXiv

Classical Liquids in Fractal Dimension

We introduce fractal liquids by generalizing classical liquids of integer dimensions $d = 1, 2, 3$ to a fractal dimension $d_f$. The particles composing the liquid are fractal objects and their configuration space is also fractal, with the same non-integer dimension. Realizations of our generic model system include microphase separated binary liquids in porous media, and highly branched liquid droplets confined to a fractal polymer backbone in a gel. Here we study the thermodynamics and pair correlations of fractal liquids by computer simulation and semi-analytical statistical mechanics. Our results are based on a model where fractal hard spheres move on a near-critical percolating lattice cluster. The predictions of the fractal Percus-Yevick liquid integral equation compare well with our simulation results.

preprint2015arXiv

Depolarized light scattering from prolate anisotropic particles: the influence of the particle shape on the field autocorrelation function

We provide a theoretical analysis for the intermediate scattering function typically measured in depolarized dynamic light scattering experiments. We calculate the field autocorrelation function $g_1^{\rm VH}(Q,t)$ in dependence on the wave vector $Q$ and the time $t$ explicitly in a vertical-horizontal scattering geometry for differently shaped solids of revolution. The shape of prolate cylinders, spherocylinders, spindles, and double cones with variable aspect ratio is expanded in rotational invariants $f_{lm}(r)$. By Fourier transform of these expansion coefficients, a formal multipole expansion of the scattering function is obtained, which is used to calculate the weighting coefficients appearing in the depolarized scattering function. In addition to translational and rotational diffusion, especially the translational-rotational coupling of shape-anisotropic objects is considered. From the short-time behavior of the intermediate scattering function, the first cumulants $Γ(Q)$ are calculated. In a depolarized scattering experiment, they deviate from the simple proportionality to $Q^2$. The coefficients $f_{lm}(Q)$ strongly depend on the geometry and aspect ratio of the particles. The time dependence, in addition, is governed by the translational and rotational diffusion tensors, which are calculated by means of bead models for differently shaped particles in dependence on their aspect ratio. Therefore, our analysis shows how details of the particle shape---beyond their aspect ratio---can be determined by a precise scattering experiment. This is of high relevance in understanding smart materials which involve suspensions of anisotropic colloidal particles.

preprint2015arXiv

Dynamical mean-field theory and weakly non-linear analysis for the phase separation of active Brownian particles

Recently, we have derived an effective Cahn-Hilliard equation for the phase separation dynamics of active Brownian particles by performing a weakly non-linear analysis of the effective hydrodynamic equations for density and polarization [Phys. Rev. Lett. 112, 218304 (2014)]. Here we develop and explore this strategy in more detail and show explicitly how to get to such a large-scale, mean-field description starting from the microscopic dynamics. The effective free energy emerging from this approach has the form of a conventional Ginzburg-Landau function. On the coarsest scale, our results thus agree with the mapping of active phase separation onto that of passive fluids with attractive interactions through a global effective free energy (mobility-induced phase transition). Particular attention is paid to the square-gradient term necessary for the dynamics. We finally discuss results from numerical simulations corroborating the analytical results.

preprint2015arXiv

Dynamics of a linear magnetic "microswimmer molecule"

In analogy to nanoscopic molecules that are composed of individual atoms, we consider an active "microswimmer molecule". It is built up from three individual magnetic colloidal microswimmers that are connected by harmonic springs and hydrodynamically interact with each other. In the ground state, they form a linear straight molecule. We analyze the relaxation dynamics for perturbations of this straight configuration. As a central result, with increasing self-propulsion, we observe an oscillatory instability in accord with a subcritical Hopf bifurcation scenario. It is accompanied by a corkscrew-like swimming trajectory of increasing radius. Our results can be tested experimentally, using for instance magnetic self-propelled Janus particles, supposably linked by DNA molecules.

preprint2015arXiv

Effect of self-propulsion on equilibrium clustering

In equilibrium, colloidal suspensions governed by short-range attractive and long-range repulsive interactions form thermodynamically stable clusters. Using Brownian dynamics computer simulations, we investigate how this equilibrium clustering is affected when such particles are self-propelled. We find that the clustering process is stable under self-propulsion. For the range of interaction parameters studied and at low particle density, the cluster size increases with the speed of self-propulsion (activity) and for higher activity the cluster size decreases, showing a non-monotonic variation of cluster size with activity. This clustering behaviour is distinct from the pure kinetic (or motility-induced) clustering of self-propelling particles which is observed at significantly higher activities and densities. We present an equilibrium model incorporating the effect of activity as activity-induced attraction and repulsion by imposing that the strength of these interactions depend on activity superlinearly. The model explains the cluster size dependence of activity obtained from simulations semi-quantitatively. Our predictions are verifiable in experiments on interacting synthetic colloidal microswimmers

preprint2015arXiv

Flow of colloidal solids and fluids through constrictions: dynamical density functional theory versus simulation

Using both dynamical density functional theory and particle-resolved Brownian dynamics simulations, we explore the flow of two-dimensional colloidal solids and fluids driven through a linear channel with a geometric constriction. The flow is generated by a constant external force acting on all colloids. The initial configuration is equilibrated in the absence of flow and then the external force is switched on instantaneously. Upon starting the flow, we observe four different scenarios: a complete blockade, a monotonic decay to a constant particle flux (typical for a fluid), a damped oscillatory behaviour in the particle flux, and a long-lived stop-and-go behaviour in the flow (typical for a solid). The dynamical density functional theory describes all four situations but predicts infinitely long undamped oscillations in the flow which are always damped in the simulations. We attribute the mechanisms of the underlying stop-and-go flow to symmetry conditions on the flowing solid. Our predictions are verifiable in real-space experiments on magnetic colloidal monolayers which are driven through structured microchannels and can be exploited to steer the flow throughput in microfluidics.

preprint2015arXiv

Getting drowned in a swirl: deformable bead-spring model microswimmers in external flow fields

Deformability is a central feature of many types of microswimmers, e.g. for artificially generated self-propelled droplets. Here, we analyze deformable bead-spring microswimmers in an externally imposed solvent flow field as simple theoretical model systems. We focus on their behavior in a circular swirl flow in two spatial dimensions. Linear (straight) two-bead swimmers are found to circle around the swirl with a slight drift to the outside with increasing activity. In contrast to that, we observe for triangular three-bead or square-like four-bead swimmers a tendency of being drawn into the swirl and finally getting drowned, although a radial inward component is absent in the flow field. During one cycle around the swirl, the self-propulsion direction of an active triangular or square-like swimmer remains almost constant, while their orbits become deformed exhibiting an ``egg-like'' shape. Over time, the swirl flow induces slight net rotations of these swimmer types, which leads to net rotations of the egg-shaped orbits. Interestingly, in certain cases, the orbital rotation changes sense when the swimmer approaches the flow singularity. Our predictions can be verified in real-space experiments on artificial microswimmers.

preprint2015arXiv

Glass transition of charged particles in two-dimensional confinement

The glass transition of mesoscopic charged particles in two-dimensional confinement is studied by mode-coupling theory. We consider two types of effective interactions between the particles, corresponding to two different models for the distribution of surrounding ions that are integrated out in coarse-grained descriptions. In the first model, a planar monolayer of charged particles is immersed in an unbounded isotropic bath of ions, giving rise to an isotropically screened Debye-Hückel- (Yukawa-) type effective interaction. The second, experimentally more relevant system is a monolayer of negatively charged particles that levitate atop a flat horizontal electrode, as frequently encountered in laboratory experiments with complex (dusty) plasmas. A steady plasma current towards the electrode gives rise to an anisotropic effective interaction potential between the particles, with an algebraically long-ranged in-plane decay. In a comprehensive parameter scan that covers the typical range of experimentally accessible plasma conditions, we calculate and compare the mode-coupling predictions for the glass transition in both kinds of systems.

preprint2015arXiv

Gravitaxis of asymmetric self-propelled colloidal particles

Many motile microorganisms adjust their swimming motion relative to the gravitational field and thus counteract sedimentation to the ground. This gravitactic behavior is often the result of an inhomogeneous mass distribution which aligns the microorganism similar to a buoy. However, it has been suggested that gravitaxis can also result from a geometric fore-rear asymmetry, typical for many self-propelling organisms. Despite several attempts, no conclusive evidence for such an asymmetry-induced gravitactic motion exists. Here, we study the motion of asymmetric self-propelled colloidal particles which have a homogeneous mass density and a well-defined shape. In experiments and by theoretical modeling we demonstrate that a shape anisotropy alone is sufficient to induce gravitactic motion with either preferential upward or downward swimming. In addition, also trochoid-like trajectories transversal to the direction of gravity are observed.

preprint2015arXiv

How does a flexible chain of active particles swell?

We study the swelling of a flexible linear chain composed of active particles by analytical theory and computer simulation. Three different situations are considered: a free chain, a chain confined to an external harmonic trap, and a chain dragged at one end. First we consider an ideal chain with harmonic springs and no excluded volume between the monomers. The Rouse model of polymers is generalized to the case of self-propelled monomers and solved analytically. The swelling, as characterized by the spatial extension of the chain, scales with the monomer number defining a Flory exponent $ν$ which is $ν=1/2, 0, 1$ in the three different situations. As a result, we find that activity does not change the Flory exponent but affects the prefactor of the scaling law. This can be quantitatively understood by mapping the system onto an equilibrium chain with a higher effective temperature such that the chain swells under an increase of the self-propulsion strength. We then use computer simulations to study the effect of self-avoidance on active polymer swelling. In the three different situations, the Flory exponent is now $ν= 3/4, 1/4, 1$ and again unchanged under self-propulsion. However, the chain extension behaves non-monotonic in the self-propulsion strength.

preprint2015arXiv

Motion of two micro-wedges in a turbulent bacterial bath

The motion of a pair of micro-wedges ("carriers") in a turbulent bacterial bath is explored using computer simulations with explicit modeling of the bacteria and experiments. The orientation of the two micro-wedges is parallel and fixed by an applied magnetic field but the translational coordinates can move freely as induced by the bacterial bath. As a result, two carriers of same orientation move such that their mutual distance decreases. Eventually the two carriers stack on each other with no intervening bacteria exhibiting a stable dynamical mode where the two micro-wedges follow each other with the same velocity. These findings are in qualitative agreement with experiments on two micro-wedges in a bacterial bath. Our results provide insight into understanding the self-assembly of many micro-wedges in an active bath.

preprint2015arXiv

Structural correlations in diffusiophoretic colloidal mixtures with nonreciprocal interactions

Nonreciprocal effective interaction forces can occur between mesoscopic particles in colloidal suspensions that are driven out of equilibrium. These forces violate Newton's third law actio=reactio on coarse-grained length and time scales. Here we explore the statistical mechanics of Brownian particles with nonreciprocal effective interactions. Our model system is a binary fluid mixture of spherically symmetric, diffusiophoretic mesoscopic particles, and we focus on the time-averaged particle pair- and triplet-correlation functions. Based on the many-body Smoluchowski equation we develop a microscopic statistical theory for the particle correlations and test it by computer simulations. For model systems in two and three spatial dimensions, we show that nonreciprocity induces distinct nonequilibrium pair correlations. Our predictions can be tested in experiments with chemotactic colloidal suspensions.

preprint2015arXiv

Swim pressure on walls with curves and corners

The concept of swim pressure quantifies the average force exerted by microswimmers on confining walls in non-equilibrium. Here we explore how the swim pressure depends on the wall curvature and on the presence of sharp corners in the wall. For active Brownian particles at high dilution, we present a coherent framework which describes the force and torque on passive particles of arbitrary shape, in the limit of large particles compared to the persistence length of the swimmer trajectories. The resulting forces can be used to derive, for example, the activity-induced depletion interaction between two disks, as well as to optimize the shape of a tracer particle for high swimming velocity. Our predictions are verifiable in experiments on passive obstacles exposed to a bath of bacteria or artificial microswimmers.

preprint2015arXiv

Symmetry-breaking in clogging for oppositely driven particles

The clogging behavior of a symmetric binary mixture of particles that are driven in opposite directions through constrictions is explored by Brownian dynamics simulations and theory. A dynamical state with a spontaneously broken symmetry occurs where one species is flowing and the other is blocked for a long time which can be tailored by the size of the constrictions. Moreover, we find self-organized oscillations in clogging and unclogging of the two species. Apart from statistical physics, our results are of relevance for fields like biology, chemistry, and crowd management, where ions, microparticles, pedestrians, or other particles are driven in opposite directions through constrictions.

preprint2015arXiv

Tailoring superelasticity of soft magnetic materials

Embedding magnetic colloidal particles in an elastic polymer matrix leads to smart soft materials that can reversibly be addressed from outside by external magnetic fields. We discover a pronounced nonlinear superelastic stress-strain behavior of such materials using numerical simulations. This behavior results from a combination of two stress-induced mechanisms: a detachment mechanism of embedded particle aggregates as well as a reorientation mechanism of magnetic moments. The superelastic regime can be reversibly tuned or even be switched on and off by external magnetic fields and thus be tailored during operation. Similarities to the superelastic behavior of shape-memory alloys suggest analogous applications, with the additional benefit of reversible switchability and a higher biocompatibility of soft materials.

preprint2015arXiv

Towards a scale-bridging description of ferrogels and magnetic elastomers

Ferrogels and magnetic elastomers differentiate themselves from other materials by their unique capability of reversibly changing shape and mechanical properties under the influence of an external magnetic field. A crucial issue in the study of these outstanding materials is the interaction between the mesoscopic magnetic particles and the polymer matrix in which they are embedded. Here we analyze interactions between two such particles connected by a polymer chain, a situation representative for particle-crosslinked magnetic gels. To make a first step towards a scale-bridging description of the materials, effective potentials for mesoscopic configurational changes are specified using microscopic input obtained from simulations. Furthermore, the impact of the presence of magnetic interactions on the probability distributions and thermodynamic quantities of the system is considered. The resulting mesoscopic model potentials can be used to economically model the system on the particle length scales. This first coarse-graining step is important to realize simplified but realistic scale-bridging models for these promising materials.

preprint2014arXiv

A deformable microswimmer in a swirl: capturing and scattering dynamics

Inspired by the classical Kepler and Rutherford problem, we investigate an analogous set-up in the context of active microswimmers: the behavior of a deformable microswimmer in a swirl flow. First we identify new steady bound states in the swirl flow and analyze their stability. Second we study the dynamics of a self-propelled swimmer heading towards the vortex center, and we observe the subsequent capturing and scattering dynamics. We distinguish between two major types of swimmers, those that tend to elongate perpendicularly to the propulsion direction and those that pursue a parallel elongation. While the first ones can get caught by the swirl, the second ones were always observed to be scattered, which proposes a promising escape strategy. This offers a route to design artificial microswimmers that show the desired behavior in complicated flow fields. It should be straightforward to verify our results in a corresponding quasi-two-dimensional experiment using self-propelled droplets on water surfaces.

preprint2014arXiv

Active colloidal suspensions: Clustering and phase behavior

We review recent experimental, numerical, and analytical results on active suspensions of self-propelled colloidal beads moving in (quasi) two dimensions. Active colloids form part of the larger theme of active matter, which is noted for the emergence of collective dynamic phenomena away from thermal equilibrium. Both in experiments and computer simulations, a separation into dense aggregates, i.e., clusters, and a dilute gas phase has been reported even when attractive interactions and an alignment mechanism are absent. Here, we describe three experimental setups, discuss the different propelling mechanisms, and summarize the evidence for phase separation. We then compare experimental observations with numerical studies based on a minimal model of colloidal swimmers. Finally, we review a mean-field approach derived from first principles, which provides a theoretical framework for the density instability causing the phase separation in active colloids.

preprint2014arXiv

Active crystals and their stability

A recently introduced active phase field crystal model describes the formation of ordered resting and traveling crystals in systems of self-propelled particles. Increasing the active drive, a resting crystal can be forced to perform collectively ordered migration as a single traveling object. We demonstrate here that these ordered migrating structures are linearly stable. In other words, during migration, the single crystalline texture together with the globally ordered collective motion is preserved even on large length scales. Furthermore, we consider self-propelled particles on a substrate that are surrounded by a thin fluid film. We find that in this case the resulting hydrodynamic interactions can destabilize the order.

preprint2014arXiv

Brownian motion and the hydrodynamic friction tensor for colloidal particles of complex shape

We synthesize colloidal particles with various anisotropic shapes and track their orientationally resolved Brownian trajectories using confocal microscopy. An analysis of appropriate short-time correlation functions provides direct access to the hydrodynamic friction tensor of the particles revealing nontrivial couplings between the translational and rotational degrees of freedom. The results are consistent with calculations of the hydrodynamic friction tensor in the low-Reynolds-number regime for the experimentally determined particle shapes.

preprint2014arXiv

Can the self-propulsion of anisotropic microswimmers be described by using forces and torques?

The self-propulsion of artificial and biological microswimmers (i.e., active colloidal particles) has often been modelled by using a force and a torque entering into the overdamped equations for the Brownian motion of passive particles. This seemingly contradicts the fact that a swimmer is force-free and torque-free, i.e., that the net force and torque on the particle vanish. Using different models for mechanical and diffusiophoretic self-propulsion, we demonstrate here that the equations of motion of microswimmers can be mapped onto those of passive particles with the shape-dependent grand resistance matrix and formally external effective forces and torques. This is consistent with experimental findings on the circular motion of artificial asymmetric microswimmers driven by self-diffusiophoresis. The concept of effective self-propulsion forces and torques significantly facilitates the understanding of the swimming paths, e.g., for a microswimmer under gravity. However, this concept has its limitations when the self-propulsion mechanism of a swimmer is disturbed either by another particle in its close vicinity or by interactions with obstacles, such as a wall.

preprint2014arXiv

Coupling between bulk- and surface chemistry in suspensions of charged colloids

The ionic composition and pair correlations in fluid phases of realistically salt-free charged colloidal sphere suspensions are calculated in the primitive model. We obtain the number densities of all ionic species in suspension, including low-molecular weight microions, and colloidal macroions with acidic surface groups, from a self-consistent solution of a coupled physicochemical set of nonlinear algebraic equations and non-mean-field liquid integral equations. Here, we study suspensions of colloidal spheres with sulfonate or silanol surface groups, suspended in demineralized water that is saturated with carbon dioxide under standard atmosphere. The only input required for our theoretical scheme are the acidic dissociation constants pKa, and effective sphere diameters of all involved ions. Our method allows for an ab initio calculation of colloidal bare and effective charges, at high numerical efficiency.

preprint2014arXiv

Glass-Transition Properties from Hard Spheres to Charged Point Particles

The glass transition is investigated in three dimensions for single and double Yukawa potentials for the full range of control parameters. For vanishing screening parameter, the limit of the one-component plasma is obtained; for large screening parameters and high coupling strengths, the glass-transition properties crossover to the hard-sphere system. Between the two limits, the entire transition diagram can be described by analytical functions. Different from other potentials, the glass-transition and melting lines for Yukawa potentials are found to follow shifted but otherwise identical curves in control-parameter space.

preprint2014arXiv

How does a thermal binary crystal break under shear?

When exposed to strong shearing, the particles in a crystal will rearrange and ultimately, the crystal will break by forming large nonaffine defects. Even for the initial stage of this process, only little effort has been devoted to the understanding of the breaking process on the scale of the individual particle size for thermalized mixed crystals. Here, we explore the shear-induced breaking for an equimolar two-dimensional binary model crystal with a high interaction asymmetry between the two different species such that the initial crystal has an intersecting square sublattice of the two constituents. Using Brownian dynamics computer simulations, we show that the combination of shear and thermal fluctuations leads to a characteristic hierarchical breaking scenario where initially, the more strongly coupled particles are thermally distorted, paving the way for the weakly coupled particles to escape from their cage. This in turn leads to mobile defects which may finally merge, proliferating a cascade of defects, which triggers the final breakage of the crystal. This scenario is in marked contrast to the breakage of one-component crystals close to melting. Moreover, we explore the orientational dependence of the initial shear direction relative to the crystal orientation and compare this to the usual melting scenario without shear. Our results are verifiable in real-space experiments of superparamagnetic colloidal mixtures at a pending air-water interface in an external magnetic field where the shear can be induced by an external laser field.

preprint2014arXiv

Magnetomechanical response of bilayered magnetic elastomers

Magnetic elastomers are appealing materials from an application point of view: they combine the mechanical softness and deformability of polymeric substances with the addressability by external magnetic fields. In this way, mechanical deformations can be reversibly induced and elastic moduli can be reversibly adjusted from outside. So far, mainly the behavior of single-component magnetic elastomers and ferrogels has been studied. Here, we go one step further and analyze the magnetoelastic response of a bilayered material composed of two different magnetic elastomers. It turns out that, under appropriate conditions, the bilayered magnetic elastomer can show a strongly amplified deformational response in comparison to a single-component material. Furthermore, a qualitatively opposite response can be obtained, i.e.\ a contraction along the magnetic field direction (as opposed to an elongation in the single-component case). We hope that our results will further stimulate experimental and theoretical investigations directly on bilayered magnetic elastomers, or, in a further hierarchical step, on bilayered units embedded in yet another polymeric matrix.

preprint2014arXiv

Mechanisms of carrier transport induced by a microswimmer bath

Recently, it was found that a wedgelike microparticle (referred to as "carrier") which is only allowed to translate but not to rotate exhibits a directed translational motion along the wedge cusp if it is exposed to a bath of microswimmers. Here we model this effect in detail by resolving the microswimmers explicitly using interaction models with different degrees of mutual alignment. Using computer simulations we study the impact of these interactions on the transport efficiency of V-shaped carrier. We show that the transport mechanisms itself strongly depends on the degree of alignment embodied in the modelling of the individual swimmer dynamics. For weak alignment, optimal carrier transport occurs in the turbulent microswimmer state and is induced by swirl depletion inside the carrier. For strong aligning interactions, optimal transport occurs already in the dilute regime and is mediated by a polar cloud of swimmers in the carrier wake pushing the wedge-particle forward. We also demonstrate that the optimal shape of the carrier leading to maximal transport speed depends on the kind of interaction model used.

preprint2014arXiv

Negative interfacial tension in phase-separated active suspensions

We study numerically a model for active suspensions of self-propelled repulsive particles, for which a stable phase separation into a dilute and a dense phase is observed. We exploit that for non-square boxes a stable "slab" configuration is reached, in which interfaces align with the shorter box edge. Evaluating a recent proposal for an intensive active swimming pressure, we demonstrate that the excess stress within the interface separating both phases is negative. The occurrence of a negative tension together with stable phase separation is a genuine non-equilibrium effect that is rationalized in terms of a positive stiffness, the estimate of which agrees excellently with the numerical data. Our results challenge effective thermodynamic descriptions and mappings of active suspensions onto passive pair potentials with attractions.

preprint2014arXiv

Structural control of elastic moduli in ferrogels and the importance of non-affine deformations

One of the central appealing properties of magnetic gels and elastomers is that their elastic moduli can reversibly be adjusted from outside by applying magnetic fields. The impact of the internal magnetic particle distribution on this effect has been outlined and analyzed theoretically. In most cases, however, affine sample deformations are studied and often regular particle arrangements are considered. Here we challenge these two major simplifications by a systematic approach using a minimal dipole-spring model. Starting from different regular lattices, we take into account increasingly randomized structures, until we finally investigate an irregular texture taken from a real experimental sample. On the one hand, we find that the elastic tunability qualitatively depends on the structural properties, here in two spatial dimensions. On the other hand, we demonstrate that the assumption of affine deformations leads to increasingly erroneous results the more realistic the particle distribution becomes. Understanding the consequences of the assumptions made in the modeling process is important on our way to support an improved design of these fascinating materials.

preprint2014arXiv

Swimming path statistics of an active Brownian particle with time-dependent self-propulsion

Typically, in the description of active Brownian particles, a constant effective propulsion force is assumed, which is then subjected to fluctuations in orientation and translation leading to a persistent random walk with an enlarged long-time diffusion coefficient. Here, we generalize previous results for the swimming path statistics to a time-dependent and thus in many situations more realistic propulsion which is a prescribed input. We analytically calculate both the noise-free and the noise-averaged trajectories for time-periodic propulsion under the action of an additional torque. In the deterministic case, such an oscillatory microswimmer moves on closed paths that can be highly more complicated than the commonly observed straight lines and circles. When exposed to random fluctuations, the mean trajectories turn out to be self-similar curves which bear the characteristics of their noise-free counterparts. Furthermore, we consider a propulsion force which scales in time $t$ as $\propto \! t^α$ (with $α=0,1,2,\ldots$) and analyze the resulting superdiffusive behaviour. Our predictions are verifiable for diffusiophoretic artificial microswimmers with prescribed propulsion protocols.

preprint2014arXiv

Transport powered by bacterial turbulence

We demonstrate that collective turbulent-like motion in a bacterial bath can power and steer directed transport of mesoscopic carriers through the suspension. In our experiments and simulations, a microwedge-like "bulldozer" draws energy from a bacterial bath of varied density. We obtain that a maximal transport speed is achieved in the turbulent state of the bacterial suspension. This apparent rectification of random motion of bacteria is caused by polar ordered bacteria inside the cusp region of the carrier, which is shielded from the outside turbulent fluctuations.

preprint2014arXiv

Tunable dynamic response of magnetic gels: impact of structural properties and magnetic fields

Ferrogels and magnetic elastomers feature mechanical properties that can be reversibly tuned from outside through magnetic fields. Here we concentrate on the question how their dynamic response can be adjusted. The influence of three factors on the dynamic behavior is demonstrated using appropriate minimal models: first, the orientational memory imprinted into one class of the materials during their synthesis; second, the structural arrangement of the magnetic particles in the materials; and third, the strength of an external magnetic field. To illustrate the latter point, structural data are extracted from a real experimental sample and analyzed. Understanding how internal structural properties and external influences impact the dominant dynamical properties helps to design materials that optimize the requested behavior.

preprint2014arXiv

Unusual swelling of a polymer in a bacterial bath

The equilibrium structure and dynamics of a single polymer chain in a thermal solvent is by now well-understood in terms of scaling laws. Here we consider a polymer in a bacterial bath, i.e. in a solvent consisting of active particles which bring in nonequilibrium fluctuations. Using computer simulations of a self-avoiding polymer chain in two dimensions which is exposed to a dilute bath of active particles, we show that the Flory-scaling exponent is unaffected by the bath activity provided the chain is very long. Conversely, for shorter chains, there is a nontrivial coupling between the bacteria intruding into the chain which may stiffen and expand the chain in a nonuniversal way. As a function of the molecular weight, the swelling first scales faster than described by the Flory exponent, then an unusual plateau-like behaviour is reached and finally a crossover to the universal Flory behaviour is observed. As a function of bacterial activity, the chain end-to-end distance exhibits a pronounced non-monotonicity. Moreover, the mean-square displacement of the center of mass of the chain shows a ballistic behaviour at intermediate times as induced by the active solvent. Our predictions are verifiable in two-dimensional bacterial suspensions and for colloidal model chains exposed to artificial colloidal microswimmers.

preprint2013arXiv

Circular motion of asymmetric self-propelling particles

Micron-sized self-propelled (active) particles can be considered as model systems for characterizing more complex biological organisms like swimming bacteria or motile cells. We produce asymmetric microswimmers by soft lithography and study their circular motion on a substrate and near channel boundaries. Our experimental observations are in full agreement with a theory of Brownian dynamics for asymmetric self-propelled particles, which couples their translational and orientational motion.

preprint2013arXiv

Compatibility waves drive crystal growth on patterned substrates

We explore the crystallization in a colloidal monolayer on a structured template starting from a few-particle nucleus. The competition between the substrate structure and that of the growing crystal induces a new crystal growth scenario. Unlike with the crystal growth in the bulk where a well-defined and connected crystal-fluid interface grows into the fluid, we identify a mechanism where a "compatibility wave" of the prescribed nucleus with the underlying substrate structure dictates the growth direction and efficiency. The growth process is strongly anisotropic and proceeds via transient island formation in front of an initial solid-fluid interface. We demonstrate the validity of this compatibility wave concept for a large class of substrate structures including a square-lattice and a quasicrystalline pattern. Dynamical density functional theory which provides a microscopic approach to the crystallization process is employed for colloidal hard spheres. Our predictions can be verified in experiments on confined colloids and also bear consequences for molecular crystal growth on structured substrates.

preprint2013arXiv

Differently Shaped Hard Body Colloids in Confinement: From passive to active particles

We review recent progress in the theoretical description of anisotropic hard colloidal particles. The shapes considered range from rods and dumbbells to rounded cubes, polyhedra and to biaxial particles with arbitrary shape. Our focus is on both static and dynamical density functional theory and on computer simulations. We describe recent results for the structure, dynamics and phase behaviour in the bulk and in various confining geometries, e.g. established by two parallel walls which reduce the dimensionality of the system to two dimensions. We also include recent theoretical modelling for active particles, which are autonomously driven by some intrinsic motor, and highlight their fascinating nonequilibrium dynamics and collective behaviour.

preprint2013arXiv

Dynamical clustering and phase separation in suspensions of self-propelled colloidal particles

We study experimentally and numerically a (quasi) two dimensional colloidal suspension of self-propelled spherical particles. The particles are carbon-coated Janus particles, which are propelled due to diffusiophoresis in a near-critical water-lutidine mixture. At low densities, we find that the driving stabilizes small clusters. At higher densities, the suspension undergoes a phase separation into large clusters and a dilute gas phase. The same qualitative behavior is observed in simulations of a minimal model for repulsive self-propelled particles lacking any alignment interactions. The observed behavior is rationalized in terms of a dynamical instability due to the self-trapping of self-propelled particles.

preprint2013arXiv

Dynamics of a deformable active particle under shear flow

The motion of a deformable active particle in linear shear flow is explored theoretically. Based on symmetry considerations, in two spatial dimensions, we propose coupled nonlinear dynamical equations for the particle position, velocity, deformation, and rotation. In our model, both, passive rotations induced by the shear flow as well as active spinning motions, are taken into account. Our equations reduce to known models in the two limits of vanishing shear flow and vanishing particle deformability. For varied shear rate and particle propulsion speed, we solve the equations numerically and obtain a manifold of different dynamical modes including active straight motion, periodic motions, motions on undulated cycloids, winding motions, as well as quasi-periodic and chaotic motions induced at high shear rates. The types of motion are distinguished by different characteristics in the real-space trajectories and in the dynamical behavior of the particle orientation and its deformation. Our predictions can be verified in experiments on self-propelled droplets exposed to a linear shear flow.

preprint2013arXiv

Effective Cahn-Hilliard equation for the phase separation of active Brownian particles

The kinetic separation of repulsive active Brownian particles into a dense and a dilute phase is analyzed using a systematic coarse-graining strategy. We derive an effective Cahn-Hilliard equation on large length and time scales, which implies that the separation process can be mapped onto that of passive particles. A lower density threshold for clustering is found, and using our approach we demonstrate that clustering first proceeds via a hysteretic nucleation scenario and above a higher threshold changes into a spinodal-like instability. Our results are in agreement with particle-resolved computer simulations and can be verified in experiments of artificial or biological microswimmers.

preprint2013arXiv

Highly asymmetric electrolytes in the primitive model: Hypernetted chain solution in arbitrary spatial dimensions

The pair-correlation functions for fluid ionic mixtures in arbitrary spatial dimensions are computed in hypernetted chain (HNC) approximation. In the primitive model, all ions are approximated as non-overlapping hyperspheres with Coulomb interactions. Our spectral HNC solver is based on a Fourier-Bessel transform introduced by Talman [J. Comput. Phys., 29, 35 (1978)], with logarithmically spaced computational grids. Numeric efficiency for arbitrary spatial dimensions is a commonly exploited virtue of this transform method. Here, we highlight another advantage of logarithmic grids, consisting in efficient sampling of pair-correlation functions for highly asymmetric ionic mixtures. For three-dimensional fluids, ion size- and charge-ratios larger than one thousand can be treated, corresponding to hitherto computationally not accessed micrometer-sized colloidal spheres in 1-1 electrolyte. Effective colloidal charge numbers are extracted from our primitive model results. For moderately large ion size- and charge-asymmetries, we present Molecular Dynamics simulation results that agree well with the approximate HNC pair correlations.

preprint2013arXiv

Introduction to colloidal dispersions in external fields

Progress in the research area of colloidal dispersions in external fields within the last years is reviewed. Colloidal dispersions play a pivotal role as model systems for phase transitions in classical statistical mechanics. In recent years the leading role of colloids to realize model systems has become evident not only for equilibrium situations but also far away from equilibrium. By using external fields (such as shear flow, electric, magnetic or laseroptical fields as well as confinement), a colloidal suspension can be brought into nonequilibrium in a controlled way. Various kinds of equilibrium and nonequilibrium phenomena explored by colloidal dispersions are described providing also a guide and summary to this special issue. Particular emphasis is put on the comparison of real-space experiments, computer simulations and statistical theories.

preprint2013arXiv

Microscopic theory for the phase separation of self-propelled repulsive disks

Motivated by recent experiments on colloidal suspensions, we study analytically and numerically a microscopic model for self-propelled particles lacking alignment interactions. In this model, even for purely repulsive interactions, a dynamical instability leading to phase separation has been reported. Starting from the many-body Smoluchowski equation, we develop a mean-field description based on a novel closure scheme and derive the effective hydrodynamic equations. We demonstrate that the microscopic origin of the instability is a force imbalance due to an anisotropic pair distribution leading to self-trapping. The phase diagram can be understood in terms of two quantities: a minimal drive and the force imbalance. At sufficiently high propulsion speeds there is a reentrance into the disordered fluid.

preprint2013arXiv

Non-Gaussian statistics for the motion of self-propelled Janus particles: experiment versus theory

Spherical Janus particles are one of the most prominent examples for active Brownian objects. Here, we study the diffusiophoretic motion of such microswimmers in experiment and in theory. Three stages are found: simple Brownian motion at short times, super-diffusion at intermediate times, and finally diffusive behavior again at long times. These three regimes observed in the experiments are compared with a theoretical model for the Langevin dynamics of self-propelled particles with coupled translational and rotational motion. Besides the mean square displacement also higher displacement moments are addressed. In particular, theoretical predictions regarding the non-Gaussian behavior of self-propelled particles are verified in the experiments. Furthermore, the full displacement probability distribution is analyzed, where in agreement with Brownian dynamics simulations either an extremely broadened peak or a pronounced double-peak structure is found depending on the experimental conditions.

preprint2013arXiv

Structure and dynamics of interfaces between two coexisting liquid crystalline phases

The phase-field-crystal model is used to access the structure and thermodynamics of interfaces between two coexisting liquid crystalline phases in two spatial dimensions. Depending on the model parameters there is a variety of possible coexistences between two liquid crystalline phases including a plastic triangular crystal (PTC). Here, we calculate numerically the profiles for the mean density and the nematic order tensor across the interface for isotropic-PTC and columnar-PTC respectively smectic A-PTC coexistence. As a general finding, the width of the interface with respect to the nematic order parameter characterizing the orientational order is larger than the width of the mean density interface. In approaching the interface from the PTC side, at first the mean density goes down and then the nematic order parameter follows. The relative shift of the two profiles can be larger than a full lattice constant of the plastic crystal. Finally, we also present numerical results for the dynamic relaxation of an initial order-parameter profile towards its equilibrium interfacial profile. Our predictions for the interfacial profiles can in principle be verified in real-space experiments of colloidal dispersions.

preprint2013arXiv

The effect of boundary adaptivity on hexagonal ordering and bistability in circularly confined quasi hard discs

The behaviour of materials under spatial confinement is sensitively dependent on the nature of the confining boundaries. In two dimensions, confinement within a hard circular boundary inhibits the hexagonal ordering observed in bulk systems at high density. Using colloidal experiments and Monte Carlo simulations, we investigate two model systems of quasi hard discs under circularly symmetric confinement. The first system employs an adaptive circular boundary, defined experimentally using holographic optical tweezers. We show that deformation of this boundary allows, and indeed is required for, hexagonal ordering in the confined system. The second system employs a circularly symmetric optical potential to confine particles without a physical boundary. We show that, in the absence of a curved wall, near perfect hexagonal ordering is possible. We propose that the degree to which hexagonal ordering is suppressed by a curved boundary is determined by the `strictness' of that wall.

preprint2013arXiv

Two-dimensional colloidal mixtures in magnetic and gravitational fields

This mini-review is concerned with two-dimensional colloidal mixtures exposed to various kinds of external fields. By a magnetic field perpendicular to the plane, dipole moments are induced in paramagnetic particles which give rise to repulsive interactions leading to complex crystalline alloys in the composition-asymmetry diagram. A quench in the magnetic field induces complex crystal nucleation scenarios. If exposed to a gravitational field, these mixtures exhibit a brazil--nut effect and show a boundary layering which is explained in terms of a depletion bubble picture. The latter persists for time-dependent gravity ("colloidal shaking"). Finally, we summarize crystallization effects when the second species is frozen in a disordered matrix which provides obstacles for the crystallizing component.

preprint2013arXiv

Two-dimensional melting under quenched disorder

We study the influence of quenched disorder on the two-dimensional melting behavior by using both video-microscopy of superparamagnetic colloidal particles and computer simulations of repulsive parallel dipoles. Quenched disorder is provided by pinning a fraction of the particles. We confirm the occurrence of the Kosterlitz-Thouless-Halperin-Nelson-Young scenario with an intermediate hexatic phase. While the fluid-hexatic transition remains largely unaffected by disorder, the hexatic-solid transition shifts towards lower temperatures for increasing disorder resulting in a significantly broadened stability range of the hexatic phase. In addition, we observe spatio-temporal critical(-like) fluctuations consistent with the continuous character of the phase transitions.

preprint2013arXiv

Vacancy diffusion in colloidal crystals as determined by dynamical density-functional theory and the phase-field-crystal model

A two-dimensional crystal of repulsive dipolar particles is studied in the vicinity of its melting transition by using Brownian dynamics computer simulation, dynamical density functional theory and phase-field crystal modelling. A vacancy is created by taking out a particle from an equilibrated crystal and the relaxation dynamics of the vacancy is followed by monitoring the time-dependent one-particle density. We find that the vacancy is quickly filled up by diffusive hopping of neighbouring particles towards the vacancy center. We examine the temperature dependence of the diffusion constant and find that it decreases with decreasing temperature in the simulations. This trend is reproduced by the dynamical density functional theory. Conversely, the phase field crystal calculations predict the opposite trend. Therefore, the phase-field model needs a temperature-dependent expression for the mobility to predict trends correctly.

preprint2012arXiv

Extended dynamical density functional theory for colloidal mixtures with temperature gradients

In the past decade, classical dynamical density functional theory (DDFT) has been developed and widely applied to the Brownian dynamics of interacting colloidal particles. One of the possible derivation routes of DDFT from the microscopic dynamics is via the Mori-Zwanzig-Forster projection operator technique with slowly varying variables such as the one-particle density. Here, we use the projection operator approach to extend DDFT into various directions: first, we generalize DDFT toward mixtures of $n$ different species of spherical colloidal particles. We show that there are in general nontrivial cross-coupling terms between the concentration fields and specify them explicitly for colloidal mixtures with pairwise hydrodynamic interactions. Secondly, we treat the energy density as an additional slow variable and derive formal expressions for an extended DDFT containing also the energy density. The latter approach can in principle be applied to colloidal dynamics in a nonzero temperature gradient. For the case without hydrodynamic interactions the diffusion tensor is diagonal, while thermodiffusion -- the dissipative cross-coupling term between energy density and concentration -- is nonzero in this limit. With finite hydrodynamic interactions also cross-diffusion coefficients assume a finite value. We demonstrate that our results for the extended DDFT contain the transport coefficients in the hydrodynamic limit (long wavelengths, low frequencies) as a special case.

preprint2012arXiv

Freezing of parallel hard cubes with rounded edges

The freezing transition in a classical three-dimensional system of parallel hard cubes with rounded edges is studied by computer simulation and fundamental-measure density functional theory. By switching the rounding parameter s from zero to one, one can smoothly interpolate between cubes with sharp edges and hard spheres. The equilibrium phase diagram of rounded parallel hard cubes is computed as a function of their volume fraction and the rounding parameter s. The second order freezing transition known for oriented cubes at s = 0 is found to be persistent up to s = 0.65. The fluid freezes into a simple-cubic crystal which exhibits a large vacancy concentration. Upon a further increase of s, the continuous freezing is replaced by a first-order transition into either a sheared simple cubic lattice or a deformed face-centered cubic lattice with two possible unit cells: body-centered orthorhombic or base-centered monoclinic. In principle, a system of parallel cubes could be realized in experiments on colloids using advanced synthesis techniques and a combination of external fields.

preprint2012arXiv

Meso-scale turbulence in living fluids

Turbulence is ubiquitous, from oceanic currents to small-scale biological and quantum systems. Self-sustained turbulent motion in microbial suspensions presents an intriguing example of collective dynamical behavior amongst the simplest forms of life, and is important for fluid mixing and molecular transport on the microscale. The mathematical characterization of turbulence phenomena in active non-equilibrium fluids proves even more difficult than for conventional liquids or gases. It is not known which features of turbulent phases in living matter are universal or system-specific, or which generalizations of the Navier-Stokes equations are able to describe them adequately. Here, we combine experiments, particle simulations, and continuum theory to identify the statistical properties of self-sustained meso-scale turbulence in active systems. To study how dimensionality and boundary conditions affect collective bacterial dynamics, we measured energy spectra and structure functions in dense Bacillus subtilis suspensions in quasi-2D and 3D geometries. Our experimental results for the bacterial flow statistics agree well with predictions from a minimal model for self-propelled rods, suggesting that at high concentrations the collective motion of the bacteria is dominated by short-range interactions. To provide a basis for future theoretical studies, we propose a minimal continuum model for incompressible bacterial flow. A detailed numerical analysis of the 2D case shows that this theory can reproduce many of the experimentally observed features of self-sustained active turbulence.

preprint2012arXiv

Microscopic approach to entropy production

It is a great challenge of nonequilibrium statistical mechanics to calculate entropy production within a microscopic theory. In the framework of linear irreversible thermodynamics, we combine the Mori-Zwanzig-Forster projection operator technique with the first and second law of thermodynamics to obtain microscopic expressions for the entropy production as well as for the transport equations of the entropy density and its time correlation function. We further present a microscopic derivation of a dissipation functional from which the dissipative dynamics of an extended dynamical density functional theory can be obtained in a formally elegant way.

preprint2012arXiv

Packing Confined Hard Spheres Denser with Adaptive Prism Phases

We show that hard spheres confined between two parallel hard plates pack denser with periodic adaptive prismatic structures which are composed of alternating prisms of spheres. The internal structure of the prisms adapts to the slit height which results in close packings for a range of plate separations, just above the distance where three intersecting square layers fit exactly between the plates. The adaptive prism phases are also observed in real-space experiments on confined sterically stabilized colloids and in Monte Carlo simulations at finite pressure.

preprint2012arXiv

Phase-field-crystal models for condensed matter dynamics on atomic length and diffusive time scales: an overview

Here, we review the basic concepts and applications of the phase-field-crystal (PFC) method, which is one of the latest simulation methodologies in materials science for problems, where atomic- and microscales are tightly coupled. The PFC method operates on atomic length and diffusive time scales, and thus constitutes a computationally efficient alternative to molecular simulation methods. Its intense development in materials science started fairly recently following the work by Elder et al. [Phys. Rev. Lett. 88 (2002), p. 245701]. Since these initial studies, dynamical density functional theory and thermodynamic concepts have been linked to the PFC approach to serve as further theoretical fundaments for the latter. In this review, we summarize these methodological development steps as well as the most important applications of the PFC method with a special focus on the interaction of development steps taken in hard and soft matter physics, respectively. Doing so, we hope to present today's state of the art in PFC modelling as well as the potential, which might still arise from this method in physics and materials science in the nearby future.

preprint2012arXiv

Tension and stiffness of the hard sphere crystal-fluid interface

A combination of fundamental measure density functional theory and Monte Carlo computer simulation is used to determine the orientation-resolved interfacial tension and stiffness for the equilibrium hard-sphere crystal-fluid interface. Microscopic density functional theory is in quantitative agreement with simulations and predicts a tension of 0.66 kT/σ^2 with a small anisotropy of about 0.025 kT and stiffnesses with e.g. 0.53 kT/σ^2 for the (001) orientation and 1.03 kT/σ^2 for the (111) orientation. Here kT is denoting the thermal energy and σthe hard sphere diameter. We compare our results with existing experimental findings.

preprint2012arXiv

Traveling and resting crystals in active systems

A microscopic field theory for crystallization in active systems is proposed which unifies the phase-field-crystal model of freezing with the Toner-Tu theory for self-propelled particles. A wealth of different active crystalline states are predicted and characterized. In particular, for increasing strength of self-propulsion, a transition from a resting crystal to a traveling crystalline state is found where the particles migrate collectively while keeping their crystalline order. Our predictions, which are verifiable in experiments and in particle-resolved computer simulations, provide a starting point for the design of new active materials.

preprint2011arXiv

Brownian dynamics of a self-propelled particle in shear flow

Brownian dynamics of a self-propelled particle in linear shear flow is studied analytically by solving the Langevin equation and in simulation. The particle has a constant propagation speed along a fluctuating orientation and is additionally subjected to a constant torque. In two spatial dimensions, the mean trajectory and the mean square displacement (MSD) are calculated as functions of time t analytically. In general, the mean trajectories are cycloids that are modified by finite temperature effects. With regard to the MSD different regimes are identified where the MSD scales with t^a with a = 0,1,2,3,4. In particular, an accelerated (a = 4) motion emerges if the particle is self-propelled along the gradient direction of the shear flow.

preprint2011arXiv

DDFT calibration and investigation of an anisotropic phase-field crystal model

The anisotropic phase-field crystal model recently proposed and used by Prieler et al. [J. Phys.: Condens. Matter 21, 464110 (2009)] is derived from microscopic density functional theory for anisotropic particles with fixed orientation. Further its morphology diagram is explored. In particular we investigated the influence of anisotropy and undercooling on the process of nucleation and microstructure formation from atomic to the microscale. To that end numerical simulations were performed varying those dimensionless parameters which represent anisotropy and undercooling in our anisotropic phase-field crystal (APFC) model. The results from these numerical simulations are summarized in terms of a morphology diagram of the stable state phase. These stable phases are also investigated with respect to their kinetics and characteristic morphological features.

preprint2011arXiv

Dynamical density functional theory for colloidal particles with arbitrary shape

Starting from the many-particle Smoluchowski equation, we derive dynamical density functional theory for Brownian particles with an arbitrary shape. Both passive and active (self-propelled) particles are considered. The resulting theory constitutes a microscopic framework to explore the collective dynamical behavior of biaxial particles in nonequilibrium. For spherical and uniaxial particles, earlier derived dynamical density functional theories are recovered as special cases. Our study is motivated by recent experimental progress in preparing colloidal particles with many different biaxial shapes.

preprint2011arXiv

Huge broadening of the crystal-fluid interface for sedimenting colloids

For sedimenting colloidal hard spheres, the propagation and broadening of the crystal-fluid interface is studied by Brownian dynamics computer simulations of an initially homogeneous sample. Two different types of interface broadenings are observed: the first occurs during growth and is correlated with the interface velocity, the second is concomitant with the splitting of the crystal-fluid interface into the crystal-amorphous and amorphous-liquid interfaces. The latter width is strongly peaked as a function of the gravitational driving strength with a huge amplitude relative to its equilibrium counterpart.

preprint2011arXiv

Microscopic and macroscopic theories for the dynamics of polar liquid crystals

We derive and analyze the dynamic equations for polar liquid crystals in two spatial dimensions in the framework of classical dynamical density functional theory (DDFT). Translational density variations, polarization, and quadrupolar order are used as order-parameter fields. The results are critically compared with those obtained using the macroscopic approach of time-dependent Ginzburg-Landau (GL) equations for the analogous order-parameter fields. We demonstrate that for both the microscopic DDFT and the macroscopic GL approach the resulting dissipative dynamics can be derived from a dissipation function. We obtain microscopic expressions for all diagonal contributions and for many of the cross-coupling terms emerging from a GL approach. Thus we establish a bridge between molecular correlations and macroscopic modeling for the dissipative dynamics of polar liquid crystals.

preprint2011arXiv

Polar liquid crystals in two spatial dimensions: the bridge from microscopic to macroscopic modeling

Two-dimensional polar liquid crystals have been discovered recently in monolayers of anisotropic molecules. Here, we provide a systematic theoretical description of liquid-crystalline phases for polar particles in two spatial dimensions. Starting from microscopic density functional theory, we derive a phase-field-crystal expression for the free-energy density which involves three local order-parameter fields, namely the translational density, the polarization, and the nematic order parameter. Various coupling terms between the order-parameter fields are obtained which are in line with macroscopic considerations. Since the coupling constants are brought into connection with the molecular correlations, we establish a bridge from microscopic to macroscopic modeling. Our theory provides a starting point for further numerical calculations of the stability of polar liquid-crystalline phases and is also relevant for modeling of microswimmers which are intrinsically polar.

preprint2011arXiv

Soft repulsive mixtures under gravity: brazil-nut effect, depletion bubbles, boundary layering, nonequilibrium shaking

A binary mixture of particles interacting via long-ranged repulsive forces is studied in gravity by computer simulation and theory. The more repulsive A-particles create a depletion zone of less repulsive B-particles around them reminiscent to a bubble. Applying Archimedes' principle effectively to this bubble, an A-particle can be lifted in a fluid background of B-particles. This "depletion bubble" mechanism explains and predicts a brazil-nut effect where the heavier A-particles float on top of the lighter B-particles. It also implies an effective attraction of an A-particle towards a hard container bottom wall which leads to boundary layering of A-particles. Additionally, we have studied a periodic inversion of gravity causing perpetual mutual penetration of the mixture in a slit geometry. In this nonequilibrium case of time-dependent gravity, the boundary layering persists. Our results are based on computer simulations and density functional theory of a two-dimensional binary mixture of colloidal repulsive dipoles. The predicted effects also occur for other long-ranged repulsive interactions and in three spatial dimensions. They are therefore verifiable in settling experiments on dipolar or charged colloidal mixtures as well as in charged granulates and dusty plasmas.

preprint2011arXiv

Stability of liquid crystalline phases in the phase-field-crystal model

The phase-field-crystal model for liquid crystals is solved numerically in two spatial dimensions. This model is formulated with three position-dependent order parameters, namely the reduced translational density, the local nematic order parameter, and the mean local direction of the orientations. The equilibrium free-energy functional involves local powers of the order parameters up to fourth order, gradients of the order parameters up to fourth order, and different couplings between the order parameters. The stable phases of the equilibrium free-energy functional are calculated for various coupling parameters. Among the stable liquid-crystalline states are the isotropic, nematic, columnar, smectic A, and plastic crystalline phases. The plastic crystals can have triangular, square, and honeycomb lattices and exhibit orientational patterns with a complex topology involving a sublattice with topological defects. Phase diagrams were obtained by numerical minimization of the free-energy functional. Their main features are qualitatively in line with much simpler one-mode approximations for the order parameters.

preprint2011arXiv

The self-propelled Brownian spinning top: dynamics of a biaxial swimmer at low Reynolds numbers

Recently, the Brownian dynamics of self-propelled (active) rod-like particles was explored to model the motion of colloidal microswimmers, catalytically-driven nanorods, and bacteria. Here, we generalize this description to biaxial particles with arbitrary shape and derive the corresponding Langevin equation for a self-propelled Brownian spinning top. The biaxial swimmer is exposed to a hydrodynamic Stokes friction force at low Reynolds numbers, to fluctuating random forces and torques as well as to an external and an internal (effective) force and torque. The latter quantities control its self-propulsion. Due to biaxiality and hydrodynamic translational-rotational coupling, the Langevin equation can only be solved numerically. In the special case of an orthotropic particle in the absence of external forces and torques, the noise-free (zero-temperature) trajectory is analytically found to be a circular helix. This trajectory is confirmed numerically to be more complex in the general case involving a transient irregular motion before ending up in a simple periodic motion. By contrast, if the external force vanishes, no transient regime is found and the particle moves on a superhelical trajectory. For orthotropic particles, the noise-averaged trajectory is a generalized concho-spiral. We furthermore study the reduction of the model to two spatial dimensions and classify the noise-free trajectories completely finding circles, straight lines with and without transients, as well as cycloids and arbitrary periodic trajectories.

preprint2010arXiv

Brownian motion of a self-propelled particle

Overdamped Brownian motion of a self-propelled particle is studied by solving the Langevin equation analytically. On top of translational and rotational diffusion, in the context of the presented model, the "active" particle is driven along its internal orientation axis. We calculate the first four moments of the probability distribution function for displacements as a function of time for a spherical particle with isotropic translational diffusion as well as for an anisotropic ellipsoidal particle. In both cases the translational and rotational motion is either unconfined or confined to one or two dimensions. A significant non-Gaussian behavior at finite times t is signalled by a non-vanishing kurtosis. To delimit the super-diffusive regime, which occurs at intermediate times, two time scales are identified. For certain model situations a characteristic t^3 behavior of the mean square displacement is observed. Comparing the dynamics of real and artificial microswimmers like bacteria or catalytically driven Janus particles to our analytical expressions reveals whether their motion is Brownian or not.

preprint2010arXiv

Chemotactic predator-prey dynamics

A discrete chemotactic predator-prey model is proposed in which the prey secrets a diffusing chemical which is sensed by the predator and vice versa. Two dynamical states corresponding to catching and escaping are identified and it is shown that steady hunting is unstable. For the escape process, the predator-prey distance is diffusive for short times but exhibits a transient subdiffusive behavior which scales as a power law $t^{1/3}$ with time $t$ and ultimately crosses over to diffusion again. This allows to classify the motility and dynamics of various predatory bacteria and phagocytes. In particular, there is a distinct region in the parameter space where they prove to be infallible predators.

preprint2010arXiv

Derivation of a three-dimensional phase-field-crystal model for liquid crystals from density functional theory

Using a generalized order parameter gradient expansion within density functional theory, we derive a phase-field-crystal model for liquid crystals composed by apolar particles in three spatial dimensions. Both the translational density and the orientational direction and ordering are included as order parameters. Different terms involving gradients in the order parameters in the resulting free energy functional are compared to the macroscopic Ginzburg-Landau approach as well as to the hydrodynamic description for liquid crystals. Our approach provides microscopic expressions for all prefactors in terms of the particle interactions. Our phase-field-crystal model generalizes the conventional phase-field-crystal model of spherical particles to orientational degrees of freedom and can be used as a starting point to explore phase transitions and interfaces for various liquid-crystalline phases.

preprint2010arXiv

Phase-field-crystal model for liquid crystals

Based on static and dynamical density functional theory, a phase-field-crystal model is derived which involves both the translational density and the orientational degree of ordering as well as a local director field. The model exhibits stable isotropic, nematic, smectic A, columnar, plastic crystalline and orientationally ordered crystalline phases. As far as the dynamics is concerned, the translational density is a conserved order parameter while the orientational ordering is non-conserved. The derived phase-field-crystal model can serve for efficient numerical investigations of various nonequilibrium situations in liquid crystals.

preprint2009arXiv

Dynamics of a microorganism moving by chemotaxis in its own secretion

The Brownian dynamics of a single microorganism coupled by chemotaxis to a diffusing concentration field which is secreted by the microorganism itself is studied by computer simulations in spatial dimensions $d=1,2,3$. Both cases of a chemoattractant and a chemorepellent are discussed. For a chemoattractant, we find a transient dynamical arrest until the microorganism diffuses for long times. For a chemorepellent, there is a transient ballistic motion in all dimensions and a long-time diffusion. These results are interpreted with the help of a theoretical analysis.

preprint2009arXiv

Non-Gaussian behaviour of a self-propelled particle on a substrate

The overdamped Brownian motion of a self-propelled particle which is driven by a projected internal force is studied by solving the Langevin equation analytically. The "active" particle under study is restricted to move along a linear channel. The direction of its internal force is orientationally diffusing on a unit circle in a plane perpendicular to the substrate. An additional time-dependent torque is acting on the internal force orientation. The model is relevant for active particles like catalytically driven Janus particles and bacteria moving on a substrate. Analytical results for the first four time-dependent displacement moments are presented and analysed for several special situations. For vanishing torque, there is a significant dynamical non-Gaussian behaviour at finite times t as signalled by a non-vanishing normalized kurtosis in the particle displacement which approaches zero for long time with a 1/t long-time tail.

preprint2009arXiv

Simulation Study of Sulfonate Cluster Swelling in Ionomers

We have performed simulations to study how increasing humidity affects the structure of Nafion-like ionomers under conditions of low sulfonate concentration and low humidity. At the onset of membrane hydration, the clusters split into smaller parts. These subsequently swell, but then maintain constant the number of sulfonates per cluster. We find that the distribution of water in low-sulfonate membranes depends strongly on the sulfonate concentration. For a relatively low sulfonate concentration, nearly all the side-chain terminal groups are within cluster formations, and the average water loading per cluster matches the water content of membrane. However, for a relatively higher sulfonate concentration the water-to-sulfonate ratio becomes non-uniform. The clusters become wetter, while the inter-cluster bridges become drier. We note the formation of unusual shells of water-rich material that surround the sulfonate clusters.

preprint2009arXiv

Statistics of polymer adsorption under shear flow

Using non-equilibrium Brownian dynamics computer simulations, we have investigated the steady state statistics of a polymer chain under three different shear environments: i) linear shear flow in the bulk (no walls), ii) shear vorticity normal to the adsorbing wall, iii) shear gradient normal to the adsorbing wall. The statistical distribution of the chain end-to-end distance and its orientational angles are calculated within our monomer-resolved computer simulations. Over a wide range of shear rates, this distribution can be mapped onto a simple theoretical finite-extensible-nonlinear-elastic dumbbell model with fitted anisotropic effective spring constants. The tails of the angular distribution functions are consistent with scaling predictions borrowed from the bulk dumbbell model. Finally, the frequency of the characteristic periodic tumbling motion has been investigated by simulation as well and was found to be sublinear with the shear rate for the three set-ups, which extends earlier results done in experiments and simulations for free and tethered polymer molecules without adsorption.

preprint2009arXiv

Ultra-fast quenching of binary colloidal suspensions in an external magnetic field

An ultra-fast quench is applied to binary mixtures of superparamagnetic colloidal particles confined at a two-dimensional water-air interface by a sudden increase of an external magnetic field. This quench realizes a virtually instantaneous cooling which is impossible in molecular systems. Using real-space experiments, the relaxation behavior after the quench is explored. Local crystallites with triangular and square symmetry are formed on different time scales and the correlation peak amplitude of the small particles evolves nonmonotonically in time in agreement with Brownian dynamics computer simulations.

preprint2008arXiv

Non-equilibrium melting of colloidal crystals in confinement

We report on a novel and flexible experiment to investigate the non-equilibrium melting behaviour of model crystals made from charged colloidal spheres. In a slit geometry polycrystalline material formed in a low salt region is driven by hydrostatic pressure up an evolving gradient in salt concentration and melts at large salt concentration. Depending on particle and initial salt concentration, driving velocity and the local salt concentration complex morphologic evolution is observed. Crystal-melt interface positions and the melting velocity are obtained quantitatively from time resolved Bragg- and polarization microscopic measurements. A simple theoretical model predicts the interface to first advance, then for balanced drift and melting velocities to become stationary at a salt concentration larger than the equilibrium melting concentration. It also describes the relaxation of the interface to its equilibrium position in a stationary gradient after stopping the drive in different manners. We further discuss the influence of the gradient strength on the resulting interface morphology and a shear induced morphologic transition from polycrystalline to oriented single crystalline material before melting.

preprint2008arXiv

Oscillatory driven colloidal binary mixtures: axial segregation versus laning

Using Brownian dynamics computer simulations we show that binary mixtures of colloids driven in opposite directions by an oscillating external field exhibit axial segregation in sheets perpendicular to the drive direction. The segregation effect is stable only in a finite window of oscillation frequencies and driving strengths and is taken over by lane formation in direction of the driving field if the driving force is increased. In the crossover regime, bands tilted relative to the drive direction are observed. Possible experiments to verify the axial segregation are discussed.