Continuous symmetry analysis and systematic identification of candidate order parameters for interacting fermion models

This paper presents a systematic framework that utilizes Majorana representation and Lie algebra theory to analyze continuous symmetries and exhaustively identify candidate order parameters in interacting fermion models with multiple internal degrees of freedom.

The Gist

Imagine you are trying to understand a massive, chaotic dance floor filled with thousands of dancers (the electrons). Some dancers are spinning, some are jumping, some are holding hands in pairs, and others are moving in perfect sync. In the world of physics, this dance floor is a material, and the "dance moves" are the laws of physics that govern how these particles interact.

The big question physicists ask is: What is the underlying rhythm of this dance? And more importantly, what happens if the dancers suddenly stop dancing in sync and start forming a new pattern? This new pattern is called an "ordered phase," and the thing that describes the change is called an "order parameter."

For simple dances, it's easy to spot the rhythm. But for complex materials with many types of dancers (spin, layers, valleys), the rhythm is hidden, and finding the new patterns is like trying to find a specific needle in a haystack made of other needles.

This paper presents a new, systematic "searchlight" method to find these hidden rhythms and predict all possible new patterns the dancers could form.

The Three-Step Magic Trick

The authors developed a three-step algorithm to solve this puzzle, which they tested on two specific types of "dance floors" (models of electrons on a honeycomb lattice, like graphene).

1. Translating the Language (The Majorana Switch)

Imagine the dancers are speaking a complex, confusing language (standard quantum mechanics). To understand them, the authors first translate everything into a simpler, more universal language called Majorana fermions.

• The Analogy: Think of this as translating a complex poem into a simple code of binary (0s and 1s). In this new language, the "dance moves" (symmetries) become simple geometric rotations. If the dancers can rotate in a certain way without changing the overall vibe of the party, that rotation is a "symmetry."

2. Mapping the Rhythm (Finding the Lie Algebra)

Once the language is translated, the authors look for the "generators" of the dance. These are the basic moves that, when combined, create all the possible symmetries.

• The Analogy: Imagine you have a giant box of Lego bricks. You want to know what structures you can build. The authors don't just guess; they use a mathematical rulebook (Lie Algebra theory) to list every single unique way the bricks can snap together.
• They found that for the first model (the Hubbard model), the hidden rhythm was SO(4) (a complex 4D rotation).
• For the second, more complex model (a double-layer system), they discovered a hidden rhythm called Spin(5) × U(1). This was a surprise! It's like finding out that a group of people who thought they were just dancing in a circle were actually performing a secret, synchronized 5-dimensional acrobatic routine.

3. Predicting the New Patterns (Order Parameters)

Now that they know the rules of the dance (the symmetries), they need to predict what happens if the music stops or changes. What new formation will the dancers take?

• The Analogy: Imagine you know the rules of a game of chess. You can now predict all the possible winning moves. The authors used a computer algorithm to systematically list every single possible way the electrons could rearrange themselves to break the symmetry.
• For the simple model, they found 7 possible new patterns.
• For the complex double-layer model, they found 18 distinct patterns.

Why Does This Matter?

In the past, physicists often had to "guess" what the new patterns might be based on intuition. If they guessed wrong, they might miss a new state of matter entirely.

This paper is like giving physicists a complete catalog of every possible outcome.

• For the "Double-Layer" Model: The authors found that between the normal state and a "symmetric mass generation" (a weird state where particles get heavy without breaking symmetry), there is an intermediate "excitonic" phase.
• The Analogy: Think of water freezing. Usually, it goes from Liquid → Ice. But this paper predicted that for this specific material, there is a "slush" phase in between. Because they had the full list of 18 patterns, they could point directly to the specific pattern (the "slush") and say, "This is what happens next."

The Takeaway

This work is a systematic toolkit. Instead of relying on a "gut feeling" to find hidden symmetries in complex quantum materials, scientists can now plug their model into this framework, and it will mathematically spit out:

1. The exact hidden symmetries.
2. A complete list of every possible new phase of matter the material could enter.

It turns the search for new quantum states from a game of "hide and seek" into a game of "checklist," ensuring that no hidden treasure is left undiscovered.

Technical Summary

1. Problem Statement

In condensed matter physics, understanding phase transitions in interacting fermion systems relies heavily on identifying the system's continuous symmetry group and the corresponding order parameters that characterize spontaneous symmetry breaking (SSB).

• The Challenge: While simple models allow for symmetry identification by inspection, complex many-body systems with multiple internal degrees of freedom (spin, orbital, valley, layer) present formidable difficulties. These systems often exhibit enlarged or "hidden" emergent symmetries that are not apparent in the standard complex fermion basis.
• The Gap: There is a lack of a systematic, algorithmic framework to exhaustively enumerate all possible candidate order parameters (bilinear and higher-order) for such complex Hamiltonians. Without a complete classification, it is difficult to rule out competing phases or establish phenomena like Symmetric Mass Generation (SMG), which requires proving the absence of long-range order in all symmetry-breaking channels.

2. Methodology

The authors propose a two-stage algorithmic framework that transforms the physical problem of symmetry analysis into a rigorous algebraic problem.

Stage I: Continuous Symmetry Analysis via Majorana Representation

1. Majorana Mapping: The fermionic Hamiltonian $H$ is mapped to a representation using Majorana fermion operators ($\gamma$). This treats particle-hole degrees of freedom on equal footing with spin and orbital indices.
2. Lie Algebra Identification: Continuous symmetries in the Majorana basis correspond to orthogonal transformations. The infinitesimal generators of these symmetries are antisymmetric matrices that commute with the Hamiltonian tensor.
   • The problem reduces to solving the linear system $[H, X] = 0$ for $X \in \mathfrak{so}(2N_f)$, where $N_f$ is the number of local fermion modes.
3. Structural Classification: Once the Lie algebra $\mathfrak{g}$ is found, the authors employ the theory of semisimple Lie algebras to classify it:
   • Decompose $\mathfrak{g}$ into ideals.
   • Identify the Cartan subalgebra and compute root systems.
   • Construct Dynkin diagrams to determine the algebraic structure (e.g., $su(2)$, $so(5)$).
   • Determine the faithful symmetry group $G$ by analyzing the action of the center of the universal covering group on the physical Hilbert space (accounting for quotient groups like $Z_2$).

Stage II: Systematic Identification of Order Parameters

1. Exterior-Power Representations: Candidate order parameters (operators) are viewed as vectors in the exterior power spaces of the Majorana vector space (e.g., bilinear operators correspond to the 2nd exterior power).
2. Intertwiner Decomposition: The authors use an algorithmic approach to decompose these induced representations into irreducible representations (irreps) of the Lie algebra $\mathfrak{g}$.
   • They compute the space of intertwiners (linear maps commuting with the algebra action).
   • If the intertwiner space is non-trivial, the representation is reducible. By diagonalizing a random linear combination of intertwiners, they extract invariant subspaces.
   • This process is applied recursively until all subspaces are irreducible.
3. Discrete Symmetry Integration: The continuous symmetry classification is refined by incorporating discrete lattice symmetries (e.g., sublattice exchange, layer exchange, time-reversal).
   • The algorithm checks if distinct invariant subspaces of the continuous group are isomorphic under the full symmetry group.
   • If isomorphic, they are recombined; if not, they form distinct irreps of the full group.

3. Key Contributions and Results

A. Benchmark: Honeycomb Hubbard Model

• Symmetry Recovery: The framework successfully recovers the well-known $SO(4)$ symmetry with algebra $\mathfrak{su}(2) \oplus \mathfrak{su}(2)$ (total spin and $\eta$-pairing pseudospin).
• Order Parameter Classification: The method exhaustively identifies 7 distinct candidate bilinear order parameters.
  • For repulsive interaction ($U>0$), it correctly identifies the antiferromagnetic order parameter.
  • For attractive interaction ($U<0$), it maps the antiferromagnetic order to superconducting and charge-density-wave orders, consistent with known physics.

B. Novel Application: Bilayer Spin-1/2 Model

The authors apply the framework to a bilayer honeycomb model with Heisenberg exchange and density-density interlayer couplings.

• New Symmetry Discovery: The model is found to possess a $Spin(5) \times U(1)/Z_2$ symmetry with algebra $\mathfrak{so}(5) \oplus \mathfrak{u}(1)$. This is a non-trivial emergent symmetry not immediately obvious in the complex fermion basis.
• Comprehensive Classification: The framework systematically classifies 18 independent local candidate bilinear order parameters. These are categorized by the specific symmetries they break (continuous $Spin(5) \times U(1)$, layer exchange $Z_2^l$, and sublattice exchange $Z_2^s$).
• Physical Implications: This classification reveals a rich landscape of potentially competing quantum phases. In a companion paper, Quantum Monte Carlo simulations on this model revealed an intermediate excitonic phase between the Dirac semimetal and the Symmetric Mass Generation (SMG) phase. The order parameter for this excitonic phase was directly identified from the classification provided in this work.

C. Extension to Pure Heisenberg Coupling

The framework was also applied to a related bilayer model with pure Heisenberg interlayer coupling, identifying an $SU(2)^3/Z_2$ symmetry ($\mathfrak{su}(2) \oplus \mathfrak{su}(2) \oplus \mathfrak{su}(2)$). This classification was crucial for proving the existence of SMG in that model by demonstrating the absence of long-range order in all 19 symmetry-breaking channels.

4. Significance

• Algorithmic Rigor: The paper replaces heuristic "guesswork" with a controlled, automated algebraic procedure. This is essential for systems with high-dimensional internal spaces where intuition often fails.
• Hidden Symmetry Detection: The Majorana-based approach is uniquely suited to uncovering emergent symmetries that are obscured in standard bases.
• Foundation for SMG: The systematic enumeration of order parameters is a prerequisite for establishing Symmetric Mass Generation (SMG), a mechanism for generating mass gaps without symmetry breaking. By proving that no order parameter condenses, one can confirm the existence of a gapped, symmetric phase.
• Generalizability: The method is readily extendable to multi-orbital, valley, and moiré systems (e.g., twisted bilayer graphene) and can be generalized to higher-order operators (composite and multipolar orders).

In summary, this work provides a powerful, general-purpose toolkit for the theoretical analysis of strongly correlated fermion systems, enabling the precise mapping of symmetry structures and the exhaustive identification of all possible phases of matter.