Quantum fluctuation-induced first-order breaking of time-reversal symmetry in unconventional superconductors

Basic Information

• Paper ID: 2510.09234
• Title: Quantum fluctuation-induced first-order breaking of time-reversal symmetry in unconventional superconductors
• Author: Yin Shi (Beijing National Laboratory for Condensed Matter Physics and Institute of Physics, Chinese Academy of Sciences)
• Classification: cond-mat.supr-con (Condensed Matter Physics - Superconductivity)
• Publication Date: October 13, 2025
• Paper Link: https://arxiv.org/abs/2510.09234

Abstract

This paper investigates the quantum phase transition of spontaneous time-reversal symmetry breaking in unconventional superconductors. In superconductors with competing non-degenerate pairing channels, such phase transitions can generate robust topological superconductivity and anomalous magnetism. Using the hole-doped square lattice t-J model, the authors derive the free energy with phase fluctuation corrections and find that quantum phase fluctuations cause the time-reversal asymmetric s+id phase region to split into a dome structure. A first-order boundary exists between this dome and the d phase, suggesting possible phase separation between time-reversal symmetric and asymmetric phases.

Research Background and Motivation

Problem Definition

1. Core Issue: Current understanding of spontaneously time-reversal symmetry-breaking superconducting states remains at the mean-field theory level, neglecting the important effects of order parameter phase fluctuations
2. Significance: Superconducting states with broken time-reversal symmetry themselves may possess topological non-triviality, offering potential for robust quantum computing protected by quantum decoherence
3. Existing Limitations:
   • Previous theoretical treatments are based on mean-field approximation
   • Neglects dynamical order parameter phase fluctuations, which are particularly important in two-dimensional systems
   • Lacks deep understanding of how quantum phase fluctuations affect phase transition properties

Research Motivation

• Time-reversal breaking superconducting phases are characterized by non-zero relative phases between two order parameters, making them particularly susceptible to order parameter phase fluctuations
• Quantum order parameter phase fluctuations have significant effects on the superconducting ground state of quasi-two-dimensional superconductors
• Need to establish a theoretical framework incorporating long-range order parameter phase fluctuations

Core Contributions

1. Theoretical Framework Innovation: Establishes a self-consistent theoretical framework incorporating long-range order parameter phase fluctuations coupled with charge density fluctuations
2. Phase Diagram Correction: Discovers that quantum phase fluctuations cause the s+id phase region to split into a dome structure, forming a first-order boundary with the d phase
3. Phase Transition Property Modification: Demonstrates that quantum fluctuations can transform a second-order phase transition into a first-order one
4. Experimental Guidance: Provides theoretical basis for understanding time-reversal symmetry breaking in twisted cuprate oxide Josephson junctions

Methodology Details

Task Definition

Study the quantum phase transition between competing s-wave and d-wave superconducting pairing channels in the square lattice t-J model, with particular focus on how order parameter phase fluctuations affect the time-reversal symmetry breaking phase transition.

Model Architecture

Hamiltonian Construction

Employs the t-J model with long-range Coulomb interactions: $H = -\sum_{rr'\sigma} \tilde{t}_{r-r'}c^\dagger_{r\sigma}c_{r'\sigma} - \frac{J}{2}\sum_{\langle rr'\rangle} b^\dagger_{rr'}b_{rr'} + \sum_q \frac{V_q}{2N}\rho^\dagger_q\rho_q$

where:

• $c_{r\sigma}$: electron annihilation operator with spin σ
• $b_{rr'} = c_{r\downarrow}c_{r'\uparrow} - c_{r\uparrow}c_{r'\downarrow}$: singlet pair annihilation operator
• $V_q = e^2/(2\epsilon_b\epsilon_0 aq)$: two-dimensional Coulomb potential

Path Integral Formulation

After decoupling interactions via Hubbard-Stratonovich transformation, the action becomes: $S = \int_0^{1/T} d\tau \left\{\sum_{rr'} \psi^\dagger_r(\tau)[\delta_{rr'}(\partial_\tau + \phi_r(\tau)\tau_3) - \tilde{t}_{r-r'}\tau_3 + \Delta_{rr'}(\tau)\tau_+ + \Delta^*_{rr'}(\tau)\tau_-]\psi_{r'}(\tau) + \frac{2}{J}\sum_{\langle rr'\rangle}|\Delta_{rr'}(\tau)|^2 - \sum_q \frac{|\phi_q(\tau)|^2}{2NV_q}\right\}$

Phase Fluctuation Treatment

Employs gauge transformation $\psi_r(\tau) \to \psi_r(\tau)e^{i\theta_r(\tau)\tau_3}$ and $\Delta_{rr'}(\tau) \to \Delta_{r-r'}e^{i[\theta_r(\tau)+\theta_{r'}(\tau)]}$, expanding the action in terms of small spacetime gradients of $\theta_r(\tau)$ and $\phi_r(\tau)$ up to quadratic order.

Technical Innovations

1. Self-Consistent Treatment: Unlike previous theories, this method achieves complete self-consistent determination of order parameters by minimizing the free energy including phase fluctuation corrections
2. Long-Range Fluctuation Coupling: Accounts for coupling between long-range order parameter phase fluctuations and charge density fluctuations
3. Beyond Gaussian Approximation: Incorporates higher-order effects of phase fluctuations through self-consistent methods

Free Energy Expression

The resulting free energy per unit cell is: $F(\Delta_x,\Delta_y) = -\frac{1}{N}\sum_k [2T\ln(1+e^{-E_k/T}) + E_k] + \frac{2}{J}(|\Delta_x|^2 + |\Delta_y|^2) + \frac{1}{N}\sum_{q<q_c}[2T\ln(1-e^{-\Omega_q/T}) + \Omega_q]$

where the Nambu-Goldstone mode spectrum is: $\Omega_q^2 = \frac{\chi V_q + 1}{\chi}\sum_{\alpha,\beta=x}^z q_\alpha q_\beta D_{\alpha\beta}$

Experimental Setup

Model Parameters

• Nearest-neighbor hopping: $t \approx 0.25$ eV
• Next-nearest-neighbor hopping: $-0.2t$
• Next-next-nearest-neighbor hopping: $0.1t$
• Exchange interaction: $J = 1.2t$
• Lattice constant: $a \approx 5.4$ Å
• Background dielectric constant: $\epsilon_b \approx 4.5$

Order Parameter Parameterization

Utilizing square lattice symmetry, the order parameters are rewritten as:

• $\Delta_x = (\Delta_s + i\Delta_d)/2$
• $\Delta_y = (\Delta_s - i\Delta_d)/2$

where $\Delta_s$ and $\Delta_d$ correspond to real superconducting gaps in the s-wave and d-wave channels, respectively.

Experimental Results

Main Findings

Phase Diagram Evolution

1. Without Phase Fluctuations: The s-d overlapping superconducting region forms an s+id phase, with both phase transitions to pure phases being second-order
2. Including Phase Fluctuations:
   • The s+id phase region splits into a small dome structure near the quantum critical point
   • The dome forms a first-order phase boundary with the d phase
   • The s+id phase region shrinks significantly (36% reduction at T=0.002t)

First-Order Phase Transition Evidence

At low temperatures ($T < 0.002t$), the superconducting gaps $\Delta_s$ and $\Delta_d$ exhibit discontinuous jumps, whereas mean-field results remain continuous. The free energy landscape reveals two minima:

• Lower minimum: stable s+id phase (finite $\Delta_s$ and $\Delta_d$)
• Higher minimum: metastable d phase (zero $\Delta_s$, finite $\Delta_d$)

Gap Modification Effects

Phase fluctuation effects on superconducting gaps:

• Generally suppress both $\Delta_s$ and $\Delta_d$
• Inside the dome: $\Delta_d$ is suppressed while $\Delta_s$ is enhanced (maximum enhancement of 31% at zero temperature)
• Bring the $\Delta_s/\Delta_d$ ratio closer to unity, counterintuitively enhancing time-reversal breaking characteristics

Phase Stiffness Suppression

Phase fluctuations suppress phase stiffness even at zero temperature:

• Maximum suppression of 15% near the quantum critical point
• Negligible suppression in regions far from the quantum critical point
• Demonstrates the important role of the time-reversal breaking quantum critical point in suppressing ground state phase stiffness

Related Work

Time-Reversal Breaking Superconductivity

1. Degenerate Channels: Such as $p_x + ip_y$ states in square lattices and $d_{x^2-y^2} + id_{xy}$ states in triangular lattices
2. Competing Non-Degenerate Channels: Magnetic field-induced phase transitions between $d_{x^2-y^2}$ and $d_{xy}$ channels in cuprates
3. Twisted Josephson Junctions: Topological $d_{x^2-y^2} + e^{i\vartheta}d_{xy}$ states in twisted cuprate bilayers

Phase Fluctuation Theory

• Emery-Kivelson mechanism: importance of phase fluctuations for superconductors with small superfluid density
• Significant suppression of zero-temperature superfluid density by order parameter phase fluctuations in quasi-two-dimensional systems

Conclusions and Discussion

Main Conclusions

1. Phase Transition Property Modification: Quantum phase fluctuations transform the time-reversal breaking quantum phase transition from second-order to first-order
2. Phase Separation Possibility: The first-order transition suggests potential spontaneous phase separation between d and s+id phases in certain doping regions
3. Phase Region Shrinkage: Phase fluctuations significantly reduce the stable region of the time-reversal breaking phase

Experimental Significance

Provides a theoretical framework for understanding recently observed phenomena in twisted cuprate oxide Josephson junction experiments:

• Explains the sharp disappearance of the Josephson diode effect at certain twist angles
• Supports the view that phase fluctuation-induced first-order transitions and phase region splitting occur

Limitations

1. Model Simplification: Based on the t-J model, may not fully capture the complexity of real materials
2. Approximate Treatment: Employs self-consistent Gaussian approximation, neglecting higher-order fluctuation effects
3. Parameter Sensitivity: The choice of cutoff wave vector $q_c$ has certain influence on results

Future Directions

1. Multi-Band Systems: Extend the theory to complex systems such as twisted cuprate oxide Josephson junctions
2. Experimental Verification: Seek direct experimental evidence for phase separation and first-order transitions
3. Topological Properties: Investigate the effects of phase fluctuations on topological superconducting properties

In-Depth Evaluation

Strengths

1. Theoretical Innovation: First systematic study of phase fluctuation effects on time-reversal breaking superconducting phase transitions
2. Advanced Methodology: Establishes a self-consistent phase fluctuation theoretical framework, surpassing previous mean-field treatments
3. Physical Insights: Reveals that quantum fluctuations can fundamentally alter the nature and order of phase transitions
4. Experimental Relevance: Provides theoretical basis for understanding current hot experimental phenomena

Weaknesses

1. Model Limitations: Only considers the t-J model; applicability to more complex systems requires verification
2. Approximation Level: Although self-consistent, still based on Gaussian approximation, potentially missing important non-Gaussian effects
3. Quantitative Predictions: Lacks quantitative comparison with specific experiments

Impact

1. Academic Value: Provides new theoretical tools for understanding phase fluctuation effects in quantum many-body systems
2. Application Prospects: Offers guidance for designing and understanding topological superconducting devices
3. Method Generalization: The developed theoretical methods can be extended to other strongly correlated electron systems

Applicable Scenarios

1. Quasi-Two-Dimensional Superconductors: Particularly systems with competing order parameters
2. Artificial Superconducting Structures: Such as twisted bilayer systems
3. Topological Superconductivity Research: Design of time-reversal breaking topological superconducting states

References

The paper cites 38 important references covering classical and recent works in time-reversal breaking superconductivity, phase fluctuation theory, t-J model research, and related fields, providing a solid theoretical foundation for the research.

This paper holds significant importance in theoretical condensed matter physics, advancing understanding of quantum phase transitions in unconventional superconductors while providing valuable theoretical explanations for related experimental phenomena. The self-consistent phase fluctuation theoretical method developed has broad application prospects.