Majorana Flat Bands in the Vortex Line of Superconducting Weyl Semimetals

Imagine you have a magical, three-dimensional crystal that acts like a highway for electrons. In this crystal, the electrons don't just flow; they behave like massless particles called Weyl fermions, zipping around in specific directions. This is a "Weyl Semimetal."

Now, imagine you take this crystal and cool it down until it becomes a superconductor (a material where electricity flows with zero resistance). When you poke a hole in this superconductor with a magnetic field, you create a tiny, swirling tunnel of energy called a vortex line.

This paper is about what happens to the electrons inside that tiny tunnel. The authors discovered something amazing: Majorana Flat Bands.

Here is the story of how they found it, explained with simple analogies.

1. The "Stack of Pancakes" Trick

The Weyl Semimetal is a 3D object, which is hard to analyze all at once. The authors had a clever idea: they treated the 3D crystal like a stack of 2D pancakes.

The Analogy: Imagine a loaf of bread. If you slice it, each slice is a 2D piece of bread. In this crystal, if you slice it at different heights (let's call the height $k_z$ ), each slice behaves like a different kind of 2D material.
Some slices are "normal" (boring).
Some slices are "topological" (special, with a twist).
The "twist" in these special slices is measured by something called a Chern number. Think of the Chern number as a "knot count." A knot count of 1 means the slice is topologically special.

2. The Magic Tunnel (The Vortex)

When the crystal becomes a superconductor and you create a vortex line (the tunnel), you are essentially poking a hole through the center of this stack of pancakes.

The Rule: Physics tells us that if you poke a hole through a "knotted" pancake (one with an odd Chern number, like 1), a special ghost-like particle appears in the hole. This is a Majorana Zero Mode (MZM). It's a particle that is its own antiparticle.
The Discovery: The authors found that in their crystal, there is a whole range of pancake slices (a range of heights) that are "knotted." Because every slice in that range has a knot, a Majorana ghost appears in the tunnel at every height in that range.
The Result: Instead of just one ghost particle, you get a continuous flat band of them. Imagine a flat, invisible highway running up the entire length of the tunnel where these ghost particles can live. This is the Majorana Flat Band (MFB).

3. Tuning the Radio

The authors showed that you can control where this "ghost highway" appears.

The Analogy: Think of the crystal like a radio. You can turn two knobs:
1. Chemical Potential ( $\mu$ ): How many electrons are in the system.
2. Pairing Strength ( $\Delta$ ): How strongly the electrons want to pair up (superconduct).
By turning these knobs, you can make the "ghost highway" appear along the entire length of the tunnel, or disappear completely. It's like tuning the radio to find the perfect station where the magic happens.

4. The "Edge vs. Core" Mix-Up

Here is a tricky part the authors solved.

The Expectation: They calculated exactly where the "ghost highway" should start and stop based on their math (the "phase diagram").
The Reality: When they ran computer simulations, the highway didn't quite reach the exact edges they predicted. It was slightly shorter.
The Explanation: They realized this was due to hybridization (mixing).
- Imagine the tunnel has a Core (the very center) and an Edge (the walls of the tunnel).
- There is a ghost particle in the Core and a ghost particle on the Edge.
- When the tunnel is very narrow (or the simulation is small), these two ghosts can "see" each other and mix up, which pushes them slightly apart in energy, making the "flat" band look a little bumpy or short.
- The Fix: The authors showed that if you make the tunnel wider (increase the lattice size), the Core and Edge ghosts get far apart, stop mixing, and the highway becomes perfectly flat and reaches the exact mathematical boundaries.

5. Making it Real (The Recipe)

Finally, the paper asks: "How do we actually build this?" You can't just wave a wand.

The authors proposed a recipe using a specific type of interaction between electrons (called Hubbard interaction).
They showed that if you take a Weyl Semimetal and apply the right amount of "attraction" between electrons (using a mean-field approach), nature will spontaneously create the superconducting state with the right properties.
They mapped out a "phase diagram" (a map of ingredients) showing exactly how much attraction and how many electrons you need to get this magical state.

Summary

In simple terms, this paper says:

We found a way to create a superconducting crystal with a special internal structure.
When we poke a magnetic tunnel through it, a continuous line of exotic "ghost" particles (Majorana fermions) appears inside the tunnel.
We can control exactly where this line appears by adjusting the electron density and superconducting strength.
We explained why the line sometimes looks a bit fuzzy at the edges (due to particles mixing) and how to fix it.
We provided a theoretical "recipe" to create this state in real materials.

This is exciting because these Majorana particles are the holy grail of quantum computing. They are incredibly stable and could be used to build computers that don't crash easily (fault-tolerant quantum computers). Finding a way to create a whole band of them, rather than just one, makes them much easier to find and use.

Here is a detailed technical summary of the paper "Majorana Flat Bands in the Vortex Line of Superconducting Weyl Semimetals" by Zhicheng Zhang and Kou-Han Ma.

1. Problem Statement

The paper addresses the nature of vortex bound states (VBSs) in superconducting (SC) time-reversal-symmetry-breaking Weyl semimetals. While previous studies have identified Majorana zero modes (MZMs) at the ends of vortex lines in time-reversal-symmetric systems or gapless chiral Majorana modes in specific pair-density-wave states, the behavior of VBSs in time-reversal-breaking SC Weyl semimetals remains underexplored. Specifically, the authors investigate whether Majorana Flat Bands (MFBs)—zero-energy flat bands in the Bogoliubov-de Gennes (BdG) spectrum—emerge within the vortex line of these materials and seek to determine their exact boundaries and topological characterization.

2. Methodology

The authors employ a multi-faceted approach combining analytical theory, numerical simulations, and mean-field approximations:

Dimensional Decomposition: They treat the 3D Weyl semimetal as a stack of 2D Chern insulators with varying Chern numbers along the $k_z$ direction. By fixing $k_z$ , the 3D problem is decomposed into a series of 2D SC Chern insulators.
Topological Phase Diagram Analysis: They analytically and numerically calculate the topological phase diagram of the constituent 2D SC Chern insulators. This involves determining BdG Chern numbers ( $C$ ) and identifying phase boundaries where bulk gaps close.
Vortex Line Modeling: To simulate the vortex line, they introduce a $\pi$ -flux (phase $e^{i\theta}$ ) into the pairing potential along the $z$ -axis. They apply open boundary conditions in the $x$ - $y$ plane (to capture edge states) and periodic boundary conditions in $z$ .
Microscopic Realization: To move beyond phenomenological pairing, they introduce an attractive Hubbard interaction to the Weyl semimetal model. Using mean-field theory and a folded Brillouin zone (FBZ) approach (to handle coexisting BCS and Fulde-Ferrell (FF) pairing momenta), they derive self-consistent gap equations to verify the emergence of BCS pairing under appropriate parameters.
Topological Invariants: They propose a $k_z$ -dependent $\mathbb{Z}_2$ Chern-Simons invariant to characterize the presence of MZMs in the vortex core for each $k_z$ slice.

3. Key Contributions

A. Mechanism of MFB Emergence

The paper establishes that MFBs arise when the BdG Chern number of the 2D SC Chern insulator slice (at a fixed $k_z$ ) is odd (specifically $C=1$ ).

In a chiral topological superconductor with an odd Chern number, a vortex core hosts a single MZM.
As $k_z$ varies, a continuous range of $k_z$ values exists where the system maintains $C=1$ . The collection of MZMs across this range forms a Majorana Flat Band along the vortex line.
The boundaries of the MFB are determined by the gap-closing conditions of the bulk spectrum, specifically where $\cos k_z = \pm \Delta$ (for $\mu=0$ ).

B. Discrepancy Between Analytical and Numerical Boundaries

A critical finding is the discrepancy between the theoretical phase boundaries (derived from the bulk topological phase diagram) and the numerically observed MFB regions.

Cause: The MFBs do not extend exactly to the analytical boundaries due to hybridization between the vortex-core MZM and the edge-localized MZM.
Mechanism: Near the phase transition points, the edge states become delocalized, increasing their overlap with the vortex core. This hybridization lifts the degeneracy of the two MZMs, opening a finite energy gap and causing the MFB to "shrink" slightly away from the exact phase boundaries.
Verification: The authors demonstrate that increasing the lattice size reduces this hybridization exponentially, causing the numerically observed MFB region to converge toward the analytical boundaries.

C. Control Parameters ( $\mu$ and $\Delta$ )

The study shows that MFBs can be tuned to exist along the entire $k_z$ axis by adjusting:

Pairing Strength ( $\Delta$ ): Increasing $\Delta$ expands the $C=1$ region in the phase diagram.
Chemical Potential ( $\mu$ ): Tuning $\mu$ deforms the phase boundaries (hyperbolas in the $\Delta$ - $m$ plane), allowing the system to traverse the $C=1$ region across the full $k_z$ range.

D. Microscopic Realization via Hubbard Model

The authors prove that the required BCS pairing is not merely an ad-hoc addition but can emerge naturally from an attractive Hubbard interaction.

They analyze the competition between BCS pairing ( $q=0$ ) and Fulde-Ferrell (FF) pairing ( $q=2Q$ ).
They identify a phase diagram in the $\mu-U$ space where the BCS phase dominates (essential for the MFBs discussed) at higher chemical potentials, while FF pairing dominates at lower potentials or very strong interactions.

E. New Topological Indicator

They propose a $k_z$ -resolved $\mathbb{Z}_2$ Chern-Simons invariant defined as $\nu(k_z) = C(k_z) \mod 2$ . This invariant serves as a direct topological indicator for the existence of MZMs in the vortex core at a specific $k_z$ , distinguishing the MFB phase from trivial phases.

4. Results

Existence: MFBs are confirmed to exist in the vortex lines of SC time-reversal-breaking Weyl semimetals.
Boundaries: The exact boundaries are analytically derived as solutions to $\cos k_z = \pm \Delta$ (for $\mu=0$ ), but numerically observed to be slightly narrower due to edge-core hybridization.
Tunability: The MFBs can be extended to cover the entire Brillouin zone in the $k_z$ direction by tuning $\mu$ or $\Delta$ .
Phase Diagram: A comprehensive phase diagram for the SC Chern insulator is constructed, mapping regions of $C=0, \pm 1, \pm 2$ .
Realization: Mean-field calculations confirm that a Weyl semimetal with attractive Hubbard interactions naturally enters a BCS superconducting phase suitable for hosting these MFBs under specific parameter regimes.

5. Significance

Theoretical Insight: The work provides a clear physical picture of how 3D topological superconductivity emerges from stacked 2D topological layers, linking bulk topological invariants to vortex line properties.
Experimental Relevance: By identifying the conditions (tuning $\mu$ and $\Delta$ ) to realize MFBs along the entire vortex line, the paper offers a roadmap for experimental detection in candidate materials (e.g., via point contact or ion sputtering techniques mentioned in the introduction).
Resolution of Hybridization Effects: The detailed analysis of the edge-core hybridization explains why numerical simulations often show gaps near phase boundaries, resolving a common confusion in topological superconductor studies.
Topological Classification: The introduction of the $k_z$ -dependent Chern-Simons invariant provides a robust tool for characterizing 1D topological defects (vortex lines) in 3D systems, distinguishing them from the 2D surface states or 0D end modes studied previously.

In conclusion, this paper establishes a rigorous framework for understanding and engineering Majorana flat bands in superconducting Weyl semimetals, bridging the gap between abstract topological phase diagrams and realistic microscopic models.