1 Introduction
Tools from topology developed over recent decades in condensed matter physics have recently shed new light on our understanding of fluid waves, from the design of microfluidic devices (Souslov et al.
Reference Souslov, van Zuiden, Bartolo and Vitelli2017), to acoustic waves (Yang et al.
Reference Yang, Gao, Shi, Lin, Gao, Chong and Zhang2015), to planetary atmospheres (Delplace, Marston & Venaille Reference Delplace, Marston and Venaille2017; Perrot, Delplace & Venaille Reference Perrot, Delplace and Venaille2018), as well as in active matter flows (Shankar, Bowick & Marchetti Reference Shankar, Bowick and Marchetti2017). It has been realized that important and robust information on the spectrum of a linear operator is encoded in the eigenmodes of this operator in an unbounded geometry, with constant coefficients. This information is revealed by the winding of the eigenmodes parameterized over a closed surface. This winding is a topological invariant, the Chern number, that can be explicitly computed. For instance, Delplace et al. (Reference Delplace, Marston and Venaille2017) showed that inertia-gravity waves in the rotating shallow water model exhibit such a topological property in
$(k_{x},k_{y},f)$
-space, with
$f$
the Coriolis parameter and
$(k_{x},k_{y})$
the wavenumber. More precisely, a Chern number of 2 was found for the positive-frequency inertia-gravity waves as they enclose the origin in parameter space, while the zero-frequency (geostrophic) modes carry a vanishing Chern number.
One spectacular manifestation of these singularities occurs when one computes the spectrum of the same operator, now assuming that one of the parameters varies spatially. In the shallow water model, such computations have for instance been performed by Matsuno (Reference Matsuno1966) on the equatorial beta plane, assuming linear variations of the Coriolis parameter in the
$y$
(meridional) direction. He discovered the existence of two branches in the dispersion relation that transit between different wavebands when the wavenumber in the
$x$
(zonal) direction is varied: the equatorial Kelvin wave and the mixed Rossby-gravity wave, now known as the Yanai wave. These modes are more localized along the equator than the others, and are unidirectional. The gradient of planetary vorticity also supports the propagation of low-frequency planetary Rossby waves that lift part of the degeneracy of geostrophic modes. However, in contrast to Kelvin and Yanai waves, Rossby waves remain in the geostrophic band. Thus, in the equatorial beta plane, when the zonal wavenumber is varied from negative to positive values, the positive-frequency inertia-gravity waveband has a net gain of two modes (Delplace et al.
Reference Delplace, Marston and Venaille2017). The correspondence between a topological invariant, the first Chern number, that describes degeneracy points for bulk waves in parameter space, on the one hand, and the number of modes that transit from one band to another in the equatorial beta plane, on the other hand, is reminiscent of the Atiyah–Singer index theorem (Bal Reference Bal2018; Faure Reference Faure2019) and could be referred to as a topological spectral flow correspondence. This correspondence has proved useful in interpreting molecular spectra (Faure & Zhilinskii Reference Faure and Zhilinskii2000), or Lamb-like waves trapped along an interface of density stratification (Perrot et al.
Reference Perrot, Delplace and Venaille2018), among other applications in physics (Nakahara Reference Nakahara2003).
There exists, however, the possibility for a stronger bulk-interface correspondence, when a topological index can be assigned to each wave band on each side of the interface, rather than to a band-crossing point in parameter space. The bulk-interface correspondence then predicts the number of unidirectional edge states trapped along the interface between the two regions as a consequence of the mismatch between the topological indices of two wavebands across the interface. In condensed matter, thanks to an underlying lattice structure, this bulk-interface correspondence naturally follows from a bulk-boundary correspondence that relates a topological index of the bulk to the number of unidirectional edge modes propagating along a boundary (Hatsugai Reference Hatsugai1993; Graf & Porta Reference Graf and Porta2013). In this context the interface is recovered by gluing together the boundaries of two topologically distinct systems. It is thus natural to ask whether equatorial Kelvin and Yanai waves can also be understood as topological interface states between two distinct topological systems. In other words, can a single hemisphere be interpreted as a topological media on its own? Here we use odd-viscosity terms to assign a well-defined topological invariant to each hemisphere, building on previous work by Volovik (Reference Volovik1988) and Souslov et al. (Reference Souslov, Dasbiswas, Fruchart, Vaikuntanathan and Vitelli2019), and show that a bulk-interface correspondence is satisfied at the equator.
This example of the linearized rotating shallow water model constitutes an exactly solvable and physically relevant example well suited to clarifying ongoing issues related to bulk-interface correspondence in continuous media (Bal Reference Bal2018; Faure Reference Faure2019), in a case where the interface is sharp. It also offers a novel interpretation of the equatorial waves as two edge states propagating along a sharp equator separating two flat hemispheres that are topologically distinct, with an explicit computation of the spectrum that is complementary to the beta plane case considered by Matsuno (Reference Matsuno1966), and to the interpretation in terms of transition modes by Iga (Reference Iga1995), who used arguments based on the conservation of zeros in eigenfunctions.
Finally, this case of a sharp equator for the Coriolis parameter provides a first step towards an understanding of the more complicated case of a sharp boundary for the flow domain, that includes the case of coastal Kelvin waves. In the presence of a boundary, Iga (Reference Iga1995) found that coastal Kelvin waves can be removed from the spectrum just by changing the boundary conditions along the wall. This seems to contradict the expected topological robustness of the boundary modes. It is then natural to ask whether this conclusion is robust to the addition of odd viscosity (Souslov et al. Reference Souslov, Dasbiswas, Fruchart, Vaikuntanathan and Vitelli2019). We will indeed confirm Iga’s result in that case, which raises an apparent paradox on the applicability of bulk-edge correspondence in fluids.
2 Linearized rotating shallow water equations with odd viscosity
We consider the rotating shallow water equations linearized around a state of rest in the Cartesian geometry, with an additional odd-viscosity term of amplitude
$\unicode[STIX]{x1D716}$
(Avron Reference Avron1998)
$$\begin{eqnarray}\displaystyle & \unicode[STIX]{x2202}_{t}\unicode[STIX]{x1D702}=-\unicode[STIX]{x2202}_{x}u-\unicode[STIX]{x2202}_{y}v & \displaystyle\end{eqnarray}$$
$$\begin{eqnarray}\displaystyle & \unicode[STIX]{x2202}_{t}u=-\unicode[STIX]{x2202}_{x}\unicode[STIX]{x1D702}+(f+\unicode[STIX]{x1D716}\unicode[STIX]{x1D6FB}^{2})v & \displaystyle\end{eqnarray}$$
$$\begin{eqnarray}\displaystyle & \unicode[STIX]{x2202}_{t}v=-\unicode[STIX]{x2202}_{y}\unicode[STIX]{x1D702}-(f+\unicode[STIX]{x1D716}\unicode[STIX]{x1D6FB}^{2})u, & \displaystyle\end{eqnarray}$$
where
$(u,v)$
are the two (depth-independent) velocity components of the flow in the plane
$(x,y)$
,
$\unicode[STIX]{x1D702}$
the interface elevation relative to the mean depth
$H$
, and
$f$
the Coriolis parameter, which may depend on spatial coordinates. The horizontal variations of interface height correspond to pressure gradients, as the thin layer of fluid satisfies hydrostatic balance. Time units have been chosen such that the shallow water phase speed is
$\sqrt{gH}=1$
, with
$g$
being the standard acceleration due to gravity.
It was realized two decades ago that a two-dimensional system with broken time-reversal symmetry at a microscopic level must include an ‘odd-viscosity’ term Avron (Reference Avron1998). This term has the same effect on the flow evolution as the Coriolis force, except that its strength depends on the wavenumber. It may be thought of as the first-order correction (in
$k$
) to the inviscid dynamics that (i) breaks time-reversal symmetry, (ii) preserves isotropy and (iii) does not dissipate energy. The existence of odd-viscosity terms is prevented by Onsager reciprocity relations only when the microscopic dynamics is time reversible. This is the case for most classical fluids. However, rotating flows are not time-reversal symmetric. Thus, if subgrid-scales dynamics of a rotating systems are modelled by analogy with molecular effects, an odd-viscosity term must be included. Such odd-viscosity terms are actually reminiscent of skew diffusion operators, which have already proved very useful in modelling the effect of baroclinic instability in coarse-resolution ocean models (Vallis Reference Vallis2017).
Here, the odd-viscosity term will be essential to assign a topological invariant to the flow model when the Coriolis parameter
$f$
is prescribed, by regularizing pathological features in bundles of eigenmodes at large wavenumbers. The addition of such terms was also proposed recently by Wiegmann (Reference Wiegmann2013) to describe a fluid of point vortices, and by Banerjee et al. (Reference Banerjee, Souslov, Abanov and Vitelli2017) and Souslov et al. (Reference Souslov, Dasbiswas, Fruchart, Vaikuntanathan and Vitelli2019) for a flow model similar to shallow water equations, motivated by plasma and active matter applications. In the latter work, regularization is also discussed. Moreover, this regularization procedure is well known in a condensed matter context, and probably goes back to Volovik (Reference Volovik1988) for the Dirac Hamiltonian in two dimensions (see also Bal (Reference Bal2018) for a more general description of this problem). In our case we consider
$\unicode[STIX]{x1D716}$
as a (arbitrarily small) constant, and will check that known results are recovered in the limit
$\unicode[STIX]{x1D716}\rightarrow 0$
. This contrasts with the effect of usual viscosity in three-dimensional turbulence, with the occurrence of anomalous dissipation in the limit of weak viscosity.
Another (complementary) way of removing singularities at large wavenumbers is to consider lattice models, such as the discrete models used for numerical simulations (Delplace et al.
Reference Delplace, Marston and Venaille2017). In that case, the admissible wavenumbers are defined on a torus (called the Brillouin zone in condensed matter), and it is possible to compute the Chern number for the bundle of discretized
$f$
-plane shallow water eigenmodes parameterized on this torus. We leave the study of such discrete models to future work, here putting emphasis on continuous models.
3 Bulk waves in the
$f$
-plane
The bulk problem is defined by the case where the flow takes place in a horizontal unbounded plane
$(x,y)$
with a given
$f$
. This is the standard
$f$
-plane approximation. We look at normal modes of the form
$(\unicode[STIX]{x1D702},u,v)=\text{e}^{\text{i}(\unicode[STIX]{x1D714}t-k_{x}x-k_{y}y)}(\hat{\unicode[STIX]{x1D702}},\hat{u} ,\hat{v})$
. The previous system becomes
$$\begin{eqnarray}\unicode[STIX]{x1D714}\left(\begin{array}{@{}c@{}}\hat{\unicode[STIX]{x1D702}}\\ \hat{u} \\ \hat{v}\end{array}\right)=\left(\begin{array}{@{}ccc@{}}0 & k_{x} & k_{y}\\ k_{x} & 0 & -\text{i}(f-\unicode[STIX]{x1D716}k^{2})\\ k_{y} & \text{i}(f-\unicode[STIX]{x1D716}k^{2}) & 0\end{array}\right)\left(\begin{array}{@{}c@{}}\hat{\unicode[STIX]{x1D702}}\\ \hat{u} \\ \hat{v}\end{array}\right),\end{eqnarray}$$
where
$k^{2}:=k_{x}^{2}+k_{y}^{2}$
. The band eigenvalues are (see figure 1
a)
$$\begin{eqnarray}\unicode[STIX]{x1D714}_{\pm }(k_{x},k_{y})=\pm \sqrt{k^{2}+(f-\unicode[STIX]{x1D716}k^{2})^{2}},\quad \unicode[STIX]{x1D714}_{0}(k_{x},k_{y})=0.\end{eqnarray}$$
The middle band
$\unicode[STIX]{x1D714}_{0}$
is flat, and corresponds to odd-geostrophic modes (pressure terms are balanced by Coriolis and odd-viscosity terms). The upper band is an interpolation between two parabolas:
$\unicode[STIX]{x1D714}_{+}\sim |\,f|+((1-2f\unicode[STIX]{x1D716})/2|\,f|)k^{2}$
when
$k\rightarrow 0$
and
$\unicode[STIX]{x1D714}_{+}\sim |\unicode[STIX]{x1D716}|k^{2}$
when
$k\rightarrow \infty$
, and similarly for the lower band. Those modes correspond to odd-Poincaré (or odd-inertia gravity) waves. Importantly, those bands are separated by a range of forbidden frequencies as long as
$f\neq 0$
: the system is gapped. As we shall see it is possible to define a topological invariant, the first Chern number, for each band only in the presence of odd viscosity:
$\unicode[STIX]{x1D716}\neq 0$
.
(a) Dispersion relation of the bulk bands
$(\unicode[STIX]{x1D714}_{-},\unicode[STIX]{x1D714}_{0},\unicode[STIX]{x1D714}_{+})$
, (b) second component
$\hat{u} _{+}\in \mathbb{C}$
of the positive-frequency inertia-gravity eigenmodes
$\unicode[STIX]{x1D6F9}_{+}$
in the
$(k_{x},k_{y})$
-plane, for
$\unicode[STIX]{x1D716}=0.2$
and (c) for
$\unicode[STIX]{x1D716}=0$
. At large
$k$
the components
$\hat{u} _{+}$
are in both cases multivalued. In the case
$\unicode[STIX]{x1D716}\neq 0$
the amplitude of
$\hat{u} _{+}$
is constant, so it can become single-valued up to a phase multiplication, allowing for a compactification at
$k\rightarrow \infty$
. In the case
$\unicode[STIX]{x1D716}=0$
the amplitude at large
$k$
is non-uniform, so compactification is impossible.
3.1 Topology of the upper band
Each eigenmode of (3.1) is defined up to a phase that one can arbitrarily choose locally (gauge freedom), namely for each point of the plane
$(k_{x},k_{y})$
. When such modes are defined on a closed surface, the existence of a smooth phase everywhere is not guaranteed. In that case, it is always possible to remove a phase singularity at a given point by a suitable gauge choice, but this singularity has to appear somewhere else over the closed surface. The impossibility to define globally a continuous phase is a topological property of the mode, captured by an integer-valued index – the first Chern number. In our case, the difficulty in studying these phase singularities originates from the fact that the parameter space is not a closed surface, but the plane
$(k_{x},k_{y})$
. Defining a meaningful topological property for the modes thus requires a compactification of the problem. Even though the compactification of a plane to a (Riemann) sphere is a standard procedure, it is not guaranteed that the eigenmodes defined on the plane can be smoothly mapped on the sphere. In particular, we will show that another kind of singularity of the eigenmodes, appearing at infinity in the plane, may prevent the compactification, and thus the definition of the topological property. However, this issue can be overcome thanks to the odd-viscosity term.
At this step, it is judicious to switch to cylindrical coordinates
$k_{x}=k\cos \unicode[STIX]{x1D719}$
,
$k_{y}=k\sin \unicode[STIX]{x1D719}$
. In particular the bands are
$\unicode[STIX]{x1D719}$
-invariant. A normalized (hence non-vanishing) family of eigenvectors associated with
$\unicode[STIX]{x1D714}_{+}$
is
$$\begin{eqnarray}\unicode[STIX]{x1D6F9}_{+}(k,\unicode[STIX]{x1D719})={\displaystyle \frac{1}{\sqrt{2}}}\left(\begin{array}{@{}c@{}}k/\unicode[STIX]{x1D714}_{+}(k)\\ \cos \unicode[STIX]{x1D719}-\text{i}\sin \unicode[STIX]{x1D719}(f-\unicode[STIX]{x1D716}k^{2})/\unicode[STIX]{x1D714}_{+}(k)\\ \sin \unicode[STIX]{x1D719}+\text{i}\cos \unicode[STIX]{x1D719}(f-\unicode[STIX]{x1D716}k^{2})/\unicode[STIX]{x1D714}_{+}(k)\end{array}\right)\end{eqnarray}$$
up to a phase that can be chosen arbitrarily (gauge freedom). These eigenmodes are regular (single-valued) on the punctured plane
$(k_{x},k_{y})\in \mathbb{R}^{2}\setminus \{0\}$
. We now ask the following questions: (1) can we extend the regularity property for
$k\rightarrow 0$
and
$k\rightarrow \infty$
? (2) If yes, can we do it simultaneously? To answer these questions, first notice that
$$\begin{eqnarray}\lim _{k\rightarrow 0}\unicode[STIX]{x1D6F9}_{+}(k,\unicode[STIX]{x1D719})={\displaystyle \frac{1}{\sqrt{2}}}\left(\begin{array}{@{}c@{}}0\\ \cos \unicode[STIX]{x1D719}-\text{i}\,\text{sign}(f)\sin \unicode[STIX]{x1D719}\\ \sin \unicode[STIX]{x1D719}+\text{i}\,\text{sign}(f)\cos \unicode[STIX]{x1D719}\end{array}\right)=\text{e}^{-\text{i}\,\text{sign}(f)\unicode[STIX]{x1D719}}{\displaystyle \frac{1}{\sqrt{2}}}\left(\begin{array}{@{}c@{}}0\\ 1\\ \text{sign}(f)\text{i}\end{array}\right),\end{eqnarray}$$
where
$\text{sign}(f)=f/|\,f|$
is the sign of
$f$
.
$\unicode[STIX]{x1D6F9}_{+}$
is then multivalued, or singular, at 0. But this singularity can be removed using gauge freedom. We define
$\unicode[STIX]{x1D6F9}_{+}^{A}(k,\unicode[STIX]{x1D719}):=\text{e}^{\text{i}\,\text{sign}(f)\unicode[STIX]{x1D719}}\unicode[STIX]{x1D6F9}_{+}(k,\unicode[STIX]{x1D719})$
, implying that
$\unicode[STIX]{x1D6F9}_{+}^{A}$
is single-valued, or regular on
$\mathbb{R}^{2}$
. The problem has been compactified at 0. Similarly,
$$\begin{eqnarray}\lim _{k\rightarrow \infty }\unicode[STIX]{x1D6F9}_{+}(k,\unicode[STIX]{x1D719})={\displaystyle \frac{1}{\sqrt{2}}}\left(\begin{array}{@{}c@{}}0\\ \cos \unicode[STIX]{x1D719}+\text{i}\,\text{sign}(\unicode[STIX]{x1D716})\sin \unicode[STIX]{x1D719}\\ \sin \unicode[STIX]{x1D719}-\text{i}\,\text{sign}(\unicode[STIX]{x1D716})\cos \unicode[STIX]{x1D719}\end{array}\right)=\text{e}^{\text{i}\,\text{sign}(\unicode[STIX]{x1D716})\unicode[STIX]{x1D719}}{\displaystyle \frac{1}{\sqrt{2}}}\left(\begin{array}{@{}c@{}}0\\ 1\\ -\text{sign}(\unicode[STIX]{x1D716})\text{i}\end{array}\right),\end{eqnarray}$$
where
$\text{sign}(\unicode[STIX]{x1D716})=\unicode[STIX]{x1D716}/|\unicode[STIX]{x1D716}|$
. The singularity of
$\unicode[STIX]{x1D6F9}_{+}$
at
$\infty$
can also be cured using gauge freedom. We define
$\unicode[STIX]{x1D6F9}_{+}^{B}(k,\unicode[STIX]{x1D719}):=\text{e}^{-\text{i}\,\text{sign}(\unicode[STIX]{x1D716})\unicode[STIX]{x1D719}}\unicode[STIX]{x1D6F9}_{+}(k,\unicode[STIX]{x1D719})$
, implying
$\lim _{k\rightarrow 0}\unicode[STIX]{x1D6F9}_{+}^{A}(k,\unicode[STIX]{x1D719})=(1/\sqrt{2})(0,1,-\text{sign}(\unicode[STIX]{x1D716})\text{i})$
so that
$\unicode[STIX]{x1D6F9}_{+}^{B}$
is regular on the punctured
$(k_{x},k_{y})$
plane where the origin is removed
$(\mathbb{R}\setminus \{0\})\cup \{\infty \}$
. The problem has been compactified at
$\infty$
, in the sense that all the infinite directions are equivalent so that we can consider them as a limiting single point. Importantly, this is not possible without odd viscosity (
$\unicode[STIX]{x1D716}=0$
). Indeed
$$\begin{eqnarray}\lim _{k\rightarrow \infty }\unicode[STIX]{x1D6F9}_{+}(k,\unicode[STIX]{x1D719})={\displaystyle \frac{1}{\sqrt{2}}}\left(\begin{array}{@{}c@{}}1\\ \cos \unicode[STIX]{x1D719}\\ \sin \unicode[STIX]{x1D719}\end{array}\right)\end{eqnarray}$$
so that the singularity is impossible to remove by the gauge freedom: the problem is not compactifiable at
$\infty$
. This is illustrated in figure 1.
So far the answer to question 1 is positive, but we defined two families
$\unicode[STIX]{x1D6F9}_{+}^{A}$
and
$\unicode[STIX]{x1D6F9}_{+}^{B}$
that differ by the choice of a phase, illustrated in figure 2. The answer to question 2 is captured by the Chern number. The latter is well defined as a topological invariant for compact manifolds only. Here we identify
$\mathbb{R}^{2}\cup \{\infty \}$
with the two-sphere
$S^{2}$
(see figure 2), for example through the stereographic projection, although we do not need any explicit transformation. The computation of the topological Chern number is now very standard through the integral of the geometrical Berry curvature (Nakahara Reference Nakahara2003). The Berry curvature is the curl of a Berry connection, which allows one to compare two adjacent normalized eigenmodes in parameter space. For
$\unicode[STIX]{x1D6FC}=A,B$
the Berry connection is
$\boldsymbol{A}_{+}^{\unicode[STIX]{x1D6FC}}=-\text{i}\langle \unicode[STIX]{x1D713}_{+}^{\unicode[STIX]{x1D6FC}},\unicode[STIX]{x1D735}\unicode[STIX]{x1D713}_{+}^{\unicode[STIX]{x1D6FC}}\rangle$
and the Berry curvature is
$\boldsymbol{B}_{+}=\unicode[STIX]{x1D735}\times \boldsymbol{A}_{+}^{\unicode[STIX]{x1D6FC}}$
, independent of
$\unicode[STIX]{x1D6FC}$
. The connection depends on the phase of the eigenmodes, but the curvature is gauge-independent.
Normalized second component
$\hat{u} _{+}/|\hat{u} _{+}|\in \mathbb{C}$
of positive-frequency inertia-gravity eigenmodes
$\unicode[STIX]{x1D6F9}+$
for
$f=1$
and
$\unicode[STIX]{x1D716}=0.2$
in the
$(k_{x},k_{y})$
plane, or equivalently on the sphere where the parallels correspond to circles of fixed radius
$k$
. (a)
$\unicode[STIX]{x1D6F9}_{+}$
is singular at 0 and at
$\infty$
, but each singularity can be cured up to a phase: (b) at 0 with
$\unicode[STIX]{x1D6F9}_{+}^{A}$
or (c) at
$\infty$
with
$\unicode[STIX]{x1D6F9}_{+}^{B}$
. The problem is compactified in the sense that each singularity can be considered as a single point since there exists an eigenvector
$\unicode[STIX]{x1D6F9}$
that is regular (for all the components) in its neighbourhood. A non-vanishing Chern number
$C_{+}$
then captures the impossibility of finding a
$\unicode[STIX]{x1D6F9}$
that is regular everywhere. In the case where
$\unicode[STIX]{x1D716}=0$
, neither
$\unicode[STIX]{x1D6F9}_{+}^{B}$
nor (c) exist, so the Chern number is not even well defined.
On
$\mathbb{R}^{2}\setminus \{0\}$
one has
$\unicode[STIX]{x1D713}_{+}^{B}(k,\unicode[STIX]{x1D711})=\text{e}^{-\text{i}(\text{sign}(f)+\text{sign}(\unicode[STIX]{x1D716}))\unicode[STIX]{x1D719}}\unicode[STIX]{x1D713}_{+}^{A}(k,\unicode[STIX]{x1D711})$
and thus
$\boldsymbol{A}_{+}^{B}=\boldsymbol{A}_{+}^{A}-(\text{sign}(f)+\text{sign}(\unicode[STIX]{x1D716}))\unicode[STIX]{x1D735}\unicode[STIX]{x1D719}$
. The Chern number is the flux of the Berry curvature through the whole (compactified) plane, namely
$$\begin{eqnarray}C_{+}={\displaystyle \frac{1}{2\unicode[STIX]{x03C0}}}\int _{0}^{\infty }\text{d}k\int _{0}^{2\unicode[STIX]{x03C0}}k\,\text{d}\unicode[STIX]{x1D719}\,\boldsymbol{B}_{+}\boldsymbol{\cdot }\boldsymbol{e}_{z}.\end{eqnarray}$$
Splitting
$\mathbb{R}^{2}\cup \{\infty \}=D_{{<}}\cup D_{{>}}$
, respectively the disk
$\{k\leqslant 1\}$
and its complementary, we apply Stokes theorem on each part and are left with a contribution at the border
$k=1$
. Explicitly
$$\begin{eqnarray}C_{+}={\displaystyle \frac{1}{2\unicode[STIX]{x03C0}}}\int _{0}^{2\unicode[STIX]{x03C0}}\text{d}\unicode[STIX]{x1D719}(\boldsymbol{A}_{+}^{A}-\boldsymbol{A}_{+}^{B})\boldsymbol{\cdot }\boldsymbol{e}_{\unicode[STIX]{x1D719}}=\text{sign}(f)+\text{sign}(\unicode[STIX]{x1D716}).\end{eqnarray}$$
Souslov et al. (Reference Souslov, Dasbiswas, Fruchart, Vaikuntanathan and Vitelli2019) reached a similar conclusion for the bulk topology when computing the integral of the Berry curvature over the whole plane
$(k_{x},k_{y})$
, and this result can be recovered for the Dirac Hamiltonian in arbitrary dimensions (Bal Reference Bal2018). Note that if we start with
$\unicode[STIX]{x1D716}=0$
the latter derivation leads to
$C_{+}=\text{sign}(f)$
, but this is no longer a Chern number, even if the Berry curvature integrated on the non-compact manifold
$\mathbb{R}^{2}$
is finite.
3.2 Topology of the lower and middle bands
The Chern number of the lower band is by the symmetry of the system
$\unicode[STIX]{x1D6F9}_{-}(k_{x},k_{y},f,\unicode[STIX]{x1D716})=\unicode[STIX]{x1D6F9}_{+}(-k_{x},-k_{y},-f,-\unicode[STIX]{x1D716})$
, leading immediately to
$C_{-}=-(\text{sign}(f)+\text{sign}(\unicode[STIX]{x1D716}))=-C_{+}$
. Finally, the normalized family of eigenvectors associated with
$\unicode[STIX]{x1D714}=0$
is
$$\begin{eqnarray}\unicode[STIX]{x1D713}_{0}(k_{x},k_{y})={\displaystyle \frac{1}{\sqrt{k^{2}+(f-\unicode[STIX]{x1D716}k^{2})^{2}}}}\left(\begin{array}{@{}c@{}}f-\unicode[STIX]{x1D716}k^{2}\\ \text{i}k\sin \unicode[STIX]{x1D719}\\ -\text{i}k\cos \unicode[STIX]{x1D719}\end{array}\right).\end{eqnarray}$$
In particular,
$\lim _{k\rightarrow 0}\unicode[STIX]{x1D6F9}_{0}=(\text{sign}(f),0,0)$
and
$\lim _{k\rightarrow \infty }\unicode[STIX]{x1D6F9}_{0}=(-\text{sign}(\unicode[STIX]{x1D716}),0,0)$
, so
$\unicode[STIX]{x1D6F9}_{0}$
is regular on the whole compactified plane
$\mathbb{R}\cup \{\infty \}$
: it is a global continuous section on a compact manifold. Thus
$C_{0}=0$
. Again, notice that when
$\unicode[STIX]{x1D716}=0$
one has
$\lim _{k\rightarrow \infty }\unicode[STIX]{x1D6F9}_{0}=(0,\text{i}\sin \unicode[STIX]{x1D719},-\text{i}\cos \unicode[STIX]{x1D719})$
, so the problem is not compactifiable at
$\infty$
.
3.3 Summarizing
If
$f$
and
$\unicode[STIX]{x1D716}$
have the same sign
$s=\pm$
, then the Chern numbers for the three wavebands are
$C_{+}=-C_{-}=2s$
and
$C_{0}=0$
. If
$f$
and
$\unicode[STIX]{x1D716}$
have opposite sign, then
$C_{+}=C_{-}=C_{0}=0$
. We now compute the wave spectrum when two hemispheres are glued together, with a given value of odd viscosity
$\unicode[STIX]{x1D716}$
. According to the bulk-interface correspondence, and whatever the sign of
$\unicode[STIX]{x1D716}$
, we expect that two unidirectional edge states should fill each frequency gap between the flat band of geostrophic modes and the inertia-gravity wavebands.
4 Interface: an equator between two flat hemispheres
We consider a sharp interface at
$y=0$
between two hemispheres that can now be interpreted as two distinct topological phases (see figure 3
a). On the upper half-plane
$f$
is constant and positive, while on the lower half-plane
$-f<0$
. The parameter
$\unicode[STIX]{x1D716}$
is constant and positive on the whole plane. To simplify the computations below we assume
$f\unicode[STIX]{x1D716}<1/4$
.
(a) Sharp interface between two hemispheres interpreted as two distinct topological phases. (b) Dispersion relation of the modes confined at the interface (a Yanai-like wave in orange and a Kelvin-like wave in red) for
$f=1$
and
$\unicode[STIX]{x1D716}=0.2$
: in each gap one has two topological solutions, in agreement with the bulk-edge correspondence. The black curve is the limit of the bulk inertia-gravity wavebands area. The vertical dashed lines correspond to
$|k_{x}|=k_{0}$
where the Kelvin mode merges into the bulk by compactification. This cutoff tends to infinity when the odd viscosity
$\unicode[STIX]{x1D716}$
tends to zero.
The translation invariance is preserved in the longitudinal direction, and we look for normal modes
$(\unicode[STIX]{x1D702},u,v)=\text{e}^{\text{i}(\unicode[STIX]{x1D714}t-k_{x}x)}(\tilde{\unicode[STIX]{x1D702}},\tilde{u} ,\tilde{v})$
, the latter vector being a function of variable
$y$
and parameters
$k_{x}$
and
$\unicode[STIX]{x1D714}$
. Dropping the tilde, it is governed by
$$\begin{eqnarray}\displaystyle & \text{i}\unicode[STIX]{x1D714}\unicode[STIX]{x1D702}=\text{i}k_{x}u-\unicode[STIX]{x2202}_{y}v & \displaystyle\end{eqnarray}$$
$$\begin{eqnarray}\displaystyle & \text{i}\unicode[STIX]{x1D714}u=\text{i}k_{x}\unicode[STIX]{x1D702}+(\pm f-\unicode[STIX]{x1D716}k_{x}^{2})v+\unicode[STIX]{x1D716}\unicode[STIX]{x2202}_{yy}v & \displaystyle\end{eqnarray}$$
$$\begin{eqnarray}\displaystyle & \text{i}\unicode[STIX]{x1D714}v=-\unicode[STIX]{x2202}_{y}\unicode[STIX]{x1D702}-(\pm f-\unicode[STIX]{x1D716}k_{x}^{2})u-\unicode[STIX]{x1D716}\unicode[STIX]{x2202}_{yy}u. & \displaystyle\end{eqnarray}$$
There is a redundancy in the system as
$\unicode[STIX]{x1D702}$
is completely determined by
$u$
and
$v$
from (4.1). This can be eliminated, leading to a two-dimensional problem of order two:
$$\begin{eqnarray}\displaystyle & \displaystyle \left(\unicode[STIX]{x1D716}\unicode[STIX]{x2202}_{yy}-{\displaystyle \frac{k_{x}}{\unicode[STIX]{x1D714}}}\unicode[STIX]{x2202}_{y}+(\pm f-\unicode[STIX]{x1D716}k_{x}^{2})\right)v={\displaystyle \frac{\text{i}}{\unicode[STIX]{x1D714}}}(\unicode[STIX]{x1D714}^{2}-k_{x}^{2})u & \displaystyle\end{eqnarray}$$
$$\begin{eqnarray}\displaystyle & \displaystyle \left(\unicode[STIX]{x1D716}\unicode[STIX]{x2202}_{yy}+{\displaystyle \frac{k_{x}}{\unicode[STIX]{x1D714}}}\unicode[STIX]{x2202}_{y}+(\pm f-\unicode[STIX]{x1D716}k_{x}^{2})\right)u=-{\displaystyle \frac{\text{i}}{\unicode[STIX]{x1D714}}}(\unicode[STIX]{x2202}_{yy}+\unicode[STIX]{x1D714}^{2})v. & \displaystyle\end{eqnarray}$$
At the interface we impose the continuity of
$u$
and
$v$
as well as their first derivative
$$\begin{eqnarray}u|_{y=0^{-}}=u|_{y=0^{+}},\quad v|_{y=0^{-}}=v|_{y=0^{-}},\quad \unicode[STIX]{x2202}_{y}u|_{y=0^{-}}=\unicode[STIX]{x2202}_{y}u|_{y=0^{+}},\quad \unicode[STIX]{x2202}_{y}v|_{y=0^{-}}=\unicode[STIX]{x2202}_{y}v|_{y=0^{-}},\end{eqnarray}$$
which implies the same for
$\unicode[STIX]{x1D702}$
. We look for the solutions that are localized near the interface, in the sense that they vanish away from it, when
$y\rightarrow \pm \infty$
. We first solve the problem in each half-plane by decomposing
$$\begin{eqnarray}u(y)=\left\{\begin{array}{@{}ll@{}}u_{\uparrow }(y), & y>0,\\ u_{\downarrow }(y), & y<0,\end{array}\right.\end{eqnarray}$$
and similarly for
$v$
, then glue the solutions at the interface through (4.6) and count the remaining degrees of freedom, leading to a number of interface modes. This number is algebraic, its sign being determined by
$\unicode[STIX]{x2202}\unicode[STIX]{x1D714}/\unicode[STIX]{x2202}k_{x}$
. The result is given in figure 3 with
$n=2$
modes in each gap. We now explain how to reconstruct it.
4.1 The compactified Kelvin wave
It is easy to see that
$u\equiv 0$
in the whole plane implies
$v\equiv 0$
. In this section we only assume
$v\equiv 0$
and look for non-trivial
$u$
. By (4.4) we infer
$k_{x}^{2}=\unicode[STIX]{x1D714}^{2}$
, and
$u$
is governed by (4.5), a homogeneous equation of order 2. The solution depends on the relative signs of
$k_{x}$
and
$\unicode[STIX]{x1D714}$
.
Case 1.
$k_{x}=\unicode[STIX]{x1D714}$
. On the upper half-plane the solution is of the form
$$\begin{eqnarray}u_{\uparrow }(y)=A_{\uparrow }\text{e}^{q_{\uparrow +}y}+B_{\uparrow }\text{e}^{q_{\uparrow -}y},\quad \text{where}~q_{\uparrow \pm }=-{\displaystyle \frac{1}{2\unicode[STIX]{x1D716}}}\left({\displaystyle \frac{k_{x}}{\unicode[STIX]{x1D714}}}\pm \sqrt{1+4\unicode[STIX]{x1D716}(\unicode[STIX]{x1D716}k_{x}^{2}-f)}\right),\end{eqnarray}$$
which is always well defined as long as
$f\unicode[STIX]{x1D716}\leqslant 1/4$
. Notice that
$q_{\uparrow \pm }$
(as well as
$A_{\uparrow }$
and
$B_{\uparrow }$
) depends on
$k_{x}$
and
$\unicode[STIX]{x1D714}$
. For
$u_{\uparrow }$
to vanish at
$y\rightarrow \infty$
we have either
$q_{\uparrow \pm }<0$
or
$A_{\uparrow }/B_{\uparrow }=0$
. Here
$q_{\uparrow +}<0$
for all
$k_{x}$
and
$q_{\uparrow -}<0$
only for
$|k_{x}|<k_{0}:=\sqrt{f/\unicode[STIX]{x1D716}}$
, so that
$$\begin{eqnarray}u_{\uparrow }(y)=\left\{\begin{array}{@{}ll@{}}A_{\uparrow }\text{e}^{q_{\uparrow +}y}+B_{\uparrow }\text{e}^{q_{\uparrow -}y}, & |k_{x}|<k_{0},\\ A_{\uparrow }\text{e}^{q_{\uparrow +}y}, & |k_{x}|\geqslant k_{0}.\end{array}\right.\end{eqnarray}$$
In the lower half-plane, one has similarly
$$\begin{eqnarray}q_{\downarrow \pm }=-{\displaystyle \frac{1}{2\unicode[STIX]{x1D716}}}\left({\displaystyle \frac{k_{x}}{\unicode[STIX]{x1D714}}}\pm \sqrt{1+4\unicode[STIX]{x1D716}(\unicode[STIX]{x1D716}k_{x}^{2}+f)}\right),\end{eqnarray}$$
but this time we select the positive roots, so that
$u_{\downarrow }$
vanishes when
$y\rightarrow -\infty$
. Here
$q_{\downarrow +}<0$
and
$q_{\downarrow -}>0$
for all
$k_{x}$
, so that
$u_{\downarrow }(y)=B_{\downarrow }\text{e}^{q_{\downarrow -}y}$
. From the interface condition (4.6) we infer two relations between
$A_{\uparrow }$
,
$B_{\uparrow }$
and
$B_{\downarrow }$
for
$|k_{x}|<k_{0}$
whereas
$A_{\uparrow }=B_{\downarrow }=0$
for
$|k_{x}|\geqslant k_{0}$
. More precisely
$$\begin{eqnarray}u(y)=\left\{\begin{array}{@{}ll@{}}A_{\uparrow }\left(\text{e}^{q_{\uparrow +}y}-{\displaystyle \frac{q_{\uparrow +}-q_{\downarrow -}}{q_{\uparrow -}-q_{\downarrow -}}}\text{e}^{q_{\uparrow -}y}\right) & y>0,\,|k_{x}|<k_{0}\\ A_{\uparrow }{\displaystyle \frac{q_{\uparrow -}-q_{\uparrow +}}{q_{\uparrow -}-q_{\downarrow -}}}\text{e}^{q_{\downarrow -}y} & y<0,\,|k_{x}|<k_{0}\\ 0 & |k_{x}|\geqslant k_{0}.\end{array}\right.\end{eqnarray}$$
We are left with one degree of freedom
$A_{\uparrow }$
, a positive dispersion relation
$\unicode[STIX]{x1D714}=k_{x}$
that is moreover compactified: the solution exists only for
$|k_{x}|<k_{0}$
. Note that
$k_{0}\rightarrow \infty$
as
$\unicode[STIX]{x1D716}\rightarrow 0$
. Moreover, in that limit
$q_{\uparrow +}\sim -1/\unicode[STIX]{x1D716}$
,
$q_{\uparrow -}\rightarrow -f$
and
$q_{\downarrow -}\rightarrow f$
. For
$y>0$
,
$\text{e}^{q_{+}y}\rightarrow 0$
so that (4.11) becomes, after renormalizing
$A_{\uparrow }=\unicode[STIX]{x1D716}\,\tilde{A}_{\uparrow }$
,
$$\begin{eqnarray}u(y)\mathop{\longrightarrow }_{\unicode[STIX]{x1D716}\rightarrow 0}\tilde{A}_{\uparrow }(-2f)^{-1}\,\text{e}^{-f|y|},\quad y\in \mathbb{R},\,k_{x}\in \mathbb{R}.\end{eqnarray}$$
It is remarkable that the limit
$\unicode[STIX]{x1D716}\rightarrow 0$
coincides with the classical (non-compactified) Kelvin wave solution obtained when
$\unicode[STIX]{x1D716}=0$
, in which case the order of the partial differential equation is lowered.
Case 2.
$k_{x}=-\unicode[STIX]{x1D714}$
. In this case
$q_{\uparrow /\downarrow \pm }$
have the same expression but their signs change. A similar inspection leads to
$u_{\uparrow }(y)=A_{\uparrow }\text{e}^{q_{\uparrow +}y}$
for
$|k_{x}|\geqslant k_{0}$
and vanishes otherwise, and
$u_{\downarrow }(y)=B_{\downarrow }\text{e}^{q_{\downarrow -}y}$
for
$k_{x}\in \mathbb{R}$
. The gluing condition implies
$A_{\uparrow }=B_{\downarrow }=0$
, so that
$u\equiv 0$
.
4.2 The compactified Yanai wave
In this section we assume
$k_{x}^{2}\neq \unicode[STIX]{x1D714}^{2}$
and
$v\neq 0$
. In that case
$u$
is entirely fixed by
$v$
through (4.4), and one can moreover combine (4.4) and (4.5) to get a fourth-order homogeneous equation for
$v$
. In the upper half-plane it reads
$$\begin{eqnarray}(\unicode[STIX]{x1D716}^{2}\unicode[STIX]{x2202}_{y}^{(4)}+(2\unicode[STIX]{x1D716}(f-\unicode[STIX]{x1D716}k_{x}^{2})-1)\unicode[STIX]{x2202}_{y}^{(2)}+(f-\unicode[STIX]{x1D716}k_{x}^{2})^{2}-(\unicode[STIX]{x1D714}^{2}-k_{x}^{2}))v=0.\end{eqnarray}$$
The corresponding algebraic equation always admits real solutions as long as
$f\unicode[STIX]{x1D716}\leqslant 1/4$
, given by
$s^{2}=S_{\pm }$
with
$$\begin{eqnarray}S_{\pm }={\displaystyle \frac{1}{2\unicode[STIX]{x1D716}^{2}}}\left(1+2\unicode[STIX]{x1D716}(\unicode[STIX]{x1D716}k_{x}^{2}-f)\pm \sqrt{1+4\unicode[STIX]{x1D716}(\unicode[STIX]{x1D716}\unicode[STIX]{x1D714}^{2}-f)}\right).\end{eqnarray}$$
In order to get a non-trivial mode at the interface, we need both
$S_{+}>0$
and
$S_{-}>0$
. This implies
$$\begin{eqnarray}\unicode[STIX]{x1D6E5}(k_{x},\unicode[STIX]{x1D714}):=k_{x}^{2}-\unicode[STIX]{x1D714}^{2}+(f-\unicode[STIX]{x1D708}k_{x}^{2})^{2}>0.\end{eqnarray}$$
In the region where
$\unicode[STIX]{x1D6E5}(k_{x},\unicode[STIX]{x1D714})\leqslant 0$
, oscillatory solutions exist: they are the normal modes from the bulk which are still allowed in the interface setting. This region corresponds to the projection of the bulk bands
$\unicode[STIX]{x1D714}_{\pm }(k_{x},k_{y})$
from (3.2) for all values of
$k_{y}\in \mathbb{R}$
. Thus (4.15) delimits the gapped region in the interface setting. In this region one has four real solutions to (4.13)
$$\begin{eqnarray}s_{\uparrow 1}=\sqrt{S_{+}},\quad s_{\uparrow 2}=\sqrt{S_{-}},\quad s_{\uparrow 3}=-\sqrt{S_{+}},\quad s_{\uparrow 4}=-\sqrt{S_{-}}\end{eqnarray}$$
and similarly for
$s_{\downarrow i}$
,
$i=1,\ldots ,4$
, where we replace
$f$
by
$-f$
in the expression of (4.14). By construction
$s_{\uparrow /\downarrow 1/2}>0$
and
$s_{\uparrow /\downarrow 3/4}<0$
regardless of
$k_{x},\unicode[STIX]{x1D714}$
or
$f$
. Consequently,
$$\begin{eqnarray}v(y)=\left\{\begin{array}{@{}ll@{}}V_{\uparrow 3}\text{e}^{s_{\uparrow 3}y}+V_{\uparrow 4}\text{e}^{s_{\uparrow 4}y}, & y>0\\ V_{\downarrow 1}\text{e}^{s_{\downarrow 1}y}+V_{\downarrow 2}\text{e}^{s_{\downarrow 2}y}, & y<0.\end{array}\right.\end{eqnarray}$$
Moreover,
$u_{\uparrow }(y)=\unicode[STIX]{x1D706}_{\uparrow 3}V_{\uparrow 3}\text{e}^{s_{\uparrow 3}y}+\unicode[STIX]{x1D706}_{\uparrow 4}V_{\uparrow 4}\text{e}^{s_{\uparrow 4}y}$
and
$u_{\downarrow }(y)=\unicode[STIX]{x1D706}_{\downarrow 1}V_{\downarrow 1}\text{e}^{s_{\downarrow 1}y}+\unicode[STIX]{x1D706}_{\downarrow 2}V_{\downarrow 2}\text{e}^{s_{\downarrow 2}y}$
where
$$\begin{eqnarray}\unicode[STIX]{x1D706}_{\uparrow /\downarrow i}={\displaystyle \frac{\unicode[STIX]{x1D714}}{\text{i}(\unicode[STIX]{x1D714}^{2}-k_{x}^{2})}}\left(\unicode[STIX]{x1D716}s_{\uparrow /\downarrow i}^{2}-{\displaystyle \frac{k_{x}}{\unicode[STIX]{x1D714}}}s_{\uparrow /\downarrow i}\pm f-\unicode[STIX]{x1D716}k_{x}^{2}\right),\end{eqnarray}$$
inferred from (4.4). The four free parameters
$V_{\uparrow 3},V_{\uparrow 4},V_{\downarrow 1}$
and
$V_{\downarrow 2}$
are constrained by the four gluing conditions (4.6), so that a non-trivial solution exists only if
$\det (\unicode[STIX]{x1D648})=0$
, where
$$\begin{eqnarray}\unicode[STIX]{x1D648}=\left(\begin{array}{@{}cccc@{}}1 & 1 & -1 & -1\\ s_{\uparrow 3} & s_{\uparrow 4} & -s_{\downarrow 1} & -s_{\downarrow 2}\\ \unicode[STIX]{x1D706}_{\uparrow 3} & \unicode[STIX]{x1D706}_{\uparrow 4} & -\unicode[STIX]{x1D706}_{\downarrow 1} & -\unicode[STIX]{x1D706}_{\downarrow 2}\\ \unicode[STIX]{x1D706}_{\uparrow 3}s_{\uparrow 3} & \unicode[STIX]{x1D706}_{\uparrow 4}s_{\uparrow 4} & -\unicode[STIX]{x1D706}_{\downarrow 1}s_{\downarrow 1} & -\unicode[STIX]{x1D706}_{\downarrow 2}s_{\downarrow 2}\end{array}\right).\end{eqnarray}$$
There is no simple expression for the dispersion relation
$\unicode[STIX]{x1D714}=f(k_{x})$
such that
$\det (\unicode[STIX]{x1D648})=0$
, however the latter gives an implicit relation between
$k_{x}$
and
$\unicode[STIX]{x1D714}$
, namely
$$\begin{eqnarray}\unicode[STIX]{x1D716}^{2}(s_{\downarrow 1}-s_{\uparrow 3})(s_{\downarrow 2}-s_{\uparrow 3})(s_{\downarrow 1}-s_{\uparrow 4})(s_{\downarrow 2}-s_{\uparrow 4})+2f{\displaystyle \frac{k_{x}}{\unicode[STIX]{x1D714}}}(s_{\downarrow 1}+s_{\downarrow 2}-s_{\uparrow 3}-s_{\uparrow 4})-4f^{2}=0,\end{eqnarray}$$
that can be computed numerically. In this way we obtain the Yanai wave of figure 3: one in each gap. An expansion around
$k_{x}\rightarrow \pm \infty$
shows the behaviour
$\unicode[STIX]{x1D714}\sim \mp (f/(2(1+k_{x}^{2}\unicode[STIX]{x1D716}^{2})))$
, so that the dispersion relation connects (asymptotically) the middle band
$\unicode[STIX]{x1D714}=0$
to the upper band at
$k_{x}=0$
. On the other hand for a finite value of
$k_{x}$
, an expansion around
$\unicode[STIX]{x1D716}\rightarrow 0$
leads to
$\unicode[STIX]{x1D714}\sim \mp f$
: the compactified Yanai waves is an inertial wave in this limit, as in Iga (Reference Iga1995), where
$\unicode[STIX]{x1D716}=0$
. For this mode the rank of
$\unicode[STIX]{x1D648}$
is 3, so that we have only one free parameter among
$V_{\uparrow 3},V_{\uparrow 4},V_{\downarrow 1}$
and
$V_{\downarrow 2}$
: we have an
$n=1$
Yanai mode in each gap.
The last case
$k_{x}^{2}=\unicode[STIX]{x1D714}^{2}$
with
$v\neq 0$
is rather tedious, but it can be checked that no extra mode appears in this case. This is done in the supplementary material, available at https://doi.org/10.1017/jfm.2019.233.
5 Bulk-boundary correspondence and boundary conditions
To emphasize the relevance of the interface studied above we briefly discuss the sharp boundary problem. Figure 4 shows the spectrum of odd-shallow water waves computed numerically with Dedalus (2016) on a channel (or infinite strip) with sharp walls:
$(x,y)\in \mathbb{R}\times [0,L]$
, where
$\unicode[STIX]{x1D716}$
and
$f$
have the same sign. The boundary conditions are either:
$v=0$
,
$u=0$
(no slip) or
$v=0$
,
$\unicode[STIX]{x2202}_{x}u-\unicode[STIX]{x2202}_{y}v=0$
(stress-free).
Dispersion relation of the shallow water model in a cylinder geometry with
$f=1$
,
$g=H=1$
, using a numerical code (Dedalus 2016). (a) No-slip boundary condition with a large odd viscosity (
$\unicode[STIX]{x1D716}=4$
,
$L_{y}=120$
), (b) no-slip boundary condition with weak odd viscosity (
$\unicode[STIX]{x1D716}=0.2$
,
$L_{y}=30$
), (c) stress-free boundary condition with the same parameters as in case (b). The colour code indicates localization of the modes in the
$y$
direction (from 0 at
$y=0$
to 1 at
$y=L_{y}$
).
Like for the interface, we recover the region of the (projected) bulk bands as well as edge modes in the gapped region, this time localized on each wall. In the following we focus on one of them (e.g.
$y=0$
). In figure 4(a) the number of modes crossing a fixed frequency line
$\unicode[STIX]{x1D714}$
in the gap is 2, in agreement with the Chern number of the upper band (Souslov et al.
Reference Souslov, Dasbiswas, Fruchart, Vaikuntanathan and Vitelli2019). However, this number becomes 1 either when
$\unicode[STIX]{x1D714}_{0}$
is too close to the middle band
$\unicode[STIX]{x1D714}_{0}=0$
or when
$\unicode[STIX]{x1D716}$
is smaller: in that case the other mode never crosses the spectral gap window
$0<\unicode[STIX]{x1D714}<f$
(figure 4
b). Moreover, when changing the boundary condition, the total number of modes surprisingly jumps from 2 to 3 in figure 4(c).
Thus it seems that the bulk-edge correspondence is not always satisfied for the sharp boundary problem since the number of edge modes depends on the choice of parameters and boundary conditions. On the other hand, our interface setting with a canonical choice of gluing condition provides a remarkable case where the bulk-interface correspondence is fully satisfied: the number of modes localized at the interface inside the spectral region
$0<\unicode[STIX]{x1D714}<f$
matches with the Chern number of the upper band.
6 Discussion and conclusion
The most representative difference of our study from the original description of equatorial waves by Matsuno (Reference Matsuno1966) is that we consider a profile
$f(y)$
with a sharp gap at the equator
$(y=0)$
and constant in each hemisphere, rather than linear variations of the Coriolis parameter with latitude (
$f=\unicode[STIX]{x1D6FD}y$
). In this celebrated equatorial beta plane configuration, all the modes are trapped along the equator, and described by Hermite functions (Matsuno Reference Matsuno1966). Here, in the case of a sharp interface, only the Yanai-like and Kelvin-like waves are exponentially trapped modes; the other modes are all delocalized. In addition, there is no Rossby wave. These properties are caused by the hypothesis of a constant value of
$f$
in each hemisphere. If we keep the hypothesis of a sharp equator but add back a constant gradient of Coriolis parameter into the problem, topological properties are preserved, but (i) the degeneracy of geostrophic modes is lifted, with the emergence of Rossby waves, and (ii) all the modes are localized at the equator (like for Hermite functions), even if the Kelvin and Yanai waves remain more localized than the others, which can be understood by counting the zeros of eigenfunctions (Iga Reference Iga1995). Using arguments based on the conservation of these zeros, Iga (Reference Iga1995) explained the global shape of equatorial wave spectra computed by Matsuno. In particular, he found that the Yanai wave can be interpreted as an inertial wave. Our study brings a complementary point view, showing the topological origin of these properties as the outcome of gluing two hemispheres with different topological indices. Our analysis also suggests that Yanai waves should be qualified as mixed geostrophic-gravity waves rather than mixed Rossby-gravity waves, as they exist even in the absence of Rossby waves. This result relied on the introduction of a regularization parameter, but we recover the spectrum of the original problem when the regularization coefficient tends to zero: there is no singular limit in the sharp interface case.
To summarize, (i) it has been possible to assign a topological index to each hemisphere through the introduction of odd-viscosity terms, as in Souslov et al. (Reference Souslov, Dasbiswas, Fruchart, Vaikuntanathan and Vitelli2019). (ii) Souslov et al. (Reference Souslov, Dasbiswas, Fruchart, Vaikuntanathan and Vitelli2019) showed a range of parameters and boundary conditions for which the correspondence between the number of unidirectional edge states filling the frequency gap and the bulk Chern number were satisfied. Building on previous work by Iga (Reference Iga1995), we noticed that the number of edge states that transit from one band to another in fact depends on the choice of the boundary condition, even in the presence of odd viscosity. This new observation raises the question of the existence of a bulk-boundary correspondence for fluids in particular, and for continuous media in general. (iii) We have considered the simpler case of a sharp interface separating the two hemispheres. This interface case bypasses the need to discuss boundary conditions, as the natural physical choice is to impose continuity of the fields and their derivatives. We have found in that case explicit analytical solutions exhibiting exactly two unidirectional modes in each gap confined along the equator, accordingly with the bulk-interface correspondence. Remarkably, this simple case provides a solvable example of a bulk-interface correspondence in fluids, in a configuration where the interface in infinitely sharp. This paves the way towards an understanding of the bulk-edge correspondence in continuous media with boundaries. We will explain in a companion paper that the apparent paradox between the number of edge states and the bulk Chern number can indeed be explained, beyond the particular case of the shallow water model.
Acknowledgements
C.T. is grateful to G. M. Graf and H. Jud for many insightful discussions. P.D. and A.V. were partly funded by ANR-18-CE30-0002-01 during this work, and thank L.-A. Couston for help with Dedalus.
Supplementary material
Supplementary material is available at https://doi.org/10.1017/jfm.2019.233.
