Surface States in Strain-Induced Nodal-Line Topological Semiconductors

This paper utilizes a minimalistic Luttinger Hamiltonian model to map the topological phase transitions of strained inverted-band-gap semiconductors among 3D topological insulators, Dirac, nodal-line, and Weyl semimetals, while deriving analytical surface state solutions that reveal their continuous evolution and a non-analytic dispersion feature at the projected nodal line.

The Gist

Imagine a crystal as a bustling city made of atoms. In most cities (standard semiconductors), the "traffic" of electrons flows smoothly, but there are strict rules about where they can and cannot go. However, in special materials like Mercury Telluride (HgTe), the city layout is "inverted." The usual rules are flipped, creating a unique environment where electrons behave like they are in a different dimension.

This paper explores what happens to the "surface traffic" (electrons living on the skin of the material) when we squeeze or stretch this crystal (apply strain) and introduce a specific type of magnetic twist (spin-orbit coupling).

Here is the story of their journey, explained through simple analogies:

1. The Stretchy City: Strain and Topology

Think of the material as a piece of rubber.

• Pulling it apart (Tensile Strain): When you stretch the rubber, you create a gap in the city. Electrons can't flow through the middle anymore. This turns the material into a Topological Insulator. It's like a city with a massive, empty moat in the center. However, the "surface" of the city has a special highway that runs right along the edge of the moat. Electrons can zip along this edge without getting stuck.
• Squeezing it together (Compressive Strain): When you squash the rubber, the moat disappears, and the city becomes a Dirac Semimetal. Now, the traffic flows freely through the center, but it does so in a very specific, cone-shaped way, like two ice cream cones touching at their tips.

2. The Magic Twist: Spin-Orbit Coupling

Now, imagine adding a "twist" to the city's rules. In the real world, this is called spin-orbit coupling (specifically from the crystal's lack of perfect symmetry).

• The Transformation: When this twist is added to the squeezed (compressed) city, the two touching ice cream cones (Dirac points) don't just stay as points. They stretch out into rings.
• The Nodal Line: These rings are called "nodal lines." Imagine a hula hoop floating in the middle of the city. Inside and outside the hoop, the rules are different. The hoop itself is a special boundary where the energy levels of the electrons cross over each other.

3. The Surface Highway: What Happens to the Edge?

The paper focuses on the "highways" that exist only on the surface of this material.

• The Smooth Ride: Without the "twist," these surface highways are smooth and predictable. They look like two lanes of traffic moving in opposite directions.
• The Kink in the Road: When the "twist" (spin-orbit coupling) is introduced, something strange happens to the surface highway as it crosses the projection of that floating hula hoop (the nodal line).
  • The road doesn't just bend; it jumps.
  • Imagine driving down a highway, and suddenly, at a specific point, the road doesn't just curve; it teleports to a slightly different elevation or changes its direction instantly. The paper calls this a non-analyticity. It's a mathematical "kink" where the rules of the road change abruptly.

4. The Patchwork Quilt: Spin Textures

The paper explains that this "kink" isn't just a glitch; it's a fundamental feature of the material's topology.

• The Mismatch: As the electron travels across this nodal line, its internal "spin" (think of it as a tiny compass needle attached to the electron) has to reorient itself.
• The Patchwork: Because of this reorientation, the surface state isn't one continuous, smooth ribbon. Instead, it's like a patchwork quilt. The electrons on one side of the nodal line belong to one "patch" with a specific spin pattern, and on the other side, they belong to a different patch.
• The Connection: The paper shows that these two patches are connected, but not in a simple way. They are linked through the nodal line like two different fabrics stitched together by a special, complex knot. You can't smoothly transition from one to the other without hitting that knot.

5. The Hierarchy of Scales: A Russian Nesting Doll

The authors also discovered that these different phases (Dirac, Nodal-Line, and Weyl) exist at different levels of energy, like a set of Russian nesting dolls:

1. The Big Doll (Dirac): You need a certain amount of energy to see the basic "ice cream cone" shape.
2. The Middle Doll (Nodal-Line): Inside that, you need to look closer (lower energy) to see the "hula hoop" rings form.
3. The Tiny Doll (Weyl): If you look even closer, the hoop breaks into tiny points (Weyl monopoles).
   The paper calculates that the "Tiny Doll" is so small that it might be very hard to see in a real experiment, but the "Middle Doll" (the Nodal Line) is clearly visible.

Summary

In short, this paper maps out the "traffic rules" for electrons on the surface of a special, strained crystal. It shows that when you twist the crystal's symmetry, the smooth surface highways develop a sudden, sharp "kink" exactly where they cross a special ring in the material's interior. This kink forces the electrons to switch their internal "compass" direction abruptly, creating a patchwork of different electron behaviors on the surface. The authors provide the exact mathematical formulas to predict exactly where these kinks happen and how the electron waves behave, unifying previous theories into one clear picture.

Technical Summary

Technical Summary: Surface States in Strain-Induced Nodal-Line Topological Semiconductors

Problem Statement
This work addresses the topological phase diagram of inverted-band-gap semiconductors, specifically focusing on HgTe and $\alpha$-Sn, under the influence of strain and spin-orbit coupling (SOC). While these materials are known to host topological surface states, the continuous evolution of these states across different phases—ranging from 3D topological insulators (TI) to Dirac semimetals, nodal-line semimetals, and Weyl semimetals—remains to be fully synthesized. A specific challenge lies in understanding how bulk inversion asymmetry (BIA), arising from the zinc-blende crystal structure of HgTe, modifies the surface state spectrum when the material is driven into a nodal-line phase by compressive strain. Previous literature often treats Dyakonov-Khaetskii (DK) and Volkov-Pankratov (VP) states within different models (Luttinger vs. Kane); this work seeks a unified description that captures the transition between these regimes and the specific non-analytic behaviors induced by the bulk nodal line.

Methodology
The authors employ a minimalistic Luttinger Hamiltonian model extended to include:

1. Strain effects: Modeled via the Bir-Pikus Hamiltonian, accounting for both tensile and compressive in-plane strain.
2. Spin-orbit coupling: Specifically, linear-in-momentum Dresselhaus terms ($H^{(1)}_{BIA}$) to account for bulk inversion asymmetry in the $T_d$ crystal class, and cubic-in-momentum terms ($H^{(3)}_{BIA}$) to explore the transition to Weyl semimetals.
3. Boundary Conditions: The study focuses on the (001) crystallographic surface. The authors derive exact analytical solutions for the surface state wave functions, utilizing a product ansatz of a spinor and an orbital wave function that satisfies the boundary condition $\psi(z=0)=0$.

The analysis is performed along high-symmetry directions of the in-plane momentum ($\mathbf{q} \parallel [110]$ and $\mathbf{q} \parallel [1\bar{1}0]$), allowing for the derivation of exact energy dispersions and spinor structures. The authors also calculate the projected density of states (PDOS) and identify van Hove singularities to characterize the bulk spectrum's influence on surface states.

Key Contributions and Results

• Unified Phase Diagram and Energy Hierarchy: The paper establishes a hierarchy of energy scales governing the sequential transitions between topological phases.

  • Tensile Strain: Opens a bulk gap, creating a 3D TI with surface states evolving from DK to VP types.
  • Compressive Strain: Creates a 3D Dirac semimetal with two Dirac points at finite $k_z$.
  • Inclusion of Linear SOC: Lifts the spin degeneracy of the Dirac points, transforming them into two nearly circular nodal lines. The energy scale for this transition is $\varepsilon_\alpha \approx \alpha\sqrt{3|\Delta|/\gamma_2}$.
  • Inclusion of Cubic SOC: Lifts the nodal line degeneracy everywhere except at four protected points, creating Weyl monopoles. The associated energy scale is $\varepsilon_\beta \approx \sqrt{3\alpha\beta\Delta}/(16\gamma_2\gamma_3)$, which is extremely small ($\sim \mu$eV) for HgTe, suggesting the nodal-line and Weyl phases may be difficult to distinguish experimentally in this material.

• Exact Analytical Solutions for Surface States: The authors derive closed-form expressions for the energy ($E$), spinor ($\chi$), and decay parameters ($\kappa, \Gamma$) of surface states along high-symmetry directions. They identify two distinct branches:

  • Branch $E_1$: An analytic branch that does not cross the nodal line.
  • Branch $E_2$: A branch that passes through the nodal line, requiring a distinction between an "inner" region ($q < q_{NL}$) and an "outer" region ($q > q_{NL}$).

• Non-Analyticity at the Nodal Line: A central finding is the non-analytic behavior of the surface states at the projected nodal line ($q_{NL} = \alpha/2\gamma_3$).

  • At this momentum, the surface state branch $E_2$ does not smoothly connect to its outer counterpart $\bar{E}_2$. Instead, the derivative $\partial E/\partial q$ is discontinuous.
  • More significantly, the spinor undergoes an abrupt reorientation. The overlap between the inner and outer spinors is $|\langle \chi_2 | \bar{\chi}_2 \rangle|^2 = \gamma_1^2 / (4\gamma_2^2)$. Since $\gamma_1/2\gamma_2 \approx 0.8$ for HgTe, the branches are neither fully connected nor fully disconnected.
  • This behavior is described as a non-abelian connection mediated by the nodal line, which acts as a degenerate subspace. The surface states form a "patch structure" in momentum space, organized into separate regions connected only through this singular point.

• Emergence of Secondary Dirac Points: Due to the spin-orbit induced shifts and the non-analytic connection, the surface state spectrum along the $[1\bar{1}0]$ direction exhibits a secondary 2D Dirac point at finite momentum outside the nodal line, in addition to the primary Dirac point at $q=0$.

• Van Hove Singularities: The paper details two types of van Hove singularities in the projected density of states ($vH_0$ and $vH_1$) and shows how they evolve with SOC, forming an envelope for the bulk spectrum extrema.

Significance and Claims
The authors claim that this work provides a unified picture synthesizing previous literature on DK and VP states, demonstrating their continuous evolution across phase boundaries. The primary significance lies in revealing that the bulk nodal line imposes a topological constraint on the surface states, manifesting as a geometric singularity (non-analyticity) and a patch structure in momentum space.

The paper asserts that while the surface spectrum resembles spin-split branches of a conventional 2D system, the underlying bulk topology enforces a distinct topological character via the non-abelian connection at the nodal line. Furthermore, the established hierarchy of energy scales ($\Delta$, $\varepsilon_\alpha$, $\varepsilon_\beta$) defines the experimental resolution required to distinguish between Dirac, nodal-line, and Weyl phases. The authors conclude that these findings advance the fundamental understanding of how bulk topology manifests in the surface electronic structure of gapless topological materials, offering specific signatures (such as the secondary Dirac points and non-analytic dispersion) for experimental detection in strained heterostructures.