Thermodynamic Constraints on Perfect Equilibrium Superconducting Diodes

Technical Summary: Thermodynamic Constraints on Perfect Equilibrium Superconducting Diodes

Problem Statement
Superconducting diodes (SDs), devices exhibiting unequal critical current magnitudes ($|I_c^+| \neq |I_c^-|$) in opposite directions, offer potential for dissipationless rectification and nonreciprocal signal processing. While the SD effect has been observed across diverse platforms—from Van der Waals heterostructures to asymmetric Josephson junctions (JJs)—achieving high efficiency (where the diode asymmetry $\epsilon = |I_c^-/I_c^+| \ll 1$) remains a challenge. In unconstrained bulk superconductors without engineered junctions, efficiencies are typically modest ($\eta \sim 10^{-2}$), whereas Josephson platforms often report higher efficiencies, sometimes approaching $\eta \approx 1$. The central question addressed is whether a "perfect" equilibrium superconducting diode ($\epsilon = 0$) is theoretically possible and what thermodynamic principles govern the limits of diode efficiency across different platforms.

Methodology and Theoretical Framework
The author employs a general thermodynamic approach based on the Gibbs free energy ($G$) of quasi-1D superconducting systems. The transport is governed by a global phase variable $\Theta$, while internal condensate degrees of freedom (e.g., relative phases in multiband systems) relax at fixed $\Theta$. The effective free energy $F(\Theta)$ is defined by minimizing the Helmholtz free energy functional over these internal variables.

Key theoretical tools include:

1. Continuity and Analyticity: The analysis assumes conventional local solids where the free energy $F(\Theta)$ must be continuous. It investigates whether $F(\Theta)$ can be a single smooth analytic branch or if it requires singularities (non-analyticities).
2. Critical Current Definitions: Critical currents are defined as the extrema of the supercurrent $I(\Theta) = 2e \partial_\Theta F(\Theta)$. The diode asymmetry is $\epsilon = |I_c^-/I_c^+|$.
3. Mathematical Bounds: For systems described by single smooth branches, the author utilizes polynomial approximations (Landau theory) and Fourier series (for Josephson junctions). The Remez inequality and properties of Chebyshev polynomials are used to derive lower bounds on $\epsilon$ for finite-order expansions.
4. Exactly Solvable Model: An Ising-coupled superconductor model is constructed to explicitly demonstrate how piecewise free energy landscapes can realize near-perfect diode behavior.

Key Contributions and Results

• Impossibility of Exact Perfect Diodes ($\epsilon = 0$): The paper proves that an exact perfect equilibrium diode ($\epsilon = 0$) is forbidden in conventional local systems. If $I_c^- = 0$ and $I_c^+ \neq 0$, the supercurrent $I(\Theta)$ would be non-negative everywhere within the superconducting interval. Integrating this positive current would result in a net change in free energy $F(\Theta_c^+) - F(\Theta_c^-) > 0$, contradicting the requirement that the free energy must be continuous and equal at the boundaries of the superconducting domain (Eq. 4). Thus, $\epsilon = 0$ would require a discontinuity in the free energy, which is excluded in conventional thermodynamics.

• Necessity of Singularities for Near-Perfect Behavior: While $\epsilon = 0$ is forbidden, the limit $\epsilon \to 0$ is thermodynamically allowed but requires specific conditions. The author demonstrates that achieving $\epsilon \ll 1$ necessitates singularities (non-analyticities) in the free energy landscape $F(\Theta)$. Specifically, the system must undergo a "softening" of the unfavorable current-carrying branch. This implies a piecewise structure where the free energy consists of broad, nearly flat regions (where $I(\Theta) < 0$ is small) alongside regions with finite slopes. These transitions correspond to branch-switching between competing superconducting states or the formation of domain walls.

• Bounds for Smooth Single-Branch Systems: For superconductors described by a single smooth analytic branch (e.g., standard Landau theory with a finite-degree polynomial), the paper derives a rigorous lower bound on $\epsilon$.

  • In Landau theory (polynomial of degree $n$), $\epsilon \geq 1/T_n(\frac{1+r}{1-r})$, where $T_n$ is the Chebyshev polynomial and $r$ is the fraction of the phase window where current is positive.
  • This implies that achieving a small $\epsilon$ requires a polynomial of parametrically large degree ($n \sim \log(1/\epsilon)$), corresponding to a complex landscape with many competing superconducting phases.
  • Similarly, for Josephson junctions, a finite truncation of the current-phase relation (finite number of harmonics) imposes a lower bound on $\epsilon$.

• Ising-Coupled Model Demonstration: The paper explicitly constructs an exactly solvable Ising-coupled superconductor model. By coupling the superconductor to an Ising variable $\sigma = \pm 1$, the free energy becomes piecewise. The model shows that by tuning parameters (stiffness $K$ and magnetoelectric coupling $\gamma$) such that one branch becomes extremely soft (nearly $q$-independent) while the other remains stiff, the system can achieve $\epsilon \ll 1$. This confirms that non-analyticities (cusps in $F(\Theta)$) enable the necessary softening of the negative current branch.

• Reconciliation of Experimental Observations: The framework explains why high efficiencies are common in Josephson platforms but rare in intrinsic bulk superconductors. Josephson devices often operate under:

  1. Phase constraints: Embedded in larger environments preventing global minimization.
  2. Metastability: Switching occurs from local minima rather than the global equilibrium state.
  3. External drives: Nonequilibrium protocols.
     These conditions allow them to evade the strict unconstrained equilibrium limits derived in the paper. Conversely, high efficiency in nominally intrinsic platforms (like twisted trilayer graphene) suggests the presence of hidden internal structures (competing orders) or metastability that creates the required non-analytic free energy landscape.

Significance and Claims
The paper provides a unified thermodynamic framework for interpreting superconducting diode efficiencies. Its primary claim is that exact perfect equilibrium diodes are forbidden by the continuity of the free energy, and near-perfect behavior in unconstrained equilibrium systems requires singular restructuring of the superconducting landscape (e.g., branch switching).

The work serves as a "no-go theorem" for simple, smooth superconductors, explaining why $\epsilon = O(1)$ is the norm in bulk materials. It posits that unusually small $\epsilon$ in equilibrium platforms is a diagnostic signature of nearby internal structure, such as competing superconducting orders or metastable states. The results offer a design principle: to achieve high efficiency in equilibrium, one must engineer systems with specific non-analytic free energy features, while high efficiencies in Josephson devices are often a consequence of operating outside strict global equilibrium constraints.