Quasiclassical theory of vortex states in locally non-centrosymmetric superconductors: application to CeRh$_{2}$As$_{2}$

Basic Information

• Paper ID: 2510.14452
• Title: Quasiclassical theory of vortex states in locally non-centrosymmetric superconductors: application to CeRh$_{2}$As$_{2}$
• Authors: Akihiro Minamide, Youichi Yanase (Department of Physics, Kyoto University)
• Classification: cond-mat.supr-con (Condensed Matter Physics - Superconductivity)
• Publication Date: October 17, 2025
• Paper Link: https://arxiv.org/abs/2510.14452

Abstract

CeRh$_{2}$As$_{2}$ is a heavy fermion superconductor discovered in 2021 that exhibits two distinct superconducting phases under magnetic field applied along the c-axis. This unconventional phase diagram is attributed to the local inversion symmetry breaking at Ce sites. At low magnetic fields, a conventional even-parity spin-singlet superconducting state is realized, while at high magnetic fields, an odd-parity spin-singlet superconducting state becomes stable, where the order parameter exhibits sign alternation between adjacent Ce layers. This study employs a quasiclassical approach to investigate vortex states in bilayer superconductors with locally broken inversion symmetry. We calculate the local density of states (LDOS) in the vortex lattice state and find that the pairing symmetry of different superconducting states is clearly manifested in the LDOS peak structure at the vortex core. Since LDOS is experimentally observable, our work provides a pathway for experimental verification of the superconducting parity transition in CeRh$_{2}$As$_{2}$.

Research Background and Motivation

Core Problem

This study addresses the microscopic mechanism of multiple superconducting phases in CeRh$_{2}$As$_{2}$, particularly the properties of vortex states and methods for their experimental identification.

Research Significance

1. Novel Physical Phenomena: CeRh$_{2}$As$_{2}$ is the first material observed to exhibit a superconducting parity transition, possessing important fundamental physics significance
2. Theoretical Challenge: The theory of vortex states in locally non-centrosymmetric systems remains incomplete, requiring development of new theoretical frameworks
3. Experimental Verification: Direct experimental methods for identifying the order parameter structure of different superconducting phases are lacking

Limitations of Existing Methods

1. Uniform State Assumption: Previous theories primarily focused on spatially uniform states, neglecting vortex degrees of freedom
2. Phenomenological Approach: Ginzburg-Landau theory lacks microscopic foundation and is only valid near the transition temperature
3. Single Vortex Approximation: Neglects inter-vortex interactions and thermodynamic stability

Research Motivation

Establish a microscopic model-based quasiclassical theory of vortex states, calculate experimentally observable physical quantities, and provide experimental verification pathways for the superconducting parity transition in CeRh$_{2}$As$_{2}$.

Core Contributions

1. Theoretical Framework: Developed a multi-band quasiclassical theory of vortex states in locally non-centrosymmetric superconductors
2. Phase Diagram Calculation: Determined the complete H-T phase diagram including orbital effects, extended to low temperature regime
3. Vortex Properties: Revealed different behaviors of vortex core radius between BCS and PDW states
4. LDOS Calculation: Predicted characteristic LDOS spectra at vortex cores for different superconducting states
5. Experimental Scheme: Proposed experimental approaches using STM-STS and other techniques to verify the superconducting parity transition
6. Novel State Prediction: Analyzed LDOS characteristics of superconducting meron states

Methodology Details

Task Definition

Study vortex lattice states in bilayer Rashba model-described locally non-centrosymmetric superconductors under magnetic field, calculate observable quantities such as local density of states, and identify characteristics of different superconducting phases.

Model Architecture

1. Bilayer Rashba Model

The Hamiltonian is:

H = (1/2)∫d²r₁d²r₂ C⃗†(r₁)Ȟ_BdG(r₁,r₂)C⃗(r₂)

where the normal state Hamiltonian is:

Ĥ^N(k) = ξ(k)σ̃₀⊗τ̃₀ + αg(k)·σ̃⊗τ̃_z + t_⊥σ̃₀⊗τ̃_x - h·σ̃⊗τ̃₀

Key parameters:

• α: Antisymmetric spin-orbit coupling strength
• t_⊥: Interlayer hopping
• g(k) = k×ẑ/k_F: Staggered Rashba-type spin-orbit coupling

2. Quasiclassical Theory

Through Wigner transformation and energy integration, the effective low-energy Eilenberger equation is obtained:

2ωₙf̃^(ν) + i(-1)^ν h cosχ[υ̃_z, f̃^(ν)] + iv_F·Πf̃^(ν) = Δ̃^(ν)g̃^(ν) + g̃^(ν)Δ̃^(ν)

where χ is defined by e^(iχ) = (t_⊥ + iα⟨|g|²⟩^(1/2))/√(t_⊥² + α²⟨|g|²⟩).

3. Landau Level Expansion

The order parameter is expanded in Landau levels:

Δⱼ(k,r) = w(k)∑_N dⱼ,ₙΨₙ(r)  (j = e,o)

where Ψₙ(r) are eigenfunctions of the triangular vortex lattice.

Technical Innovations

1. Multi-band Treatment: Accounts for contributions from two quasi-degenerate bands, handles inter-band and intra-band pairing
2. Adachi Approximation: Employs improved Pesch approximation for enhanced computational accuracy
3. Self-consistent Solution: Develops self-consistent numerical methods for vortex lattice states
4. Analytical Expressions: Obtains analytical solutions for Green's functions in pure BCS and PDW states

Experimental Setup

Model Parameters

• α/t_⊥ = 15.0: Ratio of spin-orbit coupling to interlayer hopping
• αₘ = 4.8: Maki parameter
• T^o_c0/T^e_c0 = 0.9: Ratio of transition temperatures for two channels

Computational Details

• Spatial grid: 41×41 discrete points
• Landau level cutoff: Up to 12th order
• Matsubara frequency cutoff: εc = 40T^e_c0
• Broadening parameter: η = 0.03T^e_c0

Numerical Methods

• Broyden method for self-consistent iteration
• Triangular vortex lattice assumption
• Extreme type-II superconductor approximation (κ ≫ 1)

Experimental Results

Main Results

1. Phase Diagram Determination

• Low-field BCS state: Sublattice symmetric order parameter (Δₑ,Δₒ) = (Δ,0)
• High-field PDW state: Sublattice antisymmetric order parameter (Δₑ,Δₒ) = (0,Δ)
• Parity transition line: Located where free energies of both phases are equal in the intermediate field region

2. Vortex Core Radius

Key findings:

• In the BCS state, core radius reaches minimum at intermediate field then increases again
• In the PDW state, core radius decreases monotonically
• Core radius undergoes sudden contraction at the parity transition

Physical mechanism: Paramagnetic pair-breaking effect in the BCS state leads to reduced effective superconducting gap and increased coherence length.

3. Local Density of States Characteristics

BCS State:

• Double-peak structure with up-down splitting appears at vortex core
• Energy splitting corresponds to effective Zeeman field h_eff

PDW State:

• Large zero-energy peak at vortex core
• Pseudospin-singlet pairing completely suppresses paramagnetic pair-breaking effect

4. Superconducting Meron State

• Interlayer vortex positions arranged in staggered pattern
• LDOS exhibits ridge-like structure connecting vortex cores of different layers
• Reflects quasiparticle transfer induced by interlayer coupling

Quantitative Results

• Parity transition field: μ₀H* ≈ 4T (consistent with experiment)
• Upper critical field ratio: H^o_c2/H^e_c2 ≈ 2-3
• Vortex core radius: ξ_vc ~ 10-50 R₀ (R₀ = v_F/T^e_c0)

Related Work

Theoretical Background

1. Locally Non-centrosymmetric Theory: Basic framework established by Fischer, Sigrist, and others
2. CeRh₂As₂ Theory: Uniform state theory by Yoshida, Yanase, and colleagues
3. Vortex State Theory: Single vortex studies by Higashi and others, GL theory by Möckli and colleagues

Advantages of This Work

1. Microscopic Foundation: Based on microscopic Hamiltonian, avoiding phenomenological parameters
2. Complete Treatment: Considers vortex lattice effects and thermodynamic stability
3. Low-Temperature Extension: Overcomes temperature limitations of GL theory
4. Experimental Connection: Calculates observable quantities, provides experimental verification schemes

Conclusions and Discussion

Main Conclusions

1. Confirms the microscopic mechanism of BCS-PDW parity transition in CeRh₂As₂
2. Vortex core radius and LDOS spectrum can serve as experimental probes for identifying superconducting phases
3. Predicts novel superconducting meron states and their characteristics

Limitations

1. s-wave Approximation: Neglects anisotropy of d-wave pairing
2. Weak Coupling: Does not consider strong correlation effects
3. Other Orders: Ignores coexisting magnetic and quadrupole orders
4. Temperature Limitation: Convergence difficulties in low-temperature regime

Future Directions

1. Consider anisotropic Fermi surface and d-wave pairing
2. Include strong correlation effects
3. Study coupling with other orders
4. Investigate vortex state properties under high pressure

In-Depth Evaluation

Strengths

1. Theoretical Rigor: Complete quasiclassical theory based on microscopic model
2. Computational Accuracy: Reliable numerical methods with good convergence
3. Physical Clarity: Reveals microscopic mechanisms of vortex states
4. Experimental Guidance: Proposes clear experimental verification schemes
5. Predictive Power: Predicts novel physical phenomena

Weaknesses

1. Model Simplification: s-wave isotropic assumption is overly simplified
2. Parameter Dependence: Results are sensitive to model parameters
3. Experimental Comparison: Lacks direct comparison with existing experimental data
4. Universality: Primarily focused on CeRh₂As₂, limited generalizability

Impact

1. Theoretical Contribution: Provides complete vortex state theory for locally non-centrosymmetric superconductors
2. Experimental Advancement: Provides theoretical guidance for STM-STS and other experimental techniques
3. Field Development: Promotes research on materials with multiple superconducting phases
4. Methodological Value: Quasiclassical methods can be extended to other complex superconductors

Applicable Scenarios

1. Materials Research: CeRh₂As₂ and similar locally non-centrosymmetric superconductors
2. Experimental Design: STM-STS, SANS, NMR and other vortex state measurements
3. Theoretical Development: Vortex state theory for multi-component superconductors
4. Device Applications: Quantum device design based on vortex states

References

The paper cites 119 references, primarily including:

• Experimental discovery and characterization of CeRh₂As₂ 1-4
• Locally non-centrosymmetric superconductor theory 9-15
• Quasiclassical theoretical methods 52-53, 57-61
• Vortex state-related theory 49-50, 65-77

Overall Assessment: This is a high-quality theoretical physics paper that provides important theoretical support for understanding multiple superconducting phases in CeRh₂As₂. Although some model simplifications exist, the predicted experimental phenomena have important verification value and contribute positively to advancing research in this field.