Nonhermitian topological zero modes at smooth domain walls: Exact solutions

This paper derives exact analytical solutions for nonhermitian topological zero modes at smooth domain walls, revealing a universal relation between scalar fields and the decay rates and oscillation wavelengths of these boundary modes that quantifies the bulk-boundary correspondence in line-gapped systems.

The Gist

Imagine you are walking through a vast, magical landscape made of two different types of terrain. On the left, you have a "Forest of Order" (a topologically trivial phase), and on the right, you have a "Forest of Chaos" (a topologically non-trivial phase). Usually, in physics, we think of the boundary between these two forests as a sharp cliff. But in this paper, the authors are looking at what happens when the transition is a smooth, rolling hill instead of a cliff.

Furthermore, this world isn't perfectly quiet and stable; it's a bit "noisy." Some parts of the forest are gaining energy (like a microphone picking up too much sound), while others are losing it (like a battery dying). In physics, we call this non-Hermitian physics.

Here is the story of what the authors discovered, broken down into simple concepts:

1. The "Ghost" at the Border

In the world of quantum physics, when two different materials meet, something special often happens right at the edge. A "ghost" particle appears. It's called a Zero Mode. It's a wave that gets stuck right at the border, unable to wander off into the deep forest on either side.

• The Old Way: Scientists used to know how to find these ghosts if the border was a sharp cliff and the world was perfectly stable (Hermitian).
• The New Discovery: This paper asks: "What if the border is a smooth hill? And what if the world is noisy (gaining and losing energy)?"

2. The "Jackiw-Rebbi" Equation: The Map

To find these ghosts, the authors used a special mathematical map called the Jackiw-Rebbi equation. Think of this equation as a GPS that tells you exactly where a wave will go.

• In the old days, this GPS was simple and only worked for quiet, stable worlds.
• The authors upgraded the GPS. They made it complex (allowing for imaginary numbers, which represent the gain and loss of energy) and flexible (allowing the terrain to change smoothly).

3. The "Hair" on the Ghosts

This is the most creative part of the paper. The authors classify these ghost particles based on how "fuzzy" they are. They use a funny analogy from black holes: "Hair."

• Black Holes: In physics, there's a saying that "Black holes have no hair." This means you can describe a black hole completely with just three numbers: its mass, its spin, and its charge. You don't need to know the details of what fell into it.
• The "No Hair" Ghosts: If the border between the two terrains is a sharp, instant cliff, the ghost particle is simple. It's "bald." You can describe it perfectly with just two numbers: how fast it fades away (decay rate) and how fast it wiggles (oscillation). It doesn't care about the details of the border; it only cares about the two sides.
• The "Hairy" Ghosts: But if the border is a smooth hill (a smooth domain wall), the ghost gets "hair." It becomes complex. Its shape depends on the exact details of the hill in the middle.
  • Short Hair: If the ghost is small compared to the hill, it looks simple from far away, but up close, it has a complex shape.
  • Long Hair: If the ghost is huge compared to the hill, it is complex everywhere. You can't describe it with just two numbers; you need to know the whole shape of the hill.

4. The Universal Rule (The "Recipe")

The authors found a magical recipe that connects the ghost to the terrain.
They discovered a universal relationship between:

1. The Terrain: How much energy is being gained or lost in the bulk material (the "scalar fields").
2. The Ghost: How fast the ghost fades away and how fast it wiggles.

Why is this cool?
Imagine you are a detective. You can't see the ghost directly because it's too small. But you can measure how fast the ghost fades away and how fast it wiggles. Using this new "recipe," you can work backward to figure out the exact properties of the terrain (the gain and loss) without ever seeing the ghost itself. This turns a theoretical math problem into something scientists can actually measure in a lab.

5. Why Does This Matter?

This isn't just about abstract math. This helps us understand:

• New Materials: How to build materials that conduct electricity only on the edges, even if they are losing energy or are very "noisy."
• Lasers and Light: Many modern lasers use "gain and loss" to work. This paper helps predict how light will behave in these complex systems.
• The Future of Computing: These "ghost" particles are candidates for building quantum computers that are resistant to errors.

Summary

In short, this paper is like upgrading a weather forecast.

• Before: We could only predict the weather if the storm was a sudden, sharp front and the air was calm.
• Now: The authors have figured out how to predict the weather even when the storm is a slow-moving, rolling front and the air is turbulent. They found that the "storms" (the zero modes) can be simple and bald, or complex and hairy, depending on how the terrain changes. Most importantly, they found a rule that lets us measure the storm's behavior to understand the atmosphere itself.

Technical Summary

1. Problem Statement

The bulk-boundary correspondence (BBC) is a fundamental principle in condensed matter physics stating that topological invariants of the bulk determine the existence of localized boundary modes. While well-established for Hermitian systems, extending this to non-Hermitian systems (which describe open systems with gain and loss) presents challenges.

• Limitations of Previous Work: Existing non-Hermitian BBC theories often rely on sharp domain walls (step-function transitions between phases). However, real physical systems typically exhibit smooth domain walls where scalar fields (mass and velocity) vary continuously over a finite width $w$.
• The Gap: There was a lack of analytical solutions for zero-energy modes at smooth domain walls in non-Hermitian regimes. Specifically, the relationship between the spatial variation of fields, the decay rates, oscillation wavelengths, and the topological invariants in these smooth, non-Hermitian systems was not fully understood.
• Key Question: Can one derive exact analytical wavefunctions for non-Hermitian zero modes at smooth domain walls, and can a universal relation be established between bulk physical quantities and the observable properties of these boundary modes?

2. Methodology

The authors employ an analytical approach based on a generalized Jackiw-Rebbi equation (a modified Dirac equation) to model the system.

• Model Hamiltonian: They consider a 1D system described by a second-order differential equation with space-dependent mass $m(x)$ and Dirac velocity $v(x)$, allowing for complex coefficients (non-Hermitian):
  $$\left[ (\eta p^2 + m(x))\sigma_z + 2v(x)p \sigma_y \right] \Psi(x) = E\Psi(x)$$
  where $p = -i\partial_x$. The focus is on zero-energy modes ($E=0$).
• Smooth Domain Wall Assumption: The fields $m(x)$ and $v(x)$ are assumed to vary smoothly over a characteristic width $w$, approaching constant asymptotic values ($m_{L,R}, v_{L,R}$) as $x \to \pm \infty$.
• Mathematical Transformation:
  • The authors map the infinite real line $(-\infty, \infty)$ to a finite segment $(0, 1)$ using the substitution $y(x) = \frac{1}{2}(1 + \tanh(x/2w))$.
  • They assume the fields can be expanded as polynomials in $y(x)$ (up to second order for mass, first order for velocity).
  • This transformation reduces the differential equation to a hypergeometric differential equation.
• Exact Solutions: The wavefunctions are derived analytically in terms of Gaussian hypergeometric functions (${}_2F_1$). The general solution is a linear combination of two independent solutions determined by the asymptotic exponents.

3. Key Contributions

A. Derivation of Exact Wavefunctions

The paper provides the first exact analytical form of zero-mode wavefunctions for non-Hermitian systems with smooth domain walls. The solutions are expressed as:
$$\Psi(x) \propto \begin{pmatrix} 1 \\ s \end{pmatrix} \left[ A_1 \phi_1(x) + A_2 \phi_2(x) \right]$$
where the spatial dependence involves hypergeometric functions ${}_2F_1$ with arguments dependent on the domain wall width $w$ and the asymptotic field values.

B. Generalization of Bulk-Boundary Correspondence (BBC)

The authors define a non-Hermitian topological invariant $W$ and a non-Hermitian topological mass $M$ based on the asymptotic values of the fields:
$$M = |\text{Re}(\sqrt{v^2 + m})| - |\text{Re}(v)|$$
$$W = \begin{cases} \text{sgn}(\text{Re}(v)) & \text{if } \text{sgn}(M) \le 0 \\ 0 & \text{if } \text{sgn}(M) > 0 \end{cases}$$
They prove that the number of zero modes localized at the boundary is determined by $|W_L - W_R|$, extending the BBC to non-Hermitian line-gapped systems with smooth transitions.

C. Universal Relation between Bulk and Boundary

A major contribution is the derivation of a universal relation connecting the bulk properties to the boundary mode characteristics:
$$\mu_{L,R} + i\kappa_{L,R} = \pm v_{L,R} \pm \sqrt{v_{L,R}^2 + m_{L,R}}$$
Here, $\mu$ represents the decay rate and $\kappa$ the oscillation momentum. This equation holds regardless of the specific spatial profile of the domain wall (as long as it is smooth), linking the asymptotic decay/oscillation of the wavefunction directly to the bulk parameters.

D. Classification of Zero Modes ("Hair")

The authors introduce a novel classification for zero modes based on the relationship between the domain wall width $w$ and the localization length $\xi$:

• "No Hair" (Featureless): Occurs in the sharp limit ($w \to 0$). The mode is fully described by decay rates and oscillation wavelengths on either side.
• "Short Hair": $w < \xi$. The mode appears featureless at large distances but depends on field details near the wall.
• "Long Hair": $w > \xi$. The mode is non-featureless at all length scales, with wavefunctions heavily dependent on the specific spatial variation of the fields.

4. Results

• Existence Conditions: The existence of zero modes is strictly determined by the asymptotic values of the fields. For an infinite system, zero modes exist if $W_L \neq W_R$. For a semi-infinite system, they exist if $W_R \neq 0$.
• Oscillatory vs. Exponential Decay: The behavior of the wavefunction (pure exponential decay vs. exponentially damped oscillations) is determined by the quantity $K = |\text{Im}(\sqrt{v^2+m})| + |\text{Im}(v)|$.
  • If $K=0$, the mode decays exponentially without oscillation.
  • If $K \neq 0$, the mode exhibits exponentially damped oscillations.
• Complex Wavefunctions: Unlike Hermitian Majorana modes which can be real, non-Hermitian zero modes generally possess complex wavefunctions with non-uniform phases, even in the absence of oscillations.
• Local Topological Invariant: The authors define a local topological invariant $W(x)$ that varies across the domain wall. Zero modes localize near points where $W(x)$ changes sign, providing a mechanism for localization even in complex field profiles.

5. Significance

• Experimental Relevance: The derived universal relation ($\mu, \kappa$ vs. $m, v$) offers a direct experimental protocol. By measuring the spatial decay and oscillation of boundary modes (via spatially resolved spectroscopy), one can infer the bulk topological properties and verify the non-Hermitian BBC without needing to probe the bulk directly.
• Beyond Sharp Walls: This work moves beyond the idealized "sharp wall" approximation, providing a more realistic framework for analyzing topological phases in photonic lattices, superconductors with gain/loss, and metamaterials where interfaces are naturally smooth.
• Theoretical Generalization: It demonstrates that the bulk-boundary correspondence is robust even when the system is non-Hermitian and the transition is smooth, provided the system possesses a line gap (which can be adiabatically deformed to a Hermitian gap).
• New Physical Insights: The "hair" classification highlights that the internal structure of the domain wall matters for the wavefunction's local properties, a nuance often missed in topological theories that focus solely on asymptotic invariants.

In summary, this paper provides a rigorous analytical foundation for understanding non-Hermitian topological zero modes in realistic, smooth interfaces, bridging the gap between abstract topological invariants and measurable physical quantities like decay lengths and oscillation wavelengths.