SciPost logo

SciPost Submission Page

Variational solution of the superconducting Anderson impurity model and the band-edge singularity phenomena

by Teodor Iličin, Rok Žitko

This Submission thread is now published as

Submission summary

Authors (as registered SciPost users): Rok Žitko
Submission information
Preprint Link: https://arxiv.org/abs/2503.18902v2  (pdf)
Data repository: https://zenodo.org/records/15083680
Date accepted: June 24, 2025
Date submitted: June 19, 2025, 6:46 a.m.
Submitted by: Rok Žitko
Submitted to: SciPost Physics
Ontological classification
Academic field: Physics
Specialties:
  • Condensed Matter Physics - Theory
Approaches: Theoretical, Computational

Abstract

We propose a set of variational wavefunctions for the sub-gap spin-doublet and spin-singlet eigenstates of the particle-hole symmetric superconducting Anderson impurity model. The wavefunctions include up to two Bogoliubov quasiparticles in the continuum which is necessary to correctly capture the weak-coupling asymptotics in all parameter regimes. The eigenvalue problems reduce to solving transcendental equations. We investigate how the lowest singlet state evolves with increasing charge repulsion $U$, transitioning from a proximitized state (a superposition of empty and doubly occupied impurity orbitals, corresponding to an Andreev bound state) to a local moment that is Kondo screened by Bogoliubov quasiparticles (Yu-Shiba-Rusinov state). This change occurs for $U = 2\Delta$, where $\Delta$ is the BCS gap. At this point, the band-edge effects make the eigenenergy scale in a singular way as $\Gamma^{2/3}$, where $\Gamma$ is the hybridization strength. Away from this special point, regular $\Gamma$-linear behavior is recovered, but only for $\Gamma \lesssim (U/2-\Delta)^2/\Delta$. The singular behavior thus extends over a broad range of parameters, including those relevant for some quantum devices in current use. The singular state is an equal-superposition state with maximal fluctuations between the local impurity charge configurations. Accurately capturing the band-edge singularity requires a continuum model, and it cannot be correctly described by discrete (truncated) models such as the zero-bandwidth approximation or the superconducting atomic limit. We determine the region of parameter space where the second spin-singlet state exists: in addition to the whole $U<2\Delta$ ABS region, it also includes a small part of the $U>2\Delta$ YSR region for finite values of $\Gamma$, as long as some ABS wavefunction component is admixed.

Author indications on fulfilling journal expectations

  • Provide a novel and synergetic link between different research areas.
  • Open a new pathway in an existing or a new research direction, with clear potential for multi-pronged follow-up work
  • Detail a groundbreaking theoretical/experimental/computational discovery
  • Present a breakthrough on a previously-identified and long-standing research stumbling block

Author comments upon resubmission

The revised manuscript has improved presentation of result and the reference list has been updated with DOIs. The SciPost template is now being used.

List of changes

  • Figures improved for readability (color schemes).
  • Sec IIB made more compact by shifting the discussion of symmetries to the introduction and moving the lengthy expression for the energy expectation value to an Appendix.
  • Provenience of data (NRG or variational method) indicated in all captions, details about postprocessing to obtain the heat map explained.
  • Bibliography improvements (DOIs added, missing volume and article number added).
  • Conversion to SciPost template, formating and sizing of figures updated.

Published as SciPost Phys. 19, 006 (2025)

Login to report or comment