# Stokesian dynamics simulations of a magnetotactic bacterium

Sarah Mohammadinejad, Damien Faivre, Stefan Klumpp

► To cite this version:

Sarah Mohammadinejad, Damien Faivre, Stefan Klumpp. Stokesian dynamics simulations of a magnetotactic bacterium. European Physical Journal E: Soft matter and biological physics, 2021, 44 (3), 10.1140/epje/s10189-021-00038-5. [DOI: https://doi.org/10.1140/epje/s10189-021-00038-5](https://doi.org/10.1140/epje/s10189-021-00038-5)

HAL Id: hal-03296278

https://amu.hal.science/hal-03296278

Submitted on 23 Jul 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Stokesian dynamics simulations of a magnetotactic bacterium

Sarah Mohammadinejad<sup>1,2,3,a</sup>, Damien Faivre<sup>4,5</sup>, and Stefan Klumpp<sup>1,2</sup>

Received 29 September 2020 / Accepted 15 February 2021 / Published online 23 March 2021

© The Author(s) 2021

**Abstract** The swimming of bacteria provides insight into propulsion and steering under the conditions of low-Reynolds number hydrodynamics. Here we address the magnetically steered swimming of magnetotactic bacteria. We use Stokesian dynamics simulations to study the swimming of single-flagellated magnetotactic bacteria (MTB) in an external magnetic field. Our model MTB consists of a spherical cell body equipped with a magnetic dipole moment and a helical flagellum rotated by a rotary motor. The elasticity of the flagellum as well as magnetic and hydrodynamic interactions is taken into account in this model. We characterized how the swimming velocity is dependent on parameters of the model. We then studied the U-turn motion after a field reversal and found two regimes for weak and strong fields and, correspondingly, two characteristic time scales. In the two regimes, the U-turn time is dominated by the turning of the cell body and its magnetic moment or the turning of the flagellum, respectively. In the regime for weak fields, where turning is dominated by the magnetic relaxation, the U-turn time is approximately in agreement with a theoretical model based on torque balance. In the strong-field regime, strong deformations of the flagellum are observed. We further simulated the swimming of a bacterium with a magnetic moment that is inclined relative to the flagellar axis. This scenario leads to intriguing double helical trajectories that we characterize as functions of the magnetic moment inclination and the magnetic field. For small inclination angles (≤ 20°) and typical field strengths, the inclination of the magnetic moment has only a minor effect on the swimming of MTB in an external magnetic field. Large inclination angles result in a strong reduction in the velocity in direction of the magnetic field, consistent with recent observations that bacteria with large inclination angles use a different propulsion mechanism.

## 1 Introduction

The swimming of motile microorganisms is of great scientific interest, for understanding motility in natural systems, as instances of active matter with all its intriguing physical properties and as inspiration for the development of synthetic microswimmers and micro-robots [1–3]. Microorganisms use diverse propulsion mechanisms and behavioral strategies in their swimming. While some of them such as the run-and-tumble motion of *Escherichia coli* [4,5] or the swimming sperm [6] (and in both cases, their use in chemosensing) are well studied and understood in great detail, other mechanisms of swimming and the corresponding strategies for sensing and steering have only been described recently [7–9]. Among microorganisms, magnetotactic bacteria (MTB) are of particular interest because they

can be steered by a magnetic field [10,11], which provides a promising biocompatible option for the remote control of swimming. MTB are a diverse group of microorganisms with the ability to orient and migrate along (geo-)magnetic field lines [10,12,13]. This ability is based on the presence of an intracellular structure called magnetosome chain, which acts as a magnetic compass for their navigation. The chain is composed of iron oxide or iron sulfide nanoparticles enclosed in membranous vesicles, called magnetosomes, which imparts a magnetic moment to MTB [14]. Following the Earth’s magnetic field lines is hypothesized to help them find their preferred habitat near the oxic-anoxic transition zone in stratified aquatic environment and sediments [10]. This type of aerotaxis, which is assisted by magnetic field is called magnetoaerotaxis. From a physical point of view, magnetotaxis can either provide a direction or only an axis, but in either case reduces a three-dimensional search process to a one-dimensional search along the direction defined by the field [15,16].

<footnote n="1">¹ Institute for the Dynamics of Complex Systems, University of Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany</footnote>

<footnote n="2">² Department Theory and Bio-Systems, Max Planck Institute of Colloids and Interfaces, 14424 Potsdam, Germany</footnote>

<footnote n="3">³ Department of Biological Sciences, Institute for Advanced Studies in Basic Sciences (IASBS), Zanjan 45137-66731, Iran</footnote>

<footnote n="4">⁴ Department of Biomaterials, Max Planck Institute of Colloids and Interfaces, 14424 Potsdam, Germany</footnote>

<footnote n="5">⁵ Aix-Marseille Université, CEA, CNRS, BIAM, 13108 Saint-Paul-lez-Durance, France</footnote>

<footnote>a e-mail: sarah.mohammadinejad@phys.uni-goettingen.de
 corresp odin g author)</footnote>
In recent years, MTB have attracted the attention of researchers from many different disciplines [11], including from an application point of view. MTB have been functionalized for their use as a magnetically targeted drug delivery system [17]. They also provide inspiration for synthetic and biohybrid magnetic transport: for example, magnetic guidance has been demonstrated after incorporation or attachment of magnetic particles to naturally non-magnetic organisms [18, 19] and synthetic magnetically steered swimmers were constructed with the purpose of finding the best structure for rapid swimming while being steering by magnetic field [20–22].

Despite the promising proof-of-concept applications of MTB, many aspects of the dynamics of MTB are only beginning to be known such as details of their swimming mechanisms, their strategies for responding to environmental stimuli and changing the direction of motion. Since MTB are rather diverse with respect to cell shapes, magnetic dipole orientation and number and organization of their flagella, a variety of mechanisms and strategies are expected to be observed [15].

Few experimental studies have addressed the detailed mechanisms of propulsion of MTB by their flagella. Notable exceptions are studies on the coordination of two flagella at opposite poles in *Magnetospirillum magneticum* (AMB-1) [23], of two velocity modes in *Magnetospirillum gryphiswaldense* (MSR-1) [24] and of the cooperation of two flagellar bundles almost perpendicular to the magnetic moment in *Magnetococcus marinus* (MC-1) [9]. These studies are complemented by theoretical approaches. Here simple models describe the bacteria either as self-propelled dipolar particles [25–27] or as a cell body with a rigid rotating helix [28–30] and did not include the hydrodynamics and/or mechanics of the flagellum. In this paper, we use Stokesian dynamics, which provides an efficient framework to carry out numerical simulations for particles suspended in a fluid that interact through hydrodynamic and other forces [31]. Stokesian dynamics has been used extensively for the study of the swimming of non-magnetotactic bacteria, in particular for the swimming of *Escherichia coli* including different dynamical states of its flagellum and the hydrodynamic interaction between its multiple flagella [32–34]. The same model was used to investigate the origin of screw formation in the flagellum in single-flagellated *Shewanella putrefaciens* attempting to escape from the traps [8]. Likewise, in our recent work on the *Magnetococcus marinus*, we used the same approach as a complement to experimental three-dimensional tracking of swimming bacteria to distinguish different swimming modes and to propose that the two flagella of this bacterium rotate in opposite direction [9].

Here, we report on simulations of the swimming behavior of a model MTB with spherical cell body and a single flagellum. We characterized the influence of various parameters such as the flagellum length, the orientation of the magnetosome chain relative to the flagellar axis and the strength of the magnetic field. Specifically, we present two scenarios that are relevant

from an experimental point of view: we investigated the U-turns that are observed when the magnetic field is reversed and asked whether an inclination of the magnetic moment relative to the flagellar propulsion axis is detrimental to the swimming motion.

## 2 Method

In this paper, we used Stokesian dynamics simulations to study the swimming behavior of a single-flagellated magnetotactic bacterium. We present a semi-realistic model for MTB taking into account its magnetic orientation and its propulsion by the flagellar rotation. In the model, the cell body is spherical and has a single flagellum discretized as a chain of beads (shown in Fig. 1a, with yellow and blue colors, respectively). The flagellum, in turn, consists in a helical filament, a rotary motor embedded in the cell membrane, which rotates the filament, and a short hinge for transmitting the motor torque to the filament [33]. The description of the flagellum includes its elasticity, as well as hydrodynamic and excluded volume interactions. The hydrodynamic interactions are taken into account by calculating the mutual mobility coefficients between particles representing segments of the flagellum and then solving the equations of motion for all particles using these mobility coefficients [35].

### 2.1 Model of the flagellum

The model we used for the flagellum of our MTB is the model that has been used for the *E. coli* flagellum by Vogel et al. [20, 32]. In this model, the helical filament of the flagellum is described as a discretized space curve composed of $N$ particles of diameter $a = 0.02$ µm and equilibrium separation distance $L_e = 0.2$ µm (Fig. 1a). In order to explore the bacterial behavior as a function of the flagellum length, $N$ was varied between 20 and 150. As shown in Fig. 1b, the position of each flagellum particle is represented by $\mathbf{r}_i$ ($i = 1, ..., N$). A set of three orthogonal unit vectors $\mathbf{e}_i^\alpha = \{\mathbf{e}_i^1, \mathbf{e}_i^2, \mathbf{e}_i^3\}$ (called the $i$-th triad with $\alpha = 1, 2, 3$ representing the three orthogonal unit vectors) is assigned to each particle $i$ to describe the orientation and elastic deformation of the bond between two successive particles $i$ and $i-1$ (Fig. 1b). The first triad $\mathbf{e}_1^\alpha$ is associated with the bond between the cell body and the first flagellar particle.

A regular helical filament like the one we have for the flagellum can be characterized completely by its constant curvature ($\kappa_e$) and torsion ($\tau_e$) which are related to the helix radius ($R$) and pitch ($P$) via $R = \frac{\kappa_e}{\kappa_e^2 + \tau_e^2}$ and $P = \frac{2\pi\tau_e}{\kappa_e^2 + \tau_e^2}$ [36].

Since the flagellum is an elastic filament, any deviation from its equilibrium configuration costs energy. The elastic energy for the bending and twisting of the flagellum ($H^K$) can be calculated by regarding the conformational deviation from the helical equilibrium configuration using a discretized version of Kirchhoff's the-
<figure><img src="image_2.png" /><figcaption>**Fig. 1** a Model of a magnetotactic bacterium with one flagellum. The flagellum and the cell body are shown in blue and yellow, respectively. The magnetic moment and the motor rotation axis are shown with purple and black arrows on the cell body. The flagellum is modeled as a discretized space curve composed of N particles. A part of the flagellum is magnified in **b** to indicate the triad vectors defined between two successive particles of the flagellum. **c** Representation of the bending and twist angles ($\theta_i$ and $\phi_i$) and the corresponding deformation of the triads along the flagellum length. The primed triad represents the intermediate state after the twist is applied, but before bending</figcaption></figure>

ory for elastic rods [37],

$$H^K = L_e \sum_{i=1}^{N-1} \{A[(\Omega_i^1 - \Omega_e^1)^2 + (\Omega_i^2 - \Omega_e^2)^2] + C(\Omega_i^3 - \Omega_e^3)^2\}, \quad (1)$$

in which $\Omega_i^\alpha = \{\Omega_i^1, \Omega_i^2, \Omega_i^3\}$ is the instantaneous local angular strain vector. Its components are related to the local curvature and torsion angles and the local triads via $\Omega_i^1 = -\frac{\theta_i}{\sin\theta_i}e_i^2 \cdot e_{i+1}^3$, $\Omega_i^2 = \frac{\theta_i}{\sin\theta_i}e_i^1 \cdot e_{i+1}^3$, $\Omega_i^3 = \phi_i$, where $\theta_i$ and $\phi_i$ (as shown in Fig. 1c) are the bending and twist angles between successive triads, respectively. A and C are the bending and torsional rigidities, for which the values of $A = C = 3.5$ pN µm² are used.

$\Omega_e^\alpha = \{\Omega_e^1, \Omega_e^2, \Omega_e^3\} = \{\kappa_e \cos(L_e \tau_e), \kappa_e \sin(L_e \tau_e), \tau_e\}$ is the equilibrium angular vector for the flagellum. Throughout this paper, we use $\kappa_e = 1.3 \mu m^{-1}$, $\tau_e =$

$-2.1 \mu m^{-1}$, except in the section where we study the effect of the flagellum helical geometry on the swimming velocity (Fig. 3b) where we run simulation for different $\Omega_e^\alpha$.

The corresponding elastic forces on i-th particle ($F_i^K$) and torques about i-th triad ($T_i^K$) can be obtained by calculating the numerical derivatives of $H^K$ with respect to the position of the particle ($r_i$) and the twist angle of the triads ($\phi_i$) as follows,

$$F_i^K = -\frac{\partial H^K}{\partial r_i}, \quad (2)$$

$$T_i^K = -\frac{\partial H^K}{\partial \phi_i}. \quad (3)$$

The stretching elasticity of the flagellum is taken into account using a Hookian free energy ($H^S$) of the form

$$H^S = \frac{K_s}{2L_e} \sum_{i=2}^{N} (L_i - L_e)^2, \quad L_i = |\mathbf{r}_i - \mathbf{r}_{i-1}| \quad (4)$$

with stretching rigidity of $K_s = 1000$ pN [33]. The stretching force ($F_i^S$) can be obtained by analytical differentiation of $H^S$ with respect to $r_i$. The resulting $F_i^S$ is

$$F_i^S = -\frac{K_s}{L_e} \{(L_i - L_e)e_i^3 - (L_{i+1} - L_e)e_{i+1}^3\}. \quad (5)$$

The repulsive part of a Lennard-Jones force (a WCA force) is used to model excluded volume interactions between all particles including the cell body ($i = 0$),

$$\mathbf{F}_i^{LJ} = \begin{cases} \frac{F_0}{\sigma} \left[ 2 \left( \frac{\sigma}{r_{ij}} \right)^{14} - \left( \frac{\sigma}{r_{ij}} \right)^8 \right] \mathbf{r}_{ij} & \text{if } r_{ij} < \sqrt[6]{2}\sigma, \\ 0 & \text{if } r_{ij} \ge \sqrt[6]{2}\sigma, \end{cases} \quad (6)$$

with $F_0 = 1$ pN. We set $\sigma = 2a$ for the interaction between two (non-bonded) beads of the flagellum, so that the distance between two beads on which the steric interactions become relevant is $2a$. For the interaction between flagellum beads and the cell body, we use $\sigma = a + R_b$. $\mathbf{r}_{ij}$ is the distance vector between i-th and j-th particles (including the cell body ($i = 0$)) and $r_{ij}$ its absolute value. The steric force on the i-th particle ($F_i^{LJ}$) is calculated by summation of the steric forces from all other particles ($j$) on the i-th particle. We note that despite the relatively large gaps between subsequent flagellar beads ($L_e = 10a$), we do not observe beads passing through the flagellum due to its rather high bending and twist rigidities.

The filament is attached to the surface of the bacterial cell body via its first particle ($i = 1$). We defined the first triad ($e_i^\alpha$) in such a way that its third component ($e_i^3$) is along the vector drawn from body center
to the first flagellar particle attached to body ($i$ = 1-th particle). The effect of the rotary motor was modeled by the motor torque $T_m \mathbf{e}_1^3$, for which the value $T_m = 3.4$ pN µm was assumed, acting on the first ($i$ = 1) bead of the flagellum. This torque is transmitted to the main part of the filament via the hinge, represented by the connection to the second bead. For the latter we used the rigidities $A' = A$ and $C' = 3C$. The high torsional rigidity provides the coupling of the filament to the motor rotation.

## 2.2 Dynamics of the flagellum

At low Reynolds numbers, because of the linearity of Stokes' equation, the translational and angular velocity vectors ($\mathbf{v}$ and $\omega$, respectively) of the particles have a linear relation with the forces ($\mathbf{F}$) and torques ($\mathbf{T}$) acting on those particles as follows

$$ \begin{bmatrix} \mathbf{v} \\ \omega \end{bmatrix} = \begin{bmatrix} \boldsymbol{\mu}^{tt} & \boldsymbol{\mu}^{tr} \\ \boldsymbol{\mu}^{rt} & \boldsymbol{\mu}^{rr} \end{bmatrix} \begin{bmatrix} \mathbf{F} \\ \mathbf{T} \end{bmatrix}. \quad (7) $$

Here $\mathbf{v}$ and $\omega$ are 3N-component linear and angular velocity vectors. The matrix is called the mobility matrix and $\boldsymbol{\mu}^{tt}$, $\boldsymbol{\mu}^{rr}$, $\boldsymbol{\mu}^{rt}$ and $\boldsymbol{\mu}^{tr}$ are each $3N \times 3N$ mobility sub-matrices representing the translational–translational, rotational-rotational and rotational–translational coupling between velocities and forces of all particles. In the following we refer to the mobility terms between force and velocity of the same particle as self-mobilities and those between forces and velocities of different particles as cross-mobilities.

In this study, the rotational–translational mobility sub-matrices ($\boldsymbol{\mu}^{tr}$ and $\boldsymbol{\mu}^{rt}$) and rotational-rotational cross-mobility terms ($\boldsymbol{\mu}_{i \neq j}^{rr}$) were neglected. This is mostly done for technical reasons, since these couplings may result in the rotation of the triads and cause inconsistencies with the definition of the bond direction $\mathbf{e}_i^3 = r_i - r_{i-1}$. Artifacts were reported in earlier work, specially for elastic structure [38]. From the physical point of view, these contributions can be neglected because the $\boldsymbol{\mu}_{ij}^{rr}$ and $\boldsymbol{\mu}_{i \neq j}^{tr}$ terms in the Rotne–Prager expansion are of higher order in $\frac{a}{r_{ij}}$ compared to $\boldsymbol{\mu}_{ii}^{rr}$. The remaining mobility terms for flagellar particles are the translational–translational cross-mobility ($\boldsymbol{\mu}_{i \neq j}^{tt}$), the translational–translational self-mobility ($\boldsymbol{\mu}_{ii}^{tt}$) and the rotational-rotational self-mobility ($\boldsymbol{\mu}_{ii}^{rr}$). The translational–translational cross-mobilities ($\boldsymbol{\mu}_{i \neq j}^{tt}$) between flagellar particles were calculated using the Rotne–Prager approximation [35,39],

$$ \boldsymbol{\mu}_{i \neq j}^{tt} = \mu^t \left[ \frac{3}{4} \frac{a}{r_{ij}} (\mathbf{1} + \hat{\mathbf{r}}_{ij} \otimes \hat{\mathbf{r}}_{ij}) + \frac{1}{2} (\frac{a}{r_{ij}})^3 (\mathbf{1} - 3\hat{\mathbf{r}}_{ij} \otimes \hat{\mathbf{r}}_{ij}) \right]. \quad (8) $$

Here $r_{ij}$ represents again the distance between the i-th and j-th particles. For the translational–translational self-mobility ($\boldsymbol{\mu}_{ii}^{tt}$, a 3 × 3-submatrix of $\boldsymbol{\mu}^{tt}$) of the flagellar particles we used the mobility relation for a straight

rod that is

$$ \boldsymbol{\mu}_{ii}^{tt} = \mathbf{e}_i^3 \otimes \mathbf{e}_i^3 / \gamma_{\parallel} + (\mathbf{1} - \mathbf{e}_i^3 \otimes \mathbf{e}_i^3) / \gamma_{\perp}, \quad (9) $$

in which $\gamma_{\parallel} = 3.2 \times 10^{-4}$ pNs/µm and $\gamma_{\perp} = 5.6 \times 10^{-4}$ pNs/µm were used as the anisotropic friction coefficients of the flagellum following the values used for the E. coli flagellum [33]. The rotational-rotational self-mobility of flagellar particles ($\boldsymbol{\mu}_{ii}^{rr}$) is the inverse of the rotational friction $\boldsymbol{\mu}_{ii}^{rr} = 1/\gamma_r$ with $\gamma_r = 2.52 \times 10^{-7}$ pNs µm.

Taking into account the forces and torques discussed above and including them in Eq. (7), the dynamics of the flagellar particles' positions ($\mathbf{r}_i$) and triads ($\mathbf{e}_i^\alpha$ with $\alpha = 1,2,3$ can be obtained by solving the following equations of motion for their position and orientation,

$$ \frac{\partial \mathbf{r}_i}{\partial t} = \mu_{ij}^{tt} (\mathbf{F}_j^K + \mathbf{F}_j^S + \mathbf{F}_j^{LJ}), \quad i = 2, \dots, N \quad (10) $$

$$ \frac{\partial \mathbf{e}_i^\alpha}{\partial t} = \mu_{ii}^{rr} (T_i^K \mathbf{e}_i^3 \times \mathbf{e}_i^\alpha), \quad i = 2, \dots, N \quad (11) $$

$$ \frac{\partial \mathbf{e}_{\mathbf{f}_1}^{\alpha}}{\partial t} = \mu_{11}^{rr} (T_1^K + T_m)(\mathbf{e}_1^3 \times \mathbf{e}_1^{\alpha}), \quad (12) $$

in which $\frac{\partial \mathbf{e}_{\mathbf{r}_1}^{\alpha}}{\partial t}$ is the contribution of the flagellum torques in the dynamics of the first triad. Another contribution $\frac{\partial \mathbf{e}_{\mathbf{b}_1}^{\alpha}}{\partial t}$ comes from the torques on the cell body (Sect. 2.3 Eq. (15)). These two terms add up to the total change of $\mathbf{e}_1^\alpha$ as follows

$$ \frac{\partial \mathbf{e}_1^\alpha}{\partial t} = \frac{\partial \mathbf{e}_{\mathbf{f}_1}^\alpha}{\partial t} + \frac{\partial \mathbf{e}_{\mathbf{b}_1}^\alpha}{\partial t}. \quad (13) $$

## 2.3 Dynamics of the cell body

To obtain the dynamics of the whole bacterium, the equations of motion for the cell body are also needed and have to be solved simultaneously. For simplicity, we modeled the cell body with a sphere of radius $R_b = 1$ µm. Since MTB contain a magnetic chain, which acts as a compass needle, we defined a magnetic moment vector ($\mathbf{m}$) for the cell body (shown in Fig. 1a) with a magnitude of $m = 6.2 \times 10^{-16}$ Am$^2$ [40]. In each simulation the direction of $\mathbf{m}$ with respect to the body axis is fixed and its inclination is represented by an angle $\theta_m$. We should note that $\mathbf{m}$ is oriented along the body axis ($\theta_m = 0^\circ$) except for the section where we study the effect of the inclination of $\mathbf{m}$ with respect to body axis ($\theta_m \neq 0^\circ$) explicitly. In an external magnetic field $\mathbf{B}$, the equations of motion for the cell body are

$$ \frac{\partial \mathbf{r}_b}{\partial t} = \mu_b^{tt} [\mathbf{F}_1^K + \mathbf{F}_1^S + \mathbf{F}_0^{LJ}], \quad (14) $$

$$ \frac{\partial \mathbf{e}_{\mathbf{b}_1}^{\alpha}}{\partial t} = \mu_b^{rr} [\mathbf{m} \times \mathbf{B} + (\mathbf{r}_1 - \mathbf{r}_b) \times (\mathbf{F}_1^K + \mathbf{F}_1^S)] \times \mathbf{e}_1^{\alpha}, \quad (15) $$
$$ \frac{\partial \mathbf{m}}{\partial t} = \mu_b^{rr} [\mathbf{m} \times \mathbf{B} + (\mathbf{r}_1 - \mathbf{r}_b) \times (\mathbf{F}_1^K + \mathbf{F}_1^S) - T_m \mathbf{e}_1^3] \times \mathbf{m}. \quad (16) $$

Here $\mu_b^{tt} = \frac{1}{6\pi\eta R_b}$ and $\mu_b^{rr} = \frac{1}{8\pi\eta R_b^3}$ are the translational–rotational and rotational-rotational mobility coefficients of the cell body in which $\eta$ is the viscosity of water. Cross-mobility terms between the cell body and the flagellar particles were neglected.

Equation (14) describes the translational motion of the cell body. The force on the first flagellar particle ($\mathbf{F}_1^K + \mathbf{F}_1^S$) is transmitted directly to the cell body, together with the Lennard-Jones force ($\mathbf{F}_0^{LJ}$). Eq. (15) describes the contribution of the cell body rotation on the evolution of $\mathbf{e}_1^\alpha$, the triad associated with the bond between the cell body and the first flagellar particle. Its terms are the magnetic torque and the torque resulting from the forces acting on the first bead, respectively. As mentioned earlier, Eq. (12) and (15) add up to the total change of the first triad as given by Eq. (13).

The magnetosome chain has a fixed orientation in the cell body, so the magnetic moment $\mathbf{m}$ rotates rigidly with the cell body. To determine the orientation of the cell body, we need to include one additional equation of motion in the model that describes the dynamics of an arbitrary vector fixed on the cell body. We used the equation for the magnetic moment as the equation for finding the orientation of the cell body at each time-step. Eq. (16) represents the time evolution of $\mathbf{m}$ and the torque terms in its right-hand side are the magnetic torque, the torque arising from the force on the first particle and the counter torque to motor rotation ($-T_m \mathbf{e}_1^3$).

Since the first flagellum bead ($i=1$) is rigidly attached to the cell body, its position at each time-step is determined by $\mathbf{r}_1 = \mathbf{r}_b + R_b \mathbf{e}_1^3$.

## 2.4 Implementation of the model

The model was implemented using a self-written C++ program. The equations of motion (Eqs. (10) to (16)) were solved numerically using second-order Runge-Kutta algorithm [41, 42], specifically Heun method [43], to obtain the successive configurations of beads position and triads at each time-step knowing their previous configuration. For simplicity, we chose the magnetic field to be in z-direction and its magnitude in our simulation was varied from zero to hundred times the Earth magnetic field ($B_E = 50 \mu\text{T}$).

We used a time-step of $\Delta t = 4 \times 10^{-9}$s for all simulations and an initial configuration (helical shape) of the flagellum close to its equilibrium state. Nevertheless, equilibration takes of the order of $\times 10^7$ time steps. Therefore, the simulations were run for several $\times 10^7$ steps, with the total duration adjusted to the specific parameters. In particular, in the simulations with an inclined magnetic moment (Sect. 3.4), at least one full helical turn of the trajectories was simulated; thus, longer simulations were needed for large inclination angles, where the helix radius is large. Likewise for

<figure><img src="image_3.png" /><figcaption>Fig. 2 Alignment of the model magnetotactic bacterium in a magnetic field $B = 50 B_E$. The initial configuration of the bacterium is perpendicular to the magnetic field</figcaption></figure>

U-turn simulations (Sect. 3.3), longer runs were needed for low magnetic fields. In the U-turn scenario, the simulation was run to equilibrate for $5 \times 10^6$ time steps before the field was reversed.

# 3 Results and discussion

## 3.1 Swimming and alignment with magnetic field

The Stokesian Dynamics simulation framework allows us to address the mechanisms of swimming while a MTB interacts with a magnetic field. In particular, the model incorporates the magnetic moment of the cell body, which interacts with the field, the helical structure of the flagellum as well as its elasticity and the torque due to the flagellar motor. That the model describes both propulsion and magnetic steering is shown by the following simulation: Our bacterium is put in a uniform magnetic field, aligned initially perpendicular to that field. Running the simulation shows that the model bacterium aligns in the direction of the field by turning while swimming and then continues to swim along that direction. An exemplary reorientation trajectory is depicted in Fig. 2 by illustrating successive snapshots. The corresponding movie is also available in supplementary information (Video 1). Alignment is observed for arbitrary initial orientation and for different field strengths. The time needed for the complete alignment depends on the field strengths. We will discuss the dynamics of alignment in more detail in Sect. 3.3 below.

## 3.2 Impact of flagellar geometry on swimming

To characterize our model, we first studied the dependence of the swimming velocity of our model bacterium on geometrical features of the flagellum in the absence of a magnetic field. The geometry of the flagellar helix (as shown in Fig. 1) is defined by its radius $R$, pitch $P$, pitch angle $\alpha$ (with $\tan \alpha = \frac{2\pi R}{P}$) and contour length $L$. We simulated swimming for a wide range of $L$ and $\frac{\alpha}{R}$ and determined the average swimming velocity. The swimming velocity as a function of $L$ and $\frac{\alpha}{R}$ is plotted in Fig. 3a, b, respectively.

In Fig. 3a, the equilibrium curvature and torsion (and thus $R$ and $\alpha$) of the flagellum were kept constant and
<figure><img src="image_4.png" /><figcaption>Fig. 3 The swimming velocity of the bacterium as a function of geometrical features of flagellum. a Double-logarithmic plot of the velocity as a function of the flagellum</figcaption></figure>

the length of flagellum (L) was modulated by the number of beads N. The velocity is plotted as a function of L on a double-logarithmic scale. Both the average swimming velocity ($V_z$) and the tangential velocity ($V_t$) were calculated in each simulation. Figure 3a shows that the swimming velocity decreases with L which is in qualitative agreement with a theoretical prediction [1] for a rigid flagellum,

$$V \approx \alpha \frac{\xi_{\perp} - \xi_{\parallel}}{\xi_{\parallel}} \frac{\xi_r}{\xi_{\perp}} \frac{R_b^3}{RL} \Omega_m. \quad (17)$$

A linear fit to the log-log data (indicated by the dotted lines; the data for short contour lengths were omitted in the fit) shows a decrease as a power law for long flagella, however not with an exponent (slope in the log-log plot) of -1 as predicted by the theoretical expression, but with a slope of -0.9. The corresponding slope for the tangential velocity is -0.6 (the precise value of the slope depends on the time window of the fit).

This smaller exponent may be attributed to the approximations used in the theory. This formula is derived for the limit of a very long and rigid flagellum and for very small pitch angles ($\alpha \ll 1$). We also observe in our results that the slope in the log-log plot increases for longer flagella, so it is possible that it approaches -1 for even longer flagella (> 30 µm), which are not very realistic as well as computationally expensive. Moreover, for very long flagella, the assumption of a rigid flagellum becomes more and more unrealistic, as the flagellum length exceeds its persistence length. Indeed some simulations for long flagella show very small periodic modulations of the propulsion, likely due to the elasticity of the flagellum. We also note that Eq. (17), which is derived for the limit of long flagella, is not applicable to realistic flagellar lengths (~ 10 µm).

The simulation results for different sets of the parameters {$R, \alpha$} (Fig. 3a) do not show significant differences in this slope. The velocity as a function of $\alpha/R$

<figure><img src="image_5.png" /></figure>

contour length (L) and a linear fit to it, b Velocity as a function of the helix parameter $\alpha/R$ for different helix radii (R) and pitch angles (α) with a constant contour length of 6 µm

is shown in Fig. 3b for different helix radii. Here the contour length of the flagellum was kept constant at 6 µm and its pitch angle (α) was varied. The simulations were carried out for different helix radii in the range of R = 0.02–0.12 µm. The data show that the velocity increases with increasing $\alpha/R$, as predicted by Eq. (17). The increase is, however, not exactly linear as in that expression. Our simulations indicate that the dominant dependence on the parameters α and R is through the ratio $\alpha/R$. Therefore, we varied this parameter in the simulations shown in Fig. 3. However, not all dependence on these two parameters is captured by the $\alpha/R$-dependence, as seen by comparing the data for different R with the same ratio $\alpha/R$.

### 3.3 U-turn in response to magnetic field reversal

If the direction of the magnetic field is rapidly reversed, a swimming MTB performs a U-turn to realign with the magnetic field and swim in the opposite direction. The U-turn trajectory of MTB upon magnetic field reversal is used as a method to measure the magnetic moment of magnetotactic bacteria [44]. The magnetic moment of the bacterium can be calculated from the diameter of the U-turn or the time needed for it. Esquivel et. al [45] have investigated the dynamics of a magnetotactic cell following a sudden reversal of the magnetic field by solving its torque balance equation based on Bean model [46]. Since cells swim at low Reynolds number, inertial terms can be neglected and on average, random forces will create a zero net torque. Therefore, the torque exerted on the dipole by the magnetic field should be balanced by the viscous drag torque, i.e.

$$mB \sin \theta = \Gamma \frac{d\theta}{dt}, \quad (18)$$

where m is the dipole moment, B is the external magnetic field, $\Gamma$ is the rotational drag coefficient and $\theta$ is
the angle between the direction of the magnetic field and the magnetic dipole moment.

This differential equation can be solved for the time it takes for the dipole direction to change from $\theta_i$ to $\theta$ as $\ln(\tan \frac{\theta}{2}) = \frac{mB}{\Gamma}t + \ln(\tan \frac{\theta_i}{2})$. Due to thermal fluctuations the initial angle, $\theta_i$, is considered to be very small but not zero (for the onset of field reversal). Therefore, by considering $\theta_f = \frac{\pi}{2}$ (for the middle of U-turn), the U-turn time can be obtained as

$$ \tau_u = -2 \frac{\Gamma}{mB} \ln \left( \tan \frac{\theta_i}{2} \right). \qquad (19) $$

Moreover, by integrating the component of velocity along U-turn opening ($\int_0^\infty v \sin \theta dt$, where $\sin \theta = \frac{\Gamma}{mB} \frac{d\theta}{dt}$) over time, the U-turn diameter can also be calculated as

$$ D_u = \frac{\pi v \Gamma}{mB}. \qquad (20) $$

For our model MTB, we simulated the U-turn experiment for a wide range of magnetic field strength, starting with a bacterium swimming along the direction of the field and instantaneously reversing the field direction. Figure 4a and b shows trajectories and snapshots for two different magnetic field strengths, of $10B_E$ and $500B_E$, respectively. Here, the time of the magnetic field reversal is defined as $t = 0$s.

With small fields (Fig. 4a), the bacterium aligns with the magnetic field gradually, and thus, an obvious and wide U-turn is observed in its trajectory. However, with large fields (Fig. 4b), the U-turn is very sharp and tight. Moreover, we observe that with large fields, the flagellum undergoes a pronounced bend or kink during the U-turn, as shown in Fig. 4b at $t = 9$ms and $t = 30$ms. This kink formation can be attributed to the very fast reorientation of the magnetic moment for alignment with the field. Since the flagellum is long and because of its elasticity, its distant segments can not immediately follow the very fast reorientation of the cell body while the close segments reorient due to their direct attachment to the cell body.

To make this reasoning more quantitative, we determined two characteristic time scales related to the U-turn, the time $\tau_{um}$ after which the magnetic moment is re-aligned with the magnetic field and the time $\tau_{ut}$ after which the swimming direction (and thus the trajectory) is re-aligned with the field direction. The two time scales as obtained from the simulations are plotted as functions of the magnetic field in Fig. 5a.

Figure 5a shows that U-turn time decreases with magnetic field. For fields up to $120B_E$, both characteristic times have almost the same value, but for stronger fields $\tau_{ut}$ exceeds $\tau_{um}$ remarkably. This confirms our previous interpretation that in these strong fields, the magnetic moment aligns with the field very rapidly, while the flagellum requires more time to align with the cell body's new direction. In this range of the magnetic field ($B > 120B_E$), $\tau_{ut}$ plateaus, in agreement with the expectation for a reorientation of the flagellum and

thus the trajectory driven by the flagellum's mechanics rather than by the magnetic torque. Correspondingly, we see that this plateau depends on the length of the flagellum as shown in Fig. 5b. The longer the flagellum, the slower the reorientation (increasing $\tau_{ut}$). $\tau_{um}$, by contrast is not affected by flagellum length. Since the slower of the two processes determines the overall U-turn time, we can define a crossover field strength ($B^*$) between the regimes where magnetic reorientation determines the U-turn time $\tau_u^*$ and a regime where the flagellar reorientation is dominant. This crossover field and the corresponding U-turn time are plotted as functions of flagellum length in the inset of Fig. 5b. Figure 5c shows the U-turn time for strong fields as a function of the flagellum length in more detail. In addition to the overall increase in the $\tau_{ut}$ with the length, we observe an oscillatory pattern with a period that approximately corresponds to the contour length of one turn of the helix (2.54 µm for the equilibrium curvature and torsion of $\{\kappa_e, \tau_e\} = \{1.3, -2.1\} \mu m^{-1}$).

The analytical results for the U-turn time and U-turn radius discussed above consider only the magnetic reorientation and are thus expected to be applicable to regime for weak fields. To apply the analytical expressions to our simulation data, we use the rotational friction coefficient $\Gamma$ of the bacterium as a fit parameter. We know the friction coefficient of the cell body, but this might be modulated by the presence of the flagellum. Fitting the U-turn time with the expression from eq.(19), we find good agreement with the theoretical prediction for $\Gamma/(8\pi\eta R^3) = 1.12$, i.e., for a friction coefficient of the whole bacterium that is increased compared to that of the cell body alone ($8\pi\eta R^3$).

To test this result, we also determined the friction coefficient directly but rotating the bacterium with an external torque. For this simulations, we turned off the motor torque and rotated the bacterium with an external torque applied to the cell body in a direction perpendicular to the motor axis. We then measured the angular velocity of the cell body to determine the friction coefficient. The inset in Fig. 5a shows the resulting estimate of the rotational friction coefficient of the whole bacterium as a function of its flagellum length. As expected, the drag coefficient is increased in the presence of a flagellum. The dependence on the flagellum length is relatively weak. The minimum in the friction coefficient coincides with the length where the flagellum has just one full pitch. The numerical value for a length of $L = 12 \mu m$ is 1.19, close to what we have obtained from the fit to the U-turn time.

Figure 5d shows the diameter of the U-turn (scaled by the velocity) as a function of magnetic field for flagellum lengths of $L = 4, 6, 8 \mu m$. The normalization by velocity is important because in populations of MTBs, the velocity can vary from cell to cell. Moreover, from an experimental point of view, the U-turn diameter and the velocity are more easily accessible than the U-turn time. According to Fig. 5d, the U-turn diameter also exhibits the two regimes discussed for the U-turn time. Fitting the data with Eq.(20) results in a friction coefficient of $\Gamma/(8\pi\eta R^3) = 1.23$, slightly larger than the
<figure><img src="image_6.png" /><figcaption>Fig. 4 U-turn trajectory illustrated with successive snapshots in case of **a** a weak magnetic field ($10B_E$) and **b** a strong magnetic field ($500B_E$) for a bacterium with a flagellum of length $L = 12~\mu\text{m}$</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>Fig. 5 **a** Characteristic time scales of U-turn after field reversal as a function of the magnetic field strength. The times $\tau_{um}$ (blue) and $\tau_{ut}$ (red) measure the time it takes the bacterium to turn its magnetic moments (and thus the cell body) and the time it takes to reverse the swimming direction, respectively. The solid line is a fit to the theoretical expression obtained from torque balance with $\frac{\Gamma}{8\pi\eta R^3} = 1.12$ and $\theta_i = 0.36$. The flagellum length is $L = 12~\mu\text{m}$. Inset: The friction coefficient of the whole bacterium as a function of flagellum length obtained from direct simulations in</figcaption></figure>

which the bacterium was rotated with an external torque. **b** $\tau_{ut}$ and $\tau_{um}$ as a function of magnetic field for different flagellum lengths of $L = 4, 6, 8~\mu\text{m}$. Inset: The cross-over field strength ($B^*$) and U-turn time ($\tau_u^*$) between the two regimes dominated by different characteristic times. **c** $\tau_{ut}$ as a function of flagellum length for different strengths of the magnetic field. **d** The U-turn diameter (scaled by swimming velocity) as a function of the magnetic field strength for flagellum lengths of $L = 4, 6, 8~\mu\text{m}$. The solid line is a fit to the theoretical expression with $\frac{\Gamma}{8\pi\eta R^3} = 1.23$
result from fitting the U-turn time, but in good agree-
ment with the direct measurement of the friction coeffi-
cient. The discrepancy appears to be related to the fact
that the theory assumes entirely planar motion, while
in the simulations, the bacterium makes small excur-
sions in the direction perpendicular to the main plane
of motion during the U-turn.

3.4 The effect of magnetic moment inclination on
3D motility pattern

In this section, we investigated how the swimming tra-
jectory of our model bacterium is modified by an incli-
nation of the magnetic moment with respect to the
propulsion axis. To that end, we modulated the angle
θm in our simulations (0° ≤ θm < 90°). Small incli-
nation angles are expected even for bacteria where the
nominal direction of the magnetic moment is along the
propulsion axis, as the precision of control is necessarily
finite. However, some species including Magnetococcus
marinus (MC-1) and Magnetococcus massalia (MO-1)
do show rather large angles. These do, however, not
have single flagellum as our model bacterium here and
for MC-1, we proposed recently that it swims with a dif-
ferent body orientation that reduces this angle to some
extent [9]. Nevertheless, it is interesting to quantita-
tively analyze the effect of such inclination.

When the magnetic moment is parallel to the flagel-
lum axis (θm = 0°), the case studied so far, the overall
motion is observed to be approximately on a straight
line. Closer inspection shows a wiggling of the cell body
due to a weak rotational motion around that line, so
that the trajectory is actually helical with a very small
helix radius, smaller than the size of the cell body, as
shown in Fig. 6a, where the radius is about 0.1 µm.
The incidence of such wiggling has been noted in ear-
lier work and reflects the fact that the axis of propulsion
is not exactly parallel to the axis of the cell body or, in
our case, the magnetic moment [47–49]: Due to the non-
integer number of turns in the helix of the flagellum, the
thrust force always has a component normal to the axis
of the flagellum and leads to an off-axis torque. Propul-
sion together with this off-axis torque result in wiggling.
A comprehensive study that characterized the wiggling
in microorganism trajectories was presented by Hyon
et al. [49].

When the magnetic moment is not parallel to the
flagellum (θm ≠ 0°), we observe double helical trajecto-
ries in the presence of a magnetic field. The small helix
due to the wiggling is now bent and follows the path
of a larger helix due to the precession of the magnetic
moment with a radius in the order of several times that
of the body size. An exemplary double helical trajectory
is depicted in Fig. 6a for θm = 45° and B = 500B<sub>E</sub>.
The radii of the small and large helices are represented
with r<sub>H</sub> and R<sub>H</sub>, respectively. The axis of the large
helix is parallel to the direction of the magnetic field.

We characterized the double helical trajectories by
the swimming velocity Vz along the direction of the
magnetic field (which we take to be in the z direction),

and by the radius, pitch, and period of the large helix
(R<sub>H</sub>, P<sub>H</sub> and T<sub>H</sub>, respectively). We extracted these
parameters from the trajectories obtained for different
inclination angles in the range of θ<sub>m</sub> = 0–88° and dif-
ferent strengths of the magnetic field (B = 0–1000B<sub>E</sub>).
The results for the helix radius (R<sub>H</sub>) and the swimming
velocity (V<sub>z</sub>) are plotted in Fig. 6b, c as functions of the
inclination angle θ<sub>m</sub> for different strengths of the mag-
netic field. The corresponding plots for the pitch (P<sub>H</sub>)
and the period (T<sub>H</sub>) are shown as insets in Fig. 6b, c.

Our results show that the swimming velocity decreases
with increasing inclination angle θm (Fig. 6c), while the
radius of the helix increases (Fig. 6b). Two regimes for
low and high magnetic fields can be distinguished. At
low magnetic fields (< 50B<sub>E</sub>), the large helix is not
easily detectable since it is only very slightly greater
than the small helix (only if θm approach 90°). In this
regime, the magnetic torque is not strong enough to
cause considerable precession of the bacterium. As a
consequence, the velocity is almost not affected by the
inclination angle. For small inclinations (θm < 20°),
the swimming velocity is even slightly increased by the
presence of the magnetic field, likely because the wig-
gling motion is slightly suppressed. We expect that this
effect is more pronounced in the presence of thermal
fluctuations, which provide stronger perturbations of
the alignment with the field, which in turn restrains
the bacterium to swim only in one dimension [16].

The large helix is clearly detectable for B > 120B<sub>E</sub>
and θ<sub>m</sub> > 45°. In this case, the radius of large helix
is about several times the body size and the veloc-
ity decreases strongly with the inclination angle. The
strong magnetic torque forces the magnetic moment
into a strong rotation toward the magnetic field and
hence a greater deviation from the propulsion axis.

These observations indicate that small inclination
angles (< 30°) due to imprecise orientation of the mag-
netic moments with respect to the propulsion axis do
not affect the swimming of MTB in naturally occurring
fields much. Even inclination angles of 45° would not
strongly hinder the motion. However, inclination angles
exceeding 45° as observed in strains MC-1 and MO-1
would put those cells at a disadvantage. This obser-
vation indicates that a different swimming mode with
different orientation of the magnetic moments as pro-
posed for MC-1 [9] is beneficial for these strains.

4 Conclusion

In this work, we developed a detailed model for the
propulsion and orientation of magnetotactic bacteria
(MTB) based on Stokesian dynamics and investigated
their swimming in the presence of an external mag-
netic field. There is a great variety of MTB in nature
that differ in the structure of their motility apparatus
(e.g., single- or multiflagellated with different localiza-
tions of the flagella on the cell surface) and the organi-
zation of their magnetosomes, the organelles providing
them with a magnetic moment (single chain, multiple
<figure><img src="image_8.png" /><figcaption>**Fig. 6** a Side view of the (double-)helical trajectories observed for the swimming of a model bacterium without and with inclined magnetic moment (at an angle $\theta_m = 45^\circ$ relative to the flagellar axis and for $B = 500B_E$). The direction of the magnetic field is upward. **b** Radius of the large</figcaption></figure>

chains, different chain orientation) [15]. Here, we studied a model MTB with the simplest motility apparatus, a single flagellum, and varied its parameter as well as the orientation of the magnetic moment and the field strength (we note that a variation in the strength of the magnetic moment has the same effect as the corresponding change in the magnetic field strength, as only their product enters the orientational dynamics). In the simulation, we have the possibility to investigate a wide range of values for these parameters, some of which are not biologically relevant but worth studying as they might be applicable for predicting swimming behavior in the case of artificial microswimmers.

After an initial characterization of our minimal model bacterium, we considered two scenarios that are of interest for the experimental study of magnetotactic bacteria. On the one hand, we simulated U-turns upon reversal of the magnetic field direction, a scenario often used to measure the magnetic moment of the bacterium [44]. On the other hand, we studied the three-dimensional shape of swimming trajectories in a magnetic field, a quantity that has become accessible with the introduction of 3d tracking methods [50] and that provides information about swimming mechanisms as shown recently for the *Magnetococcus marinus* MC-1 [9].

With respect to U-turns, our simulations qualitatively agree with a theoretical description based on a moving point dipole (at least for the weak fields relevant to experiments). However, they also point toward a difficulty of using U-turns to determine the magnetic moment of a bacterium, namely that the rotational friction coefficient must be known. Typically, a friction coefficient is estimated for the geometry of the cell body, but this may be a cause of inaccuracy as indicated by our simulations. For strong fields, the U-turns in our simulation do not follow the prediction of the point dipole theory, because the relaxation of the flagellum (which determined the direction of propulsion) is slower than the relaxation of the magnetic moment. Thus, our results suggest that using stronger fields than typically used in experiments may result in novel behaviors such

helix, and **c** net swimming velocity in the direction of the magnetic field as functions of the angle between the magnetic moment and the flagellar propulsion axis ($\theta_m$) for different field strengths. The insets in **b** and **c** show the corresponding pitch and period of the large helix, respectively

as formation of kink in flagellum during a U-turn as well as the observation of two different characteristic U-turn times. These should provide insights into the elastic response of the flagellum. However, we also note that with these strong deformations of the flagellum, other flagellar properties may come into play, in particular polymorphic transitions [20, 34, 51].

With respect to swimming trajectories, we observed bacteria with a magnetic moment at a non-zero angle relative to the flagellar axis to swim on double helical trajectories, spiralling in the direction defined by the magnetic field. Helical trajectories are seen in various microorganisms and are of interest for different swimming and sensing strategies [9, 18, 30, 52–54]. When we characterized the swimming velocity on these trajectories, we found that small inclination angles ($\theta_m < 20^\circ$) are not a significant disadvantage for swimming in realistic field strength. Large inclinations, as seen in some MTB such as strain MC-1, would be a problem, as the net velocity in the direction of the field is considerably reduced by the helical motion. However, our recent study [9] suggests that strain MC-1 avoids this problem by using a different swimming mechanism (cooperative pushing and pulling flagella), which reduces the inclination of the magnetic moment with respect to the propulsion axis.

The simulation method presented here is quite general and has indeed been used for different types of swimmers, specifically *Escherichia coli* with multiple flagella all over the cell body [33], *Shewanella putrefaciens* with a single flagellum [8] and *Magnetococcus marinus* with two flagella and a magnetic moment [9]. Since the model includes both the elasticity of the flagellum and hydrodynamic interactions, it is particularly suited to address questions of re-orientation as well as flagellar synchronization in multi-flagellated microorganisms [33, 55, 56].

Our model can also be used to characterize the swimming of MTB for their potential applications in MTB-inspired microrobots. For example, having a precise control on the helical trajectory might be useful in enhancing the bacterium sensing of external stimuli
[49]. One can also extend the model to include a cargo to be transported by the bacterium and study its effect on the swimming bacterium. In general, our results here and those in our previous work [9] show that Stokesian dynamics provides a useful method to address the diverse swimming mechanisms of magnetotactic bacteria.

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1140/epje/s10189-021-00038-5.

**Acknowledgements** The authors would like to thank Holger Stark and Tapan Chandra Adhyapak for helpful discussions and their support with setting up the simulations. The research leading to these results was supported by the Max Planck Society and by Deutsche Forschungsgemeinschaft (DFG) within the priority program on microswimmers SPP 1726 (Grants No. KL 818/2-2 and FA 835/7-2 to SK and DF). Further, SM was supported by Deutscher Akademiischer Austauschdienst, DAAD (Grant No. 57314018). Part of the simulations were run on the Goegrid cluster at the University of Göttingen, which is supported by the DFG (Grant INST 186/1353-1 FUGG) and MWK Niedersachsen (Grant No. 45-10-19-F-02).

**Funding** Open Access funding enabled and organized by Projekt DEAL.

**Author contribution statement**

SM, DF and SK designed research. SM wrote the code and performed simulations. SM and SK analyzed the data, and all authors wrote the manuscript.

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

**References**

1. E. Lauga, T.R. Powers, The hydrodynamics of swimming microorganisms. Rep. Prog. Phys. **72**(9), 096601 (2009)

2. J. Elgeti, R.G. Winkler, G. Gompper, Physics of microswimmers–single particle motion and collective

behavior: a review. Rep. Prog. Phys. **78**(5), 056601 (2015)

3. C. Bechinger, R. Di Leonardo, H. Löwen, C. Reichhardt, G. Volpe, G. Volpe, Active particles in complex and crowded environments. Rev. Mod. Phys. **88**, 045006 (2016)

4. H.C. Berg, *E. Coli in Motion in Motion* (Springer, Berlin, 2008)

5. N.C. Darnton, L. Turner, S. Rojevsky, H.C. Berg, On torque and tumbling in swimming *Escherichia coli*. J. Bacteriol. **189**(5), 1756–1764 (2007)

6. B.M. Friedrich, F. Jülicher, Chemotaxis of sperm cells. Proc. Natl. Acad. Sci. **104**(33), 13256–13261 (2007)

7. K. Son, J.S. Guasto, R. Stocker, Bacteria can exploit a flagellar buckling instability to change direction. Nat. Phys. **9**(8), 494 (2013)

8. M.J. Kühn, F.K. Schmidt, B. Eckhardt, K.M. Thormann, Bacteria exploit a polymorphic instability of the flagellar filament to escape from traps. Proc. Natl. Acad. Sci. **114**(24), 6340–6345 (2017)

9. K. Bente, S. Mohammadinejad, M.A. Charsooghi, F. Bachmann, A. Codutti, C.T. Lefèvre, S. Klumpp, D. Faivre, High-speed motility originates from cooperatively pushing and pulling flagella bundles in bilophotrichous bacteria. Elife **9**, e47551 (2020)

10. D.A. Bazylinski, R.B. Frankel, Magnetosome formation in prokaryotes. Nat. Rev. Microbiol. **2**(3), 217–230 (2004)

11. S. Klumpp, C.T. Lefèvre, M. Bennet, D. Faivre, Swimming with magnets: from biological organisms to synthetic devices. Phys. Rep. **789**, 1–54 (2019)

12. D. Faivre, D. Schuler, Magnetotactic bacteria and magnetosomes. Chem. Rev. **108**(11), 4875–4898 (2008)

13. A. Komeili, Molecular mechanisms of compartmentalization and biomineralization in magnetotactic bacteria. FEMS Microbiol. Rev. **36**(1), 232–255 (2012)

14. R. Blakemore, Magnetotactic bacteria. Science **190**(4212), 377–379 (1975)

15. C.T. Lefèvre, M. Bennet, L. Landau, P. Vach, D. Pignol, D.A. Bazylinski, R.B. Frankel, S. Klumpp, D. Faivre, Diversity of magneto-aerotactic behaviors and oxygen sensing mechanisms in cultured magnetotactic bacteria. Biophys. J. **107**(2), 527–538 (2014)

16. F. Popp, J.P. Armitage, D. Schüler, Polarity of bacterial magnetotaxis is controlled by aerotaxis through a common sensory pathway. Nat. Commun. **5**, 5398 (2014)

17. O. Felfoul, M. Mohammadi, S. Taherkhani, D.D. Lanauze, Y.Z. Xu, D. Loghin, S. Essa, S. Jancik, D. Houle, M. Lafleur et al., Magneto-aerotactic bacteria deliver drug-containing nanoliposomes to tumour hypoxic regions. Nat. Nanotechnol. **11**(11), 941 (2016)

18. B.-W. Park, J. Zhuang, O. Yasa, M. Sitti, Multi-functional bacteria-driven microswimmers for targeted active drug delivery. ACS Nano **11**(9), 8910–8923 (2017)

19. G. Santomauro, A.V. Singh, B.-W. Park, M. Mohammadrahimi, P. Erkoc, E. Goering, G. Schütz, M. Sitti, J. Bill, Incorporation of terbium into a microalga leads to magnetotactic swimmers. Adv. Biosyst. **2**(12), 1800039 (2018)

20. R. Vogel, H. Stark, Rotation-induced polymorphic transitions in bacterial flagella. Phys. Rev. Lett. **110**(15), 158104 (2013)
21. P.J. Vach, P. Fratzl, S. Klumpp, D. Faivre, Fast magnetic micropropellers with random shapes. Nano Lett. 15(10), 7064–7070 (2015)

22. A. Codutti, F. Bachmann, D. Faivre, S. Klumpp, Bead-based hydrodynamic simulations of rigid magnetic micropropellers. Frontiers Robot. AI 5, 109 (2018)

23. D. Murat, M. Hérisse, L. Espinosa, A. Bossa, F. Alberto, W. Long-Fei, Opposite and coordinated rotation of amphitrichous flagella governs oriented swimming and reversals in a magnetotactic spirillum. J. Bacteriol. 197, 3275–3282 (2015)

24. M. Reufer, R. Besseling, J. Schwarz-Linek, V.A. Martínez, A.N. Morozov, J. Arlt, D. Trubitsyn, F.B. Ward, W.C.K. Poon, Switching of swimming modes in magnetospirillum gryphiswaldense. Biophys. J. 106(1), 37–46 (2014)

25. K. Erglis, Q. Wen, V. Ose, A. Zeltins, A. Sharipo, P.A. Janmey, A. Cēbers, Dynamics of magnetotactic bacteria in a rotating magnetic field. Biophys. J. 93(4), 1402–1412 (2007)

26. G.-J. Liao, C.K. Hall, S.H.L. Klapp, Dynamical self-assembly of dipolar active Brownian particles in two dimensions. Soft Matter 16(9), 2208–2223 (2020)

27. A. Codutti, K. Bente, D. Faivre, S. Klumpp, Chemotaxis in external fields: simulations for active magnetic biological matter. PLoS Comput. Biol. 15(12), e1007548 (2019)

28. D. Kong, W. Lin, Y. Pan, K. Zhang, Swimming motion of rod-shaped magnetotactic bacteria: the effects of shape and growing magnetic moment. Frontiers Microbiol. 5, 8 (2014)

29. F.S. Nogueiral, H.G.P.L. de Barros, Study of the motion of magnetotactic bacteria. Eur. Biophys. J. 24(1), 13–21 (1995)

30. C. Yang, C. Chen, Q. Ma, W. Longfei, T. Song, Dynamic model and motion mechanism of magnetotactic bacteria with two lateral flagellar bundles. J. Bionic Eng. 9(2), 200–210 (2012)

31. J.F. Brady, G. Bossis, Stokesian dynamics. Ann. Rev. Fluid Mech. 20(1), 111–157 (1988)

32. R. Vogel, H. Stark, Motor-driven bacterial flagella and buckling instabilities. Eur. Phys. J. E 35(2), 15 (2012)

33. T.C. Adhyapak, H. Stark, Zipping and entanglement in flagellar bundle of *E. coli*: role of motile cell body. Phys. Rev. E 92(5), 052701 (2015)

34. T.C. Adhyapak, H. Stark, Dynamics of a bacterial flagellum under reverse rotation. Soft Matter 12(25), 5621–5629 (2016)

35. J.K.G. Dhont, *An Introduction to Dynamics of Colloids*, vol. 2 (Elsevier, Amsterdam, 1996).

36. H. Yamakawa, T. Yoshizaki, *Helical Wormlike Chains in Polymer Solutions*, vol. 1 (Springer, Berlin, 1997).

37. L.D. Landau, L.P. Pitaevskii, A.M. Kosevich, E.M. Lifshitz, *Theory of Elasticity*, vol. 7, 7th edn. (Elsevier, Amsterdam, 2012).

38. M. Reichert, *Hydrodynamic interactions in colloidal and biological systems*. PhD thesis, University of Konstanz, Faculty of Sciences, Department of Physics (2006)

39. J. Rotne, S. Prager, Variational treatment of hydrodynamic interaction in polymers. J. Chem. Phys. 50(11), 4831–4837 (1969)

40. S. Klumpp, D. Faivre, Magnetotactic bacteria. Eur. Phys. J. Spec. Top. 225(11–12), 2173–2188 (2016)

41. G. Sewell, *The Numerical Solution of Ordinary and Partial Differential Equations*, vol. 75 (Wiley, New York, 2005).

42. W.H. Press, B.P. Flannery, S.A. Teukolsky, W.T. Vetterling, *Numerical Recipes in C: The Art of Scientific Computing* (Cambridge University Press, New York, 1988).

43. A. Greiner, W. Strittmatter, J. Honerkamp, Numerical integration of stochastic differential equations. J. Stat. Phys. 51(1–2), 95–108 (1988)

44. R. Nadkarni, S. Barkley, C. Fradin, A comparison of methods to measure the magnetic moment of magnetotactic bacteria through analysis of their trajectories in external magnetic fields. PLoS ONE 8(12), e82064 (2013)

45. D.M.S. Esquivel, H.G.P.L. De Barros, Motion of magnetotactic microorganisms. J. Exp. Biol. 121(1), 153–163 (1986)

46. R.B. Frankel, Magnetic guidance of organisms. Ann. Rev. Biophys. Bioeng. 13(1), 85–103 (1984)

47. J.B. Keller, S.I. Rubinow, Swimming of flagellated microorganisms. Biophys. J. 16(2 Pt 1), 151 (1976)

48. H.C. Berg, D.A. Brown et al., Chemotaxis in *Escherichia coli* analysed by three-dimensional tracking. Nature 239(5374), 500–504 (1972)

49. Y. Hyon, T.R. Powers, R. Stocker, H.C. Fu et al., The wiggling trajectories of bacteria. J. Fluid Mech. 705, 58–76 (2012)

50. K.M. Taute, S. Gude, S.J. Tans, T.S. Shimizu, High-throughput 3d tracking of bacteria on a standard phase contrast microscope. Nat. Commun. 6(1), 1–9 (2015)

51. L. Turner, W.S. Ryu, H.C. Berg, Real-time imaging of fluorescent flagellar filaments. J. Bacteriol. 182(10), 2793–2801 (2000)

52. T. Fenchel, N. Blackburn, Motile chemosensory behaviour of phagotrophic protists: mechanisms for and efficiency in congregating at food patches. Protist 150(3), 325–336 (1999)

53. H.C. Crenshaw, Orientation by helical motion—iii. Microorganisms can orient to stimuli by changing the direction of their rotational velocity. Bull. Math. Biol. 55(1), 231–255 (1993)

54. Y. Pan, W. Lin, J. Li, W. Wenfang, L. Tian, C. Deng, Q. Liu, R. Zhu, M. Winklhofer, N. Petersen, Reduced efficiency of magnetotaxis in magnetotactic coccoid bacteria in higher than geomagnetic fields. Biophys. J. 97(4), 986–991 (2009)

55. S.Y. Reigh, R.G. Winkler, G. Gompper, Synchronization and bundling of anchored bacterial flagella. Soft Matter 8(16), 4363–4372 (2012)

56. D.R. Brumley, K.Y. Wan, M. Polin, R.E. Goldstein, Flagellar synchronization through direct hydrodynamic interactions. Elife 3, e02750 (2014)