Experimentally, the concentration of quasiparticles in gapped superconductors always largely exceeds the equilibrium one at low temperatures. Since these quasiparticles are detrimental for many applications, it is important to understand theoretically the origin of the excess. We demonstrate in detail that the dynamics of quasiparticles localized at spatial fluctuations of the gap edge becomes exponentially slow. This gives rise to the observed excess in the presence of a vanishingly weak nonequilibrium agent.
Experimentally, the concentration of quasiparticles in gapped superconductors always largely exceeds the equilibrium one at low temperatures. Since these quasiparticles are detrimental for many applications, it is important to understand theoretically the origin of the excess. We demonstrate in detail that the dynamics of quasiparticles localized at spatial fluctuations of the gap edge becomes exponentially slow. This gives rise to the observed excess in the presence of a vanishingly weak nonequilibrium agent.
At low temperature, the concentration of quasiparticles observed in superconducting circuits far exceeds the predictions of microscopic BCS theory at equilibrium. As a source of dissipation, these excess quasiparticles degrade the performance of various devices. Therefore, understanding their dynamics, especially their recombination into Cooper pairs, is an active topic of current research. In disordered superconductors, spatial fluctuations in the superconducting gap can trap quasiparticles and modify their eigenspectrum. Since this spectrum plays a key role in quasiparticle dynamics, it must be carefully investigated. To this end, we introduce a toy model of a single trap. Specifically, we consider a shallow disk-shaped gap variation in a clean superconductor. Using a semiclassical approximation, we demonstrate the existence of multiple bound states and give the dependence of their number on the size and depth of the gap suppression. Extending our analysis beyond the semiclassical regime, in dimensions larger than one, we observe an infinite number of bound states very close to the gap edge, even for an arbitrarily small trap. These results deepen our understanding of trapped quasiparticles and may have important implications for their recombination in disordered superconductors.
An excess density of quasiparticles is widely observed in superconducting films. This excess causes performance degradation in a variety of superconducting devices, including decoherence in qubits. In this Letter, we evaluate the hypothesis of [1] that the quasiparticle excess is caused by anomalously slow recombination at low quasiparticle densities due to localization in sub-gap states. We probe the density of states in aluminum and niobium films using current-voltage measurements of tunnel junctions and extract upper bounds on the energy scales of the sub-gap states and gap smearing. With these parameters, we evaluate the recombination times predicted by [1] and find that slow recombination is not predicted to occur at observed quasiparticle densities in aluminumand niobium-based superconducting devices. These results suggest that the quasiparticle excess in ordered superconductors is primarily due to non-thermal sources of quasiparticle generation and not slow recombination.
Superconducting materials are useful for a wide range of applications, ranging from particle detectors to quantum computing. The presence of excitations (quasiparticles) in the superconductor can be a resource in some devices and a disturbance in others - in kinetic inductance detectors, the generation of quasiparticles by an incoming photon can be read out and used to calculate the energy of the photon, while the tunnelling of quasiparticles is a cause of decoherence in quantum computing. While early research mostly focused on superconductors close to the critical temperature, temperatures that are far below the critical temperatures are often used for current applications. At these low temperatures, the energy exchange between quasiparticles and phonons becomes slow, and even small disturbances may drive the energy distribution of the quasiparticles out of thermal equilibrium. In fact, experiments usually show a residual number of quasiparticles even at milikelvin temperatures, for which the thermal density of quasiparticles would be negligible. This thesis consists of two parts: In the first part, I model the energy distribution of quasiparticles in superconducting resonators: Quasiparticles absorb energy by absorbing photons and release it as phonons, which drives the superconductor out of thermal equilibrium. I solve a set of rate equations, both numerical and analytical, to calculate the energy distribution of the quasiparticles. For parameters typically used in kinetic inductance detectors made of aluminum, the microwave signal used to read out the resonator dominates the energy distribution of the quasiparticles, while their density is determined by other processes. At not too low temperatures, quasiparticles are mainly produced by thermal phonons, and my calculations reproduce recently measured quality factors of resonators at temperatures above 250 mK without free parameters. I show that even a very small number of photons of higher energy, which I ascribe to noise from the microwave generator, would dominate the quasiparticle density at lower temperatures. Under this assumption, my calculations also reproduce the quality factors at lower temperatures. In the second part of this thesis, I investigate the effects of spatial fluctuations of the superconducting energy gap on the dynamics of quasiparticles. Recently, it has been proposed that such fluctuations may lead to the localization of quasiparticles in states below the energy gap, which could cause a finite quasiparticle density even at very low temperatures. Using the Keldysh-Usadel formalism, I derive rates for the different processes (recombination/(de-) localization/generation) of quasiparticles in such systems. Using these rate equations I derive new upper limits for the quasiparticle density that could result from these localized states. These upper limits show that localization caused by weak magnetic impurities in typical aluminum systems does not lead to a relevant quasiparticle density. This contradicts previous theories, which only consider the recombination between two localized quasiparticles without examining their dynamics in detail.
Material inhomogeneities in a superconductor generically lead to broadening of the density of states and to subgap states. The latter are associated with spatial fluctuations of the gap in which quasiparticles can be trapped. Recombination between such localized quasiparticles is hindered by their spatial separation and hence their density could be higher than expectations based on the recombination between mobile quasiparticles. We show here that the recombination between localized and mobile excitations can be efficient at limiting the quasiparticle density. We comment on the significance of our findings for devices such as superconducting resonators and qubits. We find that for typical aluminum devices, the subgap states do not significantly influence the quasiparticle density.
Many superconducting devices rely on the finite gap in the excitation spectrum of a superconductor: thanks to this gap, at temperatures much smaller than the critical one the number of excitations (quasiparticles) that can impact the device’s behavior is exponentially small. Nevertheless, experiments at low temperature usually find a finite, non-negligible density of quasiparticles whose origin has been attributed to various non-equilibrium phenomena. Here, we investigate the role of photons with energy exceeding the pair-breaking threshold 2\Delta 2 Δ as a possible source for these quasiparticles in superconducting resonators. Modeling the interacting system of quasiparticles, phonons, sub-gap and pair-breaking photons using a kinetic equation approach, we find analytical expressions for the quasiparticles’ density and their energy distribution. Applying our theory to measurements of quality factor as function of temperature and for various readout powers, we find they could be explained by assuming a small number of photons above the pair-breaking threshold. We also show that frequency shift data can give evidence of quasiparticle heating.
In the superconducting state, the presence of a finite gap in the excitation spectrum implies that the number of excitations (quasiparticles) is exponentially small at temperatures well below the critical temperature. Conversely, minute perturbations can significantly impact both the distribution in energy and the number of quasiparticles. Typically, the interaction with the electromagnetic environment is the main perturbation source driving quasiparticles out of thermal equilibrium, while a phonon bath is responsible for restoration of equilibrium. Here we derive approximate analytical solutions for the quasiparticle distribution function in superconducting resonators and explore the impact of nonequilibrium on two measurable quantities: the resonator's quality factor and its resonance frequency. Applying our results to experimental data, we conclude that while at intermediate temperatures there is clear evidence for the nonequilibrium effects due to heating of the quasiparticles by photons, the low-temperature measurements are not explained by this mechanism.