Puiseux series about exceptional singularities dictated by symmetry-allowed Hessenberg forms of perturbation matrices

This paper establishes a systematic framework linking the Hessenberg structure of symmetry-preserving perturbations to the order of Puiseux series singularities at exceptional points in non-Hermitian systems, demonstrating that while parity and charge-conjugation symmetries restrict third-order exceptional points to square-root branch points, parity-time symmetry generically allows for the strongest cubic-root singularities.

The Gist

Imagine you are standing in a vast, foggy landscape made of invisible hills and valleys. In the world of physics, this landscape represents the energy levels of a material. Usually, these energy levels are like smooth, rolling hills. But in a special kind of material called a Non-Hermitian system (think of it as a system that loses or gains energy, like a guitar string that keeps vibrating forever or stops abruptly), the landscape can get very strange.

Sometimes, two or more hills crash into each other at a single, sharp point. Physicists call these crash points Exceptional Points (EPs). At these points, the rules of normal physics break down: the energy levels merge, and the "directions" the particles can move (called eigenvectors) also merge into one.

This paper is like a map and a rulebook for understanding exactly how these crashes happen and what happens when you nudge the system slightly away from the crash.

Here is the breakdown of the paper's big ideas using simple analogies:

1. The "Hessenberg" Shape: The Traffic Light of Chaos

The authors discovered that the "shape" of the rules governing these systems acts like a traffic light for how the energy splits when you nudge it.

• The Analogy: Imagine you have a stack of cards representing energy levels. If you push the stack (add a small perturbation), how do the cards separate?
• The Rule: The paper says the answer depends on a specific geometric pattern in the math, called a Hessenberg form. Think of this as a "staircase" of allowed connections.
  • If the staircase has 2 steps, the energy splits like a square root (like $\sqrt{\epsilon}$). This is a "gentle" split.
  • If the staircase has 3 steps, the energy splits like a cube root (like $\sqrt[3]{\epsilon}$). This is a "sharper," more dramatic split.
  • If it has 4 steps, it splits like a fourth root ($\sqrt[4]{\epsilon}$).

The authors call this a Puiseux series, which is just a fancy math way of saying: "The answer isn't a whole number; it's a fraction."

2. The Three Symmetry Guards

The paper looks at three different "guards" (symmetries) that protect these systems: Parity (P), Charge (C), and Parity-Time (PT). These guards decide how many "steps" the staircase can have.

• The P and C Guards (The Strict Parents):
  Imagine P and C are strict parents who say, "You can only have a 2-step staircase." Even if you try to make a 3-step crash (an EP3), these guards force it to behave like a 2-step crash.

  • Result: The energy splits slowly (square root). It's like a door opening slowly.
  • The Catch: In these systems, there is always one "flat" energy level (a flat band) that doesn't move. It's like a spectator sitting on the sidelines while the other two crash. Because of this spectator, the crash can never be as wild as it could be.

• The PT Guard (The Wild Card):
  The PT guard is more flexible. It allows for a 3-step staircase.

  • Result: The energy splits much faster and more dramatically (cube root). It's like a door being blasted open.
  • Significance: This is the "strongest" possible crash for a 3-band system. The paper shows that if you want the most extreme, sensitive reaction from a material, you need the PT symmetry.

3. Tuning the Sensitivity: The "Knob" Trick

One of the coolest findings is that you can fine-tune these crashes.

• The Analogy: Imagine a radio dial. Usually, turning the dial changes the volume smoothly. But at an Exceptional Point, turning the dial changes the volume wildly (a tiny turn makes a huge noise).
• The Discovery: The authors show that by adjusting specific parameters (turning specific "knobs" in the math), you can change the nature of the crash.
  • You can make a "wild" crash (cube root) behave like a "gentle" one (square root) or even a "normal" one (linear).
  • Why does this matter? This is crucial for building sensors. Imagine a sensor that detects a tiny change in the air. If you tune it to the "wild" crash, it will scream when a tiny change happens. If you tune it to the "gentle" crash, it will whisper. You can even make it so the sensor is super sensitive to wind coming from the North, but ignores wind from the East. This is called direction-dependent sensing.

4. Real-World Examples: The "Hopf" and "Pseudospin" Models

The paper doesn't just stay in theory; it builds 3D models (like virtual Lego structures) to prove this works.

• They created models where these crashes form lines (like a wire) or surfaces (like a sheet) in space.
• They showed that in some models, these lines can get knotted (like a pretzel), which is a very topological and "cool" feature.

The Big Picture Takeaway

This paper is a user manual for chaos.

Before this, scientists knew that these "crash points" (Exceptional Points) existed and were weird. This paper tells us exactly how weird they will be based on the symmetry of the material.

• If you want a mild reaction: Use P or C symmetry.
• If you want the most extreme, sensitive reaction: Use PT symmetry.
• If you want to control the direction of the sensitivity: You can "fine-tune" the knobs to switch between these behaviors.

This is a huge step forward for designing super-sensitive sensors (for medical devices, environmental monitoring, or quantum computers) that can be programmed to react differently depending on where a signal comes from. It turns the abstract math of "fractional powers" into a practical toolkit for engineering the future.

Technical Summary

1. Problem Statement

Non-Hermitian (NH) systems exhibit Exceptional Points (EPs), which are spectral singularities where eigenvalues and eigenvectors coalesce. While second-order EPs (EP2s) are well-understood in 2D systems, the nature of higher-order EPs (EP$_n$ with $n \geq 3$) in higher-dimensional systems remains complex.

• The Challenge: Determining the precise algebraic nature of the singularity (i.e., the fractional power of the perturbation parameter $\epsilon$) around an EP$_n$.
• The Gap: General counting arguments suggest EP$_n$ requires tuning $2n-2$ parameters, implying they should only exist in high-dimensional spaces. However, symmetries can stabilize them in lower dimensions. The specific relationship between the symmetry class of a system, the structure of allowed perturbations, and the resulting Puiseux series expansion (the splitting behavior $\propto \epsilon^{1/k}$) has not been systematically categorized.

2. Methodology

The author develops a systematic framework based on linear algebra and perturbation theory:

• Jordan-Normal Basis Transformation: The analysis begins by transforming the defective Hamiltonian at an EP into its Jordan-normal form ($J_s(\lambda)$).
• Hessenberg Matrix Structure: Perturbations are analyzed in this basis. The paper establishes that symmetry-preserving perturbations manifest as upper-$k$ Hessenberg matrices ($B_{s,k}$).
  • An upper-$k$ Hessenberg matrix has non-zero elements only on the main diagonal and the first $k-1$ super-diagonals (or sub-diagonals depending on convention, here defined by the non-zero lower diagonals in the perturbation relative to the Jordan block).
• Puiseux Series Expansion: The eigenvalue splitting is derived using the Puiseux series, a generalized power series allowing fractional exponents.
  • Key Theorem: If a perturbation in the Jordan basis forms an upper-$k$ Hessenberg matrix, the leading-order eigenvalue splitting scales as $\epsilon^{1/k}$.
• Symmetry Constraints: The framework is applied to 3-band (3$\times$3) and 4-band (4$\times$4) models constrained by three specific symmetries:
  1. Parity (P)
  2. Charge-Conjugation (C)
  3. Parity-Time Reversal (PT)

3. Key Contributions

1. General Framework: A rigorous mathematical link between the Hessenberg index ($k$) of a symmetry-allowed perturbation and the singularity order ($1/k$) of the EP.
2. Symmetry Dictation: Demonstrating that specific symmetries restrict the allowed perturbation matrices to specific Hessenberg forms, thereby dictating the maximum possible singularity strength for a given system.
3. Fine-Tuning Mechanism: Showing that specific parameter choices can suppress the leading-order singular term (e.g., eliminating the $\epsilon^{1/2}$ term), resulting in a less singular behavior (e.g., $\epsilon^1$). This implies direction-dependent sensitivity.

4. Detailed Results

A. Two-Band Models (Review)

• In 2$\times$2 systems, perturbations are inherently upper-2 Hessenberg.
• Result: EP2s always exhibit square-root singularities ($\epsilon^{1/2}$).

B. Three-Band Models (3$\times$3)

The paper analyzes three symmetry classes:

1. P-Symmetric Systems:

   • Constraint: The symmetry forces one eigenvalue to be zero (flat band) and the other two to be negatives of each other.
   • Perturbation Structure: Allowed perturbations are restricted to upper-2 Hessenberg forms ($B_{3,2}$).
   • Singularity: The leading-order splitting is $\epsilon^{1/2}$ (square-root).
   • Note: An upper-3 Hessenberg perturbation would break P-symmetry.
   • Physical Example: 3D Hopf semimetals with BC-dipoles exhibit Exceptional Curves (EC3) with square-root singularities.

2. C-Symmetric Systems:

   • Constraint: Similar to P-symmetry, enforces a zero eigenvalue and symmetric pairs.
   • Perturbation Structure: Restricted to upper-2 Hessenberg forms ($B_{3,2}$).
   • Singularity: Leading-order splitting is $\epsilon^{1/2}$.
   • Physical Example: Pseudospin-1 nodes in 3D exhibit EC3s with square-root singularities.

3. PT-Symmetric Systems:

   • Constraint: Does not enforce a flat band or symmetric pairs in the same way; allows for a full 3-fold degeneracy.
   • Perturbation Structure: Allows for upper-3 Hessenberg forms ($B_{3,3}$).
   • Singularity: Supports the strongest possible singularity for 3-band systems: $\epsilon^{1/3}$ (cube-root).
   • Physical Example: A 3D model with $H_{PT}$ exhibits knotted Exceptional Curves (EC3) and surfaces (ES2) where the splitting follows the cube-root law.

C. Four-Band Models (Appendix)

• P- and C-Symmetric 4$\times$4 Matrices:
  • The analysis extends to 4-band models.
  • Result: Symmetry constraints allow for upper-4 Hessenberg perturbations.
  • Singularity: These systems host EP4s with $\epsilon^{1/4}$ singularities.

5. Significance and Implications

• Design of Sensors: The most significant practical implication is for EP-based sensors. The sensitivity of an EP sensor scales with the order of the singularity ($\epsilon^{1/k}$).
  • PT-symmetric systems offer the highest sensitivity ($\epsilon^{1/3}$) compared to P/C-symmetric systems ($\epsilon^{1/2}$).
  • Directional Sensitivity: By fine-tuning perturbations to suppress the leading singular term, one can create sensors where the response is highly sensitive in specific directions (strong singularity) but linear or weak in others. This enables multi-functional, direction-dependent sensing.
• Topological Classification: The work clarifies the topological landscape of NH systems, showing how symmetries stabilize higher-order EPs (ECs and ESs) in lower dimensions than generic counting arguments would suggest.
• Future Directions: The paper suggests exploring the connection between the degree of EP singularity and the Non-Hermitian Skin Effect (NHSE), as well as the role of these singularities in topological pumping and charge transport.

Conclusion

The paper provides a unifying mathematical language (Hessenberg forms and Puiseux series) to predict the behavior of exceptional singularities. It establishes that symmetry is the primary dictator of singularity strength: P and C symmetries limit 3-band systems to square-root singularities, while PT symmetry unlocks the stronger cube-root singularities. This insight is crucial for engineering next-generation non-Hermitian devices with tailored sensitivity profiles.