A Composite Phononic Crystal Design for Quasiparticle Lifetime Enhancement in Kinetic Inductance Detectors

A nanoscale phononic crystal filter (reflector) is designed for a kinetic inductance detector where the reflection band is matched to the quasiparticle recombination phonons with the aim to increase quasiparticle lifetime in the superconducting resonator. The inductor is enclosed by a 1-μm-wide phononic crystal membrane section with two simple hole patterns that each contain a partial spectral gap for various high-frequency phonon modes. The phononic crystal is narrow enough for low-frequency thermal phonons to propagate unimpeded. With 3D phonon scattering simulation over a 40 dB attenuation in transmitted power is found for the crystal, which is estimated to give a lifetime enhancement of nearly two orders of magnitude.


Introduction
Phononic crystals (PnCs) are periodic structures designed to modify properties of elastic and acoustic waves by the Bragg scattering mechanism or by localized resonances [1]. Applications have been proposed in many engineering and physics areas, including vibration and noise absorption, acoustic cloaking and metamaterials research [2]. Recently, with the ability to fabricate microscopic phononic crystals, thermal applications have also appeared benefiting from the wave interference of phonons at thermal wavelengths. This could be utilized for instance in improving thermoelectric devices by reduction in thermal conductivity, guiding and focusing heat in nanostructures and using heat in phononic information processing [3][4][5]. At sub-Kelvin temperatures, phononic crystals fabricated in a suspended membrane have been shown to reduce thermal conductance by several orders of magnitude [6,7], which is a key to improve sensitivity in many low-temperature bolometric detector technologies.
Kinetic inductance detectors (KID), on the other hand, are extremely sensitive non-bolometric photon detectors based on superconducting microwave resonators [8]. The operation principle of a KID is to break Cooper pairs into quasiparticles in a superconducting microwave resonator with signal photons, thus changing the impedance of the resonator, which can be measured with a coupled microwave readout line. Quasiparticles recombine quickly back into Cooper pairs by emitting phonons at frequencies ≥ 2 ∕h , corresponding to the energy gap of the superconducting material. The recombination phonons can break other Cooper pairs, but phonons escaping the system result in a lost signal limiting the responsivity of the detector. Simply removing material around the resonator does not help because a good thermal contact is desired so that low-frequency phonons due to dissipation of the microwave readout signal can escape. Recently, Rostem et al. presented calculations for an effective recombination lifetime enhancement of quasiparticles in a phononic crystal isolated kinetic inductance detector (PhiKID) [9]. They showed how a PnC notch filter with a 0.3 ∕h wide reflection band and a 40 dB attenuation positioned slightly off from the gap frequency 2 ∕h of a superconductor can increase the lifetime by two orders of magnitude.
A number of publications in the field of ultrasound physics discuss the formation of band gaps or scattering of elastic waves in plates with obstacles [10][11][12], which are useful starting points in finding suitable PnC geometries with full band gaps for the recombination phonons. However, it will be difficult to fabricate such PnCs for the most commonly used low-T c superconductors, such as aluminum, which has 2 ∕h ∼ 89 GHz. For example, creating a matched reflection band for aluminum with a thru-hole pattern in a suspended silicon nitride membrane would steer the fabrication toward very thin membranes ( ∼ 33 nm) and the smallest features in a few nanometer length scale [6]. Even for thin-film hafnium, which has 2 ∕h ∼ 32 GHz, a matching PnC design with a thru-hole pattern would require a 95-nm-thick membrane with a distance between holes only 5 nm.
In this paper, we propose a simulation-based design for a composite PnC consisting of two simple hole patterns with minimum feature size 50 nm, and a relatively low combined width of 1 μm. Each of the patterns possess partial spectral gaps for various modes, which effectively create a full reflection band for highfrequency athermal recombination phonons, while thermal low-frequency phonons pass through the structure unimpeded. Interestingly, a PnC will also reject athermal phonons released during the energy downconversion of quasiparticles excited to energies significantly greater than . The downconversion process releases phonons at the scale of the Debye frequency of the metal. The loss of these so-called hot phonons from the superconductor to the bath is known to limit the performance of many cryogenic detectors, including STJs [13], TESs [14] and KIDs [15]. A PnC can effectively recycle hot short wavelength phonons through a purely geometric scattering mechanism from the hole walls, due to the lack of material that is etched away. At the same time, due to the scale of the PnC pattern, large wavelength thermal phonons with frequencies below 10 GHz have unity transmission. Thus, a PnC may be effective at improving the performance of KIDs that utilize aluminum, or other metals with similar T c values, to detect visible and near-infrared photons.

Scattering Simulations
Phonon propagation and scattering in PnC membranes at sub-Kelvin temperatures is described with a 3D continuum linear elasticity theory [16]. For smooth membranes, i.e., where the surface roughness is small compared to the phonon wavelength, phonon mean free paths have been measured to exceed 100 μm [17], thus phonons can be treated as ballistic waves that can only scatter from the PnC structure. The membrane material can be assumed to be isotropic and homogeneous, and it is thus characterized by a constant density and two independent elastic material constants, often denoted by the Lame parameters and . The displacement field U = (u, v, w) in such system obeys where is the domain describing the infinite PnC membrane geometry (see Fig. 1 left) and F = F(x, y, z, t) is the driving force applied in the domain. On the domain boundaries (the top and bottom surfaces, and hole boundaries), we apply the stressfree boundary condition n ⋅ (U) = 0 . Here, (⋅) = ( ij (⋅)) , i, j = x, y, z is the stress tensor and n is the surface normal vector. Assuming time harmonic force F and solution U, the derivative 2 ∕ t 2 can be replaced by − 2 , where is the angular frequency of the force [16].
The time harmonic problem of Eq. (1) is solved in the frequency domain using the finite element method and COMSOL Multiphysics v5.3 software. When the incident wave strikes the PnC from a 90° angle ( Fig. 1 left), the scattered field is periodic in in the y-direction, and thus the infinite domain can be truncated as shown in Fig. 1 (right). Periodic boundary conditions for the solution U are set on the x-z boundaries so that unphysical coupling or conversion of the modes does not occur there. Finally, perfectly matched layers pml are created at the ends to absorb the scattered phonons. We remark that it is in principle possible to simulate the system in a full-scale 3D model without periodicity allowing variable incident angles, but this would greatly increase the computational cost.
To study reflection and mode conversion properties of PnCs, we aim to run the scattering simulations with pure shear and Lamb modes [16], which are the solutions for Eq. (1) in the regions far away from the PnC. However, a typical problem with excitation of pure modes is that the scattered field re-scatters from the emitter and thus interferes with the original incident phonon. Neglecting the interference can lead to incorrect evaluation of transmission and reflection coefficients. We solved this problem by taking advantage of the biorthogonality property of the modes in the force term F applied in the 3D subdomain ex [18,19]. If the solution U on ex is expressed as U = U + inc + U scat , where U scat = ∑ j j U − j and plus (minus) denotes waves propagating in the positive (negative) x-direction, then by the biorthogonality property for all mode indices j. Here, x (⋅) ∶= ( xx (⋅), xy (⋅), xz (⋅)) T is a vector of the stress tensor components in the x-direction, and the overline denotes the complex conjugate. In the weak formulation of the elasticity equation (1) (see e.g., [20]), the domain force term appears in the same form as the first term of the biorthogonality property (2), if the force term is set to F = x (U + inc ) and the functions U − j are replaced by the finite element basis functions. We then add another integral term in the weak formulation corresponding to the second term of the biorthogonality property, which together with the first term guarantees that all scattered wave components U − j in the solution U result to zero in these integrals and thus cannot backscatter from the subdomain ex . Phonon modes U ± j in the membrane can be calculated either numerically, or semi-analytically by solving the Rayleigh-Lamb dispersion equations [16]. After normalization, functions U ± j are inserted into the weak formulation integrals for a numerical solution of the scattering problem.
Similarly, the phonon transmission power can be evaluated by the biorthogonality relation [19]. Total mechanical power of the elastic field U through a boundary A can be expressed as the sum of the single mode U + j contributions P j , where

Results
To find a practical design for the PnC membrane with good reflective properties in the 32 GHz region (corresponding to the energy gap of thin-film hafnium), we numerically calculated phonon dispersion relations for several simple infinite cylindrical hole patterns in crystalline silicon membranes varying the lattice constant, hole filling fraction, membrane thickness and the lattice geometry. The designs with a full spectral gap at 32 GHz found this way were not feasible for fabrication with current lithography processes because of too small feature sizes going below 5-nm length scale. A full spectral gap, however, is not mandatory, because a similar effect can be achieved by combining several simpler structures in series each possessing a partial spectral gap for some Lamb modes. With the scattering simulation, we can then estimate the effect of the composite structure for each phonon mode. The most practical geometry found, shown in Fig. 2, is a composite thru-hole PnC design, which functions as a phononic multi-pole filter [21]. The geometry contains cylindrical holes in hexagonal and square lattices in a thin silicon membrane with thickness of 100 nm. The lattice constant for both lattices is 110 nm and the hole diameter is 60 nm. The square lattice is 5 cells and the hexagonal lattice 5 cells wide making the composite crystal only about 1 μm wide with the smallest feature (neck width) 50 nm. The corresponding dispersion spectra are also shown for each of the components. This geometry can be fabricated with a single exposure etch process.
It is immediately seen from the dispersion relations that at the low-frequency (sub-10 GHz) range, the spectra resemble a full-membrane spectrum [16], thus the PnC should have minimal effect on thermal conductance in the membrane. In addition, neither of the lattices seem to have a full spectral gap at 32 GHz. However, further analysis of the modes reveals partial spectral gaps that can be utilized for the same effect. We remark that pure Lamb modes do not exist in the PnCs, but in these lattices (anti-)symmetrically polarized incident waves are likely to remain (anti-) symmetric also in the PnC. In other words, if there are no propagating antisymmetric out-of-plane modes in the PnC, then it can be expected that incident antisymmetric out-of-plane Lamb modes are reflected from the PnC. To locate the partial  Fig. 2, we classified the modes based on a component analysis of the fields U. In Fig. 2, red corresponds to the (mainly) antisymmetric out-of-plane mode, green is the in-plane shear, blue is the symmetric out-of-plane mode, and black is the second-order antisymmetric out-of-plane mode. The gray highlighted area shows the direction perpendicular to the PnC, which can be simulated with the scattering model. Following the colored curves, partial spectral gaps can be identified for each of the mode types in this region.
Finally, the scattering simulation is used to calculate the attenuation level of the composite structure. In Fig. 3, the transmitted power of each of the propagating Lamb and horizontal shear modes is shown around the 32-GHz frequency region. Colors of the various curves denote different mode polarizations, and the incident mode color is shown at the top-left corner. For all modes, a fairly wide reflection band is created at 32 GHz, which is highlighted by blue. Cut-off bands for the second-order modes are highlighted by red, which means that these modes do not exist at low-frequency range. Incident modes in most cases forward-scatter to the same mode, but a small amount of coupling can be identified for the first-and the second-order antisymmetric out-of-plane modes as revealed by the secondary curves in those plots. We remark that changing the incident angle slightly will not cause significant qualitative changes to the transmitted power in Fig. 3, and similar attenuation levels and band locations are likely to exist for a range of incident angles. Effective quasiparticle lifetime enhancement for each phonon mode can be calculated by inserting reflection band widths and locations for each of curves in Fig. 3 into the equations in Rostem et al. [9]. A single effective reflection band location and width parameters for the total power cannot be formulated from our data, because the emitted phonon angular distribution and population in various mode types were not fully analyzed. In a PhiKID detector, however, the antisymmetric modes are likely to be dominant, because the metal thin-film emitter is located only on the top-side of the suspended membrane, thus making the source geometry antisymmetric relative to the membrane. Thus, for an estimate of the total effective lifetime enhancement, the first antisymmetric mode band gap position 32 GHz ( = 2 ∕h ) and width 5 GHz ( = 0.3 ∕h ) are shown in Fig. 3, and then the lifetime enhancement factor can be read from Rostem et al. [9] as shown in Fig. 5. With these numbers, the value of approaches 10 2 indicating roughly two orders of magnitude increase in quasiparticle lifetime.

Conclusion
Nanoscale phononic crystals (PnCs) appear as promising functional components in many low temperature detector technologies. Kinetic inductance detectors have been suggested to benefit from notch filter-type PnCs when the filter is tuned to the quasiparticle recombination phonon frequencies of the KID superconducting resonator. However, for many superconductor materials, fabricating such a filter is challenging even with state-of-the-art nanofabrication techniques, because of extremely small features required for the PnC. Using 3D elastic wave scattering simulations, we demonstrated how a composite PnC assembled from two simple hole patterns can reach at least 40 dB power attenuation level for the recombination phonons incident from the perpendicular direction, which has been estimated to result in effective quasiparticle lifetime enhancement of two orders of magnitude. The smallest feature in this design was only 50 nm and the structure was only 1 μm wide making it relatively simple to fabricate with modern electron beam lithography techniques. Based on these results, a phononic crystal isolated kinetic inductance detector could be realized in the near future.