PT-symmetric time delay oscillator modelling beyond the weak coupling limit via a scattering matrix formulation

Technical Summary: 𝒫𝒯-symmetric time delay oscillator modelling beyond the weak coupling limit via a scattering matrix formulation

Problem Statement
Parity–time (𝒫𝒯) symmetry in time-delay oscillators (TDOs), such as lasers and optoelectronic oscillators (OEOs), offers a pathway to enhanced spectral purity, including reduced phase noise and improved sidemode suppression. However, existing theoretical descriptions rely on coupled-mode formulations derived under slowly varying envelope and near-degeneracy assumptions. These perturbative approaches restrict validity to weak coupling, small gain/loss contrast, and small detuning. In large-delay systems (e.g., those employing kilometer-scale fiber loops), the propagation delay is comparable to or exceeds the system's characteristic timescales, rendering the assumption of spatially uniform internal fields invalid. Consequently, there is a need for a non-perturbative formulation that treats propagation delay explicitly without relying on modal truncation.

Methodology
The authors develop a non-perturbative formulation based on a delay–difference equation and a scattering matrix representation of the coupling network. The model considers a pair of identical delay loops with round-trip delay $\tau_D$, interconnected by a linear time-invariant coupler characterized by a transmission matrix $\mathbf{K}$.

Key methodological steps include:

1. Formulation: The system is described by the equation $\mathbf{u} = \kappa \mathbf{K} \mathcal{D}_{\tau_D} \mathbf{u}$, where $\mathbf{u}$ is the state vector, $\kappa$ is the round-trip gain/phase, and $\mathcal{D}_{\tau_D}$ is the delay operator.
2. Symmetry Conditions: The conditions for 𝒫𝒯-invariance are expressed directly in terms of the coupler's transmission matrix. The matrix $\mathbf{K}$ must satisfy $\mathbf{K} = \mathbf{J}(\mathbf{K}^{-1})^\dagger \mathbf{J}$ (where $\mathbf{J}$ is the exchange matrix), ensuring the network preserves the pseudo-unitary metric.
3. Matrix Decomposition: The transmission matrix is decomposed into a product of a lossless coupler matrix $\mathbf{U}(\theta)$ and a balanced gain/loss matrix $\mathbf{\Gamma}(\gamma)$. This allows the system to be parameterized by a coupling phase $\theta$ and a gain/loss contrast $\gamma$.
4. Eigenvalue Analysis: The exact eigenvalue structure is derived in closed form by solving the characteristic polynomial of the transmission matrix, avoiding the approximations inherent in coupled-mode theory.

Key Contributions and Results

• Exact Eigenvalue Structure: The authors derive the exact eigenvalues of the system, showing they depend on a dimensionless order parameter $\xi = \cosh(\gamma)\cos(\theta)$. This parameter governs the transition between symmetry regimes across the full parameter space.
• Regime Characterization:
  • Unbroken Regime ($\xi^2 < 1$): The eigenvalues form a complex conjugate pair on the unit circle, corresponding to neutral mode doublets with frequency splitting $\nu$.
  • Broken Regime ($\xi^2 > 1$): The eigenvalues form a reciprocal real pair, corresponding to modes with reciprocal amplification and attenuation (gain-loss splitting $\mu$).
  • Exceptional Points ($\xi^2 = 1$): These points delineate the phase transition where the matrix becomes non-diagonalizable.
• Global Symmetry Transition: The analysis demonstrates that the symmetry transition is a global property of the coupled system, determined solely by coupling and gain/loss parameters, and is independent of the overall loop gain or frequency-dependent amplification.
• Recovery of Coupled-Mode Theory: The paper proves that conventional coupled-mode theory emerges as the asymptotic limit of this exact formulation for small parameters (small $\gamma$, $\theta$, and $\ln(\kappa)$). This establishes a direct mathematical connection between perturbative and non-perturbative descriptions.
• Visualizations: The study provides level curves for the order parameter, frequency splitting, and gain-loss splitting, contrasting the exact results with the dotted curves of the conventional coupled-mode approximation, highlighting the divergence in strong coupling regimes.

Significance
The paper claims to provide a unified and physically transparent framework for analyzing 𝒫𝒯-symmetric delay systems beyond the weak-coupling limit. By explicitly treating propagation delay and avoiding modal truncation, the formulation remains valid for arbitrary coupling strength, gain/loss imbalance, and resonance detuning.

The significance lies in:

1. Design Optimization: The framework offers a rigorous basis for the design and optimization of low-noise oscillators and photonic systems, particularly those operating in strong coupling or large-delay regimes inaccessible to perturbative models.
2. Theoretical Clarity: It clarifies the range of validity of widely used coupled-mode models by showing them as asymptotic limits of a more general exact theory.
3. Distinguishing Effects: The identification of a global symmetry transition independent of loop gain provides a basis for distinguishing intrinsic 𝒫𝒯-symmetry effects from frequency-selective behaviors arising from external filtering or gain dispersion.

The authors note that while the current work focuses on the linear delay-difference model, future extensions could include dispersion, bidirectional propagation, and nonlinear saturation mechanisms to further apply these findings to the prediction of phase noise and mode selection in advanced oscillator architectures.