## Quantum Hall effect induced by electron-phonon interaction

Andreas Sinner and Klaus Ziegler

Institut f¨r u Physik, Universit¨t a Augsburg, D-86135 Augsburg, Germany

When phonons couple to fermions in 2D semimetals, the interaction may turn the system into an insulator. There are several insulating phases in which the time reversal and the sublattice symmetries are spontaneously broken. Examples are many-body states commensurate to Haldane's staggered flux model or to lattice models with periodically modulated strain. We find that the effective field theories of these phases exhibit characteristic Chern-Simons terms, whose coefficients are related to the topological invariants of the microscopic model. This implies that the corresponding quantized Hall conductivities characterize these insulating states.

## I. INTRODUCTION

The ground state of a many-body Hamiltonian may pass through a transition between phases with very different properties under a change of model parameters. This process is usually accompanied by the symmetry change of ground states on both sides of the transition, the phenomenon known as the spontaneous symmetry breaking. The phase with the lower symmetry is characterized by the emergent macroscopic ordering of elementary building blocks of the system, technically captured by an order parameter. Phenomenological models describing the physics of the symmetry broken phase in terms of order parameters have proven very successful in the description of various thermal and quantum phase transitions. In these models, critical parameters at which phase transitions occur correspond to emergent non-trivial potential minima. They can be approached by variational techniques, which lead to a set of mean-field equations for all order parameters. Sometimes, the ground states of a system is not uniquely determined by the order parameter alone, as they may have additional topological ordering [1-3]. An important observation is that different phases of a degenerated ground state may be described by different topological field theories. Usually in condensed matter physics the emergence of the topological ChernSimons excitations has not been considered in the broader context of symmetry breaking and phase transitions. To some extend, this is the case for superfluid Helium, where 3 He -A phase exhibits also a spontaneous time reversal symmetry breaking among other spontaneously broken symmetries [4-6].

The phase of the wave function can play a crucial role for the properties of the related physical system [7, 8]. In particular, if this phase is associated with a topological invariant through the winding number, it may describe robust macroscopic properties. Examples in condensed matter physics are the quantized Hall conductivity of electrons in 2d [1, 3, 9-13] and in quasi 1d structures [14]. Another example is the quantized transverse conductance of the supercurrent in the 3 He -A-phase of superfluid liquid Helium [4-6], and quantized conductivities of different anomalous Hall effects [15-17]. In general, topological structures in 2+1-dimensional gauge field theories are related to Chern-Simons terms, where the coefficient in front of the latter is related to a topological invariant [18-23].

The idea that dynamical lattice degrees of freedom can be treated as effective gauge fields is not new [24, 25]. Here we study a system of monochromatic phonons minimally coupled to Dirac fermions in 2D. As the electron-phonon coupling strength increases the system should approach a transition into the structurally different phase. Our intention is to identify possible gapped phases, the corresponding quasiparticles and observables. The connection with tight-binding models will help us to link these phases to different lattice distortions, which lead to the formation of spectral gaps with broken time reversal and sublattice symmetries. To emphasize the similarity of phonons and gauge fields we consider the E 1 -optical phonon mode of a honeycomb lattice which represents individual vibrations of the sublattices in opposite directions [26, 27]. The number of modes is therefore naturally restricted to 2 and there is no time-like component of the effective gauge field. This resembles the popular in electrodynamics 'light cone' gauge fixing ˆ n · /vector A = 0, with ˆ denoting a unit vector pointing into the time direction [2]. n It was shown in [28, 29] that in such systems phonon Chern-Simons excitations can be generated.

## II. THE MODEL

The properties of tight-binding Hamiltonians which describe the motion of electrons on the honeycomb lattice are well known and investigated in details in the existing literature [30, 31]. In sublattice representation it reads

$$H _ { \text{el} } = \sum _ { r, r ^ { \prime } } h _ { 1 ; \text{rr} ^ { \prime } } c _ { \text{r} } ^ { \dagger$$

with the Pauli matrices σ µ reflecting the sublattice degree of freedom. The hopping matrix elements are

$$& h _ { 1 ; \text{rr} ^ { \prime } } \ = \ - \frac { t } { 2 } \sum _ { i = 1, 2, 3 } \left ( \delta _ { \text{rr$$

$$h _ { 2 ; \text{rr} ^ { \prime } } \, = \, - \frac { t } { 2 i } \sum _ { i = 1, 2, 3 } ^ { \,, - \,, \$$

where t denotes the hopping amplitude and { a i } i =1 2 3 , , are the nearest neighbor basis vectors on a lattice. The Fourier components of the hopping matrices read h 1 ( p ) = -t ∑ i =1 ··· 3 cos( p · a i ) and h 2 ( p ) = -t ∑ i =1 ··· 3 sin( p · a i ), and which gives for the kernel of the Hamiltonian H 0 ( p ) = h 1 ( p ) σ 1 + h 2 ( p ) σ . 2 (4)

The tight-binding Hamiltonian Eq. (1) is symmetric under time reversal H T el = H el , which translates into the Fourier space as

$$H _ { 0 } ^ { * } ( - \mathbf p ) = H _ { 0 } ( \mathbf p ), \quad \begin{array} {$$

with the complex conjugation acting only on the sublattice space. A spectral gap appears in the Hamiltonian in the form of the Dirac mass parameter [31] m ∑ r c σ c † r 3 r , and does not affect the property (5).

The dynamics of two-component (in-plane) monochromatic phonons is governed by

$$H _ { \text{ph} } & = \sum _ { \mu = 1, 2 } \sum _ { \text{r} } \frac { 1 } { 2 } ( P _ { \mu ; \text{r} }$$

with dispersionless frequency ω and operators of canonical momentum ( P ) and position ( A )

$$P _ { \mu ; \mathfrak { r } } \, & = \, i \sqrt { \frac { \omega } { 2 } } \sum _ { q } \left [ b _ { \mu, q } e ^ { i q \cdot \mathfrak { r } } - b _ { \mu, q } ^ { \dagger } e ^ { - i q \cdot \mathfrak { r } } \right ], \\ A \quad & = 1 \quad \sum _ { \Gamma _ { L } \quad \mathfrak { r } _ { L } i q \cdot \mathfrak { r } _ { L } \dots i q \cdot \mathfrak { r } _ { L } } \quad \mathfrak { r } _ { L } \circ$$

$$A _ { \mu ; r } & = \, \frac { 1 } { \sqrt { 2 \omega } } \sum _ { q } ^ { \psi } \left [ b _ { \mu, q } e ^ { i q \cdot r } + b _ { \mu, q } ^ { \dagger } e ^ { - i q \cdot r } \right ], \\ \intertext { c a t h o m i o n a t h } \tilde { \sim } \, \tilde { h } \, \text{which anillate} \, \text{thematrix} \, \text{and} \, \text{when} \, \tilde { t } \, \text{notification} \,.$$

expressed in terms of bosonic operators b µ, r ( b † µ, r ) which annihilate (create) a phonon at position r . Phonons interact with electrons via attachment to the lattice bonds

$$H _ { \text{I} } & = U \sum _ { \mu = 1, 2 } \sum _ { \text{r} } c _ { \text{r} } ^ { \dagger } \cdot A _ { \mu ; \text{r} } \sigma _ { \mu } c _ { \text{r} }, \\ \text{d. \, d. \, d } \, \text{r} \, \dots \, \text{d. \, d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d.. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \ \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \ } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d } \, \text{d. \, d }$$

whereas the coupling parameter U can in principle be time dependent and describe a periodic driving. In doing this one can open up new directions toward physics of emergent dynamical Floquet topological phases [32, 33] or non-equilibrium superconductivity [34-36].

## A. Functional integral representation

The coherent state functional integral representation of the zero temperature partition function corresponding to the model Hamiltonian

$$H = H _ { \text{el} } + H _ { \text{ph} } + H _ { \text{l} }$$

reads

$$\mathcal { Z } = \int \mathcal { D } [ \psi ^ { \dagger }, \psi, A ] \exp \{ - \mathcal { S } [ \psi ^ { \dagger }, \psi, A ] \}, \text{attice action}$$

with the Euclidean lattice action

$$\mathcal { S } [ \psi ^ { \dagger }, \psi, A ] = \mathcal { S } _ { \text{el} } [ \psi ^ { \dagger }, \psi ] + \mathcal { S } _ { \text{ph} } [ A ] + \mathcal { S } _ { \text{i} } [ \psi ^ { \dagger }, \psi, A ]. & & ( 1 2 )$$

The elements of the action are

$$\mathcal { S } _ { \text{el} } [ \psi ^ { \dagger }, \psi ] & \, = \, \sum _ { r, r ^ { \prime } } \int d \tau \, \psi _ { r } ^ { \dagger } ( \tau ) [ \partial _ { \tau } \sigma _ { 0 } + H _ { 0 ; \text{rr} } ] \psi _ { r ^ { \prime } } ( \tau ), \\ & \, \sim \, \tau \,. & \, 1 \, \sqrt { \dagger } \, \sqrt { \dagger } \,. \, \tau \,. \, \tau \,. \, \tau \,. \, \tau \,. \, \tau \,. \, \tau \,. \, \tau \,. \, \tau \,. \, \tau \,. \, \tau \,. \, \tau \,. \, \tau \,. \, \tau \,. \, \tau \,. \, \tau \,. \, \tau \,. \, \tau \,. \, \tau \,. \, \tau \,. \, \tau \,. \, \tau \,. \, \tau \,. \, \tau \,. \, \tau \,. \, \tau \,. \, \tuple \,.$$

$$\mathcal { S } _ { 1 } [ \psi ^ { \dagger }, \psi, A ] \, = \, U ^ { \prime } \sum _ { \mu = 1, 2 } ^ { \prime } \sum _ { \mathbf r } \int d \tau \, A _ { \mu ; \mathbf r } ( \tau ) \psi _ { \mathbf r } ^ { \dagger } ( \tau ) \sigma _ { \mu } \psi _ { \mathbf r } ( \tau ). \\ \mathbf r _ { \dots } \, \mathcal { H } \, - \, ( \mathcal { H } \, - \, \mathcal { H } \, \dots \, \mathcal { H } \, - \, \mathcal { H } \, \dots \, \mathcal { H } \, - \, \mathcal { H } \, \dots \, \mathcal { H } \, - \, \mathcal { H } \, \dots \, \mathcal { H } \, - \, \mathcal { H } \, \dots \, \mathcal { H } \, - \, \mathcal { H } \, \dots \, \mathcal { H } \, - \, \mathcal { H } \, \dots \, \mathcal { H } \, - \, \mathcal { H } \, \dots \, \mathcal { H } \, - \, \mathcal { H } \, \dots \, \mathcal { H } \dots \, \mathcal { H } \, - \, \mathcal { H } \, \dots \, \mathcal { H } \, - \, \mathcal { H } \, \dots \, \mathcal { H } \, - \, \mathcal { H } \, \dots \, \mathcal { H } \, - \, \mathcal { H } \, \dots \, \mathcal { H } \, - \, \mathcal { H } \, \dots \, \mathcal { H } \, - \, \mathcal { H } \, \dots \, \mathcal { H } \, - \, \mathcal { H } \, \dots \, \mathcal { H } \, - \, \mathcal { H } \, \dots \, \mathfrak { H } \, - \, \mathcal { H } \, \dots \, \mathfrak { H } \, - \, \mathcal { H } \, \dots \, \mathfrak { H } \, - \, \mathcal { H } \, \dots \, \mathfrak { H } \, - \, \mathcal { H } \, \dots \, \mathfrak { H } \, - \, \mathfrak { H } \, \dots \, \mathfrak { H } \, - \, \mathfrak { H } \, \dots \, \mathfrak { H } \, - \, \mathfrak { H } \, \dots \, \mathfrak { H } \, - \, \mathfrak { H } \, \dots \, \mathfrak { H } \, - \, \mathfrak { H } \, \dots \, \mathfrak { H } \, - \, \mathfrak { H } \, \dots \, \mathfrak { H } \, - \, \mathfrak { H } \, \dots \, \mathfrak { H } \, - \, \mathfrak { H } \, \dots \, \mathfrak { H }$$

$$\mathcal { S } _ { \text{ph} } [ A ] & \, = \, \frac { 1 } { 2 } \sum _ { \mu = 1, 2 } \sum _ { \text{r} } \int d \tau \, \left \{ [ \partial _ { \tau } A _ { \mu ; \text{r} } ( \tau ) ] ^ { 2 } + [ \nabla _ { \text{r} } A _ { \mu ; \text{r} } ( \tau ) ] ^ { 2 } + \omega ^ { 2 } A _ { \mu ; \text{r} } ^ { 2 } ( \tau ) \right \}, \quad & & ( 1 4 ) \\ \mathfrak { c } \, \left ( \int \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \ dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \\ \mathfrak { c } \, \left ( \int \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \ dots \dots \dots \dots \dots \ dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots$$

Here ψ = ( ψ , ψ 1 2 ) represents a two-component time and position dependent Grassmann spinor with the component index referring to the sublattice. In this representation, the model resembles the standard twodimensional Su-Schrieffer-Heeger (SSH) model of electron-phonon interaction [37] with phonons coupling to the electron currents [38, 39]. In contrast to the conventional SSH model our phonons are dispersionless, which allows us to make the crucial approximation by dropping the kinetic term ( ∂ A τ µ ) 2 +( ∇ r A µ ) 2 . The physical picture of this approximation consists in neglecting the slow dynamics of heavy lattice ions due to thermal fluctuations close to the absolute zero of temperature in comparison to the induced dynamics due to the interaction with band electrons. This approximation is justified only if phonons have a relatively large spectral gap, which is the case for optical phonon branches of honeycomb lattices [27, 40-42].

## B. Effective low-energy model

The low-energy part of the spectrum of H 0 defined in Eq. (4) is characterized by two Dirac cones at K ± , which represent Fermi quasiparticles with different chiralities and introduces an additional discrete degree of freedom into the model. This requires operating in the space of 4 × 4 complex matrices Σ ij = σ i ⊗ σ , j i, j = 0 1 2 3. , , , In total there are 15 traceless matrices corresponding to a particular representation of the group of SU (4)-transformations and a 4-dimensional unit matrix Σ 00 . For any complex 4 × 4-matrix Q there is a unique decomposition Q = Q ij Σ . Using this notation the effective ij low-energy Hamiltonian reads

$$H _ { \text{W} } ( \mathbf p ) = p _ { 1 } \Sigma _ { 0 1 } + p _ { 2 } \Sigma _ { 3 2 }.$$

A reasonable translation of the property (5) for the low-energy Hamiltonian is Σ 02 H T W ( -p )Σ 02 = H W ( p ) for the time reversal and Σ 33 H W ( p )Σ 33 = -H W ( p ) for the sublattice (or chiral) symmetry [43].

While in the lattice model defined in the previous paragraph, the phonon field is attached to the lattice bonds, in low-energy approximation it is attached to each sublattice and accounts for the scattering between both Dirac nodes. This is dictated by the properties of the C 6 v group [26, 29] and models the E 1 -in-plane optical modes. The spectrum of these modes reveals a weak alteration over the entire Brillouin zone [27, 40, 41], which makes it possible to model them in the form of dispersionless monochromatic lattice vibrations. The effective low-energy Hamiltonian which describes the Dirac fermions coupled to the in-plane phonons in 2D reads

$$\mathcal { H } = i \vec { \partial } \cdot \vec { j } + \vec { A } \cdot \vec { \Pi }, \quad \ \ \ \ ( 1 7 )$$

where we assembled two particular subsets with two Σ-matrices each to vectors /vector j = { Σ 01 , Σ 32 } and /vector Π = { Σ 13 , Σ 20 } . The set /vector j couples to the kinetic energy of fermions, while Π to the two-component /vector phonon field /vector A . The corresponding Euclidean action becomes:

$$\mathcal { S } [ A, \psi ^ { \dagger }, \psi ] = \frac { 1 } { 2 g } \vec { A } \cdot \vec { A } + \psi ^ { \dagger } \cdot [ \partial _ { \tau } \Sigma _ { 0 0 } + \mathcal { H } ] \, \psi,$$

where ψ = ( ψ 11 , ψ 12 , ψ 21 , ψ 22 ) T are the complex Grassmann fields with the first index referring to the sublattice and the second to the respective Dirac point, and H is given in Eq. (17). The coupling constant g ∼ U /ω 2 2 is related to the inverse frequency of monochromatic phonons [26, 28]. The correlation functions of currents can be obtained from the generating functional

$$\mathcal { W } [ \mathcal { A } ] = \log \langle \exp \{ - \mathcal { S } [ \mathcal { A } ] \} \rangle _ { \psi ^ { \dagger } \psi, A }$$

by repeating variations with respect to the fields A = { /vector α, /vector β } , S A [ ] is the action (18) augmented by an auxiliary source term

$$\psi ^ { \dagger } [ \vec { \alpha } \cdot \vec { j } + \vec { \beta } \cdot \vec { \Pi } ] \psi = \vec { \alpha } \cdot \vec { J } ^ { \text{intra} } + \vec { \beta } \cdot \vec { J } ^ { \text{inter} }$$

with introduced intranodal and internodal current operators

$$\vec { J } ^ { \text{intra} } = \vec { j } = \{ \Sigma _ { 0 1 }, \Sigma _ { 3 2 } \}, \ \vec { J } ^ { \text{inter} } = \vec { \Pi } = \{ \Sigma _ { 1 3 }, \Sigma _ { 2 0 } \}.$$

The averaging operator 〈· · · 〉 ψ ψ,A † denotes the functional integration over all degrees of freedom. Field α is the external gauge field while field β could be an external mechanical modulated strain field. Particularly interesting for symmetry broken phases are non-vanishing average currents ( µ = 1 2) ,

$$\bar { J } _ { \mu, r } ^ { \text{intra} } \ & = \ & - \frac { \delta } { \delta \alpha _ { \mu, r } } \Big | _ { \mathcal { A } = 0 } \mathcal { W } [ \mathcal { A } ], & & ( 2 2 ) \\ \bar { J } _ { \mu, r } ^ { \text{inter} } \ & = \ & - \frac { \delta } { \delta \mathcal { A } } \Big | _ { \mathcal { W } } \quad \mathcal { W } [ \mathcal { A } ]. & & ( 2 3 )$$

∣ Following the procedure developed in Ref. [29], we integrate out the bosonic field /vector A , which creates a four-fermion interaction term. The latter can be decoupled anew by 4 × 4 matrix fields Q :

$$J ^ { \text{mra} } _ { \mu, r } \, = \, - \, \frac { } { \delta \alpha _ { \mu, r } } \Big | _ { \mathcal { A } = 0 } \mathcal { W } [ \mathcal { A } ], & ( 2 2 ) \\ \bar { J } ^ { \text{inter} } _ { \mu, r } \, = \, - \, \frac { \delta } { \delta \beta _ { \mu, r } } \Big | _ { \mathcal { A } = 0 } \mathcal { W } [ \mathcal { A } ]. & ( 2 3 ) \\ \text{developed in Ref. } [ 2 9 ], \, \text{we integrate out the bosonic field } \vec { A }, \, \text{which creates a}$$

$$\frac { 1 } { 2 g } A _ { \mu } ^ { 2 } + A _ { \mu } \psi ^ { \dagger } \Pi _ { \mu } \psi \to - \frac { g } { 2 } \left ( \psi ^ { \dagger } \Pi _ { \mu } \psi \right ) ^ { 2 } = \frac { g } { 2 } \text{tr} \left ( \Pi _ { \mu } \psi ^ { \dagger } \right ) ^ { 2 } \to \frac { 1 } { 2 g } \text{tr} \mathcal { Q } _ { \mu } ^ { 2 } + i \psi ^ { \dagger } \mathcal { Q } _ { \mu } \Pi _ { \mu } \psi, \\ \text{where the summation over } \mu = 1, 2 \text{ is understood and the decoupling field decomposes as } \mathcal { Q } _ { \cdot } = \mathcal { Q } _ { \cdot } ^ { i j } \Sigma _ { i i }.$$

where the summation over µ = 1 2 is understood and the decoupling field decomposes as , Q µ = Q ij µ Σ . ij The full action in this representation reads

$$\mathcal { S } [ \mathcal { Q }, \psi ^ { \dagger }, \psi ] = \frac { 1 } { 2 g } \text{tr} \mathcal { Q } _ { \mu } ^ { 2 } + \psi ^ { \dagger } \cdot \left [ G _ { 0 } ^ { - 1 } + i \mathcal { Q } _ { \mu } \Pi _ { \mu } \right ] \psi, \ \ G _ { 0 } ^ { - 1 } = \partial _ { \tau } \Sigma _ { 0 0 } + i \vec { \partial } \cdot \vec { j }. \\ \text{low we can integrate out the fermions to arrive at the bosonic action}$$

Now we can integrate out the fermions to arrive at the bosonic action

$$\mathcal { S } [ \mathcal { Q } ] = \frac { 1 } { 2 g } \text{tr} \mathcal { Q } _ { \mu } ^ { 2 } - \text{tr} \log \left [ G _ { 0 } ^ { - 1 } + i \mathcal { Q } _ { \mu } \Pi _ { \mu } \right ].$$

## III. SADDLE-POINT ANALYSIS

## A. Effective potential

In simplest mean field approximation we replace quantum fields Q µ in action (26) by corresponding spatially uniform classical background fields M µ

$$V _ { \text{eff} } = \frac { 1 } { 2 g } \text{tr} M _ { \mu } ^ { 2 } - \log \det \left [ G _ { 0 } ^ { - 1 } + i M _ { \mu } \Pi _ { \mu } \right ]. & & ( 2 7 ) \\ \text{ral gap, the gap parameter must couple to a matrix which does not commute with the }$$

To provide a spectral gap, the gap parameter must couple to a matrix which does not commute with the low-energy Hamiltonian H W in Eq. (16). Technically this prevents singularities on the real axis in the Green's functions. This restricts our freedom to M 1 = (∆Σ 01 + m Σ 20 ) / 2 and M 2 = (∆Σ 32 -m Σ 13 ) / 2, such that

$$i M _ { \mu } \Pi _ { \mu } = m \Sigma _ { 3 3 } + \Delta \Sigma _ { 1 2 },$$

FIG. 1: Left panel: Low-energy landscape of the mean-field free energy for g/g c &lt; 1 with the only minimum at m = 0 = ∆. Right panel: Low-energy landscape of the mean-field free energy with visible barriers between stable phases. The model parameter are chosen to give γ = g/g c ∼ 3 5. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00442v2-figure-52989f5d-c195-4349-acce-f09c953f2f3d.jpg)

which indeed anticommutes with H W . There are two further matrices which anticommute with H W : Σ 03 and Σ 22 . It is easy to see though, that the self-consistent equations for the respective order-parameters are unstable. For Σ 03 this is demonstrated in Ref. [29], while the case of Σ 22 is ruled out in the Appendix A.

The minus sign between both terms in Eq. (27) suggests a competition between them. For small g the first term dominates and the effective potential has the shape of a convex hull as it is shown in the left panel of Figure 1. As the interaction becomes larger, both terms become comparable in size until the second term destroys the convexity of the potential. The determinant can be readily evaluated giving

$$\log \det \left [ G _ { 0 } ^ { - 1 } + i M _ { \mu } \Pi _ { \mu } \right ] = \int \frac { d ^ { 3 } Q } { ( 2 \pi ) ^ { 3 } } \ \log \left [ \Delta ^ { 4 } + 2 \Delta ^ { 2 } ( Q ^ { 2 } - m ^ { 2 } ) + ( Q ^ { 2 } + m ^ { 2 } ) ^ { 2 } \right ], \quad \left ( 2 9 \right ) \\ \text{where} \, Q ^ { 2 } = a ^ { 2 } + \mathfrak { n } ^ { 2 }. \ \text{The argument of the loan/thm is not indifferent to the interchange of } m \text{ and } \Delta. \ \text{which}$$

where Q 2 = q 2 0 + q 2 . The argument of the logarithm is not indifferent to the interchange of m and ∆, which suggests the anisotropy of the potential. All integrals diverge and we need to perform the integrations up to a spherical cutoff. The potential landscape plotted in the right panel of Figure 1 reveals a high, yet discrete symmetry. Visually one observes four angular minima, corresponding to two stable phases with degenerated vacua separated from each other by potential walls. On the ridge of the potential one recognizes a fourfold degenerated saddle point like structure, which corresponds to the unstable states with strong tendency to decay in one of both stable phases. The ∆- and the m -phase are associated with extra terms in the tight-binding Hamiltonian (1) in real space representation:

$$H _ { \text{TB} ; \text{rr} ^ { \prime } } + \Delta \sigma _ { 2 } \cos ( \text{G} \cdot \text{r} ) \delta _ { \text{r}, \text{r} ^ { \prime } },$$

where G = K + -K -, and

$$H _ { \text{TB} ; \text{rr} ^ { \prime } } + i m \sigma _ { 3 } \sum _ { j = 1 } ^ { 3 } ( \delta _ { r ^ { \prime }, \text{r} + \text{c} _ { j } } - \delta _ { r ^ { \prime }, \text{r} - \text{c} _ { j } } ), & & ( 3 1 ) \\ \dots & \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \.$$

{ c j } j =1 2 3 , , denoting nearest neighbor basis vectors on one sublattice [11, 45]. The ∆-phase is characterized by a strain field which modulates with the nodal wave vector G (modulated strain phase), while the m -phase is a flux phase that is commensurate with one sublattice (Haldane phase). Both phases are subject to time reversal and sublattice symmetry breaking.

The minima corresponding to both phases are equally deep, i.e. the phases are evenly likely to emerge. A natural question to ask is therefore, if this observation points to a possible symmetry between both phases. Then there would be a unitary transformation which maps one vacuum on another and vice versa. It is indeed possible to construct a unitary transformation which transforms Σ 12 on Σ 33 , i.e.

$$T ^ { \dagger } \Sigma _ { 1 2 } T = \Sigma _ { 3 3 }.$$

FIG. 2: Two examples for broken symmetry between both ground states. Upper row, left panel: Two perturbations which anticommute with the bare electronic Hamiltonian make the ∆-phase less pronounced compared to the m -phase. They couple to matrices Σ 03 (staggered chemical potential) and Σ 22 . Upper row, right panel: The situation which favors the ∆-phase can be realized with two further perturbations which anticommute with the bare electronic Hamiltonian. The couple to Σ 33 and Σ 12 phases. Depending on the sign of the perturbation, the degeneracy between the two symmetric minima is lifted as well.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00442v2-figure-7cb2c176-d923-4a1e-9f0f-6101a22ea4eb.jpg)

Explicitly we get

$$T ^ { \dagger } = \frac { 1 } { 2 } \left ( \begin{array} { c c c c } 1 & - i & 1 & - i \\ 1 & i & 1 & i \\ 1 & - i & - 1 & i \\ 1 & i & - 1 & - i \end{array} \right ), \ T = \frac { 1 } { 2 } \left ( \begin{array} { c c c c } 1 & 1 & 1 & 1 \\ i & - i & i & - i \\ 1 & 1 & - 1 & - 1 \\ i & - i & - i & i \end{array} \right ). \\ \text{information} \, T ^ { \dagger } \text{ is not a symmetrv of the bare electronic Hamiltonian. nor of the mean-field}$$

The transformation T † is not a symmetry of the bare electronic Hamiltonian, nor of the mean-field Hamiltonian, but because of its unitarity leaves the eigenvalues unchanged.

Technically, the simplest way to break the symmetry between both stable ground states is by a perturbation which anticommutes with the bare electronic Hamiltonian Eq. (16). This property guaranties, that the determinant of the perturbed Hamiltonian remains rotationally invariant which greatly simplifies the calculations. If the perturbation mimics the order parameter, i.e. it couples either to Σ 12 or Σ 33 matrix, playing the role of the remanent density of respective fields, then the strain phase (order parameter ∆) is energetically more preferable and even the Z 2 -symmetry between the both minima is violated, as it can be seen in the right panel of Figure 2. The two other possible perturbations couple to either Σ 03 or Σ 22 and favor the flux phase (order parameter m ). In this case the Z 2 -degeneracy of the ground state remains present, cf. left panel of Figure 2. On the lattice the corresponding contributions to the Hamiltonian are: if Σ 03 is chosen

$$\delta H _ { \mathbf p } ^ { ( 1 ) } = u _ { 1 } \sigma _ { 3 } \delta _ { \mathbf r, \mathbf r ^ { \prime } },$$

and if Σ 22

where G = K + -K -.

## B. Solution of the saddle-point equations

A more detailed picture of the symmetry breaking follows from the solutions of mean field equations. They are obtained if potential (27) is varied with respect to M α :

$$M _ { \alpha } = i g \Pi _ { \alpha } \left [ G _ { 0 } ^ { - 1 } + i M _ { \mu } \Pi _ { \mu } \right ] _ { r r } ^ { - 1 }.$$

$$\delta H _ { \mathrm p } ^ { ( 2 ) } = i u _ { 2 } \sigma _ { 2 } \sin ( \mathbf G \cdot \mathbf r ) \delta _ { \mathrm r, \mathbf r ^ { \prime } },$$

Inserting the order parameter matrix in (36), performing the frequency integration from -∞ to + ∞ and the radial integral up to the cutoff Λ, related to the band width we get

$$m \pm \Delta = \gamma ( m \pm \Delta ) \left [ \sqrt { 1 + ( m \pm \Delta ) ^ { 2 } } - | m \pm \Delta | \right ], & & ( 3 7 ) \\ \L / 2 \pi \text{ and } m \text{ and } \Delta \text{ are rescaled in units of } \Lambda. \text{ First we notice that } m = \Delta \text{ might be a}$$

where γ = g Λ 2 / π and m and ∆ are rescaled in units of Λ. First we notice that m = ∆ might be a solution for both +, as well as trivially for -sign. This also includes m = 0 = ∆ case. For + sign, the non-trivial solutions are extracted from

$$1 = \gamma [ \sqrt { 1 + 4 m ^ { 2 } } - 2 | m | ], & & ( 3 8 ) \\ \text{ration strength exceeds a critical value of the order of inverse band width}$$

which exists only if the interaction strength exceeds a critical value of the order of inverse band width γ = g/g c /greaterorequalslant 1, g c = 2 π/ Λ

$$| m | = | \Delta | = \frac { \gamma ^ { 2 } - 1 } { 4 \gamma } \Theta ( \gamma - 1 ).$$

/negationslash

This solution corresponds to the state which is visible on the ridge of the phase separating walls in Figure 1. This state is stable only classically (mean-field) and is destroyed by quantum fluctuations as it will be shown below. Further solutions are there for m = 0, ∆ = 0 and for m = 0, ∆ = 0. For ∆-phase we get

/negationslash

$$| \Delta | = \frac { \gamma ^ { 2 } - 1 } { 2 \gamma } \Theta ( \gamma - 1 ),$$

and the same for m -phase. These solutions correspond to the four angular minima visible in Figure 1. They turn out to be stable against quantum fluctuations as well.

## IV. QUANTUM FLUCTUATIONS

Expanding Eq. (26) to the second order in the fluctuations around ground states of respective phases we get

$$\mathcal { S } _ { G } [ \mathcal { Q } ] = \text{tr} \left \{ \frac { 1 } { 2 g } \mathcal { Q } _ { \mu } ^ { 2 } - \frac { 1 } { 2 } G \mathcal { Q } _ { \alpha } \Pi _ { \alpha } G \mathcal { Q } _ { \beta } \Pi _ { \beta } \right \}, \\ \dots \dots \dots \dots \dots.$$

which can be brought further into the conventional quadratic form

$$\mathcal { S } _ { G } [ \mathcal { Q } ] = \int \frac { d ^ { 3 } P } { ( 2 \pi ) ^ { 3 } } \, \left ( \begin{array} { c } \vec { \mathcal { Q } } _ { 1 } \\ \vec { \mathcal { Q } } _ { 2 } \end{array} \right ) _ { P } ^ { T } \left ( \begin{array} { c } \mathbf A & \mathbf B \\ \mathbf C & \mathbf D \end{array} \right ) _ { P } \left ( \begin{array} { c } \vec { \mathcal { Q } } _ { 1 } \\ \vec { \mathcal { Q } } _ { 2 } \end{array} \right ) _ { - P }$$

with the vector fields

$$\vec { \mathcal { Q } } _ { i } = ( \mathcal { Q } _ { i } ^ { 0 0 }, \mathcal { Q } _ { i } ^ { 0 1 }, \mathcal { Q } _ { i } ^ { 0 2 }, \mathcal { Q } _ { i } ^ { 0 3 }, \mathcal { Q } _ { i } ^ { 1 0 }, \mathcal { Q } _ { i } ^ { 1 1 }, \mathcal { Q } _ { i } ^ { 1 2 }, \mathcal { Q } _ { i } ^ { 1 3 }, \mathcal { Q } _ { i } ^ { 2 0 }, \mathcal { Q } _ { i } ^ { 2 1 }, \mathcal { Q } _ { i } ^ { 2 2 }, \mathcal { Q } _ { i } ^ { 3 0 }, \mathcal { Q } _ { i } ^ { 3 1 }, \mathcal { Q } _ { i } ^ { 3 2 }, \mathcal { Q } _ { i } ^ { 3 3 } ) \,, \ i = 1, 2. \quad ( 4 3 ) \\ \text{The kernel matrix of the quadratic form, (i.e. inverse bosonic propagator matrix)}$$

The kernel matrix of the quadratic form, (i.e. inverse bosonic propagator matrix)

$$\mathbf P ^ { - 1 } ( P ) = \left ( \begin{array} { c c } \mathbf A & \mathbf B \\ \mathbf C & \mathbf D \end{array} \right ) _ { P }$$

/negationslash is a 32 × 32 matrix, i.e., each matrix block itself is a 16 × 16 matrix. In accord with the intuitive picture drawn from the analysis of the effective potential, the momentum-independent terms in Gaussian action (42) (which are commonly called the Proca matrices) of both m -phase and ∆-phase are indeed strictly positive, while in the mixed phase between them with m = ∆ = 0 found in Eq. (39) there are negative eigenvalues, rendering this phase unstable, cf. Appendix B. Because the Proca matrix has no zero eigenvalues the kernel matrix in (42) can be inverted even for zero momentum. An astonishing fact about the Proca matrices is that some of their elements sum up to multiples of γ -1 , cf. Appendix B. In

calculating the Hall conductivities, Proca matrices guarantee the convergence of the functional integral but their elements do not appear in the ultimate result.

In exploring the structure of the gradient expansion we restrict our analysis to the linear terms only. In the usual three-dimensional electrodynamics it is referred to as the induced Chern-Simons terms [1821]. These terms have several important properties: they are odd under time-reversal and because they are not associated with any length scale they represent pristine topological excitations with both scale and conformal invariance. Because of these qualities they are believed to represent the effective field theoretical description of the quantum Hall effect [12]. In our formalism things are getting more involved. At first glance the frequency dependent part of what is going to become the Chern-Simons tensor reads

$$\mathcal { S } _ { C S } ^ { ( 0 ) } & = \int \frac { d ^ { 3 } P } { ( 2 \pi ) ^ { 3 } } \ p _ { 0 } \mathcal { Q } _ { \mu, P } ^ { a b } \mathcal { Q } _ { \nu, - P } ^ { \alpha \beta } \ \mathcal { P } _ { a b | \alpha \beta } ^ { \mu \nu },$$

where

$$\mathcal { P } _ { a b | \alpha \beta } ^ { \mu \nu } \, = \, \frac { i } { 2 } \text{Tr} \int \frac { d ^ { 3 } Q } { ( 2 \pi ) ^ { 3 } } \, G ( Q ) \Sigma _ { a b } \Pi _ { \mu } G ( Q ) G ( Q ) \Sigma _ { \alpha \beta } \Pi _ { \nu } \\ i \, \infty \, \left ( \hat { \sim } \, \infty \, \pi \, \infty \, \pi \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inyeset \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \. \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \}$$

with Σ = sgn( ˆ m )Σ 33 in m -phase and Σ = sgn(∆)Σ ˆ 12 in ∆-phase. For spatial components we get at first

$$\varkappa = \frac { \mu } { 3 2 \pi } & \text{$\nu$} \quad \cdots \quad / \\ = \, \frac { i } { 3 2 \pi } & \text{Tr} \left \{ \hat { \Sigma } \left ( \Sigma _ { a b } \Pi _ { \mu } \Sigma _ { \alpha \beta } \Pi _ { \nu } - \Sigma _ { \alpha \beta } \Pi _ { \nu } \Sigma _ { a b } \Pi _ { \mu } \right ) \right \}, \\ \varkappa \text{ in $m$-phase and $\hat{\Sigma} = \text{sgn}(\Delta)\Sigma_{\alpha}$ in $\Delta$-phase. For spatial components we get at first}$$

$$\mathcal { S } _ { C S } & = \int \frac { d ^ { 3 } Q } { ( 2 \pi ) ^ { 3 } } \ \mathcal { Q } _ { \mu, P } ^ { a b } \mathcal { Q } _ { \nu, - P } ^ { \alpha \beta } \ p _ { i } \cdot \mathcal { T } _ { a b | \alpha \beta } ^ { i, \mu \nu },$$

$$= \begin{array} { c c c } \bar { \mu } ^ { \downarrow } \quad \bar { \nu } ^ { \dots \, \bar { \nu } } \\ - \frac { 1 } { 3 2 \pi } \text{Tr} \left \{ \hat { \Sigma } \left ( \Sigma _ { a b } \Pi _ { \mu } j _ { i } \Sigma _ { \alpha \beta } \Pi _ { \nu } - \Sigma _ { \alpha \beta } \Pi _ { \nu } j _ { i } \Sigma _ { a b } \Pi _ { \mu } \right ) \right \}. \\ \text{ification of the above expressions becomes possible with the knowledge of certain distinct} \end{array}$$

where

$$\mathcal { T } _ { a b | \alpha \beta } ^ { i, \mu \nu } \, = \, \frac { 1 } { 2 } \text{Tr} \int \frac { d ^ { 3 } Q } { ( 2 \pi ) ^ { 3 } } \, G ( Q ) \Sigma _ { a b } \Pi _ { \mu } G ( Q ) j _ { i } G ( Q ) \Sigma _ { \alpha \beta } \Pi _ { \nu } \\ 1 \, \sum \, f \hat { \sim } \, \nu \, \cdots \, \pi \, \cdots \, \pi \, \cdots \, \pi \, \cdots \, \pi \, \cdots \, \pi \, \cdots \, \pi \, \cdots \, \pi \, \cdots \, \pi \, \cdots \, \pi \, \cdots \, \pi \, \cdots \, \pi \, \cdots \, \pi \, \cdots \, \pi \, \cdots \, \pi \, \cdots \, \pi \, \cdots \, \pi \, \cdots \, \pi \, \cdots \, \pi \, \cdots \, \pi \, \cdots \, \pi \, \cdots \, \pi \, \cdots \, \pi \, \cdots \, \pi \end{matrix}$$

A decisive simplification of the above expressions becomes possible with the knowledge of certain distinct features of Σ-matrices. The full set of 16 Σ-matrices subdivides into four subsets of four matrices each, with three of them anticommuting among each other and the fourth commuting with all. Calling them Λ a i with a = 1 2 3 4 and , , , i = 0 1 2 3 and imposing the normalization condition Tr[ , , , Λ Λ Λ Λ a 0 a 1 a 2 a 3 ] = 4 for all a , the subsets become

$$\Lambda ^ { 1 } & = \{ i \Sigma _ { 0 0 }, - \Sigma _ { 0 1 }, - \Sigma _ { 3 2 }, - \Sigma _ { 3 3 } \}, \quad \Lambda ^ { 2 } = \{ - i \Sigma _ { 2 1 }, \Sigma _ { 2 0 }, - \Sigma _ { 1 3 }, \Sigma _ { 1 2 } \}, \\ \Lambda ^ { 3 } & = \{ \Sigma _ { \dots } \, - i \Sigma _ { \dots } \, i \Sigma _ { \dots } \, - i \Sigma _ { \dots } \} \quad \Lambda ^ { 4 } = \{ - i \Sigma _ { \dots } \, \Sigma _ { \dots } \, \Sigma _ { \dots } \, - \Sigma _ { \dots } \}$$

$$A ^ { 3 } = \{ \Sigma _ { 3 0 }, - i \Sigma _ { 3 1 }, i \Sigma _ { 0 2 }, - i \Sigma _ { 0 3 } \}, \quad A ^ { 4 } = \{ - i \Sigma _ { 1 1 }, \Sigma _ { 1 0 }, \Sigma _ { 2 3 }, - \Sigma _ { 2 2 } \}.$$

One recognizes that Λ 1 and Λ 2 include vectors /vector j and Π of the model defined in (17). /vector Linear terms of gradient expansion of the action (41) in both phases can be brought into the form structurally similar to the conventional Chern-Simons terms. The crucial difference is that rather than in terms of single components of Q -fields, the Chern-Simons-like structures appear in terms of

$$\Lambda _ { i } ^ { a } = \text{Tr} [ \Lambda _ { i } ^ { a } ( \mathcal { Q } _ { \mu } \Pi _ { \mu } ) ]$$

for each set Λ a . In m -phase it resembles the conventional Chern-Simons terms, giving one Chern-Simons term for each set

$$\mathcal { S } _ { C S } ^ { m } = S _ { m } \epsilon _ { i j k } \Lambda _ { i } ^ { a } \cdot i \partial _ { j } \Lambda _ { k } ^ { a },$$

where the dot-product implies the integration over the entire 3d space-time, i, j, k = 0 1 2, and , , S m = sgn( m / ) 64 π . The fact that in Eq. (54) the derivatives with respect to spatial degrees of freedom appear alongside with the temporal can be interpreted as the restoration of the full gauge invariance, which formally lacks in action Eq. (18) [28].

For the ∆-phase we get two morphologically different terms. The first one resembles the generalized multi-field Chern-Simons term [23]

$$\mathcal { S } _ { C S } ^ { \Delta } = S _ { \Delta } \epsilon _ { i j k } K _ { a b } \Lambda _ { i } ^ { a } \cdot i \partial _ { j } \Lambda _ { k } ^ { b }, \quad \ \.$$

where upper indices run only over a, b = 1 2 and diagonal elements of symmetric tensor , K are zero. Effectively it connects the intranodal and internodal currents with each other. The two remaining matrix sets Λ 3 and Λ 4 interact with each other in an unexpected way

$$\tilde { \mathcal { S } } _ { C S } ^ { \Delta } & = S _ { \Delta } K _ { a b } \left [ \Lambda _ { 3 } ^ { a } \cdot i \partial _ { i } \Lambda _ { i } ^ { b } - \Lambda _ { i } ^ { b } \cdot i \partial _ { i } \Lambda _ { 3 } ^ { a } \right ], & ( 5 6 ) \\ \intertext { d a, b = 3, 4. \, \begin{array} { c } \left ( \begin{array} { c } \left ( \end{array} { c } \left ( \end{array} { c } \left ( \begin{array} { c } \left ( \begin{array} { c } \left ( \end{array} { c } \left ( \begin{array} { c } \left ( \end{array} { c } \left ( \begin{array} { c } \left ( \end{array} { c } \left ( \begin{array} { c } \left ( \end{array} { c } \left ( \begin{array} { c } \left ( \end{array} { c } \left ( \begin{array} { c } \left ( \end{array} { c } \left ( \begin{array} { c } \left ( \end{array} { c } \left ( \begin{array} { c } \left ( \end{array} { c } \left ( \begin{array} { c } \left ( \end{array} { c } \left ( \begin{array} { c } \left ( \end{ array} { c } \left ( \begin{array} { c } \left ( \end{array} { c } \left ( \begin{array} { c } \left ( \end{array} { c } \left ( \begin{array} { c } \left ( \end{array} { c } \left ( \begin{array} { c } \left ( \end{array} { c } \left ( \begin{array} { c } \left ( \end{array} { c } \left ( \begin{array} { c } \left ( \end{array} { c } \left ( \begin{array} { c } \left ( \end{array} { c } \left ( \begin{array} { c } \left ( \end{array} { c } \left ( \begin{array} { c } \left ( \end{array} \\ \left ( \begin{array} { c } \left ( \end{array} { c } \left ( \begin{array} { c } \left ( \end{array} { c } \left ( \begin{array} { c } \left ( \end{array} { c } \left ( \begin{array} { c } \left ( \end{array} { c } \left ( \begin{array} { c } \left ( \end{array} { c } \left ( \begin{array} { c } \left ( \end{array} { c } \left ( \begin{array} { c } \left ( \end{array} { c } \left ( \begin{array} { c } \left ( \end{array} { c } \left ( \begin{array} { c } \left ( \end{array} { c } \right ) } \end{array} } \end{array} } \end{array} }$$

where i = 0 1 2 and , , a, b = 3 4. The two terms (55) and (56) are hardly discussed in context of induced , Chern-Simons terms [12, 21]. From Eq. (54) and (55) it is possible to construct a generalized ChernSimons-like tensor for both phases. For this we retain only two Pauli sets Λ 1 and Λ 2 which are connected to the effective model (17) in the sense that they contain both intra- and internode currents. With α, β = 1 2 we can combine them in one single expression ,

$$\mathcal { S } _ { C S } = \epsilon _ { i j k } T _ { c d } \mathbf \Lambda _ { i } ^ { c } \cdot i \partial _ { j } \mathbf \Lambda _ { k } ^ { d },$$

where the diagonal elements of the symmetric tensor T cc = sgn( m / ) 64 π and non-diagonal ones T cd = sgn(∆) 64 / π . Together with the Proca matrices of respective phases given in Appendix B, Eq. (57) is used to calculate the topological response function specific to each identified phase.

## A. Generalized Hall response

We introduce a generalized current-current correlation function

$$\langle J _ { x, \mu } ^ { c } J _ { y, \nu } ^ { d } \rangle _ { \psi ^ { \dagger } \psi, A } = \frac { \delta } { \delta A _ { y, \nu } ^ { d } } \frac { \delta } { \delta A _ { x, \mu } ^ { c } } \Big | _ { \mathcal { A } = 0 } \mathcal { W } [ \mathcal { A } ], \\ \text{lenotes the normalized functional integration with respect to the action in Eq. (18) and}$$

∣ where 〈· · · 〉 ψ ψ,A † denotes the normalized functional integration with respect to the action in Eq. (18) and W A [ ] the generating functional of current-current correlator functions defined in Eq. (19). The lower current indices refer to the position and time coordinates of the local current on the lattice ( x, y ) and to their spatial direction ( µ, ν ). Upper indices ( a, b ) refer to the type of the current and distinguishes between inter- and intranode currents. The index combination c = d = 1 corresponds to the correlator of usual intranodal currents,

$$\langle ( \psi ^ { \dagger } j _ { \mu } \psi ) _ { x } ( \psi ^ { \dagger } j _ { \nu } \psi ) _ { y } \rangle _ { \psi ^ { \dagger } \psi, A } \to \frac { 1 } { 4 \gamma ^ { 2 } } \langle \text{Tr} \left ( j _ { \mu } \Pi _ { a } \mathcal { Q } _ { a } \right ) _ { x } \text{Tr} \left ( j _ { \nu } \Pi _ { b } \mathcal { Q } _ { b } \right ) _ { y } \rangle \to \frac { 1 } { 4 \gamma ^ { 2 } } \langle \Lambda _ { \mu, x } ^ { 1 } \Lambda _ { \nu, y } ^ { 1 } \rangle _ { g }, \\ \text{where in the second term we brought it into the $\mathcal{Q}$-matrix representation and introduced in the third}$$

where in the second term we brought it into the Q -matrix representation and introduced in the third Λ -notation according to Eq. (53). The operator 〈· · · 〉 Q denotes the average over the non-linear action Eq. (26) while 〈· · · 〉 G over its Gaussian approximation Eq. (41). In the same manner we recognize that c = d = 2-case represents the correlator internodal currents. The response in the ∆-phase with c = 1, d = 2 and vice versa represents correlators over products of an intra- and an internodal current flowing in perpendicular directions. With all four terms we get for Eq. (58)

$$\langle J _ { x, \mu } ^ { c } J _ { y, \nu } ^ { d } \rangle _ { \psi ^ { \dagger } \psi, A } & \approx \frac { 1 } { 4 \gamma ^ { 2 } } \langle \mathbf A _ { \mu, x } ^ { c } \mathbf A _ { \nu, y } ^ { d } \rangle _ { \mathcal { G } }.$$

The physical content of this correlator is simple: We create a local current density of the sort c at position x in direction µ . Then we measure the local current of type d at position y in direction ν . This suggests that not only correlators between pairs of intra-(inter-)node currents are possible, but also mixed correlators of the type 〈 J intra x,µ J inter y,ν 〉 . With explicit knowledge of Proca matrices for each phase and of the generalized Chern-Simons term Eq. (57), the correlator of each two currents can be evaluated giving to linear order in gradient expansion, c.f. Appendix C

$$\frac { 1 } { 4 \gamma ^ { 2 } } \langle \Lambda _ { i, r } ^ { c } \Lambda _ { j, r ^ { \prime } } ^ { d } \rangle _ { \mathcal { G } } = 3 2 T _ { c d } i \epsilon _ { i j k } \partial _ { k } \delta ( r - r ^ { \prime } ).$$

In deriving this expression we used the fact that some elements of the Proca matrices sum up to multiples of γ -1 . This means that while the presence of the Proca matrix is absolutely indispensable in order to carry out the functional integration over the Q -fields, its elements do not appear in the correlator.

The current-current correlator is related to the Hall conductivity via the Kubo formula [47]

$$\bar { \sigma } _ { \mu \nu } ^ { c d } & = \lim _ { \omega \to 0 } \frac { 2 \pi } { \omega } \int d ^ { 3 } x \left [ e ^ { - i \omega ( x _ { 0 } - y _ { 0 } ) } - 1 \right ] \langle J _ { x, \mu } ^ { c } J _ { y, \nu } ^ { d } \rangle _ { \psi ^ { \dagger } \psi, A }. \\ \text{rough the above two limitations} \, f i n d$$

Going through the algebra we ultimately find

$$\bar { \sigma } _ { \mu \nu } ^ { c d } = 6 4 \pi T _ { c d } \epsilon _ { \mu \nu }.$$

Diagonal elements of (63) give the quantized response sgn( m /epsilon1 ) µν in the m -phase and the two off-diagonal components give the quantized response sgn(∆) /epsilon1 µν of the ∆-phase. c = d = 1 corresponds to the usual intranodal Hall conductivity of two Dirac cones, while c = d = 2 represents the average of two internodal currents, which gives the internodal Hall conductivity. The response in the ∆-phase with c = 1, d = 2 and vice versa represents the averages over products of an intra- and an internodal current flowing in perpendicular directions, which too has a quantized response.

## V. DISCUSSIONS

The low-energy behavior of free fermions on the honeycomb lattice is determined by separated spectral nodes. Because of these one can distinguish between currents directly linked to each node, i.e. the intranodal currents, and those between the nodes, i.e. the internodal currents. The system of such fermions coupled to the in-plane phonons can be formally interpreted as a gauge field theory. At sufficiently large electron-phonon interaction strength a transition into a structurally different lattice phase takes place, as visualized in Figure 1. We identify two distinct ground states with order parameters which break both the time reversal and the sublattice symmetries. The effective phonon model exhibits the Chern-Simons-like terms in both phases with coefficients related to quantized Hall conductivities. A similar effect of spontaneous time-reversal symmetry breaking was previously discussed for a time-periodic driven quantum system, where also a topological state with non-zero Chern numbers was observed [32]. The similarity of the two cases is reflected by the mean-field Hamiltonian (30), which agrees with the effective Floquet Hamiltonian in low energy approximation [32, 33]. This indicates that the time-dependent phonon field plays a similar role as a time-periodic driving field.

For the experimental observation of the proposed effects we suggest the following setup: the measurement could be performed on a suspended tightened graphene sample which must be sufficiently flat in order to rule out disorder due to ripples and out-of-plane phonon modes. The in-plane phonon modes could be generated mechanically by applying modulated strain with variable frequency on the sample and by adjusting carefully the strain frequency to the energy of the E 1 -mode. Estimations based on ab-initio calculations vary between 0 15 . -0 2eV [26, 41], which is well within the range of infrared Raman scattering . technique [42]. In order to measure the Hall conductivities a weak electric gradient should be applied to the sample. In the case of the ∆-phase, both electric and strain fields should be applied in the same direction to create the proposed quantized response.

## ACKNOWLEDGMENTS

Illuminating discussions with John Chalker, Bertrand Halperin, Eugene Demler, and Alexei Tsvelik are gratefully acknowledged. We thank Aditi Mitra for drawing our attention to connections between timeperiodically driven quantum systems and systems with electron-phonon interaction. This research was supported by the grants of the Julian Schwinger Foundation for Physics Research.

[1] E. Fradkin, Field theories of condensed matter physics , Cambridge University Press, Cambridge, UK (2013).

- [2] X. G. Wen, Vacuum degeneracy of chiral spin states in compactified space , Phys. Rev. B 40 , 7387(R) (1989).
- [3] D. J. Thouless. M. Kohmoto, M. P. Nightingale, and M. den Nijs, Quantized Hall Conductance in a TwoDimensional Periodic Potential , Phys. Rev. Lett. 49 , 405 (1982).
- [4] G. E. Volovik, Quantized Hall effect in superfluid Helium-3 film , Phys. Lett. A 128 , 277 (1988).
- [5] G. E. Volovik, An analog of the quantum Hall effect in a superfluid He 3 film , Zh. Eksp. Teor. Fiz. 94 , 123 (1988) [Sov. Phys. JETP 67, 1804 (1988)].
- [6] D. Vollhardt and P. W¨lfle, o The Superfluid Phases of Helium 3 , Taylor &amp; Francis, London, (1990).
- [7] Y. Aharonov and D. Bohm, Significance of electromagnetic potentials in quantum theory , Phys. Rev. 115 , 485 (1959).
- [8] M. V. Berry, Quantal Phase Factors Accompanying Adiabatic Changes , Proc. R. Soc. A 392 , 45 (1984).
- [9] B. I. Halperin, Quantized Hall conductance, current-carrying edge states, and the existence of extended states in a two-dimensional disordered potential , Phys. Rev. B 25 , 2185 (1982).
- [10] V. M. Yakovenko, Chern-Simons terms and n field in Haldane's model for the quantum Hall effect without Landau levels , Phys. Rev. Lett. 65 , 251 (1990).
- [11] F. D. M. Haldane, Model for a quantum Hall effect without Landau levels: Condensed-matter realization of the 'parity anomaly' , Phys. Rev. Lett. 61 , 2015 (1988).
- [12] A. Zee, Quantum hall fluids , edited by H. B. Geyer Field Theory, Topology and Condensed Matter Physics. Lecture Notes in Physics (Springer, Berlin, Heidelberg, 1995), Vol. 456 , 99-153.
- [13] J. Fr¨hlich and A. Zee, o Large scale physics of the quantum Hall fluid , Nucl. Phys. B, 364 , 517 (1991).
- [14] H. Hettmansperger, F. Duerr, J. B. Oostinga, C. Gould, B. Trauzettel, and L. W. Molenkamp, Quantum Hall effect in narrow graphene ribbons , Phys. Rev. B 86 , 195417 (2012).
- [15] E. M. Hankiewicz, J. Li, T. Jungwirth, Q. Niu, S.-Q. Shen, and J. Sinova, Charge Hall effect driven by spin-dependent chemical potential gradients and Onsager relations in mesoscopic systems , Phys. Rev. B 72 , 155305 (2005).
- [16] P. Jacquod, R. S. Whitney, J. Meair, and M. B¨ttiker, u Onsager relations in coupled electric, thermoelectric, and spin transport: The tenfold way , Phys. Rev. B 86 , 155118 (2012).
- [17] N. Nagaosa, J. Sinova, S. Onoda, A. H. MacDonald, and N. P. Ong, Anomalous Hall effect , Rev. Mod. Phys. 82 , 1539 (2010).
- [18] A. N. Redlich, Parity violation and gauge noninvariance of the effective gauge field action in three dimensions , Phys. Rev. D 29 , 2366 (1984).
- [19] R. Jackiw, Fractional charge and zero modes for planar systems in a magnetic field , Phys. Rev. D 29 , 2375 (1984).
- [20] E. Fradkin and F. Schaposnik, The fermion-boson mapping in three-dimensional quantum field theory , Phys. Lett. B 338 , 253 (1994).
- [21] G. Dunne, Aspects of Chern-Simons theory , in Topological aspects of low dimensional systems , edited by A. Comtet, T. Jolicoeur, S. Ouvry, and F. David, Les Houches Summer School, (Springer, Berlin, Heidelberg, 1999), Vol. 69 , 177-263.
- [22] X.-L. Qi, T. L. Hughes, and S. C. Zhang, Topological field theory of time-reversal invariant insulators , Phys. Rev. B 78 , 195424 (2008).
- [23] Y.-M. Lu and A. Vishwanath, Theory and classification of interacting integer topological phases in two dimensions: A Chern-Simons approach , Phys. Rev. B 86 , 125119 (2012).
- [24] E. Zohar, J. I. Cirac, and B. Reznik, Quantum simulations of lattice gauge theories using ultracold atoms in optical lattices , Rep. Prog. Phys. 79 , 014401 (2016).
- [25] Sh. Heidari, A. Cortijo, and R. Asgari, Hall viscosity for optical phonons , Phys. Rev. B 100 , 165427 (2019).

[26]

D. M. Basko and I. L. Aleiner,

Interplay of

Coulomb and electron-phonon interactions in graphene

,

Phys.

Rev.

77

,

041409(R) (2008).

- [27] B. Amorim, A. Cortijo, F. de Juan, A. G. Grushin, F. Guinea, A. Gutirrez-Rubio, H. Ochoa, V. Parente, R. Rold´ an, P. San-Jos´, J. Schiefele, M. Sturla, and M. A. H. Vozmediano, e Novel effects of strains in graphene and other two dimensional materials , Phys. Rep. 617 , 1 (2016).
- [28] A. Sinner and K. Ziegler, Emergent Chern-Simons excitations due to electron-phonon interaction , Phys. Rev. B 93 , 125112 (2016).
- [29] A. Sinner and K. Ziegler, Spontaneous mass generation due to phonons in a two-dimensional Dirac fermion system , Ann. Phys. 400 , 262 (2019).
- [30] P. K. Wallace, The band theory of graphite , Phys. Rev. 71 , 622 (1947).
- [31] G. W. Semenoff, Condensed-matter simulation of a three-dimensional anomaly , Phys. Rev. Lett. 53 , 2449 (1984).
- [32] N.H. Lindner, G. Rafael, and V. Galitski, Floquet topological insulator in semiconductor quantum wells , Nature Phys. 7 , 490 (2011).
- [33] J. Cayssol, B. D´ra, F. Simon and R. Moessner, o Floquet topological insulators , Phys. Status Solidi RRL 7 , 101-108 (2013).
- [34] M. Knap, M. Babadi, G. Refael, I. Martin, and E. Demler, Dynamical Cooper pairing in nonequilibrium electron-phonon systems , Phys. Rev. B 94 , 214504 (2016).

- [35] M. Babadi, M. Knap, G. Refael, I. Martin, and E. Demler, Theory of parametrically amplified electron-phonon superconductivity , Phys. Rev. B 96 , 014512 (2017).
- [36] M. Buzzi, G. Jotzu, A. Cavalleri, J. I. Cirac, E. A. Demler, B. I. Halperin, M. D. Lukin, T. Shi, Y. Wang, and D. Podolsky, Higgs-mediated optical amplification in a non-equilibrium superconductor , arXiv:1908.10879.
- [37] W. P. Su, J. R. Schrieffer, and A. J. Heeger, Solitons in polyacetylene , Phys. Rev. Lett. 42 , 1698 (1979).
- [38] M. Zoli, Path integral computation of phonon anharmonicity , Eur. Phys. J. B 40 , 79 (2004).
- [39] M. Zoli, Path integral of the two-dimensional Su-Schrieffer-Heeger model , Phys. Rev. B 71 , 205111 (2005).
- [40] A. C. Ferrari, Raman spectroscopy of graphene and graphite: Disorder, electron-phonon coupling, doping and nonadiabatic effects , Sol. St. Comm. 143 , 47 (2007).
- [41] G. D. Sanders, A. R. T. Nugraha, K. Sato, J.-H. Kim, J. Kono, R. Saito and C. J. Stanton, Theory of coherent phonons in carbon nanotubes and graphene nanoribbons , J. Phys.: Condens. Matter 25 , 144201 (2013).
- [42] A. Jorio, Raman Spectroscopy in Graphene-Based Systems: Prototypes for Nanoscience and Nanometrology , ISRN Nanotechnology, Vol. 2012 , pp. 1-12 (2012).
- [43] A.W.W. Ludwig, M.P.A. Fisher, R. Shankar, and G. Grinstein, Integer quantum Hall transition: An alternative approach and exact results , Phys. Rev. B 50 , 7526 (1994).
- [44] Y. Hatsugai and M. Kohmoto, Energy spectrum and the quantum Hall effect on the square lattice with nextnearest-neighbor hopping , Phys. Rev. B 42 , 8282 (1990).
- [45] A. Hill, A. Sinner and K. Ziegler, Valley symmetry breaking and gap tuning in graphene by spin doping , N. J. Ph. 13 , 035023 (2011).
- [46] C.-Y. Hou, C. Chamon, and C. Mudry, Electron fractionalization in two-dimensional graphene like structures , Phys. Rev. Lett. 98 , 186809 (2007).
- [47] P.B. Allen, Electron Transport , in S.G. Louie, M.L. Cohen (Editors), Conceptual Foundations of Materials , Chapter 6, Elsevier B.V. (2006).

## Appendix A: Exclusion of ∆ 2 parameter

Using the ansatz with which

$$M = i M _ { \mu } \Pi _ { \mu } = m \Sigma _ { 3 3 } + \Delta _ { 1 } \Sigma _ { 1 2 } + \Delta _ { 2 } \Sigma _ { 2 2 }.$$

we get to the set of mean-field equations

Projecting both sides on the subspaces of the composite order parameter we get for each α

$$M _ { \alpha } = i g \Pi _ { \alpha } \left [ G _ { 0 } ^ { - 1 } + M \right ] _ { r r } ^ { - 1 }. & & ( A 3 ) \\ \i e \text{ subspaces of the composite order parameter we get for each $\alpha$} \\ & \alpha. &.$$

$$\ m & = \ \frac { g } { 2 } \text{Tr} \Sigma _ { 3 3 } \left [ G _ { 0 } ^ { - 1 } + M \right ] _ { r r } ^ { - 1 }, & ( \text{A4} ) \\ \Delta _ { 1 } & = \ \frac { g } { 2 } \text{Tr} \Sigma _ { 1 2 } \left [ G _ { 0 } ^ { - 1 } + M \right ] _ { r r } ^ { - 1 }, & ( \text{A5} )$$

$$\begin{array} { c c c c c c c } \cdot & 2 & \dashv \cdot \lambda & \lambda r r \cdot & \cdot \cdot \\ \Delta _ { 2 } & = & - \frac { g } { 2 } \text{Tr} \Sigma _ { 2 2 } \left [ G _ { 0 } ^ { - 1 } + M \right ] _ { r r } ^ { - 1 }, & & ( A 6 ) \\ \text{erting the inverse propagator matrix} \end{array}$$

$$\Delta _ { 1 } \, & = \, \begin{array} { c c c c } 2 & \infty \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \\ \frac { g } { 2 } \text{Tr} \Sigma _ { 1 2 } \left [ G _ { 0 } ^ { - 1 } + M \right ] _ { r r } ^ { - 1 }, & \cdot \\ \Delta _ { 2 } \, & = \, - \frac { g } { 2 } \text{Tr} \Sigma _ { 2 2 } \left [ G _ { 0 } ^ { - 1 } + M \right ] _ { r r } ^ { - 1 }, & \cdot \\ \end{array}$$

which becomes after inverting the inverse propagator matrix

$$m \ = \ 2 m g \int \frac { \stackrel { \mu } { d ^ { 3 } Q } } { ( 2 \pi ) ^ { 3 } } \ \frac { \stackrel { \mu } { q _ { 0 } ^ { 2 } + q ^ { 2 } + m ^ { 2 } - \Delta _ { 1 } ^ { 2 } - \Delta _ { 2 } ^ { 2 } } { m ^ { 4 } + 2 m ^ { 2 } ( q _ { 0 } ^ { 2 } + q ^ { 2 } - \Delta _ { 1 } ^ { 2 } - \Delta _ { 2 } ^ { 2 } ) + ( q _ { 0 } ^ { 2 } + q ^ { 2 } + \Delta _ { 1 } ^ { 2 } + \Delta _ { 2 } ^ { 2 } ) ^ { 2 } } } { \stackrel { \mu } { r } \ \stackrel { d ^ { 3 } Q } { r } \ \stackrel { a _ { 0 } ^ { 2 } + q ^ { 2 } - m ^ { 2 } + \Delta _ { 2 } ^ { 2 } + \Delta _ { 2 } ^ { 2 } } { a _ { 0 } ^ { 2 } + q ^ { 2 } - m ^ { 2 } + \Delta _ { 2 } ^ { 2 } + \Delta _ { 2 } ^ { 2 } } }, \quad \ \ ( A 7 ) \\. \quad \ \ \ \. \quad \ r \ \, d ^ { 3 } O \ \ \ \ \ a _ { 0 } ^ { 2 } + q ^ { 2 } - m ^ { 2 } + \Delta _ { 2 } ^ { 2 } + \Delta _ { 2 } ^ { 2 } }$$

$$\Delta _ { 2 } \ = \ - 2 \Delta _ { 2 } g \int \frac { \lambda ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 \prime } } { ( 2 \pi ) ^ { 3 } } \ \frac { \lambda ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \ddots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 2 } \cdots ^ { 3 } } { \quad } \quad \quad \quad$$

$$\Delta _ { 1 } \ = \ 2 \Delta _ { 1 } g \int \frac { \lambda ^ { \prime } } { ( 2 \pi ) ^ { 3 } } \ \frac { \lambda ^ { \prime } _ { 0 } \, \cdots \, \lambda ^ { \prime } _ { 0 } \, \cdots \, \lambda ^ { \prime } _ { 1 } \, \cdots \, \lambda ^ { \prime } _ { 2 } \, \cdots \, \lambda ^ { \prime } _ { 1 } \, \cdots \, \lambda ^ { \prime } _ { 2 } } { m ^ { 4 } + 2 m ^ { 2 } ( q _ { 0 } ^ { 2 } + q ^ { 2 } - \Delta _ { 1 } ^ { 2 } - \Delta _ { 2 } ^ { 2 } ) + ( q _ { 0 } ^ { 2 } + q ^ { 2 } + \Delta _ { 1 } ^ { 2 } + \Delta _ { 2 } ^ { 2 } ) ^ { 2 } } } { \sum _ { 0 } ^ { 2 } + q ^ { 2 } - m ^ { 2 } + \Delta _ { 1 } ^ { 2 } + \Delta _ { 2 } ^ { 2 } }, \quad \ \ ( A 8 ) \\ \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

If there are any mutually independent non-trivial solutions for ∆ 1 and ∆ 2 then they follow from

$$1 \ = \ 2 g \int \frac { d ^ { 3 } Q } { ( 2 \pi ) ^ { 3 } } \ \frac { q _ { 0 } ^ { 2 } + q ^ { 2 } - m ^ { 2 } + \Delta _ { 1 } ^ { 2 } + \Delta _ { 2 } ^ { 2 } } { m ^ { 4 } + 2 m ^ { 2 } ( q _ { 0 } ^ { 2 } + q ^ { 2 } - \Delta _ { 1 } ^ { 2 } - \Delta _ { 2 } ^ { 2 } ) + ( q _ { 0 } ^ { 2 } + q ^ { 2 } + \Delta _ { 1 } ^ { 2 } + \Delta _ { 2 } ^ { 2 } ) ^ { 2 } }, \\ \quad \ f \ d ^ { 3 } Q \quad \ a _ { 0 } ^ { 2 } + a ^ { 2 } - m ^ { 2 } + \Delta _ { 1 } ^ { 2 } + \Delta _ { 2 } ^ { 2 } \quad \dots.$$

$$1 \, = \, - 2 g \int \frac { \lambda ^ { \prime } Q } { ( 2 \pi ) ^ { 3 } } \, \frac { \lambda ^ { \prime } \, \cdots \, \lambda ^ { \prime } \, \cdots \, \lambda ^ { \prime } \, \cdots \, \lambda ^ { \prime } \, \cdots \, \lambda ^ { \prime } \, \cdots \, \lambda ^ { \prime } \, \cdots \, \lambda ^ { \prime } \, \cdots \, \lambda ^ { \prime } \, \cdots \, \lambda ^ { \prime } \, \cdots \, \lambda ^ { \prime } } { m ^ { 4 } + 2 m ^ { 2 } ( q _ { 0 } ^ { 2 } + q ^ { 2 } - \Delta _ { 1 } ^ { 2 } - \Delta _ { 2 } ^ { 2 } ) + ( q _ { 0 } ^ { 2 } + q ^ { 2 } + \Delta _ { 1 } ^ { 2 } + \Delta _ { 2 } ^ { 2 } ) ^ { 2 } }, \quad \quad ( A 1 1 )$$

$$M _ { 1 } = \frac { \Delta _ { 1 } } { 2 } \Sigma _ { 0 1 } - i \frac { \Delta _ { 2 } } { 2 } \Sigma _ { 3 1 } + \frac { m } { 2 } \Sigma _ { 2 0 }, \ M _ { 2 } = \frac { \Delta _ { 1 } } { 2 } \Sigma _ { 3 2 } - i \frac { \Delta _ { 2 } } { 2 } \Sigma _ { 0 2 } - \frac { m } { 2 } \Sigma _ { 1 3 },$$

i. e. only one of them can be true and they cannot be fulfilled together. The two possibilities are

$$\Delta _ { 1 } = 0, \ \Delta _ { 2 } \neq 0 ;$$

$$\Delta _ { 2 } = 0, \ \Delta _ { 1 } \neq 0.$$

/negationslash

/negationslash

The scenario of Eq. (A13) is discussed in the main text. Here we consider the scenario of Eq. (A12), i. e. ∆ = 0 and ∆ 1 2 = ∆ is finite. The system of equations (A7)-(A9) reduces to

$$m \ = \ 2 m g \int \frac { d ^ { 3 } Q } { ( 2 \pi ) ^ { 3 } } \, \frac { q _ { 0 } ^ { 2 } + q ^ { 2 } + m ^ { 2 } - \Delta ^ { 2 } } { m ^ { 4 } + 2 m ^ { 2 } ( q _ { 0 } ^ { 2 } + q ^ { 2 } - \Delta ^ { 2 } ) + ( q _ { 0 } ^ { 2 } + q ^ { 2 } + \Delta ^ { 2 } ) ^ { 2 } }, \\. \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
 \end{cases}$$

$$\Delta \ = \ - 2 \Delta g \int \frac { \ell ^ { 2 } } { ( 2 \pi ) ^ { 3 } } \ \frac { \ell ^ { 2 } _ { 0 } + q ^ { 2 } - m ^ { 2 } + \Delta ^ { 2 } } { m ^ { 4 } + 2 m ^ { 2 } ( q _ { 0 } ^ { 2 } + q ^ { 2 } - \Delta ^ { 2 } ) + ( q _ { 0 } ^ { 2 } + q ^ { 2 } + \Delta ^ { 2 } ) ^ { 2 } }, \\. \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot$$

which upon integrating the frequency and angles becomes

$$m \ & = \ \frac { g } { 8 \pi } \int _ { 0 } ^ { \Lambda ^ { 2 } } \ d q ^ { 2 } \ \left [ \frac { m - \Delta } { \sqrt { q ^ { 2 } + ( m - \Delta ) ^ { 2 } } } + \frac { m + \Delta } { \sqrt { q ^ { 2 } + ( m + \Delta ) ^ { 2 } } } \right ], \\ \Delta \ & = \ \frac { g } { \circ _ { - } } \int _ { 0 } ^ { \Lambda ^ { 2 } } \ d q ^ { 2 } \ \left [ \frac { m - \Delta } { \sqrt { \rho _ { 0 } + 1 } } - \frac { m + \Delta } { \sqrt { \rho _ { 0 } + 1 } } \right ].$$

Summing and subtracting them we get

$$\Delta \ = \ \frac { \mu } { 8 \pi } \int _ { 0 } ^ { \Lambda ^ { 2 } } \, d q ^ { 2 } \, \left [ \frac { m - \Delta } { \sqrt { q ^ { 2 } + ( m - \Delta ) ^ { 2 } } } - \frac { m + \Delta } { \sqrt { q ^ { 2 } + ( m + \Delta ) ^ { 2 } } } \right ]. \quad \text{(A17)} \\ \text{ing and subtracting them we get}$$

$$m + \Delta \ = \ \frac { g } { 4 \pi } \int _ { 0 } ^ { \Lambda ^ { 2 } } \ d q ^ { 2 } \frac { m - \Delta } { \sqrt { q ^ { 2 } + ( m - \Delta ) ^ { 2 } } }, \\ \quad \ a \quad \ r ^ { \Lambda ^ { 2 } } \quad \cdot \quad \ m \pm \Lambda$$

/negationslash

$$m + \Delta \, & = \, \frac { g } { 4 \pi } \int _ { 0 } ^ { \cdot \cdot } \, d q ^ { 2 } \frac { m - \Delta } { \sqrt { q ^ { 2 } + ( m - \Delta ) ^ { 2$$

While the limit ∆ = 0 , m = 0 is well defined, another limit m = 0 , ∆ = 0 is not (left hand side is negative, right hand side is positive) which suggests that the ∆-phase cannot exist alone. It remains to check the case of coexistence of both parameters, which arises as the solution of

/negationslash

$$m \pm \Delta = \gamma ( m \mp \Delta ) \left [ \sqrt { 1 + ( m \mp \Delta ) ^ { 2 } } - | m \mp \Delta | \right ], & & ( A 2 0 ) \\ \text{s after integrating the momenta and rescaling all quantities$$

which appears after integrating the momenta and rescaling all quantities as m → Λ m, ∆ → Λ∆ , γ = g Λ 2 / π . The expression in the squared brackets is positive, i. e. the full equation has any sense only if the sign of m ± ∆ and m ∓ ∆ are the same, which for positive m and ∆ suggests | m &gt; | | ∆ . | Below we assume both sides to be positive . Rewriting it as

$$\frac { m \pm \Delta } { \gamma } + ( m \mp \Delta ) ^ { 2 } = ( m \mp \Delta ) \sqrt { 1 + ( m \mp \Delta ) ^ { 2 } }, & & ( A 2 1 ) \\ \text{oth sides we get}$$

and squaring both sides we get

$$\frac { 1 } { \gamma ^ { 2 } } ( m \pm \Delta ) ^ { 2 } = ( m \mp \Delta ) ^ { 2 } \left [ 1 - \frac { 2 } { \gamma } ( m \pm \Delta ) \right ]. \quad \quad \quad$$

Separating terms which are unique in the sign from those appearing with ± we find

$$\frac { m ^ { 2 } + \Delta ^ { 2 } } { \gamma ^ { 2 } } - \left ( 1 - \frac { 2 m } { \gamma } \right ) \left ( m ^ { 2 } + \Delta ^ { 2 } \right ) - \frac { 4 m \Delta$$

which can only make sense if both sides vanish individually. Left hand side yields

$$\Delta ^ { 2 } & = - m ^ { 2 } \frac { 1 - \gamma ^ { 2 } + 2 \gamma m } { 1 - \gamma ^ { 2 } - 2 \gamma m },$$

while right hand side gives

Equating both right hand sides yields

$$| m | = \frac { 1 - \gamma ^ { 4 } } { 4 \gamma },$$

i. e. it is semi-positive if γ /lessorequalslant 1 and predicts a non-physical phase for small γ . Moreover, plugging it back yields for ∆

$$\Delta ^ { 2 } = - \frac { 1 - \gamma ^ { 4 } } { 1 6 \gamma ^ { 2 } } \left [ \gamma ^ { 4 } + 4 \gamma ^ { 2 } + 3 \right ], & & ( A 2 7 ) \\ \lesssim 1. \text{ But this means that } \Delta \text{ is imaginary which conflicts with the requirement}$$

which is negative for γ /lessorequalslant 1. But this means that ∆ is imaginary which conflicts with the requirement for the effective Hamiltonian to be hermitian. Hence the order parameter ∆ 2 must be excluded from the further consideration.

## Appendix B: Structure of the Proca term

The Gaussian action is defined by Eq. (19) in the main part

$$\mathcal { S } _ { G } [ \mathcal { Q } ] = \vec { \mathcal { Q } } _ { r } \cdot \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \mathbf \Psi } _ { r r ^ { \prime } },$$

where the real vector field /vector Q r = ( /vector Q 1 , /vector Q 2 ) T r

$$\vec { \mathcal { Q } } _ { i } = ( \mathcal { Q } _ { i } ^ { 0 0 }, \mathcal { Q } _ { i } ^ { 0 1 }, \mathcal { Q } _ { i } ^ { 0 2 }, \mathcal { Q } _ { i } ^ { 0 3 }, \mathcal { Q } _ { i } ^ { 1 0 }, \mathcal { Q } _ { i } ^ { 1 1 }, \mathcal { Q } _ { i } ^ { 1 2 }, \mathcal { Q } _ { i } ^ { 1 3 }, \mathcal { Q } _ { i } ^ { 2 0 }, \mathcal { Q } _ { i } ^ { 2 1 }, \mathcal { Q } _ { i } ^ { 2 2 }, \mathcal { Q } _ { i } ^ { 2 3 }, \mathcal { Q } _ { i } ^ { 3 0 }, \mathcal { Q } _ { i } ^ { 3 1 }, \mathcal { Q } _ { i } ^ { 3 2 }, \mathcal { Q } _ { i } ^ { 3 3 } ) \,, \ i = 1, 2. \quad ( B 2 ) \\ \text{The matrix $\mathcal{ }D_{rrr}^{-1}$ is different for each phase.  Below we give Proca matrices for all phases. }$$

The matrix D -1 rr ′ is different for each phase. Below we give Proca matrices for all phases.

## 1. m -phase:

$$1. \ m { -phase: } \\ \Pi ^ { - 1 } ( Q = 0 )$$

Explicitly, the matrix elements are

$$\begin{array} {$$

The difference in comparison to representation in Ref. [29] is explained by the fact that some of the matrix elements given there happen to vanish by integration. The matrix elements fulfill following

$$\Delta ^ { 2 } = - \frac { m } { \gamma } [ 1 + \$$

equalities:

$$M _ { 1 } + M _ { 5 } = \frac { 4 } { \gamma }$$

$$M _ { 2 } + M _ { 4 } = \frac { \cdot } { \gam$$

The eigenvalues of the Gauss matrix for this case are

$$E _ { 1 \cdots 2 8 } ^ { m } = \frac { 2 } { \gamma }, \ \ E _ { 2 9, 3 0 } ^ { m } = \frac { 2 ( 3 + \gamma ^ { 2 } ) } { \gamma + \gamma ^ { 3 } }, \ \ E _ { 3 1, 3 2 } ^ { m } = \frac { 2 ( \gamma ^ { 2 } - 1 ) } { \gamma + \gamma ^ { 3 } }.$$

## 2. ∆-phase:

$$2. \Delta - \text{phase} \colon & \quad \begin{array} { c c c } \Pi ^ { - 1 } ( Q = 0 ) & = & M _ { 1 } ( e _ { 1, 1 } + e _ { 1 0, 1 0 } + e _ { 1 7, 1 7 } + e _ { 2 6, 2 6 } ) + M _ { 2 } ( e _ { 2, 2 } + e _ { 9, 9 } + e _ { 2 4, 2 4 } + e _ { 3 1, 3 1 } ) \\ & & + & M _ { 3 } ( e _ { 3, 3 } + e _ { 1 2, 1 2 } + e _ { 2 1, 2 1 } + e _ { 3 0, 3 0 } ) + M _ { 4 } ( e _ { 4, 4 } + e _ { 1 1, 1 1 } + e _ { 2 0, 2 0 } + e _ { 2 7, 2 7 } ) \\ & & + & M _ { 5 } ( e _ { 5, 5 } + e _ { 1 4, 1 4 } + e _ { 1 9, 1 9 } + e _ { 2 8, 2 8 } ) + M _ { 6 } ( e _ { 6, 6 } + e _ { 1 3, 1 3 } + e _ { 2 2, 2 2 } + e _ { 2 9, 2 9 } ) \\ & & + & M _ { 7 } ( e _ { 7, 7 } + e _ { 1 6, 1 6 } + e _ { 2 3, 2 3 } + e _ { 3 2, 3 2 } ) + M _ { 8 } ( e _ { 8, 8 } + e _ { 1 5, 1 5 } + e _ { 1 8, 1 8 } + e _ { 2 5, 2 5 } ) \\ & & + & i D _ { 1 } ( e _ { 1, 3 2 } - e _ { 7, 2 6 } + e _ { 1 0, 2 3 } - e _ { 1 6, 1 7 } - e _ { 1 7, 1 6 } + e _ { 2 3, 1 0 } - e _ { 2 6, 7 } + e _ { 3 2, 1 } ) \\ & & + & D _ { 2 } ( e _ { 2, 3 1 } - e _ { 9, 2 4 } - e _ { 2 4, 9 } + e _ { 3 1, 2 } ) + D _ { 3 } ( e _ { 3, 3 0 } + e _ { 1 2, 2 1 } + e _ { 2 1, 1 2 } + e _ { 3 0, 3 } ) \\ & & - & i D _ { 4 } ( e _ { 4, 2 9 } + e _ { 6, 2 7 } - e _ { 1 1, 2 2 } - e _ { 1 3, 2 0 } - e _ { 2 0, 1 3 } - e _ { 2 2, 1 1 } + e _ { 2 7, 6 } + e _ { 2 9, 4 } ) \\ & & + & D _ { 5 } ( e _ { 5, 2 8 } + e _ { 1 4, 1 9 } + e _ { 1 9, 1 4 } + e _ { 2 8, 5 } ) \. \\ \text{The matrix elements are} \\ M, = & \frac { 2 } { } + \frac { 1 } { } \frac { 1 } { } \cdots \quad M _ { 2 } = \frac { 2 } { } + 2 \Delta - \frac { 1 + 2 \Delta ^ { 2 } } { } \cdots \quad M _ { 2 } = \frac { 2 } { } - \Delta + \frac { \Delta ^ { 2 } } { } \cdots. \end{array}$$

The matrix elements are

$$& \text{The matrix elements are} \\ & M _ { 1 } = \frac { 2 } { \gamma } + \frac { 1 } { 2 \sqrt { 1 } + \Delta ^ { 2 } }, \quad M _ { 2 } = \frac { 2 } { \gamma } + 2 \Delta - \frac { 1 + 2 \Delta ^ { 2 } } { \sqrt { 1 } + \Delta ^ { 2 } }, \quad M _ { 3 } = \frac { 2 } { \gamma } - \Delta + \frac { \Delta ^ { 2 } } { \sqrt { 1 } + \Delta ^ { 2 } }, \\ & M _ { 4 } = \frac { 2 } { \gamma } - \Delta + \frac { 1 + 2 \Delta ^ { 2 } } { 2 \sqrt { 1 } + \Delta ^ { 2 } }, \quad M _ { 5 } = \frac { 2 } { \gamma } - \Delta + \sqrt { 1 + \Delta ^ { 2 } }, \quad M _ { 6 } = \frac { 2 } { \gamma } + \Delta - \frac { 1 + 2 \Delta ^ { 2 } } { 2 \sqrt { 1 } + 2 }, \\ & M _ { 7 } = \frac { 2 } { \gamma } - \frac { 1 } { 2 \sqrt { 1 } + \Delta ^ { 2 } }, \quad M _ { 8 } = \frac { 2 } { \gamma }, \quad \\ & D _ { 1 } = - \frac { 1 } { 2 \sqrt { 1 } + \Delta ^ { 2 } }, \quad D _ { 2 } = 2 \Delta - \frac { 1 + 2 \Delta ^ { 2 } } { 2 \sqrt { 1 } + \Delta ^ { 2 } }, \quad D _ { 3 } = \Delta - \frac { \Delta ^ { 2 } } { \sqrt { 1 } + \Delta ^ { 2 } }, \\ & D _ { 4 } = - \Delta + \frac { 1 + 2 \Delta ^ { 2 } } { 2 \sqrt { 1 } + \Delta ^ { 2 } }, \quad D _ { 5 } = \sqrt { 1 + \Delta ^ { 2 } } - \Delta, \\ & \text{We notice following most important equalities which are fulfilled by the matrix elements:} \\ & \quad,$$

We notice following most important equalities which are fulfilled by the matrix elements:

$$M _ { 1 } + M _ { 7 } = \frac { 4 } { \gamma }, \quad M _ { 1 } + D _ { 1 } = \frac { 2 } { \gamma } = M _ { 7 } - D _ { 1 }$$

$$M _ { 4 } + M _ { 6 } = \frac { \overset { \prime } { 4 } } { \gamma }, \ \ M _ { 2 } - D _ { 2 } = \frac { \overset { \prime } { 2 } } { \gamma } = M _ { 4 } - D _ { 4 }.$$

$$( B 9 )$$

The stability matrix has the following eigenvalues :

$$E _ { 1 \cdots 2 6 } ^ { \Delta } = \frac { 2 } { \gamma }, \ \ E _ { 2 7, 2 8 } ^ { \Delta } = \frac { 4 } { \gamma }, \ \ E _ { 2 9, 3 0 } ^ { \Delta } = \frac { 4 } { \gamma + \gamma ^ { 3 } }, \ \ E _ { 3 1, 3 2 } ^ { \Delta } = \frac { 2 ( \gamma ^ { 2 } - 1 ) } { \gamma + \gamma ^ { 3 } }. \quad \quad ( B 1 0 )$$

Last two eigenvalues are related to the order parameter of the phase ∆ in Eq. (17) in the main part.

- 3. Coexisting phase m = ∆:

$$3. \, \text{Coexisting phase } m = \Delta \colon \\ \Pi ^ { - 1 } ( Q = 0 ) \, = \, M _ { 1 } \left ( e _ { 1, 1 } + e _ { 1 0, 1 0 } + e _ { 1 7, 1 7 } + e _ { 2 6, 2 6 } \right ) + M _ { 2 } \left ( e _ { 2, 2 } + e _ { 9, 9 } + e _ { 2 4, 2 4 } + e _ { 3 1, 3 1 } \right ) \\ & + \, M _ { 3 } \left ( e _ { 3, 3 } + e _ { 1 2, 1 2 } + e _ { 2 1, 2 1 } + e _ { 3 0, 3 0 } \right ) + M _ { 4 } \left ( e _ { 4, 4 } + e _ { 1 1, 1 1 } + e _ { 2 0, 2 0 } + e _ { 2 7, 2 7 } \right ) \\ & + \, M _ { 5 } \left ( e _ { 5, 5 } + e _ { 1 4, 1 4 } + e _ { 1 9, 1 9 } + e _ { 2 8, 2 8 } \right ) + M _ { 6 } \left ( e _ { 6, 6 } + e _ { 1 3, 1 3 } + e _ { 2 2, 2 2 } + e _ { 2 9, 2 9 } \right ) \\ & + \, M _ { 7 } \left ( e _ { 7, 7 } + e _ { 1 6, 1 6 } + e _ { 2 3, 2 3 } + e _ { 3 2, 3 2 } \right ) + M _ { 8 } \left ( e _ { 8, 8 } + e _ { 1 5, 1 5 } + e _ { 1 8, 1 8 } + e _ { 2 5, 2 5 } \right ) \\ & + \, A _ { 1 } \left ( e _ { 1, 1 0 } - e _ { 7, 1 6 } + e _ { 1 0, 1 } - e _ { 1 6, 7 } - e _ { 1 e, 1 2 } \right ) + e _ { 7, 1 7 } - e _ { 1 0, 3 2 } + e _ { 1 6, 2 6 } \\ & + \, e _ { 1 7, 7 } - e _ { 2 3, 1 } + e _ { 2 6, 1 6 } - e _ { 3 2, 1 0 } + e _ { 1 7, 2 6 } - e _ { 2 3, 3 2 } + e _ { 2 6, 1 7 } - e _ { 3 2, 2 3 } \right ) \\ & + \, A _ { 2 } \left ( e _ { 2, 9 } + e _ { 9, 2 } - e _ { 2, 2 4 } + e _ { 9, 3 1 } - e _ { 2 4, 2 } + e _ { 3 1, 9 } - e _ { 2 4, 3 1 } - e _ { 3 1, 2 4 } \right ) \\ & + \, A _ { 3 } \left ( e _ { 3, 3 0 } + e _ { 1 2, 2 1 } + e _ { 2 1, 1 2 } + e _ { 3 0, 3 } \right ) \\ & + \, i D _ { 1 } \left ( e _ { 1, 3 2 } - e _ { 7, 2 6 } + e _ { 1 0, 2 3 } - e _ { 1 6, 1 7 } - e _ { 1 7, 1 6 } + e _ { 2 3, 1 0 } - e _ { 2 6, 7 } + e _ { 3 2, 1 } \right ) \\ & + \, D _ { 2 } \left ( e _ { 2, 3 1 } - e _ { 9, 2 4 } - e _ { 2 4, 9 } + e _ { 3 1, 2 } \right ) \\ & - \, i D _ { 3 } \left ( e _ { 4, 2 9 } + e _ { 6, 2 7 } - e _ { 1 1, 2 2 } - e _ { 1 3, 2 0 } - e _ { 2 0, 1 3 } - e _ { 2 2, 1 1 } + e _ { 2 7, 6 } + e _ { 2 9, 4 } \right ) \\ & + \, D _ { 4 } \left ( e _ { 5, 2 8 } + e _ { 1 4, 1 9 } + e _ { 1 9, 1 4 } + e _ { 2 8, 5 } \right ) \,. \\ \text{Each quantity in the Proca matrix means} \\ & A _ { 1 } \, = \, \frac { 1 - \sqrt { 1 + 4 m ^ { 2 } } } { 4, \frac { 1 + 4 m ^ { 2 } } { 4 m ^ { 2 } } } \,,$$

Each quantity in the Proca matrix means

$$\text{tity in the Proca matrix means} \\ A _ { 1 } \, & = \, \frac { 1 - \sqrt { 1 + 4 m ^ { 2 } } } { 4 \sqrt { 1 + 4 m ^ { 2 } } }, \\ A _ { 2 } \, & = \, \frac { 1 } { 2 } \left [ 1 + 4 m - 2 \sqrt { 1 + 4 m ^ { 2 } } + \frac { 1 } { \sqrt { 1 + 4 m ^ { 2 } } } \right ], \\ A _ { 3 } \, & = \, - \frac { 1 } { 6 m ^ { 2 } } \left [ 1 - \sqrt { 1 + 4 m ^ { 2 } } + 2 m ^ { 2 } \left ( \sqrt { 1 + 4 m ^ { 2 } } - 2 m \right ) \right ]. \\ \text{ments } A _ { 1 \dots 3 } \text{ appear only as consequence of coexisting order parameters } m \text{ and}$$

Matrix elements A 1 ··· 3 appear only as consequence of coexisting order parameters m and ∆.

$$\begin{array} { c c c c c c } & & \cdots & \cdots & \cdots & \cdots \\ \colon \text{elements } A _ { 1 } \cdots _ { 3 } \text{ appear only as consequence of coexisting order parameters } m \text{ and } \Delta \\ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &  & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & \cdots & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & && & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & \\ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &   & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &  \cdots & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &  & &  & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &  & &
 & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &
 & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &  & & \\ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &  & & \ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & - & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & _ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &. & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & ; & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & } & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &, & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &$$

where (and further on) m is defined in Eq. (17) in the main part.

$$D _ { 1 } \ & = \ - \frac { 1 + \sqrt { 1 + 4 m ^ { 2 } } } { 4 \sqrt { 1 + 4 m ^ { 2 } } }, \\ D _ { 2 } \ & = \ \frac { 1 } { 2 } \left [ 4 m - 1 - 2 \sqrt { 1 + 4 m ^ { 2 } } + \frac { 1 } { \sqrt { 1 + 4 m ^ { 2 } } } \right ], \\ D _ { 3 } \ & = \ - \frac { 1 } { 1 2 m ^ { 2 } } \left [ 1 - \sqrt { 1 + 4 m ^ { 2 } } - 4 m ^ { 2 } \left ( \sqrt { 1 + 4 m ^ { 2 } } - 2 m \right ) \right ], \\ D _ { 4 } \ & = \ \sqrt { 1 + 4 m ^ { 2 } } - 2 m. \\ \text{ility matrix has following eigenvalues} \\ 2 \quad 4 \quad 8 \quad 1 + 2 \varepsilon \quad 2 ( \varepsilon ^ { 2 } - 1 )$$

The stability matrix has following eigenvalues

$$E _ { 1 \cdots 2 6 } = \frac { 2 } { \gamma }, \ \ E _ { 2 7, 2 8 } = \frac { 4 } { \gamma }, \ \ E _ { 2 9, 3 0 } = \frac { 8 } { 3 \gamma } \frac { 1 + 2 \gamma } { ( 1 + \gamma ) ^ { 2 } }, \ \ E _ { 3 1 } = \frac { 2 ( \gamma ^ { 2 } - 1 ) } { \gamma + \gamma ^ { 3 } } \quad \quad ( B 1 1 )$$

the latter related to the order parameter defined in Eq. (17) in the main text. However, the last remaining eigenvalue

$$E _ { 3 2 } = - 2 + \frac { 2 } { \gamma }$$

is negative for all γ &gt; 1.

## Appendix C: Evaluation of the current-current correlator

Here we evaluate Eq. (61) for the case of usual electronic currents. This correlator is zero in the ∆-phase but finite in the m -phase, i.e. we have to chose c = d = 1 and i = 1 , j = 2

$$\frac { 1 } { 4 \gamma ^ { 2 } } \left ( \Lambda _ { 1, r } ^ { 1 } \Lambda _ { 2, r ^ { \prime } } ^ { 1 } \right ) = \frac { 4 } { \gamma ^ { 2 } } \left ( i \mathcal { Q } _ { 1 r } ^ { 1 2 } \mathcal { Q } _ { 1 r ^ { \prime } } ^ { 2 1 } - i \mathcal { Q } _ { 2 r } ^ { 2 1 } \mathcal { Q } _ { 2 r ^ { \prime } } ^ { 1 2 } - \mathcal { Q } _ { 1 r } ^ { 1 2 } \mathcal { Q } _ { 2 r ^ { \prime } } ^ { 1 2 } - \mathcal { Q } _ { 2 r } ^ { 2 1 } \mathcal { Q } _ { 1 r ^ { \prime } } ^ { 2 1 } \right ). \\ \colon & \text{forming the functional integration we get}$$

Performing the functional integration we get

$$\frac { 1 } { 4 \gamma ^ { 2 } } \langle \Lambda ^ { 1 } _ { 1, r } \Lambda ^ { 1 } _ { 2, r ^ { \prime } } \rangle _ { G } = \frac { 4 } { \gamma ^ { 2 } } \int \frac { d ^ { 3 } Q } { ( 2 \pi ) ^ { 2 } } \ e ^ { i Q ( r - r ^ { \prime } ) } \frac { 1 } { 2 } \left [ \Pi _ { 2 6, 1 0 } ( q _ { 0 } ) + \Pi _ { 7, 2 3 } ( q _ { 0 } ) + i \Pi _ { 2 6, 2 3 } ( q _ { 0 } ) - i \Pi _ { 7, 1 0 } ( q _ { 0 } ) \right ], \ \, \\ \text{where the factor $1/2$ appears because of the Gauss integral. and the propagator is defined as inverse}$$

where the factor 1 / 2 appears because of the Gauss integral, and the propagator is defined as inverse kernel matrix from Eq. (42)

$$\mathbf P ( Q ) = \left ( \begin{array} { c c } \mathbf A & \mathbf B \\ C & \mathbf D \end{array} \right ) _ { Q } ^ { - 1 }.$$

In m -phase we get to linear order in q 0 via a direct evaluation

$$\frac { 2 } { \gamma ^ { 2 } } \left [ \Pi _ { 2 6, 1 0 } ( q _ { 0 } ) + \Pi _ { 7, 2 3 } ( q _ { 0 } ) + i \Pi _ { 2 6, 2 3 } ( q _ { 0 } ) - i \Pi _ { 7, 1 0 } ( q _ { 0 } ) \right ] & \ \approx \ - \text{sgn} ( m ) \frac { 2 } { \gamma ^ { 2 } } \frac { q _ { 0 } } { 4 \pi } \frac { ( M _ { 1 } + M _ { 5 } ) ^ { 2 } } { ( M _ { 1 } M _ { 2 } + \Delta ^ { 2 } ) ^ { 2 } } \\ & = \ - \text{sgn} ( m ) \frac { q _ { 0 } } { 2 \pi },$$

where the equalities Eq. (B4) were used. Hence we finally obtain for the correlator

$$\frac { 1 } { 4 \gamma ^ { 2 } } \langle \Lambda _ { 1, r } ^ { 1 } \Lambda _ { 2, r ^ { \prime } } ^ { 1 } \rangle _ { \mathcal { G } } \approx - \frac { \text{sgn} ( m ) } { 2 \pi } \int \frac { d ^ { 3 } Q } { ( 2 \pi ) ^ { 2 } } \, e ^ { i Q ( r - r ^ { \prime } ) } \ q _ { 0 } = \frac { \text{sgn} ( m ) } { 2 \pi } \ i \partial _ { \tau } \delta ( r - r ^ { \prime } ). \quad \quad ( C 6 ) \\ \sim \dots \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \colon \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \dots$$

Calculating correlators for all index combinations and bringing them into a compact form we get Eq. (61).
Andreas Sinner Klaus Ziegler
Quantum Hall effect induced by electron-phonon interaction
When phonons couple to fermions in 2D semimetals, the interaction may turn
the system into an insulator. There are several insulating phases in which the
time reversal and the sublattice symmetries are spontaneously broken. Examples
are many-body states commensurate to Haldane's staggered flux model or to
lattice models with periodically modulated strain. We find that the effective
field theories of these phases exhibit characteristic Chern-Simons terms, whose
coefficients are related to the topological invariants of the microscopic
model. This implies that the corresponding quantized Hall conductivities
characterize these insulating states.