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

Resumo Técnico: Análise de sincronização resolvida no tempo de junções de Josephson intrínsecas empilhadas de um supercondutor de cuprato com espectros de radiação terahertz modulados em frequência

Problema
As junções de Josephson intrínsecas (IJJs) em supercondutores de cuprato, especificamente o Bi$_2$Sr$_2$CaCu$_2$O$_{8+\delta}$ (Bi2212), servem como sistemas quânticos macroscópicos para o estudo de osciladores não lineares interagentes. Embora esses sistemas possam gerar radiação terahertz (THz) via o efeito Josephson AC, uma compreensão quantitativa abrangente de seus mecanismos de sincronização dinâmica — particularmente sob modulação de frequência (FM) — permaneceu elusiva. Trabalhos anteriores dos autores demonstraram espectros FM e atribuíram os fenômenos observados a um tempo de relaxação de sincronização finito, porém uma estimativa precisa deste tempo e uma discussão detalhada do mecanismo subjacente eram ausentes. A complexidade surge da coexistência de duas interações distintas: o acoplamento indutivo através de centenas de junções (caracterizado pelo comprimento de penetração $\lambda_{ab}$) e o acoplamento capacitivo entre vizinhos imediatos (caracterizado pelo comprimento de Debye $\mu$).

Metodologia
O estudo utilizou um emissor de plasma de Josephson (JPE) fabricado a partir de um cristal único de Bi2212, apresentando antenas de patch triangulares. O dispositivo foi polarizado com uma tensão DC $V_0$ e um sinal sobreposto de modulação de frequência $V_m \sin(2\pi f_m t)$, onde a frequência de modulação $f_m$ foi variada de 3 GHz até 100 kHz.

Os espectros de radiação foram medidos usando um espectrômetro de infravermelho por transformada de Fourier (FTIR) do tipo Martin-Puplet equipado com um bolômetro de silício. Um avanço metodológico chave neste estudo envolveu o processamento rigoroso dos dados do interferograma. Diferentemente de iterações anteriores, os autores simetrizaram os dados de interferograma assimétricos (usando uma técnica de aliasing "OR" em relação ao surto central) e aplicaram zero-padding para um comprimento de $2^{17}$ antes de realizar uma Transformada Rápida de Fourier (FFT) discreta com apodização Hanning. Isso melhorou significativamente a resolução de frequência e a precisão espectral.

A análise focou em:

1. Espectros Não Modulados ($I_{UM}(\omega)$): Caracterização da intensidade de radiação sob polarização DC para estabelecer uma linha de base.
2. Espectros FM: Análise da evolução espectral conforme $f_m$ era decrementada, especificamente observando a largura de banda ($B$) e a razão da intensidade espectral abaixo da frequência portadora ($\eta$).
3. Modelagem Numérica: Desenvolvimento de um modelo onde a amplitude da intensidade de radiação é atrasada por um tempo de relaxação de sincronização $\tau_s$ em relação à tensão instantânea, enquanto a frequência responde instantaneamente.

Principais Resultados

• Estrutura de Dupla Ressonância na Radiação Não Modulada: A intensidade de radiação não modulada $I_{UM}(\omega)$ exibe uma estrutura distinta de pico duplo-Gaussiana. Um pico está centrado em 849 GHz (estreito, FWHM 5,0 GHz) e o outro em 866 GHz (largo, FWHM 30 GHz). Os autores atribuem o pico largo à estrutura da antena de patch e o pico estreito à ressonância da cavidade da mesa (especificamente o modo $TM_{1,3}$). A observação simultânea desses modos sugere transições temporais ou bifurcações rápidas entre modos de ressonância em uma escala de tempo sub-nanossegundo, em vez de inhomogeneidade espacial ou flutuações de tensão de polarização.
• Assimetria Espectral e Corte: Nos espectros FM, uma assimetria significativa foi observada dependendo da frequência de modulação. Em altas $f_m$ (ex: 3 GHz), os subbandas inferiores foram mais proeminentes do que as subbandas superiores, mesmo quando a frequência portadora estava posicionada onde a intensidade não modulada $I_{UM}(\omega)$ estava próxima de zero. Isso contradiz a expectativa padrão onde o espectro FM é simplesmente o produto da intensidade não modulada e um pente de frequência.
• Estimativa Quantitativa do Tempo de Relaxação: Ao analisar a largura de banda espectral e a razão de intensidade $\eta$ como função de $f_m$, os autores identificaram uma região de transição crítica em torno de $f_m \sim 1$ GHz. Um modelo simplificado assumindo uma dependência de função degrau da intensidade de radiação em relação à frequência, combinado com um atraso de sincronização finito $\tau_s$, reproduziu com sucesso os dados experimentais. Esta análise resultou em um tempo de relaxação de sincronização de $\tau_s \simeq 0,28$ ns a 21,5 K.
• Reprodução Numérica: Simulações numéricas incorporando este atraso ($\tau_s = 0,28$ ns) reproduziram com sucesso a transformação espectral observada, especificamente o aumento da intensidade na borda de frequência inferior quando a frequência portadora estava abaixo do pico principal de ressonância. Isso confirma que, embora a frequência de radiação responda instantaneamente à modulação de tensão, a intensidade da radiação requer um tempo finito para se sincronizar.

Significância e Alegações
O artigo afirma fornecer a primeira estimativa precisa e quantitativa do tempo de relaxação de sincronização em IJJs empilhadas. A contribuição primária é a demonstração de que a sincronização macroscópica das oscilações de plasma de Josephson não é instantânea, mas governada por uma dinâmica de não-equilíbrio sub-nanossegundo ($\tau_s \approx 0,28$ ns).

Os autores afirmam que este tempo de relaxação corresponde ao tempo de vida do plasma macroscópico excitado dentro da estrutura da mesa, distinto do tempo de vida microscópico dos plasmas de Josephson observados na ressonância de milimétricas. O estudo destaca que a interação entre a mesa da IJJ e as antenas de patch cria um cenário de ressonância complexo onde as transições de modo ocorrem rapidamente. Ao introduzir a análise resolvida no tempo via espectros FM, os autores demonstram um método para extrair parâmetros dinâmicos de sistemas de junções de Josephson acopladas, oferecendo uma compreensão mais profunda dos mecanismos de sincronização nesses sistemas quânticos macroscópicos. O trabalho sugere que futuros designs de dispositivos JPE devem considerar essas interações antena-mesa e o tempo de sincronização finito para otimizar as propriedades de radiação.