Source author record

Sabine H. L. Klapp

Sabine H. L. Klapp 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

45works
7topics
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

45 published item(s)

preprint2023arXiv

Clustering and flocking of repulsive chiral active particles with non-reciprocal couplings

Recently, non-reciprocal systems have become a focus of growing interest. Examples occur in soft and active matter, but also in engineered quantum materials and neural (brain) networks. Here, we investigate the impact of non-reciprocity on the collective behavior of a system of (dry) chiral active matter. Specifically, we consider a mixture of "circle swimmers" with steric interactions and non-reciprocal alignment couplings. Based on hydrodynamic equations which we derive from a set of Langevin equations, we explore the interplay of non-reciprocity, finite size, and chirality. We first consider, as a reference, one-species systems with reciprocal couplings. Based on a linear stability analysis and numerical simulations, we here observe three different types of collective behavior, that is, flocking, motility-induced phase separation, and a combination of both. Turning then to a non-reciprocal system, we find that non-reciprocity can turn otherwise stationary instabilities into oscillatory ones, affect the relative orientation of flocks, and, crucially, change the general type of instability. This illustrates the drastic impact of non-reciprocity and chirality on the emergent collective dynamics of active matter systems, with potentially far-reaching biological implications.

preprint2022arXiv

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

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

preprint2022arXiv

Ising-like critical behavior of vortex lattices in an active fluid

Turbulent vortex structures emerging in bacterial active fluids can be organized into regular vortex lattices by weak geometrical constraints such as obstacles. Here we show, using a continuum-theoretical approach, that the formation and destruction of these patterns exhibit features of a continuous second-order equilibrium phase transition, including long-range correlations, divergent susceptibility, and critical slowing down. The emerging vorticity field can be mapped onto a two-dimensional (2D) Ising model with antiferromagnetic nearest-neighbor interactions by coarse-graining. The resulting effective temperature is found to be proportional to the strength of the nonlinear advection in the continuum model.

preprint2022arXiv

Optimal Turbulent Transport in Microswimmer Suspensions

Microswimmer suspensions self-organize into complex spatio-temporal flow patterns, including vortex lattices and mesoscale turbulence. Here we explore the consequences for the motion of passive tracers, based on a continuum model for the microswimmer velocity field. We observe two qualitatively different regimes distinguished via the dimensionless Kubo number $K$. At advection strengths right above the transition to turbulence, the flow field evolves very slowly ($K \gg 1$) and the spatial vortex structures lead to dominant trapping effects. In contrast, deep in the turbulent state, much faster dynamics ($K \ll 1$) consistent with the so-called sweeping hypothesis leads to transport properties completely determined by the temporal correlations. In between ($K \approx 1$), we observe a regime of optimal transport, signaled by a maximum of the diffusion coefficient.

preprint2021arXiv

Correlation functions of non-Markovian systems out of equilibrium: Analytical expressions beyond single-exponential memory

This paper is concerned with correlation functions of stochastic systems with memory, a prominent example being a molecule or colloid moving through a complex (e.g., viscoelastic) fluid environment. Analytical investigations of such systems based on non-Markovian stochastic equations are notoriously difficult. A common approximation is that of a single-exponential memory, corresponding to the introduction of one auxiliary variable coupled to the Markovian dynamics of the main variable. As a generalization, we here investigate a class of "toy" models with altogether three degrees of freedom, giving rise to more complex forms of memory. Specifically, we consider, mainly on an analytical basis, the under- and overdamped motion of a colloidal particle coupled linearly to two auxiliary variables, where the coupling between variables can be either reciprocal or non-reciprocal. Projecting out the auxiliary variables, we obtain non-Markovian Langevin equations with friction kernels and colored noise, whose structure is similar to that of a generalized Langevin equation. For the present systems, however, the non-Markovian equations may violate the fluctuation-dissipation relation as well as detailed balance, indicating that the systems are out of equilibrium. We then study systematically the connection between the coupling topology of the underlying Markovian system and various autocorrelation functions.We demonstrate that already two auxiliary variables can generate surprisingly complex (e.g., non-monotonic or oscillatory) memory and correlation functions. Finally, we show that a minimal overdamped model with two auxiliary variables and suitable non-reciprocal coupling yields correlation functions resembling those describing hydrodynamic backflow in an optical trap.

preprint2020arXiv

Dynamical modes of sheared confined microscale matter

Based on (overdamped) Stokesian dynamics simulations and video microscopy experiments, we study the non equilibrium dynamics of a sheared colloidal cluster, which is confined to a two-dimensional disk. The experimental system is composed of a mixture of paramagnetic and non magnetic polystyrene particles, which are held in the disk by time shared optical tweezers. The paramagnetic particles are located at the center of the disk and are actuated by an external, rotating magnetic field that induces a magnetic torque. We identify two different steady states by monitoring the mean angular velocities per ring. The first one is characterized by rare slip events, where the inner rings momentarily depin from the outer ring, which is kept static by the set of optical traps. For the second state, we find a bistability of the mean angular velocities, which can be understood from the analysis of the slip events in the particle trajectories. We calculate the particle waiting- and jumping time distributions and estimate a time scale between slips, which is also reflected by a plateau in the mean squared azimuthal displacement. The dynamical transition is further reflected by the components of the stress tensor, revealing a shear-thinning behavior as well as shear stress overshoots. Finally, we briefly discuss the observed transition in the context of stochastic thermodynamics and how it may open future directions in this field.

preprint2020arXiv

Freezing and re-entrant melting of hard discs in a one-dimensional potential:\\ Predictions based on a pressure-balance equation

We investigate theoretically the freezing behaviour of a two-dimensional (2D) system of hard discs on a one-dimensional (1D) external potential (typically called laser-induced freezing). As shown by earlier theoretical and numerical studies, one observes freezing of the modulated liquid upon increase of the substrate potential amplitude, and re-entrant melting back into the modulated liquid when the substrate potential amplitude is increased even further. The purpose of our present work is to calculate the freezing and re-entrant melting phase diagram based on information from the bulk system. To this end, we employ an integrated pressure-balance equation derived from density functional theory [Phys. Rev. E \textbf{101}, 012609 (2020)]. Furthermore, we define a measure to quantify the influence of registration effects that qualitatively explain re-entrant melting. Despite severe approximations, the calculated phase diagram shows good agreement with the known phase diagram obtained by Monte Carlo simulations.

preprint2019arXiv

Freezing of a soft-core fluid in a one-dimensional potential: Predictions based on a pressure-balance equation

Using concepts from classical density functional theory (DFT) we investigate the freezing of a two-dimensional (2D) system of ultra-soft particles in a one-dimensional (1D) external potential; a phenomenon often called laser-induced freezing (LIF). In the first part of the paper, we present numerical results from free minimization of a mean-field density functional for a system of particles interacting via the GEM-4 potential. We show that the system does indeed display a LIF transition, although the interaction potential is markedly different from the cases studied before. We also show that one may consider the (suitably defined) effective density within the potential wells, $\barρ_{\text{eff}}$, as a control parameter of LIF, rather than the amplitude of the external potential as in the common LIF scenario. In the second part, we suggest a new theoretical description of the onset of LIF which bases on the pressure balance equation relating the pressure tensor and the external potential. Evaluating this equation for the modulated liquid phase at effective density $\barρ_{\text{eff}}$ and combining it with the (known) stability threshold of the corresponding bulk fluid, we can predict the critical effective density or, equivalently, the potential amplitude related to the onset of LIF. Our approach yields very good results for the model at hand, and it is transferable, in principle, to other model systems.

preprint2016arXiv

A minimal model for short-time diffusion in periodic potentials

We investigate the dynamics of a single, overdamped colloidal particle, which is driven by a constant force through a one-dimensional periodic potential. We focus on systems with large barrier heights where the lowest-order cumulants of the density field, that is, average position and the mean-squared displacement, show nontrivial (non-diffusive) short-time behavior characterized by the appearance of plateaus. We demonstrate that this "cage-like" dynamics can be well described by a discretized master equation model involving two states (related to two positions) within each potential valley. Non-trivial predictions of our approach include analytic expressions for the plateau heights and an estimate of the "de-caging time" obtained from the study of deviations from Gaussian behaviour. The simplicity of our approach means that it offers a minimal model to describe the short-time behavior of systems with hindered dynamics.

preprint2016arXiv

Angle-resolved effective potentials for disk-shaped molecules

We present an approach for calculating coarse-grained angle-resolved effective pair potentials for uniaxial molecules. For integrating out the intramolecular degrees of freedom we apply umbrella sampling and steered dynamics techniques in atomistically-resolved molecular dynamics (MD) computer simulations. Throughout this study we focus on disk-like molecules such as coronene. To develop the methods we focus on integrating out the van-der-Waals and intramolecular interactions, while electrostatic charge contributions are neglected. The resulting coarse-grained pair potential reveals a strong temperature and angle dependence. In the next step we fit the numerical data with various Gay-Berne-like potentials to be used in more efficient simulations on larger scales. The quality of the resulting coarse-grained results is evaluated by comparing their pair and many-body structure as well as some thermodynamic quantities self-consistently to the outcome of atomistic MD simulations of many-particle systems. We find that angle-resolved potentials are essential not only to accurately describe crystal structures but also for fluid systems where simple isotropic potentials start to fail already for low to moderate packing fractions. Further, in describing these states it is crucial to take into account the pronounced temperature dependence arising in selected pair configurations due to bending fluctuations.

preprint2016arXiv

Coarse-grained electrostatic interactions of coronene: Towards the crystalline phase

In this article we present and compare two different, coarse-grained approaches to model electrostatic interactions of disc-shaped aromatic molecules, specifically coronene. Our study builds on previous work [J. Chem. Phys. 141, 214110 (2014)] where we proposed, based on a systematic coarse-graining procedure starting from the atomistic level, an anisotropic effective (Gay-Berne-like) potential capable of describing van-der-Waals contributions to the interaction energy. To take into account electrostatics, we introduce, first, a linear quadrupole moment along the symmetry axis of the coronene disc. The second approach takes into account the fact that the partial charges within the molecules are distributed in a ring-like fashion. We then reparametrize the effective Gay-Berne-like potential such that it matches, at short distances, the ring-ring potential. To investigate the validity of these two approaches, we perform many-particle Molecular Dynamics (MD) simulations, focusing on the crystalline phase (karpatite) where electrostatic interaction effects are expected to be particularly relevant for the formation of tilted stacked columns. Specifically, we investigate various structural parameters as well as the melting transition. We find that the second approach yields consistent results with those from experiments despite the fact that the underlying potential decays with the wrong distance dependence at large molecule separations. Our strategy can be transferred to a broader class of molecules, such as benzene or hexabenzocoronene.

preprint2016arXiv

Collective dynamics of dipolar and multipolar colloids: from passive to active systems

This article reviews recent research on the collective dynamical behavior of colloids with dipolar or multipolar interactions. Indeed, whereas equilibrium structures and static self-assembly of such systems are now rather well understood, the past years have seen an explosion of interest in understanding dynamicals aspects, from the relaxation dynamics of strongly correlated dipolar networks over systems driven by time-dependent, electric or magnetic fields, to pattern formation and dynamical control of active, self-propelled systems. Unraveling the underlying mechanisms is crucial for a deeper understanding of self-assembly in and out of equilibrium and the use of such particles as functional devices. At the same time, the complex dynamics of dipolar colloids poses challenging physical questions and puts forward their role as model systems for nonlinear behavior in condensed matter physics. Here we attempt to give an overview of these developments, with an emphasis on theoretical and simulation studies.

preprint2016arXiv

Depinning and heterogeneous dynamics of colloidal crystal layers under shear flow

Using Brownian dynamics (BD) simulations and an analytical approach we investigate the shear-induced, nonequilibrium dynamics of dense colloidal suspensions confined to a narrow slit-pore. Focusing on situations where the colloids arrange in well-defined layers with solidlike in-plane structure, the confined films display complex, nonlinear behavior such as collective depinning and local transport via density excitations. These phenomena are reminiscent of colloidal monolayers driven over a periodic substrate potential. In order to deepen this connection, we present an effective model which maps the dynamics of the shear-driven colloidal layers to the motion of a single particle driven over an effective substrate potential. This model allows to estimate the critical shear rate of the depinning transition based on the equilibrium configuration, revealing the impact of important parameters such as the slit-pore width and the interaction strength. We then turn to heterogeneous systems where a layer of small colloids is sheared with respect to bottom layers of large particles. For these incommensurate systems we find that the particle transport is dominated by density excitations resembling the so-called "kink" solutions of the Frenkel-Kontorova (FK) model. In contrast to the FK model, however, the corresponding "antikinks" do not move.

preprint2016arXiv

Equilibrium structures of anisometric, quadrupolar particles confined to a monolayer

We investigate the structural properties of a two-dimensional system of ellipsoidal particles carrying a linear quadrupole moment in their center. These particles represent a simple model for a variety of uncharged, non-polar conjugated organic molecules. Using optimization tools based on ideas of Evolutionary Algorithms, we first examine the ground state structures as we vary the aspect ratio of the particles and the pressure. Interestingly, we find, besides the intuitively expected T-like configurations, a variety of complex structures, characterized with up to three different particle orientations. In an effort to explore the impact of thermal fluctuations, we perform constant-pressure Molecular Dynamics simulations within a range of rather low temperatures. We observe that ground state structures formed by particles with a large aspect ratio are in particular suited to withstand fluctuations up to rather high temperatures. Our comprehensive investigations allow for a deeper understanding of molecular or colloidal monolayer arrangements under the influence of a typical electrostatic interaction on a coarse-grained level.

preprint2016arXiv

Non-equilibrium surface growth in a hybrid inorganic-organic system

Using kinetic Monte Carlo simulations, we show that molecular morphologies found in non-equilibrium growth can be strongly different from those at equilibrium. We study the prototypical hybrid inorganic-organic system 6P on ZnO(10-10) during thin film adsorption, and find a wealth of phenomena including re-entrant growth, a critical adsorption rate and observables that are non-monotonous with the adsorption rate. We identify the transition from lying to standing molecules with a critical cluster size and discuss the competition of time scales during growth in terms of a rate equation approach. Our results form a basis for understanding and predicting collective orientational ordering during growth in hybrid material systems.

preprint2016arXiv

Shear banding in nematogenic fluids with oscillating orientational dynamics

We investigate the occurrence of shear banding in nematogenic fluids under planar Couette flow, based on mesoscopic dynamical equations for the orientational order parameter and the shear stress. We focus on parameter values where the sheared homogeneous system exhibits regular oscillatory orientational dynamics, whereas the equilibrium system is either isotropic (albeit close to the isotropic--nematic transition) or deep in its nematic phase. The numerical calculations are restricted to spatial variations in shear gradient direction. We find several new types of shear banded states characterized by regions with regular oscillatory orientational dynamics. In all cases shear banding is accompanied by a non--monotonicity of the flow curve of the homogeneous system; however, only in the case of the initially isotropic system this curve has the typical $S$--like shape. We also analyze the influence of different orientational boundary conditions and of the spatial correlation length.

preprint2015arXiv

A scale-bridging modeling approach for anisotropic organic molecules at patterned semiconductor surfaces

Hybrid systems consisting of organic molecules at inorganic semiconductor surfaces are gaining increasing importance as thin film devices for optoelectronics. The efficiency of such devices strongly depends on the collective behavior of the adsorbed molecules. In the present paper we propose a novel, coarse-grained model addressing the condensed phases of a representative hybrid system, that is, para-sexiphenyl (6P) at zinc-oxide (ZnO). Within our model, intermolecular interactions are repre- sented via a Gay-Berne potential (describing steric and van-der-Waals interactions) combined with the electrostatic potential between two linear quadrupoles. Similarly, the molecule-substrate interactions include a coupling between a linear molecular quadrupole to the electric field generated by the line charges characterizing ZnO(10-10). To validate our approach, we perform equilibrium Monte Carlo simulations, where the lateral positions are fixed to a 2D lattice, while the rotational degrees of freedom are continuous. We use these simulations to investigate orientational ordering in the condensed state. We reproduce various experimentally observed features such as the alignment of individual molecules with the line charges on the surface, the formation of a standing uniaxial phase with a herringbone structure, as well as the formation of a lying nematic phase.

preprint2015arXiv

Binary mixtures of rod-like colloids under shear: microscopically-based equilibrium theory and order--parameter dynamics

This paper is concerned with the dynamics of a binary mixture of rod--like, repulsive colloidal particles driven out of equilibrium by means of a steady shear flow (Couette geometry). To this end we first derive, starting from a microscopic density functional in Parsons--Lee approximation, a mesoscopic free energy functional whose main variables are the orientational order parameter tensors. Based on this mesoscopic functional we then explore the stability of isotropic and nematic equilibrium phases in terms of composition and rod lengths. Second, by combining the equilibrium theory with the Doi--Hess approach for the order parameter dynamics under shear, we investigate the orientational dynamics of binary mixtures for a range of shear rates and coupling parameters. We find a variety of dynamical states, including synchronized oscillatory states of the two components, but also symmetry breaking behavior where the components display different in--plane oscillatory states.

preprint2015arXiv

Enhancement of mobility in an interacting colloidal system under feedback control

Feedback control schemes are a promising way to manipulate transport properties of driven colloidal suspensions. In the present article we suggest a feedback scheme to enhance the collective transport of colloidal particles with repulsive interactions through a one-dimensional tilted washboard potential. The control is modelled by a harmonic confining potential, mimicking an optical "trap", with the center of this trap moving with the (instantaneous) mean particle position. Our theoretical analysis is based on the Smoluchowski equation combined with Dynamical Density Functional Theory (DDFT) for systems with hard-core or ultra-soft (Gaussian) interactions. For either type of interaction we find that the feedback control can lead to an enhancement of the mobility by several orders of magnitude relative to the uncontrolled case. The largest effects occur for intermediate stiffness of the trap and large particle numbers. Moreover, in some regions of the parameter space the feedback control induces oscillations of the mean velocity. Finally, we show that the enhancement of mobility is robust against a small time delay in implementing the feedback control.

preprint2015arXiv

Generic model for tunable colloidal aggregation in multidirectional fields

Based on Brownian Dynamics computer simulations in two dimensions we investigate aggregation scenarios of colloidal particles with directional interactions induced by multiple external fields. To this end we propose a model which allows continuous change in the particle interactions from point-dipole-like to patchy-like (with four patches). We show that, as a result of this change, the non-equilibrium aggregation occurring at low densities and temperatures transforms from conventional diffusion-limited cluster aggregation (DLCA) to slippery DLCA involving rotating bonds; this is accompanied by a pronounced change of the underlying lattice structure of the aggregates from square-like to hexagonal ordering. Increasing the temperature we find a transformation to a fluid phase, consistent with results of a simple mean-field density functional theory.

preprint2015arXiv

Hydrodynamic length-scale selection and effective viscosity in microswimmer suspensions

A universal characteristic of mesoscale turbulence in active suspensions is the emergence of a typical vortex length scale, distinctly different from the scale-invariance of turbulent high-Reynolds number flows. Collective length-scale selection has been observed in bacterial fluids, endothelial tissue and active colloides, yet the physical origins of this phenomenon remain elusive. Here, we systematically derive an effective fourth-order field theory from a generic microscopic model that allows us to predict the typical vortex size in microswimmer suspensions. Building on a self-consistent closure condition, the derivation shows that the vortex length scale is determined by the competition between local alignment forces and intermediate-range hydrodynamic interactions. Vortex structures found in simulations of the theory agree with recent measurements in Bacillus subtilis suspensions. Moreover, our approach correctly predicts an effective viscosity enhancement (reduction), as reported experimentally for puller (pusher) microorganisms.

preprint2015arXiv

Lane formation in a system of dipolar microswimmers

Using Brownian Dynamics (BD) simulations we investigate the non-equilibrium structure formation of a two-dimensional (2D) binary system of dipolar colloids propelling in opposite directions. Despite of a pronounced tendency for chain formation, the system displays a transition towards a laned state reminiscent of lane formation in systems with isotropic repulsive interactions. However, the anisotropic dipolar interactions induce novel features: First, the lanes have themselves a complex internal structure characterized by chains or clusters. Second, laning occurs only in a window of interaction strengths. We interprete our findings by a phase separation process and simple force balance arguments.

preprint2015arXiv

Particle-resolved dynamics during multilayer growth of C$_{60}$

Using large-scale kinetic Monte-Carlo (KMC) simulations, we investigate the non-equilibrium surface growth of the fullerene C$_{60}$. Recently, we have presented a self-consistent set of energy barriers that describes the nucleation and multilayer growth of C$_{60}$ for different temperatures and adsorption rates in quantitative agreement with experiments [Bommel et al., Nat. Comm. 5, 5388 (2014)]. We found that C$_{60}$ displays lateral diffusion resembling colloidal systems, however it has to overcome an atom-like energetic step-edge barrier for interlayer diffusion. Here, we focus on the particle-resolved dynamics, and the interplay between surface morphology and particle dynamics during growth. Comparing C$_{60}$ growth with an atom-like system, we find significant differences in the evolution of the surface morphology, as well as the single-particle dynamics on the growing material landscape. By correlating the mean-squared-displacement of particles with their current neighborhood, we can identify the influence of the different time scales that compete during growth and can pinpoint the differences between the two systems.

preprint2015arXiv

Self-assembly of three-dimensional ensembles of magnetic particles with laterally shifted dipoles

We consider a model of colloidal spherical particles carrying a permanent dipole moment which is laterally shifted out of the particles' geometrical centres, i.e. the dipole vector is oriented perpendicular to the radius vector of the particles. Varying the shift $δ$ from the centre, we analyze ground state structures for two, three and four hard spheres, using a simulated annealing procedure. We also compare to earlier ground state results. We then consider a bulk system at finite temperatures and different densities. Using Molecular Dynamics simulations, we examine the equilibrium self-assembly properties for several shifts. Our results show that the shift of the dipole moment has a crucial impact on both, the ground state configurations as well as the self-assembled structures at finite temperatures.

preprint2015arXiv

Spontaneous ordering of magnetic particles in liquid crystals: From chains to biaxial lamellae

Using Monte Carlo (MC) computer simulations we explore the self-assembly and ordering behavior of a hybrid, soft magnetic system consisting of small magnetic nano-spheres in a liquid-crystalline (LC) matrix. Inspired by recent experiments with colloidal rod-matrices we focus on conditions where the sphere and rod diameters are comparable. Already in the absence of a magnetic field, the nematic ordering of the LC can stabilize formation of magnetic chains along the nematic or smectic director, yielding a state with local (yet no macroscopic) magnetic order. The chains, in turn, increase the overall nematic order, reflecting the complex interplay of the structure formation of the two components. Increasing the sphere diameter the spontaneous uniaxial ordering is replaced by biaxial lamellar morphologies characterized by alternating layers of rods and magnetic chains oriented perpendicular to the rod's director. These novel ordering scenarios at zero field suggest complex response of the resulting hybrid to external stimuli such as magnetic fields and shear forces.

preprint2015arXiv

Tunable structures of mixtures of magnetic particles in liquid-crystalline matrices

We investigate the self-organization of a binary mixture of similar sized rods and dipolar soft spheres by means of Monte-Carlo simulations. We model the interparticle interactions by employing anisotropic Gay-Berne, dipolar and soft-sphere interactions. In the limit of vanishing magnetic moments we obtain a variety of fully miscible liquid crystalline phases including nematic, smectic and lamellar phases. For the magnetic mixture, we find that the liquid crystalline matrix supports the formation of orientationally ordered ferromagnetic chains. Depending on the relative size of the species the chains align parallel or perpendicular to the director of the rods forming uniaxial or biaxial nematic, smectic and lamellar phases. As an exemplary external perturbation we apply a homogeneous magnetic field causing uniaxial or biaxial ordering to an otherwise isotropic state.

preprint2014arXiv

Delay-induced transport in a rocking ratchet under feedback control

Based on the Fokker-Planck equation we investigate the transport of an overdamped colloidal particle in a static, asymmetric periodic potential supplemented by a time-dependent, delayed feedback force, $F_{\mathrm{fc}}$. For a given time $t$, $F_\mathrm{fc}$ depends on the status of the system at a previous time $t-τ_\mathrm{D}$, with $τ_\mathrm{D}$ being a delay time, specifically on the delayed mean particle displacement (relative to some "switching position"). For non-zero delay times $F_{\mathrm{fc}}(t)$ develops nearly regular oscillations generating a net current in the system. Depending on the switching position, this current is nearly as large or even larger than that in a conventional open-loop rocking ratchet. We also investigate thermodynamic properties of the delayed non-equilibrium system and we suggest an underlying Langevin equation which reproduces the Fokker-Planck results.

preprint2014arXiv

Manipulating shear-induced non-equilibrium transitions by feedback control

Using Brownian Dynamics (BD) simulations we investigate non-equilibrium transitions of sheared colloidal films under controlled shear stress $σ_{\mathrm{xz}}$. In our approach the shear rate $\dotγ$ is a dynamical variable, which relaxes on a timescale $τ_c$ such that the instantaneous, configuration-dependent stress $σ_{\mathrm{xz}}(t)$ approaches a pre-imposed value. Investigating the dynamics under this "feedback-control" scheme we find unique behavior in regions where the flow curve $σ_{\mathrm{xz}}(\dotγ)$ of the uncontrolled system is monotonic. However, in non-monotonic regions our method allows to {\em select} between dynamical states characterized by different in-plane structure and viscosities. Indeed, the final state strongly depends on $τ_c$ relative to an {\em intrinsic} relaxation time of the uncontrolled system. The critical values of $τ_c$ are estimated on the basis of a simple model.

preprint2014arXiv

Novel structure formation of a phase separating colloidal fluid in a ratchet potential

Based on Dynamical Density Functional Theory (DDFT) we investigate a binary mixture of interacting Brownian particles driven over a substrate via a one-dimensional ratchet potential. The particles are modeled as soft spheres where one component carries a classical Heisenberg spin. In the absence of a substrate field, the system undergoes a first-order fluid-fluid demixing transition driven by the spin-spin interaction. We demonstrate that the interplay between the intrinsic spinodal decomposition and time-dependent external forces leads to a novel dynamical instability where stripes against the symmetry of the external potential form. This structural transition is observed for a broad range of parameters related to the ratchet potential. Moreover, we find intriguing effects for the particle transport.

preprint2014arXiv

On the waiting time distribution for continuous stochastic systems

The waiting time distribution (WTD) is a common tool for analysing discrete stochastic processes in classical and quantum systems. However, there are many physical examples where the dynamics is continuous and only approximately discrete, or where it is favourable to discuss the dynamics on a discretized and a continuous level in parallel. An example is the hindered motion of particles through potential landscapes with barriers. In the present paper we propose a consistent generalisation of the WTD from the discrete case to situations where the particles perform continuous barrier-crossing characterised by a finite duration. To this end, we introduce a recipe to calculate the WTD from the Fokker-Planck/Smoluchowski equation. In contrast to the closely related first passage time distribution (FPTD), which is frequently used to describe continuous processes, the WTD contains information about the direction of motion. As an application, we consider the paradigmatic example of an overdamped particle diffusing through a washboard potential. To verify the approach and to elucidate its numerical implications, we compare the WTD defined via the Smoluchowski equation with data from direct simulation of the underlying Langevin equation and find full consistency provided that the jumps in the Langevin approach are defined properly. Moreover, for sufficiently large energy barriers, the WTD defined via the Smoluchowski equation becomes consistent with that resulting from the analytical solution of a (two-state) master equation model for the short-time dynamics developed previous by us [PRE 86, 061135 (2012)]. Thus, our approach "interpolates" between these two types of stochastic motion. We illustrate our approach for both symmetric systems and systems under constant force.

preprint2014arXiv

Spinodal decomposition of a binary magnetic fluid confined to a surface

In our previous work [J. Chem. Phys. \textbf{136}, 024502 (2012)], we reported a demixing phase transition of a two-dimensional, binary Heisenberg fluid mixture driven by the ferromagnetic interactions of the magnetic species. Here, we present a theoretical study for the \textit{time-dependent} coarsening occuring within the two-phase region in the density-concentration plane, also known as spinodal decomposition. Our investigations are based on Dynamical Density Functional Theory (DDFT). The particles in the mixture are modelled as Gaussian soft spheres on a two-dimensional surface, where one component carries a classical spin of Heisenberg type. To investigate the two-phase region, we first present a linear stability analysis with respect to small, harmonic density perturbations. Second, to capture non-linear effects, we calculate time-dependent structure factors by combining DDFT with Percus' test particle method. For the growth of the average domain size $l$ during spinodal decomposition with time $t$, we observe a power-law behavior $l\propto t^{δ_α}$ with $δ_m\simeq 0.333$ for the magnetic species and $δ_n\simeq 0.323$ for the non-magnetic species.

preprint2013arXiv

Feedback control of flow alignment in sheared liquid crystals

Based on a continuum theory, we investigate the manipulation of the non-equilibrium behavior of a sheared liquid crystal via closed-loop feedback control. Our goal is to stabilize a specific dynamical state, that is, the stationary "flow-alignment", under conditions where the uncontrolled system displays oscillatory director dynamics with in-plane symmetry. To this end we employ time-delayed feedback control (TDFC), where the equation of motion for the ith component, q_i(t), of the order parameter tensor is supplemented by a control term involving the difference q_i(t)-q_i(t-τ). In this diagonal scheme, τis the delay time. We demonstrate that the TDFC method successfully stabilizes flow alignment for suitable values of the control strength, K, and τ; these values are determined by solving an exact eigenvalue equation. Moreover, our results show that only small values of K are needed when the system is sheared from an isotropic equilibrium state, contrary to the case where the equilibrium state is nematic.

preprint2013arXiv

Non-equilibrium dynamics of a confined colloidal bilayer in planar shear flow

Using Brownian dynamics (BD) simulations we investigate the impact of shear flow on structural and dynamical properties of a system of charged colloids confined to a narrow slit pore. Our model consists of spherical microions interacting through a Derjaguin-Landau-Verwey-Overbeek (DLVO) and a soft-sphere potential. The DLVO parameters were chosen according to a system of moderately charged silica particles (with valence Z$\sim$35) in a solvent of low ionic strength. At the confinement conditions considered, the colloids form two well-pronounced layers. In the present study we investigate shear-induced transitions of the translational order and dynamics in the layers, including a discussion of the translational diffusion. In particular, we show that diffusion in the shear-melted state can be described by an analytical model involving a single shear-driven particle in a harmonic trap. We also explore the emergence of zig-zag motion characterized by spatio-temporal oscillations of the particles in the layers in the vorticity direction. Similar behavior has been recently observed in experiments of much thicker colloidal films.

preprint2013arXiv

Oscillatory motion of sheared nanorods beyond the nematic phase

We study the role of the control parameter triggering nematic order (temperature or concentration) on the dynamical behavior of a system of nanorods under shear. Our study is based on a set of mesoscopic equations of motion for the components of the tensorial orientational order parameter. We investigating these equations via a systematic bifurcation analysis based on a numerical continuation technique, focusing on spatially homogeneous states. Exploring a wide range of parameters we find, unexpectedly, that states with oscillatory motion can exist even under conditions where the equilibrium system is isotropic. These oscillatory states are characterized by wagging motion of the paranematic director, and they occur if the tumbling parameter is sufficiently small. We also present full non-equilibrium phase diagrams, in the plane spanned by the concentration and the shear rate.

preprint2012arXiv

Critical exponents of colloid particles in bulk and confinement

Using grand canonical Monte Carlo simulations, we investigate the percolation behavior of a square well fluid with an ultra-short range of attraction in three dimension (3D) and in confined geometry. The latter is defined through two parallel and structureless walls (slit-pore). We focus on temperatures above the critical temperature of the (metastable) condensation transition of the 3D system. Investigating a broad range of systems sizes, we first determine the percolation thresholds, i. e., the critical packing fraction for percolation $η_{c}$. For the slit-pore systems, $η_{c}$ is found to vary with the wall separation $L_{z}$ in a continuous but non-monotonic way, $η_{c}(L_{z}\rightarrow\infty)=η_{c}^{\text{3D}}$. We also report results for critical exponents of the percolation transition, specifically, the exponent $ν$ of the correlation length $ξ$ and the two fisher exponents $τ$ and $σ$ of the cluster-size distribution. These exponents are obtained from a finite-size analysis involving the cluster-size distribution and the radii of gyration distribution at the percolation threshold. Within the accuracy of our simulations, the values of the critical exponents of our 3D system are comparable to those of 3D random percolation theory. For narrow slit-pores, the estimated exponents are found to be close to those obtained from the random percolation theory in two dimensions.

preprint2012arXiv

Feedback-induced oscillations in one-dimensional colloidal transport

We investigate a driven, one-dimensional system of colloidal particles in a periodically currogated narrow channel subject to a time-delayed feedback control. Our goal is to identify conditions under which the control induces oscillatory, time-periodic states. The investigations are based on the Fokker-Planck equation involving the density distribution of the system. First, by using the numerical continuation technique, we determine the linear stability of a stationary density. Second, the nonlinear regimes are analyzed by studying numerically the temporal evolution of the first moment of the density distribution. In this way we construct a bifurcation diagram revealing the nature of the instability. Apart from the case of a system with periodic boundary conditions, we also consider a microchannel of finite length. Finally, we study the influence of (repulsive) particle interactions based on Dynamical Density Functional Theory (DDFT).

preprint2012arXiv

Non-equilibrium condensation and coarsening of field-driven dipolar colloids

In colloidal suspensions, self-organization processes can be easily fueled by external fields. One particularly interesting class of phenomena occurs in monolayers of dipolar particles that are driven by rotating external fields. Here we report results from a computer simulation study of such systems focusing on the clustering behavior also observed in recent experiments. The key result of this paper is a novel interpretation of this pattern formation phenomenon: We show the clustering to be a by-product of a vapor-liquid first order phase transition. In fact, the observed dynamic coarsening process corresponds to the spindodal demixing that occurs during such a transition

preprint2012arXiv

Percolation and orientational ordering in systems of magnetic nanorods

Based on Monte Carlo (MC) computer simulations we study the structure formation of a system of magnetic nanorods. Our model particles consist of fused spheres with permanent magnetic dipole moments, as inspired by recent experiments. The resulting system behaves significantly different from a system of hard (non-magnetic) rods or magnetic rods with a single longitudinal dipole. In particular, we observe for the magnetic nanorods a significant decrease of the percolation threshold (as compared to non-magnetic rods) at low densities, and a stabilization of the high-density nematic phase. Moreover, the percolation threshold is tunable by an external magnetic field.

preprint2012arXiv

Rotational ratchets with dipolar interactions

We report results from a computer simulation study on the rotational ratchet effect in systems of magnetic particles interacting via dipolar interactions. The ratchet effect consists of directed rotations of the particles in an oscillating magnetic field, which lacks a net rotating component. Our investigations are based on Brownian dynamics simulations of such many-particle systems. We investigate the influence of both, the random and deterministic contributions to the equations of motion on the ratchet effect. As a main result, we show that dipolar interactions can have an enhancing as well as a dampening effect on the ratchet behavior depending on the dipolar coupling strength of the system under consideration. The enhancement is shown to be caused by an increase in the effective field on a particle generated by neighboring magnetic particles, while the dampening is due to restricted rotational motion in the effective field. Moreover, we find a non-trivial influence of the short-range, repulsive interaction between the particles.

preprint2011arXiv

Pattern formation of dipolar colloids in rotating fields: Layering and synchronization

We report Brownian dynamics (BD) simulation and theoretical results for a system of spherical colloidal particles with permanent dipole moments in a rotating magnetic field. Performing simulations at a fixed packing fraction and dipole coupling parameter, we construct a full non-equilibrium phase diagram as function of the driving frequency ($ω_0$) and field strength ($B_0$). This diagram contains both synchronized states, where the individual particles follow the field with (on average) constant phase difference, and asynchronous states. The synchronization is accompanied by layer formation, i.e. by spatial symmetry-breaking, similar to systems of induced dipoles in rotating fields. In the permanent-dipole case, however, too large $ω_0$ yield a breakdown of layering, supplemented by complex changes of the single-particle rotational dynamics from synchronous to asynchronous behavior. We show that the limit frequencies $ω_c$ can be well described as a bifurcation in the nonlinear equation of motion of a single particle rotating in a viscous medium. Finally, we present a simple density functional theory, which describes the emergence of layers in perfectly synchronized states as an equilibrium phase transition.

preprint2010arXiv

Crossover from normal to anomalous diffusion in field-aligned dipolar systems

Using molecular dynamics simulations we investigate the translational dynamics of particles with dipolar interactions in homogenous external fields. For a broad range of concentrations, we find that the anisotropic, yet normal diffusive behavior characterizing weakly coupled systems becomes anomalous both parallel and perpendicular to the field at sufficiently high dipolar coupling and field strength. After the ballistic regime, chain formation first yields cage-like motion in all directions, followed by transient, mixed diffusive-superdiffusive behavior resulting from cooperative motion of the chains. The enhanced dynamics disappears only at higher densities close to crystallization.

preprint2010arXiv

Feedback-controlled transport in an interacting colloidal system

Based on dynamical density functional theory (DDFT) we consider a non-equilibrium system of interacting colloidal particles driven by a constant tilting force through a periodic, symmetric "washboard" potential. We demonstrate that, despite of pronounced spatio-temporal correlations, the particle current can be reversed by adding suitable feedback control terms to the DDFT equation of motion. We explore two distinct control protocols with time delay, focussing on either the particle positions or the density profile. Our study shows that the DDFT is an appropriate framework to implement time-delayed feedback control strategies widely used in other fields of nonlinear physics

preprint2010arXiv

Non-linear rheology of active particle suspensions: Insights from an analytical approach

We consider active suspensions in the isotropic phase subjected to a shear flow. Using a set of extended hydrodynamic equations we derive a variety of {\em analytical} expressions for rheological quantities such as shear viscosity and normal stress differences. In agreement to full-blown numerical calculations and experiments we find a shear thickening or -thinning behaviour depending on whether the particles are contractile or extensile. Moreover, our analytical approach predicts that the normal stress differences can change their sign in contrast to passive suspensions.

preprint2010arXiv

Self-assembly of amphiphilic Janus particles at planar walls: A density functional study

We investigate the structure formation of amphiphilic molecules at planar walls using density functional theory. The molecules are modeled as (hard) spheres composed of a hydrophilic and hydrophobic part. The orientation of the resulting Janus-particles is described as a vector representing an internal degree of freedom. Our density functional approach involves Fundamental Measure Theory combined with a mean-field approximation for the anisotropic interaction. Considering neutral, hydrophilic and hydrophobic walls, we study the adsorption of the particles, focussing on the competition between the surface field and interact ion-induced ordering phenomena. Finally, we consider systems confined between two planar walls. It is shown that the anisotropic Janus interaction yields pronounced frustration effects at low temperatures.

preprint2010arXiv

Shear-stress controlled dynamics of nematic complex fluids

Based on a mesoscopic theory we investigate the non-equilibrium dynamics of a sheared nematic liquid, with the control parameter being the shear stress $σ_{\mathrm{xy}}$ (rather than the usual shear rate, $\dotγ$). To this end we supplement the equations of motion for the orientational order parameters by an equation for $\dotγ$, which then becomes time-dependent. Shearing the system from an isotropic state, the stress- controlled flow properties turn out to be essentially identical to those at fixed $\dotγ$. Pronounced differences when the equilibrium state is nematic. Here, shearing at controlled $\dotγ$ yields several non-equilibrium transitions between different dynamic states, including chaotic regimes. The corresponding stress-controlled system has only one transition from a regular periodic into a stationary (shear-aligned) state. The position of this transition in the $σ_{\mathrm{xy}}$-$\dotγ$ plane turns out to be tunable by the delay time entering our control scheme for $σ_{\mathrm{xy}}$. Moreover, a sudden change of the control method can {\it stabilize} the chaotic states appearing at fixed $\dotγ$.