Time-resolved synchronization analysis of stacked intrinsic Josephson junctions of a cuprate superconductor with frequency-modulated terahertz radiation spectra

Résumé technique : Analyse de la synchronisation résolue en temps des jonctions Josephson intrinsèques empilées d'un supraconducteur de type cuprate avec des spectres de rayonnement terahertz à modulation de fréquence

Énoncé du problème
Les jonctions Josephson intrinsèques (IJJ) dans les supraconducteurs de type cuprate, spécifiquement le Bi$_2$Sr$_2$CaCu$_2$O$_{8+\delta}$ (Bi2212), servent de systèmes quantiques macroscopiques pour l'étude d'oscillateurs non linéaires en interaction. Bien que ces systèmes puissent générer un rayonnement térahertz (THz) via l'effet Josephson alternatif (AC), une compréhension quantitative complète de leurs mécanismes de synchronisation dynamique — particulièrement sous modulation de fréquence (FM) — est restée élusive. Des travaux antérieurs par les auteurs ont démontré des spectres FM et ont attribué les phénomènes observés à un temps de relaxation de synchronisation fini, mais une estimation précise de ce temps ainsi qu'une discussion détaillée du mécanisme sous-jacent faisaient défaut. La complexité provient de la coexistence de deux interactions distinctes : le couplage inductif à travers des centaines de jonctions (caractérisé par la longueur de pénétration $\lambda_{ab}$) et le couplage capacitif entre voisins proches (caractérisé par la longueur de Debye $\mu$).

Méthodologie
L'étude a utilisé un émetteur de plasma Josephson (JPE) fabriqué à partir d'un monocristal de Bi2212, doté d'antennes à patch triangulaires. Le dispositif a été polarisé par une tension continue $V_0$ et un signal superposé à modulation de fréquence $V_m \sin(2\pi f_m t)$, où la fréquence de modulation $f_m$ a été variée de 3 GHz jusqu'à 100 kHz.

Les spectres de rayonnement ont été mesurés à l'aide d'un spectromètre infrarouge à transformée de Fourier (FTIR) de type Martin-Puplet équipé d'un bolomètre au silicium. Une avancée méthodologique clé de cette étude a consisté en un traitement rigoureux des données d'interférogramme. Contra irement aux itérations précédentes, les auteurs ont symétrisé les données d'interférogramme asymétriques (en utilisant une technique de repliement "OR" par rapport au pic central) et ont appliqué un remplissage nul (zero-padding) à une longueur de $2^{17}$ avant d'effectuer une transformée de Fourier rapide (FFT) discrète avec apodisation de Hanning. Cela a considérablement amélioré la résolution en fréquence et la précision spectrale.

L'analyse s'est concentrée sur :

1. Spectres non modulés ($I_{UM}(\omega)$) : Caractérisation de l'intensité du rayonnement sous polarisation DC pour établir une ligne de base.
2. Spectres FM : Analyse de l'évolution spectrale à mesure que $f_m$ diminuait, spécifiquement en examinant la largeur de bande ($B$) et le rapport d'intensité spectrale sous la fréquence porteuse ($\eta$).
3. Modélisation numérique : Développement d'un modèle où l'amplitude de l'intensité du rayonnement est retardée par un temps de relaxation de synchronisation $\tau_s$ par rapport à la tension instantanée, tandis que la fréquence répond instantanément.

Résultats clés

• Structure de double résonance dans le rayonnement non modulé : L'intensité du rayonnement non modulé $I_{UM}(\omega)$ présente une structure distincte à double pic de type Gaussienne. Un pic est centré à 849 GHz (étroit, FWHM 5,0 GHz) et l'autre à 866 GHz (large, FWHM 30 GHz). Les auteurs attribuent le pic large à la structure de l'antenne patch et le pic étroit à la résonance de la cavité de la mesa (spécifiquement le mode $TM_{1,3}$). L'observation simultanée de ces modes suggère des transitions temporelles ou des bifurcations rapides entre les modes de résonance sur une échelle de temps sub-nanoseconde, plutôt qu'une inhomogénéité spatiale ou des fluctuations de la tension de polarisation.
• Asymétrie spectrale et coupure : Dans les spectres FM, une asymétrie significative a été observée selon la fréquence de modulation. À haute $f_m$ (par exemple, 3 GHz), les bandes latérales inférieures étaient plus proéminentes que les bandes latérales supérieures, même lorsque la fréquence porteuse était positionnée là où l'intensité non modulée $I_{UM}(\omega)$ était proche de zéro. Cela contredit l'attente standard où le spectre FM est simplement le produit de l'intensité non modulée et d'un peigne de fréquences.
• Estimation quantitative du temps de relaxation : En analysant la largeur de bande spectrale et le rapport d'intensité $\eta$ en fonction de $f_m$, les auteurs ont identifié une région de transition critique autour de $f_m \sim 1$ GHz. Un modèle simplifié supposant une dépendance en fonction de saut (step-function) de l'intensité du rayonnement sur la fréquence, combiné à un délai de synchronisation fini $\tau_s$, a réussi à reproduire les données expérimentales. Cette analyse a permis d'obtenir un temps de relaxation de synchronisation de $\tau_s \simeq 0,28$ ns à 21,5 K.
• Reproduction numérique : Des simulations numériques intégrant ce délai ($\tau_s = 0,28$ ns) ont réussi à reproduire la transformation spectrale observée, spécifiquement l'augmentation de l'intensité au bord de la fréquence inférieure lorsque la fréquence porteuse était située en dessous du pic de résonance principal. Cela confirme que, bien que la fréquence du rayonnement réponde instantanément à la modulation de tension, l'intensité du rayonnement nécessite un temps fini pour se synchroniser.

Signification et revendications
L'article affirme fournir la première estimation précise et quantitative du temps de relaxation de synchronisation dans les IJJ empilées. La contribution primaire est la démonstration que la synchronisation macroscopique des oscillations de plasma Josephson n'est pas instantanée mais régie par une dynamique hors équilibre sub-nanoseconde ($\tau_s \approx 0,28$ ns).

Les auteurs soutiennent que ce temps de relaxation correspond à la durée de vie du plasmon macroscopique excité au sein de la structure de la mesa, distincte de la durée de vie microscopique des plasmas Josephson observée dans la résonance millimétrique. L'étude souligne que l'interaction entre la mesa de l'IJJ et les antennes patch crée un paysage de résonance complexe où les transitions de mode se produisent rapidement. En introduisant une analyse résolue en temps via les spectres FM, ils démontrent une méthode pour extraire les paramètres dynamiques des systèmes de jonctions Josephson couplés, offrant une compréhension plus profonde des mécanismes de synchronisation dans ces systèmes quantiques macroscopiques. Le travail suggère que les conceptions futures de dispositifs JPE doivent tenir compte de ces interactions entre l'antenne et la mesa ainsi que du temps de synchronisation fini pour optimiser les propriétés de rayonnement.