Topological phases of the Bogoliubov de Gennes Hamiltonian

This paper investigates a two-dimensional superconducting system with a periodically modulated order parameter, demonstrating that the modulation period determines the topological winding number and establishing an analytic bulk-edge connection that explains the emergence of localized edge states from the bulk spectrum.

The Gist

The Big Picture: Dancing Electrons on a Twisted Stage

Imagine a superconductor not as a boring, flat sheet of metal, but as a giant, magical dance floor. In a normal superconductor, the electrons (the dancers) move in perfect, synchronized pairs. They glide without friction, creating a supercurrent.

Usually, this dance floor is flat and uniform. But in this paper, the author, Klaus Ziegler, asks: What happens if we twist the dance floor?

He imagines a superconductor shaped like a ring (like a wedding band or a tire). If you push a supercurrent around this ring, it creates a "vortex"—a giant whirlpool of energy. This twist forces the "dance steps" (the quantum phase of the electrons) to rotate as you go around the ring.

The paper investigates how this twisting affects the topology of the system. In physics, "topology" is like the shape of a donut versus a coffee mug. You can stretch a mug into a donut without tearing it, but you can't turn a donut into a sphere without making a hole. The "winding number" is a count of how many times the dancers twist around the ring before they return to their starting point.

The Main Characters: The Spinors and the Bloch Vector

To understand the dancers, the author uses two main tools:

1. The Spinor (The Dancer's Pose): This describes the exact state of an electron pair. Think of it as a dancer's pose. If the dancer spins around the ring, their pose changes.
2. The Bloch Vector (The Shadow on the Wall): This is a 3D arrow that represents the dancer's pose. Imagine shining a light on the dancer; the shadow cast on a wall is the Bloch vector.
   • If the dancer spins once around the ring, the shadow on the wall might trace a circle.
   • If the dancer spins twice, the shadow traces the circle twice.

The Winding Number is simply counting how many full circles the shadow makes on that wall as you walk once around the ring.

The Discovery: Twisting Creates New Paths

The paper finds two very different types of "dancers" (solutions to the equations) depending on where they are standing:

1. The Bulk Dancers (The Crowd in the Middle)
These are the electrons in the middle of the superconductor.

• The Analogy: Imagine a crowd of people walking in a circle on a flat track. If the track is twisted (the order parameter has a phase), the crowd's "shadow" (the Bloch vector) winds around a specific number of times.
• The Result: The number of twists (the winding number) is directly determined by how much the superconductor is twisted. If you increase the supercurrent, you increase the twist, and the winding number goes up. It's like winding a rubber band tighter; the more you twist it, the higher the number of coils.

2. The Edge Dancers (The Ghosts on the Wall)
These are the special electrons that only exist at the very edge or boundary of the material.

• The Analogy: Imagine the edge of the dance floor is a cliff. Some dancers get stuck there, unable to leave. These are "edge modes."
• The Surprise: The paper shows that these edge dancers behave differently than the crowd in the middle. Their "winding number" depends on the specific rules of the edge (boundary conditions).
• The Connection: The author uses a clever mathematical trick (analytic continuation) to show how the "bulk" dancers (in the middle) can morph into "edge" dancers (at the boundary) when the energy conditions change. It's like realizing that a person walking in the middle of a room can suddenly become a ghost haunting the corner if the lighting changes just right.

The "Magic" Transformation

The author uses a mathematical "magic trick" called a Unitary Transformation.

• The Metaphor: Imagine you are watching a movie of the dancers. You can either watch them dance with a spinning background (the twisted phase), or you can change the camera angle so the background looks flat, but the dancers themselves are now spinning.
• The Point: The paper proves that the "twist" in the superconductor (the phase) is mathematically equivalent to a "magnetic field" or a "gauge field" pushing the electrons. This allows the author to calculate the winding number easily by looking at the shadow (the Bloch vector) rather than the complex dance moves.

Why Does This Matter? (The "So What?")

1. Robustness: Topological numbers (like the winding number) are like a knot in a rope. You can wiggle the rope, but you can't untie the knot without cutting the rope. This means the properties of these edge electrons are very stable. They won't disappear just because the material is a little bit dirty or imperfect.
2. Control: The paper shows that we can control these topological states by simply changing the supercurrent or applying a magnetic field. It's like a dial: turn the dial (change the current), and you change the number of twists (the winding number).
3. Future Tech: Because these edge states are so stable, they are perfect candidates for quantum computers. If you can store information in these "twisted" edge states, the computer won't crash easily due to noise or interference.

Summary in a Nutshell

• The Setup: A superconducting ring with a twisted phase (like a giant vortex).
• The Tool: A 3D arrow (Bloch vector) that tracks the "twist" of the electrons.
• The Finding: The amount of twist in the ring determines a "winding number."
• The Twist: There are two types of electrons: those in the middle (bulk) and those on the edge. The edge electrons appear when the bulk electrons hit a wall, and their behavior is dictated by the same winding number.
• The Takeaway: We can control these stable, topological states by adjusting the supercurrent, opening the door to more robust quantum technologies.

In short, the paper explains how spinning a superconductor creates a specific, countable "knot" in the quantum world, and how this knot protects special electrons living on the edge of the material.

Technical Summary

1. Problem Statement

The paper investigates the topological properties of two-dimensional superconducting systems characterized by a spatially periodic order parameter. While standard superconductors have a uniform order parameter phase, certain physical scenarios (such as Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phases, rotating superfluids, or superconducting rings with macroscopic vortices) exhibit a phase that varies periodically in space.

The core problem addressed is:

• How does a spatially modulated order parameter phase affect the topological characterization (specifically the winding number) of the quasiparticle eigenstates described by the Bogoliubov–de Gennes (BdG) Hamiltonian?
• What is the relationship between the periodicity of the order parameter, the Bloch vector, and the emergence of edge states?
• How do boundary conditions (bulk vs. edge) influence the topological invariants in such modulated systems?

2. Methodology

The author employs a theoretical framework based on the SU(2) Hamiltonian formalism, mapping the BdG Hamiltonian to a spin-1/2 system coupled to a pseudomagnetic field.

• Hamiltonian Formulation: The BdG Hamiltonian is expressed as $H_{BdG} = \vec{h} \cdot \vec{\sigma}$, where $\vec{\sigma}$ are Pauli matrices and $\vec{h} = (\Delta', \Delta'', D)^T$. Here, $\Delta = |\Delta|e^{ix/L}$ represents a complex, spatially periodic order parameter, and $D$ is a differential (continuum) or difference (lattice) operator.
• Analytical Solutions: The author solves the eigenvalue problem $\vec{h} \cdot \vec{\sigma} \Psi = E \Psi$ by assuming spinor solutions of the form $\Psi(x) = (a_1 e^{\gamma_1 x}, a_2 e^{\gamma_2 x})^T$. This leads to a quadratic equation for the energy $E$ in terms of the complex wave number $\gamma$.
• Bloch Vector Analysis: The spinor solutions are mapped to a three-dimensional Bloch vector $\vec{s}(x)$, defined as the expectation value of the Pauli vector. The trajectory of this vector on the Bloch sphere is analyzed to determine the winding number ($w$).
• Models Studied:
  1. Continuum Model: $D = -\partial_x^2 - \partial_y^2$.
  2. Tight-Binding Model: $D$ acts as a difference operator on a square lattice ($D\psi_x = \psi_{x+1} + \psi_{x-1} - 2\psi_x$).
• Bulk-Edge Connection: The paper utilizes analytic continuation of the wave number $k$ (where $\gamma = ik$) into the complex plane to identify exponential solutions corresponding to edge modes.
• Unitary Transformations: The author employs unitary transformations to relate the phase of the order parameter to gauge fields and to demonstrate how the winding number is preserved or transformed under specific operations.

3. Key Contributions

• Direct Link Between Order Parameter and Topology: The paper establishes that the periodicity of the order parameter phase ($L$) directly determines the winding number of the eigenfunctions. Specifically, for a system of length $l$ with periodic boundary conditions, the winding number is quantized as $w = n$, where $L = l/2\pi n$.
• Identification of a Topological Invariant: A topological invariant is identified that links the winding number directly to the Bloch vector. The winding number is calculated as the integral of the Berry connection (or the solid angle swept by the Bloch vector) around the north-south axis of the Bloch sphere.
• Distinction Between Bulk and Edge Modes: The study reveals a crucial distinction:
  • Bulk (Plane-wave) Solutions: Exist for purely imaginary $\gamma$ ($\gamma = ik$) and have real energy dispersions.
  • Edge (Exponential) Solutions: Exist for complex $\gamma$ ($\gamma = \kappa + ik$ with $\kappa \neq 0$). These solutions are exponentially localized at the boundaries.
  • Boundary Condition Dependence: The winding numbers for bulk and edge modes can differ because they are subject to different boundary conditions. The paper highlights that the winding number is not an intrinsic property of the Hamiltonian alone but depends on the specific boundary conditions imposed on the system.
• Comparison with $\pi$-Flux Hamiltonian: The author contrasts the BdG system with the $\pi$-flux Hamiltonian. In the BdG system, increasing the order parameter phase gradient increases the winding number. In contrast, for the $\pi$-flux model, the winding number is robust ($w=0, \pm 1$) and depends on the trajectory on the torus rather than a simple phase gradient.

4. Key Results

• Energy Dispersion:
  • For the continuum model, the dispersion relation is derived as $E_{\pm}(\gamma) = \frac{i\gamma}{L} - \frac{1}{2L^2} \pm \frac{1}{2}\sqrt{[\gamma^2 + (\gamma + i/L)^2]^2 + 4|\Delta|^2}$.
  • Real energy solutions exist for both plane waves (bulk) and specific exponential decays (edge).
• Bloch Sphere Trajectories:
  • The Bloch vector traces a closed trajectory on the Bloch sphere with a constant latitude determined by the parameter $b = |E - d(\gamma)|/|\Delta|$.
  • The winding number $w$ corresponds to the number of times this trajectory winds around the north-south axis.
  • Crucial Finding: The spinor periodicity is $4\pi L$, while the Bloch vector periodicity is $2\pi L$. Consequently, the winding number of the Bloch vector is an observable integer, whereas the spinor itself may require a double loop to return to its original state.
• Edge State Emergence:
  • Edge states emerge from the bulk spectrum when the analytic continuation of the wave number leads to complex $\gamma$ values that satisfy the condition of real energy $E$.
  • The existence of these edge modes is strictly tied to the boundary conditions (e.g., a superconducting ring with a Josephson barrier or a macroscopic vortex).
• Gauge Field Interpretation: The spatial modulation of the phase $\phi = x/L$ is shown to be equivalent to the introduction of a gauge field (supercurrent) in the Hamiltonian. The unitary transformation $e^{i\phi/2}$ shifts the derivative operator, effectively creating a Peierls phase in the tight-binding model or a gauge field in the continuum model.

5. Significance

• Control of Topological Phases: The work demonstrates that topological invariants (winding numbers) in superconductors can be controlled macroscopically by manipulating the supercurrent or the order parameter phase (e.g., via magnetic fields or rotation). This offers a potential pathway for engineering topological states without changing the material's microscopic structure.
• Experimental Accessibility: The Bloch vector is identified as an observable quantity. The paper suggests that the first two components of the Bloch vector ($s_1, s_2$) can be measured via linear response to changes in the order parameter ($\partial_{\Delta'} E, \partial_{\Delta''} E$), providing an experimental route to detect these topological phases.
• Fundamental Understanding of Edge States: By clarifying the relationship between bulk and edge winding numbers and the role of boundary conditions, the paper provides a rigorous framework for understanding the emergence of Majorana-like or chiral edge modes in periodically modulated superconductors.
• Future Directions: The paper opens avenues for studying non-Hermitian extensions (non-Bloch BdG Hamiltonians) and systems with random or piecewise linear phase modulations, which could mimic disordered systems or create new types of localized wave functions.

In summary, Ziegler's paper provides a rigorous analytical derivation showing that the spatial modulation of the superconducting order parameter acts as a tunable topological knob, directly dictating the winding number of quasiparticle states and the existence of edge modes, with the Bloch vector serving as the primary observable for this topological characterization.