Conventional readout of a superconducting nanowire single-photon detector (SNSPD) sets an upper bound on the output voltage to be the product of the bias current and the load impedance, IB × Zload, where Zload is limited to 50 Ω in standard r.f. electronics. Here, we break this limit by interfacing the 50 Ω load and the SNSPD using an integrated superconducting transmission line taper. The taper is a transformer that effectively loads the SNSPD with high impedance without latching. At the expense of reduced maximum counting rate, it increases the amplitude of the detector output while preserving the fast rising edge. Using a taper with a starting width of 500 nm, we experimentally observed a 3.6× higher pulse amplitude, 3.7× faster slew rate, and 25.1 ps smaller timing jitter. The results match our numerical simulation, which incorporates both the hotspot dynamics in the SNSPD and the distributed nature in the transmission line taper. The taper studied here may become a useful tool to interface high-impedance superconducting nanowire devices to conventional low-impedance circuits.

Superconducting nanowire single-photon detectors (SNSPDs) are the leading single-photon detection technology at infrared wavelengths.1,2 With exceptional performance,3–7 they have played essential roles in various applications, especially quantum information science8,9 and deep-space optical communication.10 

A common problem with SNSPDs is their low output voltage and signal-to-noise ratio (SNR), which has been a limiting factor in detector timing jitter.11 A simple lumped-circuit model dictates that the output voltage from the nanowire cannot exceed IB × Zload, where IB is the bias current and Zload is the load impedance.12IB is usually limited by the nanowire's switching current in the μA range. It can be increased by connecting multiple wires in parallel and utilizing the cascade switchings among them.13,14Zload is set by the input impedance of the coaxial cable and r.f. electronics, which is conventionally 50 Ω. To improve readout SNR, significant progress has been made in developing cryogenic amplifiers with low noise, dissipation, and cost, e.g., using silicon germanium and gallium arsenide transistors.15–18 Digital readout circuits built directly from superconducting electronics, such as nanocryotrons19 and single flux quantum (SFQ) circuits,20,21 have also been demonstrated. These integrated superconducting circuits have low noise and are scalable but usually require additional biasing and suffer from leakage current and crosstalk.

An alternative approach to increase the output signal is to increase Zload. Compared to a standard 50 Ω load, a high-impedance load is often more desirable—it not only increases the detector output but also enables direct mapping of hotspot resistance and photon number/energy resolution.22–24 However, high-impedance loading is difficult to achieve in practice. The lack of high-impedance coaxial cables makes it necessary to place the high-impedance amplifiers close to the detectors (at the low-temperature stage), which imposes a more stringent power budget. More importantly, even if a high-impedance amplifier is available,17 loading a standard SNSPD directly with high impedance can lead to latching.12 

In this work, without the need of high-impedance cryogenic amplifiers or any active circuit elements, we break the IB × 50 Ω limit by using an integrated superconducting transmission line taper. The taper gradually transforms its characteristic impedance from kΩ to 50 Ω, which effectively loads the SNSPD with a kΩ impedance without latching. We designed the taper to be a co-planar waveguide (CPW) and fabricated it from the same superconducting thin film as the SNSPD. Using a taper with a starting width of 500 nm and a nominal passband from 200 MHz (the taper is a high-pass filter), we experimentally observed 3.6× higher output voltage and no added noise compared to the non-tapered reference device. This voltage gain is equivalent to an 11 dB passive, dissipation-free cryogenic amplifier. The large inductance of the taper reduces the detector's maximum counting rate from 40.3 MHz to 9.1 MHz (estimated from reset time) but does not slow down its fast rising edge, resulting in an increased slew rate and reduced timing jitter (from 48.9 ps to 23.8 ps). The integrated impedance taper demonstrated here is useful for interfacing high-impedance nanowire-based devices to conventional low-impedance components, such as memory and electrical or optical modulators.

Figure 1(a) shows a circuit model of a conventional SNSPD readout circuit, where the detector is modeled as a kinetic inductor LK in series with a time-dependent variable resistor RN. When an incident photon triggers the detector, RN switches from 0 to ≈kΩ within ≈100 s of ps and diverts the bias current to the load. The evolution of RN is determined by the non-linear electrothermal feedback in the detector.12,25 The currents from the bias source (IB), in the nanowire (ID), and to the load (IL) simply follow Kirchhoff's law IL = IBID. The maximum IL is therefore limited to IB, corresponding to the case where RN pushes all the current out of the nanowire (ID = 0). The output voltage on the load thus cannot exceed IB × 50 Ω. In practice, due to the electro-thermal feedback,12ID usually has some remainder, which depends on the bias current, kinetic inductance, and thermal constants of the materials.

FIG. 1.

Circuit diagram and micrographs of the tapered SNSPD readout. (a) A circuit diagram of a conventional SNSPD readout. (b) A circuit diagram of a tapered readout. The taper loads the SNSPD with high impedance while interfacing at its other end to the readout electronics at 50 Ω, resulting in a larger output voltage. (c) Schematic diagram of a co-planar waveguide transmission line taper. When an electrical pulse is launched from the high-impedance end, its voltage drops but current increases while traveling towards the low-impedance end. (d) An optical micrograph of the integrated transmission line taper. Light area: NbN; violet outlines: substrate. (e) A scanning electron micrograph of the SNSPD. Dark area: NbN; light area: substrate.

FIG. 1.

Circuit diagram and micrographs of the tapered SNSPD readout. (a) A circuit diagram of a conventional SNSPD readout. (b) A circuit diagram of a tapered readout. The taper loads the SNSPD with high impedance while interfacing at its other end to the readout electronics at 50 Ω, resulting in a larger output voltage. (c) Schematic diagram of a co-planar waveguide transmission line taper. When an electrical pulse is launched from the high-impedance end, its voltage drops but current increases while traveling towards the low-impedance end. (d) An optical micrograph of the integrated transmission line taper. Light area: NbN; violet outlines: substrate. (e) A scanning electron micrograph of the SNSPD. Dark area: NbN; light area: substrate.

Close modal

Figure 1(b) shows a simplified circuit diagram for the tapered readout. The taper is inserted between the SNSPD and the load, with a low impedance of ZL = 50 Ω on the load end and a high impedance of ZH on the detector end. In our implementation, the taper consists of a continuous nanowire transmission line without any dissipative elements. The taper is high-pass—it works as a transformer at high frequency but acts as an inductor at low frequency. When an incident photon triggers the SNSPD, RN switches on and pushes the current away from the nanowire. Instead of diverting the current directly to the 50 Ω load as in the conventional readout, the SNSPD injects current to the taper at ZH. As the electrical pulse travels towards the low impedance end, its current amplitude increases while the voltage amplitude drops, with a ratio that satisfies the change of impedance [Fig. 1(c)]. Assuming an ideal broadband transformer with perfect impedance matching and power transmission, the current leaving the low-impedance end (to the load) ΔIDL is related to the current injected to the high-impedance end ΔIDH by ΔIDL2ZL=ΔIDH2ZH. In our transmission line taper, this relation is valid only at high frequency (passband of the taper), which dominates the rising edge of the detector pulse. In the extreme case where the SNSPD pushes all the bias current out, i.e., ΔIDH = IB, the current diverted to the load can be as large as IL=ΔIDL=IBZHZL, corresponding to an output voltage of VLtaper=IB×50ΩZHZL and an effective voltage gain of ZHZL with respect to the conventional readout. In practice, when terminated with the high-impedance taper, the SNSPD latches, leaving the residual current at the hotspot current Iss, and then resets through reflection from the taper. As we will show later, the actual voltage gain is always less than ZHZL due to the electro-thermal feedback and the limited taper bandwidth.

Figure 1(d) shows an optical micrograph of a fabricated SNSPD with a meandered transmission-line taper. The bright area is NbN, and the violet area is the substrate, where the NbN was etched away. The NbN was sputtered at room temperature on a silicon substrate with a 300 nm thick thermal oxide layer.26 The film had a critical temperature of 8.1 K and a room-temperature sheet resistance of 342 Ω/sq. The sheet inductance was estimated to be 80 pH/sq by fitting the falling edge of the output pulse from a reference detector. The nanowire fabrication process is described in Refs. 4 and 27. The taper was made from a CPW with a fixed gap size of 3 μm and a varying center conductor width from 135 μm (50 Ω) to 500 nm (1.7 kΩ; see Fig. S1). Transmission lines made from thin-film superconductors operate in the kinetic-inductive limit and have been studied previously.28,29 The taper's left/wide end was wire bonded to an external circuit board, and the right/narrow end was connected to the SNSPD through a 1 μm-long hyperbolic taper. The short hyperbolic taper is not for r.f. impedance matching, but to allow width transition without causing d.c. current crowding. The SNSPD was 100 nm wide, densely packed with a 50% fill factor and spanned a rectangular area of 11 μm × 10 μm [see Fig. 1(e)]. A 200 nm gap surrounded the detector region to reduce the proximity effect in fabrication. On the same chip, we also fabricated non-tapered detectors as references.

The taper was designed to be a 5672-section cascaded transformer with a lower cut-off frequency of 200 MHz and a total electrical length of 851 mm, following the Klopfenstein taper profile.30 The physical length was 77.9 mm due to the slow phase velocity of the superconducting transmission line, and the total inductance was 1.410 μH (see Fig. S2 for the simulated S parameters). This length was chosen so that the maximum reflection in the passband (≥200 MHz) would not exceed −20 dB. The total electrical length is calculated as le = Ac/2πfco, where fco is the nominal cut-off frequency and c is the speed of light in vacuum. A is a factor that determines the maximum reflection in the passband and is calculated as cosh(A)=ρ0/ρpb, where ρ0 is the initial reflection coefficient (i.e., without taper) and ρpb is the maximum passband reflection coefficient (taken to be 0.1 here). For design convenience, we followed Klopfenstein's original approach and took ρ0 = 0.5ln(ZH/ZL) instead of (ZHZL)/(ZH + ZL).30 

The detectors were operated at 1.3 K in a closed-cycle cryostat. Both the bias circuit and the readout electronics were at room temperature. The output signal of the detectors was amplified using a 2.5 GHz, 25 dB gain low-noise amplifier (RF BAY LNA-2500), and a 3 dB attenuator was inserted before the amplifier to reduce reflection and prevent latching. The output pulses from the amplifier were then acquired using a 6 GHz real-time oscilloscope (Lecroy 760Zi). In fact, the output of the tapered SNSPD is large enough to operate without amplifiers. However, adding a room-temperature amplifier increases the SNR and reduces the timing jitter. The detector chip was flood illuminated using attenuated sub-ps pulsed lasers at 1550 nm (FPL-02CCF) through an optical fiber (SMF-28e). The laser pulses were split into two arms, one to a variable attenuator then to the cryostat and the other to a fast photodiode (Thorlabs DET08CFC) as timing references. Since the distance between the non-tapered detector and the tapered detector (≈5 mm) was much less than the distance between the detector chip and the fiber tip (≈10 cm), we expect the difference in photon arrival time to be <1 ps. Both the tapered and non-tapered detectors had a switching current of 30 μA and were biased at 27.5 μA throughout the measurement.

Figure 2(a) shows the measured pulse shapes from the reference and tapered detectors. The amplifier gain was removed to better compare with simulations. To avoid phase distortion in reconstructing the unamplified pulses, we used a weighted gain, G¯=dfPSD(f)G(f)/dfPSD(f), where PSD(f) is the power spectral density of the pulse, and G(f) is the measured system gain spectrum (see Figs. S3 and S4). G¯ was calculated to be 20.5 dB. As shown in Fig. 2(a), we observed a voltage gain of 3.6 and a delay of 2.8 ns from the tapered device compared to the reference device (by aligning the electrical pulses to the optical references). This voltage enhancement is equivalent to a passive, excessive-noise-free 11 dB amplifier.

FIG. 2.

Measured detector pulses and comparison with SPICE simulation. (a) Measured voltage with the amplifier gain removed. The black traces are single-shot waveforms, and the superimposed colored lines are averaged waveforms. (b) Simulated output voltages for both the tapered and non-tapered detectors. (c) Simulated current distributions. ID is the current in the SNSPD for a non-tapered readout. IDL and IDH are the currents on the low-impedance and high-impedance ends of the tapers in the tapered detector, respectively.

FIG. 2.

Measured detector pulses and comparison with SPICE simulation. (a) Measured voltage with the amplifier gain removed. The black traces are single-shot waveforms, and the superimposed colored lines are averaged waveforms. (b) Simulated output voltages for both the tapered and non-tapered detectors. (c) Simulated current distributions. ID is the current in the SNSPD for a non-tapered readout. IDL and IDH are the currents on the low-impedance and high-impedance ends of the tapers in the tapered detector, respectively.

Close modal

We simulated the tapered detector using a SPICE model that incorporates both the hotspot dynamics in the SNSPD and the distributed nature of the transmission line taper.31,32 The simulation was implemented in LTspice, a free electrical circuit modeling software. The SPICE model for SNSPD was developed by Berggren et al.31 based on the phenomenological hotspot velocity model by Kerman et al.12 The taper was simulated as cascaded lossless transmission lines (down-sampled to 300 sections),32 and each section was implemented using the Lossy Transmission Line Model (LTRA) in LTspice with different length, inductance, and capacitance settings.

Figure 2(b) plots the simulated load voltages, showing a voltage gain of 3.5 and a delay of 2.8 ns, as compared to the measured gain of 3.6 and delay of 2.8 ns. The subsequent peaks in the output voltage are spaced by ≈4.2 ns, which should correspond to the round trip time in the taper. The single-trip delay of the taper calculated from the reflection peaks (2.1 ns) is shorter than the delay between the tapered detector and the reference detector (2.82 ns) because the hotspot grows for a longer period of time and to a larger resistance in the tapered detector, as can be seen in the simulated currents.

Figure 2(c) shows the simulated currents. For the tapered detector, the current in the nanowire (IDH) first drops at a similar rate as the non-tapered case (ID) and then enters a intermediate plateau due to latching.12 Similar current plateau and latching behaviors are often observed when loading an SNSPD with a kΩ resistor. However, at ≈4.6 ns, IDL drops again and kicks the detector out of the latching state. The drop in current is from the reflection in the transmission line taper. Alternatively, it can be interpreted as the distributed capacitors in the transmission line drawing current from the SNSPD. After a few oscillations (high frequency), the current in the detector recovers with a τ = L/R exponential time constant (low frequency). Here, L is the total inductance of the SNSPD and the taper (at low frequency, the taper behaves as an inductor), and R is 50 Ω. The simulated currents at the high- and low-impedance ends of the taper roughly follow our intuitive understanding of how a transformer works. IDH drops from the original bias current 27.5 μA to ≈4.6 μA at the plateau (ΔIDH = −22.9 μA), while IDL drops from 27.5 μA to −54.7 μA (ΔIDL = −81.5 μA), giving ΔIDLDH ≈ 3.6. This is 39% less than the theoretical limit for a matched source and ideal transformer (1.7kΩ/50Ω=5.8), due to the complex electrothermal feedback and finite taper bandwidth. In this particular detector and taper design, the maximum counting rate (estimated as 1/3τ) decreases from 40.3 MHz (L =414 nH) for the non-tapered detector to 9.1 MHz (L =1.824 μH) for the tapered detector. We would like to point out that the added inductance and consequently the slower maximum counting rate are currently the major drawbacks of the tapered readout and may limit its usage for certain applications, such as high-data-rate communication. A study of the trade-off between the gain factor and taper inductance based on current CPW design can be found in Fig. S8 in the supplementary material. To reduce the taper inductance, one may use microstrip or CPWs with smaller gaps or with additional top/bottom ground. With larger distributed capacitances, these transmission lines will have shorter effective wavelengths, and the tapers will thus have smaller inductances.

The impedance taper amplifies the output pulse without sacrificing the fast rising edge, resulting in a faster slew rate. Figure 3(a) compares the averaged rising edges of the detector pulses from the reference and tapered detectors (with the amplifier gain). The sampling rate was 40 GS/s. As shown in Fig. 3(b), the maximum slew rates (dV/dt) were 39 μV/ps for the reference detector, but 143 μV/ps (3.7 times faster) for the tapered detector.

FIG. 3.

Experimental observation on the reduction of timing jitter as a result of the faster slew rate. (a) Averaged rising edges of the detector pulses from the tapered and non-tapered detectors (amplifier gain not removed); (b) corresponding slew rate calculated as dV/dt. (c) The measured FWHM timing jitter reduced from 48.9 ps to 23.8 ps with the tapered readout at 1550 nm.

FIG. 3.

Experimental observation on the reduction of timing jitter as a result of the faster slew rate. (a) Averaged rising edges of the detector pulses from the tapered and non-tapered detectors (amplifier gain not removed); (b) corresponding slew rate calculated as dV/dt. (c) The measured FWHM timing jitter reduced from 48.9 ps to 23.8 ps with the tapered readout at 1550 nm.

Close modal

The slew rate directly impacts the electrical noise contribution to the timing jitter, usually referred to as noise jitter, σnoise.11,33 We sampled the background electrical noise on the oscilloscope for both detectors by measuring the voltage at 400 ps before the rising edge of the pulses. The noise followed a Gaussian distribution and had standard deviations of 559 μV and 547 μV for the reference and tapered detector, respectively (see Fig. S5). Taking their respective fastest slew rates, we calculated that the reference detector would have a standard deviation σnoise of 14.3 ps and the tapered detector would have a σnoise value of 3.8 ps.

We measured the jitter of the detectors following the procedure described in Ref. 34. The discrimination levels for time tagging were set to voltages with the fastest slew rates. Figure 3(c) shows the instrument response function (IRF) of the reference and tapered detectors at an illumination wavelength of 1550 nm. With the impedance taper, the full-width half-maximum (FWHM) jitter reduced from 48.9 ps to 23.8 ps. We fitted the IRF using an exponentially modified Gaussian distribution4 and found σ = 16.8 ps and 1/λ = 17.4 ps for the reference detector and σ = 6.5 ps and 1/λ = 13.6 ps for the tapered detector. Here, σ is the standard deviation of the normal distribution and λ is the exponential decay rate. The detectors showed similar jitter reduction at 1064 nm, where both detectors operated on the saturation plateau (see Fig. S6). The FWHM jitter reduced from 47.0 ps (σ = 16.4 ps, 1/λ = 15.9 ps) to 22.4 ps (σ = 6.2 ps, 1/λ = 12.5 ps). We observed a leading edge tail in IRF for the tapered detector. It is likely due to the counting events from the taper or the transition region between the taper and the detector. This effect could be reduced if the optical mode is focused in the center of the detector active area, through self-aligned fiber coupling or focusing lenses.35,36

As a final remark, we have treated the SNSPD as a lumped element in this paper because the nanowire was closely meandered and had a dispersion similar to an ideal inductor at the frequency of interest.37 Despite this choice, multi-photon absorption would generate a different hotspot resistance compared to the single-photon events.22,38 The impedance taper provides an effective kΩ load and may thus allow direct discrimination of hotspot resistance and hence photon numbers. In another scheme, where the nanowire is sparse or designed to a transmission line,32 the taper provides impedance-matched readout and has been used to resolve the photon location and photon numbers.27,39 We expect the integrated taper to become a widely used tool for matching high-impedance nanowire-based devices to low-impedance systems.

See supplementary material for more information on the taper design, device measurement, and discussion on trade-off between the gain factor and taper inductance.

Part of this work was performed at the Jet Propulsion Laboratory, California Institute of Technology, under Contract with the National Aeronautics and Space Administration. Support for this work was provided in part by the JPL Strategic University Research Partnerships program, DARPA Defense Sciences Office through the DETECT program, and National Science Foundation under Contract Nos. ECCS-1509486 (MIT) and ECCS-1509253 (UNF). D.Z. is supported by the National Science Scholarship from A*STAR, Singapore.

1.
R.
Hadfield
,
Nat. Photonics
3
,
696
(
2009
).
2.
C. M.
Natarajan
,
M. G.
Tanner
, and
R. H.
Hadfield
,
Supercond. Sci. Technol.
25
,
063001
(
2012
).
3.
F.
Marsili
,
V. B.
Verma
,
J. A.
Stern
,
S.
Harrington
,
A. E.
Lita
,
T.
Gerrits
,
I.
Vayshenker
,
B.
Baek
,
M. D.
Shaw
,
R. P.
Mirin
, and
S. W.
Nam
,
Nat. Photonics
7
,
210
(
2013
).
4.
B. A.
Korzh
,
Q.-Y.
Zhao
,
S.
Frasca
,
J. P.
Allmaras
,
T. M.
Autry
,
E. A.
Bersin
,
M.
Colangelo
,
G. M.
Crouch
,
A. E.
Dane
,
T.
Gerrits
,
F.
Marsili
,
G.
Moody
,
E.
Ramirez
,
J. D.
Rezac
,
M. J.
Stevens
,
E. E.
Wollman
,
D.
Zhu
,
P. D.
Hale
,
K. L.
Silverman
,
R. P.
Mirin
,
S. W.
Nam
,
M. D.
Shaw
, and
K. K.
Berggren
, preprint arXiv:1804.06839 (
2018
).
5.
E. E.
Wollman
,
V. B.
Verma
,
A. D.
Beyer
,
R. M.
Briggs
,
B.
Korzh
,
J. P.
Allmaras
,
F.
Marsili
,
A. E.
Lita
,
R. P.
Mirin
,
S. W.
Nam
, and
M. D.
Shaw
,
Opt. Express
25
,
26792
(
2017
).
6.
A.
Vetter
,
S.
Ferrari
,
P.
Rath
,
R.
Alaee
,
O.
Kahl
,
V.
Kovalyuk
,
S.
Diewald
,
G. N.
Goltsman
,
A.
Korneev
,
C.
Rockstuhl
, and
W. H. P.
Pernice
,
Nano Lett.
16
,
7085
(
2016
).
7.
F.
Marsili
,
F.
Bellei
,
F.
Najafi
,
A. E.
Dane
,
E. A.
Dauler
,
R. J.
Molnar
, and
K. K.
Berggren
,
Nano Lett.
12
,
4799
(
2012
).
8.
H.-L.
Yin
,
T.-Y.
Chen
,
Z.-W.
Yu
,
H.
Liu
,
L.-X.
You
,
Y.-H.
Zhou
,
S.-J.
Chen
,
Y.
Mao
,
M.-Q.
Huang
,
W.-J.
Zhang
,
H.
Chen
,
M. J.
Li
,
D.
Nolan
,
F.
Zhou
,
X.
Jiang
,
Z.
Wang
,
Q.
Zhang
,
X.-B.
Wang
, and
J.-W.
Pan
,
Phys. Rev. Lett.
117
,
190501
(
2016
).
9.
X.
Qiang
,
X.
Zhou
,
J.
Wang
,
C. M.
Wilkes
,
T.
Loke
,
S.
O'Gara
,
L.
Kling
,
G. D.
Marshall
,
R.
Santagati
,
T. C.
Ralph
,
J. B.
Wang
,
J. L.
O'Brien
,
M. G.
Thompson
, and
J. C. F.
Matthews
,
Nat. Photonics
12
,
534
(
2018
).
10.
M. E.
Grein
,
A. J.
Kerman
,
E. A.
Dauler
,
M. M.
Willis
,
B.
Romkey
,
R. J.
Molnar
,
B. S.
Robinson
,
D. V.
Murphy
, and
D. M.
Boroson
,
Proc. SPIE
9492
,
949208
(
2015
).
11.
J.
Wu
,
L.
You
,
S.
Chen
,
H.
Li
,
Y.
He
,
C.
Lv
,
Z.
Wang
, and
X.
Xie
,
Appl. Opt.
56
,
2195
(
2017
).
12.
A. J.
Kerman
,
J. K. W.
Yang
,
R. J.
Molnar
,
E. A.
Dauler
, and
K. K.
Berggren
,
Phys. Rev. B
79
,
100509
(
2009
).
13.
M.
Ejrnaes
,
R.
Cristiano
,
O.
Quaranta
,
S.
Pagano
,
A.
Gaggero
,
F.
Mattioli
,
R.
Leoni
,
B.
Voronov
, and
G.
Gol'tsman
,
Appl. Phys. Lett.
91
,
262509
(
2007
).
14.
Y.
Cheng
,
H.
Liu
,
C.
Gu
,
X.
Zhu
, and
X.
Hu
, in
Conference Lasers Electro-Optics
(OSA, Washington, D.C.,
2017
), p.
JW2A.120
.
15.
J. C.
Bardin
,
P.
Ravindran
,
S.-W.
Chang
,
C.
Mohamed
,
R.
Kumar
,
J. A.
Stern
,
M. D.
Shaw
,
D.
Russell
,
F.
Marsili
,
G.
Resta
, and
W. H.
Farr
,
Proc. SPIE
9114
,
911404
(
2014
).
16.
H.
Mani
and
P.
Mauskopf
, in
2014 11th International Workshop on Low Temperature Electronics
(IEEE,
2014
), pp.
17
20
.
17.
C.
Wan
,
Z.
Jiang
,
L.
Kang
, and
P.
Wu
,
IEEE Trans. Appl. Supercond.
27
,
1502605
(
2017
).
18.
C.
Cahall
,
D. J.
Gauthier
, and
J.
Kim
,
Rev. Sci. Instrum.
89
,
063117
(
2018
).
19.
A. N.
McCaughan
and
K. K.
Berggren
,
Nano Lett.
14
,
5748
(
2014
).
20.
H.
Terai
,
S.
Miki
,
T.
Yamashita
,
K.
Makise
, and
Z.
Wang
,
Appl. Phys. Lett.
97
,
112510
(
2010
).
21.
T.
Ortlepp
,
M.
Hofherr
,
L.
Fritzsch
,
S.
Engert
,
K.
Ilin
,
D.
Rall
,
H.
Toepfer
,
H.-G.
Meyer
, and
M.
Siegel
,
Opt. Express
19
,
18593
(
2011
).
22.
M.
Bell
,
A.
Antipov
,
B.
Karasik
,
A.
Sergeev
,
V.
Mitin
, and
A.
Verevkin
,
IEEE Trans. Appl. Supercond.
17
,
289
(
2007
).
23.
J.
Kitaygorsky
,
S.
Dorenbos
,
E.
Reiger
,
R.
Schouten
,
V.
Zwiller
, and
R.
Sobolewski
,
IEEE Trans. Appl. Supercond.
19
,
346
(
2009
).
24.
S.
Jahanmirinejad
and
A.
Fiore
,
Opt. Express
20
,
5017
(
2012
).
25.
J. K. W.
Yang
,
A. J.
Kerman
,
E. A.
Dauler
,
V.
Anant
,
K. M.
Rosfjord
, and
K. K.
Berggren
, “
IEEE Trans
,”
Appl. Supercond.
17
,
581
(
2007
).
26.
A. E.
Dane
,
A. N.
McCaughan
,
D.
Zhu
,
Q.
Zhao
,
C.-S.
Kim
,
N.
Calandri
,
A.
Agarwal
,
F.
Bellei
, and
K. K.
Berggren
,
Appl. Phys. Lett.
111
,
122601
(
2017
).
27.
Q.-Y.
Zhao
,
D.
Zhu
,
N.
Calandri
,
A. E.
Dane
,
A. N.
McCaughan
,
F.
Bellei
,
H.-Z.
Wang
,
D. F.
Santavicca
, and
K. K.
Berggren
,
Nat. Photonics
11
,
247
(
2017
).
28.
J. C.
Swihart
,
J. Appl. Phys.
32
,
461
(
1961
).
29.
K.
Yoshida
,
M. S.
Hossain
,
T.
Kisu
,
K.
Enpuku
, and
K.
Yamafuji
,
Jpn. J. Appl. Phys., Part 1
31
,
3844
(
1992
).
30.
R. W.
Klopfenstein
,
Proc. IRE
44
,
31
(
1956
).
31.
K. K.
Berggren
,
Q.-Y.
Zhao
,
N.
Abebe
,
M.
Chen
,
P.
Ravindran
,
A.
McCaughan
, and
J. C.
Bardin
,
Supercond. Sci. Technol.
31
,
055010
(
2018
).
32.
Q.-Y.
Zhao
,
D. F.
Santavicca
,
D.
Zhu
,
B.
Noble
, and
K. K.
Berggren
,
Appl. Phys. Lett.
113
,
082601
(
2018
).
33.
Q.
Zhao
,
L.
Zhang
,
T.
Jia
,
L.
Kang
,
W.
Xu
,
J.
Chen
, and
P.
Wu
,
Appl. Phys. B
104
,
673
(
2011
).
34.
F.
Najafi
,
F.
Marsili
,
E.
Dauler
,
R. J.
Molnar
, and
K. K.
Berggren
,
Appl. Phys. Lett.
100
,
152602
(
2012
).
35.
A. J.
Miller
,
A. E.
Lita
,
B.
Calkins
,
I.
Vayshenker
,
S. M.
Gruber
, and
S. W.
Nam
,
Opt. Express
19
,
9102
(
2011
).
36.
F.
Bellei
,
A. P.
Cartwright
,
A. N.
McCaughan
,
A. E.
Dane
,
F.
Najafi
,
Q.
Zhao
, and
K. K.
Berggren
,
Opt. Express
24
,
3248
(
2016
).
37.
D. F.
Santavicca
,
J. K.
Adams
,
L. E.
Grant
,
A. N.
Mccaughan
, and
K. K.
Berggren
,
J. Appl. Phys.
119
,
234302
(
2016
).
38.
C.
Cahall
,
K. L.
Nicolich
,
N. T.
Islam
,
G. P.
Lafyatis
,
A. J.
Miller
,
D. J.
Gauthier
, and
J.
Kim
,
Optica
4
,
1534
(
2017
).
39.
D.
Zhu
,
Q.-Y.
Zhao
,
H.
Choi
,
T.-J.
Lu
,
A. E.
Dane
,
D.
Englund
, and
K. K.
Berggren
,
Nat. Nanotechnol.
13
,
596
(
2018
).

Supplementary Material