A One-Particle Density Matrix Framework for Mode-Shell Correspondence: Characterizing Topology in Amorphous Higher-Order Topological Insulators

Basic Information

• Paper ID: 2509.03632
• Title: A One-Particle Density Matrix Framework for Mode-Shell Correspondence: Characterizing Topology in Amorphous Higher-Order Topological Insulators
• Authors: Miguel F. Martínez, Lucien Jezequel, Jens H. Bardarson, Thomas Klein Kvorning, Julia D. Hannukainen
• Classification: cond-mat.mes-hall cond-mat.dis-nn
• Submission Date: September 2024 (arXiv)
• Paper Link: https://arxiv.org/abs/2509.03632

Abstract

This paper proposes a framework for directly characterizing higher-order topological phases from the one-particle density matrix without reference to the underlying Hamiltonian. The method extends the mode-shell correspondence, originally formulated for single-particle Hamiltonians, to Gaussian states subject to chiral constraints. In this correspondence, the mode index counts topological boundary modes, while the shell index quantifies the bulk topology surrounding a mode, providing a bulk-boundary diagnostic. In one-dimensional topological insulators, the shell index reduces to a local chiral marker that recovers the winding number in the translation-invariant limit. The authors apply the mode-shell correspondence to $C_4$-symmetric higher-order topological insulators with chiral constraints and demonstrate that fractional shell indices imply the higher-order phase is intrinsic.

Research Background and Motivation

Problem Definition

The characterization of topological quantum states traditionally relies on momentum-space analysis of the Hamiltonian, but this approach faces fundamental difficulties when dealing with disordered systems and amorphous materials. Specifically:

1. Absence of Translational Invariance: In amorphous materials, the conventional concept of the Brillouin zone fails, and topological invariants based on momentum space cannot be defined
2. Complexity of Higher-Order Topological Phases: Corner modes and hinge modes in higher-order topological insulators require novel characterization methods
3. Distinguishing Intrinsic from Extrinsic Phases: It is necessary to determine whether the topological phase is guaranteed by bulk topology (intrinsic) or depends on specific boundary termination (extrinsic)

Research Significance

This research addresses important problems in condensed matter physics:

• Provides characterization tools for amorphous topological materials
• Establishes methods for analyzing topological properties directly from quantum states
• Extends the applicability of topological physics to disordered systems

Limitations of Existing Methods

Traditional approaches such as nested Wilson loops and quantized multipole moments all depend on translational symmetry and cannot be applied to:

• Amorphous materials
• Strongly disordered systems
• Certain cases of interacting systems

Core Contributions

1. Theoretical Framework Innovation: Generalizes the mode-shell correspondence from Hamiltonian formulation to one-particle density matrix formulation
2. Real-Space Topological Characterization: Provides a method for characterizing topological phases that does not depend on translational invariance
3. Intrinsicness Criterion: Establishes criteria based on fractional shell indices for determining whether higher-order topological phases are intrinsic
4. Application to Amorphous Systems: First verification of the method in amorphous higher-order topological insulators with $C_4$ symmetry
5. Extension to Interacting Systems: Demonstrates that the method can be extended to interacting systems with gapped one-particle density matrices

Detailed Methodology

One-Particle Density Matrix Framework

For a many-body state $|\Psi\rangle$, the Bogoliubov-de Gennes one-particle density matrix has the block form:

$\rho = \begin{pmatrix} \tilde{\rho} & \kappa \\ \kappa^\dagger & 1-\tilde{\rho}^* \end{pmatrix}$

where $\tilde{\rho}_{ij} = \langle\Psi|c_i^\dagger c_j|\Psi\rangle$ and $\kappa_{ij} = \langle\Psi|c_i^\dagger c_j^\dagger|\Psi\rangle$.

Mode-Shell Correspondence

Mode Index Definition

The mode index is defined as: $I^{\text{mode}} = \sum_\lambda 4\lambda(1-\lambda)\langle\lambda|\theta C|\lambda\rangle$

where:

• $\lambda$ are eigenvalues of the restricted one-particle density matrix $\rho_A$
• $|\lambda\rangle$ are the corresponding eigenstates
• $\theta$ is the projector onto the subregion $A_{\text{shell}} \subset A$
• $C$ is the chiral symmetry operator

Shell Index Derivation

Through the cyclic property of the trace and the chiral constraint $\{\rho_A, C\} = C$, the mode index can be rewritten as: $I^{\text{shell}} = -2\text{Tr}(C\rho_A[\rho_A, \theta])$

Under gradient expansion: $I^{\text{shell}} = 2\sum_i \text{Tr}(C\rho_A[\rho_A, X_i]\partial_{x_i}\theta)$

Equivalence in One Dimension

In one-dimensional systems, the shell index is equivalent to the local chiral marker: $\nu(x) = -2\text{tr}\langle x|\rho CX\rho|x\rangle$

In the bulk region: $I^{\text{shell}} = \mp\nu$

Application to Higher-Order Topological Insulators

BBH Model

Consider the real-space Hamiltonian of the Benalcazar-Bernevig-Hughes (BBH) model: $h_{rr'} = -\gamma(\sigma_z+\sigma_y\tau_y)\delta_{rr'} - \frac{\eta}{2}(\sigma_z + i\sigma_y\tau_z)\delta_{r',r+a_x} - \frac{\eta}{2}(\sigma_y\tau_y + i\sigma_y\tau_x)\delta_{r',r+a_y}$

Intrinsicness Criterion

For $C_4$-symmetric higher-order topological insulators:

• Odd Mode Index: Each edge-shell crossing point contributes a half-integer value, indicating an intrinsic topological phase
• Even Mode Index: Integer-valued contributions cannot distinguish between intrinsic and extrinsic phases

Experimental Setup

Crystalline System Parameters

• Square lattice linear size: $L = 20$
• BBH model parameters: $\gamma = 0.5$, $\eta = 1$
• Restricted region $A$: lower-left quadrant of the lattice ($L_A = 10$)

Amorphous System Generation

1. Generate small square lattice ($\tilde{L} = L/2$)
2. Introduce disorder degree $w$ (amorphousness parameter) via normal distribution
3. Create four copies and apply $C_4$ rotational symmetry
4. Set parameters: $w = 0.1$, cutoff distance $r_0 = 1.3$

Shell Function Design

$\theta(x,y) = 1 - [1 + e^{-c(\sqrt{x^2+y^2}-\omega)}]^{-1}$

where $\omega = L_A/2$ determines the shell position, and $c = L_A/2$ (crystalline) or $c = 2L_A/3$ (amorphous) controls the transition sharpness.

Experimental Results

Crystalline BBH Model Results

• Mode Index: $I^{\text{mode}} = 0.998 \approx 1$
• Shell Index: $I^{\text{shell}} = 0.998 \approx 1$
• Four $\lambda = 1/2$ Eigenvalues: Corresponding to four chiral corner modes
• Shell Index Distribution: Contributions only at edge-shell crossing points, each contributing $1/2$

Amorphous BBH Model Results

• Mode Index: $I^{\text{mode}} = 0.996 \approx 1$
• Shell Index: $I^{\text{shell}} = 0.996 \approx 1$
• Structural Disorder Effects: More eigenvalues appear in the $(0,1)$ interval
• Bulk Contributions: Shell index shows positive and negative contributions in the bulk region, but averages to zero

Key Findings

1. Quantization Accuracy: The mode-shell correspondence maintains good quantization even in amorphous systems
2. Intrinsicness Confirmation: Fractional shell indices ($1/2$) confirm the intrinsic higher-order topological nature of the BBH model
3. Robustness: The method is robust against structural disorder

Related Work

Real-Space Topological Markers

• Local Chiral Marker: One-dimensional real-space topological invariant proposed by Hannukainen et al.
• Bott Index: Real-space topological marker for disordered systems
• Other Real-Space Methods: Including spectral localization and local Chern markers

Higher-Order Topological Insulators

• Nested Wilson Loops: Traditional method for characterizing higher-order topology
• Multipole Moments: Topological invariants based on Wannier centers
• Symmetry Indicators: Topological classification based on band representations

Amorphous Topological Materials

• Amorphous Chern Insulators: Earliest studies of amorphous topological materials
• Amorphous Higher-Order Topological Insulators: Recently developed research direction

Conclusions and Discussion

Main Conclusions

1. Theoretical Breakthrough: Successfully generalizes the mode-shell correspondence to the one-particle density matrix formulation of quantum states
2. Practical Verification: Validates the method's effectiveness in both crystalline and amorphous higher-order topological insulators
3. Intrinsicness Criterion: Establishes criteria for intrinsic higher-order topological phases based on fractional shell indices
4. Broad Applicability: The method can be extended to interacting systems and arbitrary dimensions

Limitations

1. Chiral Constraint Requirement: The method only applies to systems with chiral symmetry
2. Gap Requirement: Requires the one-particle density matrix to be gapped
3. Shell Selection: The choice of shell region requires careful consideration of relevant length scales
4. Computational Complexity: Diagonalization of the one-particle density matrix may be computationally expensive for large systems

Future Directions

1. Other Symmetry Classes: Extension to other Altland-Zirnbauer symmetry classes
2. Three-Dimensional Systems: Application to three-dimensional higher-order topological insulators
3. Strongly Interacting Systems: Further exploration of applications in strongly interacting systems
4. Experimental Verification: Verification of the method in real materials

In-Depth Evaluation

Strengths

1. Strong Theoretical Innovation: First generalization of mode-shell correspondence to density matrix formulation, opening new avenues for amorphous topological material research
2. Good Method Universality: Does not depend on translational invariance and can be applied to a wide range of system types
3. Clear Physical Picture: The correspondence between mode and shell indices provides intuitive understanding of bulk-boundary relationships
4. Sufficient Numerical Verification: Detailed numerical verification in both crystalline and amorphous systems
5. Novel Intrinsicness Criterion: Proposes a new method for determining whether higher-order topological phases are intrinsic

Weaknesses

1. Symmetry Restrictions: Method is limited to systems with chiral symmetry, restricting its applicability
2. Computational Cost: High computational complexity for density matrix diagonalization in large systems
3. Parameter Sensitivity: The choice of shell function $\theta$ may affect the accuracy of results
4. Limited Experimental Connection: Lacks direct connection to measurements in real materials

Impact and Significance

1. Academic Contribution: Provides important tools for theoretical development in topological physics
2. Application Prospects: Provides theoretical foundation for design and characterization of amorphous topological materials
3. Method Generalization: May inspire development of other real-space topological characterization methods
4. Interdisciplinary Value: Method may have applications in quantum information and materials science

Applicable Scenarios

1. Amorphous Topological Materials: Systems with high disorder but retaining certain symmetries
2. Finite-Size Systems: Small topological systems where boundary effects are important
3. Interacting Systems: Strongly correlated systems with gapped one-particle density matrices
4. Numerical Studies: Theoretical calculations requiring real-space topological characterization

References

This paper cites 98 references, covering core literature in multiple important research fields including topological insulators, real-space topological markers, amorphous materials, and higher-order topological phases, providing a solid theoretical foundation for this research.