Systematic Schrieffer-Wolff-transformation approach to Josephson junctions: quasiparticle effects and Josephson harmonics

This paper employs a systematic Schrieffer-Wolff transformation to derive an effective Hamiltonian for Josephson junctions that recovers the conventional cosine potential while simultaneously revealing how Bogoliubov quasiparticles induce correlated dynamics and how higher-order terms naturally generate Josephson harmonics linked to microscopic junction properties.

The Gist

Imagine a Josephson junction as a very special, super-conducting bridge between two islands. In the ideal world, this bridge is crossed only by "Cooper pairs"—these are like perfectly synchronized dance couples (two electrons holding hands) that glide across the tunnel without any friction. This smooth, synchronized crossing is what makes superconducting quantum computers work.

However, sometimes the dance couples break up. Individual electrons, now called "quasiparticles," get left behind. These lone dancers are messy; they don't follow the rhythm, and when they try to cross the bridge, they disrupt the perfect flow of the couples. This is known as "quasiparticle poisoning," and it's a headache for scientists trying to build stable quantum devices.

This paper introduces a new mathematical tool called the Schrieffer-Wolff Transformation (SWT). Think of this tool as a sophisticated translator or a "lens" that allows physicists to look at the messy, complex reality of individual electrons and translate it into a simpler, effective story about the whole system.

Here is what the authors discovered using this lens:

1. Recovering the Classic Story (The Baseline)
First, the authors used their tool on a "clean" bridge where no lone dancers (quasiparticles) are present. By starting with the complex rules of individual electron tunneling and applying their transformation, they successfully recreated the famous, simple equation used by everyone in the field: $H = -E_J \cos(\phi)$.

• The Analogy: It's like taking a chaotic crowd of people moving randomly and mathematically proving that, on average, they move like a single, smooth wave. This confirmed their tool works and connects the microscopic chaos to the macroscopic order.

2. The Messy Reality: When Lone Dancers Join the Party
Next, they relaxed the rules and allowed a single "lone dancer" (a quasiparticle) to exist on the bridge.

• The Discovery: They found that the lone dancer doesn't just sit there; they interact with the dancing couples. The movement of the lone dancer becomes "entangled" with the movement of the couples.
• The Result: This interaction changes the energy landscape of the bridge. In their simple "toy model" (a simplified version of the bridge), they showed that the presence of a lone dancer shifts the "sweet spot" where the system is most stable and changes the "stiffness" (curvature) of the energy curve.
• Why it matters: In a real quantum computer, this means the presence of these lone dancers changes the frequency at which the qubit (the computer bit) vibrates. It's like if a single person walking on a trampoline changed the bounce frequency for everyone else jumping on it.

3. Discovering Hidden Harmonics
Finally, the authors used their tool to look even deeper, going beyond the standard second-level math to a fourth-level calculation.

• The Discovery: They found that the bridge doesn't just have one simple rhythm (the main cosine wave). It has "harmonics"—subtle, higher-frequency ripples in the energy landscape.
• The Connection: The size of these ripples isn't random; it is directly tied to the microscopic details of the materials used to build the bridge.
• The Benefit: Their math provides a recipe to calculate exactly how strong these ripples are based on the specific properties of the superconducting leads. This could help engineers tune their devices to control these harmonics if they choose to.

In Summary
The paper doesn't propose a new device or a medical cure. Instead, it provides a better map.

• It confirms that the standard map (the simple cosine equation) is a valid approximation of the complex reality.
• It draws a new, more detailed map that shows exactly how "messy" lone electrons distort the path of the "clean" dancing couples.
• It reveals hidden "ripples" (harmonics) in the path and explains how to calculate their size based on the materials used.

Essentially, the authors have built a systematic way to translate the complex, messy language of individual electrons into the clean, effective language used to design and understand superconducting quantum circuits.

Technical Summary

Technical Summary: Systematic Schrieffer-Wolff-transformation approach to Josephson junctions: quasiparticle effects and Josephson harmonics

Problem Statement
Josephson junctions (JJs) are fundamental to superconducting quantum devices, typically modeled by the effective Hamiltonian $H = -E_J \cos \hat{\phi}$, which relies on the coherent tunneling of Cooper pairs (CPs). However, in realistic scenarios, the system is subject to "quasiparticle (QP) poisoning," where elementary excitations (QPs) tunnel across the junction, disrupting coherent behavior. While experimental efforts focus on mitigating QP populations, theoretical challenges remain in describing the dynamics of these inevitable QPs within effective models for devices like transmons and flux qubits. Conventional approaches often rely on master-equation methods that treat QPs and CPs separately or lack a direct operator-level connection to the underlying microscopic theory. Furthermore, standard perturbation theory yields corrections to energies and wavefunctions, requiring a challenging reconstruction of the effective Hamiltonian.

Methodology
The authors employ the Schrieffer-Wolff transformation (SWT), a systematic perturbative method involving a unitary transformation followed by a projection onto a relevant low-energy subspace. The starting point is the charge-conserving Bardeen-Cooper-Schrieffer (BCS) theory, where the order parameter is defined via a condensate creation operator $S^+$, ensuring the total particle number is conserved. The system consists of two superconducting leads separated by a tunnel barrier, described by a total Hamiltonian $H = H_0 + V$, where $H_0$ represents the leads and $V$ represents single-electron tunneling terms expressed in terms of Bogoliubov QP operators.

The authors perform the SWT under two distinct regimes:

1. Second-order expansion with QP-free subspace: Restricting the low-energy subspace to states without QPs to recover the conventional effective model.
2. Second-order expansion with single-QP subspace: Relaxing the restriction to include states containing exactly one QP, thereby deriving an effective model for a junction populated by a single QP.
3. Fourth-order expansion with QP-free subspace: Returning to the QP-free subspace but extending the perturbative expansion to fourth order to capture higher-order tunneling processes.

Key Contributions and Results

• Recovery of the Conventional Model: By performing a second-order SWT on the charge-conserving BCS theory restricted to QP-free states, the authors directly recover the standard effective Hamiltonian $H_{\text{eff}} = -E_J \cos \hat{\phi}$. Crucially, this derivation provides an explicit operator-valued phase bias $\hat{\phi}$ and demonstrates the unitary transformation connecting the microscopic single-electron tunneling terms to the macroscopic Josephson term. This contrasts with conventional perturbation theory, which typically yields energy corrections rather than operator expressions.

• Quasiparticle-Coupled Dynamics: When the low-energy subspace is extended to include single-QP states, the effective Hamiltonian reveals new terms describing the coupling between QP motion and CP tunneling. The derived effective Hamiltonian (Eq. 18) includes:

  • First-order terms representing direct QP tunneling between leads.
  • Second-order terms describing intralead QP scattering, which can be assisted by CP tunneling (e.g., a QP scattering from mode $k$ to $k'$ accompanied by a CP transfer).
  • These terms introduce correlated dynamics that reshape the energy spectrum.

• Impact on Qubit Excitation Frequencies: Using a minimal toy model with a single energy level at the gap edge, the authors analyze the spectrum of the junction with a single QP. The results show that the presence of a QP:

  • Introduces a vertical energy shift (on the order of the gap $\Delta$).
  • Shifts the minimum of the energy landscape in phase space ($\phi$) away from zero.
  • Modifies the curvature of the potential well near the minimum.
    Since the curvature determines the qubit excitation frequency, the presence of QPs alters this frequency, a result that would require a full circuit QED model to quantify precisely but is qualitatively established here.

• Josephson Harmonics: By performing a fourth-order SWT within the QP-free subspace, the authors derive expressions for Josephson harmonics. These higher-order terms naturally arise from the perturbative expansion and describe amplitudes directly related to the microscopic properties of the superconducting leads and the junction. The authors provide exact expressions for calculating these effective-model parameters from microscopic quantities, offering a framework for tuning the ratio between different harmonics.

Significance and Claims
The paper claims that the SWT approach offers a systematic and rigorous framework for deriving effective JJ Hamiltonians from microscopic theories. Its primary advantages are:

1. Direct Operator Derivation: It provides operator expressions for the effective Hamiltonian directly, avoiding the reconstruction step required by standard perturbation theory.
2. Systematic Extensibility: The method can be systematically extended to more complex scenarios, such as including QPs or calculating higher-order harmonics, by simply adjusting the low-energy subspace or the order of the expansion.
3. Microscopic-Macroscopic Connection: It explicitly highlights the unitary transformation between the microscopic BCS theory and the resulting phase-dependent Hamiltonian.

The authors emphasize that their work extends the conventional $H = -E_J \cos \phi$ model to account for QP presence and higher-order harmonics, providing a theoretical basis for understanding how QP poisoning and junction microstructure influence the energy spectrum and dynamics of superconducting quantum devices. They note that while the voltage dependence of the Josephson energy appears in their derivation, it is typically a minor correction in practical applications.