🔬 mesoscale physics

Entropy transport through a superfluid quantum point contact: A Keldysh field-theory approach

Using the Keldysh formalism, this study derives particle and entropy current-bias characteristics for a superfluid quantum point contact connecting two Fermi-gas reservoirs, revealing an oscillatory entropy current at low voltages in the ballistic limit and comparing these theoretical findings with experimental data from unitary cold atomic gases.

Original authors: Davide Bertolusso, C. J. Bolech, Thierry Giamarchi

Published 2026-05-04

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Davide Bertolusso, C. J. Bolech, Thierry Giamarchi

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine two large, calm lakes filled with a special kind of "super-fluid" water. In this fluid, the tiny particles (atoms) move in perfect harmony, like a synchronized dance troupe, rather than bumping into each other chaotically. Now, imagine connecting these two lakes with a very narrow, single-lane bridge. This bridge is our "Quantum Point Contact."

The scientists in this paper are studying what happens when they push the water from one lake to the other through this bridge. They aren't just looking at how many water droplets (particles) move across; they are also measuring something more abstract called "entropy," which you can think of as the disorder or messiness of the flow.

Here is the breakdown of their discovery using simple analogies:

1. The Setup: The Dance Floor and the Bridge

The two lakes are kept at slightly different "pressure" levels (chemical potential). This pressure difference acts like a slope, encouraging the water to flow from the high-pressure lake to the low-pressure one.

The Particles: These are the water droplets trying to cross the bridge.
The Entropy: This is the "chaos" or "heat" carried along with the droplets.

2. The Special Rules of the Game (Superfluids)

In normal water, if you push a droplet across a bridge, it just goes straight through. But in this super-fluid, the particles are "entangled" in pairs (like dance partners holding hands).

The Barrier: There is a "dance floor" rule (the superconducting gap) that makes it hard for single dancers to cross unless they have enough energy.
The Trick (Andreev Reflection): If a single dancer tries to cross but hits the rule, they don't just bounce back. Instead, they grab a partner from the other side, turn into a "hole" (a missing dancer), and bounce back. This is called Andreev Reflection.
The Multi-Step Dance (MAR): If the pressure difference is just right, the dancer can do a complex routine: cross, bounce back, grab another partner, cross again, and so on. This is called Multiple Andreev Reflection (MAR). It's like a dancer doing a series of backflips and spins to get across the bridge.

3. The Big Discovery: The Oscillating Entropy

The scientists calculated two things:

Particle Current: How many droplets cross.
Entropy Current: How much "messiness" or heat crosses.

The Particle Result:
The number of droplets crossing behaves exactly as physicists expected. As they increase the pressure, more droplets flow. It's a smooth, predictable curve.

The Entropy Result (The Surprise):
The flow of "messiness" (entropy) does not behave smoothly. Instead, it oscillates (wiggles up and down) like a heartbeat as they increase the pressure.

Why? The paper explains this is a tug-of-war between two types of "dance moves":
- The "Reflection" Move: A dancer bounces back and forth within their own lake, carrying a lot of heat.
- The "Tunneling" Move: A dancer successfully crosses to the other lake, carrying less net heat.
As the pressure increases, these two moves turn on and off at different specific thresholds. When the "Reflection" move is strong, entropy goes up. When the "Tunneling" move takes over, entropy dips. This switching back and forth creates the wiggly, oscillating pattern.

4. The "Perfect" Bridge vs. The "Leaky" Bridge

The team tested the bridge at different levels of "transparency" (how easy it is to cross).

Low Transparency (A Leaky Bridge): The flow is weak, and the wiggles are small.
High Transparency (A Perfect, Ballistic Bridge): When the bridge is perfect, the wiggles in the entropy flow become very clear and pronounced. The scientists found that in this perfect state, the entropy flow is surprisingly small compared to what experiments with real cold gases have seen.

5. The Takeaway

The paper concludes that while their mathematical model (BCS theory) perfectly predicts how many particles move, it underestimates the entropy flow seen in real experiments.

This suggests that the real world is more complex than their "perfect dance floor" model. The real atoms might be doing things the model didn't account for, such as extra "fluctuations" or interactions that aren't part of the standard synchronized dance. The oscillating entropy is a signature of these complex quantum dance moves, but the fact that the model doesn't match the real data perfectly tells scientists they need to look for new physics beyond their current equations.

In short: They built a mathematical model of a super-fluid bridge, found that the "messiness" of the flow wiggles up and down in a complex pattern due to quantum dance moves, and realized that real-world experiments show even more chaos than their model predicted.

1. Problem Statement

The paper addresses the theoretical challenge of understanding entropy transport in ultra-cold neutral Fermi gases coupled via a Quantum Point Contact (QPC) under a chemical potential gradient. While particle transport in such systems is well-understood and matches experimental data, the behavior of entropy (and heat) transport, particularly in the superconducting (BCS) regime and under high-transparency (ballistic) conditions, remains theoretically underexplored.

Key specific problems addressed include:

The failure of linear dispersion approximations to capture heat and entropy currents due to inherent particle-hole symmetry.
The complex interplay of Multiple Andreev Reflections (MARs) in high-transparency junctions, where standard perturbative treatments fail.
The discrepancy between mean-field BCS predictions and experimental observations in unitary Fermi gases regarding entropy transport.

2. Methodology

The authors employ a rigorous Keldysh non-equilibrium field-theory approach combined with a tunneling-Hamiltonian formalism.

Model System: Two superfluid reservoirs (Left $c$ and Right $d$ ) described by a standard BCS mean-field model, connected by a single-point QPC. The system is driven out of equilibrium by a chemical potential difference ( $\Delta\mu$ ) and potentially a temperature difference.
Hamiltonian: The system is modeled using Nambu spinors to handle superconducting pairing. Crucially, the authors move beyond the linear dispersion approximation near the Fermi energy, adopting a quadratic dispersion relation ( $\epsilon_k = k^2/2m - E_F$ ). This is essential to break particle-hole symmetry and allow for non-zero heat/entropy currents.
Keldysh Formalism:
- The authors utilize the closed-time-path integral to compute non-equilibrium steady-state currents.
- They perform a Keldysh rotation to separate classical ( $\psi_{cl}$ ) and quantum ( $\psi_q$ ) fields.
- The particle and heat currents are derived from the time evolution of number and energy operators, expressed as correlation functions between the left and right reservoirs.
Numerical Implementation:
- The Green's functions for the coupled system are calculated in frequency-momentum space.
- Due to the non-diagonal nature of the action in frequency space (caused by MARs), the authors truncate the matrix representation after a finite number of "co-tunneling events" (up to $n=100$ for particle current and $n=50$ for entropy current).
- They analyze regimes of High Transparency ( $T \to 1$ , ballistic limit) and Low Transparency (perturbative regime).

3. Key Contributions

Derivation of Entropy Current in BCS Regime: The paper provides the first comprehensive derivation of the entropy current for a superfluid-superfluid (SS) junction using a microscopic BCS model with quadratic dispersion.
Identification of Oscillatory Behavior: The authors discover that the entropy current exhibits a non-monotonic, oscillatory behavior at low voltages in the ballistic limit, a feature absent in the particle current.
Mechanism of Oscillations: They attribute these oscillations to the competition between two distinct microscopic transport mechanisms:
- Andreev Composite Reflection (oMAR): An odd number of MARs where a particle is reflected as a hole within the same reservoir. This process yields a large additive heat contribution.
- Andreev Composite Tunneling (eMAR): An even number of MARs where a particle tunnels to the opposite reservoir. This process yields a subtractive heat contribution (difference between outgoing and incoming particle heat).
- As the bias voltage ( $\Delta\mu$ ) increases, different MAR channels (even vs. odd) are activated sequentially. The alternating dominance of these processes, which have opposite signs for their heat contributions, drives the oscillations.
Role of Dispersion: The study demonstrates that the quadratic dispersion is critical. In a linear dispersion model, particle-hole symmetry would cause the heat current to vanish. The quadratic term introduces the necessary asymmetry to generate the observed entropy transport.

4. Key Results

Particle Current: The calculated particle current in the ballistic limit ( $T=0.99$ ) reproduces the well-known non-linear behavior and MAR steps observed in previous literature and experiments, validating the model.
Entropy Current Oscillations:
- In the ballistic limit, the entropy current shows distinct oscillations within the superconducting gap ( $\Delta\mu < 2\Delta$ ).
- The oscillations correspond to the activation thresholds of MAR processes: $\Delta\mu = 2\Delta/n$ .
- The amplitude of these oscillations is significantly reduced as transparency decreases (low $T$ ), where the system becomes perturbative and MAR effects are suppressed.
Comparison with Unitary Regime:
- The authors compare their BCS mean-field results with recent experiments on unitary Fermi gases (strongly interacting).
- Discrepancy: While the BCS model accurately predicts the particle current, it underestimates the entropy current by approximately two orders of magnitude compared to experimental data in the unitary regime.
- Implication: This suggests that entropy transport is highly sensitive to correlation effects beyond mean-field theory (e.g., pairing fluctuations), which are significant in the unitary regime but neglected in the standard BCS approach.
SN Junctions: In Superconductor-Normal (SN) junctions, the entropy current shows a quadratic-like behavior inside the gap, distinct from the linear behavior outside, driven by single Andreev reflection.

5. Significance and Outlook

Theoretical Insight: The work clarifies the microscopic origin of entropy transport in superconductors, distinguishing between tunneling and reflection processes and highlighting the necessity of including particle-hole asymmetry (quadratic dispersion) for accurate thermodynamic predictions.
Experimental Relevance: The findings provide a benchmark for interpreting cold-atom experiments where entropy transport is being measured. The large deviation between the BCS prediction and unitary gas experiments highlights the need for theoretical models that go beyond mean-field approximations to capture strong correlation effects.
Future Directions: The authors suggest that investigating systems with unequal superconducting gaps ( $\Delta_c \neq \Delta_d$ ) or incorporating beyond-mean-field fluctuations could resolve the quantitative discrepancies with unitary gas experiments and further elucidate the role of pairing fluctuations in entropy transport.

In summary, this paper establishes a robust theoretical framework for entropy transport in superfluid QPCs, revealing complex oscillatory dynamics driven by the interplay of Andreev processes and identifying the limitations of mean-field theory in describing entropy in strongly interacting quantum systems.