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

This study analyzes the spectral dynamics of frequency-modulated terahertz radiation from stacked intrinsic Josephson junctions in a cuprate superconductor, revealing a double-peak structure driven by electromagnetic coupling and quantifying a sub-nanosecond synchronization time that governs the system's non-equilibrium behavior.

The Gist

Imagine a superconductor not as a solid block of metal, but as a towering skyscraper made of thousands of tiny, identical floors. Each floor is a "Josephson junction," a microscopic sandwich that allows electricity to flow without resistance. When you push a specific amount of voltage through this skyscraper, all these floors start to "dance" together, vibrating in perfect unison. This collective dance emits a beam of invisible light called Terahertz radiation, which sits between microwaves and infrared light on the spectrum.

This paper is like a high-speed camera study of that dance, specifically looking at what happens when you try to change the rhythm of the music very quickly.

The Setup: A Dance Floor with Two Favorite Spots

The researchers built their "skyscraper" (a device made of a material called Bi2212) and attached special triangular antennas to the top, like satellite dishes.

When they let the device dance at a steady pace (using a constant voltage), they expected to see one single, clear peak in the frequency of the light it emitted. Instead, they found two distinct peaks in the intensity of the light.

• The Analogy: Imagine a singer who has two favorite notes they can hit perfectly. Even when they try to sing just one note, their voice naturally splits into a duet of two very close pitches.
• The Cause: The paper suggests this happens because the device has two different ways to resonate: one way is determined by the shape of the "skyscraper" itself (the mesa), and the other is determined by the "satellite dishes" (the antennas) attached to it. Sometimes the device dances to the rhythm of the building; other times, it dances to the rhythm of the antennas. At certain voltages, it's hard to tell which one it's doing, so you see both.

The Experiment: Changing the Beat

To study how quickly these dancers can change their rhythm, the researchers applied a "frequency modulation." Think of this as taking a steady drumbeat and wiggling the tempo up and down very fast.

• They wiggled the tempo at different speeds: from very slow (100,000 wiggles per second) to very fast (3 billion wiggles per second).
• The Expectation: If the dancers were perfect and instantaneous, the light they emitted should simply spread out into a neat pattern of "comb teeth" (a series of evenly spaced frequencies) that perfectly matches the shape of the original two peaks.

The Surprise: The "Lag" in the Dance

Here is where the paper gets interesting. When they wiggled the tempo at a medium speed (around 1 billion wiggles per second), the pattern of light broke the rules.

• The light didn't just spread out evenly. Instead, the "comb teeth" on the low-frequency side became much brighter than the ones on the high-frequency side, even though the original dance floor preferred the high frequencies.
• The Analogy: Imagine a group of dancers who are told to speed up and slow down. If they are perfect, they change speed instantly. But if they are human, they have a reaction time. If you tell them to slow down, they keep speeding up for a tiny fraction of a second before they actually stop.
• The Discovery: The researchers found that the intensity of the light (how bright the dance is) has a "reaction time" or a lag. The frequency of the light changes instantly with the voltage, but the brightness takes a tiny moment to catch up.

The Result: Measuring the Lag

By analyzing exactly how the light pattern distorted at different wiggle speeds, the team calculated this lag time.

• They found the synchronization relaxation time (the time it takes for the whole skyscraper to agree on a new brightness level) is about 0.28 nanoseconds.
• To put that in perspective: A nanosecond is a billionth of a second. 0.28 nanoseconds is so fast that light only travels about 8 centimeters (3 inches) in that time. Yet, for these quantum dancers, that tiny delay is huge and completely changes the shape of the light they emit.

Why It Matters (According to the Paper)

The paper claims that by understanding this tiny delay, they can now mathematically reproduce the strange, distorted light patterns they saw. They showed that if you assume the dancers have this specific 0.28-nanosecond "thinking time," your computer model perfectly matches the real-world experiment.

In short, this paper didn't just watch the dancers; it figured out exactly how long it takes for the entire building to agree on a new move, revealing that even in a world of super-fast quantum physics, there is still a tiny, measurable moment of hesitation.

Technical Summary

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

Problem Statement
Intrinsic Josephson junctions (IJJs) in cuprate superconductors, specifically Bi$_2$Sr$_2$CaCu$_2$O$_{8+\delta}$ (Bi2212), serve as macroscopic quantum systems for studying interacting nonlinear oscillators. While these systems can generate terahertz (THz) radiation via the AC Josephson effect, a comprehensive quantitative understanding of their dynamic synchronization mechanisms—particularly under frequency modulation (FM)—has remained elusive. Previous work by the authors demonstrated FM spectra and attributed observed phenomena to a finite synchronization relaxation time, yet a precise estimation of this time and a detailed discussion of the underlying mechanism were lacking. The complexity arises from the coexistence of two distinct interactions: inductive coupling across hundreds of junctions (characterized by the penetration length $\lambda_{ab}$) and capacitive coupling between nearest neighbors (characterized by the Debye length $\mu$).

Methodology
The study utilized a Josephson plasma emitter (JPE) fabricated from a Bi2212 single crystal, featuring triangular patch antennas. The device was biased with a DC voltage $V_0$ and a superimposed frequency-modulated signal $V_m \sin(2\pi f_m t)$, where the modulation frequency $f_m$ was varied from 3 GHz down to 100 kHz.

Radiation spectra were measured using a Martin-Puplet-type Fourier transform infrared (FTIR) spectrometer equipped with a silicon bolometer. A key methodological advancement in this study involved rigorous data processing of the interferogram. Unlike previous iterations, the authors symmetrized the asymmetric interferogram data (using an "OR" aliasing technique relative to the center burst) and applied zero-padding to a length of $2^{17}$ before performing a discrete Fast Fourier Transform (FFT) with Hanning apodization. This significantly improved frequency resolution and spectral precision.

The analysis focused on:

1. Unmodulated Spectra ($I_{UM}(\omega)$): Characterizing the radiation intensity under DC bias to establish a baseline.
2. FM Spectra: Analyzing the spectral evolution as $f_m$ was decremented, specifically looking at bandwidth ($B$) and the ratio of spectral intensity below the carrier frequency ($\eta$).
3. Numerical Modeling: Developing a model where the radiation intensity amplitude is delayed by a synchronization relaxation time $\tau_s$ relative to the instantaneous voltage, while the frequency responds instantaneously.

Key Results

• Double-Resonance Structure in Unmodulated Radiation: The unmodulated radiation intensity $I_{UM}(\omega)$ exhibits a distinct double-Gaussian peak structure. One peak is centered at 849 GHz (narrow, FWHM 5.0 GHz) and the other at 866 GHz (broad, FWHM 30 GHz). The authors attribute the broad peak to the patch antenna structure and the narrow peak to the mesa cavity resonance (specifically the $TM_{1,3}$ mode). The simultaneous observation of these modes suggests rapid temporal transitions or bifurcations between resonance modes on a sub-nanosecond timescale, rather than spatial inhomogeneity or bias voltage fluctuations.
• Spectral Asymmetry and Cutoff: In the FM spectra, a significant asymmetry was observed depending on the modulation frequency. At high $f_m$ (e.g., 3 GHz), the lower sidebands were more prominent than the upper sidebands, even when the carrier frequency was positioned where the unmodulated intensity $I_{UM}(\omega)$ was near zero. This contradicts the standard expectation where the FM spectrum is simply the product of the unmodulated intensity and a frequency comb.
• Quantitative Estimation of Relaxation Time: By analyzing the spectral bandwidth and the intensity ratio $\eta$ as a function of $f_m$, the authors identified a critical transition region around $f_m \sim 1$ GHz. A simplified model assuming a step-function dependence of radiation intensity on frequency, combined with a finite synchronization delay $\tau_s$, successfully reproduced the experimental data. This analysis yielded a synchronized relaxation time of $\tau_s \simeq 0.28$ ns at 21.5 K.
• Numerical Reproduction: Numerical simulations incorporating this delay ($\tau_s = 0.28$ ns) successfully reproduced the observed spectral transformation, specifically the enhanced intensity at the lower frequency edge when the carrier frequency was below the main resonance peak. This confirms that while the radiation frequency responds instantaneously to voltage modulation, the radiation intensity requires a finite time to synchronize.

Significance and Claims
The paper claims to provide the first precise, quantitative estimation of the synchronization relaxation time in stacked IJJs. The primary contribution is the demonstration that the macroscopic synchronization of Josephson plasma oscillations is not instantaneous but governed by a sub-nanosecond non-equilibrium dynamic ($\tau_s \approx 0.28$ ns).

The authors assert that this relaxation time corresponds to the lifetime of the macroscopic plasmon excited within the mesa structure, distinct from the microscopic lifetime of Josephson plasmons observed in millimeter-wave resonance. The study highlights that the interaction between the IJJ mesa and the patch antennas creates a complex resonance landscape where mode transitions occur rapidly. By introducing time-resolved analysis via FM spectra, the authors demonstrate a method to extract dynamical parameters of coupled Josephson junction systems, offering a deeper understanding of the synchronization mechanisms in these macroscopic quantum systems. The work suggests that future designs of JPE devices must account for these antenna-mesa interactions and the finite synchronization time to optimize radiation properties.