🔬 optics

Intensity-Based Criterion for Determining Exceptional Point in Parity-Time (PT) Symmetric Coupled Array of Optical Waveguides

This paper proposes a novel, efficient intensity-based criterion for determining the exceptional point in PT-symmetric optical waveguide arrays, which avoids the computational complexity of Hamiltonian diagonalization while offering results consistent with eigenvalue analysis and general applicability to tight-binding systems.

Original authors: Mahla Bahar, Mojtaba Golshani, Mostafa Motamedifar, Khatereh Jafari

Published 2026-03-31

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Mahla Bahar, Mojtaba Golshani, Mostafa Motamedifar, Khatereh Jafari

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Big Picture: A Light Show with a "Tipping Point"

Imagine you have a long row of 50 connected water pipes (these are the optical waveguides). You can pour water (light) into one of them. In a normal system, the water would flow from pipe to pipe, spreading out evenly like a ripple in a pond.

But in this experiment, the scientists set up a special "Parity-Time" (PT) symmetric system. Think of this as a row of pipes where some are leaking (losing water) and others are being refilled (gaining water) at the exact same rate.

The Goal: They wanted to find the exact moment when this delicate balance breaks. They call this moment the Exceptional Point (EP).
The Problem: Usually, to find this breaking point, you have to do incredibly difficult math (solving a giant puzzle called "diagonalizing the Hamiltonian"). It's like trying to predict exactly how a complex machine will break by calculating the stress on every single screw inside it. It takes a supercomputer and a lot of time.

The New Idea: A Simple "Thermometer"

The authors of this paper said, "Wait a minute. We don't need to solve the whole puzzle. We just need to watch the water flow."

They developed a new, simple rule (an intensity-based criterion) to spot the breaking point just by looking at how bright the light is in the pipes over time.

Here is how their "thermometer" works:

Before the Breaking Point (The Calm Phase):
Imagine the water is sloshing back and forth between the pipes. It goes up, it goes down, it spreads out. The amount of water in the first half of the journey is roughly the same as the amount in the second half.
- The Math: The system is stable. The "score" (called $I_r$ ) stays at 1.
After the Breaking Point (The Chaos Phase):
Suddenly, the balance tips. The pipes that were "refilling" start to explode with water, while the "leaking" pipes go dry. The light doesn't just spread; it explodes exponentially. The second half of the journey has way more water than the first half.
- The Math: The system is unstable. The "score" ( $I_r$ ) drops to 0.

The Discovery: The exact moment the score drops from 1 to 0 is the Exceptional Point.

Why This Matters: The "Magic" of the Tipping Point

Why do we care about this tipping point? Because systems right at the edge of breaking are super-sensitive.

Analogy: Think of a pencil balanced perfectly on its tip. It's in a "critical" state. If you blow the tiniest breath of air on it, it falls.
Real World Use: If you build a sensor using this "tipping point" light system, it can detect the tiniest changes in the environment (like a tiny virus or a slight change in temperature) that normal sensors would miss. The system reacts massively to tiny nudges.

The Quantum Twist: The "Ghost" Particles

The paper also looked at what happens when you send in quantum light (photons) instead of just a steady beam. They used different "types" of light:

Number States: Like dropping exactly 2 marbles into the pipe.
Coherent States: Like a steady stream of rain.
Entangled States: Like two magical marbles that are linked; if you find one here, the other is definitely there, even if they are far apart.

The Surprise: They found that no matter what kind of light you send in (marbles, rain, or magic links), the tipping point (the EP) stays exactly the same. It's a property of the pipes themselves, not the water you pour in.

They also looked at how the light particles "correlated" (how likely they were to be found together).

Before the break: The particles spread out nicely, sometimes staying close, sometimes moving apart.
After the break: The particles get "scared" and huddle together right next to the input pipe, because that's where the "gain" (the refilling) is strongest.

The "Cheat Code" for Scientists

The biggest takeaway from this paper is a computational shortcut.

The Old Way: To find the tipping point, you have to calculate the energy levels of the system. If you have many waveguides, the math becomes impossible (the "Hilbert space" grows exponentially, like trying to count every possible arrangement of a deck of cards as you add more cards).
The New Way: Just watch the light intensity. If it starts growing wildly in the second half of the journey, you've passed the tipping point.

The Analogy:

Old Way: Trying to predict a traffic jam by calculating the speed, weight, and fuel efficiency of every single car on the highway.
New Way: Just looking out the window. If you see a massive pile-up forming, you know traffic has broken down. You don't need to know the details of every car to know the jam is there.

Summary

This paper introduces a simple, visual way to find the "breaking point" of a special light system. Instead of doing impossible math, scientists can just watch how the light intensity changes. This makes it much easier to design super-sensitive sensors and understand how quantum light behaves in complex networks.

1. Problem Statement

Parity-Time (PT) symmetric optical systems, characterized by balanced gain and loss, exhibit a phase transition at an Exceptional Point (EP). Below the EP, the system has real eigenvalues (unbroken PT symmetry); above it, eigenvalues become complex, leading to exponential growth or decay (broken PT symmetry).

The Challenge: Identifying the EP is crucial for applications like enhanced sensing and zero group velocity. Traditionally, the EP is found by diagonalizing the system's Hamiltonian to find where eigenvalues transition from real to complex.
The Bottleneck: In quantum optical systems (waveguide arrays), the Hamiltonian acts on an infinite-dimensional Hilbert space. Even when truncated to a specific photon number, the dimension grows exponentially with the number of waveguides ( $N$ ) and photons ( $m$ ). This makes direct Hamiltonian diagonalization computationally prohibitive for large arrays or complex quantum input states (e.g., coherent or squeezed states).
The Goal: The authors aim to develop a computationally efficient, intensity-based criterion to determine the EP without requiring full Hamiltonian diagonalization, applicable to various quantum input states.

2. Methodology

The study investigates a 1D array of $N$ coupled optical waveguides with alternating gain and loss (PT-symmetric). The system is modeled using a non-Hermitian Hamiltonian within the tight-binding approximation.

A. Theoretical Framework

Hamiltonian: Described using creation ( $\hat{a}^\dagger$ ) and annihilation ( $\hat{a}$ ) operators. The propagation constant includes an imaginary term ( $\pm ig$ ) representing gain and loss.
Dynamics: The evolution of field operators is governed by the Heisenberg equation of motion. The solution is expressed via a Green's function ( $G_{jm}(z)$ ), which relates the field at waveguide $j$ and distance $z$ to the initial input at waveguide $m$ .
Input States: The study simulates various quantum inputs:
- Separable States: Number states ( $|n\rangle$ ), Coherent states ( $|\alpha\rangle$ ), and Squeezed states ( $|\xi\rangle$ ).
- Entangled States: NOON states ( $|n\rangle_p|0\rangle_q + |0\rangle_p|n\rangle_q$ ) injected into specific waveguide pairs.
Observables:
1. Intensity ( $I_j$ ): The expectation value of the photon number, representing the probability of detecting a photon.
2. Site-to-Site Correlation ( $g^{(2)}_{ji}$ ): The second-order correlation function, representing the joint probability of detecting two photons at waveguides $i$ and $j$ .

B. Two Approaches for EP Determination

Standard Method (Eigenvalue Analysis):
- The Hamiltonian is constructed in a truncated $m$ -photon subspace.
- Eigenvalues are computed numerically.
- The EP ( $g_c$ ) is identified as the gain/loss strength where the maximum imaginary part of the eigenvalues becomes non-zero.
Proposed Method (Intensity-Based Criterion):
- Instead of diagonalization, the authors solve the coupled linear differential equations for the Green's function to simulate light propagation.
- They define a metric, $I_r$ , based on the spatial intensity distribution:
  $I_r = \exp\left(-\gamma \left| \frac{\bar{I}_2 - \bar{I}_1}{\bar{I}_2 + \bar{I}_1} \right|^\beta \right)$
  Where $\bar{I}_1$ and $\bar{I}_2$ are the mean intensities in the first and second halves of the propagation length, respectively.
- Logic:
  - Unbroken Phase ( $g < g_c$ ): Energy oscillates; $\bar{I}_1 \approx \bar{I}_2$ , so $I_r \approx 1$ .
  - Broken Phase ( $g > g_c$ ): Energy grows exponentially; $\bar{I}_2 \gg \bar{I}_1$ , so the ratio approaches 1, and $I_r \approx 0$ .
- The EP is identified as the transition point where $I_r$ drops from 1 to 0.

3. Key Results

Propagation Patterns and Correlations

Hermitian/Unbroken PT Regime ( $g < g_c$ ):
- Propagation patterns exhibit discrete diffraction.
- Total intensity oscillates but remains bounded.
- Site-to-site correlations show distinct spatial features (e.g., spatial anti-bunching for central inputs, bunching for boundary inputs) but remain stable and similar to Hermitian systems.
Broken PT Regime ( $g > g_c$ ):
- Exponential Growth: Intensity grows exponentially, dominated by waveguides with gain near the input.
- Correlation Shift: The site-to-site correlation function concentrates heavily around the input waveguide and nearby gain regions, indicating that photons are most likely to be detected near the source due to gain amplification.

Validation of the Intensity-Based Criterion

Quantitative Agreement: The EP ( $g_c$ $g_{c}$ ) calculated using the $I_r$ $I_{r}$ criterion was compared against the eigenvalue method.
- Example: For $N=10$ waveguides, the eigenvalue method yielded $g_c = 0.284$ , while the $I_r$ method yielded $g_c \approx 0.289$ . The difference is negligible ( $\sim 0.005$ ).
Input Independence: The EP value was found to be invariant regardless of the input quantum state (number, coherent, squeezed, or entangled). This confirms that $g_c$ is an intrinsic property of the system's Hamiltonian, not the initial state.
Photon Number Independence: The onset of complex eigenvalues (the EP) does not depend on the number of photons ( $m$ ) in the subspace, validating the use of single-photon subspaces to determine the EP for the full system.

Dependence on System Size

Waveguide Count ( $N$ ): As $N$ increases, the critical gain/loss strength ( $g_c$ ) decreases. This means larger arrays enter the broken PT-symmetry phase at lower gain/loss thresholds.
Parity Effect: Arrays with an odd number of waveguides (featuring a central neutral waveguide) exhibit a higher $g_c$ compared to even arrays (which alternate strictly between gain and loss).

4. Key Contributions

Novel Criterion: Introduction of a straightforward, intensity-based metric ( $I_r$ ) to detect the PT-symmetry breaking point without diagonalizing the Hamiltonian.
Computational Efficiency: The proposed method avoids the exponential computational cost associated with diagonalizing high-dimensional Hilbert spaces, making it scalable for large waveguide arrays and complex quantum states (like coherent/squeezed states where photon number is undefined).
Universality: Demonstrated that the criterion works across a wide class of input states (separable and entangled) and is applicable to any system describable by a tight-binding framework.
Physical Insight: Provided a detailed analysis of how quantum correlations (bunching/anti-bunching) evolve across the PT phase transition, showing a distinct shift from delocalized interference patterns to localized gain-dominated behavior in the broken phase.

5. Significance

Scalability: This approach solves a major bottleneck in simulating non-Hermitian quantum systems. It allows researchers to study PT-symmetry breaking in large-scale arrays where traditional eigenvalue methods are computationally infeasible.
Experimental Relevance: Since the criterion relies on intensity measurements (which are directly observable in experiments), it offers a practical, real-time method for identifying the EP in experimental setups without needing full state tomography or complex theoretical reconstruction.
Sensor Design: By accurately and efficiently locating the EP, this method aids in the design of ultra-sensitive optical sensors that operate near the exceptional point, where system response to perturbations is maximized.
General Applicability: The framework extends beyond optics to other Schrödinger-type dynamics, including spin chains and population transfer in multi-level systems.