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.
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.
Low temperature microwave superconducting resonators (SRs) are attractive candidates for producing quantum-sensitive, arrayable energy or power detectors for astrophysical and other precision measurement applications. Their readout uses a microwave probe signal with quanta of energy well below the threshold for pair-breaking in the superconductor. We have calculated the non-equilibrium quasiparticle and phonon distributions generated by the photons of the probe signal of a resonator operating well below its superconducting transition temperature T c as the absorbed probe power was changed using the coupled kinetic equations described by Chang and Scalapino. The calculations give insight into a rate equation estimate which suggests that the quasiparticle distributions can be driven far from their thermal equilibrium value for typical readout powers. From the driven quasiparticle distribution functions, the driven quasiparticle number densities and lifetimes were calculated. An effective temperature to describe the driven quasiparticles was defined. The non-equilibrium lifetimes were compared to the distribution-averaged thermal lifetimes at the effective temperature and good agreement was found typically within a few per cent. We used the non-equilibrium quasiparticle distribution to model a representative SR. The complex conductivity and hence the frequency dependence of the experimentally measured forward scattering parameter S 21 of the SR as a function of absorbed power were found. The non-equilibrium S 21 cannot be accurately modeled by a thermal distribution even at its own elevated temperature, having a higher quality factor in all cases studied, although for low absorbed powers the two effective temperatures are similar. From the non-equilibrium lifetimes and number densities we determined the achievable noise equivalent power (NEP) of the resonator used as a power detector as a function of absorbed microwave power. Simpler expressions to evaluate the effective quasiparticle temperature as a function of absorbed power have also been derived. We conclude that multiple photon absorption from the microwave probe increases the quasiparticle number above the thermal background and ultimately limits the achievable NEP of the resonator at temperatures well below T c .
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 this thesis I describe work on developing theoretical and numerical models of supercon- ducting thin-film microwave resonators. Superconducting resonators are used in a variety of applications, one of which is as kinetic inductance detectors (KIDs). KIDs are ultra-low noise, highly sensitive, multiplexable detectors, with uses in a wide variety of fields including astrophysics, medical imaging, and particle physics. Resonators are also crucial for supercon- ducting qubit readout, superconducting mixers and parametric amplifiers, and as multiplexers for other devices. The results described in this thesis apply to all thin film resonators, but are primarily described in the context of KIDs. I develop models of how constant absorbed power affects superconducting thin films, driving them out of thermal equilibrium with the bath. Using the Chang & Scalapino equations, I calculate numerically the steady-state quasiparticle and phonon distributions for various sources of power (sub-gap and above-gap photons and phonons, single frequency and broadband). Many new results emerge, a few of which are: the quasiparticle heating effects of microwave power, explaining the experimentally measured saturation of resonator performance; the frequency dependence of the quasiparticle generation efficiency in the sub-mm wavelength range; and the importance of phonon trapping in thin films at low temperatures. I then use these nonequilibrium results in higher-level device models. Starting from a generic framework for resonators, I present and implement a model of a KID with variable thermal isolation. I calculate the effects of quasiparticle heating on both the large-signal and small-signal behaviour, highlighting the effects of electrothermal feedback caused by readout power. Electrothermal feedback is shown to be able to increase or decrease the magnitude and bandwidth of the responsivity and noise, depending on the operating point. Finally, I propose an experimental measurement of quasiparticle heating effects in a new device – a four-port ring resonator. In such a device, unlike in the standard two-port resonator, the device can be heated and read out independently. I develop detailed models for the thermal and electrical behaviour of these devices, and suggest a scheme by which the quasiparticle heating effects of microwave power, predicted in this thesis, can be measured.
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.
Nonequilibrium quasiparticle and phonon dynamics are central to the operation of superconducting devices. Superconducting detectors, such as superconducting nanowire single-photon detectors, transition-edge sensors, or microwave kinetic inductance detectors, perform best when a large quasiparticle population is generated in response to small perturbations. Conversely, for superconducting qubits and topologically protected Majorana fermions, even relatively small quasiparticle densities can lead to significant performance degradation. Hence, ideal materials for these devices would be less susceptible to quasiparticle poisoning. However, existing models of these devices often rely on approximations and phenomenology. Therefore, they lack a rigorous description of the underlying quasiparticle and phonon dynamics that are responsible for device performance. In this article, we combine kinetic equations with density functional theory to model the nonequilibrium quasiparticle and phonon dynamics of a thin superconducting film ab initio. To demonstrate the universality of our model, we illustrate two independent example applications: (1) we develop a theoretical model describing the detection of single photons in superconducting nanowires, and (2) we calculate the energy-relaxation rate of a transmon qubit due to the presence of excess quasiparticles. Our examples demonstrate from first principles that niobium nitride is well suited to be used for single-photon detection and that tantalum transmon qubits possess reduced sensitivity to quasiparticle poisoning relative to other materials, which is likely in part responsible for their longer coherence times. In contrast to previous models of superconducting devices, our ab initio approach makes predictions of device performance without experimental input and thus can be used to accelerate progress in device development. Moreover, by considering the full-bandwidth electron-phonon coupling, our approach can incorporate strong-coupling effects. Our methods effectively integrate ab initio materials modeling with nonequilibrium theory of superconductivity to perform practical modeling of superconducting devices, providing a comprehensive approach that connects fundamental theory with device-level applications.
This dissertation describes energy dissipation and microwave loss due to non-equilibrium quasiparticles in superconducting transmon qubits and titanium nitride coplanar waveguide resonators. During the measurements of transmon T1 relaxation time and resonator quality factor QI, I observed reduced microwave loss as the temperature increased from 20 mK to approximately Tc/10 at which the loss takes on a minimum value. I argue that this effect is due to non-equilibrium quasiparticles. I measured the temperature dependence of the relaxation time T1 of the excited state of an Al/AlOx/Al transmon and found that, in some cases, T1 increased by almost a factor of two as the temperature increased from 30 mK to 100 mK with a best T1 of 0.2 ms. I present an argument showing this unexpected temperature dependence occurs due to the behavior of non-equilibrium quasiparticles in devices in which one electrode in the tunnel junction has a smaller volume, and slightly smaller superconducting energy gap, than the other electrode. At sufficiently low temperatures, non-equilibrium quasiparticles accumulate in the electrode with the smaller gap, leading to a relatively high density of quasiparticles at the junction and a short T1. Increasing the temperature gives the quasiparticles enough thermal energy to occupy the higher gap electrode, reducing the density at the junction and increasing T1. I present a model of this effect, extract the density of quasiparticles and the two superconducting energy gaps, and discuss implications for increasing the relaxation time of transmons. I also observed a similar phenomenon in low temperature microwave studies of titanium nitride coplanar waveguide resonators. I report on loss in a resonator at temperatures from 20 mK up to 1.1 K and with the application of infrared pair breaking radiation (λ=1.55 μm). With no applied IR light, the internal quality factor increased from QI = 800,000 at T < 70 mK up to QI=(2×10^6 ) at 600 mK. The resonant frequency f0 increased by 2 parts per million over the same temperature range. Above 600 mK both QI and f0 decreased rapidly, consistent with the increase in the density of thermally generated quasiparticles. With the application of IR light and for intensities below 1 aW μm^(-2) and T < 400 mK, QI increased in a similar way to increasing the temperature before beginning to decrease with larger intensities. I show that a model involving non-equilibrium quasiparticles and two regions of different superconducting gaps can explain this unexpected behavior.