Abstract
Characterizing superconducting microwave resonators with highly dissipative elements is a technical challenge, but a requirement for implementing and understanding the operation of hybrid quantum devices involving dissipative elements, e.g., for thermal engineering and detection. We present experiments on λ/4 superconducting niobium coplanar waveguide resonators, terminating at the antinode by a dissipative copper microstrip via aluminum leads, such that the resonator response is difficult to measure in a typical microwave environment. By measuring the transmission both above and below the superconducting transition of aluminum, we are able to isolate the resonance. We then experimentally verify this method with copper microstrips of increasing thicknesses, from 50 nm to 150 nm, and measure quality factors in the range of 10-67 in a consistent way.
| Original language | English |
|---|---|
| Article number | 022601 |
| Journal | Applied Physics Letters |
| Volume | 115 |
| Issue number | 2 |
| DOIs | |
| Publication status | Published - 8 Jul 2019 |
| MoE publication type | A1 Journal article-refereed |
Funding
https://orcid.org/0000-0001-9101-0973 Chang Yu-Cheng 1,2,3 a) Karimi Bayan 3 Senior Jorden 3 Ronzani Alberto 3 https://orcid.org/0000-0002-3903-7623 Peltonen Joonas T. 3 Goan Hsi-Sheng 1 Chen Chii-Dong 2,3 Pekola Jukka P. 3 1 Department of Physics, National Taiwan University , Taipei 106, Taiwan 2 Institute of Physics, Academia Sinica , Taipei 115, Taiwan 3 QTF Centre of Excellence, Department of Applied Physics, Aalto University , P.O. Box 13500, FI-00076 Aalto, Finland a) Electronic mail: [email protected] 08 07 2019 115 2 022601 01 04 2019 18 06 2019 08 07 2019 © 2019 Author(s). 2019 Author(s) 0003-6951/2019/115(2)/022601/5 All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license ( http://creativecommons.org/licenses/by/4.0/ ). Characterizing superconducting microwave resonators with highly dissipative elements is a technical challenge, but a requirement for implementing and understanding the operation of hybrid quantum devices involving dissipative elements, e.g., for thermal engineering and detection. We present experiments on λ /4 superconducting niobium coplanar waveguide resonators, terminating at the antinode by a dissipative copper microstrip via aluminum leads, such that the resonator response is difficult to measure in a typical microwave environment. By measuring the transmission both above and below the superconducting transition of aluminum, we are able to isolate the resonance. We then experimentally verify this method with copper microstrips of increasing thicknesses, from 50 nm to 150 nm, and measure quality factors in the range of 10–67 in a consistent way. Academy of Finland 297240 H2020 Marie Skłodowska-Curie Actions http://doi.org/10.13039/100010665 742559 766025 Academy of Finland http://doi.org/10.13039/501100002341 303677 312057 crossmark Superconducting microwave resonators are the cornerstone of much of the current state-of-the art superconducting quantum technologies. Their intrinsic electromagnetic properties routinely enable high quality factors, typically over 10 000, 1 allowing for extremely sensitive measurements, for example, the multiplexed readout of the dispersive shift of a weakly coupled quantum bit, 2–4 a fundamentally important tool for quantum information, or photon absorption of a microwave kinetic inductance detectors (MKIDs), 5–8 a tool widely used in astronomy. On the other hand, dissipative elements added to superconducting microwave circuits are finding increasing applications in the fields of microwave amplification, 9,10 circuit quantum environmental engineering, 11–13 quantum information, 14 and circuit quantum thermodynamics. 15–17 In superconducting circuits, normal metal elements can be easily integrated into the existing chip architecture and provide channels for intentional decoherence, such as for quantum bit initialization, 11–13 precise temperature control (electronic heating and cooling), 18,19 and ultrasensitive calorimetry. 20 Notably, these devices are in use as high-speed thermometers and microwave photon sources. 21,22 Experimentally characterizing such dissipative resonators is challenging by conventional single tone spectroscopy, as the dissipation results in a broad and shallow resonance, often unresolvable from the electronic noise of a cryogenic high electron mobility transistor (HEMT) amplifier and the frequency dependent variations of transmission in the experimental setup. In this report, we present a method for this task utilizing an intermediate superconductor with a lower critical temperature, enabling us to isolate the resonance by performing characterization at differing bath temperatures, with measured quality factors as low as 10. The method presented is demonstrated here using a structure consisting of two λ /4 niobium coplanar waveguide (CPW) resonators, both inductively coupled to a common CPW transmission line, shown in Fig. 1(a) , used for multiplexed readout. On the left is a fully superconducting λ /4 CPW resonator to act as a reference. On the right resonator, the voltage node is terminated by an aluminum-copper-aluminum constriction, forming a Nb/Al/Cu/Al/Nb heterostructure, with the copper acting as the dissipative normal metal shown in the micrograph in Fig. 1(b) . FIG. 1. Device and measurement setup. (a) The image shows two λ /4 resonators inductively coupled to a transmission line. The left resonator is a standard λ /4 superconducting resonator with one end of the center conductor directly shorted to the ground plane on chip, and the other end is an open circuit to the ground. This resonator acts as a reference. The right side resonator is a copper-terminated (Al-Cu-Al) λ /4 resonator. At right and left bottom, four test junctions are made for characterizing the DC electronic properties. (b) Scanning-electron microscopy image of the yellow outlined region of (a), highlighting the Al/Cu/Al junction in contact with the center conductor of the resonator at the right side and with ground at the left side. The image is colored with orange for copper, cyan for aluminum, and purple for niobium. (c) Equivalent lumped circuit for the copper terminated λ /4 resonator. The coupled-inductor symbol represents the dominant coupling mechanism between the in/output transmission line and the resonator. The resistance R represents the copper termination. (d) Measurement setup. The microwave signal is introduced from port 1 by a vector network analyzer (VNA) at room temperature and passes through several attenuators distributed at different temperatures to the input of the device at 10 mK. The output microwave is measured at the port 2 of VNA through two isolators and two amplifiers. Samples are fabricated on a 330 μ m thick c-plane sapphire substrate using a process described in Ref. 17 . The sapphire surface is initially cleaned with an argon ion plasma milling before a 200 nm thick niobium film is deposited by DC magnetron sputtering. The coplanar waveguide resonator patterns are written by electron-beam lithography (EBL) and transferred to the niobium film using an SF 6 + O 2 reactive ion etching process. As the niobium layer is thick, subsequently comparatively thin evaporated aluminum and copper films are fragile and can become discontinuous at the intersection. Therefore, during the EBL exposure, the dose at the interface is incrementally changed, forming a ramp when etched to remove the discontinuity and increase the surface area. The Al/Cu/Al terminations are written by electron-beam lithography onto a bilayer resist and grown using double-angle deposition in an electron-beam evaporator, with the galvanic contact between lithographic layers facilitated by an in situ argon ion plasma milling process removing native oxides on the niobium. At the beginning of deposition, 10 nm of aluminum is evaporated to improve adhesion to the sapphire surface, followed by copper (of variable thickness; 50 nm, 100 nm, and 150 nm for different devices). Finally, two 110 nm thick aluminum contacts to the niobium CPW are deposited. After processing, the substrates are diced with a diamond-embedded resin blade, wire-bonded to the sample stage assembly described in Ref. 23 . Finally, the sample is loaded into a cryogen-free dilution refrigerator with a base temperature of 10 mK. All spectroscopic measurements are performed using a vector network analyzer (VNA) at room temperature with the signal reaching the sample via an attenuated microwave line. The attenuation is distributed at various temperature stages of the dilution refrigerator, as shown in Fig. 1(d) . This setup reduces Johnson-Nyquist noise from the attenuators at higher temperature stages. 24 The signal leaving the sample is then passed through two circulators at base temperature, to a low noise HEMT amplifier mounted at 4 K, providing 80 dB isolation from the amplifier input. Outside of the dilution refrigerator, the signal is passed through an additional 30 dB amplifier, before it is received by the VNA, capturing both I and Q quadratures, from which we reconstitute the transmission S 21 . The transmission through the system reads S 21 notch = a e i α e − 2 π if τ [ 1 − ( Q l / | Q c | ) e i ϕ 1 + 2 i Q l ( f / f r − 1 ) ] . (1) Here, a is the overall amplitude, α is the phase shift contributed by various components in the circuit, τ is the electronic delay caused by the length of the cable and the finite speed of light, and ϕ quantifies the impedance mismatch. 25,26 Parameters f and f r denote the probe frequency and the resonance frequency of the resonator, and Q l and Q c are the loaded quality factor and the coupling quality factor, respectively. 27 The quality factor is the ratio of energy stored in a resonator to average energy lost per cycle. The photon loss rate, inversely proportional to the quality factor, is a linear combination of the internal losses of the resonator and corresponds to the losses arising in the coupling from the resonator to the transmission line. The inverse of the loaded quality factor reads 1 / Q l = 1 / Q i + 1 / Q c . The Q l / Q c ratio determines the depth of the notched transmission and is maximized at Q l ≈ Q c at the resonance frequency. In order to increase the measured signal for a lossy resonator, we require Q c ≪ Q i , meaning that there are more photons leaving the resonator than dissipated ones. In the samples measured, a quality factor of 20 was desired corresponding to favorable operation of the envisioned device, the quantum heat engine. 16,17 The quality factor of our samples can be estimated by a simple model which terminates the center conductor of the CPW transmission line and the ground through a lumped resistance R at the voltage node and capacitively to the ground at the current node. The impedance of our λ /4 transmission line, shown in Fig. 1(c) , is given near this resonance by Z = − i π Z ∞ ( f / f r − f r / f ) / 4 , where Z ∞ = 50 Ω is the characteristic impedance of an infinite transmission line. Power injected to the termination at the half-power points f = f r ± Δ f / 2 can be written as P ≃ V 2 / ( R ( 1 + ( π Z ∞ / 4 R ) 2 ( Δ f / f r ) 2 ) ) = V 2 / 2 R , where Δ f is the width of the resonance peak. By definition, the internal quality factor Q i of the resistively terminated resonator can be written as Q i = π Z ∞ 4 R . (2) To measure the low Q i , coupling should be increased ( Q c ≪ Q i ), 28 in order to maximize the measurable signal. For us, however, even though the coupler covers half of the cavity length to achieve the strong coupling regime, Q l is still dominated by Q i . Due to the low overall quality factor, the depth of the notch in S 21 is of order 1 dB, which typically cannot be resolved within the background of the microwave setup. To isolate the resonator, the method presented here is based on measuring the spectrum above and below the critical temperature of aluminum T c , Al . The lower quality factor Q i of the resonator above T c , Al suppresses the resonance. Under these high temperature conditions, the transmittance can be regarded as a reference measurement for the background. The transmission through the CPW line is shown in Fig. 2(b) , with the red and blue traces measured at different temperatures T H and T L , above and below the critical temperature of the aluminum T c , Al , respectively, yet staying sufficiently below the critical temperature T c , Nb of the niobium transmission line. Here, T c , Al and T c , Nb can be determined by the resistance-temperature characteristic of a coprocess sample shown in the inset of Fig. 2(a) . A narrow band of S 21 ( T L ) is shown in the right inset of Fig. 2(b) , emphasizing the narrow bandwidth of the standard reference resonator. Here, f 0 is extracted as 7.246 GHz with a corresponding 5 MHz shift with temperature, as the kinetic inductance increases with decreasing temperature. 29,30 The dissipative resonator is not clearly visible in either trace. We observe, however, that as expected, due to lower operating temperature and thus more ideal superconducting characteristics of the CPW transmission line, the blue trace is consistently higher than the red trace, except in the region highlighted in yellow, centered around the design frequency of the dissipative resonator. By taking the ratio of the two S 21 measured at T L and T H , one can isolate the resonance of the dissipative resonator. This ratio is shown in Fig. 2(b) . Fitting a standard notched resonator model 27 to the trace, the parameters f r , Q c , and Q i can be extracted, with data for different samples (with variable thicknesses of copper) presented in Table I . FIG. 2. (a) Transmittance spectra of the device with 100 nm thick copper termination. The red and blue curves are measured at T H = 2 K and T L = 10 mK, respectively, as indicated by arrows in the inset. Here, T L and T H satisfy T L ≪ T c , Al < T H < T c , Nb . The inset shows resistance-temperature characteristics of a nominally identical coprocess dissipative element, demonstrating three plateaus corresponding to the superconducting transitions of aluminum and niobium, respectively. The yellow region demonstrates the frequency range of the resistively terminated superconducting resonator. We present in the right inset of (b) the quality factor of 2.8 × 10 4 at resonant frequency f 0 = 7.246 GHz of the reference resonator. (b) The amplitude of S 21 ( T L ) / S 21 ( T H ) and the phase in the left inset. Within the yellow highlighted frequency region, we see the resonance of the dissipative resonator. Here, we extracted f r and Q i as 6.71 GHz and 45 corresponding to the green dashed fitting line. TABLE I. Parameters of resistively terminated superconducting resonators extracted based on the measured S 21 . Thickness 50 nm 100 nm 150 nm 100 nm a Sample A-1 B-1 B-2 B-3 C-1 D-1 D-2 f r (GHz) 6.54 6.74 6.69 6.70 6.53 6.67 6.67 Q c 330 340 200 250 320 150 460 Q i 10.3 45.3 53.2 44.7 66.3 20.4 26.5 Q i , err 0.5 0.5 0.9 0.3 1.3 0.2 0.8 a Samples with increased contact resistance between the niobium and aluminum layers. In order to investigate how the Q i of the resistively terminated superconducting resonator depends on its resistance R , the coprocess samples with identical termination elements [see the red dashed frame in Fig. 1(a) ] are measured in a current biased four probe configuration at a 50 mK bath temperature. A measured IV curve is shown for the 150 nm thick copper termination in Fig. 3(a) . By ramping the bias current up through the termination, we observe first the superconducting energy gap and then two resistive branches at higher bias. The superconducting state exists by virtue of the proximity effect in the superconductor-normal metal-superconductor (SNS) Josephson junction induced by the Al/Cu/Al structure, 31,32 and the termination switches to the first dissipative branch at current I sw 1 = 5.6 μ A with normal state resistance R = 0.87 Ω. The normal resistance R due to the copper wire and the imperfect contact between the aluminum and the copper layers determines the quality factor of the resonator in the microwave measurement at T L . Continually increasing the current, the resistance increases at I sw 2 = 43.3 μ A with resistance R > = 4.83 Ω when the current exceeds the critical current of aluminum. R > consists of R , the normal state resistance of aluminum wire and the Al/Nb contact. All the measured parameters I sw 1 , I sw 2 , R , and R > extracted from IV measurement are given in Table II . FIG. 3. (a) and (b) VI measurement for the termination element of Nb/Al/Cu/Al/Nb of the 150 nm thick copper. I sw 1 is the critical current of the SNS junction and I sw 2 denotes the critical current of aluminum. Applying current below I sw 2 and over I sw 1 leads to the slope marked with a light blue dashed line originating from the resistance R of copper wire and the Cu/Al contact resistance. At a bias current exceeding the critical current of the aluminum wire, the copper, aluminum, Cu/Al contact resistance, and Nb/Al contact resistance all contribute to the total resistance R > . (c) Q i extracted by fitting two models, notch S 21 and divided notch S 21 / S 21 ( R > ) vs the reciprocal measured termination resistance R in blue dots and green hollow circles, respectively. The orange solid line denotes theoretical Q i as a function of reciprocal termination resistance. TABLE II. DC properties of Nb/Al/Cu/Al/Nb junctions. Thickness 50 nm 100 nm 150 nm 100 nm a Sample A-2 A-3 B-4 B-5 C-2 C-3 D-3 D-4 I sw 1 ( μ A) 0.2 0.3 2.1 1.8 5.6 2.6 0.74 1.0 I sw 2 ( μ A) 19.6 22.8 39.0 40.0 43.3 42.0 21.0 21.1 R (Ω) 5.3 4.3 1.7 1.6 0.9 1.3 3.3 2.8 R > (Ω) 21.5 16.0 8.2 7.9 4.8 5.2 19.4 19.5 a Samples with increased contact resistance between the niobium and aluminum layers. The presence of supercurrent in SNS junctions is due to the well-known proximity effect. 33–36 Based on the ratio of the Thouless energy E T = ℏ D / L 2 to the superconducting gap of aluminum Δ, the junction is in the long junction regime, when E T Δ ≪ 1 and the zero temperature eRI c is found to be proportional to E T in this limit. Here, L is the length of the junction, D = v F l e / 3 is the diffusion constant of the N metal, v F is the Fermi velocity, I c is the critical current of the junction, and l e is the elastic mean free path of electrons. In Table II , eR I sw 1 is smaller than 5 μ eV, indicating a long junction limit. The N lead between S superconductors of the resonator is an ideal element to localize heat. The small volume of normal metal enhances the temperature of the element at fixed transferred power, whereas the S forms a Cooper pair while injecting an electron from N and reflects a hole known as Andreev reflection that drops the temperature at the interface and localizes the heat efficiently in the N. 37 Depending on the application and operation regime, the quality factor Q i of the resistively terminated resonator can be designed merely by the resistance of the N element in the long junction regime. In Ref. 17 , f r / Q i is identified as the coupling between the copper heat bath and the resonator. In Fig. 3(c) , we plot the measured Q i from Table I as blue dots. The theoretically expected Q i is shown as an orange line. The expected slope by Eq. (2) is shown with Z ∞ set to 50 Ω. For a more general analysis of the internal quality factor, a model that takes into account the nonvanishing quality factor above the aluminum superconducting transition temperature can be used to simulate the measurement technique by inserting R > into Eq. (2) to obtain Q i ( T H ) . Q i ( T L ) and Q i ( T H ) are independent dissipation channels with 1 / Q l = 1 / Q i + 1 / Q c , which are then inserted into Eq. (1) to obtain the general two temperature model, S 21 ( T L ) / S 21 ( T H ) . Here, Q i ( T L ) is extracted from the model and plotted as green hollow circles. The slopes of extracted Q i vs 1/ R from the two models differ by less than 12%, which places the values within the uncertainty of the fitting algorithm. According to Table II , R > decreases while the thickness of the copper increases. The lower R > in the thicker samples yields a broader resonance when aluminum undergoes the transition to normal state, which explains the difference of extracted Q i for the thicker Cu films. However, the measured Q i is approximately double the expected value. The plausible origin of the discrepancy with respect to Eq. (2) is that the admittance of the SNS junction is composed of parallel dissipative (real) and reactive (imaginary) components, Y R and Y I , respectively. 38 In this case, the actual quality factor Q i is enhanced with respect to Eq. (2) by a factor ( Y R 2 + Y I 2 ) / Y R 2 , and while the contribution of the phase is challenging to measure, based on the I sw 1 and R values presented in Table II , we estimate that this can lead to an enhancement of approximately factor 2, which would explain the discrepancy. In conclusion, we have presented a technique for isolating the resonance of low-quality dissipative superconducting resonators based on the transition temperature of an intermediate superconductor with a lower energy gap as compared to that of the superconducting resonator. We have verified this method by characterizing a series of resistively terminated λ /4 superconducting CPW resonators. From microwave measurements, we have extracted extremely low quality factors of 10–67 for resonators with different resistances of the copper termination element. While similar characterization exploiting the superconducting transition could be achieved by supplying sufficient power to the resonator in excess of the critical current of aluminum or bypassed completely by using a microwave switch to measure a background reference, we believe the technique presented here is more versatile, applicable also in highly attenuated microwave input lines, such as those being used in circuit quantum electrodynamics and thermodynamics experiments and more accurately measuring the microwave background closer to the device being characterized. We acknowledge the Micronova Nanofabrication Centre of OtaNano infrastructure for providing the processing facilities and the European MicroKelvin Platform for the use of the dilution refrigerator. Y.-C.C. was supported by the visiting project of the Centre for Quantum Engineering, Finland and Institute of Physics, Academia Sinica, Taiwan (R.O.C.) during this experiment. B.K. acknowledges the grant by Marie Sklodowska-Curie actions (Grant Agreement Nos. 742559 and 766025). 1. A. Megrant , C. Neill , R. Barends , B. Chiaro , Y. Chen , L. Feigl , J. Kelly , E. Lucero , M. Mariantoni , P. J. O'Malley et al. , Appl. Phys. Lett. 100 , 113510 ( 2012 ). 10.1063/1.3693409 2. J. Majer , J. Chow , J. Gambetta , J. Koch , B. Johnson , J. Schreier , L. Frunzio , D. Schuster , A. Houck , A. Wallraff et al. , Nature 449 , 443 ( 2007 ). 10.1038/nature06184 3. Y. Chen , D. Sank , P. O'Malley , T. White , R. Barends , B. Chiaro , J. Kelly , E. Lucero , M. Mariantoni , A. Megrant et al. , Appl. Phys. Lett. 101 , 182601 ( 2012 ). 10.1063/1.4764940 4. M. Jerger , S. Poletto , P. Macha , U. Hübner , E. Ilichev , and A. V. Ustinov , Appl. Phys. Lett. 101 , 042604 ( 2012 ). 10.1063/1.4739454 5. L. Swenson , A. Cruciani , A. Benoit , M. Roesch , C. Yung , A. Bideaud , and A. Monfardini , Appl. Phys. Lett. 96 , 263511 ( 2010 ). 10.1063/1.3459142 6. S. McHugh , B. A. Mazin , B. Serfass , S. Meeker , K. O'Brien , R. Duan , R. Raffanti , and D. Werthimer , Rev. Sci. Instrum. 83 , 044702 ( 2012 ). 10.1063/1.3700812 7. R. Janssen , J. Baselmans , A. Endo , L. Ferrari , S. Yates , A. Baryshev , and T. Klapwijk , Appl. Phys. Lett. 103 , 203503 ( 2013 ). 10.1063/1.4829657 8. J. van Rantwijk , M. Grim , D. van Loon , S. Yates , A. Baryshev , and J. Baselmans , IEEE Trans. Microwave Theory Tech. 64 , 1876 ( 2016 ). 10.1109/TMTT.2016.2544303 9. D. Hover , S. Zhu , T. Thorbeck , G. Ribeill , D. Sank , J. Kelly , R. Barends , J. M. Martinis , and R. McDermott , Appl. Phys. Lett. 104 , 152601 ( 2014 ). 10.1063/1.4871088 10. M. Mück , B. Schmidt , and J. Clarke , Appl. Phys. Lett. 111 , 042604 ( 2017 ). 10.1063/1.4985384 11. P. J. Jones , J. A. M. Huhtamäki , J. Salmilehto , K. Y. Tan , and M. Möttönen , Sci. Rep. 3 , 1987 ( 2013 ). 10.1038/srep01987 12. J. Tuorila , M. Partanen , T. Ala-Nissila , and M. Möttönen , npj Quantum Inf. 3 , 27 ( 2017 ). 10.1038/s41534-017-0027-1 13. M. Partanen , K. Y. Tan , S. Masuda , J. Govenius , R. E. Lake , M. Jenei , L. Grönberg , J. Hassel , S. Simbierowicz , V. Vesterinen , J. Tuorila , T. Ala-Nissila , and M. Möttönen , Sci. Rep. 8 , 6325 ( 2018 ). 10.1038/s41598-018-24449-1 14. F. Verstraete , M. M. Wolf , and J. I. Cirac , Nat. Phys. 5 , 633 ( 2009 ). 10.1038/nphys1342 15. A. O. Niskanen , Y. Nakamura , and J. P. Pekola , Phys. Rev. B 76 , 174523 ( 2007 ). 10.1103/PhysRevB.76.174523 16. B. Karimi and J. P. Pekola , Phys. Rev. B 94 , 184503 ( 2016 ). 10.1103/PhysRevB.94.184503 17. A. Ronzani , B. Karimi , J. Senior , Y.-C. Chang , J. T. Peltonen , C. Chen , and J. P. Pekola , Nat. Phys. 14 , 991 ( 2018 ). 10.1038/s41567-018-0199-4 18. M. Nahum , T. M. Eiles , and J. M. Martinis , Appl. Phys. Lett. 65 , 3123 ( 1994 ). 10.1063/1.112456 19. F. Giazotto , T. T. Heikkilä , A. Luukanen , A. M. Savin , and J. P. Pekola , Rev. Mod. Phys. 78 , 217 ( 2006 ). 10.1103/RevModPhys.78.217 20. S. Gasparinetti , K. L. Viisanen , O.-P. Saira , T. Faivre , M. Arzeo , M. Meschke , and J. P. Pekola , Phys. Rev. Appl. 3 , 014007 ( 2015 ). 10.1103/PhysRevApplied.3.014007 21. O.-P. Saira , M. Zgirski , K. L. Viisanen , D. S. Golubev , and J. P. Pekola , Phys. Rev. Appl. 6 , 024005 ( 2016 ). 10.1103/PhysRevApplied.6.024005 22. S. Masuda , K. Y. Tan , M. Partanen , R. E. Lake , J. Govenius , M. Silveri , H. Grabert , and M. Möttönen , Sci. Rep. 8 , 3966 ( 2018 ). 10.1038/s41598-018-21772-5 23. R. E. George , J. Senior , O.-P. Saira , J. P. Pekola , S. E. de Graaf , T. Lindström , and Y. A. Pashkin , J. Low Temp. Phys. 189 , 60 ( 2017 ). 10.1007/s10909-017-1787-x 24. V. J. Urick , K. J. Williams , and J. D. McKinney , Fundamentals of Microwave Photonics ( John Wiley & Sons , 2015 ), Vol. 1 . 25. M. S. Khalil , M. J. A. Stoutimore , F. C. Wellstood , and K. D. Osborn , J. Appl. Phys. 111 , 054510 ( 2012 ). 10.1063/1.3692073 26. C. Deng , M. Otto , and A. Lupascu , J. Appl. Phys. 114 , 054504 ( 2013 ). 10.1063/1.4817512 27. S. Probst , F. B. Song , P. A. Bushev , A. V. Ustinov , and M. Weides , Rev. Sci. Instrum. 86 , 024706 ( 2015 ). 10.1063/1.4907935 28. M. Göppl , A. Fragner , M. Baur , R. Bianchetti , S. Filipp , J. M. Fink , P. J. Leek , G. Puebla , L. Steffen , and A. Wallraff , J. Appl. Phys. 104 , 113904 ( 2008 ). 10.1063/1.3010859 29. R. Meservey and P. M. Tedrow , J. Appl. Phys. 40 , 2028 ( 1969 ). 10.1063/1.1657905 30. K. Inomata , T. Yamamoto , M. Watanabe , K. Matsuba , and J.-S. Tsai , J. Vac. Sci. Technol., B 27 , 2286 ( 2009 ). 10.1116/1.3232301 31. J. Wei , P. Cadden-Zimansky , V. Chandrasekhar , and P. Virtanen , Phys. Rev. B 84 , 224519 ( 2011 ). 10.1103/PhysRevB.84.224519 32. L. Angers , F. Chiodi , G. Montambaux , M. Ferrier , S. Guéron , H. Bouchiat , and J. C. Cuevas , Phys. Rev. B 77 , 165408 ( 2008 ). 10.1103/PhysRevB.77.165408 33. J. Clarke , Proc. R. Soc. London, Ser. A 308 , 447 ( 1969 ). 10.1098/rspa.1969.0020 34. J. G. Shepherd , Proc. R. Soc. London, Ser. A 326 , 421 ( 1972 ). 10.1098/rspa.1972.0018 35. P. G. de Gennes , Rev. Mod. Phys. 36 , 225 ( 1964 ). 10.1103/RevModPhys.36.225 36. P. Dubos , H. Courtois , B. Pannetier , F. K. Wilhelm , A. D. Zaikin , and G. Schön , Phys. Rev. B 63 , 064502 ( 2001 ). 10.1103/PhysRevB.63.064502 37. A. F. Andreev , Sov. Phys. JETP 20 , 1490 ( 1965 ). 38. P. Virtanen , F. S. Bergeret , J. C. Cuevas , and T. T. Heikkilä , Phys. Rev. B 83 , 144514 ( 2011 ). 10.1103/PhysRevB.83.144514