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

Resumen Técnico: Modelado de osciladores con retardo temporal con simetría $\mathcal{PT}$ más allá del límite de acoplamiento débil mediante una formulación de matriz de dispersión

Planteamiento del problema
La simetría paridad-tiempo ($\mathcal{PT}$) en osciladores con retardo temporal (TDO, por sus siglas en inglés), como láseres y osciladores optoelectrónicos (OEO), ofrece una vía para mejorar la pureza espectral, incluyendo la reducción del ruido de fase y la mejora de la supresión de modos secundarios. Sin embargo, las descripciones teóricas existentes se basan en formulaciones de modos acoplados derivadas bajo supuestos de envolvente de variación lenta y casi degeneración. Estos enfoques perturbativos restringen la validez al acoplamiento débil, al contraste de ganancia/pérdida pequeño y al desintonizamiento pequeño. En sistemas de gran retardo (por ejemplo, aquellos que emplean bucles de fibra de escala kilométrica), el retardo de propagación es comparable o superior a las escalas temporales características del sistema, lo que invalida la suposición de campos internos espacialmente uniformes. En consecuencia, existe la necesidad de una formulación no perturbativa que trate el retardo de propagación explícitamente sin depender de la truncación modal.

Metodología
Los autores desarrollan una formulación no perturbativa basada en una ecuación de diferencia de retardo y una representación de matriz de dispersión de la red de acoplamiento. El modelo considera un par de bucles de retardo idénticos con retardo de ida y vuelta $\tau_D$, interconectados por un acoplador lineal invariante en el tiempo caracterizado por una matriz de transmisión $\mathbf{K}$.

Pasos metodológicos clave:

1. Formulación: El sistema se describe mediante la ecuación $\mathbf{u} = \kappa \mathbf{K} \mathcal{D}_{\tau_D} \mathbf{u}$, donde $\mathbf{u}$ es el vector de estado, $\kappa$ es la ganancia/fase de ida y vuelta, y $\mathcal{D}_{\tau_D}$ es el operador de retardo.
2. Condiciones de simetría: Las condiciones para la invariancia $\mathcal{PT}$ se expresan directamente en términos de la matriz de transmisión del acoplador. La matriz $\mathbf{K}$ debe satisfacer $\mathbf{K} = \mathbf{J}(\mathbf{K}^{-1})^\dagger \mathbf{J}$ (donde $\mathbf{J}$ es la matriz de intercambio), asegurando que la red preserve la métrica pseudo-unitaria.
3. Descomposición de matrices: La matriz de transmisión se descompone en un producto de una matriz de acoplador sin pérdidas $\mathbf{U}(\theta)$ y una matriz de ganancia/pérdida balanceada $\mathbf{\Gamma}(\gamma)$. Esto permite que el sistema sea parametrizado por una fase de acoplamiento $\theta$ y un contraste de ganancia/pérdida $\gamma$.
4. Análisis de autovalores: Se deriva la estructura exacta de los autovalores resolviendo el polinomio característico de la matriz de transmisión, evitando las aproximaciones inherentes a la teoría de modos acoplados.

Contribuciones clave y resultados

• Estructura exacta de autovalores: Los autores derivan los autovalores exactos del sistema, mostrando que dependen de un parámetro de orden adimensional $\xi = \cosh(\gamma)\cos(\theta)$. Este parámetro gobierna la transición entre regímenes de simetría a través de todo el espacio de parámetros.
• Caracterización de regímenes:
  • Regimen no roto ($\xi^2 < 1$): Los autovalores forman un par conjugado complejo sobre el círculo unitario, lo que corresponde a dobletes de modo neutro con división de frecuencia $\nu$.
  • Regimen roto ($\xi^2 > 1$): Los autovalores forman un par real recíproco, lo que corresponde a modos con amplificación y atenuación recíprocas (división de ganancia-pérdida $\mu$).
  • Puntos excepcionales ($\xi^2 = 1$): Estos puntos delinean la transición de fase donde la matriz se vuelve no diagonalizable.
• Transición de simetría global: El análisis demuestra que la transición de simetría es una propiedad global del sistema acoplado, determinada únicamente por los parámetros de acoplamiento y ganancia/pérdida, y es independiente de la ganancia total del bucle o de la amplificación dependiente de la frecuencia.
• Recuperación de la teoría de modos acoplados: El artículo demuestra que la teoría convencional de modos acoplados emerge como el límite asintótico de esta formulación exacta para parámetros pequeños (pequeños $\gamma$, $\theta$ y $\ln(\kappa)$). Esto establece una conexión matemática directa entre las descripciones perturbativas y no perturbativas.
• Visualizaciones: El estudio proporciona curvas de nivel para el parámetro de orden, la división de frecuencia y la división de ganancia-pérdida, contrastando los resultados exactos con las curvas punteadas de la aproximación convencional de modos acoplados, resaltando la divergencia en los regímenes de acoplamiento fuerte.

Significancia
El artículo afirma proporcionar un marco unificado y físicamente transparente para analizar sistemas de retardo con simetría $\mathcal{PT}$ más allá del límite de acoplamiento débil. Al tratar explícitamente el retardo de propagación y evitar la truncación modal, la formulación sigue siendo válida para cualquier fuerza de acoplamiento, desequilibrio de ganancia/pérdida y desintonización de resonancia.

La significancia radica en:

1. Optimización del diseño: El marco ofrece una base rigurosa para el diseño y la optimización de osciladores de bajo ruido y sistemas fotónicos, particularmente aquellos que operan en regímenes de acoplamiento fuerte o de gran retardo, inaccesibles para los modelos perturbativos.
2. Claridad teórica: Clarifica el rango de validez de los modelos de modos acoplados ampliamente utilizados, mostrando que son límites asintóticos de una teoría exacta más general.
3. Distinción de efectos: La identificación de una transición de simetría global independiente de la ganancia del bucle proporciona una base para distinguir los efectos intrínsecos de la simetría $\mathcal{PT}$ de los comportamientos selectivos de frecuencia derivados de filtrado externo o dispersión de ganancia.

Los autores señalan que, si bien el presente trabajo se centra en el modelo de diferencia de retardo lineal, extensiones futuras podrían incluir dispersión, propagación bidireccional y mecanismos de saturación no lineal para aplicar aún más estos hallazgos a la predicción del ruido de fase y la selección de modos en arquitecturas avanzadas de osciladores.