Arxiv:2009.00118V1 [Physics.Flu-Dyn] 31 Aug 2020 Ary Layer, That Is, Νe → 0 Or Cs → ∞, Both the ﬂuid the Negative Sign Refers to the Left-Moving Ones

Non-linear shallow water dynamics with odd viscosity

Gustavo M. Monteiro1 and Sriram Ganeshan1 1Department of Physics, City College, City University of New York, New York, NY 10031, USA (Dated: September 2, 2020) In this letter, we derive the Korteweg-de Vries (KdV) equation corresponding to the surface dynamics of a shallow depth (h) two-dimensional ﬂuid with odd viscosity (νo) subject to gravity (g) in the long wavelength weakly nonlinear limit. In the long wavelength limit, the odd viscosity term plays the role of surface tension albeit with opposite signs for the right and left movers. We show that there exists two regimes with a sharp transition point within the applicability of the p 3 p 3 KdV dynamics, which we refer to as weak (|νo| < gh /6) and strong (|νo| > gh /6) parity- breaking regimes. While the ‘weak’ parity breaking regime results in minor qualitative diﬀerences in the soliton amplitude and velocity between the right and left movers, the ‘strong’ parity breaking regime on the contrary results in solitons of depression (negative amplitude) in one of the chiral sectors.

Introduction: Parity-breaking phenomena in two- modified pressure term at the edge. dimensional fluids such as odd viscosity effects have been In short, for an irrotational flow, that is, v = ∇θ, at the center of investigation in diverse platforms. For the effective boundary dynamics can be expressed as two dimensional isotropic and non-dissipative flows, odd a Laplace equation for the velocity potential θ in the viscosity is the only non-vanishing coefficient of the vis- bulk and Bernoulli’s equation at the boundary with the cosity tensor and has a geometrical connection to the modified pressure [40]. The variational principle for this adiabatic curvature on the space of flat background met- boundary dynamics was later formulated in terms of a rics [1, 2]. Examples of quantum systems where odd vis- geometric action which resulted in the odd viscosity in- cous effects are important include electron fluids in meso- duced effective pressure at the boundary [41]. In the scopic systems [3–5], quantum Hall fluids [6–26], and chi- limit of infinitely deep fluid the weakly non-linear dy- ral superfluids/superconductors [27]. namics within a small angle approximation was shown to Chiral-Burgers In classical fluids, odd viscosity shows up in polyatomic be governed by the novel equation [40]. gases [28–31], chiral active matter [32–34], vortex dy- In this letter, we study the shallow depth limit of namics in 2D [35–38] and chiral active fluids [32–34]. the weakly non-linear surface dynamics with odd vis- For incompressible flows, it has been shown by one of cosity and gravitational force (confining potential) (see the authors that odd viscosity effects are absent when the schematic in Fig. 1). We assume that the bound- the fluid is spread on the entire plane or confined in ary layer is the shortest length scale and the effective rigid domains with no-slip boundary conditions [39]. In dynamics is irrotational, using the hydrodynamic equa- other words, the velocity profile is independent of the tions from [41] as the starting point. We show that, for odd viscosity. Nevertheless, the signature of this parity- later times and long wavelengths, the weakly non-linear breaking coefficient is present in surface waves and in the dynamics is given by the integrable Kortweg De Vries interface between two fluids governed by kinematic and (KdV) equation with the kinematic odd viscosity νo en- no-stress boundary conditions, which explicitly depends tering the coefficient of the dispersive term, on the odd viscosity. The dynamical surface problem ! in the presence of odd viscosity results in an oscillat- ηt 3 2 1 νo ηx ± √ + η ηx + h ± p ηxxx = 0. (1) gh 2h 6 gh3 ing boundary layer where√ the vorticity is confined within some thickness of δ ∝ νe [40, 41] (where νe is the kine- −1 matic shear viscosity) for the dissipative case and δ ∝ cs Here, η(x, t) is the boundary shape profile, h is the aver- (where cs is the sound velocity) for the non-dissipative age depth of the fluid and g is the acceleration of gravity. compressible case [42]. In the limit of a very thin bound- The positive sign refers to right-moving solitons whereas arXiv:2009.00118v1 [physics.flu-dyn] 31 Aug 2020 ary layer, that is, νe → 0 or cs → ∞, both the fluid the negative sign refers to the left-moving ones. The manifestation of odd viscosity in the above equation is pressure√ at the edge and the surface vorticity diverge as 1/ νe or cs, but the quantityp ˜ = p − νoρ ω remains fi- similar to that of the surface tension although with dif- nite. We refer top ˜ as modified pressure, where νo is the ferent signs for the left and right movers. odd viscosity and the variables p, ω and ρ are the fluid The odd viscosity entering the KdV equation has ma- pressure, vorticity and constant background density re- jor consequences due to its parity breaking effects. We spectively. This cancelation of divergences allows us to show that there exists two regimes with a sharp transi- write the dynamical surface problem with odd viscosity tion point within the applicability of the KdV dynam- p 3 as an effective irrotational system where all the effects of ics, which we refer to as weak (|νo| < gh /6) and p 3 odd viscosity and boundary layer can be absorbed into a strong (|νo| > gh /6) parity-breaking regimes. In the 2

y the fluid normal to the boundary is equal to the rate of x change of the boundary shape. The pair of dynamical 2a h boundary conditions imposes that there are no normal and tangent forces acting on an element of the fluid sur- y = ⌘(t, x) face. The presence of two dynamical boundary conditions along with incompressibility requires a singular boundary layer where the vorticity is confined. The role of this singular boundary layer is to ensure that there are Figure 1. Schematics of the shallow fluid dynamics with free no tangent forces on this interface. It has been shown surface. Here a denotes the amplitude of the edge profile and that regardless of the boundary layer mechanism (dis- δ is the boundary layer thickness. The vertical direction is sipative or compressible), the normal component of the exaggerated to highlight these features. DBC is universal and geometric in nature [40–42]. As- suming that the boundary layer is stable and is confined weak parity breaking regime, the left and right mov- to short length scales, the free surface problem can be ing solitons slightly differ in amplitude and speed. In written as an effective description of the fluid with the the strong parity breaking regime, one of the sectors be- effects of boundary layer encoded in the odd viscosity comes solitonic waves of depression. At the critical value modified pressure term, | | p 3 ( νo = gh /6), one of the sectors becomes unstable and 2νo p˜ y=η(t,x) = p ∂xvn . (4) higher order derivative terms become important. The 1 + η2 parity breaking KdV dynamics discussed here is in stark x contrast to the parity preserving case of shallow water Here vn is the velocity component which is normal to the KdV dynamics without odd viscosity. boundary, taken at y = η(t, x). Incompressible fluids with odd viscosity: The hy- This effective free surface dynamics, that is, Eqs. (2, drodynamic equations for incompressible fluids with odd 4) can be expressed in the form of an action principle, viscosity consist of the Euler equation, with modified as shown in [41]. For irrotational flows, i.e., v = ∇θ, pressure, together with the incompressibility condition, this action simplifies and can be thought of as an odd that is, viscosity extension of the Luke’s variational principle [43]. Following [41], the hydrodynamic action becomes 1 ∇ · v = 0 , ∂tv + (v · ∇)v = − ∇p˜ − gyˆ. (2) ρ ZZ Z η(t,x) 1 2 2 S = − dt dx dy θt + θx + θy Here, v is the flow velocity and ρ is the constant and uni- −h 2 ZZ form fluid density. This differs from the ordinary Euler 1 2 −1 − dt dx gη − νo ηt tan (ηx) , (5) equation in the definition of the modified pressure 2

p˜ := p − νoρ ω , (3) where η(t, x) as the top surface shape function. The domain of the fluid is bounded by a finite depth at the bot- where ω = ∇ × v is the fluid vorticity. tom. The bulk and boundary equations of motion are Since the fluid density is constant, there is no equation obtained by varying the above action for a finite depth of state and the pressure is completely determined by the fluid. The bulk equation for the irrotational system is flow. The curl of the Euler equation is an equation for simply the Laplace equation for the potential ∆θ = 0 the flow vorticity, which does not depend on the pressure. defined in the domain −h < y < η(t, x). The boundary This vorticity equation together with the incompressibil- conditions at the top and bottom of the fluid domain can ity condition completely determines the fluid flow, up to be written as, boundary conditions. Therefore, the presence of odd viscosity will only change the velocity flow if the boundary θy = 0 , y = −h, (6) conditions depend on νo, otherwise, it only modifies how ηt + ηxθx = θy , y = η(t, x), (7) the pressure depends on velocity flow, which is not an 2 2 " # easily accessible quantity in experiments. θx + θy 2νo ηt θt + + gη = ∂x , y = η(t, x). p 2 p 2 Irrotational limit of the free surface dynamics 2 1 + ηx 1 + ηx and boundary layer approximation: Bulk equations (8) of motion (2) must be supplemented with boundary conditions. The fluid free surface is a dynamical interface Eq. (7) is the kinematic boundary condition and Eq. (8) is y = η(t, x) between two fluids where we impose one the odd viscosity modified dynamic boundary condition. kinematic and two dynamical boundary conditions. The Linear waves in long wavelength limit: Before kinematic boundary condition states that the velocity of we dive into the non-linear dynamics of the shallow fluid 3 regime, let us focus on the linearized free surface problem the boost of the center of mass and it only appears in the with finite depth. In this limit, we drop all the quadratic combination X − T , for right-movers, or X + T , for left- terms in Eqs. (7, 8) and evaluate the derivatives of θ at movers. On the other hand, the cubic term scales as ε3/2 y = 0. For the monochromatic surface profile η(x, t) = and give rise to a dispersive group velocity of the boosted a cos(kx − Ωt) as an input, we find the velocity potential wavepacket. This effect becomes relevant at much later to be times, in comparison to T , and introduce a second time a Ω scale, which governs the time evolution of the boosted θ(x, y, t) = cosh[k(y + h)] sin(kx − Ωt), wavepacket and we denote by τ. This means that both k sinh(kh) variables θ and η evolve according to this double time with the surface dispersion relation given by scale, such that, the time derivative becomes h i r r 3 2 p 2 4 gε gε Ω = tanh kh −νok ± ν k + gk coth(kh) . (9) ∂t = ∂T + ∂τ . o h h In the absence of gravity (g = 0) and in the deep ocean In the following, we derive the full non-linear shallow limit (h → ∞), we have that tanh kh ≈ k/|k| and we water dynamics with odd viscosity and discuss how the recover the known odd viscosity dominated dispersion parity breaking effects manifest in the non-linear dynam- 1 p 3 Ω = {0, −2νok|k|}. The weakly non-linear dynamics for ics. In particular we show that νo = 6 gh manifests as this system was discussed in Ref. [40]. a critical point that separates two qualitatively different Shallow waves, on the other hand, arise when the fluid regimes of non-linear dynamics. depth h is much smaller than the characteristic wave- Non-linear shallow depth waves: Korteweg-de lengths of the system. In other words, they are charac- Vries (KdV) equation arises in the study of shallow water terized by kh 1. In this approximation, the leading waves of long wavelengths and small amplitudes. In the terms in the dispersion (9) are given by following analysis, we show how the KdV equation corresponding to the shallow depth limit is modified by the r " ! # g νo 1 3 presence of the odd viscosity term. The counting scheme Ω ≈ ±kh − p ± (kh) . (10) h gh3 6 for the KdV equation is chosen such that small amplitudes regime corresponds to η ∼ O(ε), that is, η is of 2 The first term is the usual shallow water gravity wave the same order as ∂x. Thus, we can rescale the boundary dispersion, whereas the second one is the odd viscosity shape in terms of h as η = εhy. modified Korteweg de-Vries dispersive term. Prima facie KdV regime happens for sufficiently later times, so that it seems that the odd viscosity is qualitatively similar to right-moving and left-moving solutions are independent the surface tension effect for the shallow water surface and well-separated. Here, we restrict ourselves to only dispersion. However, the physical manifestation of odd right-moving propagation, since the analysis for the the viscosity is completely different, since the coefficient of left-movers follows similarly. Hence, let us assume θ and the cubic term can develop a relative sign change between η of the form p the left mover and right mover for |ν | > 1 gh3, what o 6 p 3 indicates a strong parity-breaking phenomena. θ(t, x, y) = εh g ϑ(τ, σ, Y ; ε) , (11) The shallow wave condition naturally introduces a η(t, x) = εh y(τ, σ; ε) , (12) power expansion in kh. Formally, it is convenient to de- √ with σ = X − T . Under these conditions, the bulk equa- fine an expansion parameter ε 1, such that kh = ε k¯ tion of motion and the boundary condition at the flat and k¯ is a dimensionless wavenumber. Since an expan- bottom become sion in powers of kh can be translated into a derivative expansion for the fluid dynamics, this rescaling is equiv- ε ϑσσ + ϑYY = 0, −1 < Y < εy, (13) √h alent to the redefinition x = ε X, where X is the di- ϑY = 0,Y = −1. (14) mensionless horizontal coordinate. In the same way, we − can define y = hY , with Y being the dimensionless ver- Let us denote ϑ(τ, σ, 1; ε) by φ(τ, σ; ε). This way, tical coordinate. In this counting scheme, we have that the solution of Eq. (13) with the condition (14) can be 1/2 written as ∂x ∼ O ε , whereas ∂y ∼ O(1). ∞ The wave dynamics dictated by Eq. (10) evolves ac- X (−ε)n(1 + Y )2n ϑ(τ, σ, Y ; ε) = ∂2nφ . (15) cording to two distinct time scales. The linear term (2n)! σ √ n=0 scales with ε and governs the splitting of an initial dis- turbance into right-moving and left-moving wavepackets. Plugging Eq. (15) into Eqs. (7, 8) and neglecting terms 2 The first time scale, which we denote by T , is the char- of O(ε ) or higher, we obtain acteristic time in which left and right-movers are so far 1 2 1 φσ = y + ε φτ + φ + φσσσ + 2¯νoyσσ , (16) apart, we can study them separately. In other words, for 2 σ 2 y − 1 y − y 2 y sufficiently later times, the only role of T is to account for σ φσσ = ε 2 σσσ τ + 3 φσσσσ + ∂σ(φσ ) . (17) 4

p 3 Here, we denotedν ¯o = νo/ gh . Eq. (16) allow us to X+ perturbatively express φ in terms of y. Substituting this σ 2.0 expression into Eq. (17), the leading order equation for the right-moving surface wave in the boosted reference frame becomes 1.5

3 1 , τ ) , τ ) yτ + 2 y yσ + 6 +ν ¯o yσσσ = 0, (18) - 1.0 + η ( X η ( X which is nothing but the well-known KdV equation. In terms of the dimensionful variables t, x and η(t, x), it 0.5 becomes Eq. (1). As previously mentioned, we could repeat the same 0.0 analysis for the left-moving solitons. For ξ = X + T , we -40 -20 0 20 40 obtain X- 3 1 yτ − 2 y yξ − 6 − ν¯o yξξξ = 0. (19) Figure 2. Left (blue) and right (yellow) moving soliton in the 1 ¯ Under reﬂection about y-axis (parity operation), y → y, weak parity breaking regime with parameters |ν¯o| = 15 , k = 1 ξ → −σ and Eq. (19) becomes and τ = 1.

3 1 X yτ + 2 y yσ + 6 − ν¯o yσσσ = 0. (20) + 4 The odd viscosity term breaks parity symmetry of the problem [44], since the left-moving soliton under reﬂec- 3 tion about the y-axis do not behave like the right-moving soliton. The odd viscosity term entering the KdV equa- 2 , τ ) , τ )

tion is similar to the presence of the surface tension. - + Within this analogy of odd viscosity as surface tension, 1 η ( X η ( X the left mover and right mover will have opposite signs 0 of surface tension due to the parity breaking effects of odd viscosity. In other words, odd viscosity in the KdV -1 regime acts as chirality dependent surface tension term. Soliton solution: In the following we analyze the role -2 -40 -20 0 20 40 of odd viscosity in the single soliton solution of Eqs. (18) and (19). Although multi-soliton solutions also show the X- same qualitative behavior, they are out of the scope of Figure 3. Left (blue) and right (yellow) moving soliton in this letter. The single soliton solution corresponding to 1 the strong parity breaking regime with parameters |ν¯o| = 3 , the left and right movers can be written as, k¯ = 1 and τ = 1. ¯2 1 2 ¯ ¯3 1 y(τ, x±) = 8k 6 ∓ ν¯o sech kx± ± 4k τ 6 ∓ ν¯o , (21) has positive amplitude, whereas the other corresponds ¯ where k is the dimensionless wavenumber and we denoted to solitonic waves of depression or depletion as shown in x+ = ξ and x− = σ in order to shorten the notation. Fig. 3. 1 Moreover, τ = 0 was chosen such that the soliton center ‘Critical’ dynamics |ν¯o| = : At this critical points, 1 6 of mass is at x± = 0. Note that 6 ∓ ν¯o enters both the dispersive term in one of the sectors vanish and we the amplitude and the wave speed. Therefore, the odd end up with the inviscid Burger’s equation for such sector viscosity modification to the KdV soliton dynamics can as shown in Fig. 4. In fact, it is known that solutions of be separated into three regimes depending on the value the inviscid Burger’s equation are subjected to a blow up 1 of 6 ∓ ν¯o time, in which the spatial derivative of y becomes infinite 1 ‘Weak’ parity breaking regime |ν¯o| < 6 : In this case, and higher order derivative terms become important. 1 6 ∓ ν¯o > 0, with left and right moving solitons only In general, non-localized solutions for the KdV are differ in the magnitude of the amplitude and velocity as given in terms of the Jacobi elliptic function, also known shown in Fig. 2 . We refer to this as weak parity breaking as cnoidal functions, regime. 1 ¯2 1 2 ¯ ¯3 1 ‘Strong’ parity breaking regime |ν¯o| > 6 : In this case, y(τ, x±) = 8k 6 ∓ ν¯o cn kx± ± 4k τ 6 ∓ ν¯o ; κ . 1 6 ∓ ν¯o have opposite signs. We call this strong par- (22) ity breaking regime, because the difference between left The parameter κ ∈ [0, 1] interpolates between the long and right moving solitons is visual, that is, one sector linear waves for κ → 0 (cn → cos) and the single soliton 5

X+ [5] A. Berdyugin, S. Xu, F. Pellegrino, R. K. Kumar, 3.0 A. Principi, I. Torre, M. B. Shalom, T. Taniguchi, 2.5 K. Watanabe, I. Grigorieva, et al. Measuring hall viscosity of graphene’s electron ﬂuid. Science, page eaau0685 2.0 (2019). [6] I. Tokatly. Magnetoelasticity theory of incompressible

, τ ) 1.5 , τ ) quantum Hall liquids. Physical Review B, 73, 205340 - + (2006).

η ( X 1.0 η ( X [7] I. Tokatly and G. Vignale. New collective mode in the 0.5 fractional quantum Hall liquid. Physical review letters, 98, 026805 (2007). 0.0 [8] I. Tokatly and G. Vignale. Erratum: Lorentz shear modulus of a two-dimensional electron gas at high magnetic -0.5 -40 -20 0 20 40 field [Phys. Rev. B 76, 161305 (R)(2007)]. Physical Re- view B, 79, 199903 (2009). X- [9] N. Read. Non-Abelian adiabatic statistics and Hall viscosity in quantum Hall states and p x+ i p y paired su- Figure 4. Right (yellow) moving soliton in the critical parity perfluids. Physical Review B, 79, 045308 (2009). 1 ¯ breaking |ν¯o| = 6 , k = 1, and τ = 1. The left mover has zero [10] F. Haldane. Geometrical description of the fractional amplitude for this case. quantum Hall effect. Physical review letters, 107, 116801 (2011). [11] F. Haldane. Self-duality and long-wavelength behavior of solution for κ → 1 (cn → sech). the Landau-level guiding-center structure function, and the shear modulus of fractional quantum Hall fluids. Discussion and Outlook: In this letter, we derived arXiv preprint arXiv:1112.0990 (2011). parity broken generalization of the Korteweg de-Vries [12] C. Hoyos and D. T. Son. Hall viscosity and electro- equation for shallow depth fluid with odd viscosity and magnetic response. Physical review letters, 108, 066805 gravity in the long wavelength weakly non-linear limit. (2012). The presence of odd viscosity manifests weak and strong [13] B. Bradlyn, M. Goldstein, and N. Read. Kubo formulas parity breaking regimes in the two chiral sectors of the for viscosity: Hall viscosity, Ward identities, and the relation with conductivity. Physical Review B, 86, 245309 KdV dynamics. The odd viscosity term plays the role (2012). of surface tension albeit with opposite signs for the right [14] B. Yang, Z. Papić,E. Rezayi, R. Bhatt, and F. Haldane. and left movers. In future work, we aim to specialize Band mass anisotropy and the intrinsic metric of frac- this result to chiral active fluids, where odd viscous ef- tional quantum Hall systems. Physical Review B, 85, fects have been observed in free surface dynamics [34]. 165318 (2012). In order to make contact with experiments, we will nu- [15] A. G. Abanov. On the effective hydrodynamics of the merically study the Cauchy initial value problem of an fractional quantum Hall effect. Journal of Physics A: Mathematical and Theoretical, 46, 292001 (2013). initial perturbation that evolves into left and right mov- [16] T. L. Hughes, R. G. Leigh, and O. Parrikar. Tor- ing solitons and quantify conditions under which weak sional anomalies, Hall viscosity, and bulk-boundary cor- and strong parity breaking KdV dynamics emerges. respondence in topological states. Physical Review D, 88, Acknowledgments. We thank Alexander Abanov 025040 (2013). and Vincenzo Vitelli for helpful discussions and sugges- [17] C. Hoyos. Hall viscosity, topological states and effective tions about to this project. This work is supported theories. International Journal of Modern Physics B, 28, 1430007 (2014). by NSF CAREER Grant No. DMR-1944967 (SG) and [18] M. Laskin, T. Can, and P. Wiegmann. Collective field partly from PSC-CUNY Award. GM was supported by theory for quantum Hall states. Physical Review B, 92, 21st century foundation startup award from CCNY. 235141 (2015). [19] T. Can, M. Laskin, and P. Wiegmann. Fractional quantum Hall effect in a curved space: gravitational anomaly and electromagnetic response. Physical review letters, 113, 046803 (2014). [1] J. Avron, R. Seiler, and P. G. Zograf. Viscosity of quan- [20] T. Can, M. Laskin, and P. B. Wiegmann. Geometry of tum Hall fluids. Physical review letters, 75, 697 (1995). quantum Hall states: Gravitational anomaly and trans- [2] J. Avron. Odd viscosity. Journal of statistical physics, port coefficients. Annals of Physics, 362, 752–794 (2015). 92, 543–557 (1998). [21] S. Klevtsov and P. Wiegmann. Geometric adiabatic [3] T. Scaffidi, N. Nandi, B. Schmidt, A. P. Mackenzie, and transport in quantum Hall states. Physical review let- J. E. Moore. Hydrodynamic electron flow and Hall vis- ters, 115, 086801 (2015). cosity. Physical review letters, 118, 226601 (2017). [22] S. Klevtsov, X. Ma, G. Marinescu, and P. Wiegmann. [4] F. M. Pellegrino, I. Torre, and M. Polini. Nonlocal trans- Quantum Hall effect and Quillen metric. Communica- port and the Hall viscosity of two-dimensional hydrody- tions in Mathematical Physics, 349, 819–855 (2017). namic electron liquids. Physical Review B, 96, 195401 [23] A. Gromov and A. G. Abanov. Density-curvature re- (2017). sponse and gravitational anomaly. Physical review let- 6

ters, 113, 266802 (2014). surface dynamics in the absence of dissipation. arXiv [24] A. Gromov, G. Y. Cho, Y. You, A. G. Abanov, and preprint arXiv:1907.11196 (2019). E. Fradkin. Framing anomaly in the effective theory of [43] J. C. Luke. A variational principle for a fluid with a free the fractional quantum hall effect. Physical review letters, surface. Journal of Fluid Mechanics, 27, 395–397 (1967). 114, 016805 (2015). [44] A right moving soliton reflected about Y-axis is same as [25] A. Gromov, K. Jensen, and A. G. Abanov. Boundary this solution played backwards in time. effective action for quantum Hall states. Physical review letters, 116, 126802 (2016). [26] P. Alekseev. Negative magnetoresistance in viscous flow of two-dimensional electrons. Physical review letters, 117, 166601 (2016). [27] C. Hoyos, S. Moroz, and D. T. Son. Effective theory of chiral two-dimensional superfluids. Physical Review B, 89, 174507 (2014). [28] J. Korving, H. Hulsman, H. Knaap, and J. Beenakker. Transverse momentum transport in viscous flow of di- atomic gases in a magnetic field. Physics Letters, 21, 5–7 (1966). [29] H. Knaap and J. Beenakker. Heat conductivity and viscosity of a gas of non-spherical molecules in a magnetic field. Physica, 33, 643–670 (1967). [30] J. Korving, H. Hulsman, G. Scoles, H. Knaap, and J. Beenakker. The influence of a magnetic field on the transport properties of gases of polyatomic molecules;: Part I, Viscosity. Physica, 36, 177–197 (1967). [31] H. Hulsman, E. Van Waasdijk, A. Burgmans, H. Knaap, and J. Beenakker. Transverse momentum transport in polyatomic gases under the influence of a magnetic field. Physica, 50, 53–76 (1970). [32] D. Banerjee, A. Souslov, A. G. Abanov, and V. Vitelli. Odd viscosity in chiral active fluids. Nature Communi- cations, 8, 1573 (2017). [33] A. Souslov, K. Dasbiswas, M. Fruchart, S. Vaikun- tanathan, and V. Vitelli. Topological waves in fluids with odd viscosity. Physical Review Letters, 122, 128001 (2019). [34] V. Soni, E. Bililign, S. Magkiriadou, S. Sacanna, D. Bar- tolo, M. J. Shelley, and W. Irvine. The free surface of a colloidal chiral fluid: waves and instabili- ties from odd stress and Hall viscosity. arXiv preprint arXiv:1812.09990 (2018). [35] P. Wiegmann and A. G. Abanov. Anomalous hydrodynamics of two-dimensional vortex fluids. Physical review letters, 113, 034501 (2014). [36] X. Yu and A. S. Bradley. Emergent non-eulerian hydrodynamics of quantum vortices in two dimensions. Physical review letters, 119, 185301 (2017). [37] A. Bogatskiy and P. Wiegmann. Edge Wave and Bound- ary Layer of Vortex Matter. Phys. Rev. Lett., 122, 214505 (2019). [38] A. Bogatskiy. Vortex flows on closed surfaces. arXiv preprint arXiv:1903.07607 (2019). [39] S. Ganeshan and A. G. Abanov. Odd viscosity in two- dimensional incompressible fluids. Physical Review Flu- ids, 2, 094101 (2017). [40] A. Abanov, T. Can, and S. Ganeshan. Odd surface waves in two-dimensional incompressible fluids. SciPost Physics, 5, 010 (2018). [41] A. G. Abanov and G. M. Monteiro. Free surface variational principle for an incompressible fluid with odd viscosity. Phys Rev Lett, 122, 154501 (2019). [42] A. G. Abanov, T. Can, S. Ganeshan, and G. M. Mon- teiro. Hydrodynamics of two-dimensional compressible fluid with broken parity: variational principle and free