Magnetic flux controlled current phase relationship in double Quantum Dot Josephson junction

This study investigates a double quantum dot Josephson junction threaded by magnetic flux, revealing how flux and interaction strength govern complex phase transitions between singlet, doublet, and triplet ground states and significantly influence the critical current.

The Gist

The Big Picture: A Quantum Traffic Jam

Imagine a superhighway made of electricity, but instead of cars, it's made of electrons. Usually, electrons bump into each other and get stuck. But in a superconductor (the "highway" in this story), electrons pair up and glide smoothly without any friction. This is the Josephson effect: a supercurrent flowing effortlessly across a bridge.

The scientists in this paper built a tiny, artificial bridge with two small "parking spots" (called Quantum Dots) in the middle. They wanted to see how these parking spots affect the flow of the supercurrent, especially when they added a "magnetic twist" (magnetic flux) to the road.

The Tools: A Digital Model

The real physics here is incredibly complex, like trying to calculate the movement of every single grain of sand on a beach. To solve this, the researchers used a clever trick called a "Surrogate Model."

Think of the superconducting highway as an infinite ocean. You can't simulate an infinite ocean on a computer. So, they replaced the infinite ocean with just three specific islands (discrete energy levels). They adjusted the "bridges" connecting these islands to the parking spots so that, mathematically, the three islands acted exactly like the infinite ocean. This allowed them to solve the equations perfectly on a computer without needing supercomputers.

The Experiments: What Happened on the Bridge?

The researchers tested three different scenarios, like changing the rules of a game:

1. The Empty Parking Spots (No Interaction)

First, they looked at the dots when the electrons didn't bother each other.

• The Magic Trick: They applied a magnetic twist (flux) to the loop.
• The Result: When the twist was just right (a specific angle), the two paths the electrons could take canceled each other out perfectly, like two waves crashing and making a flat line. The current stopped completely.
• The Analogy: Imagine two people walking around a circular track in opposite directions. If they start at the same time and the track is twisted just right, they meet exactly in the middle and stop. The traffic halts.

2. One Grumpy Spot (One Dot has Interaction)

Next, they made one parking spot "grumpy" (adding Coulomb interaction), meaning electrons there didn't like to be alone and preferred to pair up or avoid each other.

• The Dance: The system started to switch between two "moods" or states:
  • The Singlet (The Couple): Two electrons holding hands, moving together.
  • The Doublet (The Soloist): One electron moving alone.
• The Magnetic Control: By turning the magnetic twist knob, they could force the system to switch from "Couple" mode to "Soloist" mode.
• The Surprise: The magnetic field didn't just change the mood; it shifted when the switch happened. It was like a magnetic lever that pushed the transition point back and forth.

3. Two Grumpy Spots (Both Dots have Interaction)

Finally, they made both parking spots grumpy. This is where things got really interesting.

• The Three-Act Play: As they changed the phase difference (the "traffic signal"), the ground state (the most comfortable state for the electrons) went through a three-stage evolution:
  1. Doublet (Soloist)
  2. Singlet (Couple)
  3. Triplet (A weird, high-energy trio state)
• The Magnetic Brake: When they increased the magnetic twist, it acted like a brake on the "Soloist" and "Triplet" states. Eventually, the magnetic field was so strong that it forced the system to stay in the "Couple" (Singlet) state no matter what.
• The Critical Peak: Even though the magnetic field usually stops the current, they found a special "sweet spot" in the interaction strength where the current suddenly spiked to a maximum. It was like finding a specific gear in a car where the engine roars to its loudest, most efficient point.

The "Triple Point"

In the final scenario, where the connection between the dots and the highway was very strong (stronger than the superconducting gap itself), the researchers found a "Triple Point."

• The Analogy: Think of a map where three different countries meet at a single point. In this quantum map, there is a specific combination of magnetic twist and traffic signal where the "Soloist," "Triplet," and "Couple" states all merge into one. It's a rare moment where three different quantum realities collide.

Summary

The paper shows that by using a clever computer model (the surrogate), they could map out exactly how magnetic fields and electron interactions change the flow of supercurrents in a double-dot system. They discovered that magnetic fields can act as a precise switch, toggling the system between different quantum states (Soloist vs. Couple) and even creating a special "sweet spot" where the current is strongest. They also found a rare "triple point" where three different quantum states meet.

Technical Summary

1. Problem Statement

The paper investigates the transport properties and quantum phase transitions in a Josephson junction where the central region consists of two parallel-coupled Quantum Dots (QDs) connected to two superconducting (SC) leads. The system is threaded by a magnetic flux ($\Phi_B$), introducing an Aharonov-Bohm phase.

The study aims to address several open questions in this regime:

• How does magnetic flux control the ground state parity (singlet, doublet, triplet) and the resulting Josephson current?
• How do Coulomb interactions (on-site repulsion $U$) in one or both dots influence the $0-\pi$ transition and the critical current?
• What is the behavior of the system in the strong tunneling regime ($\Gamma > \Delta$), where perturbative methods fail?
• Existing studies often rely on weak-tunneling perturbation theory or complex numerical methods like Numerical Renormalization Group (NRG). The authors seek a method that balances computational efficiency with accuracy across a broad parameter space.

2. Methodology

The authors employ a combination of three theoretical approaches to tackle the problem:

A. Surrogate Hamiltonian Model (Primary Method)

To handle the strong correlations and the continuum of the superconducting leads without the high computational cost of NRG, the authors utilize a surrogate model (related to the zero-bandwidth approximation).

• Discretization: The continuous SC leads are discretized into a finite number of effective energy levels (specifically three levels).
• Renormalization: The tunneling coefficients are renormalized to preserve the spectral properties of the original continuum.
• Implementation: The continuous momentum $k$ is replaced by discrete indices $l$. The resulting Hamiltonian has a finite dimension ($4^8 = 65,536$ states for the full system), allowing for exact diagonalization.
• Advantage: This method provides direct access to the full energy spectrum, eigenstates, and physical quantities (entropy, parity, current) with high accuracy, validated against NRG results in single-dot cases.

B. Path-Integral Formalism

For the non-interacting limit ($U=0$), the authors employ a functional integral approach.

• They integrate out the superconducting leads to derive an effective action involving a self-energy term.
• The Josephson current is calculated via the derivative of the partition function with respect to the SC phase difference ($\phi$).
• This serves as a benchmark to validate the surrogate model results in the non-interacting regime.

C. Low-Energy Effective Hamiltonian

To gain physical insight, the authors derive a low-energy effective model in the limit of infinite superconducting gap ($\Delta \to \infty$).

• This model neglects quasiparticle excitations in the leads and focuses on the effective coupling between the dots.
• It allows for analytical solutions (using $3\times3$ or $4\times4$ matrices) to explain the mechanisms behind phase transitions and current profiles.

3. Key Contributions and Results

The paper categorizes its findings based on the presence and symmetry of Coulomb interactions ($U_1, U_2$):

Case 1: Non-Interacting Limit ($U_1 = U_2 = 0$)

• Symmetric Dots: When dot energies are symmetric ($\epsilon_1 = \epsilon_2$), the Josephson current vanishes at $\Phi_B = \pi$ due to destructive interference. The system remains in a singlet state (parity 0) throughout; no $0-\pi$ transition occurs, only a suppression of current.
• Asymmetric Dots: If dot energies differ, the current remains non-zero at $\Phi_B = \pi$. A $0-\pi$ transition can be induced by varying $\Phi_B$, driven by the crossing of subgap Andreev Bound States (ABS).
• Current Components: The total current is dominated by the discrete ABS contribution, while the continuum contribution is minor and less sensitive to flux.

Case 2: Single Dot Interaction ($U_1 \neq 0, U_2 = 0$)

• Phase Transitions: The system exhibits a magnetic-flux-controlled transition between a Singlet ($|S\rangle$, parity 0) and a Doublet state.
• Flux Dependence: The phase boundary between Singlet and Doublet shifts almost linearly with the magnetic flux $\Phi_B$.
• Current Behavior: Unlike the standard $0-\pi$ transition, the doublet state contributes both positive and negative current components depending on the phase and flux, leading to complex current profiles.
• Triplet State: The triplet state does not become the ground state in this configuration.

Case 3: Double Dot Interaction ($U_1 = U_2 \neq 0$)

• RKKY Interaction: The presence of interactions in both dots introduces Ruderman-Kittel-Kasuya-Yosida (RKKY) interactions, enabling transitions involving the Triplet ($|T\rangle$) state.
• Sequential Transitions: For weak interactions, the ground state evolves from Doublet $\to$ Singlet $\to$ Triplet as the SC phase difference $\phi$ increases.
• Strong Interaction Limit: As $U$ increases, the Doublet phase is suppressed. The system transitions directly from Singlet to Triplet.
• Magnetic Flux Effect: Increasing $\Phi_B$ suppresses both Doublet and Triplet phases, stabilizing the Singlet ground state even at $\Phi_B = \pi$.
• Critical Current Peak: In the regime where the ground state is Singlet, increasing $U$ does not induce a Singlet-Doublet transition. Instead, it drives a transition between upper and lower Singlet states. This results in a critical current peak at a specific interaction strength ($U^* \approx 0.6735$), where the critical phase is pinned at $\phi = \pi$. This is distinct from the $0-\pi$ transition; it is an internal $0-0$ transition.

Case 4: Strong Tunneling Regime ($\Gamma > \Delta$)

• Triple Point: When the tunneling rate exceeds the superconducting gap, single-particle transport becomes dominant. The Doublet phase region expands and touches the Triplet region.
• Convergence: In the presence of magnetic flux, the Singlet, Doublet, and Triplet phases converge at a specific point in the $(\phi, \Phi_B)$ parameter space, forming a triple point.
• State Ordering: The energy level ordering at $\phi=\pi$ changes from $|S\rangle \to \text{Doublet} \to |T\rangle$ (weak tunneling) to $|T\rangle \to \text{Doublet} \to |S\rangle$ (strong tunneling).

4. Significance

• Methodological Advancement: The paper demonstrates that the surrogate Hamiltonian approach is a powerful, computationally efficient alternative to NRG for studying strongly correlated Josephson junctions with complex geometries (like DQDs with flux). It achieves quantitative agreement with exact methods while being more tractable for parameter sweeps.
• Novel Physics: The discovery of a critical current peak driven by an internal Singlet-Singlet transition (rather than a Singlet-Doublet transition) in the presence of magnetic flux provides new insight into the interplay between Coulomb interactions and superconductivity.
• Control Mechanisms: The work establishes magnetic flux as a precise tuning knob for controlling ground state parity and critical currents in double-dot systems, suggesting potential applications in superconducting qubits and quantum phase transition engineering.
• Triple Point Prediction: The identification of a triple point in the parameter space for the strong tunneling regime offers a specific experimental signature for future verification in nanowire or carbon nanotube-based devices.

In summary, this work provides a comprehensive theoretical framework for understanding how magnetic flux and Coulomb interactions jointly dictate the quantum phase and transport characteristics of parallel double quantum dot Josephson junctions, bridging the gap between weak and strong coupling regimes.