我们为具有粒子-空穴对称性的超导Anderson杂质模型的亚阈值自旋双重态和自旋单重态本征态提出了一组变分波函数。这些波函数包括连续区域内最多两个Bogoliubov准粒子，这对于在所有参数范围内正确捕捉弱耦合渐近行为是必要的。本征值问题归结为求解超越方程。我们研究了最低的自旋单重态随电荷排斥力$U$增加时的演化过程，从一个近邻化状态（空态和双占据杂质轨道的叠加，对应于Andreev束缚态）过渡到由Bogoliubov准粒子Kondo屏蔽的局域磁矩（Yu-Shiba-Rusinov态）。这一转变发生在$U=2\Delta$处，其中$\Delta$为BCS能隙。在此点处，带边效应使本征能量尺度以奇异的方式随着$\Gamma^{2/3}$变化，其中$\Gamma$为杂化强度。远离此特殊点时，在$\Gamma \lesssim (U/2-\Delta)^2/\Delta$范围内可恢复常规的$\Gamma$-线性行为。因此，这种奇异行为扩展到了广泛的参数范围，其中包括一些当前使用的量子器件相关的参数。奇异态是一种等幅叠加态，其局部杂质电荷构型之间的涨落达到最大值。准确捕捉带边奇异现象需要连续模型，而离散（截断）模型如零带宽近似或超导原子极限无法正确描述该现象。我们确定了第二自旋单重态存在的参数空间区域：除了整个$U<2\Delta$的ABS区域外，还包括有限$\Gamma$值下$U>2\Delta$的YSR区域的一部分，只要存在一定的ABS成分混合。

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 component is admixed.