Dirac Magnon Nodal Loops in Quasi-2D Quantum Magnets

TL;DR

This paper introduces the concept of 2D Dirac magnon nodal-line loops in quasi-2D quantum magnets, demonstrating their topological protection, robustness, and potential realization in layered honeycomb ferromagnets.

Contribution

It proposes the existence of 2D Dirac magnon nodal-line loops in quasi-2D quantum magnetic systems, a novel topological feature in magnon band structures.

Findings

Magnon bands form 1D closed lines of Dirac nodes in bilayer honeycomb ferromagnets.

These nodal-line loops are topologically protected by inversion and time-reversal symmetry.

They are robust against weak Dzyaloshinskii-Moriya interactions and support chiral magnon edge modes.

Abstract

In this report, we propose a new concept of one-dimensional (1D) closed lines of Dirac magnon nodes in two-dimensional (2D) momentum space of quasi-2D quantum magnetic systems. They are termed "2D Dirac magnon nodal-line loops". We utilize the bilayer honeycomb ferromagnets with intralayer coupling $J$ and interlayer coupling $J_L$, which is realizable in the honeycomb chromium compounds CrX$_3$ (X $\equiv$ Br, Cl, and I). However, our results can also exist in other layered quasi-2D quantum magnetic systems. Here, we show that the magnon bands of the bilayer honeycomb ferromagnets overlap for $J_L\neq 0$ and form 1D closed lines of Dirac magnon nodes in 2D momentum space. The 2D Dirac magnon nodal-line loops are topologically protected by inversion and time-reversal symmetry. Furthermore, we show that they are robust against weak Dzyaloshinskii-Moriya interaction $ \Delta_{DM}< J_L$ and possess chiral magnon edge modes.}

Figures

Equations

$$\mathcal{H}$$

$$E_{\pm }^{\alpha }(\mathbb{k})=E_0\pm S[|J{f}_{\mathbb{k}}|+(-1{)}^{\alpha }|J_L|],$$

$$q_x^2+q_y^2=J_L^2/c_s^2,$$

$$q_x^2+q_y^2=[J_L^2-(3J{)}^2]/c_0^2,$$

$$\mathcal{H}(q)=E_0{\mathbb{I}}_{2\times 2}-v_0(\pm q_x{\sigma }_x+q_y{\sigma }_y),$$

$$(-1{)}^{\nu }=\prod _{i=1}^4{\delta }_i,$$

$${\mathcal{H}}_{DM}={\Delta }_{DM}\sum _{⟨⟨i,j⟩⟩;\tau }{\nu }_{ij}\hat{z}\cdot {\mathbf{S}}_{i,\tau }\times {\mathbf{S}}_{j,\tau },$$

$$E_{\pm }^{\alpha }(\mathbb{k})=E_0\pm S[\sqrt{|J{f}_{\mathbb{k}}{|}^2+m_{\mathbb{k}}^2}+(-1{)}^{\alpha }{J}_L],$$

$$q_x^2+q_y^2=(J_L^2-{\Delta }_{DM}^2)/c_s^2,$$

$${\mathcal{C}}_{\pm }^{\alpha }=\frac{1}{2\pi }{\int }_{BZ}{d}^{2}k{\Omega }_{\mathbb{k}\pm }^{\alpha }.$$

$$S_{i,\tau }^z=S-a_{i,\tau }^{†}{a}_{i,\tau },$$

$$S_{i,\tau }^y=i\sqrt{\frac{S}{2}}(a_{i,\tau }^{†}-a_{i,\tau }),$$

$$S_{i,\tau }^x=\sqrt{\frac{S}{2}}(a_{i,\tau }^{†}+a_{i,\tau }),$$

$${\mathcal{H}}_{\mathbb{k}}=E_0{I}_{\tau }\otimes I_{\sigma }-J_{L}S{\tau }_x\otimes {\sigma }_x-JS{\tau }_x\otimes [{\sigma }_{+}{f}_{\mathbb{k}}+{\sigma }_{-}{f}_{\mathbb{k}}^{\ast }],$$

$${\sigma }_{\pm }\to {\tau }_x{\sigma }_{\pm },{\tau }_{\pm }\to {\sigma }_x{\tau }_{\pm }.$$

$${\mathcal{H}}_{\mathbb{k}}=(E_0{I}_{\tau }-J_{L}S{\tau }_x)\otimes I_{\sigma }-JS{I}_{\tau }\otimes [{\sigma }_{+}{f}_{\mathbb{k}}+{\sigma }_{-}{f}_{\mathbb{k}}^{\ast }].$$

$${\mathcal{H}}_{\mathbb{k}}^{\alpha }=[E_0+(-1{)}^{\alpha }|J_{L}S|]I_{\sigma }-JS[{\sigma }_{+}{f}_{\mathbb{k}}+{\sigma }_{-}{f}_{\mathbb{k}}^{\ast }],$$

$${\psi }_{\mathbb{k}\pm }^{\alpha }=\frac{1}{\sqrt{2}}\left(\genfrac{}{}{0ex}{}{1}{(-1{)}^{\alpha }}\right)\otimes \frac{1}{\sqrt{2}}\left(\genfrac{}{}{0ex}{}{1}{\pm e^{i{\varphi }_{\mathbb{k}}}}\right),$$

$$\mathcal{P}:{\mathcal{H}}_{\mathbb{k}}\to I_{\tau }\otimes {\sigma }_x{\mathcal{H}}_{-\mathbb{k}}{I}_{\tau }\otimes {\sigma }_x.$$

$$\mathcal{H}(\mathbb{k}={\Gamma }_i)=d_0+d_1(\mathbb{k}={\Gamma }_i)\mathcal{P},$$

$$(-1{)}^{\nu }=\prod _{i=1}^4{\delta }_i.$$

Full text

The experimental observations of topological insulators top3 ; top4 and topological semimetals Xu ; lv1 ; lv2 in electronic systems have inspired a great interest in condensed matter physics. Consequently, this has led to the re-examination of band structures in bosonic systems. One of the areas of recent interest is the topological magnon bands in insulating ordered quantum magnets with inversion symmetry breaking lifa ; alex4 ; alex5a ; shin1 ; sol ; sol1 ; kkim , which allows a Dzyaloshinskii-Moriya (DM) spin-orbit interaction (SOI) dm ; dm2 . The study of topological magnetic spin excitations in quasi-two-dimensional (2D) quantum magnetic systems is currently an active research field in condensed matter physics both theoretically and experimentally. Topological magnonic systems are expected to open the next frontier of condensed matter, because they are potential candidates towards magnon spintronics and magnon thermal devices magn . Unlike electronic charge particles, magnons are charge-neutral bosonic quasiparticles and they do not experience a Lorentz force and do not have conduction and valence bands. However, a temperature gradient can induce a heat current and the Berry curvature induced by the DM SOI acts as an effective magnetic field in momentum space. This leads to a thermal version of the anomalous Hall effect characterized by a temperature dependent thermal Hall conductivity as predicted theoretically alex0 ; alex5 ; alex2 ; alex22 and subsequently observed experimentally alex1a ; alex1 ; alex6 .

Recently, Mook et al. mok1 have shown that the magnon bands in the three-dimensional (3D) anisotropic pyrochlore ferromagnets without DM SOI form a magnon analogue of 3D electronic nodal-line semimetals aab2 ; aab3 ; aab4 ; aab5 . However, the concept of the quasi-2D magnon nodal-line has never been proposed in quantum magnetic systems. In principle, there are more quantum magnetic materials in quasi-2D forms than 3D forms. In addition, quasi-2D quantum magnetic systems are simpler to study and manipulate both theoretically and experimentally. Therefore, it is interesting to seek for a 2D analogue of Dirac magnon nodal-lines. Nevertheless, the concept of 2D electronic Dirac nodal-line semimetals has remained elusive for a while. It was recently put forward and it is predicted to exist in composite lattices such as mixed honeycomb and kagomé lattices lu , mixed square lattice yang , and 2D trilayers che . In this case, the mixed energy bands have the tendency to overlap in momentum space and form 1D closed lines of Dirac nodes in 2D momentum space. A magnonic analogue of these systems should be of interest in the study of quantum magnetism both theoretically and experimentally.

In this report, we propose a 2D Dirac magnon nodal loop or 1D closed lines of Dirac magnon nodes in 2D momentum space as a direct analogue of 2D electronic Dirac nodal-line semimetals yang ; lu ; che . We utilize the honeycomb bilayer ferromagnet as an example, which can be realized in the hexagonal chromium compounds CrX3 (X $\equiv$ Br, Cl, and I) with honeycomb lattices coupled by small interlayer interaction dav0 ; dav ; dav1 ; foot ; dejon . However, the general idea of 1D closed lines of Dirac magnon nodes in 2D momentum space can be easily extended to other quasi-2D lattice structures with more than two sublattices in the unit cell. We show here that in the realistic regime with the interlayer coupling ($J_{L}$) smaller than the intralayer coupling ($J$), i.e. $J_{L}<J$, there are Dirac magnon nodal loops centred at the corners of the Brillouin zone. We show that the 2D Dirac nodal-line loops are topologically protected by the $\mathbb{Z}_{2}$ invariance of the parity eigenvalue resulting from the presence of inversion symmetry. In principle, the DM SOI with strength $\Delta_{DM}$ can be allowed on the honeycomb lattice because of inversion symmetry breaking between the bonds of second-nearest neighbour sites sol . With the inclusion of DM SOI, we show that the quasi-2D Dirac magnon nodal-line loops are robust for $\Delta_{DM}<J_{L}$. This is in stark contrast to 3D magnon nodal-lines for which the DM SOI transforms the magnon nodal-lines into magnon Weyl points mok1 .

Results.

Model. We consider the ferromagnetic spin Hamiltonian of a bilayer honeycomb lattice, given by

[TABLE]

The first term describes the intralayer nearest-neighbour (NN) interaction on the top ($\tau=T$) and bottom ($\tau=B$) layers respectively, and the second term is the interlayer interaction. The last term is the Zeeman magnetic field in units of $g\mu_{B}$ applied out of the bilayer plane (in the $z$-direction). The bilayer honeycomb lattice can be stacked in two forms: AB-stacked and AA-stacked. In the former the top layer is slightly shifted from the bottom layer, whereas they are placed right above each other in the latter. We will assume the latter as the former does not possess the Dirac magnon nodal-line loops proposed in this report. Also, in the latter the Hamiltonian (1) is a good approximation of the honeycomb ferromagnetic chromium compounds CrX3 (X $\equiv$ Br, Cl, and I) dav0 ; dav ; dav1 ; foot ; dejon . The intralayer coupling is ferromagnetic $J>0$ and the small interlayer coupling can be ferromagnetic $J_{L}>0$ as in X $\equiv$ Br and I dav0 ; dav ; dav1 or antiferromagnetic $J_{L}<0$ as in X $\equiv$ Cl foot . We note that CrCl3 and CrBr3 have an estimated value of $|J_{L}|/J\sim 0.004$ and $J_{L}/J\sim 0.06$ respectively dav0 ; dav . These values are very close to the magnetic susceptibility measurements dejon . As we will show later, the Dirac magnon-nodal lines exist for any value of $J_{L}\neq 0$. In the presence of a strong magnetic field, the antiferromagnetic coupled layers can be mapped to the ferromagnetic coupled layers at the saturation field $H_{s}=2J_{L}S$, where $S$ is the value of the spin. Therefore, we will consider only ferromagnetic coupled layers. Due to small interlayer coupling $J_{L}<J$, it is reasonable to consider the quasi-2D system.

Magnon band structures. At low-temperatures the magnetic excitations of ferromagnetic materials can be measured by inelastic neutron scattering experiment dav1 ; alex6 and they correspond to the quanta of spin waves also known as magnons. The linearized Holstein-Primakoff (HP) hp spin-boson mapping suffices in the low-temperature regime. The corresponding magnon bands of Eq. (1) are given by (see Methods)

[TABLE]

where $E_{0}=3JS+J_{L}S+H$ and $f_{\mathbb{k}}=1+e^{-ik_{a}}+e^{-i(k_{a}+k_{b})}$, $\hat{a}=\sqrt{3}\hat{x}$ and $\hat{b}=-\sqrt{3}\hat{x}/2+3\hat{y}/2$ with $k_{a}=\mathbb{k}\cdot\hat{a}$ and $k_{b}=\mathbb{k}\cdot\hat{b}$. Here, the lattice constant is denoted by $a$ and $\alpha=1$ or $2$. The magnetic field $H$ only introduces a gap at the Goldstone mode at $\mathbb{k}=0$ and therefore it has no effects on the Dirac magnon points or the Dirac magnon nodal-lines. Henceforth, we set $H=0$. The topmost band is $E_{+}^{2}$ and the bottommost band is $E_{-}^{2}$ and the two middle bands are $E_{\pm}^{1}$. For theoretical purposes, we will vary the value of $J_{L}/J$ in order to magnify the Dirac magnon nodal-lines as well as the Dirac points. For instance, the magnon bands are depicted in Fig. 2(a) along the Brillouin zone paths in Fig. 1(b) for $J_{L}/J=0.8$. The topmost band $E_{+}^{2}$ and the lower middle band $E_{+}^{1}$ form a Dirac magnon point at $\pm{\bf K}$. Similarly, the lowest band $E_{-}^{2}$ and the upper middle band $E_{-}^{1}$ form a Dirac magnon point at $\pm{\bf K}$. This is very similar to single-layer honeycomb ferromagnet mag . In other words, the Dirac magnon points do not necessarily require a bilayer form.

Dirac magnon nodal-line loops. The bilayer ferromagnets have another important features that are not possible in single-layer ferromagnetic form mag . As evident from Fig. 2(a) the middle magnon bands $E_{\pm}^{1}$ in the bilayer honeycomb ferromagnets overlap, say, along ${\bf\Gamma}$–${\bf K}$ line in the Brillouin zone and they form 1D closed lines of Dirac nodes in 2D momentum space or 2D Dirac magnon nodal-line closed loops. They are different from the Dirac magnon points formed by discrete band touching at the high symmetry $\pm{\bf K}$-points. They are the exact analogues of 2D Dirac node-line semimetals in electronic systems that also exist in mixed lattices such as the honeycomb and kagomé lattices lu , mixed square lattice yang , and 2D trilayers che . Indeed, for any value of $J_{L}\neq 0$ there exists 2D Dirac magnon nodal lines in the magnon band structures. The magnon density of states (DOS) in Fig. 2(b) linearly cross between the bands forming the Dirac magnon nodal-loops and they are distinguished from that of Dirac magnon points. As in electronic systems the 2D materials offer additional tunability which may not be available in 3D systems. If the interlayer coupling $J_{L}$ is considered as a tunable parameter we see from Fig. 3(a) that upon increasing $J_{L}$ the Dirac magnon points remain intact while the Dirac magnon nodal closed loops along ${\bf K}$–${\bf M}$ form a line node near the ${\bf M}$-point. The corresponding magnon density of states in Fig. 3(b) also show a line of degeneracy consistent with the magnon bands.

One of the advantages of magnon bands in quasi-2D quantum magnets is that they offer an explicit analytical simplification, which may not be possible in 3D systems such as the pyrochlore ferromagnets mok1 . In this respect, the existence of Dirac magnon nodal-loops and Dirac magnon points are consequences of the magnon band structures and can be directly understood from Eq. 2. In magnonic systems both the Dirac magnon points and the Dirac magnon nodal-loops occur at finite energy $E_{0}$. The solution of $E_{\pm}^{\alpha}=E_{0}$ determines the nodal-line degeneracy of the bands. There is no solution for the topmost and bottommost bands $E_{\pm}^{2}$, which is consistent with the absence of Dirac magnon closed loops. However, the zeros of the middle bands $E_{\pm}^{1}=E_{0}$ give the solution where the two bands overlap in momentum space, i.e., 2D Dirac magnon nodal-line closed loops.

The main result of this report is depicted in Fig. 4. For $J_{L}<J$ the Dirac magnon nodal loops are centred at the corners of the Brillouin zone. At the symmetric point $J_{L}=J$ there are additional line nodes connecting the Dirac magnon nodal loops. Although $J_{L}>J$ is usually unrealistic, there are still Dirac magnon nodal loops in this regime as shown in Fig. 4, but now they are centred at the ${\bf\Gamma}$-point and shrink as $J_{L}$ increases and subsequently turn into a quadratic point node at the unrealistic limit critical value $J_{L}^{c}=3J$.

To get an intuition about the development of Dirac magnon nodal loops we consider the low-energy limit. Near ${\bf k}={\bf K}$ the zeros of $E_{\pm}^{1}=E_{0}$ gives

[TABLE]

where $c_{s}=3J/2$ and $\mathbb{q}=\mathbb{k}-{\bf K}$. The radius of the Dirac magnon nodal loops $|J_{L}/c_{s}|$ increases with $J_{L}$ consistent with Fig. 4 for $J_{L}<J$. On the other hand, near ${\bf k}={\bf\Gamma}$ the zeros of $E_{\pm}^{1}=E_{0}$ gives

[TABLE]

where $c_{0}=3J/\sqrt{2}$ and $\mathbb{q}=\mathbb{k}-{\bf\Gamma}$. As noted above the Dirac magnon nodal loops shrink to a quadratic point node at the critical value $J_{L}=J_{L}^{c}$ as can be seen in Fig. 4.

By linearizing $\mathcal{H}_{\mathbb{k}}$ (see Methods) near the Dirac magnon points and the Dirac magnon nodal-line loops, and projecting onto the subspace associated with the middle bands, we obtain a 2D Dirac-like Hamiltonian

[TABLE]

where $v_{0}=c_{s}(c_{0})$ is the spin wave velocity. The linearized Hamiltonian has a Berry phase kai defined as $\gamma=\oint_{\mathcal{C}}\mathcal{A}(\mathbb{q})\cdot d{\mathbb{q}}$, where $\mathcal{A}(\mathbb{q})$ is the Berry connection given by $\mathcal{A}(\mathbb{q})=i\braket{\psi_{\mathbb{q}}}{{\vec{\nabla}}_{\mathbb{q}}\psi_{\mathbb{q}}}$. The Berry phase is defined for a closed path encircling the Dirac point nodes in momentum space with $\gamma=\pi$, whereas $\gamma=0$ otherwise. Therefore the Dirac point nodes can be interpreted as topological defects in 2D momentum space. On the other hand, the 2D Dirac magnon nodal-line loops are topologically protected through the parity eigenvalues at the time-reversal-invariant momenta (TRIM) points due to the presence of inversion symmetry. In two dimensions there are four TRIM points at $\Gamma_{i}=(\boldsymbol{\Gamma},3\mathbb{M}^{\prime}s)$ points, where the $\mathbb{M}$ points are the bisects of the Brillouin zone corners. The $\mathbb{Z}_{2}$ invariance $\nu$ is given by fu

[TABLE]

where $\delta_{i}=\prod_{n=1}^{2N}\xi_{n}(\Gamma_{i})$, and $\xi_{n}(\Gamma_{i})$ is the parity eigenvalue associated with the magnon bands forming the Dirac nodal-line loops. As shown in Methods $\delta_{i}=-\text{sgn}\left(d_{1}(\mathbb{k}=\Gamma_{i})\right)$, where $d_{1}(\mathbb{k})=-JS\text{Re}f_{\mathbb{k}}$. We find that $\nu=1$ which yields a nonzero topological invariant $\mathbb{Z}_{2}=1$. This confirms the topological protection of the odd number of Dirac nodal-line loops in the 2D Brillouin zone.

Effects of Dzyaloshinskii-Moriya interaction. The honeycomb lattice can also allow a DM SOI due to inversion symmetry breaking at the midpoint of the bonds connecting the second-nearest neighbours which form a triangular plaquette. The form of the DM SOI on the honeycomb lattice is given by sol

[TABLE]

where $\nu_{ij}=\pm$ for hopping from right to left on each of the two-layer sublattices. The DM SOI is imaginary in the spin-boson mapping and the resulting magnon band structures are given by

[TABLE]

where $m_{\mathbb{k}}=-m_{-\mathbb{k}}=2\Delta_{DM}[\sin k_{a}+\sin k_{b}-\sin(k_{a}+k_{b})]$. The combination of ferromagnetic spontaneous magnetization and DM SOI breaks time-reversal symmetry macroscopically and the Dirac points are gapped leading to a topological magnon insulator in single-layer honeycomb ferromagnets also known as the Haldane magnon insulator sol ; kkim with finite thermal Hall effect sol1 .

It is interesting to ask what happens to the Dirac magnon nodal loops and the Dirac magnon points as the DM SOI is taken into account in the bilayer system. We find that the Dirac magnon points are gapped due to broken time-reversal symmetry. On the other hand, the Dirac magnon nodal loops are not gapped by the DM SOI in the realistic limit $\Delta_{DM}<J_{L}$ as shown in Fig. 5; therefore they are robust in this regime. Nevertheless, the opposite limit $\Delta_{DM}>J_{L}$ has gap Dirac magnon loops and Dirac magnon points and the system becomes a full topological magnon insulator as depicted in Fig. 6. These interesting features can be well-understood in the low-energy limit by expanding the magnon energy bands near ${\bf k}={\bf K}$. The zeros of $E_{\pm}^{1}=E_{0}$ give

[TABLE]

where $\Delta_{DM}=3\sqrt{3}\tilde{\Delta}_{DM}$. Evidently, the presence of Dirac magnon nodal loop requires $\Delta_{DM}<J_{L}$.

Another interesting feature of this system is the existence of chiral magnon edge modes on the boundary. We have solved for the magnon edge modes using a strip geometry with open boundary conditions along the $y$-direction and infinite along $x$-direction. As shown in Fig. 7 there are four chiral magnon edge modes in both regimes $\Delta_{DM}<J_{L}$ (a) and $\Delta_{DM}>J_{L}$ (b) with each band contributing a single edge mode. The topological protection of the magnon edge modes is encoded in the Chern numbers of the magnon bulk bands. For a gap topological system the Berry curvature is given by $\Omega_{\mathbb{k}\pm}^{\alpha}=\left(\boldsymbol{\nabla}\times\boldsymbol{\mathcal{A}}_{\pm\mathbb{k}}^{\alpha}\right)_{z}$, where $\boldsymbol{\mathcal{A}}_{\pm\mathbb{k}}^{\alpha}=i\braket{\psi_{\pm\mathbb{k}}^{\alpha}}{\boldsymbol{\nabla}}{\psi_{\pm\mathbb{k}}^{\alpha}}$ is the Berry connection and $\psi_{\pm\mathbb{k}}^{\alpha}$ are the eigenvectors of $\mathcal{H}_{\mathbb{k}}$ given in Methods. The Berry phase $\gamma$ now defines a Chern number given by the integration of the Berry curvature over the Brillouin zone,

[TABLE]

The Berry curvature has its maximum at $\pm{\bf K}$. Therefore, the major contributions to the Chern numbers come from the states near the Dirac points. By direct numerical integration, we find that for the topmost and the upper middle bands $E_{+}^{2}$ and $E_{-}^{1}$ respectively the Chern numbers are $\mathcal{C}_{\mp}^{1(2)}=\pm 1$. Similarly, for the bottommost and the lower middle bands $E_{-}^{2}$ and $E_{+}^{1}$ we find $\mathcal{C}_{\mp}^{2(1)}=\pm 1$. This is consistent with the number of chiral magnon edge modes in the system.

Conclusion. We have shown that the magnon bands in the quasi-2D quantum magnetic systems offer another new topological distinction — the quasi-2D Dirac magnon nodal loop or 1D closed lines of Dirac magnon nodes in 2D momentum space. They are the direct analogues of 2D electronic Dirac nodal-line semimetals in composite lattices lu ; yang . The Dirac magnon-nodal loops occur when two magnon bands overlap in momentum space as opposed to Dirac magnon points which occur when two magnon bands touch at discrete points at the high symmetry points of the Brillouin zone. We also showed that the quasi-2D Dirac magnon-nodal loops can exist even in the presence of DM SOI as opposed to 3D counterparts on the pyrochlore ferromagnets mok1 . It would be of interest to search for the existence of Dirac magnon-nodal loops in quantum magnetic systems via the inelastic neutron scattering experiments. Recently, many different topological magnonic systems have been proposed; they include Weyl magnon points in 3D pyrochlore antiferromagnets fei and ferromagnets mok ; su , and also Dirac points in quasi-2D quantum magnetic systems with DM SOI sol4 ; oku , which can be categorized as the magnonic analogue of 2D Dirac semimetals ste . The present work has never been studied in any quasi-2D quantum magnetic system; hence it completes the magnonic analogues of electronic topological semimetals in 2D systems. The present work is applicable to the honeycomb chromium compounds CrX3 (X $\equiv$ Br, Cl, and I) dav0 ; dav ; dav1 ; foot ; dejon . Therefore, the magnon bands of these honeycomb magnetic materials should be re-examined experimentally in the context of topological magnonics. We note that the study of magnonic analogues of electronic topological semimetals is still developing. There are very few experiments alex1a ; alex1 ; alex6 and lots of theoretical proposals have not been achieved experimentally. Because of the bulk sensitivity of inelastic neutron scattering methods, the chiral magnon edge modes have not been experimentally measured in any quantum magnetic system, but the bulk topological magnon bands have been realized in the kagomé lattice ferromagnet Cu(1-3, bdc) alex5a . Therefore, edge sensitive methods should be applied in the study of topological magnonics. In this respect, it be would interesting to implement other experimental techniques such as light luuk and electronic kha scattering methods. We also note that the electronic properties of the honeycomb chromium compounds CrX3 (X $\equiv$ Br, Cl, and I) have been studied by density functional theory and other methods sem1 ; sem2 and they possess semiconducting features.

Methods The formalism we implement in this report is the linear spin wave theory. This approximation is valid at low-temperature in the ordered regime when few magnons are thermally excited. We assume that the spins are polarized along the $z$-axis, which can be achieved by applying a magnetic field along this direction. The linearized Holstein-Primakoff (HP) spin-boson mapping is given by

[TABLE]

where $a_{i,\tau}^{\dagger}(a_{i,\tau})$ are the bosonic creation (annihilation) operators. There are four sites in the unit cell, i.e. $\tau=1,2$ on the top layer and $\tau=3,4$ on the bottom layer. Next, we substitute the spin-boson transformation into Eq. 1 and Fourier transform the quadratic bosonic Hamiltonian which is similar to that of AA-stacked graphene dav2 ; dav3 ; davv , apart from the constant energy shift $E_{0}=3JS+J_{L}S+H$ that pushes the negative bands to positive bands as expected for magnonic (bosonic) systems. In other words, there is no conduction and valence bands in magnonic systems as the spin excitations are charge-neutral and obey Bose-Einstein statistics. Instead, the magnon bands are thermally populated at low temperatures. The momentum space Hamiltonian $\mathcal{H}_{\mathbb{k}}$ can be written as $\mathcal{H}=\sum_{\mathbb{k}}\psi^{\dagger}_{\mathbb{k}}\cdot\mathcal{H}_{\mathbb{k}}\cdot\psi_{\mathbb{k}},$ where $\psi_{\mathbb{k}}^{\dagger}=(a_{\mathbb{k},1}^{\dagger},a_{\mathbb{k},4}^{\dagger},a_{\mathbb{k},3}^{\dagger},a_{\mathbb{k},2}^{\dagger})$ is the eigenvectors. The momentum space Hamiltonian is given by

[TABLE]

where $f_{\mathbb{k}}=1+e^{-ik_{a}}+e^{-i(k_{a}+k_{b})}$, $\hat{a}=\sqrt{3}\hat{x}$ and $\hat{b}=-\sqrt{3}\hat{x}/2+3\hat{y}/2$ with $k_{a}=\mathbb{k}\cdot\hat{a}$ and $k_{b}=\mathbb{k}\cdot\hat{b}$. Here, $\boldsymbol{\tau}$ and $\boldsymbol{\sigma}$ are Pauli matrices acting on the layer space and sublattice space respectively. Whereas $I_{\tau}$ and $I_{\sigma}$ are $2\times 2$ identity matrix in each space and $\sigma_{\pm}=(\sigma_{x}\pm i\sigma_{y})/2$. In order to diagonalize the Hamiltonian, we perform the following canonical (i.e. commutation relation preserving) transformation of the $\boldsymbol{\sigma}$ and $\boldsymbol{\tau}$ operators

[TABLE]

The Hamiltonian Eq. 14 now takes the form

[TABLE]

Next, we diagonalize $(E_{0}I_{\tau}-J_{L}S\tau_{x})$ and the resulting Hamiltonian is given by

[TABLE]

where $\alpha=1,2$ for the top and bottom layers respectively. The eigenvalues of Eq. 17 yield the magnon bands given in Eq. (2). The corresponding normalized eigenvectors are tensor product of the two spaces, i.e.

[TABLE]

where $\phi_{\mathbb{k}}=\tan^{-1}\left({\text{Im}f_{\mathbb{k}}}/{\text{Re}f_{\mathbb{k}}}\right)$; Re and Im denote the real and imaginary parts.

The inversion operation changes the sign of the momentum and interchanges the two sublattices in each layer:

[TABLE]

It follows that Eq. 14 is invariant under inversion symmetry $\mathcal{P}$. Similar to graphene fu , at the time-reversal-invariant momenta (TRIM) points the Hamiltonian has the form

[TABLE]

where $d_{1}(\mathbb{k})=-JS\text{Re}f_{\mathbb{k}}$. Therefore the parity eigenvalues $\xi_{n}$ for the states at $\mathbb{k}=\Gamma_{i}$ are given by the eigenvalues of $\mathcal{P}$. Hence, $\delta_{i}=-\text{sgn}\left(d_{1}(\mathbb{k}=\Gamma_{i})\right)$ and the $\mathbb{Z}_{2}$ invariance $\nu$ is given by fu

[TABLE]

Indeed, the DM spin-orbit interaction has no effects on the values $\delta_{i}$. It only ensures that a topological gap exists in the system. However, in the absence of any topological gap the $\mathbb{Z}_{2}$ invariance $\nu$ defines the topological protection of the 2D nodal-line loops yang ; che .

References