Flocking from a quantum analogy: Spin-orbit coupling in an active fluid

Systems composed of strongly interacting self-propelled particles can form a spontaneously flowing polar active fluid. The study of the connection between the microscopic dynamics of a single such particle and the macroscopic dynamics of the fluid can yield insights into experimentally realizable active flows, but this connection is well understood in only a few select cases. We introduce a model of self-propelled particles based on an analogy with the motion of electrons that have strong spin-orbit coupling. We find that, within our model, self-propelled particles are subject to an analog of the Heisenberg uncertainty principle that relates translational and rotational noise. Furthermore, by coarse-graining this microscopic model, we establish expressions for the coefficients of the Toner-Tu equations---the hydrodynamic equations that describe an active fluid composed of these"active spins."The connection between self-propelled particles and quantum spins may help realize exotic phases of matter using active fluids via analogies with systems composed of strongly correlated electrons.

Active liquids exhibit striking phenomena due to the unusual nature of their hydrodynamics [1]. Such phenomena have been observed in naturally occurring collections of live animals [2][3][4] and cells [5][6][7][8][9], as well as synthetically prepared systems of granules [10,11], robots [12], colloids [13][14][15], and molecules [16][17][18][19]. Coarse-grained descriptions that capture these phenomena may be either constructed based solely on symmetry and lengthscale considerations or derived from simple particle-based models [15,[20][21][22]. A crucial advantage of the latter, microscopic, approach is that it connects the hydrodynamic coefficients (such as viscosity, diffusivity, and compressibility) to the microscopic parameters of the model. In experimental realizations of active fluids, this connection between microscopics and hydrodynamics can be used to construct design principles targeting the realization of novel materials and devices. For example, recent work has focused on the robustness of active liquids against disorder [23], the design of flow patterns in confined active fluids [24][25][26][27][28], and the use of such channel networks for the design of topological metamaterials [29] and logic gates [30].
One specific challenge is that the coarse graining of a microscopic model of self-propelled particles is, in general, technically difficult. As a result, specific counterintuitive phenomena associated with active-liquid hydrodynamics are difficult to describe in generic, modelindependent terms. In pursuit of this goal, the introduction of additional minimal models of self-propelled particles, along with their coarse-grained hydrodynamics, can serve to strengthen the connection between small-and large-scale phenomena in active systems.
In this work, we introduce a particularly simple, minimal model of self-propelled particles, and we explore their individual and collective statistical dynamics in order to arrive at a hydrodynamic description. The simplicity of the model stems from its connection to the Schrödinger equation describing a quantum particle. As such, we use basic results from quantum mechanics to develop physical intuition for active-fluid phenomena. For example, we describe active-fluid analogs of such well-known quantummechanical concepts as spin, spin-orbit coupling, and the Heisenberg uncertainty principle. We discuss how the analog of a spinor can be used to introduce a propulsion direction via spin-orbit coupling. We then construct a probabilistic, Fokker-Planck interpretation for the dynamics of a single self-propelled particle in the presence of translational noise; see Fig. 1. We show that the microscopic model we consider includes feedback between rotational and translational noise, which we interpret as an analog of the Heisenberg uncertainty relation. Crucially, we use this single-particle model to construct a hydrodynamic description of a system of many self-propelled particles. We thus obtain simple relations between the coefficients in the Toner-Tu model [31] and the microscopic parameters of the individual particles under consideration, including their interactions. We are then able to conclude that, as for any model in the Toner-Tu universality class, the many-particle system we consider exhibits long-range orientational order in two dimensions as a consequence of activity [31].

I. MODEL OF SELF-PROPELLED PARTICLES
We begin with the well-known connection between non-relativistic quantum mechanics (described via the Schrödinger equation) and classical statistical mechanics (described via the diffusion equation). Consider the Schrödinger equation: For the free-particle Hamiltonian operatorĤ =p 2 /2m = − 2 ∇ 2 /2m (in the position representation), a rotation of time into the imaginary axis via t → −it transforms this Schrödinger equation into the diffusion equation: x Translational Noise Self-propulsion FIG. 1. a) Schematic illustration of single-particle dynamics in the model we consider. The particle trajectory is composed of displacements due to self-propulsion combined with translational noise. Such noise may arise, e.g., from a fluid in which the particles are suspended. In addition to the translational noise, which alters the particle displacement, the particle is subject to orientational noise, which alters the direction of self-propulsion. In the model we consider, the two sources of noise are intimately coupled, which leads to a relation analogous to the Heisenberg uncertainty principle. b) For a system of many particles, we consider processes that can align the self-propulsion directions of two colliding particles. For example, elongated active particles may prefer to align due to the dynamics of their collisions. In the coarse-grained model, we capture the strength of the alignment interaction via the parameter g.
In the diffusion equation, Ψ can be identified with the particle density ρ and 2m with the diffusion constant D. This bridge allows us to use tools from quantum mechanics to characterize classical stochastic phenomena. However, this approach does not capture self-propulsion or spontaneous active flow, which cannot be described via the diffusion equation.
In order to capture self-propulsion, each particle ought to carry information about its direction of motion, for which we need to introduce additional degrees of freedom. On the quantum side of the analogy, these degrees of freedom capture the quantum spin state. For a system having spin, Ψ is an n-component spinor (we consider n = 2), and has additional symmetries with respect to spin rotation, which we consider in the following subsection. Significantly, we consider a two-dimensional quantum system with spin-orbit coupling, i.e., particles whose momentum operator is coupled to their spin state. The Hamiltonian (with = 1 henceforth) of the system is:Ĥ where σ ·∇ ≡ σ x ∂ x +σ y ∂ y is the spin-orbit coupling term and are the Pauli spin matrices. Of these, only (σ x , σ y ) are associated with the two-dimensional coordinate frame. In the quantum system, if κ → ∞ thenĤ becomes the two-dimensional Dirac Hamiltonian in the Weyl representation (up to a global energy shift). Dimensionality plays an important role in this description: the one-dimensional Dirac equation is given byĤ with κ → ∞ and for Ψ independent of y. Rotated into imaginary time, the one-dimensional case yields a description of the time-dependent probability distribution of a persistent random walker. The stochastic process corresponding to such a walker is a Poisson process, and the one-dimensional Dirac equation in imaginary time can be restated as the telegrapher's equations describing this process. Such walkers move along a line at constant speed and with a turning rate (i.e., a rate for changing direction) given by m. This analogy has its origins in the path-intergral formulation of the Dirac equation [32], also restated in imaginary time in Refs. [33][34][35].
Although we develop some intuition by considering the quantum side of the analogy, we mostly focus on the mathematical description of self-propelled particles by performing a rotation of time into the imaginary axis: t → it. One of our main conclusions is that the (imaginary-time Schrödinger) equation in two dimensions, withĤ given by Eq. (3), describes the time-evolution of the probability distribution Ψ(x, t) of a self-propelled particle subject to two sources of noise: translational noise (controlled by the strength of the diffusion constant 1/κ), and rotational noise in the orientation angle [controlled by the parameter m in the term m(I − σ z )t]. This latter term describes the ability of active particles to change their direction of motion, as it does for particles in the one-dimensional Dirac equation [33][34][35]. We also show that, unlike in the one-dimensional case, the two-dimensional Dirac Hamiltonian does not consistently describe the probability density of a single, self-propelled particle: to ensure a physical description, κ must be constrained to be less than 8m. We now proceed to derive these results.
A. Spin, rotation, and velocity Let us now discuss the motivation for using spin as the carrier of information about the direction of propulsion. We use the spinorial representation of the rotation group, in which a rotation by angle θ around the axisn is generated by the unitary operatorÛn ,θ ≡ exp(−iθn · σ 3 /2), and where we have used the spin vector σ 3 defined by σ 3 ≡ (σ x , σ y , σ z ). We are considering particles constrained to be in the xy-plane and, therefore, all rotations are around the z-axis: σ z generates this abelian rotation group. These rotation operators are given bŷ As the global phase e −iθ/2 does not change the physical quantum state, we redefine the operator aŝ The action of the rotationÛ θ on the spinor (a, b) transforms it into the spinor (a, be iθ ). Note that the second component is rotated in the complex plane by the angle θ. Thus, the phase of this spinor component can be interpreted as the orientation of a polar particle, i.e., a particle that carries information about its orientation. Without loss of generality, we choose a global phase such that the first component of the spinor is real. Then, the spinor describing a particle oriented alongn = (cos θ, sin θ) is given by with s 1 real. The particle orientation is then given in terms of the real ( ) and imaginary ( ) parts of s 2 bŷ n ≡ ( s 2 , s 2 )/|s 2 |. We use spinors to describe the stochastic nature of motion and orientation for a self-propelled particle. We aim to capture quantities characterizing this particle such as the (scalar) probability ρ(r) of finding the particle at position r (regardless of its orientation) and corresponding probability current j(r). We show that the spinor encodes this information via ρ = (s 1 ) and j = ( s 2 , s 2 ). We can then construct the probability density P (r, θ) for finding the particle at position r and oriented along angle θ via where v(θ) ≡ (cos θ, sin θ) and s 2 = ( s 2 , s 2 ). For a self-propelled particle, motion and orientation are coupled. Using the probability density in Eq. (8), we note that the average of the particle's orientation at position r is proportional to j. We can conclude that the eigenvectors of the x-and y-components of the operator σ = (σ x , σ y ) correspond to particles oriented along the x-and y-directions, respectively. We then note that σ ·∇ is reminiscent of the convective derivative v·∇: σ·∇ convects the probability density in the direction along which the spinor points. Alternatively, this relation between σ and the velocity operator can be gathered directly from the Heisenberg equation of motion. If we consider the simplified HamiltonianĤ = σ · ∇, we can express the time-derivative of the position operator r (i.e., the velocity operator) as As a result of the additional terms in the Hamiltonian (3) of the model we consider, the translational noise contributes along with σ to the velocity operator.

B. Probabilistic interpretation
We now proceed to demonstrate the link between (imaginary-time Schrödinger) Eq. (4) and the probability densities and currents of self-propelled particles. To do so, we decompose the spinor Ψ into real and imaginary parts: where ρ, χ, j x , and j y are real-valued functions of the position r and time t (and are independent of θ). With this parametrization, Eq. (4) becomes where we have introduced j ≡ (j x , j y ) and, for any vector a, a ⊥ ≡ (a y , −a x ). Note that the first of these equations can be interpreted as a continuity equation, with ρ taking the role of a density and with both j and ∇ρ contributing to the current. Furthermore, we can interpret χ as a gauge degree of freedom for the orientation of the local coordinate frame. We make the simplest choice of gauge: χ = 0. Substituting this condition into Eqs. (11), we find We check the consistency of our gauge choice by noting that if the initial conditions satisfy Eq. (13), the evolution given by Eqs. (12,14) remains consistent with Eq. (13).
Indeed, we find this consistency condition to hold by applying ∇ ⊥ to Eq. (14): In what follows, we parametrize the velocity of selfpropulsion via v(θ) = (cos θ, sin θ), and use Eq. (12) to show that the system of equations (12)(13)(14) is equivalent to a Fokker-Planck equation that describes the dynamics of the probability density P (r, θ) ≡ ρ + v(θ) · j. Physically, P describes the probability of having a particle near r and oriented at an angle near θ. We first decompose the current into components parallel to and perpendicular to the velocity v via j = (v · j)v + (v ⊥ · j)v ⊥ and substitute this identity into the continuity equation (12). To get the dynamics of the distribution of the orientation angle θ, we multiply Eq. (14) by v. We add this equation for the time evolution of the current to the continuity equation to find an equation for the probability density P : (16) The two terms that are not yet expressed in terms of P can be addressed in the following way. First, note that v · j encodes orientational diffusion: We also find an extra component to diffusion that couples translational noise, rotational noise, and convection. This can be obtained using the identity v ⊥ · j = −∂ θ P : By substituting all of the diffusive and convective terms in Eqs. (17)(18) into Eq. (16), we arrive at the Fokker-Planck equation for P : From this formulation, we may note that m plays the role of a diffusion constant for the orientation θ, and κ −1 plays that role for the position r. The derivation of Eq. (19) is one of our main results: we showed that the Fokker-Planck equation (19) is equivalent to the imaginary-time Schrödinger equation (4) with the Hamiltonian (3), which includes a spin-orbit coupling term.

C. Microscopic Langevin equation
We now derive the precise microscopic model that corresponds to the Fokker-Planck equation (19). Before doing so, first note that for a Fokker-Planck equation to represent a stochastic microscopic model, the associated diffusion matrix must be positive-definite. In this subsection, we derive and analyze the drift vector and diffusion matrix for the model defined by Eqs. (3)(4), and derive the conditions under which an underlying microscopic model exists. To begin this analysis, we first define z ≡ aθ, where a is a particle lengthscale, and we re-express θ in terms of z to ensure that the different entries in the diffusion matrix have the same dimensionality. We compare Eq. (19) with the usual Fokker-Planck equation, viz., in which µ is the drift vector and D is the diffusion matrix. We then read off as follows: using a vector notation in which the third component corresponds to θ and As mentioned above, in order to construct a particlebased model for Eq. (19), it is necessary that D be positive-definite. To check this, we examine its eigenvalues (λ 1 , λ + , λ − ): κ, a, and m real and positive guarantees that λ 1 and λ + are positive. From the expression for λ − we conclude that D is positive-definite if and only if 8m > κ.
Thus, in order for Eq. (19) with m > 0 to correspond to a microscopic model, translational noise must be present, because without translational noise (i.e., in the limit κ → ∞) there isn't an orientational noise parameter m that satisfies Eq. (25). This implies that, in two dimensions, the Dirac equation alone cannot describe a self-propelled particle, in contrast to the one-dimensional case. As mentioned in the introduction to this section, the one-dimensional (imaginary-time) Dirac equation is equivalent to the telegrapher's equations, which do describe a stochastic process that can be interpreted as self-propulsion in one dimension. From Eq. (19), further insight into the relationship between m and κ and their physical interpretations can be gained by deriving the connection between this Fokker-Planck equation and the underlying microscopic process, i.e., the stochastic Langevin equation In Eq. (26), R t has N components (corresponding to the random variables), µ is an N -component associated drift vector, Σ(R t , t) is an N × M matrix, and W t is an M -dimensional Wiener process interpreted in either the Itô or Stratonovich sense [36,37]. For example, in the present case of N = 3, R t consists of the position vector R = (x, y) and orientation angle θ.
Formally, this interpretation can be established for an Itô process by considering a diffusion matrix D of the form ΣΣ T . Notice that if one such Σ 0 furnishes this decomposition then, for any orthogonal matrix R, Σ = Σ 0 R also satisfies it, so there are many different Langevin equations that yield the same Fokker-Planck equation. By using a Cholesky decomposition [38], we find a particular solution for the case M = 3: which yields the following Langevin equation: Here, ξ = (ξ 1 , ξ 2 ) is the two-dimensional translational noise that acts on the position of the particle, whereas ξ 3 is a rotational noise influencing the polarization angle θ. Note that to interpret this microscopic model we have assumed that Eqs. (27,28) are Itô stochastic differential equations. Generally, this differs from a Stratonovich process by an extra, noise-induced, drift vector having components However, in the case we are considering, the corresponding term is identically zero, and thus Eqs. (27)(28) can also be seen as a Stratonovich stochastic differential equation.

D. Noise and the uncertainty principle
Let us now discuss the physical picture of the singleparticles dynamics described by Eqs. (27)(28). This microscopic model has similarities to the models of active particles used, e.g., in Refs. [1,39]. At each instant in time, a particle is oriented at an angle θ and attempts to propagate in this direction at a constant speed. However, translational noise can change the direction of propagation away from the particle polarization. As a unique feature, the model we consider has feedback between translational and rotational noise: the larger the translational noise, the weaker the rotational noise. Quantitatively, if we define α ξ to be the angle that the force from the translational noise, ξ = (ξ 1 , ξ 2 ), makes with respect to the x-axis, we can rewrite the rotational noise term From Eq. (30), we observe that particles in effect try to oppose the translational noise, and prefer to align opposite to the direction of each kick. This coupling acts as a guidance system: in the absence of translational noise, the particle does not know which way to point. This is a consequence of how κ enters Eq. (28): the translational noise strength is inversely proportional κ, whereas the rotational noise strength is proportional to κ. Thus, the particle depends on feedback from translational noise to decide where to go. (A curious analogy emerges from the physics of hair cells in the inner ear, which depend on the presence of external noise to complete their function [40].) To further examine this feedback feature, consider the extreme case in which m is only slightly bigger than the lower-bound of κ/8, i.e., m = κ/8 + δ, with δ/κ 1. If δ is sufficiently small then rotational noise becomes irrelevant, compared to the large translational noise, and Eq. (28) becomes In this regime, the noise dominates over the self-propelled aspect of particle motion. The strength of this noise, quantified by κ −1 , is not subject to any restrictions, and both large-and small-noise regimes are physically accessible. Curiously, the interplay between the strength of the translational noise and its effect on the polarization is an expression of an uncertainty principle in this model. To see this, disregard the drift, and consider the feedback on the angle as simple additive noise. One then finds [R(t)−R(0)] 2 ∼ 4t/κ and [θ(t)−θ(0)] 2 ∼ tκ/2, which suggest the relation: This is a direct analog of the Heisenberg uncertainty principle, which relates the uncertainty of the position and velocity (here captured by orientation θ) of a quantum particle.
We now compare this microscopic model with others discussed in the literature. One related example involves a system composed of self-propelled hard rods that also experience translational noise, as examined in Refs. [41][42][43]. In these models, the translational and orientational noises are assumed to be uncorrelated. A situation closer to ours is explored in Ref. [44], in which the translational noise affects both orientational and spatial diffusion. In that case, the effects of these correlations have been examined in the inertial regime, in which Fokker-Planck dynamics are not equivalent to the imaginarytime Schrödinger equation that we examine here.
To conclude this section, let us generalize this interplay between translational and rotational noise and give it an arbitrary strength. In this case, Eq. (28) acquires an additional arbitrary (real) parameter λ via: (with m > κλ 2 /8). This parameter λ controls the response of a self-propelled particle to translational noise. The sign of λ determines the type of response: for λ < 0, the particle turns in the direction of any translational kick, whereas for λ > 0, as in the case above, the particle reacts in opposition to the kick. The Fokker-Planck equation associated with the Langevin dynamics of Eq. (27) is given by Although the additional parameter λ introduces more flexibility into the model, it destroys the uncertainty principle (32) and the bridge with the Schrödinger equation describing a two-component spinor.

II. HYDRODYNAMICS OF ACTIVE SPINS
In the previous section, we concluded that the Schrödinger equation for a model with a spin-orbit coupling term can be interpreted as an equation for the probability density of a self-propelled particle. In this section, we start with the many-body version of such a model, and go on to derive the coarse-grained hydrodynamic description of this polar active fluid.
In the noninteracting limit, the N -body Schrödinger equation can be written in terms of the elementary gen-eralizationĤ N of the one-body Hamiltonian (3), given byĤ where the summation i = 1, . . . , N is performed over the N particles. For noninteracting particles, we extract the probabilistic interpretation of the many-body wavefunction by noting that the probability density for many independent processes must obey P N (r 1 , θ 1 ; . . . ; r N , θ N ) = N i=1 P i (r i , θ i ). Therefore, the one-particle quantities ρ and j can be written as On the other hand, the many-body spinor associated to P N has 2 N components and the following structure: , as in the case of many non-interacting and uncorrelated quantum particles. Notice that the probability density P can capture more details regarding the distribution of the angular variable θ than the components of the spinor Ψ can encode. Indeed, the spinorial description assumes that only the first two Fourier modes in the angle θ are relevant, and disregards all higher Fourier components. Thus, in order to reduce a description in terms of P to one in terms of Ψ, the quantities 2π 0 e inθ P (θ)dθ must be negligible for all |n| ≥ 2.
Before we go on to include many-particle interactions, we generalize the route that took us from the Schrödinger equation (4) to the Fokker-Planck equation (19) to include multiple particles. We then arrive at the following many-particle Fokker-Planck equation: In order to uncover collective phenomena, we need to consider inter-particle interactions. For example, let us consider the alignment interaction typical of, e.g., the XY model, which can be included via a potential in the many-particle Langevin equation (28): wherein ({r, θ}) is a shorthand notation for (r 1 , θ 1 ; . . . ; r N , θ N ). In Eq. (38), the inter-particle interactions are encoded in the potential V i , which is defined via and acts with interaction strength g. The inter-particle separation enters the interaction potential via the function R(r i − r j ), which includes the characteristic range of the interactions. We now add this two-particle interaction to the many-particle Fokker-Planck equation (37) Equation (40) describes the time evolution for the probability distribution of many interacting self-propelled particles.

A. Self-consistent approximation
Instead of trying to exactly solve Eq. (40), in the present section we introduce a self-consistent approximation. To do this, we rewrite the potential as where The defining assumption of the self-consistent approximation is that where the average · · · C is taken with respect to the conditional probability that one of the particles is at position r i and oriented along θ i : This approximation treats the inter-particle interaction as an external potential due to the average effect of all the other particles. An explicit computation of the conditional average in Eq. (43) leads to (46) and the integral is taken over the two-dimensional area A. For simplicity, we now consider a purely local interaction [that is to say, taking R(r i − r j ) → δ(r i − r j )], in which case the self-consistency condition reduces to the simple form: 2π 0 dθ j e iθj P j (r i , θ j ). (47) Within the self-consistent approximation, all particles are identical and experience the same forcing. This forcing is, in turn, determined by considering the effect of a particle on its neighbors. The assumption of identical particles leads to all particles having the same probability distributions for all observables. In terms of probabilities, we thus have P j (r, θ) = P (r, θ) for all j. By using Eq. (47) and taking the N 1 limit, we obtain Substituting Eq. (48) into the expression for the potential, we find that the self-consistent potential has the form: For convenience, we rewrite this expression using an external alignment field h(r), defined via In terms of h, the potential has the form The self-consistent alignment field satisfies |h(r)| = h(r) and h(r) = πN j, i.e., it is a measure of the spontaneous alignment between the particle velocities. The advantage of using this self-consistent approximation is that it reduces the many-body Fokker-Planck equation to the one-particle non-linear equation, i.e., In the present subsection we have restricted ourselves to considering a description of interacting, selfpropelled particles in terms of the probability density P N rather than in terms of the Schrödinger equation. This is done out of necessity: the external potential term g ∂ θ [(h ⊥ · v)P ] in Eq. (52) cannot be captured within the Hamiltonian (3). To demonstrate this impossibility within a concrete example, let us consider a term in the Hamiltonian of the form h ⊥ · σ as a possible candidate. Such a term presents two issues that cannot be overcome within the framework we are considering: (i) Such a term generates a nonzero value of χ (in the imaginary part of the spinor). This issue can be overcome if one considers a more general framework in which the space of quantum states includes four-spinors with the structure Ψ = (φ,φ) as well as by including additional terms in the Hamiltonian (3). (ii) More significantly, the external potential term in Eq. (52) couples the lowest two Fourier modes of the orientation to higher Fourier modes. As a result, a description based on only the first two modes does not form a closed system of equations. We thus conclude that, in general, the Schrödinger equation in imaginary time with a spin-orbit coupling term describes single-particle dynamics only.

B. Onset of alignment
Although the spinorial description works for singleparticle dynamics only, we can use the Fokker-Planck description to examine the stability of the interacting isotropic active gas. In this subsection, we explore the onset of alignment due to inter-particle interactions. We follow the standard approach based on the dynamics of Fourier modes of the distribution of orientations θ [21].
First, we expand the single-particle probability density in Fourier modes: where v n ≡ (cos[nθ], sin[nθ]) and j n are the vectors whose components are the distinct Fourier modes of P, i.e., This expansion is similar to the one introduced in Ref. [21]; while the one in Ref. [21] is a decomposition using complex numbers, the one used above is purely real. Substituting Eqs. (54a, 54b) into Eq. (52) and using the linear independence of the Fourier components leads to the following set of coupled equations describing the time-evolution of the 3 lowest Fourier modes: where for compactness we have introduced the notation for the * product of two vectors, defined via a * b ≡ a x b x − a y b y . From Eqs. (55a-55e) we explicitly see that the interaction terms (which are proportional to g) couple the higher-order Fourier modes to the lowest ones. Nevertheless, notice that the dependence of j 2 on j is of higher order in the nonlinearity: h ∝ j, and the interaction term is quadratic in the currents. Thus, we may deduce whether the isotropic phase is stable by performing a linear stability analysis in which we assume that the current density |j| is small compared to the particle density ρ. Then, all higher Fourier modes, such as j 2 , may be neglected, and we obtain a linearized theory. As this approach neglects all stabilizing nonlinear terms, it does not yield a description of the polar active phase-we leave that task to the following subsections. We thus proceed with examining the stability of the isotropic phase, while neglecting all non-linear terms. The linearized equations are In order to study the stability of the isotropic phase within Eqs. (56a, 56b), we first look for solutions of the form (ρ, j) = (ρ 0 , j 0 ) e λt . We take spatial Fourier transforms, which re-express the gradient terms through the wavevector k ≡ (k x , k y ). These steps allow us to transform the above differential equations into an eigenvalue problem, wherein ρ 0 and j 0 act as eigenvector components and λ as an eigenvalue. The stability of the solutions of this system can then be analyzed by look-ing at the sign of the eigenvalues for each value of k. Specifically, there are three eigenvalues associated with the right-hand side of Eqs. (56a, 56b): where k ≡ |k|. From Eqs. (57a, 57b), we note that for small wavenumbers (or, equivalently, long wavelengths) the eigenvalues have a real part that is negative only if m > gN/2A. Defining the areal number density of the particle system via ρ d ≡ N/A, we may rewrite this condition as Physically, this condition corresponds to the regime in which the disordered (i.e., isotropic or non-polar) fluid phase is stable with respect to small fluctuations. In fact, the same condition appears in the analysis of the XY and Kuramoto models of, respectively, two-dimensional spins and synchronizing oscillators; see, e.g., Refs. [45,46]. In the special case of maximally large orientational noise (i.e., m → κ/8), this condition implies that This counterintuitive result states that the translational noise must be large in order for the disordered phase to be stable. This is a special feature of the model we are considering, in which the coupling between translational and orientational noise (and the resulting "uncertainty relation") requires a large translational noise if the orientations of the particles are to be distributed narrowly.

C. The polar active liquid
In the following subsection, we examine the behavior of long-wavelength fluctuations in the vicinity of the transition between the polar and disordered phases of the active fluid. Let us first consider the case in which the disordered phase is highly unstable, and show that the system's preferred spatially homogeneous state has rotational-symmetry-broken polar order. As a polar fluid maintains translational symmetry, all spatial derivatives in the Fokker-Planck equation (52) vanish, and the equation reduces to In this homogeneous case, pressure and translational noise are neglected, and the only obstacle to polar order is the orientational noise. The exact solution of this equation takes the form P = 1 AZ exp − g m h · v , in which Z is a normalization constant. Although this is an exact solution, we obtain a simpler expression by only considering the case in which |j| is small. It is more illustrative, though, to operate in terms of a Landau theory: in this regime, we may consider only up to the first three Fourier modes, which brings us to the following systems of equations: Solving these equations, and recalling that h = πN j, leads us to the following self-consistent equation for the alignment vector h: where h = |h|. Solutions of Eq. (62) are either the trivial h = 0 solution which is unstable for gρ d < 2m , or the solution This solution is only physical for gρ d ≥ 2m. Thus, this inequality presents the boundary between the disordered fluid phase and the ordered, polar, phase of the active fluid. We show this boundary in Fig. 2. Note that although in this subsection we have neglected all spatial derivatives, it is crucial that we have included the contributions of the higher Fourier modes in the orientational angle, i.e., j 2 . The coupling of this mode to the lower Fourier modes introduces nonlinearity in the system, which is necessary to stabilize the ordered polar phase.
The self-consistent approximation is known not to hold for phase transitions in which fluctuations are strong. Fluctuations can shift the transition line in parameter space and, furthermore, can change the nature of the transition. We thus expect the prediction that the orderdisorder transition in polar active fluids is continuous not to hold-in physical polar active fluids, the transition may be discontinuous.

D. Hydrodynamic equations: coarse-graining to
Toner-Tu theory Let us now go beyond the assumption of spatial homogeneity and examine the effects of temporally slow, longwavelength fluctuations in the density and the current for the case of a polar active fluid. This yields the hydrodynamic theory which, as for any polar active fluid, reduces to a form of Toner-Tu theory [31,47]. We proceed via the method introduced in Ref. [21]. Starting from Eqs. (55a-55e), we go to the regime in which the order is weak and therefore P (r, θ) depends only weakly on θ. Physically, this is the regime in which the active flow is much slower than the microscopic speed of each particle. Mathematically, this regime allows us to discard all Fourier modes higher than j 2 . [NB: For the sake of generality, let us also include the parameter λ from Eq. (33) which controls the feedback between the translational noise and the angle.] Then, we consider the equations Note two important aspects of the associated approximation: (i) |j 2 | is much smaller than |j|; and (ii) as we are interested in hydrodynamics, we only consider timeand length-scales much larger than the microscopic ones. Therefore, we consider the regime characterized by ∂ t j 2,x , ∂ t j 2,y mj 2,x , mj 2,y in which both the time-derivatives ∂ t j 2 and the term 1 κ ∇ 2 j 2 are neglected. Re-expressing the current j via j = h/πN and ρ d = 2πN ρ, we can rewrite the above equations as and Finally, the unusual h ⊥ terms may be rewritten in a more familiar way via the identities where, as above, h ≡ |h|. The resulting hydrodynamic equations read: Equation (66) is the continuity equation that established the conservation of the number of self-propelled particles.
On the other hand, Eq. (69) describes the hydrodynamics of the polarization order parameter for the active fluid. Note that the order parameter h is not identical to the current of the active fluid, which in addition includes the part of the motion due to translational diffusion. Let us now discuss the result for the hydrodynamics, Eq. (69). The second term on the left-hand side is an advection term. In general, as consequence of the breaking of Galilean invariance in our system, its prefactor is not unity. This results from the presence of a special reference frame in the system: there is only one frame in which the particles can propagate with equal speed, regardless of their direction of movement. On the righthand side, the first term is the one responsible for the spontaneous breaking of the rotational symmetry of the fluid in the polar state. The second and third terms are pressure-like terms. The third term includes the effects of a non-linear compressibility.
Turning now to the fourth term, this is a viscous damping term familiar from, e.g., the Navier-Stokes equations. In the present context, this term includes the effects of both the translational noise and the coupling between the translational and orientational noises. Curiously, the parameter λ has a strong effect on the effective viscosity. E.g., if λ = 1 (the model we focus on for the singleparticle dynamics), this viscosity reduces to 1/2κ. In that case, the viscosity depends only on the translational noise. At first glance, the expression for the effective viscosity does not appear to be positive-definite. If the viscosity were to change sign, this would suggest that the homogeneous fluid state should become unstable. However, in the previous section we showed that the particle dynamics is physical only if the inequality m > κλ 2 /8 holds. We now re-express the mass in terms of a positive parameter , via m = κ[ + (λ 2 /8)], and write the effective viscosity as .
The effective viscosity as given by Eq. (70) is positive for λ real and nonnegative. To see this, note that when 2 − λ − λ 2 < 0, Eq. (70) yields the minimal possible value for the effective viscosity for each value of λ in the limit → 0. This value is always positive and converges to 1/4κ as λ → ±∞. On the other hand, when 2 − λ − λ 2 > 0, the minimum value that the effective viscosity can achieve is 1/2κ. We plot the minimum value of the effective viscosity for each value of λ in Fig. 3. Curiously, it is also worth noting that in the case λ = −1 (a noise coupling with the same strength but with the opposite sign of our original model), the advection term vanishes.
The hydrodynamic equations (69), along with the continuity equation (66), are a version of the Toner-Tu equations, first derived in Ref. [31] based on symmetry considerations. By contrast, we obtain these equations based on the microscopic single-particle model that we introduced in the previous section along with inter-particle interactions. The coefficients that we obtained thus explicitly depend on the microscopic parameters of the model, which allows not only for the form but also for the precise numerical evaluation for the hydrodynamic coefficients in the Toner-Tu equations. In the special limit (κ, λ) → (∞, 0), with λ 2 κ → 0, the hydrodynamic equations that we derive coincide exactly with the ones derived in Ref. [47].

III. CONCLUSIONS
Because the model we consider falls within the Toner-Tu universality class, we can immediately conclude that The parameter λ controls the coupling between orientational and translational noises. Note that the effective viscosity is always positive, and the minimal value has a plateau for −2 ≤ λ ≤ 1. For |λ| 1, ν min eff asymptotes to 0.25. For λ > 0, this asymptote is approached from below, whereas for λ < 0 the asymptote is approached from above. This plot shows how the coupling between orientational and translational noise of active particles can control the effective viscosity of an active fluid.
the polar active particles that we consider do exhibit long-range order despite the two-dimensional character of the system [31]. This is in contrast to equilibrium two-dimensional systems with short-ranged interactions, which cannot break a continuous symmetry, such as the rotational symmetry of the XY model, as embodied in the Mermin-Wagner theorem [48]. Therefore, the hydrodynamic theory that we obtain, described by Eqs. (69), differs from the regular Navier-Stokes equations in two crucial ways: (i) due to the breaking of Galilean invariance via both activity and momentum exchange characteristic of dry active matter, the hydrodynamic theory includes terms prohibited in the Navier-Stokes equations; and (ii) as a result of these additional terms, interacting self-propelled particles exhibit long-range polar order.
To summarize, we have introduced a model of active particles based on an analogy with a Schrödinger equation that describes the propagation of an electron subject to spin-orbit coupling. We show that this model has a standard description as a stochastic process in terms of either a Fokker-Planck or a Langevin equation, both of which we derive. We note that within this stochastic interpretation, the orientational and the translational noise of the active particles we consider are coupled via a relation reminiscent of the Heisenberg uncertainty principle. Based on this single-particle physics, we derive a description for a polar active fluid in which the particles preferentially align their velocities. Within this description, we characterize the transition from a disordered to an ordered (i.e., polar) state via a hydrodynamic Toner-Tu theory.
Previously, analogies between classical processes and quantum dynamics have found use in areas as diverse as polymer physics, liquid crystal elasticity, hydrodynamics, and financial markets [49]. Such analogies are often drawn via a path-integral formalism but they may instead be formulated by rotating the time axis into the complex plane. We have shown that this latter approach can be extended to the study active fluids composed of self-propelled particles. In order to account for the selfpropulsion, we employ concepts familiar from the study of correlated electron fluids. These connections have the potential to help uncover novel phases of active matter via analogies with electronic counterparts.

IV. ACKNOWLEDGMENTS
We thank D. Bartolo, J. Toner, V. Vitelli, and R. Hipolito for insightful discussions. B.L. would like to thank CONICYT and Becas Chile for their financial support. A.S. acknowledges funding from the Delta Institute for Theoretical Physics. P.M.G. acknowledges for its hospitality the Aspen Center for Physics, which is supported by National Science foundation grant PHY-1607611.