Braiding rule of boundary Majorana-like zero mode

Basic Information

• Paper ID: 2501.00749
• Title: Braiding rule of boundary Majorana-like zero mode
• Authors: Qiyun Ma, Hailong He, Meng Xiao, Zhengyou Liu
• Classification: cond-mat.mes-hall (Condensed Matter Physics - Mesoscopic and Nanoscale Systems)
• Publication Date: January 2025
• Paper Link: https://arxiv.org/abs/2501.00749

Abstract

The study of topological states has become an important topic in solid-state systems and artificial structures such as photonic and phononic crystals. Among these, Majorana zero modes exhibiting non-trivial braiding processes have attracted considerable research interest. The classical wave analogue of Majorana zero modes—Majorana-like zero modes (MLZMs)—has also recently received significant attention. However, the vast majority of previous work has focused on MLZMs bound to bulk vortex cores. This work reveals the tunable braiding rules of MLZMs around the boundaries of vortex-textured Kekulé-modulated graphene. The study demonstrates that the existence of these zero-dimensional boundary MLZMs is protected by the Zak phase of one-dimensional boundary states. Consequently, multiple MLZMs can be constructed and their corresponding braiding processes analyzed. Furthermore, an implementation scheme for realizing boundary MLZMs in acoustic crystals is provided.

Research Background and Motivation

Problem Definition

The core problems addressed in this research are:

1. Limitations of conventional MLZMs: Existing MLZM research primarily focuses on zero modes at bulk vortex centers, whose positions are difficult to control
2. Missing mechanisms of boundary MLZMs: Although theoretically predicted to exist as complements to vortex-core MLZMs, the detailed mechanisms, wavefunctions, and braiding rules of boundary zero modes remain unstudied
3. Demand for tunability: Development of position-tunable MLZMs is needed to extend applications of topological state braiding

Significance

1. Quantum computing applications: The non-Abelian braiding properties of MLZMs are significant for topological quantum computation
2. Extension to artificial structures: Realization of analogous functionality in classical systems such as photonic and phononic crystals
3. Theoretical refinement: Enhancement of understanding regarding topological zero modes in vortex-textured modulated systems

Limitations of Existing Methods

1. All MLZMs are bound to vortex centers with fixed, non-tunable positions
2. Lack of systematic study of boundary MLZMs (BMLZMs)
3. Braiding rules limited to interactions between vortex-core MLZMs

Core Contributions

1. First systematic study of boundary MLZMs: Reveals the existence and properties of tunable MLZMs at the boundaries of vortex-textured Kekulé-modulated graphene
2. Establishment of theoretical framework: Demonstrates that BMLZMs are zero-dimensional domain wall states protected by the Zak phase of one-dimensional boundary states
3. Derivation of braiding rules: Establishes complete braiding matrix theory encompassing both center MLZMs (CMLZMs) and boundary MLZMs
4. Proposed implementation scheme: Designs an acoustic crystal realization scheme, providing a feasible path for experimental verification
5. Extended application scenarios: Provides new degrees of freedom for topological state braiding through position tunability

Detailed Methodology

Theoretical Model

Tight-binding Hamiltonian

The research is based on a tight-binding model of graphene with vortex-textured Kekulé modulation:

H = Σ_r∈Λ_a Σ_{i=1}^3 [t + δt_{r,i}] a_r† b_{r+s_i} + H.c.

Where:

• δt_{r,i} = Δ_0 e^{i(K_+·s_i + (K_+ - K_-)·r + φ_r)/3} + c.c.
• φ_r = θ_r + α_0 contains the azimuthal angle θ_r and initial offset α_0
• The system employs armchair boundaries to avoid trivial zero-energy modes

Effective k·p Hamiltonian

The effective Hamiltonian in the low-energy limit takes a form analogous to the Bogoliubov-de Gennes (BdG) Hamiltonian of topological superconductors:

H = [u_b†, u_a†, v_a†, v_b†] × 4×4 matrix × [u_b, u_a, v_a, v_b]^T

This establishes the connection between MLZMs and genuine Majorana zero modes.

Boundary Zero Mode Mechanism

Zak Phase Protection Mechanism

Through analysis of the band structure of one-dimensional boundary states, the following is discovered:

1. Boundary states form one-dimensional Dirac points at φ = ±0.873π
2. The Zak phase of boundary bands undergoes a transition (0,0) → (π,π) upon crossing the Dirac point
3. This topological phase transition leads to the emergence of domain wall states, analogous to the Jackiw-Rebbi mechanism

Analytical Form of Zero Mode Solutions

The analytical solutions for boundary zero modes are:

u_a = exp[Δ_0 cos(φ_r-φ')/vt (x-x_0)] exp[-Δ_0 sin(φ_r-φ')/vt (y-y_0)]
u_b = exp[-Δ_0 cos(φ_r-φ')/vt (x-x_0)] exp[-Δ_0 sin(φ_r-φ')/vt (y-y_0)]

Derivation of Braiding Rules

Phase Evolution of Individual BMLZMs

As the initial offset α_0 varies from 0 to 2π:

1. The BMLZM rotates clockwise once along the boundary
2. A phase of π is acquired each time a lattice corner is crossed
3. After complete rotation, a total phase of 3π is acquired; accounting for 2π redundancy, the net acquired phase is π

Braiding Matrix for Multiple MLZMs

For n BMLZMs, the braiding matrix is:

U_n = Σ [U_n]_{i,j} |i⟩⟨j|

Where [U_n]_{1,n} = (-1)^{n_1}, [U_n]_{i,i-1} = (-1)^{n_i}, with n_i recording the number of times the i-th BMLZM crosses a corner.

Experimental Setup

Numerical Calculation Parameters

• Lattice size: Number of boundary lattice points N_b = 18 (Figure 1)
• Hopping strength: t = 1, Δ_0 = -1.5
• Boundary type: Three armchair boundaries with outward normals parallel to ΓK_+
• System geometry: Equilateral triangle with vortex core at geometric center

Computational Methods

1. Eigenvalue problem solving: Numerical diagonalization of tight-binding Hamiltonian
2. Band structure calculation: Application of Bloch periodic boundary conditions
3. Zak phase calculation: Berry phase integration to compute boundary state Zak phases
4. Wavefunction analysis: Extraction of spatial distribution of zero-energy eigenstates

Verification Scheme

• Comparison with analytical k·p model results
• Scaling behavior examination across different system sizes
• Acoustic crystal implementation scheme design (Supplementary Material Section IX)

Experimental Results

Main Findings

Zero Mode Distribution Characteristics

1. CMLZM: Always located at the vortex core, distributed only on sublattice Λ_b, position independent of α_0
2. BMLZM: Located at the boundary where φ_r = 0.873π, distributed only on sublattice Λ_a, position tunable with α_0

Topological Phase Transition of Boundary States

• One-dimensional Dirac point exists in boundary states at φ = 0.873π
• Zak phase undergoes transition at this point: (0,0) → (π,π)
• Domain wall states emerge at boundaries with different φ values

Braiding Behavior Verification

Braiding of two BMLZMs:

• Initial state: |ψ_0⟩ = 1,1^T
• Post-braiding state: |ψ_T⟩ = U_2|ψ_0⟩ = [0,1; -1,0]1,1^T = -1,1^T
• Braiding matrix: U_2 = 0,1; -1,0

Multi-BMLZM configurations:

• Number of BMLZMs on different boundaries affects braiding period and final state
• Braiding matrix form: Non-zero elements of U_n are U_n{1,n} = (-1)^{n_1} and U_n{i,i-1} = (-1)^{n_i}

Numerical Verification Results

1. Zero energy precision: Numerically computed zero mode energy at 10^{-10} level
2. Wavefunction localization: BMLZM and CMLZM exhibit exponential localization at boundary and center respectively
3. Phase evolution: Numerical results for BMLZM position variation with α_0 perfectly match theoretical predictions

Related Work

Majorana Zero Mode Research

• Solid-state systems: MZMs in p-wave superconductors Read & Green, 2000
• Artificial structures: MLZMs in photonic crystals Iadecola et al., 2016 and phononic crystals Gao et al., 2019
• Braiding theory: Theoretical foundations of non-Abelian braiding Ivanov, 2001; Nayak et al., 2008

Graphene Topological States

• Kekulé modulation: Vortex-textured modulation producing fractional charge zero modes proposed by Hou et al. (2007)
• Photonic realization: MLZMs realized in photonic crystals by multiple research groups Noh et al., 2020; Yang et al., 2022
• Acoustic realization: Topological zero modes in phononic crystals Wu et al., 2021

Novelty of This Work

Compared with existing work, this paper is the first to:

1. Systematically study the existence and properties of boundary MLZMs
2. Establish the Zak phase protection mechanism for BMLZMs
3. Derive complete braiding theory encompassing both CMLZMs and BMLZMs
4. Propose position-tunable MLZM realization schemes

Conclusions and Discussion

Main Conclusions

1. Existence of boundary MLZMs: Demonstrates the existence of zero-dimensional topological zero modes protected by Zak phase at boundaries of vortex-textured Kekulé-modulated graphene
2. Tunability: BMLZM positions can be continuously adjusted by varying the initial offset α_0
3. Braiding rules: Establishes a complete braiding theory framework incorporating corner point phase jumps
4. Implementation pathways: Provides realization schemes in artificial structures such as acoustic crystals

Limitations

1. Theoretical model: Based on tight-binding approximation, neglecting higher-order effects
2. Boundary conditions: Primarily considers armchair boundaries; other boundary types require further investigation
3. Experimental verification: Lacks direct experimental observation verification
4. Interactions: Does not account for electron-electron interactions and other many-body effects

Future Directions

1. Experimental realization: Verification of theoretical predictions in acoustic and optical systems
2. Boundary engineering: Investigation of effects of different boundary types on BMLZMs
3. Many-body effects: Consideration of interaction effects on MLZM properties
4. Quantum applications: Exploration of potential applications in quantum information processing

In-Depth Evaluation

Strengths

1. Theoretical completeness: Establishes a complete theoretical framework from tight-binding model to effective theory
2. Physical insight: Reveals the Zak phase protection mechanism of BMLZMs, providing profound physical understanding
3. Numerical verification: Detailed numerical calculations verify theoretical predictions
4. Practical value: Tunability offers new perspectives for topological quantum device design
5. Clear presentation: Well-structured paper with intuitive physical imagery

Weaknesses

1. Experimental absence: Lacks experimental verification; theoretical predictions await experimental confirmation
2. Model limitations: Limited applicability range of tight-binding model
3. Boundary dependence: Primarily addresses specific boundary types; universality requires verification
4. Application outlook: While application prospects are proposed, specific implementation pathways require further exploration

Impact

1. Academic value: Provides new theoretical framework and physical mechanisms for topological zero mode research
2. Application potential: Possesses potential applications in quantum computing and topological device design
3. Cross-disciplinary impact: Generalizable to multiple artificial structure fields including photonics and acoustics
4. Subsequent research: Provides important reference for related theoretical and experimental investigations

Applicable Scenarios

1. Artificial crystal design: Topological zero mode design in photonic and phononic crystals
2. Quantum simulation: Simulation of Majorana zero modes in cold atom systems
3. Topological devices: Development of quantum devices based on tunable zero modes
4. Fundamental research: Theoretical investigation of topological phase transitions and domain wall states

References

This paper cites 61 important references spanning multiple fields including topological physics, Majorana zero modes, graphene physics, and artificial structures, providing a solid theoretical foundation for the research.