---
abstract: 'In recent theoretical and experimental investigations, researchers have linked the low-energy field theory of a Weyl semimetal gapped with a charge-density wave (CDW) to high-energy theories with axion electrodynamics. However, it remains an open question whether a lattice regularization of the dynamical Weyl-CDW is in fact a single-particle axion insulator (AXI). In this Letter, we use analytic and numerical methods to study both lattice-commensurate and incommensurate minimal (magnetic) Weyl-CDW phases in the mean-field state. We observe that, as previously predicted from field theory, the two inversion- ($\mathcal{I}$-) symmetric Weyl-CDWs with $\phi = 0,\pi$ differ by a topological axion angle $\delta\theta_{\phi}=\pi$. However, we crucially discover that *neither* of the minimal Weyl-CDW phases at $\phi=0,\pi$ is individually an AXI; they are instead quantum anomalous Hall (QAH) and “obstructed” QAH insulators that differ by a fractional translation in the modulated cell, analogous to the two phases of the Su-Schrieffer-Heeger model of polyacetylene. Using symmetry indicators of band topology and non-abelian Berry phase, we demonstrate that our results generalize to multi-band systems with only two Weyl fermions, establishing that minimal Weyl-CDWs unavoidably carry nontrivial Chern numbers that prevent the observation of a static magnetoelectric response. We discuss the experimental implications of our findings, and briefly outline generalizations to nonmagnetic semimetals with additional Weyl points.'
author:
- 'Benjamin J. Wieder'
- 'Kuan-Sen Lin'
- Barry Bradlyn
bibliography:
- 'refs.bib'
title: 'Is the Dynamical Axion Weyl-Charge-Density Wave an Axionic Band Insulator?'
---

In condensed matter physics, one of the most important tools is low-energy field theory. From the ${\bf k}\cdot {\bf p}$ Hamiltonian of a solid-state material, an effective action can be developed to characterize the robust, and frequently topological, long-wavelength response effects [@wang2013chiral; @you2016response; @TitusCDW; @qi2008topological; @essin2009magnetoelectric; @wilczekaxion; @zhang1989effective; @halperin1993theory; @son2015composite]. However, to extrapolate from a low-energy field theory to an experimentally observable response effect, one must carefully complete the theory to short (UV) wavelengths – specifically, two field theories that are identical at the ${\bf k}\cdot {\bf p}$ level may differ at large momenta, leading to distinct physical interpretations. For example, the ${\bf k}\cdot{\bf p}$ Hamiltonian $\mathcal{H}({\bf q}) = \sigma^{x}q_{x} + \sigma^{y}q_{y}$ can characterize one of the twofold Dirac points in a graphene-like 2D semimetal [@MeleGraphene; @SemenoffGraphene; @GrapheneReview], or the isolated Dirac point on the surface of a time-reversal- ($\mathcal{T}$-) symmetric 3D topological insulator (TI) [@fukanemele; @Fu2007; @Zahid3DTI; @xia2009observation]. In both cases, $\mathcal{H}({\bf q})$ carries the response of a half-level $2+1$-D Chern-Simons theory [@qi2008topological; @essin2009magnetoelectric; @fukanemele; @Fu2007; @LapaParity]. For a 2D semimetal, the UV completion must account for additional Dirac fermions at higher momenta, which are required be present by the parity anomaly [@LapaParity; @WittenParity; @RedlichParity; @JackiwParity; @haldanemodel; @MulliganAXI; @DiracInsulator], leading to an overall integer-valued Hall conductivity [@qi2008topological; @haldanemodel; @bernevigbook; @hasan2010colloquium; @fradkin1]. However, in the case of a 3D TI, the Dirac fermion is unpaired, and a $\mathcal{T}$-preserving UV completion is impossible without accounting for the bulk of the system, which conversely contributes an anomalous *half-integer* Hall conductivity to a magnetically gapped TI surface [@fukanemele; @Fu2007; @qi2008topological; @LapaParity; @Hughes11; @DiracInsulator; @wieder2018axion; @HOTIFan; @vanderbiltaxion; @essin2009magnetoelectric; @MulliganAXI; @witten2016fermion].

Some of the most intriguing low-energy field theories involve condensed matter realizations of high-energy electrodynamics [@wilczekaxion]. Specifically, in 3D insulators, the long-wavelength linear response can be formulated as an effective action[^1]: $$S[A_\mu]=\frac{1}{4\pi^2}\int d^4x \epsilon^{\mu\nu\lambda\rho}(\theta\partial_\mu+ v_\mu)A_\nu\partial_\lambda A_\rho,
\label{eq:mainAction}$$ in which $\mathbf{A_\mu}$ is the electromagnetic gauge potential [@qi2008topological; @essin2009magnetoelectric; @you2016response]. In Eq. (\[eq:mainAction\]), ${\bf v}$ is a rotational-symmetry-breaking vector that determines the quantum anomalous Hall (QAH) response; when $v_\mu=(\mathbf{v},0)$ is constant vector: $$\sigma^H_{ij}=\frac{e^2v_k}{h}\epsilon_{kij}$$ where $\sigma^{H}_{ij}$ is the Hall conductivity tensor. For gapped periodic systems: $$\mathbf{v} = C(\mathbf{k}\cdot \mathbf{R}_i)\mathbf{G}_i,$$ where $\mathbf{R}_i$ is a primitive lattice vector, $\mathbf{G}_i$ is a primitive reciprical lattice vector, and $C(\mathbf{k}\cdot \mathbf{R}_i)\equiv v_i/|\mathbf{G}_i|=\nu_i$ is the *weak* Chern number \[Fig.\[fig:scheme\](b)\], which can be evaluated as the integral of the abelian Berry curvature over a plane in the Brillouin zone (BZ) normal to $\mathbf{R}_i$ [^2] [@halperin1987possible; @kohmoto1992diophantine; @qi2008topological; @haldane2004berry]. In Eq. (\[eq:mainAction\]), $\theta$ is known as the axion angle, and governs magnetoelectric response [@qi2008topological; @essin2009magnetoelectric]; in terms of the non-abelian Berry connection ${\bf \mathcal{A}}$: $$\theta[{\bf\mathcal{A}}]=\frac{1}{4\pi}\int d^3k\epsilon^{ijk}\mathrm{tr}\left(\mathcal{A}_i\partial_j\mathcal{A}_k-\frac{2i}{3}\mathcal{A}_i\mathcal{A}_j\mathcal{A}_k\right).
\label{eq:thetadef}$$ Because $\theta$ is the coefficient of ${\bf E}\cdot {\bf B}$ in Eq. (\[eq:mainAction\]), then spatial inversion ($\mathcal{I}$) and $\mathcal{T}$ symmetries (as well as other, more complicated crystal and magnetic symmetries) [@wieder2018axion] act to quantize $\theta$ as a $\mathbb{Z}_{2}$ topological invariant for which $\theta=0$ ($\theta=\pi$) is the trivial (topological) value [@qi2008topological; @essin2009magnetoelectric; @witten2016fermion][^3]. In particular, $\nu_{x,y,z}=0$, $\theta=\pi$ is one definition of a 3D TI when $\mathcal{T}$ symmetry quantizes $\theta$, and defines a magnetic axion insulator (AXI) when $\mathcal{T}$ symmetry is absent and $\theta$ is instead quantized by $\mathcal{I}$ [@qi2008topological; @essin2009magnetoelectric; @wieder2018axion; @Hughes11; @Turner2012; @AFMTIMoore]. In 3D TIs and AXIs, the combination of $\theta=\pi$ and $\nu_{x,y,z}=0$ leads to several unusual response properties, such as low-energy excitations resembling magnetic monopoles (the Witten effect [@WittenEffect; @rosenberg2010witten]) and quantized Faraday and Kerr rotations [@qi2008topological; @wu2016quantized]. AXIs have recently been recognized to be a form of “higher-order” TI [@wieder2018axion; @xu2020high; @MTQC; @HOTIFan; @vanderbiltaxion; @benalcazar2017quantized; @benalcazar2017electric; @Po2017; @langbehn2017reflection; @watanabe2018structure; @khalaf; @khalaf2018higher; @KoreanAXI; @NicoDavidAXI2; @TMDHOTI; @fang2019new; @JiabinZhidaAXIDirac; @MurakamiAXI1; @MurakamiAXI2; @pozo2019quantization] (HOTI) with gapped surfaces and gapless hinges, where, in particular, $\mathcal{I}$-symmetric, finite-sized AXIs host an odd number of sample-encircling chiral hinge modes \[Fig. \[fig:scheme\](a)\]. AXI phases have recently been theoretically proposed in a number of compounds [@xu2019higher; @xu2020high; @SarmaMobiusAXI] and observed in Mn-doped 3D TIs [@AXIExp1; @AXIExp2; @AXIExp3; @SlagerAshvinAXI].

![(a) An inversion- ($\mathcal{I}$-) symmetric axion insulator (AXI) [@qi2008topological; @essin2009magnetoelectric; @wieder2018axion; @Hughes11; @Turner2012; @AFMTIMoore] with $\nu_{x,y,z}=0,\ \theta=\pi$ in Eq. (\[eq:mainAction\]). (b) A weak, $\hat{\bf z}$-directed, $\mathcal{I}$-symmetric stack of Chern insulators with a nontrivial quantum anomalous Hall (QAH) effect [@halperin1987possible; @kohmoto1992diophantine] \[$\nu_{z}=-1,\ \nu_{x}=\nu_{y}=\theta=0$\]. (c) A superposition of the insulators in (a) and (b) is can be deformed into an “obstructed” QAH (oQAH) insulator \[$\nu_{z}=-1,\ \nu_{x}=\nu_{y}=0,\ \theta=\pi$\] equivalent to shifting the QAH insulator in (b) by a half-lattice translation. A finite-sized oQAH insulator either exhibits coexisting surface and hinge states or exactly one fewer (or one more) QAH surface state, depending on whether $\mathcal{I}$ symmetry is weakly broken.[]{data-label="fig:scheme"}](AXIoQAH.pdf){width="\columnwidth"}

In several recent proposals [@wang2013chiral; @you2016response; @zyuzin2012weyl; @maciejko2014weyl], it was also recognized that the low-energy field theory of a bulk-gapless topological material could realize a gapped $\theta=\pi$ phase in the presence of a correlated charge-density wave (CDW). Specifically, in Weyl semimetals (WSMs) – whose bulk Fermi pockets \[Weyl points (WPs)\] are sources and sinks of Berry curvature characterized by integer-valued topological (chiral) charges [@Wan11; @BalentsWeyl; @armitage2018weyl] \[Fig. \[fig:zonefold\](a)\] – it was shown at the $\mathbf{k}\cdot\mathbf{p}$ level that: $$\theta({\bf x},t)=\theta_0+\phi({\bf x},t),
\label{eq:CDWthetaphi}$$ where $\theta_0=\mathbf{Q}\cdot\mathbf{x}$ is the contribution to $\theta({\bf x},t)$ from two WPs separated by a momentum $\mathbf{Q}$, and $\phi({\bf x},t)$ is the (dynamical) phase of the CDW order parameter [@wang2013chiral; @you2016response; @zyuzin2012weyl; @maciejko2014weyl]. The appearance of an axionic response was attributed to the chiral anomaly in quantum field theory. In the mean-field approximation, $\phi({\bf x},t)$ is a Goldstone mode, and hence can be tuned freely; uniform shifts of $\phi({\bf x},t)$ are the current-carrying sliding mode of the CDW. In disordered or incommensurate CDWs, however, $\phi({\bf x},t)$ is typically pinned to a non-universal value [@fukuyama1978dynamics; @bardeen1980tunneling; @gruner1988dynamics]. The proposals in [@wang2013chiral; @you2016response; @zyuzin2012weyl; @maciejko2014weyl] have recently been revisited in light of the discovery of high-temperature WPs in the archetypal CDW compound (TaSe$_4$)$_2$I [@shi2019charge; @OtherARPESCDWWeyl; @Ta2Se8IPrepare; @Cava1986], where nonlinear negative magnetoresistance consistent with a gapped dynamical Weyl-CDW phase was observed in (TaSe$_4$)$_2$I below the CDW transition temperature [@gooth2019evidence].

![Band structure as a function of $k_z$ demonstrating commensurate BZ folding from a CDW with ${\bf Q}=(2\pi/Nc)\hat{\bf z}$, $N=3$. (a) A minimal $\mathcal{I}$-symmetric semimetal with two Weyl points (WPs) at $k_{z} = \pm \pi/3c$ and a BZ with ${\bf G}_{z} = (2\pi/c)\hat{\bf z}$ \[Eqs. (\[eq:positionTB\]) and (\[eq:ham0\])\]. Solid blue (red dashed) \[black dotted\] bands lie in the first (second) \[third\] reduced BZ (rBZ). (b) The rBZ, for which ${\bf G}_{z}' = {\bf Q}$. The WPs from (a) are folded in (b) into an (unstable) fourfold degeneracy at the rBZ boundary.[]{data-label="fig:zonefold"}](zonefold-final.pdf){width="0.7\columnwidth"}

Confusingly, the gapped dynamical Weyl-CDW phase is frequently labeled an AXI in the literature. However, because Eq. (\[eq:CDWthetaphi\]) was derived from a ${\bf k}\cdot {\bf p}$ approximation, the relationship between the field-theory analysis and mean-field band topology should not be taken for granted. There have remained open and urgent questions whether there exists a UV completion in which $\delta\theta_{\phi} = \theta_{\phi=\pi}-\theta_{\phi=0}\text{ mod }2\pi = \pi$ emerges due to the topology of band electrons, and whether the bulk at $\phi=0,\pi$ is a single-particle electronic AXI. In this Letter, we use analytic and numerical methods to demonstrate that the UV completion of the simplest dynamical Weyl-CDW phase is in fact *not* an AXI, but is simply an $\mathcal{I}$-symmetric QAH insulator. However, we discover that, depending on the CDW angle $\phi$, the Weyl-CDW phase can be decomposed into two *topologically distinct* QAH phases – a QAH insulator or an “obstructed” QAH (oQAH) insulator \[Fig. \[fig:scheme\](b,c)\] – that differ by fractional lattice translation. Crucially, because $\theta$ is origin-dependent in the presence of a background QAH [@qi2008topological; @vanderbiltaxion; @NicoDavidAXI2], then we find that the QAH and oQAH phases still *differ* by an *origin- (gauge-) independent, topological* axion angle $\delta\theta_{\phi}=\pi$ that reflects a difference in $\mathcal{I}$-quantized “Chern number polarization” [@NicoDavidAXI2; @zilberberg2018photonic]. This provides a direct analogy between axionic CDWs and the Su-Schrieffer-Heeger model of polyacetylene [@ssh1979], in which both phases are trivial atomic limits that differ by a fractional lattice translation corresponding to an $\mathcal{I}$-quantized topological polarization. We demonstrate that the relative axionic response of two $\mathcal{I}$-symmetric Weyl-CDWs can be attributed to the single-particle quasiparticle band structure. Additionally, we generalize our findings to multi-band systems with two WPs and to incommensurate CDWs, establishing that QAH insulators and topological phase shifts $\delta\theta_{\phi}=\pi$ are generic in minimal $\mathcal{I}$-symmetric Weyl-CDWs. We conclude by discussing the experimental implications of our findings in the larger search for AXIs, and briefly outline generalizations to $\mathcal{T}$-symmetric WSMs with additional WPs.

#### Model:

We begin by introducing a simple model of a minimal $\mathcal{T}$-broken, $\mathcal{I}$-symmetric magnetic WSM with two WPs and with orthorhombic lattice vectors of length $a,b,c$ in the $\mathbf{\hat{x}},\mathbf{\hat{y}},\mathbf{\hat{z}}$ directions, respectively[@mccormick2017minimal]: $$\begin{aligned}
    H_0&=\left(\sum_{\mathbf{R}}\left[it_xc^\dag_R\sigma^xc_{\mathbf{R}+\mathbf{\hat{x}}}+it_yc^\dag_{\mathbf{R}}\sigma^yc_{\mathbf{R}+\mathbf{\hat{y}}}+t_zc^\dag_{\mathbf{R}}\sigma^z c_{\mathbf{R}+\mathbf{\hat{z}}}\right]\right. \nonumber \\
    &+\sum_{\mathbf{R}}\frac{m}{2}\left(c^\dag_\mathbf{R}\sigma^zc_{\mathbf{R}+\mathbf{\hat{x}}} + c^\dag_\mathbf{R}\sigma^zc_{\mathbf{R}+\mathbf{\hat{y}}}  -2 c^\dag_\mathbf{R}\sigma^zc_{\mathbf{R}} \right) \nonumber \\
    &\left.-\sum_\mathbf{R} t_z\cos\frac{Qc}{2}c^\dag_\mathbf{R}\sigma^z c_\mathbf{R}\right) +\mathrm{h.c.},
\label{eq:positionTB}\end{aligned}$$ where $m/2>t_x,t_y,t_z>0$. While opaque in position space, Eq. (\[eq:positionTB\]) assumes a simpler form in ${\bf k}$-space: $$\begin{aligned}
    H(\mathbf{k})&= -2[t_x\sigma^x\sin (k_xa) +t_y\sigma^y\sin (k_yb)]\nonumber \\
    & +2t_z\sigma^{z}[\cos (k_zc) -\cos\frac{Qc}{2}] \nonumber \\
    &-m\sigma^{z}[2-\cos (k_xa) - \cos (k_yb)].
\label{eq:ham0}\end{aligned}$$ Eq. (\[eq:ham0\]) is gapped at half filling at all ${\bf k}$ points away from two WPs at $\mathbf{k}=(0,0,\pm Q/2$) with chiral charges $|C|=1$ \[Fig. \[fig:zonefold\](a)\]. The WPs in Eq. (\[eq:ham0\]) are related by $\mathcal{I}$, which is represented by: $$\mathcal{I}H(\mathbf{k})\mathcal{I}^{-1} = \sigma^z H(-\mathbf{k})\sigma^z.
\label{eq:idef}$$ As shown in [@QWZ; @Fang2012] and in the Supplementary Material (SM), the occupied parity ($\mathcal{I}$) eigenvalues imply the ${\bf k}$-space Chern numbers $C(ck_{z}=0)\text{ mod }2=-1$, $C(ck_{z}=\pi)\text{ mod }2=0$, mandating the appearance of the $|C|=1$ WPs.

We next construct a $\mathbf{k}\cdot\mathbf{p}$ expansion of Eq. (\[eq:ham0\]) about the two WPs: $$H(\mathbf{q})\approx -(2t_{x}aq_x\sigma^x+2t_{y}bq_{y}\sigma^y)\tau^0 +2t_{z}cq_{z}\sin\frac{Qc}{2}\sigma^z\tau^z, 
\label{eq:kdp}$$ where the Pauli matrices $\vec{\tau}$ act in the space of $c_{1/2,{\bf k}}$ in: $$c_\mathbf{R}\approx\sum_\mathbf{k} c_{1\mathbf{k}}e^{i{\bf R}\cdot[(Q/2)\hat{\bf z}+{\bf k}]}+c_{2\mathbf{k}}e^{-i{\bf R}\cdot[(Q/2)\hat{\bf z}-{\bf k}]}.
\label{eq:expandedop}$$ Eq. (\[eq:kdp\]) can be gapped by a CDW distortion: $$\begin{aligned}
H_{CDW}&=2\sum_{\mathbf{R}}|\Delta|\cos(QR_z+\phi)c^\dag_\mathbf{R}\sigma^z c_\mathbf{R} \\
    &=\Delta\sum_\mathbf{k} (c^\dag_{\mathbf{k}-\frac{Q}{2}\mathbf{\hat{z}}}\sigma^z c_{\mathbf{k}+\frac{Q}{2}\mathbf{\hat{z}}}e^{-i\phi} + \mathrm{h.c.}),
\label{eq:cdwmass}\end{aligned}$$ which breaks the translation symmetry of Eqs. (\[eq:positionTB\]) and (\[eq:ham0\]), coupling the two WPs. At the ${\bf k}\cdot {\bf p}$ level, Eq. (\[eq:cdwmass\]) induces a mass in Eq. (\[eq:kdp\]): $$V_\phi = |\Delta|\sigma^z(\tau^{x}\cos\phi + \tau^{y}\sin\phi).
\label{eq:kdpMass}$$ Crucially, $\mathcal{I}$ symmetry is now represented in Eqs. (\[eq:kdp\]) and (\[eq:kdpMass\]) by: $$\mathcal{I}H(\mathbf{q})\mathcal{I}^{-1} = \sigma^z \tau^x H(-\mathbf{q})\sigma^z\tau^x,
\label{eq:idefKDP}$$ such that Eq. (\[eq:kdpMass\]) only preserves $\mathcal{I}$ (centered at the origin) for $\phi=0,\pi$. For $Q\neq \pi/c$ and $|\Delta|>0$, this model is gapped for all $\phi$.

Consistent with previous works [@wang2013chiral; @you2016response; @zyuzin2012weyl; @maciejko2014weyl], we recognize that a domain wall between $\phi=0,\pi$ is equivalent to the critical point between a trivial insulator and an AXI [@wieder2018axion; @qi2008topological]. Correspondingly, because $\{H({\bf q}),V_{\phi}\}=0$ for all $\phi$, $\mathcal{I}$-breaking defects in the space $(\Delta,\phi)$ will bind chiral modes [@teo2010topological] (the axion strings of Ref. ). In [@wilczekaxion; @wang2013chiral; @you2016response; @ryu2010topological; @qi2008topological], the authors used the chiral anomaly to motivate the appearance of chiral modes, specifically identifying the relationship $\theta=[\pi/2](1-\mathrm{sgn}[\cos\phi])\text{ mod }2\pi$ for $\phi=0,\pi$. The change in axion angle $\delta \theta_{\phi} = \delta \phi$ is also consistent with magnetic symmetry-based indicators $\{\tilde{z}_4|\tilde{z}_{2x}\tilde{z}_{2y}\tilde{z}_{2z}\}$ for 3D crystals with only $\mathcal{I}$ and translation symmetries [@xu2020high; @MTQC; @Po2017; @watanabe2018structure; @khalaf; @khalaf2018higher; @wieder2018axion; @JiabinZhidaAXIDirac; @MurakamiAXI1; @MurakamiAXI2]: $$\begin{aligned}
\tilde{z}_4&=\frac{1}{2}\sum_{\mathbf{k}_a\in\mathrm{TRIMS}} (n^a_+- n^a_-)\text{ mod }4, \nonumber \\
\tilde{z}_{2,i} &= \frac{1}{2}\sum_{{\bf k_a\cdot R_{i}}=\pi\in\mathrm{TRIMS}}(n^a_+- n^a_-)\text{ mod }2,
\label{eq:z4def}\end{aligned}$$ where $n^a_{\pm}$ are the number of states with $\pm 1$ parity eigenvalues at $\mathbf{k}_a$. Specifically, $\delta \phi = \pi$ in Eq. (\[eq:kdpMass\]) implies that $|\delta \tilde{z}_{4}| = 2$. However, because weak Chern numbers are only $\mathcal{I}$-symmetry-indicated modulo $2$ [@QWZ; @Fang2012], $|\delta\tilde{z}_4|=2$ does not in itself indicate an AXI transition. Furthermore, when $\nu_z\neq 0$, defining $\theta$ uniquely requires the specification of a reference state and $\mathcal{I}$ center (*i.e.*, an origin) [@NicoDavidAXI2]. Correspondingly, if there are other bulk or surface contributions to the topological response (*e.g.*, other massive Dirac fermions at larger momenta, or a background QAH), then defects in $V_{\phi}$ will host additional states that coexist with and obscure the AXI bound states. Therefore, in order to fully determine the topology of the Weyl-CDW, we must examine the lattice-regularized UV completion \[Eqs. (\[eq:ham0\]) and (\[eq:cdwmass\])\] in detail.

#### Commensurate UV Completion:

When $Q=2\pi/Nc$, $N\in\mathbb{Z}^{+}$ in Eq. (\[eq:ham0\]), the CDW is lattice-commensurate, and Eq. (\[eq:ham0\]) remains periodic with respect to the reciprocal lattice vector ${\bf G} = Q\hat{\bf z}$ of the folded (reduced) BZ (rBZ). Folding into the rBZ translates bands with $|k_z|> \pi/(Nc)$ into the rBZ \[depicted in Fig. \[fig:zonefold\](b) for $N=3$\]. For all values of $N$, the two WPs fold into a linear fourfold (Dirac) [@Young12] degeneracy at the rBZ boundary \[Eq. (\[eq:kdp\])\]. From the bulk parity eigenvalues, we deduce that, for all $N\in\mathbb{Z}^{+}$, $C_z(|k_zc| < \pi/N)=-1\ ,C_z(|k_zc|>\pi/N)=0$ [^4]. This implies that $\nu_z=-1$ in the rBZ \[Fig. \[fig:zonefold\](b)\], *independent of whether $\phi=0,\pi$*. Combining $\nu_{z}=-1$ with the ${\bf k}\cdot{\bf p}$ analysis in the text preceding Eq. (\[eq:z4def\]) and fixing the origin to $z=0$ in the modulated cell, we find that $\phi=0$ \[$\phi=\pi$\] corresponds to a $\{2|001\}$ $\hat{\bf{z}}$-directed weak Chern (*i.e.* QAH) insulator \[$\{0|001\}$ oQAH insulator\] (see Fig. \[fig:scheme\]) with $\nu_{x,y}=0$, $\nu_{z}=-1$ and $\theta=0$ \[$\theta=\pi$\]. Despite the QAH and oQAH insulators differing by a translation $(Nc/2)\hat{\bf z}$, $\delta\theta_{\phi}=\pi$, independent of the choice of origin. In the SM, we provide further details, and analytically compute $\nu_{x,y,z}$ and $\delta\theta_{\phi}$ for $N=2$.

![Bulk ($\nu_{z}$) and slab ($\mathcal{G}_{z}$) Chern numbers for an $N=3$ ($Q=2\pi/3c$) commensurate CDW. (a) The $\hat{\bf y}$-surface spectral function at $E=0$ as a function of $k_x'$ and $k_z'$ in the surface rBZ exhibits a flat band. (b) The $\hat{\bf y}$-surface spectral function at $k_z'=0$ as a function of $E$ and $k_{x}'$ exhibits $C(3k'_zc)=-1$ spectral flow. The spectral functions in (a,b) are the same for $\phi=0,\pi$, and together imply that $\nu_{z}(\phi=0,\pi)=-1$. To measure the bulk axion angle $\theta$, we construct an $\mathcal{I}$-symmetric, $\hat{\bf z}$-directed slab of Eq. (\[eq:ham0\]) with 15 layers (5 unit cells) and $Q=2\pi/3c$. (c) The $\hat{\bf y}$-directed slab Berry phase [@wieder2018axion] exhibits $\mathcal{G}_{z}=-5$ ($\mathcal{G}_{z}=-4$) spectral flow for $\phi=0$ ($\phi=\pi$). Along with $\nu_{z}(\phi=0,\pi)=-1$ obtained from (a,b), the spectral flow in (c) implies through Eq. (\[eq:thetaconductance\]) that $\theta_{\phi=0,\pi}=0,\pi$, such that the Weyl-CDW at $\phi=0$ ($\phi=\pi$) is a QAH (oQAH) insulator \[Fig. \[fig:scheme\](b,c)\], and that $\delta\theta_{\phi}=\pi$. Further calculation details are provided in the SM.[]{data-label="fig:commensuratewilson"}](n3slab.pdf){width="0.8\columnwidth"}

To explicitly determine the bulk topology, we will employ model-agnostic numerical methods [@wieder2018axion; @vanderbiltaxion; @NicoDavidAXI2] to extract $\nu_{z}(\phi)$ and $\delta\theta_{\phi}$. To begin, we fix the origin to the $\mathcal{I}$ center at $(x,y,z)=(0,0,0)$ in the modulated cell, and form an $\mathbf{R}_i$-directed, $\mathcal{I}$-symmetric slab. The Hall *conductance* $G_{H,i}$ of the slab consists of an extensive contribution from the bulk QAH and an *intensive* contribution from $\theta$ that either reflects the bulk magnetoelectric polarizability or a QAH effect offset from the origin [@qi2008topological; @vanderbiltaxion; @NicoDavidAXI2]: $$G_{H,i}=\sigma_{H,i}L_i + (e^2\theta/h\pi),
\label{eq:thetaconductance}$$ where $\sigma_{H,i}=e^2\nu_i\mathbf{G}_i/2\pi h$ is the Hall conductivity (given by the weak Chern number $\nu_{i}$), where $L_i$ is the (lattice-regularized) thickness of the slab. Because a slab is a quasi-2D system, it carries a quantized Chern number $\mathcal{G}_{i}$ that is related to Eq. (\[eq:thetaconductance\]) by $\mathcal{G}_{i}e^2/h = G_{H,i}$. For $\mathcal{I}$-symmetric slab geometries, $\theta$ remains quantized to the bulk value, and provides an odd-integer contribution to Eq. (\[eq:thetaconductance\]) when $\theta\text{ mod }2\pi = \pi$ [@wieder2018axion] – this effect manifests in finite 3D systems with $\nu_{x,y,z}=0$ (*i.e.* AXIs) as chiral hinge modes \[Fig. \[fig:scheme\](a)\]. Therefore, given a fixed, $\mathcal{I}$-symmetric, $\mathbf{R}_i$-directed slab, and knowledge of $\nu_{i}$, changes in $\theta$ can be numerically extracted through Eq. (\[eq:thetaconductance\]). For the Weyl-CDWs in this Letter, our choice of origin corresponds to a convention in which $\theta=0$ when $\mathcal{G}_{z}=\nu_{z}L_{z}/(Nc)$ [^5].

#### $N=3$:

As an example, we analyze Eq. (\[eq:ham0\]) with commensurate $Q=2\pi/3c$. The tight-binding model in the rBZ exhibits the symmetry-based indicators $\{2|001\}$ ($\{0|001\}$) at $\phi=0$ ($\phi=\pi$) \[Eq. (\[eq:z4def\])\] (see SM). In Fig. \[fig:commensuratewilson\](a) \[Fig. \[fig:commensuratewilson\](b)\], we plot the $\mathbf{\hat{y}}$-normal surface spectral (Green’s) function as a function of $E$ \[the rBZ momentum $3k_z'c$\] and $k_x'a$, which are identical for $\phi=0,\pi$. Because surface Green’s functions do not capture hinge states [@wieder2018axion; @TMDHOTI; @HingeSM; @hotis; @schindler2018higher] (or origin-dependent changes in $\theta$), then the spectral flow in Fig. \[fig:commensuratewilson\](b) and the flat band in Fig. \[fig:commensuratewilson\](a) indicate a bulk contribution to the QAH from a weak Chern number $\nu_{z}(\phi=0,\pi)=-1$.

![(a) Band structure for the incommensurate Weyl-CDW in an $\mathcal{I}$-symmetric, $\hat{\bf x}$-directed rod geometry [@wieder2018axion; @TMDHOTI; @HingeSM; @hotis; @schindler2018higher] with 21 sites in the $\hat{\bf y}$ and $\hat{\bf z}$ directions and with $\phi=0$. The chiral states traversing the gap have degeneracy proportional to the rod thickness, indicating that they are QAH background surface states. When $\phi=\pi$, the rod spectrum is qualitatively the same as in (a), but exhibits one additional pair of chiral modes. (b) The $\hat{\bf y}$-directed Wilson loop of an $\mathcal{I}$-symmetric, $\hat{\bf z}$-directed slab [@wieder2018axion] of the model in (a) with 21 layers exhibits $\mathcal{G}_{z}=-9$ ($\mathcal{G}_{z}=-8$) spectral flow for $\phi=0$ ($\phi=\pi$). Along with the unchanged QAH contribution to $G_{H,z}$ in (a), this implies that $\delta\theta_{\phi}=\pi$.[]{data-label="fig:goldenrod"}](incomslab.pdf){width="0.8\columnwidth"}

Next, we cut the $N=3$ model in Fig. \[fig:commensuratewilson\](a,b) into an $\mathcal{I}$-symmetric, $\hat{\bf z}$-directed slab geometry with $L_{z}/3c=5$ unit cells, and calculate the $\hat{\bf y}$-directed non-abelian Berry phase (Wilson loop) $\mathcal{W}$ [@fidkowski2011bulk; @Yu11; @Alexandradinata16; @wieder2018axion] of the slab \[Fig. \[fig:commensuratewilson\](c)\], whose winding indicates that $\mathcal{G}_{z}(0)=-5=\nu_zL_{z}/3c,\ \mathcal{G}_{z}(\pi)=-4=\nu_zL_{z}/3c+1$. Along with $\nu_{z}(\phi=0,\pi)=-1$, the difference $|\delta \mathcal{G}_{z}| = 1$ indicates through Eq. (\[eq:thetaconductance\]) that $\theta_{\phi=0,\pi}=0,\pi$, such that the insulating Weyl-CDW is a QAH (oQAH) insulator for $\phi=0$ ($\phi=\pi$), and $\delta\theta_{\phi}=\pi$.

#### Incommensurate UV completion

: Having demonstrated that $\mathcal{I}$-symmetric Weyl-CDWs are either QAH or oQAH insulators in the commensurate case, we next explore the case of incommensurate modulation. Although an incommensurate CDW is not translationally-invariant, neither QAH nor oQAH phases require translation symmetry [@qi2008topological; @essin2009magnetoelectric; @wieder2018axion; @Hughes11; @Turner2012; @AFMTIMoore]. Consequently, Eq. (\[eq:thetaconductance\]) still applies, without modification, to $\hat{\bf z}$-directed slabs of Eq. (\[eq:ham0\]) with incommensurate values of $Q$.

To confirm this result, we cut $H_0+H_{CDW}$ \[Eqs. (\[eq:ham0\]) and (\[eq:cdwmass\])\] with $Q=\varphi\pi/2c$ (where $\varphi$ is the golden ratio) into an $\mathcal{I}$-symmetric rod geometry. First, we observe an extensive number of QAH surface states along the rod for $\phi=0,\pi$; \[$\phi=0$ is shown in Fig. \[fig:goldenrod\](a)\]; we additionally calculate that there is exactly one fewer chiral mode per surface at $\phi=\pi$. Next, to measure $\theta$, we cut the incommensurate Weyl-CDW into an $\mathcal{I}$-symmetric, $\hat{\bf z}$-directed slab geometry and calculate the $\hat{\bf y}$-directed Berry phase, as we did in the commensurate case \[Fig. \[fig:commensuratewilson\]\]. In the slab geometry, $|\delta \mathcal{G}_{z}|=1$ between $\phi =0,\pi$ \[Fig. \[fig:goldenrod\](b)\]. Furthermore, in incommensurate CDWs, a tuning cycle in $\phi$ changes the bulk wavefunctions, but not the bulk energy spectrum [@kraus2012topological; @de1983electrons], such that, if an incommensurate Weyl-CDW with $\phi=0$ is bulk-insulating, then an incommensurate Weyl-CDW with arbitrary $\phi$ is also bulk insulating. Additionally, the weak QAH invariant cannot change without a bulk gap closure, whereas the value of $\theta$ is free to wind between $0$ and $\pi$ at $\mathcal{I}$-breaking CDW angles away from $\phi=0,\pi$ [@qi2008topological; @essin2009magnetoelectric; @wieder2018axion]. As shown in the SM, $|\delta \mathcal{G}_{z}|=1$ in Fig. \[fig:goldenrod\](b), along with $G_{H,z}$ calculated for successive rational appoximants of an irrational $Q$, imply that, as in the commensurate case (Fig. \[fig:commensuratewilson\]), the incommensurate Weyl-CDW carries the relative axion angle $\delta\theta_{\phi}=\pi$.

#### Experimental Implications:

Our results have several implications for the search for axionic responses in Weyl-CDW materials. First, we have demonstrated that a large QAH effect is unavoidable and guaranteed by band topology for both commensurate and incommensurate Weyl-CDWs, independent of $\phi$. Second, the interplay between lattice and phase-angle defects, which both bind 1D chiral modes, is a fruitful area for future study, though one must cautiously separate contributions from $\theta$ and those from a background QAH effect. Next, we emphasize that the quantized axionic response $\delta\theta_{\phi}=\pi$ in Weyl-CDWs is measurable through the dynamical dependence of the quasi-2D QAH effect on $\phi$, rather than through the static magnetoelectric polarizability at fixed $\phi$ [@armitage2019matter; @lin2020prep]. Furthermore, soliton-like defects in $\phi$ can in principle be manipulated by exciting the CDW sliding mode. Finally, in $\mathcal{I}$-symmetric, magnetic Weyl-CDWs, our results highlight the experimental and theoretical difficulty of distinguishing QAH, oQAH, and AXI phases. However, our results do imply that, by carefully computing ${\bf k}$-space Chern numbers and then zone-folding, it is possible to predict the topology of Weyl-CDWs in real materials without performing intensive quasiperiodic calculations.

#### Outlook:

Going forward, our analysis can also be extended to CDWs and spin-density waves [@RahulSDW; @VladCDW; @BeenakkerKramersWeyl; @VTe2CDW] in $\mathcal{T}$-symmetric semimetals, including Dirac [@Young12; @ZJDirac; @ZJDirac2], Weyl, and nodal-line semimetals [@Kim2015; @fang2015topological; @Wieder2016; @LeslieLine1], which have recently been shown to exhibit signatures of higher-order topology [@TMDHOTI; @HingeSM; @TaylorToy; @HingeDiracExp]. Most interestingly, because rotation- and $\mathcal{T}$-symmetric HOTIs [@wieder2018axion; @fang2019new] can be formed from weak stacks of 2D TIs [@Po2017; @khalaf], then it is likely that rotation-symmetric CDWs in $\mathcal{T}$-symmetric WSMs, such as the CDW in (TaSe$_4$)$_2$I [@shi2019charge], may also exhibit nontrivial response effects. Specifically, a CDW can fold the four WPs in a minimal rotation- and $\mathcal{T}$-symmetric WSM into an eightfold double Dirac point for which $\mathcal{T}$-symmetric, line-like defects bind *helical* modes equivalent to HOTI hinge states [@Wieder2016b; @Bradlyn2016]. Additionally, recent experiments have demonstrated hinge-state-like step-edge helical modes and robust edge supercurrents in rotation- and $\mathcal{T}$-symmetric WSMs [@DavidMoTe2Exp; @MazWTe2Exp; @WTe2HingeStep; @PhuanOngMoTe2Hinge; @MoTe2WeylInterface]. This provides further motivation for the elucidation of quantized response effects in $\mathcal{T}$-symmetric HOTIs and Weyl-CDWs with trivial $\theta$ angles, which are currently unknown.

We acknowledge helpful conversations with Peter Abbamonte, B. Andrei Bernevig, Matthew Gilbert, as well as crucial input from Zhida Song on an early version of this Letter. BJW acknowledges support from B. Andrei Bernevig through Department of Energy Grant No. DE-SC0016239, Simons Investigator Grant No. 404513, BSF Israel US Foundation Grant No. 2018226, and ONR Grant No. N00014-20-1-2303.

[^1]: In this Letter we will use Greek indices as spacetime indices, and Roman indices as purely spatial indices

[^2]: Note that for gapless systems, such as Weyl semimetals, $C(\mathbf{k}\cdot \mathbf{R}_\mu)$ will depend on $\mathbf{k}$.

[^3]: We note that, because Eq. (\[eq:mainAction\]) is an action, then only $\theta\text{ mod }2\pi$ is a gauge-invariant quantity. This is reflected in the gauge-ambiguity of the Chern-Simons term \[Eq. (\[eq:thetadef\])\].

[^4]: Throughout this Letter, we choose parameters for which $\nu_z=-1$.

[^5]: This convention relies on $N$ being odd to avoid slabs with fractional modulated cells. In the SM, we consider a Weyl-CDW with $N=2$, which can only be cut into an $\mathcal{I}$-symmetric slab with a half-integer number of modulated cells.
