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 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 kinetic inductance detectors (KIDs) and other similar applications of superconducting microresonators, both the large and small-signal behaviour of the device may be affected by electrothermal feedback. Microwave power applied to read out the device is absorbed by and heats the superconductor quasiparticles, changing the superconductor conductivity and hence the readout power absorbed in a positive or negative feedback loop. In this work, we explore numerically the implications of an extensible theoretical model of a generic superconducting microresonator device for a typical KID, incorporating recent work on the power flow between superconductor quasiparticles and phonons. This model calculates the large-signal (changes in operating point) and small-signal behaviour of a device, allowing us to determine the effect of electrothermal feedback on device responsivity and noise characteristics under various operating conditions. We also investigate how thermally isolating the device from the bath, for example by designing the device on a membrane only connected to the bulk substrate by thin legs, affects device performance. We find that at a typical device operating point, positive electrothermal feedback reduces the effective thermal conductance from the superconductor quasiparticles to the bath, and so increases responsivity to signal (pair-breaking) power, increases noise from temperature fluctuations, and decreases the noise equivalent power (NEP). Similarly, increasing the thermal isolation of the device while keeping the quasiparticle temperature constant decreases the NEP, but also decreases the device response bandwidth.
Abstract In kinetic inductance detectors (KIDs) and other similar applications of superconducting microresonators, both the large and small-signal behaviour of the device may be affected by electrothermal feedback. Microwave power applied to read out the device is absorbed by and heats the superconductor quasiparticles, changing the superconductor conductivity and hence the readout power absorbed in a positive or negative feedback loop. In this work, we explore numerically the implications of an extensible theoretical model of a generic superconducting microresonator device for a typical KID, incorporating recent work on the power flow between superconductor quasiparticles and phonons. This model calculates the large-signal (changes in operating point) and small-signal behaviour of a device, allowing us to determine the effect of electrothermal feedback on device responsivity and noise characteristics under various operating conditions. We also investigate how thermally isolating the device from the bath, for example by designing the device on a membrane only connected to the bulk substrate by thin legs, affects device performance. We find that at a typical device operating point, positive electrothermal feedback reduces the effective thermal conductance from the superconductor quasiparticles to the bath, and so increases responsivity to signal (pair-breaking) power, increases noise from temperature fluctuations, and decreases the noise equivalent power (NEP). Similarly, increasing the thermal isolation of the device while keeping the quasiparticle temperature constant decreases the NEP, but also decreases the device response bandwidth.
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.
Understanding the detailed behaviour of superconducting pair breaking photon detectors such as Kinetic Inductance Detectors (KIDs) requires knowledge of the nonequilibrium quasiparticle energy distributions. We have previously calculated the steady state distributions resulting from uniform absorption of monochromatic sub gap and above gap frequency radiation by thin films. In this work, we use the same methods to calculate the effect of illumination by broadband sources, such as thermal radiation from astrophysical phenomena or from the readout system. Absorption of photons at multiple above gap frequencies is shown to leave unchanged the structure of the quasiparticle energy distribution close to the superconducting gap. Hence for typical absorbed powers, we find the effects of absorption of broadband pair breaking radiation can simply be considered as the sum of the effects of absorption of many monochromatic sources. Distribution averaged quantities, like quasiparticle generation efficiency η, match exactly a weighted average over the bandwidth of the source of calculations assuming a monochromatic source. For sub gap frequencies, however, distributing the absorbed power across multiple frequencies does change the low energy quasiparticle distribution. For moderate and high absorbed powers, this results in a significantly larger η–a higher number of excess quasiparticles for a broadband source compared to a monochromatic source of equal total absorbed power. Typically in KIDs the microwave power absorbed has a very narrow bandwidth, but in devices with broad resonance characteristics (low quality factors), this increase in η may be measurable.
The response of superconducting pair-breaking detectors is dependent on the details of the quasiparticle distribution. In Kinetic Inductance Detectors (KIDs), where both pair breaking and non-pair breaking photons are absorbed simultaneously, calculating the detector response therefore requires knowledge of the often nonequilibrium distributions. The quasiparticle effective temperature provides a good approximation to these nonequilibrium distributions. We compare an analytical expression relating absorbed power and the quasiparticle effective temperature in superconducting thin films to full solutions for the nonequilibrium distributions, and find good agreement for a range of materials, absorbed powers, photon frequencies and temperatures typical of KIDs. This analytical expression allows inclusion of nonequilibrium effects in device models without solving for the detailed distributions. We also show our calculations of the frequency dependence of the detector response are in agreement with recent experimental measurements of the response of Ta KIDs at THz frequencies.
Abstract Understanding the detailed behaviour of superconducting pair breaking photon detectors such as Kinetic Inductance Detectors (KIDs) requires knowledge of the nonequilibrium quasiparticle energy distributions. We have previously calculated the steady state distributions resulting from uniform absorption of monochromatic sub gap and above gap frequency radiation by thin films. In this work, we use the same methods to calculate the effect of illumination by broadband sources, such as thermal radiation from astrophysical phenomena or from the readout system. Absorption of photons at multiple above gap frequencies is shown to leave unchanged the structure of the quasiparticle energy distribution close to the superconducting gap. Hence for typical absorbed powers, we find the effects of absorption of broadband pair breaking radiation can simply be considered as the sum of the effects of absorption of many monochromatic sources. Distribution averaged quantities, like quasiparticle generation efficiency η , match exactly a weighted average over the bandwidth of the source of calculations assuming a monochromatic source. For sub gap frequencies, however, distributing the absorbed power across multiple frequencies does change the low energy quasiparticle distribution. For moderate and high absorbed powers, this results in a significantly larger η –a higher number of excess quasiparticles for a broadband source compared to a monochromatic source of equal total absorbed power. Typically in KIDs the microwave power absorbed has a very narrow bandwidth, but in devices with broad resonance characteristics (low quality factors), this increase in η may be measurable.