Lorentz forces induce inhomogeneity and fluxes in active systems

We consider the nonequilibrium dynamics of a charged active Brownian particle in the presence of a space dependent magnetic field. It has recently been shown that the Lorentz force induces a particle flux perpendicular to density gradients, thus preventing a diffusive description of the dynamics. Whereas a passive system will eventually relax to an equilibrium state, unaffected by the magnetic field, an active system subject to a spatially varying Lorentz force settles into a nonequilibrium steady state characterized by an inhomogeneous density and divergence-free bulk fluxes. A macroscopic flux of charged active particles is induced by the gradient of the magnetic field only and does not require additional symmetric breaking such as density or potential gradients. This stands in marked contrast to similar phenomena in condensed matter such as the classical Hall effect. In a confined geometry we observe circulating fluxes, which can be reversed by inverting the direction of the magnetic field. Our theoretical approach, based on coarse-graining of the FokkerPlanck equation, yields analytical results for the density, fluxes, and polarization in the steady state, all of which are validated by direct computer simulation. We demonstrate that passive tracer particles can be used to measure the essential effects of the Lorentz force on the active particle bath, and we discuss under which conditions the effects of the flux could be observed experimentally.


I. INTRODUCTION
The Lorentz force modifies the trajectory of a charged particle without performing work on it. In general the Lorentz force F L acts on a moving charge q as where v denotes the velocity vector of the charge and B is the applied magnetic field. The Lorentz force gives rise to important phenomena in solid-state physics such as the classical Hall effect, in which electrical conductors develop a voltage difference transverse to the direction of the current if a magnetic field is applied perpendicular to the direction of the flow of charge carriers [1]. Note that a macroscopic flux of charge carriers induced by Lorentz forces requires an explicitly broken symmetry that gives rise to a macroscopic velocity vector in addition to the symmetry breaking due to the magnetic field vector. Simply speaking the classical Hall effect leads to a deviation of an existing macroscopic flux. The main finding of this work is that in active systems driven by internal forces a macroscopic flux emerges from a flux-free system in the presence of a magnetic field gradient. While the effect of the Lorentz force on the properties of materials has been thoroughly studied in the context of solid-state physics [2], much less is known about its influence on soft-matter systems, which are dominated by overdamped diffusive dynamics. Part of the reason is that Lorentz forces may appear to be too weak to affect slowly moving colloidal particles significantly. However, for Lorentz force to affect the motion of the particle, it is * sharma@ipfdd.de not its absolute value that is important but how its magnitude compares with the frictional force coming from the solvent. Since both the forces depend on velocity, the strength of Lorentz force relative to frictional force can be quantified and shown to be significant for highly charged mesoscopic particles in strong magnetic fields.
Before coming back to these points it is important to understand how the Lorentz force affects diffusive dynamics and how the diffusion equation can be extended to incorporate nonconservative velocity-dependent forces. It is known that as a consequence of the Lorentz force the corresponding Smoluchowski equation for the probability distribution acquires a tensorial diffusion coefficient [3], the components of which are determined by the applied magnetic field and the friction coefficient. Recently it has been shown that this tensor has an antisymmetric part, which gives rise to additional fluxes in the direction perpendicular to density gradients [4,5]. In a passive system, these fluxes affect the dynamics of a system as it relaxes from an initial nonequilibrium state to an equilibrium state. The equilibrium properties, however, are unaffected by the Lorentz force. A known consequence of the Lorentz force is a reduction of the diffusion coefficient in the plane perpendicular to the magnetic field, which can be understood by the deflection of the diffusing particle in the direction perpendicular to its trajectory according to the Eq. (1). This does not affect the equilibrium distribution of a system of passive particles.
A completely different situation arises when a driven system is subjected to an external magnetic field. A particularly interesting example of a driven system is active matter, which is ubiquitous in biology; examples include molecular motors inside cells [6] and microscopic swimmers such as bacteria and algae [7]. Besides biological active matter, there exists synthetic active matter, most notably Janus particles, which are 'two-faced' particles with only one of them functionalized. When the functionalized half of the particle is activated, chemically or under illumination, the particle is propelled forward [8][9][10][11][12]. In contrast to externally driven systems, active matter has the hallmark feature that it is out of equilibrium without an externally broken symmetry. Much progress has been made in understanding the properties of active matter by using active Brownian particles (ABPs) as a model system, which captures the essential features of an internally driven system in a minimalistic fashion. ABPs violate time-reversal symmetry by consuming fuel to generate motion, often referred to as self-propulsion. Besides the application to biological systems, active matter serves as a paradigm to study the effect of broken timereversal symmetry and nonequilibrium steady states in general [7,[13][14][15].
The absence of time-reversal symmetry can be exploited to obtain active fluxes by breaking the spatial symmetry of the system. For instance, sliding motion of active particles along an asymmetric wall [17][18][19][20][21] can give rise to fluxes along a boundary. Similarly, in systems without walls or obstacles, fluxes can be induced by making the activity spatially inhomogeneous [22]. Explicitly including alignment interactions, such as torque on an active particle, can also give rise to fluxes in a spatially inhomogeneous activity field as has been shown experimentally in Ref. [23]. Whereas in these studies activity has been time independent, one can induce fluxes by making the swim speed space and time dependent, as has been shown in theoretical and numerical studies [24][25][26][27][28][29]. In these systems the active particles 'surf' on an activity wave resulting in active fluxes.
In this work, we show that the Lorentz force can induce macroscopic effects in a system of ABPs such as a magnetic-field dependent density distribution and fluxpatterns. Note that, in contrast to externally driven system (for example, charge carriers in conductors under an applied voltage), ABPs do not have a preferential direction of motion. Thus, inhomogeneities are solely controlled by the spatial variation of the magnetic field. Whereas interaction between particles were necessary to induce bulk fluxes in a system with a space-dependent swim speed which is constant in time [22], this is not the case for the system studied here. Moreover, with Lorentz forces, the fluxes can be reversed by reversing the direction of the magnetic field.
For slowly moving particles, Lorentz forces can be extremely small. However, the effect of the Lorentz force on the motion of a particle is determined by ratio of the magnitudes of the Lorentz force and the frictional force. This ratio, which we refer to as the diffusive Hall constant, may become significant for highly charged particles in strong magnetic fields. We note that ABPs with magnetic interaction have been studied before; however, only by including a magnetic dipole moment without considering the Lorentz force. [30][31][32][33][34].
As an example of the effect of the Lorentz force on charged active particles we show the results of Brownian dynamics simulations where the ABPs are subjected to an inhomogeneous magnetic field in Fig. 1. The Gaussian shaped magnetic field is radially symmetric. The two signatures of the Lorentz force on ABPs are clearly visible: Circulating fluxes perpendicular to the gradient of the magnetic field and an inhomogeneous density distribution of the particles. The fluxes observed in ABPs under Lorentz force are reminiscent of the Corbino effect in electrical conductors [16]. Whereas a radial flux of charges is necessary for the Corbino effect in conductors, this is not the case for ABPs which do not require explicit symmetry breaking in the system except the spatial inhomogeneity of the magnetic field.
On a microscopic level, fluxes in our system are related to the polarization of the active particles. Polar ordering can emerge in active systems without aligning interactions. This has been shown in previous studies, for instance, under gravity [35,36] and inhomogeneous activity [22,37]. As in these studies, the origin of polarization in our system is kinetic in nature. However, in contrast to previous studies, the particles in our system orient not only along the gradient of the applied field but also perpendicular to it. Related to this polar order are active fluxes in the same direction.
A common starting point of the theoretical description of ABPs is the Langevin equation. Most models of ABP ignore inertial effects and are based on overdamped equations of motion [27,[37][38][39][40][41]. Due to the nonconservative nature of the Lorentz force, the overdamped equation of motion cannot be obtained by simply setting the mass of the particles to zero [4,42]. Though there exists a nontrivial limiting procedure [43][44][45] that yields the small-mass limit of the (velocity) Langevin equation, the resulting overdamped equation is not suitable for determining velocity dependent variables such as flux and entropy [4,5,[46][47][48] despite the fact that it captures the position statistics accurately. In this work, we have avoided these problems by performing simulations using the (velocity) Langevin equation with a finite but small mass. Our theoretical approach is based on coarsegraining of the Fokker-Planck equation using a gradient expansion [39,49] , which yields an equation for the flux and density distribution as a function of time and the position variable alone. This equation also yields the space-dependent polarization in the system by identifying the corresponding contribution from the orientational degrees of freedom. We further show that one can improve the accuracy of the results by combining the coarsegraining approach with linear-response theory. The results of this method match those of Brownian dynamics simulations.
Our calculations lead to a straightforward estimation of the magnitude of this new kind of Hall effect in ABPs by comparing the amount of directed motion with the undirected motion of the active particles. We present a possible experimental realization of our findings: The motion of a large tracer particle, insensitive to Lorentz force, would become anisotropic due to the presence of y/l y/l y/l x/l x/l x/l FIG. 1. Nonequilibrium steady state of an internally driven system (active Brownian particles) of noninteracting particles subjected to a Lorentz force arising from an external magnetic field. l denotes the persistence length of the ABP, q its charge, γ the friction coefficient, v0 its self-propulsion speed, and ρ b is the bulk density. The Gaussian shaped magnetic field is radially symmetric and points in the z direction B(x, y)êz. The steady-state is characterized by (i) inhomogeneously distributed density and (ii) bulk fluxes. The direction of the flux, indicated by the arrows, is perpendicular to the gradient of the magnetic field. The flux can be reversed by reversing the magnetic field. In the absence of internal driving, there are neither fluxes nor density gradients in the steady state for any magnetic field. This also holds for spatially homogeneous driving and no magnetic field. However, when one considers the two together, homogeneously driven particles in an inhomogeneous magnetic field, the system exhibits the unusual nonequilibrium behavior. The circulating fluxes are reminiscent of the Corbino effect in conductors [16].
The data was obtained from Brownian dynamics simulations of active Brownian particles subjected to Lorentz force (App. F).
the (active) flux patterns. Using reasonable estimates for the system parameters, our calculations suggest that this anisotropy could be measured experimentally.

II. ACTIVE PARTICLES WITH LORENTZ FORCES
We consider a single spherical, self-propelled, charged Brownian particle in a magnetic field B(r) with a spatially varying magnitude but constant direction. The active motion of the particle is modeled as a force in the direction of the orientation of the particle specified by an embedded unit vector p undergoing rotational diffusion. The state of the particle is determined by the position coordinate r(t), the velocity v(t) and the orientation p(t).
The dynamics are described by the following (Stratonovich) stochastic differential equations: and where m is the particle mass, γ is the friction coefficient, f is the self-propulsion force, T is the temperature in units such that Boltzmann constant is equal to one and D r is the rotational-diffusion constant. The (dimensionless) matrix Γ(r) = 1 + qB(r) γ M , where q is the charge of the particle and M is a matrix such that B(r)M · v = B(r) × v, and B(r) ≡ |B(r)|. The stochastic vectors ξ and η are Gaussian distributed with zero mean and autocorrelation ξ(t)ξ T (t ) = η(t)η T (t ) = δ(t − t )1. Note that there is no direct coupling between the orientation and the magnetic field.
We start with Eq. (3) with the inertia term because in presence of a Lorentz force, the overdamped equation of motion cannot be used to determine observables depending on the velocity process (such as flux), even though the position-statistics can be accurately described [4,5,44]. However, the Fokker-Planck equation (FPE) for the overdamped motion corresponding to Eqs. (2) and (3) can be obtained by a different method based on an expansion in powers of mass. This expansion removes the velocity degrees of freedom (see App. A) and yields the following FPE for the probability density Q(t) ≡ Q(r, p, t): is the translational flux, is the rotational flux, R = p × ∇ p is the rotation operator, and with This parameter, which we refer to as the diffusive Hall constant, measures the strength of the Lorentz force relative to the frictional force in the system. When this parameter is comparable to unity, the Lorentz force affects the motion of the particle significantly. As we show later in Sec. V, the Lorentz force may become comparable to the frictional force for highly charged particles in strong magnetic fields. Note that the matrix Γ −1 in Eq. (8) is not symmetric, and therefore, T Γ −1 /γ cannot be interpreted as a diffusion tensor [5]. Due to the antisymmetric terms in Γ −1 , there are additional fluxes in the system that are perpendicular to the density gradient.
Equation (5) can be coarse grained further by a gradient expansion to remove the orientational degrees of freedom. We consider the ansatz where the scalar ρ, the vector σ and the matrix τ are function of r and t, and Ω contains all higher-order terms. Equations for the time evolution of the probability density ρ, the total orientation σ/3 and higher-order tensors can be obtained by inserting Eq. (10) into Eq. (5) and integrating out the orientation degrees of freedom (see App. B). In the limit where gradients are small compared to the persistence length of the ABP [39,49], one obtains a FPE for the probability density ρ(r, t) as a function of time and position alone: where where is the persistence length of the ABP, andp(r, t) is the polarization, which is defined as the average orientation per particle pδ (r − r) /ρ(r, t) where the angle brackets indicate an average with respect to the probability density of r , p and v at time t.

III. NONEQUILIBRIUM STEADY STATE
Here we focus on the steady-state bulk behavior of the system and therefore use periodic boundary conditions in all directions. Furthermore, the analysis is restricted to a magnetic field with a constant orientation, in this case the z direction, and a magnitude that depends on a single coordinate, say x. Since the motion along the direction of the magnetic field is not affected by the Lorentz force, there is neither flux nor density variation along the z direction in steady state. For the chosen magnetic field, the system is translationally invariant in the y direction; therefore, the flux and the polarization in the y direction are independent of the y coordinate. This also implies that in steady state there can be no flux in the x direction because of the periodicity and equivalence of the positive and negative x directions. Moreover, there can be no local x fluxes forming closed loops via the y direction because of the translation invariance in the y direction. The absence of fluxes in the x direction has been corroborated by simulations.
By setting the x component of the flux [Eq. (12)] to zero, we obtain the following analytical expression for density: where where ∆T = f 2 /6D r γ. The polarization can be calculated by inserting Eq. (15) in Eq. (13), which yields for the x component, and for the y component, where κ (x) = dκ(x) dx . The flux has two contributions, one coming from the polarization and the other from the gradient in the density [see Eq. (12)]. In steady state, the two contributions to the flux in the x direction are equal in magnitude, but in opposite direction, resulting in zero flux in this direction. In the y direction there can be no density gradients because of the translation symmetry in this direction. This means that the contribution to the flux in the y direction coming from a density gradient along this direction is zero; however, the contribution from the polarization is not zero, resulting in a flux in the y direction. This flux can be calculated by using the solution for density [Eq. (15)] as input to the y component of Eq. (12), yielding Note that the orientation and flux in the y direction are perpendicular to the magnetic field gradient. Before we test the predictions of the theory, a few remarks are in order. First, ∆T appearing in Eqs. (15)- (18) has units of energy (note that the Boltzmann constant is unity, so T has units of energy as well) and has previously been interpreted as an effective increase in the temperature due to activity [50,51]. However, we do not interpret ∆T as increase in temperature because this system cannot be mapped to an equilibrium system at a higher temperature. Second, the expressions for the density, flux and polarization can be improved significantly by combining the gradient expansion with linear-response calculations of the polarization. Such a calculation results in a differential equation for the polarization, which can be solved numerically. This solution can be used to replace the linear-order contribution of the activity to the polarization calculated using the gradient expansion and using this polarization as input to the flux in Eq. (15). Details of the linear-response theory, its application to the polarization and the combined solution are shown, respectively, in Appendices C, D and E. Finally, it is interesting to note that j y (x)/ρ(x) is the macroscopic velocity of the ABP in the presence of the gradient field, while f /γ gives the intrinsic microscopic velocity of the ABP. While the intrinsic velocity leads to random motion on large time scales, the field-induced motion persists (in our case in the y direction) on any time scale. Thus the polarization [Eq. (18)] gives the fraction of the swim force that is converted into macroscopic directed motion. Figure 2 shows the predictions of the theory for the density, flux and polarization for κ(x) = 4 sin(2πx/10) and γ = 1 together with the results from Brownian dynamics (BD) simulations. Consistent with the prediction of equation (15), the density is high where the magnitude of the magnetic field is large (see Fig. 2 a). Qualitatively, the accumulation of particles in the high magnetic field region can be understood as follows. A magnetic field hinders the diffusive motion of Brownian particles because there is a force perpendicular to the velocity which results in circular motion. For passive Brownian particles this means that on time scales longer than the velocity autocorrelation time the effect of a uniform magnetic field is a decrease of the diffusion constant in the plane perpendicular to the magnetic field [3]. For a spatially varying magnetic field, the diffusion coefficient of particles is small in the region where magnetic field is large. Whereas, for passive particles, the steady-state density distribution is independent of the diffusion coefficient, active particles accumulate in the region of low swim speed, which correspond to regions of high friction, or in this case high magnetic field [7,37].
Passive Brownian particles in a space dependent magnetic field have a uniform steady-state density (a Boltzmann distribution with a uniform potential). A system of ABPs with a homogeneous activity also has a uniform steady-state density. However, the combination of a space-dependent magnetic field and uniform activity results in an inhomogeneous steady-state density. An inhomogeneous density without a space dependent potential means that the density does not follow a Boltzmann distribution, a hallmark of nonequilibrium systems.
There is polarization in the plane perpendicular to the magnetic field; see Eqs. (17) (18) and Figs. 2 c and d. In steady state, the polarization and the gradient of the density in the x direction are such that the fluxes resulting from them cancel each other. The polarization in the y direction, however, cannot be canceled by a density gradient because when the magnetic field varies only along the x direction, the system is translation invariant in the y direction and therefore, the density, flux and polarization cannot depend on the y coordinate. Because there is polarization in the y direction, but no density gradient to counter the flux from the polarization, there is in steady state a flux in the y direction [see Eq. (19) and Fig. 2 b]. Due to the periodicity of the field, there is flux in both the positive as the negative y direction resulting in zero center-of-mass movement.
Note that there is no direct coupling between the magnetic field and the polarization [see Eqs. (3) and (4)], and the inhomogeneous density along the x direction is not a result of the polarization. Both the polarization and the inhomogeneous density are kinetic in nature, and result from a filtering mechanism. The magnetic field hinders the motion of the ABP and a particle with an orientation parallel to the gradient of the magnetic field will get stuck. A particle with an orientation antiparallel to the gradient of the magnetic field can more easily move towards regions of lower magnetic field. These effects result in a force imbalance between regions of small and large magnetic fields. In steady state the force imbalance is canceled by gradient in the chemical potential due to the increase of density in the high magnetic region.
The presence of steady-state bulk fluxes is a key feature of the nonequilibrium steady state induced by Lorentz forces in a system of active particles. These fluxes are unique in the sense that the underlying mechanism is independent of inter-particle interactions [22], aligning torques [23] or sliding motion along asymmetric walls [17][18][19][20][21]. Moreover, the direction of fluxes is perpendicular to the gradient of the applied magnetic field. This particular property is what gives rise to the circulating steady-state fluxes shown in Fig. 1 in which the magnetic field is radially symmetric. The direction of the fluxes can be reversed by simply reversing the direction of the magnetic field.
In case of a finite system with walls, the bulk properties of the system are similar to that of an infinite system; however, there is also flux parallel to the wall resulting in closed-loop fluxes (see Fig. 3). A flux along walls has recently been investigated in the case of chiral active particle [53]. The flux along the wall induced by Lorentz force can possibly be explained by the same mechanism.  19)]. There is flux in negative y direction around x = 0 and x = 40l, and in the positive direction around x = 20l. Note that, due to the periodicity of the system, the integral of the flux over the x coordinate is zero and there is no center-of-mass motion. From (c) and (d) it is clear that the polarization is nonzero in a narrow space interval and changes rapidly in this interval; for example, from x ≈ 18l to x ≈ 22l the polarization rotates clockwise from the negative to the positive x direction. The analytical predictions from the gradient expansion overestimate the change in density, the flux and the orientation; the prediction obtained from the combination of the gradient expansion and linear response theory, however, has good agreement with the simulation results.

IV. ACTIVE DIFFUSION UNDER LORENTZ FORCE
Here we investigate the mean squared displacement (MSD) of an ABP in the three directions. The MSD in the y direction is particularly interesting: for a magnetic field varying in the x direction, the steady-state fluxes in the y direction constitute layers of alternating flows (see inset to Fig. 4 and Fig. 2). How does the emergence of these layers of alternating flows influence the MSD in the y direction?
Free ABPs have a swim speed of f /γ [54]; therefore, on short time scales the motion of ABPs is ballistic and is the same in the three directions: where ∆x(t) is x(t) − x(0), and similarly for ∆y(t) and ∆z(t). The diffusion parallel to the magnetic field orientation is unaffected by the field. This implies that the diffusion coefficient in the z direction is that of a freely diffusing ABP, which is Since a magnetic field hinders the motion in the plane perpendicular to its direction, the diffusion in the x and y directions is reduced. For an inhomogeneous magnetic field, the spatially varying diffusion coefficient in the x and y directions can be obtained from the symmetric part of the tensor in Eq. (8) as On a length scale larger than the period of the magnetic field, the space dependence of the diffusion constant be- comes irrelevant, and one can approximate the effective diffusion constant for MSD in the x direction by averaging over space: where L is the period of the magnetic field. There are three time scales in the system. Under the assumption that the length scale of ballistic motion is much smaller than the period of the field, the shortest time scale, τ 1 , is the transition from ballistic to diffusive motion in the x direction. This time scale can be obtained from (f τ 1 /γ) 2 = 2D (ef f ) x τ 1 , and is The second time scale, τ 2 , is the transition from ballistic to diffusive motion in the z direction, which can be calculated from (f τ 2 /γ) 2 = 2D z τ 2 as The largest time scale, τ 3 , is the transition to diffusive motion in the y direction. The diffusion in the y direction is increased due to the flux lanes in the system. Because the net flux in the system is zero, the motion in the y direction becomes diffusive on the same time scale as the time it takes to diffuse from on lane to another. This time scale can be calculated from 2D The MSD in the y direction becomes subballistic on the same time scale as the x direction, namely τ 1 ; however the diffusive regime starts later at τ 3 . This means that there is a super diffusive regime that persists over several decades in time. Figure 4 shows that when t > τ 3 the MSD in the y direction has the same effective diffusion coefficient as the z direction. The magnetic field hinders diffusion in both x and y direction. Whereas the diffusion in x direction is expectedly smaller than in the z direction, the flux lanes increase the effective diffusion along y and therefore the MSD in the y direction 'catches up' with z direction.  21)]. The motion of an ABP along the x is also diffusive at long times ∆x(t) 2 ∼ 2D ef f x t (dotted line) albeit with a reduced diffusion coefficient as given in Eq. (23). This is expected since a magnetic field hinders the motion of particle in the plane perpendicular to its orientation. The motion along the y direction, however, is affected differently than along x, which is related to the flux lanes (see inset). At short times, t < τ2, the motion along y is same as x. The time scale of the transition from ballistic to diffusive motion in the z direction is τ2 [Eq. (25)]. However, at long times t > τ3, the rate of diffusion along y becomes the same as that of a freely diffusing ABP. τ3 is the time scale of the transition to diffusive motion in the y direction [Eq. (26)]. At intermediate times τ2 < t < τ3, the motion along y is super diffusive.

V. PASSIVE DIFFUSION OF TRACER PARTICLE
We showed in Sec III that the polarization [Eq. (18)] is the fraction of the swim force that is converted into macroscopic motion. Since polarization scales as lκ 1, it is small in the limit where gradients are small compared to the persistence length of the ABP. To overcome this limitation, we consider a system in which the effect of the anisotropic motion of the ABPs on a tracer particle can accumulate. Consider again the system of ABPs in steady state in a periodically varying magnetic field. As discussed above the ABPs will form a periodic array of lanes carrying an alternating macroscopic flux in y direction, see Fig. 5. We now immerse an uncharged, passive tracer particle in the system and study its MSD. On can expect that the steady-state fluxes of the ABPs affect the diffusion in the y direction of the tracer particle. Similar problems have been studied in different contexts, such as the diffusion of a passive particle in presence of random velocity fields [55][56][57] and the motion of polymer chains in random layered flows [58].
Due to the (alternating) velocity field, one can expect that the MSD of the tracer particle along y direction is enhanced. This implies that the motion of the particle is anisotropic. The enhancement of the diffusion along the y direction can be estimated as follows: The particle diffuses freely with a diffusion rate D p between the velocity lanes. Assuming that inside the lane, the particle moves with the local velocity, the distance that the particle travels along the lane is l f = v y d 2 /2D p , where v y is the local velocity and d is the width of the lane. The particle experiences these flow lanes every τ y = L 2 /8D p seconds, where L/2 is the distance between two neighboring lanes. The increase in the diffusion rate along the y direction can be estimated as l 2 y /2τ y ∼ v 2 y d 4 /D p L 2 . The anisotropy parameter, defined as the ratio of the MSDs in the x and y directions, is An experimental measurement of the anisotropy of the diffusion of a tracer particle, as described above, would provide an, albeit indirect, experimental verification of the active fluxes in the system. Since Lorentz force is relativistic in nature and ABPs are colloidal particles, one expects that v y is extremely small. As we show below, this is indeed the case for the known experimentally realizable ABPs. However, since the anisotropy parameter depends inversely on the diffusion coefficient of the tracer particle, it might still be observed experimentally for sufficiently large tracer particles on long time scales. The smaller the diffusion constant of the tracer particle the longer it will stay in one lane and thus accumulating the effect by drifting along with the velocity field. In the extreme case of a vanishing diffusion constant, D p =0, the only observable displacement on a long time scale is the drift and thus r a 1. The anisotropy parameter can be estimated as follows: The lane velocity can be obtained from Eq. (19) as v y = j y /ρ = v 0py , where v 0 = f /γ is the self-propulsion speed andp y is the y component of the polarization vector inside a lane. It can be obtained from Eq. (18) as p y ≈ lq|∇B|/3γ, where l = v 0 /2D r is the persistence length and |∇B| is the (maximum) gradient of the magnetic field. With this the lane velocity can be estimated as v y ∼ v 2 0 q|∇B|/(2D r γ). To estimate the charge, we consider that the active particle is a Janus-shaped particle with half of its surface grafted with polymer chains, and the other half serves as the propulsion engine [8][9][10][11][12]. The total charge on such an ABP is q = 2πR 2 σ b λe, where R is the radius of the ABP, σ b is the grafting density, λ is the number of charges per polymer chain, and e is the elementary charge. Using this together with the Stokes relations for the rotational diffusion coefficient D r and γ, one obtains In order to estimate the magnitude of lane velocity, we consider T = 298/k b kgm 2 s −2 , where k b is the Boltzmann constant, a dense grafting density of σ b = 1nm −2 and highly charged chains with λ = 100. For the magnetic field we assume a period of L ∼ 1 cm. For the magnetic field gradient we use 10 T /cm, which can be obtained in magnetic matrices used for filtering fine and weakly magnetic particles [59]. For the ABPs we take the values v 0 ≈ 3 µms −1 and R ≈ 2 µm from Ref. [23]. With these, we obtain the lane velocity v y ≈ 2.5 × 10 −3 µms −1 . The width of a lane can be approximated as d = L/(2π) ∼ 0.1 cm. Note that for the chosen parameters the diffusive Hall-constant is of the order, κ = qB/γ ∼ 0.1 (29) implying that Lorentz force becomes comparable to frictional force on the particle. Finally, assuming that the diffusion coefficient of the tracer particle is D p ∼ 10 −13 m 2 /s, a reasonable estimate for large tracer particles (∼ µm), we obtain from Eq. (27) r a ≈ 2. Of course long-time scales are necessary to detect the effect. However, the distance that a tracer particle drifts along a velocity lane in a day time turns out to be approximately 0.2 mm, which can be detected experimentally. These estimates suggest that Lorentz forces may be observable experimentally for charged, active particles. We note that it is not the magnitude of the magnetic field that is important for this Hall-like effect, but it is only the gradient which enters the calculation. A few remarks are necessary with respect to a possible experimental setup. First, the total charge of the polyelectrolyte brush is compensated by counter-ions arising from the dissociation of the charged groups of polymers in water. In a salt-free solution the counter-ions will be strongly localized inside the brush and thus effectively cancel the effect of the Lorentz-force [60]. However, a superior aspect of the external magnetic field is that it is not screened by electric charges. Thus, at a high salt level in the solution the counter-ions are quasi-free [61,62]. Furthermore, the exchange dynamics of them will be much faster than the Lorentz-induced drift. Therefore, brushdecorated active colloidal particles in a salted solution could be a promising experimental system.

VI. DISCUSSION AND CONCLUSIONS
Lorentz force has the unique property that it depends on the velocity of the particle and is always perpendicular to it. Although this force generates particle currents, Inside a lane, the particle drifts along with the local velocity of the lane shown as white arrows. This results in anisotropic diffusion of the particle. Since the time spent by a particle inside a lane increases with decreasing diffusion coefficient, the anisotropy can be large for small diffusion coefficients.
they are purely rotational and do no work on the system. As a consequence, the equilibrium properties of a Brownian system, for instance the steady-state density distribution, are independent of the applied magnetic field. The dynamics, however, are affected by Lorentz force: the Smoluchowski equation picks up a tensorial (diffusion) coefficient, which reflects the anisotropy of the particle's motion [3]. The diffusion rate perpendicular to the direction of the magnetic field decreases with increasing field whereas the rate along the field remains unaffected. In addition to this effect, Lorentz force gives rise to nondiffusive fluxes which are perpendicular to the density gradients [4,5].
The effects caused by the Lorentz force, however, occur only in nonequilibrium and cease to exist when the distribution of charged particles reaches equilibrium. A system of ABPs, in contrast, is continuously driven out of equilibrium. Lorentz force couples to the nonequilibrium dynamics of an ABP via its self-propulsion. In this work, we showed that when subjected to an inhomogeneous magnetic field, a system of uniformly active ABPs settles into an unusual nonequilibrium steady state. This state is characterized by (a) inhomogeneously distributed density with the same spatial symmetry as that of the magnetic field and (b) macroscopic fluxes perpendicular to the gradient of the magnetic field.
The theoretical approach is based on coarse graining of the Smoluchowski equation and yields analytical results for the density distribution, active fluxes, and polarization in the steady state. We also showed that when combined with the linear-response theory, the theoretical predictions match the results from BD simulations.
That the density distribution in our system is inhomogeneous can be understood as follows. The swim force increases the ABP's velocity in the direction of its orientation. The Lorentz force rotates the velocity vector. The effect of this rotation is a decrease in the speed of the particle in the direction of its orientation, which is similar as decreasing the effect of the swim force. From this it follows that a spatially varying magnetic field has a similar effect as ABP as subjecting them to an inhomgeneous activity field. However, the existence of nondiffusive fluxes perpendicular to the density gradient in the steady-state is a unique feature of this system and is a consequence of the antisymmetric coupling of the x and y directions due to the Lorentz force. In case of a finite system with walls, the bulk properties of the system are similar to that of an infinite system; however, there is also flux parallel to the wall resulting in closed-loop fluxes The MSD of an ABP in the presence of an inhomogeneous magnetic field has nontrivial behavior related to the existence of flux lanes in the system. Whereas the motion along the direction of the field is unaffected by the field, it is hindered in the plane perpendicular to the field. The diffusion coefficient in the z direction is, as expected, the same as that of a freely diffusing ABP. The diffusion in the x direction is expectedly smaller than in the z direction; however, the same does not hold for the y direction. The diffusion in the y direction is super diffusive at intermediate times. At long times, the motion along the y direction becomes diffusive with the same rate as along the z coordinate. The emergence of flux lanes in the y direction gives rise to this nontrivial behavior of the MSD in the y direction.
Our study naturally raises the question whether the effects induced by Lorentz force in a system of ABPs can be observed experimentally. A Lorentz force on a colloidal particle is extremely small due to the particle's low velocity. However, for such a particle to be affected by a Lorentz force, it is not the absolute value of the Lorentz force that is important but how it compares with the frictional forces. Since both frictional and Lorentz force depend on the velocity, for Lorentz force to have any significant effect on the motion of the particle the diffusive Hall constant κ = qB/γ must be of the order of unity. In Sec. V, we considered the possibility of experimental verification of this study where we have proposed to use the high total charge of a poly-electrolyte brush grafted to an ABP to obtain a value of κ 0.1. Direct experimental verification of the flux lanes might still be difficult. A possibility is to measure the effect of the flux lanes on the diffusive behavior of a passive tracer particle on which the effect of the flux lanes can accumulate. Using values for propulsion speed and particle size from experiments [23] and large but reasonable estimates for other parameters in the system, we estimated a drift of the tracer particle along the flux lines of the active system of the order of millimetre per day. This indicates that magnetic-field induced effect on a charged active system could be measured experimentally.
Finally, note some general aspects of the application of magnetic fields to soft-matter systems. First of all, magnetic fields are not screened by static charges such as electric field are. Strong electric field strength in particular combined with thin-film geometries are prone to short-cuts, irreversibly damaging the system. By contrast the application even of high magnetic fields is possible and is used in, for example, high-resolution NMR-experiments to medicine. In this work we have shown that Lorentz forces may be detectable in active soft-matter systems, which opens new possibilities to manipulate such systems. Interestingly, macroscopic effects are induced by the symmetry breaking of the magnetic field only and do not require a macroscopic flux of the particles which establishes a novel quality of the Hall-like effect, which in this case depends on the gradient of the field only.
Appendix A: small-mass limit In order to obtain the small-mass limit of the FPE corresponding to Eqs. (2)(3) and (4), one can use the method described in Chap. 10 of Ref. [63]. This was done for a single passive particle in a homogeneous magnetic field in Ref. [4]. We start with the FPE corresponding to Eqs. (2), (3) and (4) for the time evolution of the probability density where the time-evolution operator has been split up in a reversible part an irreversible part without the part coming from Eq. (4) the part representing the rotation of the orientation vector with the rotation operator R = p × ∇ p [65], and the active, nonequilibrium operator coming from the selfpropulsion where L can be any of the operators in Eq. (A1), and is the solution to L irr R(v) = 0, normalized such that the integral over v is one. The transformed operators arē The eigenfunctions of the operator b † α b α , where α is x,y or z, are and where H n are Hermite polynomials. The operators b † α and b α are the raising and lowering operators of the eigen- . Because the eigenfunctions are orthonormal, they can be used to expandP (t): (A17) where c nx,ny,nz ≡ c nx,ny,nz (r, p, t).
The probability density for the position and orientation, Q(t) ≡ Q(r, p, t), is given by the first expansion coefficient: can be evaluated by using Eqs. (A7), (A17) and J (r, p, t) = T /mc 1 (r, t). (A22) Using this, together with the definition of the flux in orientation space one obtains the following equation for the time evolution of the probability density: Next, we need an expression for the c 1 in terms of Q(t). We are interested in the case where the magnetic field is oriented perpendicular to the gradient of the field, and take this to be the z direction, so B(r) = B(r)ẑ. Equation (A6) together with the orthonormality of the eigenfunctions yields an hierarchy of equations for the functions c nx,ny,nz called a Brinkman hierarchy [66]: ∂ ∂t c nx,ny,nz = − γ m c nx,ny,nz (n x + n y + n z ) − D ·   √ n x + 1c nx+1,ny,nz n y + 1c nx,ny+1,nz √ n z + 1c nx,ny,nz+1 √ n x c nx−1,ny,nz √ n y c nx,ny−1,nz √ n z c nx,ny,nz−1   + qB(r) m n x (n y + 1)c nx−1,ny+1,nz − qB(r) m (n x + 1)n y c nx+1,ny−1,nz + D r R 2 c nx,ny,nz , where D = T m ∇ r andD = T m ∇ r − 1 √ T m f p. The equation governing the time evolution of the first expansion coefficient is where c 1 = (c 1,0,0 , c 0,1,0 , c 0,0,1 ) T . This is of course the same as Eq. (A24). The equation for c 1 (t) is again a differential equation in t, which can be solved after a Laplace transformation: where the tilde indicates a Laplace transformation from the variable t to s, D 1,2 is a matrix of which the entries are linearly proportional to D x , D y or D z , c 2 (t) is a vector of which the elements are all c nx,ny,nz such that n x + n y + n z = 2, and Note that Γ m is a matrix of which the entries are at least of order m 0 , D 2,3 is a matrix whose elements are proportional to D x , D y or D z ,D 2,1 is a matrix whose elements are proportional toD x ,D y orD z , and c 3 (t) is a vector of which the elements are all c nx,ny,nz such that n x + n y + n z = 3.
By only retaining the zeroth-order contribution in m to Eq. (A26), one obtains an equation that is independent of v and holds for sm/γ << 1 (that is, for times longer than the velocity autocorrelation time). The vector D is of order m −1/2 , so only the terms of order m 1/2 of c 1 are needed, which means that one does not need the matrices D 1,2 , D 2,3 andD 2,1 . Equation (A29) can be used to replace the c 2 in Eq. (A27), which, after inverting the Laplace transformation and only retaining terms of order m 1/2 or lower, becomes which gives for the translational flux where The result for the translational flux together with Eqs. (A20) and (A26) gives following Fokker-Planck equation: where the ∂ c Υ cab comes from the projection on to the third and higher-order harmonics. Equations (B2) (B4) (B5) are exact, but not closed since Υ is not specified. We now consider the limit in which the gradients in the system are small in comparison to the persistence length of the ABP [39,49]. The relaxation time of ρ is of order ∼ (∇) −1 , whereas all the other harmonics relax in times of order ∼ 1. This implies that in the limit of small gradients, ρ is the slowest mode. We thus assume that the time derivative of σ and terms with both a σ and a ∇ can be neglected. This yields which is related to the polarization byp = 1 3 σ/ρ. With this expansion the flux becomes In a system in steady state with periodic boundary conditions and a magnetic field of the form B(x)ê z , the flux in the x direction must be zero (see Sec. III). This condition together with Eq. (B7) yields the steady-state density: where the proportionality constant is determined by the normalization of the density, δ = ∆T T +∆T and ∆T = f 2 6Drγ . Using this density to calculate the y component of the steady-state flux gives The steady-state polarization in the x and y directions can be determined from Eq. (B6) and the density: Appendix C: Linear-response theory Linear-response theory for ABPs has been developed in Ref. [38] to study the average swim speed, and was later applied to a system of interacting ABPs with inhomogeneous activity to study the density profile [37] and flux [27]. Here we adapt the linear-response theory for ABPs and apply it to an charged ABP in a space-dependent magnetic field. where is the operator corresponding to the particle's diffusive behavior, L r = L r is the rotational time-evolution operator and is the nonequilibrium time-evolution operator. Next, we split the probability density Q(t) in an equilibrium (Q eq ) and a nonequilibrium (δQ(t)) part, where the equilibrium part is the steady-state solution to Eq. (C1) with f = 0, which is constant in space. The equation for the nonequilibrium part is because ∂ t Q eq = (L d + L r ) Q eq = 0. We are interested in nonequilibrium steady state averages, and therefore assume that the system started out at t = −∞ in the equilibrium state. The previous equation can be integrated from t = −∞ to t = 0 by treating L a Q eq as an inhomogeneity: (C5) where δQ ss = δQ(t = ∞) is the steady-state solution of the nonequilibrium part of the probability density. This equation is exact and can be used to calculate the steadystate nonequilibrium average of a generic function g = g(r, p): g r,p = dr dp Q ss (r, p)g where Q ss (r, p) = Q eq +δQ ss (r, p) is the nonequilibrium steady-state probability density, the super script r, p denotes the variables that are averaged over, the subscript eq indicates that the average is taken with respect to the equilibrium distribution, and a † indicates the adjoint of an operator. The adjoint operators acting on a function g are and where the superscript T indicates a transpose. The response of the system to the activity up to linearorder in f is obtained by removing from the full timeevolution operator the nonequilibrium operator L † a : g r,p This Green-Kubo relation defines the active transport coefficient α corresponding to an functions g by α = lim f →0 (< g > lr − < g > eq ) /f . Equation (C11) shows that only quantities that are odd functions of p have a nonzero linear response because the angular average in equilibrium any function of an odd power of p yields zero by symmetry.
Without loss of generality one can consider the case where the magnetic field is oriented along the z direction. In this case the matrices M and M 2 are so clearly there is no orientation along the direction of the magnetic field. The operator L † d acting on r in Eq. (D6) has the same effect as acting with the adjoint on r (see, for example, Chap. 4.2 of Ref. [63]); thereforē and was used. In the linear-response calculation of the orientation ρ(r ) can be replaced with ρ b . Equation (D12) can be simplified tō which is equivalent to the following differential equation: This equation is the main result of this appendix. Solutions can be obtained numerically, or by further approximations. For large D r and a magnetic field with small gradients, the contribution from the operator L d to the left hand side of Eq. (D18) is negligible compared to the contribution coming from the part with D r ; therefore, one can approximate, in this case, the solution to this equation byp (1) Note that this is the same Eqs. (B10) and (B11) without the terms of cubic order in f . Even though the density does not have a linear response to the swim speed, the linear response calculations for the orientation can be used to calculate the density and fluxes [37]. First, one integrates out the orientational degrees of freedom from Eq. (A24): ρ(r, t) = −∇ r · j(r, t), with j(r, t) = − 1 γ Γ −1 · [T ∇ r − fp(r, t)] ρ(r, t), wherep(r, t) ≡ 1 ρ(r,t) dpQ(r, p, t) is the polarization. When the magnetic field is oriented in the z direction and depends only on the x coordinate, equating the steady-state flux in the x direction to zero gives the following density: where the proportionality constant is determined by normalization, and the superscript (2) indicates that it is correct up to second order in the swim force (there is no linear-order response in the density). With this density the steady-state flux in the y direction can be calculated from Eq. (D21): where ρ b is the bulk density.

Appendix E: Improved solution
Only using the linear order contributions in f to the orientation gives poor agreement with simulations for the density; see Fig. 6. The terms of cubic order in f contribute to the diffusion in the system, and by ignoring them there has to be a larger density gradient to compensate for the flux resulting from the orientation. This problem can be overcome by combining the first order result from linear-response theory with the cubic terms from the gradient expansion.
Using the gradient expansion (App. B) we found that the flux is wherep =p(r, t) is the average polarization at position r, which is The second term on the right hand side is higher order in the swim force because a gradient in the density is at least second order. The first term on the right hand side of this equation is proportional to the swim force and is equal to the linear-response solution in the limit of small gradients. However, Eq. (D18) gives the linear-order response of the polarization without assuming that the gradients are small. We therefore replace the linear order term with the numerical solution to the linear-response equation for the polarization [Eq. (D18)] and rewrite the previous equation as When the magnetic field is oriented in the z direction and varies in the x direction, one can use this expression for the polarization together with the expression for the flux [Eq. (E1)] to obtain the steady-state density by equating the flux in the x direction to zero: where N is a normalization constant such that And finally, the density can be used in Eq. (E3) to obtain the polarization in the x and y directions in terms of the linear-response polarization as andp y (x) =p (1) y − δκ 1 + κ 2 p (1) x + κp (1) y . (E7)