### Inhalt des Dokuments

# Statistische Physik weicher Materie und biologischer Systeme

- [1]
- © AG Stark

Die Arbeitsgruppe bearbeitet ein breites Spektrum von Projekten aus den Bereichen der statistischen Physik weicher Materie und biologischer Systeme. Neben Frage- stellungen im thermischen Gleichgewicht werden vor allem Systeme untersucht, die durch externe Felder im Nichtgleichgewicht gehalten werden oder von sich heraus einen Nichtgleichgewichtszustand annehmen, wie zum Beispiel aktive Brownsche Teilchen. Dabei lassen wir uns von hochaktuellen Forschungsgebieten wie der Mikrofluidik und Biofluiddynamik leiten.

Derzeit ist ein Schwerpunkt unserer Forschung die aktive Bewegung in künstlichen und in biologischen Systemen. Dabei geht es um das Verständnis von Design- prinzipen zur Generierung aktiver Bewegung, um deren generische Eigenschaften und um kollektives Verhalten.

Unser Themenspektrum umfasst:

- Weiche Materie: insbes. Flüssigkristalle, kolloidale Dispersionen, (Bio-)Polymere
- Strukturbildung weicher Materie an Grenzflächen
- Weiche Materie und biologische Systeme im Nichtgleichgewicht, Mikrofluidik
- Biofluiddynamik und Fortbewegung von Mikroorganismen, aktive Bewegung
- Biomimetische Materialien und Biomechanik

# Aktuelles

Seminar

Die AG Stark bietet im Wintersemester 2021/2022 das Seminar

Themen und Methoden der Statistischen Mechanik weicher Materie und
biologischer Systeme [2] an.

Auf der ISIS-Seite [3]
finden Sie aktuelle Informationen und den Zoom-Link zur
Online-Veranstaltung**Bitte melden Sie sich über das
SAP-System an, um Modulpunkte erhalten zu können.**

Weitere Neuigkeiten:

**Q&A with Prof. Holger Stark, Editor-in-Chief of EPJ E**

Link to the Springer website [4]

**Prof. Stark bietet folgende Prüfungstermine an**- Dönnerstag, 16.12.2021
- Donnerstag, 20.01.2022
- Donnerstag, 17.02.2022

** Bitte
melden Sie sich direkt über das SAP-Sytem zur Prüfung
an!**

**Brennpunkt-Artikel im Physik Journal:**

Die Körperform macht's

Holger Stark

Physik Journal [5]

- [6]
- © Zantop, Stark

**Multi-Particle
Collision Dynamics with a Non-Ideal Equation of State: Part
I**

**Arne Zantop and Holger Stark**

The method of multi-particle collision dynamics (MPCD) and its different implementations are commonly used in the field of soft matter physics to simulate fluid flow at the micron scale. Typically, the coarse-grained fluid particles are described by the equation of state of an ideal gas, and the fluid is rather compressible. This is in contrast to conventional fluids, which are incompressible for velocities much below the speed of sound, and can cause inhomogeneities in density. We propose an algorithm for MPCD with a modified collision rule that results in a non-ideal equation of state and a significantly decreased compressibility. It allows simulations at less computational costs compared to conventional MPCD algorithms. We derive analytic expressions for the equation of state and the corresponding compressibility as well as shear viscosity. They show overall very good agreement with simulations, where we determine the pressure by simulating a quiet bulk fluid and the shear viscosity by simulating a linear shear flow and a Poiseuille flow.

J. Chem. Phys. 154, 024105 (2021) [7]

- [8]
- © Kneževic, Stark

**Effective Langevin
equations for a polar tracer in an active bath**

**by Milos Knezevic and Holger Stark**

We study the motion of a polar tracer, having a concave surface, immersed in a two-dimensional suspension of active particles. Using Brownian dynamics simulations, we measure the distributions and auto-correlation functions of force and torque exerted by active particles on the tracer. The tracer experiences a finite average force along its polar axis, while all the correlation functions show exponential decay in time. Using these insights we construct the full coarse-grained Langevin description for tracer position and orientation, where the active particles are subsumed into an effective self-propulsion force and exponentially correlated noise for both translations and rotations. The ensuing mesoscopic dynamics can be described in terms of five dimensionless parameters. We perform a thorough parameter study of the mean squared displacement, which illustrates how the different parameters influence the tracer dynamics, which crosses over from a ballistic to diffusive motion. We also demonstrate that the distribution of tracer displacements evolves from a non-Gaussian shape at early stages to a Gaussian behavior for sufficiently long times. Finally, for a given set of microscopic parameters, we establish a procedure to estimate the matching parameters of our effective model, and show that the resulting dynamics is in a very good quantitative agreement with the one obtained in Brownian dynamics simulations.

New J. Phys. **22**, 113025 (2020) [9]

- [10]
- © van Buel, Stark

**Active open-loop
control of elastic turbulence**

**by Reinier van Buel and Holger Stark**

We demonstrate through numerical solutions of the Oldroyd-B model in a two-dimensional Taylor–Couette geometry that the onset of elastic turbulence in a viscoelastic fluid can be controlled by imposed shear-rate modulations, one form of active open-loop control. Slow modulations display rich and complex behavior where elastic turbulence is still present, while it vanishes for fast modulations and a laminar response with the Taylor–Couette base flow is recovered. We find that the transition from the laminar to the turbulent state is supercritical and occurs at a critical Deborah number. In the state diagram of both control parameters, Weissenberg versus Deborah number, we identify the region of elastic turbulence. We also quantify the transition by the flow resistance, for which we derive an analytic expression in the laminar regime within the linear Oldroyd-B model. Finally, we provide an approximation for the transition line in the state diagram introducing an effective critical Weissenberg number in comparison to constant shear. Deviations from the numerical result indicate that the physics behind the observed laminar-to-turbulent transition is more complex under time-modulated shear flow.

Sci. Rep. **10**, 15704 (2020)

- [11]
- © EDP Sciences, SIF and Springer-Verlag

**Particle
pairs and trains in inertial microfluidics**

**by Christian Schaaf and Holger Stark**

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

Eur. Phys. J. E **43**, 50 (2020). [12]

selected for front cover picture [13]

- [14]
- © Rühle, Schaaf, Stark

**Optimal Control of
Colloidal Trajectories in Inertial Microfluidics Using the Saffman
Effect**

**by Felix Rühle, Christian Schaaf, and Holger
Stark**

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

Micromachines 11, 592 (2020) [15]

- [16]
- © Zantop, Stark

**Squirmer rods as
elongated microswimmers: flow fields and confinement**

**by Arne W. Zantop and Holger Stark**

Microswimmers or active elements, such as bacteria and active filaments, have an elongated shape, which determines their individual and collective dynamics. There is still a need to identify what role long-range hydrodynamic interactions play in their fascinating dynamic structure formation. We construct rods of different aspect ratios using several spherical squirmer model swimmers. With the help of the mesoscale simulation method of multi-particle collision dynamics we analyze the flow fields of these squirmer rods both in a bulk fluid and in Hele-Shaw geometries of different slab widths. Based on the hydrodynamic multipole expansion either for bulk or confinement between two parallel plates, we categorize the different multipole contributions of neutral as well as pusher-type squirmer rods. We demonstrate how confinement alters the radial decay of the flow fields for a given force or source multipole moment compared to the bulk fluid.

Soft Matter **16**, 6400 (2020) [17]

- [18]
- © EDP Sciences, SIF and Springer-Verlag

**Emergent
collective dynamics of bottom-heavy squirmers under
gravity**

**by Felix Rühle and Holger Stark**

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

Eur. Phys. J. E **43**, 26 (2020)

see also: *'Bottom-heavy' squirmers adopt characteristic group
behaviours*, article as EPJE Highlight [19] and at EurekAlert
[20]

- [21]
- © Europhysics Letters Association

**Capillary
condensation in an active bath**

**by Milos Knezevic and Holger Stark**

e study capillary condensation in a bath of active Brownian particles (ABPs) and the forces acting on the capillary close to the motility-induced phase separation (MIPS). The capillary is modeled as two parallel rods, which are fixed in space. We consider a bath of ABPs having a self-propulsion speed much larger than the critical speed necessary for MIPS to occur.We gradually increase the packing fraction of ABPs, starting from a dilute phase of ABPs and going towards the binodal of MIPS. In stark contrast to conventional capillary condensation, we do not observe any hysteresis in the capillary packing fraction and attribute this to strong temporal fluctuations in the capillary packing fraction . Depending on the packing fraction of ABPs and capillary width, we find that the effective force between the capillary rods can be either attractive or repulsive. In fact,with increasing width it shows damped oscillations as long as capillary condensation occurs. We analyze them in detail by studying the distribution of particle distances from the inner and outer wall of the capillary, respectively. In addition, we examine the capillary in the active bath close to the critical point. We do not observe signs of the presence of long-range Casimir interactions.

Europhys. Lett. **128**, 40008 (2020)

- [22]
- © Copyright??

**A flowing pair of particles in inertial
microfluidics**

**by Christian Schaaf, Felix Rühle and Holger
Stark**

A flowing pair of particles in inertial microfluidics gives important insights into understanding and controlling the collective dynamics of particles like cells or droplets in microfluidic devices. They are applied in medical cell analysis and engineering. We study the dynamics of a pair of solid particles flowing through a rectangular microchannel using lattice Boltzmann simulations. We determine the inertial lift force profiles as a function of the two particle positions, their axial distance, and the Reynolds number. Generally, the profiles strongly differ between particles leading and lagging in flow and the lift forces are enhanced due to the presence of a second particle. At small axial distances, they are determined by viscous forces, while inertial forces dominate at large separations. We identify cross-streamline pairs as stable fixed points in the lift force profiles and argue that same-streamline configurations are only one-sided stable. Depending on the initial conditions, the two-particle lift forces in combination with the Poiseuille flow give rise to three types of unbound particle trajectories, called moving-apart, passing, and swapping, and one type of bound trajectory, where the particles perform damped oscillations towards the cross-stream line configuration. The damping rate scales with Reynolds number squared, since inertial forces are responsible for driving the particles to their steady-state positions.

Soft Matter **15**, 1988 (2019)

- [23]
- © Copyright??

**Traveling concentration pulses of bacteria
in a generalized Keller–Segel model**

**by Maximilian Seyrich, Andrzej Palugniok, and Holger
Stark**

We formulate a Markovian response theory for the tumble rate of a bacterium moving in a chemical field and use it in the Smoluchowski equation. Based on a multipole expansion for the one-particle distribution function and a reaction-diffusion equation for the chemoattractant field, we derive a polarization extended model, which also includes the recently discovered angle bias. In the adiabatic limit we recover a generalized Keller–Segel equation with diffusion and chemotactic coefficients that depend on the microscopic swimming parameters. Requiring the tumble rate to be positive, our model introduces an upper bound for the chemotactic drift velocity, which is no longer singular as in the original Keller–Segel model. Solving the Keller–Segel equations numerically, we identify traveling bacterial concentration pulses, for which we do not need a second, signaling chemical field nor a singular chemotactic drift velocity as demanded in earlier publications. We present an extensive study of the traveling pulses and demonstrate how their speeds, widths, and heights depend on the microscopic parameters. Most importantly, we discover a maximum number of bacteria that the pulse can sustain—the maximum carrying capacity. Finally, by tuning our parameters, we are able to match the experimental realization of the traveling bacterial pulse.

New J. Phys. **21, **103001 (2019)

- [24]
- © Copyright??

**Chemotaxis in a binary mixture of active
and passive particles**

**by Julian Stürmer , Maximilian Seyrich , and Holger
Stark**

Mixtures of active and passive colloids show an intriguing dynamics
of self-assembling, which is driven by the active component.
Self-phoretic active colloids generate sinks in a chemical
concentration field that cause passive colloids to drift toward active
colloids by diffusiophoresis. The strength of this effective
attraction is governed by the diffusiophoretic parameter, which
determines the drift velocity. Simulating the Langevin dynamics of the
colloids, we determine the state diagram for increasing
diffusiophoretic strength and fixed active velocity. Three main states
are distinguished. For weak attraction, passive particles are first
scattered in the simulation box and then form a colloidal cloud around
its center. Increasing the diffusiophoretic parameter further, passive
particles oscillate between the cloud and a compact cluster, which
embeds active colloids. Ultimately, in the third state, all particles
collapse into a single stable cluster. In the collapse regime, the
clustering dynamics of the largest cluster follows a logistic function
and the mean cluster velocity vs cluster size decays with a power law.
Throughout this article, we discuss our simulation results with regard
to the experiments of Singh et al., Adv. Mater. **29**
(32), 1701328 (2017).

J. Chem. Phys. **150**, 214901 (2019)

- [25]
- © Copyright??

**Optimal steering of a smart active
particle**

**by Elias Schneider and Holger Stark**

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

EPL **127**, 64003 (2019) [26]

- [27]
- © Copyright??

**Collective Dynamics in a Monolayer of
Squirmers Confined to a Boundary by Gravity**

**by Jan-Timm Kuhr, Felix Rühle, and Holger
Stark**

We present a hydrodynamic study of a monolayer of squirmer model microswimmers confined to aboundary by strong gravity using the simulation method of multi-particle collision dynamics. Thesquirmers interact with each other via their self-generated hydrodynamic flow fields and thereby form avariety of fascinating dynamic states when density and squirmer type are varied. Weak pushers, neutralsquirmers, and pullers have an upright orientation. With their flow fields they push neighbors away andthereby form a hydrodynamic Wigner fluid at lower densities. Furthermore, states of fluctuating chainsand trimers, of kissing, and at large densities a global cluster exist. Finally, pushers at all densities can tiltagainst the wall normal and their in-plane velocities align to show swarming. It turns into chaoticswarming for strong pushers at high densities. We characterize all these states quantitatively.

Soft Matter **15**, 5685 (2019)

- [28]
- © Copyright??

**Artificial Chemotaxis of Self-Phoretic
Active Colloids: Collective Behavior**

**by Holger Stark**

Microorganisms use chemotaxis, regulated by internal complex
chemical pathways, to swim along chemical gradients to find better
living conditions. Artificial microswimmers can mimic such a strategy
by a pure physical process called diffusiophoresis, where they drift
and orient along the gradient in a chemical density field. Similarly,
for other forms of taxis in nature such as photo- or thermotaxis the
phoretic counterpart exists.

In this Account,
we concentrate on the chemotaxis of self-phoretic active colloids.
They are driven by self-electro- and diffusiophoresis at the particle
surface and thereby acquire a swimming speed. During this process,
they also produce nonuniform chemical fields in their surroundings
through which they interact with other colloids by translational and
rotational diffusiophoresis. In combination with active motion, this
gives rise to effective phoretic attraction and repulsion and thereby
to diverse emergent collective behavior. A particular appealing
example is dynamic clustering in dilute suspensions first reported by
a group from Lyon. A subtle balance of attraction and repulsion causes
very dynamic clusters, which form and resolve again. This is in stark
contrast to the relatively static clusters of motility-induced phase
separation at larger densities.

To treat
chemotaxis in active colloids confined to a plane, we formulate two
Langevin equations for position and orientation, which include
translational and rotational diffusiophoretic drift velocities. The
colloids are chemical sinks and develop their long-range chemical
profiles instantaneously. For dense packings, we include screening of
the chemical fields. We present a state diagram in the two
diffusiophoretic parameters governing translational, as well as
rotational, drift and, thereby, explore the full range of phoretic
attraction and repulsion. The identified states range from a gaslike
phase over dynamic clustering states 1 and 2, which we distinguish
through their cluster size distributions, to different types of
collapsed states. The latter include a full chemotactic collapse for
translational phoretic attraction. Turning it into an effective
repulsion, with increasing strength first the collapsed cluster starts
to fluctuate at the rim, then oscillates, and ultimately becomes a
static collapsed cloud. We also present a state diagram without
screening. Finally, we summarize how the famous Keller−Segel model
derives from our Langevin equations through a multipole expansion of
the full one-particle distribution function in position and
orientation. The Keller−Segel model gives a continuum equation for
treating chemotaxis of microorganisms on the level of their spatial
density.

Our theory is extensible to mixtures
of active and passive particles and allows to include a dipolar
correction to the chemical field resulting from the dipolar symmetry
of Janus colloids.

Acc. Chem. Res. **51**, 2681 (2018)

- [29]
- © Copyright??

**Elastic turbulence in two-dimensional
Taylor-Couette flows **

**by R. van Buel, C. Schaaf, and H. Stark**

We report the onset of elastic turbulence in a two-dimensional
Taylor-Couette geometry using numerical solutions of the Oldroyd-B
model generated with the program OpenFOAM®. Beyond a critical
Weissenberg number, an elastic instability causes a supercritical
transition from the laminar Taylor-Couette flow to a turbulent flow.
The order parameter, the time average of secondary-flow strength,
follows the scaling law Φ∝ Wi -Wi_{c}^{γ} with Wi_{c}=10 and γ=
0.45 . Additionally, the flow resistance increases beyond Wi_{c} . The temporal power spectra of the velocity
fluctuations show a power-law decay with a characteristic exponent in
the range 2< α < 4 , which strongly depends on the
radial position. The characteristic exponent β for the spatial power
spectra obeys the necessary condition β > 3 , associated with
elastic turbulence, for all Wi > Wi_{c} .

Europhys. Lett. **124,** 14001 (2018)

- [30]
- © Copyright??

**Feedback control of photoresponsive fluid
interfaces**

**by Josua Grawitter and Holger Stark**

Photoresponsive surfactants provide a unique microfluidic driving mechanism. Since they switch between two molecular shapes under illumination and thereby affect surface tension of fluid interfaces, Marangoni flow along the interface occurs. To describe the dynamics of the surfactant mixture at a planar interface, we formulate diffusion–advection–reaction equations for both surfactant densities. They also include adsorption from and desorption into the neighboring fluids and photoisomerization by light. We then study how the interface responds when illuminated by spots of light. Switching on a single light spot, the density of the switched surfactant spreads in time and assumes an exponentially decaying profile in steady state. Simultaneously, the induced radial Marangoni flow reverses its flow direction from inward to outward. We use this feature to set up specific feedback rules, which couple the advection velocities sensed at the light spots to their intensities. As a result two neighboring spots switch on and off alternately. Extending the feedback rule to light spots arranged on the vertices of regular polygons, we observe periodic switching patterns for even-sided polygons, where two sets of next-nearest neighbors alternate with each other. A triangle and pentagon also show regular oscillations, while heptagon and nonagon exhibit irregular oscillations due to frustration. While our findings are specific to the chosen set of parameters, they show how complex patterns at photoresponsive fluid interfaces emerge from simple feedback coupling.

Soft Matter **14**,
1856 (2018)

AG_Stark/Bilder_Aktuelle_Forschung/AB_emblem.png

e/seminar_statistische_physik_weicher_materie_und_biolo

gischer_systeme/parameter/de/font0/maxhilfe/

er-physics/blog2

14

AG_Stark/Bilder_Aktuelle_Forschung/2021/zantop_mpcd_202

1.png

/AG_Stark/Publikationen_PDF/2021/zantop_mpcd_2021.pdf

AG_Stark/Bilder_Aktuelle_Forschung/2020/fig1b.png

/AG_Stark/Publikationen_PDF/2020/Knezevic_2020_New_J_Ph

ys_22_113025.pdf

/AG_Stark/Bilder_Aktuelle_Forschung/2020/vanBuelSciRep2

020.png

/AG_Stark/Publikationen_PDF/2020/schaaf_2020_cover_epje

.png

d/AG_Stark/Publikationen_PDF/2020/Schaaf.pdf

d/AG_Stark/Publikationen_PDF/2020/schaaf_epje_front_cov

er2020.pdf

/AG_Stark/Publikationen_PDF/2020/micromachines-11-00592

-ag-550.jpg

d/AG_Stark/Publikationen_PDF/2020/RuehleStark2020_micro

machines.pdf

/AG_Stark/Publikationen_PDF/2020/zantop_abstract.png

d/AG_Stark/Publikationen_PDF/2020/Zantop.pdf

/AG_Stark/Publikationen_PDF/2020/10189_0043_005_2020-xe

xternal.pdf

bottom-heavy-squirmers-adopt-characteristic-group-behav

iours

sa052820.php

/AG_Stark/Bilder_Aktuelle_Forschung/2020/Knezeciv.png

/AG_Stark/Bilder_Aktuelle_Forschung/2019/SchaafSoftMatt

er.png

/AG_Stark/Bilder_Aktuelle_Forschung/2019/GraphicalAbstr

act.png

/AG_Stark/Bilder_Aktuelle_Forschung/2019/Stuermer.png

/AG_Stark/Bilder_Aktuelle_Forschung/2019/Schneider.png

d/AG_Stark/Publikationen_PDF/2019/Schneider.pdf

/AG_Stark/Bilder_Aktuelle_Forschung/2019/kuhr.png

/AG_Stark/Bilder_Aktuelle_Forschung/2018/StarkBild.pdf

/AG_Stark/Bilder_Aktuelle_Forschung/2018/VanBuel.png

/AG_Stark/Bilder_Aktuelle_Forschung/2018/Grawitter.png

ge_ag_stark/current_research_of_previous_years/paramete

r/de/font0/maxhilfe/