Paper The following article is Open access

Nexus and Dirac lines in topological materials

and

Published 11 September 2015 © 2015 IOP Publishing Ltd and Deutsche Physikalische Gesellschaft
, , Citation T T Heikkilä and G E Volovik 2015 New J. Phys. 17 093019 DOI 10.1088/1367-2630/17/9/093019

1367-2630/17/9/093019

Abstract

We consider the Z2 topology of the Dirac lines, i.e., lines of band contacts, on an example of graphite. Four lines—three with topological charge ${N}_{1}=1$ each and one with ${N}_{1}=-1$—merge together near the H-point and annihilate due to summation law $1+1+1-1=0$. The merging point is similar to the real-space nexus, an analog of the Dirac monopole at which the Z2 strings terminate.

Export citation and abstract BibTeX RIS

Content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.

1. Introduction

Dirac points in 2D systems and Dirac lines in 3D are examples of the exceptional points and lines of level crossing analyzed by von Neumann and Wigner [1]. They are typically protected by symmetry and are described by the Z2 topological invariant (see [2]). Close to the 2D Dirac point with nontrivial topological charge N1, the energy spectrum after deformation can be represented by the $2\times 2$ matrix ${H}_{\mathrm{Dirac}}=c({\sigma }_{1}{p}_{x}\pm {\sigma }_{2}{p}_{y})$, where the Pauli matrices ${\sigma }_{\mathrm{1,2}}$ describe the pseudo-spin induced in the vicinity of the level crossing. This gives rise to the conical spectrum near the Dirac point $E={E}_{0}\pm c| {{\bf{p}}}_{\perp }| $. For Dirac lines in 3D, the components px and py are in the transverse plane.

The Dirac Hamiltonian anticommutes with ${\sigma }_{3}$, which allows us to have the analytic form for the topological charge N1, see e.g. review [3]:

Equation (1)

Here C is an infinitesimal contour in momentum space around the Dirac point or the Dirac line. The topological charge N1 in equation (1) is integer: ${N}_{1}=1$ for sign + and ${N}_{1}=-1$ for sign −. However, the integer-valuedness emerges only in the vicinity of the Dirac point or line. In general, the summation rule is $1+1=2\equiv 0$. This means that the Dirac line with ${N}_{1}=2$ can be continuously deformed to the trivial configuration.

In time reversal symmetric superconductors, due to chiral symmetry the Dirac lines may have zero energy and thus correspond to the nodal lines in the spectrum. Such lines exist in cuprate superconductors. According to the bulk-boundary correspondence the nodal lines may produce a dispersionless spectrum on the boundary—the flat band with zero energy [46]. If time reversal symmetry is violated (for example, by supercurrent), the Dirac line acquires a nonzero energy, a Fermi surface is formed, and the Dirac line lives inside the Fermi surface. In cuprate superconductors, the Fermi surfaces created by supercurrent around Abrikosov vortices give rise, at zero temperature, to the finite density of states proportional to $\sqrt{B}$, where B is the magnetic field [7].

Examples of semimetals with 2D Dirac points and 3D Dirac lines are provided by graphene and graphite, correspondingly. In these materials the spin–orbit interaction can be ignored, and one can consider them as spinless materials. In graphene there are two Dirac points in the Brillouin zone. If graphene is treated as spinless, one Dirac point has ${N}_{1}=+1$ and another one has ${N}_{1}=-1$ (for spinful electrons these are ${N}_{1}=+2$ and ${N}_{1}=-2$, correspondingly).

In rhombohedral graphite (ABCABC... stacking) two Dirac points of graphene layers generate two well-separated Dirac lines. Each has a form of a spiral and is characterized by the topological charge ${N}_{1}=+1$ or ${N}_{1}=-1$ [8, 9]. These Dirac lines have nonzero energy and live inside the chains of the hole and electron Fermi surfaces, discussed by McClure (see figure 2 in [10]). The projection of the spiral to the top or bottom surface of graphite determines the boundary of the (approximate) surface flat band [8, 9]. This flattening of the spectrum has been recently observed in epitaxial rhombohedral multilayer graphene[11].

2. Bilayer graphene

In Bernal graphite (ABAB... stacking) the geometry of the Dirac lines is essentially different from rhombohedral graphite, because Bernal graphite is the 3D extension of bilayer graphene with AB stacking. Let us hence start by describing the bilayer graphene with AB stacking.

Since in each graphene layer the topological charges are ${N}_{1}=+1$ and ${N}_{1}=-1$ in the two valleys, in bilayer graphene the charges are summed up giving rise to trivial topological charges ${N}_{1}=\;+2\equiv 0$ and ${N}_{1}=-2\equiv 0$. Interaction between the layers may lead to several possible scenarios of the geometry of the fermionic spectrum in bilayer graphene:

  • (i)  
    If there is some special symmetry, such as the fundamental Lorentz invariance, one obtains a doubly degenerate conical spectrum. However, in bilayer graphene there is no symmetry which could support such a scenario.
  • (ii)  
    The topological charge ${N}_{1}=2$ gives rise to the Dirac fermions with parabolic energy spectrum near the Dirac point, $E={E}_{0}\pm {{\bf{p}}}_{\perp }^{2}/2m$ (figure 1(left)). Such a spectrum belongs to the trivial element of the group Z2, that is why the neglected hopping elements destroy the parabolic spectrum. The Dirac point disappears, and finite gap (the Dirac mass) emerges.
  • (iii)  
    The Dirac point with ${N}_{1}=+2$ splits into two Dirac points, each with ${N}_{1}=+1$ and with conical spectrum (see [12] for the relativistic 3+1 system). Such splitting, however, violates the hexagonal symmetry of graphene.
  • (iv)  
    The splitting of the ${N}_{1}=+2$ Dirac point, which is consistent with the hexagonal symmetry, is in figure 1(right). The ${N}_{1}=+2$ Dirac point splits into four Dirac conical points: the central Dirac point has ${N}_{1}=-1$, while three others with ${N}_{1}=+1$ each are connected by the C3 symmetry, see figure 1(right). This is called the trigonal warping. The topological charges ${N}_{1}=\pm 1$ belong to the nontrivial element of the group Z2. That is why each Dirac point is topologically stable and is not destroyed by addition of the neglected hopping elements, if they obey the time reversal and sublattice symmetry.

Figure 1.

Figure 1. Z2-topology of Dirac points in bilayer graphene. Left: dirac point with quadratic spectrum characterized by topological invariant ${N}_{1}=2$. Since the topological invariant is trivial in Z2 topology, such a spectrum is unstable towards the gap formation, or to splitting. Right: splitting to topologically stable Dirac points due to trigonal warping.

Standard image High-resolution image

The bilayer graphene chooses the scenario (iv) with the trigonal warping [13, 14], which prohibits the annihilation of the Dirac point due to splitting into the topologically stable Dirac points with ${N}_{1}=\pm 1$. Altogether the Brillouin zone of the bilayer graphene contains 8 Dirac conical points.

3. From topology of bilayer graphene to Bernal graphite

The trigonal warping also stabilizes the Dirac lines in the Bernal graphite (ABAB... stacking). Altogether the Brillouin zone of the Bernal graphite contains 8 Dirac lines [15, 16], which originate from 8 Dirac conical points of bilayer graphene. These Dirac lines have finite energy and live inside the Fermi surface of graphite.

However, the topologically stable Dirac lines appear only in a certain region of momenta pz. Outside that region the scenario (ii) without the Dirac lines takes place. This means that there exists an exceptional point in the spectrum, at which four topologically stable Dirac lines merge together and annihilate each other [15, 16]. Such a point is a direct analog of the nexus, a point in the real space where the Z2 or Z4 topological defects merge and annihilate each other, see figure 2(top) [18].

Figure 2.

Figure 2. Real-space and momentum space nexus. Top left: real-space nexus in the dipole locked superfluid 3He-A, where vortices are described by Z2-topology. Two vortices, with ${N}_{1}=1$ each, merge together at the hedgehog in the field of the orbital momentum vector $\hat{{\bf{l}}}$ and annihilate each other since ${N}_{1}=2\equiv 0$. Top right: real-space nexus in the dipole unlocked superfluid 3He-A, where vortices are described by Z4-topology. Four half-quantum vortices, with ${N}_{1}=1/2$ each, merge together at the hedgehog and annihilate each other. Bottom: illustration of the topology of the momentum-space nexus in graphite. Dirac lines obey the Z2-topology. Four Dirac lines, with ${N}_{1}=+1$, ${N}_{1}=+1$, ${N}_{1}=+1$ and ${N}_{1}=-1$, merge together at an exceptional point of the spectrum and annihilate each other, because ${N}_{1}=1+1+1-1=2\equiv 0$.

Standard image High-resolution image

Let us find this exceptional point. According to ${D}_{6h}$ symmetry of graphite, the Dirac lines are situated in the vertical mirror planes, and thus the exceptional point is on the H-K-H line in the Brillouin zone at ${p}_{z}={P}^{*}$.

The Hamiltonian of Bernal graphite is a 3D extension of the 2D Hamiltonian describing bilayer graphene:

Equation (2)

This corresponds to a rotated version of that written in [10] (we write the Hamiltonian in a basis spanned by the two layers ×sublattice points, whereas [10] writes it in the eigenbasis defined by the strongest interlayer couplings), but we assume that the intralayer coupling ${\gamma }_{0}$ dominates over the other terms, and therefore also make the $\sigma \cdot p$ approximation for each layer. Here $p{{\rm{e}}}^{{\rm{i}}\phi }={p}_{x}+{\rm{i}}{p}_{y}$; ${\rm{\Gamma }}({p}_{z})=\mathrm{cos}((\pi /2a){p}_{z})$, where a is the distance between the K-point and H-point in the Brillouin zone of graphite and ${p}_{z}=0$ at the K-point so that ${\rm{\Gamma }}(K)=1$ and ${\rm{\Gamma }}(H)=0$. The function ${\gamma }_{2}({p}_{z})={\gamma }_{2}{{\rm{\Gamma }}}^{2}({p}_{z})/2$ describes the coupling across two layers and Δ denotes the locally broken A-B sublattice symmetry, which still preserves the global A-B symmetry. The coefficients ${\gamma }_{\mathrm{3,4}}$ are related to the original tight-binding coefficients [10] ${\tilde{\gamma }}_{\mathrm{3,4}}$ via ${\gamma }_{\mathrm{3,4}}={\tilde{\gamma }}_{\mathrm{3,4}}/{\gamma }_{0}$. In graphite, ${\gamma }_{3}\approx 0.1$ and ${\gamma }_{4}\approx 0.01$ [17]. We neglect spin–orbit coupling and below we set ${v}_{F}=1$ for simplicity.

Hamiltonian (2) is invariant under time reversal symmetry (complex conjugation) combined with reflection from plane ${p}_{y}=0$ ($\phi \to -\phi $). This symmetry supports the Z2 topology of Dirac points in graphene (figure 1(right)) and Dirac lines in graphite.

4. Point of merging of Dirac lines

Since the merging point ${P}^{*}$ is on the H-K-H line, we consider the Hamiltonian at p = 0 which has the following form (we consider only the traceless part):

Equation (3)

This can be rewritten in terms of the Pauli matrices ${\sigma }_{i}$ and ${\tau }_{i}$:

Equation (4)

Let us introduce new matrices ${\tilde{\tau }}_{3}=-{\tau }_{3}$, ${\tilde{\tau }}_{1}=-{\tau }_{1}$, ${\tilde{\tau }}_{2}={\tau }_{2}$ and the 'total spin' ${\bf{J}}=\frac{1}{2}({\boldsymbol{\sigma }}+{\boldsymbol{\tau }})$:

Equation (5)

Equation (6)

The Hamiltonian at p = 0 can be expressed in terms of quantum numbers J and J3:

Equation (7)

The exceptional point ${p}_{z}={P}^{*}$ is at $f({p}_{z})+g({p}_{z})=0$, where three branches of the spectrum with J = 1 are degenerate. Its position is determined by equation

Equation (8)

In other notations this is the known equation for the exceptional point in graphite, see e.g. [16]. It is situated in the vicinity of the H-point. Away from this point, two branches with ${J}_{3}=\pm 1$ remain degenerate for all pz. Close to the exceptional point

Equation (9)

where $\lambda =f^{\prime} ({P}^{*})+g^{\prime} ({P}^{*})$.

The exceptional point ${P}^{*}$ is the merging point of four Dirac lines living in the symmetry planes at ${p}_{z}\lt {P}^{*}$. At ${p}_{H}\gt {p}_{z}\gt {P}^{*}$ the spectrum remains degenerate: the branches with ${J}_{3}=+1$ and ${J}_{3}=-1$ have the same energy. But these are the band contact lines rather than the Dirac lines: the branches have the same sign of energy counted from the line position, as distinct from ${p}_{z}\lt {P}^{*}$, where the touching branches have energies with opposite sign, which is the characteristics of a Dirac line (see figure 3). This means that at ${p}_{z}={P}^{*}$ a band inversion occurs. The band inversion simultaneously happens for the branch with quantum numbers J = 1 and ${J}_{3}=0$.

Figure 3.

Figure 3. (a) Positions of the Dirac lines in the energy spectrum are seen as the points in momentum space where the two eigenenergies meet (plotted for ${\rm{\Gamma }}({p}_{z})=0.5$). (b) Cut through the $\phi =0$ axis revealing the positions of two of the Dirac lines. With a large $| {\rm{\Delta }}| $ they shift towards lower values of p. (c) Same as previous but with ${\rm{\Gamma }}({p}_{z})=-2{\rm{\Delta }}/{\gamma }_{1}$, showing also the higher-energy branches. Now the finite-momentum Dirac lines are too close to the center point to be visible. (d) Dispersion at the nexus, ${\rm{\Gamma }}({p}_{z})=-{\rm{\Delta }}/(2{\gamma }_{1})$; (e) above the nexus, ${\rm{\Gamma }}({p}_{z})=-{\rm{\Delta }}/(4{\gamma }_{1})$ and (f) at the H-point ${\rm{\Gamma }}({p}_{z})=0$. The figures are computed with ${\rm{\Delta }}=-0.05{\gamma }_{1}$, ${\gamma }_{3}=0.1$, ${\gamma }_{2}={\gamma }_{4}=0$.

Standard image High-resolution image

5. 2D Dirac fermions at H-point

To see that the topology at ${p}_{H}\gt {p}_{z}\gt {P}^{*}$ is trivial, let us consider the ${p}_{x/y}$ plane containing the H-point. Since ${\rm{\Gamma }}({p}_{H})=0$, the Hamiltonian at the plane ${p}_{z}={p}_{H}$ is

Equation (10)

It represents two copies of massive 2D fermions with Dirac spectrum:

Equation (11)

and thus two doubly degenerate branches. In addition to the contact line of ${J}_{z}=\pm 1$ branches, the other two branches, J = 1, ${J}_{3}=0$ and $J={J}_{3}=0$, contact each other. The Hamiltonian (10) has trivial N1 topology. Due to continuity, the topology in the plane slightly below the H-point is also trivial. This demonstrates that at ${p}_{H}\gt {p}_{z}\gt {P}^{*}$ there are no Dirac lines. Instead there are band-contact lines, with trivial topological charge ${N}_{1}=0$. Thus the exceptional point ${P}^{*}$ is the point of merging of four Dirac lines.

6. Finite-energy Dirac lines

Let us consider the occurrence of the Dirac lines at a finite energy. For simplicity, we neglect ${\gamma }_{2}({p}_{z})$ and ${\gamma }_{4}$, as they do not change the qualitative behavior of the Dirac lines close to the nexus. In this case, the characteristic equation for the bulk spectrum reads

Equation (12)

This has four solutions for epsilon. Let us analyze the solutions for three cases: (a) $| {p}_{z}| \ll {P}^{*}$, (b) $| {p}_{z}| ={P}^{*}$, and (c), ${p}_{H}\geqslant | {p}_{z}| \gt {P}^{*}$. The relevant eigensolutions for a specific set of parameters in the three cases are plotted in figure 3.

For region (a), there are two eigenenergies $\epsilon \approx \pm {\gamma }_{1}{\rm{\Gamma }}({p}_{z})$, and two closer to zero. The latter two are plotted in figures 3(a), (b) for a fixed ${\rm{\Gamma }}({p}_{z})$. Two low-energy solutions touch at specific points in the transverse momentum direction, either at p = 0 and close to $\mathrm{cos}(3\phi )=-1$, $p={p}_{0}=4{\gamma }_{1}{\gamma }_{3}{{\rm{\Gamma }}}^{2}({p}_{z})$, which is the exact position of the Dirac line for ${\rm{\Delta }}=0$. For small Δ, the Dirac lines at finite momentum obtain a finite energy $\epsilon ={\rm{\Delta }}[1-1/(1+4{\gamma }_{3}^{2}{\rm{\Delta }}{{\rm{\Gamma }}}^{2}({p}_{z}))]$ and slightly shift the line away from $p={p}_{0}$. The Dirac line character however persists. Figure 3(c) shows all four eigenenergies in region (a). Here the splitting of the two low-energy lines is not visible as they are quite close to each other.

For ${p}_{z}={P}^{*}$, the four solutions are plotted in figure 3(d). For low Δ and p, the energies are given by

Equation (13)

There are hence two gapped and quadratic and two gapless linear branches, and three of them meet at $\epsilon =p=0$.

Above the nexus (region (c)), we may first set ${\rm{\Gamma }}({p}_{z})=0$, i.e., consider the H point. There, the finite-energy solutions are given by (shifted by a constant energy from equation (11))

Equation (14)

The energies, plotted in figure 3(f), are doubly degenerate, but the degeneracy is lifted for a finite ${\rm{\Gamma }}({p}_{z})$, yielding two pairs of quadratic branches, the gap depending on the precise value of ${\rm{\Gamma }}({p}_{z})$. The lifting of this degeneracy is shown in figure 3(e).

7. Nexus

In the real space, the nexus is a kind of a Dirac monopole that terminates cosmic strings (the so-called Z-string in the Standard Model) [19]. In chiral superfluid 3He-A the nexus is a hedgehog which terminates vortices with the Z2 or Z4 topology (figure 17.3 in the book [18] and figure 2(top)). In the dipole-locked 3He-A the homotopy group for vortices is Z2. That is why two vortices can terminate at the hedgehog in the field of the orbital momentum vector $\hat{{\bf{l}}}$, which is the analog of the Dirac magnetic monopole.

In the same manner the Dirac lines can be described by the Z2 group, and if so, they may terminate on a Dirac monopole in some field. If the invariant N1 belongs to the Z2 group, which means that ${N}_{1}=2$ and ${N}_{1}=0$ are equivalent, then by some transformation of hopping elements the line described by even ${N}_{1}=2k$ can disappear. The lines with ${N}_{1}=\pm 1$ are stable. But when four such lines meet each other, their total charge is even, and thus they may annihilate each other. This also happens, since at ${p}_{z}\gt {P}^{*}$ the band contact lines are topologically trivial, ${N}_{1}=0$.

So, the point ${p}_{z}={P}^{*}$ represents the nexus in momentum space, see figure 2(bottom). It is distinct from the 3D Weyl point in momentum space, which also represents the momentum-space analog of a Dirac monopole. Such a monopole contains the Dirac string: a singularity in the Berry phase which terminates at the monopole, see figure 11.4 in [18]. But this string is not observable.

8. Conclusion

The Dirac lines are suggested to exist in different semimetals [2026]. The materials with Dirac lines are important, because they may have an (approximate) flat band on the boundary or at the interface between materials with different topological properties. The high density of states in the flat band provides a possible route to room-temperature superconductivity [27].

However, the Bernal graphite is a very instructive example. Its spectrum demonstrates the possible interplay of several relevant topological invariants. One of them characterizes the local line element as in [2]. Another one is the global invariant, which characterizes for example the closed nodal ring as a whole: it is the same invariant which characterizes the 3D Dirac or Weyl point obtained by shrinking the nodal ring. There are also the topological invariants that characterize the Fermi surface(s). These include the local invariant of the Fermi surface [18] and the topology of the shape of the Fermi surface [28]. All this may combine to produce exotic topological patterns. The interplay of Weyl point and Fermi surface topologies with exchange of the Berry flux, when two Fermi surfaces contact each other, have been discussed in [26, 29]. The Bernal graphite provides an example of the interplay of Dirac lines in the mirror planes, an exceptional point (nexus), and Fermi surfaces with touching points between electron and hole pockets, see [15, 16].

Acknowledgments

This work was supported by the Academy of Finland through its Center of Excellence program and by the European Research Council (Grant No. 240362-Heattronics).

Please wait… references are loading.
10.1088/1367-2630/17/9/093019