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

Résumé technique : Modélisation d'oscillateurs à retard avec symétrie $\mathcal{PT}$ au-delà de la limite de couplage faible via une formulation par matrice de diffusion

Énoncé du problème
La symétrie parité-temps ($\mathcal{PT}$) dans les oscillateurs à retard temporel (TDO), tels que les lasers et les oscillateurs optoélectroniques (OEO), offre une voie vers une pureté spectrale accrue, incluant une réduction du bruit de phase et une amélioration de la suppression des modes secondaires. Cependant, les descriptions théoriques existantes reposent sur des formulations de modes couplés dérivées sous des hypothèses d'enveloppe à variation lente et de quasi-dégénérescence. Ces approches perturbatives restreignent la validité aux cas de couplage faible, de faible contraste gain/perte et de faible désaccord. Dans les systèmes à retard important (par exemple, employant des boucles de fibre à l'échelle du kilomètre), le retard de propagation est comparable ou supérieur aux échelles de temps caractéristiques du système, rendant invalide l'hypothèse de champs internes spatialement uniformes. Par conséquent, il est nécessaire d'une formulation non perturbative qui traite explicitement le retard de propagation sans recourir à la troncature modale.

Méthodologie
Les auteurs développent une formulation non perturbative basée sur une équation de différence à retard et une représentation par matrice de diffusion du réseau de couplage. Le modèle considère une paire de boucles à retard identiques avec un retard de tour $\tau_D$, interconnectées par un coupleur linéaire invariant dans le temps caractérisé par une matrice de transmission $\mathbf{K}$.

Étapes méthodologiques clés :

1. Formulation : Le système est décrit par l'équation $\mathbf{u} = \kappa \mathbf{K} \mathcal{D}_{\tau_D} \mathbf{u}$, où $\mathbf{u}$ est le vecteur d'état, $\kappa$ est le gain/phase par tour, et $\mathcal{D}_{\tau_D}$ est l'opérateur de retard.
2. Conditions de symétrie : Les conditions d'invariance $\mathcal{PT}$ sont exprimées directement en termes de la matrice de transmission du coupleur. La matrice $\mathbf{K}$ doit satisfaire $\mathbf{K} = \mathbf{J}(\mathbf{K}^{-1})^\dagger \mathbf{J}$ (où $\mathbf{J}$ est la matrice d'échange), garantissant que le réseau préserve la métrique pseudo-unitaire.
3. Décomposition matricielle : La matrice de transmission est décomposée en un produit d'une matrice de coupleur sans perte $\mathbf{U}(\theta)$ et d'une matrice de gain/perte équilibrée $\mathbf{\Gamma}(\gamma)$. Cela permet de paramétrer le système par une phase de couplage $\theta$ et un contraste gain/perte $\gamma$.
4. Analyse des valeurs propres : La structure exacte des valeurs propres est dérivée sous forme fermée en résolvant le polynôme caractéristique de la matrice de transmission, évitant ainsi les approximations inhérentes à la théorie des modes couplés.

Contributions clés et résultats

• Structure exacte des valeurs propres : Les auteurs dérivent les valeurs propres exactes du système, montrant qu'elles dépendent d'un paramètre d'ordre adimensionnel $\xi = \cosh(\gamma)\cos(\theta)$. Ce paramètre régit la transition entre les régimes de symétrie à travers tout l'espace des paramètres.
• Caractérisation des régimes :
  • Régime non brisé ($\xi^2 < 1$) : Les valeurs propres forment une paire conjuguée complexe sur le cercle unité, correspondant à des doublets de modes neutres avec un éclatement de fréquence $\nu$.
  • Régime brisé ($\xi^2 > 1$) : Les valeurs propres forment une paire réelle réciproque, correspondant à des modes avec amplification et atténuation réciproques (éclatement gain-perte $\mu$).
  • Points exceptionnels ($\xi^2 = 1$) : Ces points délimitent la transition de phase où la matrice devient non diagonalisable.
• Transition de symétrie globale : L'analyse démontre que la transition de symétrie est une propriété globale du système couplé, déterminée uniquement par les paramètres de couplage et de gain/perte, et qu'elle est indépendante du gain de boucle global ou de l'amplification dépendante de la fréquence.
• Récupération de la théorie des modes couplés : L'article prouve que la théorie conventionnelle des modes couplés émerge comme la limite asymptotique de cette formulation exacte pour de petits paramètres (petits $\gamma$, $\theta$, et $\ln(\kappa)$). Cela établit une connexion mathématique directe entre les descriptions perturbative et non perturbative.
• Visualisations : L'étude fournit des courbes de niveau pour le paramètre d'ordre, l'éclatement de fréquence et l'éclatement gain-perte, contrastant les résultats exacts avec les courbes en pointillés de l'approximation conventionnelle des modes couplés, mettant en évidence la divergence dans les régimes de couplage fort.

Signification
L'article affirme fournir un cadre unifié et physiquement transparent pour analyser les systèmes à retard $\mathcal{PT}$-symétriques au-delà de la limite de couplage faible. En traitant explicitement le retard de propagation et en évitant la troncature modale, la formulation reste valide pour un couplage, un déséquilibre gain/perte et un désaccord de résonance arbitraires.

La signification réside dans :

1. Optimisation de la conception : Le cadre offre une base rigoureuse pour la conception et l'optimisation d'oscillateurs à faible bruit et de systèmes photoniques, notamment ceux opérant dans des régimes de couplage fort ou de grand retard, inaccessibles aux modèles perturbatifs.
2. Clarté théorique : Il clarifie la portée de validité des modèles de modes couplés largement utilisés en montrant qu'ils sont des limites asymptotiques d'une théorie exacte plus générale.
3. Distinction des effets : L'identification d'une transition de symétrie globale indépendante du gain de boucle fournit une base pour distinguer les effets intrinsèques de la symétrie $\mathcal{PT}$ des comportements sélectifs en fréquence résultant de la dispersion du gain ou du filtrage externe.

Les auteurs notent que bien que le présent travail se concentre sur le modèle de différence de retard linéaire, des extensions futures pourraient inclure la dispersion, la propagation bidirectionnelle et les mécanismes de saturation non linéaires afin d'appliquer davantage ces résultats à la prédiction du bruit de phase et de la sélection de mode dans les architectures d'oscillateurs avancées.