Recently, the possibility to tune the critical current of conventional metallic superconductors via electrostatic gating was shown in wires, Josephson weak-links, and superconductor-normal metal–superconductor junctions. Here, we exploit such a technique to demonstrate a gate-controlled vanadium-based Dayem nano-bridge operated as a half-wave rectifier at 3 K. Our devices exploit the gate-driven modulation of the critical current of the Josephson junction and the resulting steep variation of its normal-state resistance, to convert an AC signal applied to the gate electrode into a DC one across the junction. All-metallic superconducting gated rectifiers could provide the enabling technology to realize tunable photon detectors and diodes useful for superconducting electronics circuitry.

In the last few years, the possibility to suppress the critical current in conventional Bardeen–Cooper–Schrieffer (BCS) superconducting wires,1 Dayem bridge (DB) Josephson junctions (JJs),2–5 and superconductor-normal metal–superconductor weak-links6 via electrostatic gating was demonstrated.7 Such a technique is at the basis of a class of innovative all-metallic superconducting Josephson transistors, where a gate electrode is used to modulate down to full suppression the critical current of a superconducting channel, in analogy to what occurs for the current in standard gated semiconductor transistors. Furthermore, the impact of the electrostatic gating on the macroscopic phase and on the phase slip dynamics was proved in superconducting quantum interference devices (SQUIDs)8 and in gated titanium DBs, respectively.9 So far, the microscopic origin of the gating effect in metallic superconductors is unclear, and the debate in the search of a satisfactory explanation for the aforementioned phenomenology has just started.4,10,11 Electrostatic-field-driven surface-orbital polarization was proposed as the knob to control both the amplitude and the phase of gated superconductors.4,11 Furthermore, it was conjectured about the role of a tiny field-emission gate-superconductor current.12,13 Yet, several technological implementations based on gated superconducting transistors for both classical and quantum computation architectures were proposed7 until now and promise new paradigms to realize innovative superconducting circuits. These are expected to benefit from operation frequencies close to those typical of rapid single flux quantum devices (up to 750 GHz),14 but in a voltage-driven fashion, which would provide a direct and convenient interface with complementary metal-oxide technology (CMOS).

Here, we demonstrate a vanadium (V) Josephson half-wave rectifier based on a scheme recently proposed,5 which exploits the strong non-linearity of the resistance (R) of a gated superconducting nanotransistor vs gate voltage (VG). Such a device might form the rectifying core of a photon sensor operating at high frequencies, relevant, for instance, for cosmological applications15 or could be exploited as a comparator16 to be used in all-superconducting digital electronic circuits.7 The potentially high operation frequency, the ease of fabrication, and the simple geometry of our nanodevices make them attractive as effective alternatives to conventional all-semiconductor17 and semiconductor–superconductor hybrid18–23 field-effect-transistor-based technologies.24–28 Below, we focus on the manipulation of the dissipative response of the V-based Josephson weak-link via the application of a voltage to a gate capacitively coupled to the junction and show its time-resolved response to periodic square-wave and sinusoidal signals applied to the gate, thereby leading to electric rectification.

Our device consists of a planar 60-nm-tick, 160-nm-long, 90-nm-wide V Dayem bridge JJ with a 70-nm-apart, 120-nm-wide side-gate aligned to the constriction. The device fabrication was performed on a silicon oxide substrate by single-step electron beam lithography followed by vanadium deposition performed at a rate of 0.36 nm/s in an ultra-high vacuum electron-beam evaporator with a base pressure of 1011 Torr.29–35Figure 1(a) shows the false color scanning electron micrograph of a representative V JJ transistor along with the 4-wire biasing scheme used for the low-temperature current vs voltage (IV) characterization of the devices performed in a filtered cryogen-free 3He–4He dilution refrigerator.

FIG. 1.

(a) False-color scanning electron micrograph of a typical gate-controlled vanadium JJ nanotransistor. The DB constriction (light blue) is current biased, and the voltage drop is measured in a 4-wire configuration. The electrostatic field is applied via the gate electrode (orange). (b) Back and forth current I vs voltage V characteristics of a representative device measured at different bath temperatures from 2.0 K to 3.8 K. The curves are horizontally offset for clarity. The inset shows the 4-wire lock-in measurement (with biasing current I15 nA) of the R vs T characteristic of the junction. The superconducting transition occurs at TC3.6K. (c) Evolution of the critical current IC (light blue dots) and retrapping current IR (yellow dots) vs bath temperature T. The black dashed line represents the best fit of IC(T) with Bardeen's theory, which yields IC0=(2.2±0.1) mA and TCfit=(3.62±0.07) K as fitting parameters.

FIG. 1.

(a) False-color scanning electron micrograph of a typical gate-controlled vanadium JJ nanotransistor. The DB constriction (light blue) is current biased, and the voltage drop is measured in a 4-wire configuration. The electrostatic field is applied via the gate electrode (orange). (b) Back and forth current I vs voltage V characteristics of a representative device measured at different bath temperatures from 2.0 K to 3.8 K. The curves are horizontally offset for clarity. The inset shows the 4-wire lock-in measurement (with biasing current I15 nA) of the R vs T characteristic of the junction. The superconducting transition occurs at TC3.6K. (c) Evolution of the critical current IC (light blue dots) and retrapping current IR (yellow dots) vs bath temperature T. The black dashed line represents the best fit of IC(T) with Bardeen's theory, which yields IC0=(2.2±0.1) mA and TCfit=(3.62±0.07) K as fitting parameters.

Close modal

For temperatures below 3.5 K, the V JJ is superconducting and exhibits, at 2 K, a critical current IC of 1.42 mA, a retrapping current IR380μA, and a normal-state resistance RN of around 110Ω. The forward and backward current I vs V characteristics for temperatures T ranging from 2.0 K to 3.8 K are shown in Fig. 1(b), where the full evolution of IC as a function of T is highlighted by the gray area. The hysteretic behavior of the curves stems from quasiparticle over-heating when the junction switches from the superconducting to the normal state.36,37 Moreover, the additional features present on the IV characteristics at higher current are likely to be related to the switching of the DB banks to the normal state. The decay of both IC and IR with the temperature is shown in Fig. 1(c), where the black dashed line represents a fit of the IC vs T characteristic performed with Bardeen's equation: IC(T)=IC0[1(TTC)2]32, where IC0 is the zero-temperature critical current and TC is the DB critical temperature. The fitting procedure yielded IC0=(2.2±0.1) mA and TCfit=(3.62±0.07) K. The latter value is in agreement with the experimental critical temperature determined from the 4-wire lock-in measurement of the junction resistance vs bath temperature T, as shown in the inset of Fig. 1(b).

The preliminary characterization of the gating effect on IC and R was carried out by acquiring the forward and backward IV characteristics of a representative vanadium DB as a function of the voltage applied to the gate electrode, VG. Figure 2(a) displays the IC(VG) characteristics extracted from IV acquired at several bath temperatures. The suppression of IC occurs for both positive and negative values of VG (VGVG), and the gating effect persists up to 3.3 K, i. e., 0.9TC. We note that the sharper suppression of IC observed for positive values of VG seems to be in contrast to a possible contribution to the critical current suppression driven by cold field-emission of electrons. Indeed, electron extraction from the gate (occurring at negative gate bias values) is expected to be more effective than that from the DB (occurring at positive gate bias values) owing to the substantial asymmetry of the device geometry.13 

FIG. 2.

(a) IC vs VG characteristics for different temperatures ranging from 2.0 to 3.3 K. Data were computed averaging 25 acquisitions of the switching current. The error bars represent the standard deviation of the ensemble of measurements. (b) Resistance R=V/IB vs gate voltage VG for different values of bias current IB from 1.4 to 150 μA. (c) Transresistance dRdVG vs gate voltage VG for the same bias currents of panel (b). Measurements of panels (b) and (c) were carried out at TB=3.0 K.

FIG. 2.

(a) IC vs VG characteristics for different temperatures ranging from 2.0 to 3.3 K. Data were computed averaging 25 acquisitions of the switching current. The error bars represent the standard deviation of the ensemble of measurements. (b) Resistance R=V/IB vs gate voltage VG for different values of bias current IB from 1.4 to 150 μA. (c) Transresistance dRdVG vs gate voltage VG for the same bias currents of panel (b). Measurements of panels (b) and (c) were carried out at TB=3.0 K.

Close modal

Let us now describe further the transport properties of the transistor by focusing on the behavior of the device resistance R under the action of the electric field. The steep modulation of the critical current yields, at constant current bias IB, an analogous modulation of the resistance of the device. Figure 2(b) displays the evolution of R as a function of gate voltage VG for several values of bias current IB, which was always set below the retrapping current so that no hysteretic behavior was expected from the device. We note that R remains almost equal to 0 for VG values such that IC(VG)>IB. Then, for higher gate voltage, the resistance jumps to finite values due to the gate-driven superconducting-to-normal state transition of the device. Finally, R increases with VG, probably owing to the enlargement of the normal-state region within the junction.2 In particular, for IB>30μA, R exhibits a linear dependence on VG. Moreover, the steepness of the resistance jump also depends on IB, as customarily described by the figure of merit of the resistance derivative with respect to the gate voltage (dR/dVG), which is shown in Fig. 2(c). The amplitude of the differential resistance peaks, corresponding to the transition to the normal state, increases with IB reaching the maximum value of 30Ω/ V at IB = 150 μA. Furthermore, the peak center shifts toward lower VG values as the bias current is increased.

The strong non-linearity typical of the R(VG) characteristics, provided by the transition to the dissipative regime, could be conveniently exploited to implement a superconducting half-wave rectifier converting an AC signal applied to the gate into a DC one across the DB. Such a mechanism has been recently proposed for Nb gate-controlled transistors,5 where the rectification properties of a DB device have been analyzed by a standard lock-in technique, which does not provide information neither on the time-resolved behavior of the devices nor on the shape of the output signal. These features could allow us to realize, for instance, a radiation detector in which the AC electric field collected by an antenna coupled to the gate electrode would be rectified, thereby leading to an amplified non-zero average signal. The gain of such a device can be defined as the ratio g=Vout/Vin, where Vin and Vout are the peak-to-peak amplitudes of the signal applied to the gate and the voltage drop across the DB, respectively. The quantity g can be calculated through the relation g=R(VG)IBτ/σ, where R(VG)IB=Vout is the product of the gate-dependent DB resistance R and the current bias IB, while σ/τ=Vin is given by the ratio of the typical width (σ) of the switching current probability distribution for gate-driven superconducting field-effect transistors9 and the transconductance τ=dIC/dVG.7 For V-based devices, g is expected to reach values as high as 7 with a typical power dissipation of 40 nW. Our V DB nanotransistor shows gain performance similar to cryogenic semiconductor devices38,39 with a reduction in power consumption of about three orders of magnitude. Noteworthily, a series of N rectifiers can be realized by feeding the gate electrode of the nth rectifier with the output voltage of the (n1)th one, resulting in a total gain equal to gN.

In the following, we shall focus on the time-resolved investigation of the response of V Josephson transistors to AC square-wave and sinusoidal gate signals. The measurement setup, comprising an input/output analog-to-digital/digital-to-analog converter (ADC/DAC) board for the generation and the acquisition of time-resolved alternate current and voltage signals, is depicted in Fig. 3(a). The voltage drop V(t) across the current-biased DB was measured as a function of a time-dependent gate voltage VG(t). In this configuration, when IC[VG(t)]<IB, the junction is in the normal state and a finite voltage drop is built across the DB (high-state). By contrast, when IC(VG(t))>IB, the DB is superconducting and V0 (low-state). The response of the device to a transistor–transistor-logic (TTL)-like square-wave excitation was probed by feeding the gate with a signal consisting of a VDC = 10 V DC bias added to a VAC = 5 V square-wave signal with frequencies up to 50 Hz, as shown in Fig. 3(c). The voltage drop in the high and low states for some representative values of IB is shown in Fig. 3(d) as pairs of points of the same color superimposed on the IV characteristics. Such curves are taken at selected values of VG. As displayed in Fig. 3(d), the voltage drop across the junction preserves the shape of the input signal. In this regard, we emphasize that in our experiments, the maximum operation frequency is limited by the low-pass filtering stages of the electrical lines of our setup. A lower bound to the operation frequency of the V gate-controllable rectifiers might arise in the 1081010 Hz range due to the typical electron–phonon relaxation rates at these bath temperatures.40 Furthermore, the frequency defined by 2Δ/h265 GHz, where Δ is the superconducting gap energy and h Planck's constant, is likely to provide the ultimate operation limit.41 The latter value is comparable to state-of-the-art hybrid and CMOS field-effect transistors;17,42,43 nonetheless, we stress that at the present time, the timescale of gate-driven superconducting phase transition is still totally unknown and is worth to be investigated through specifically designed experiments.

FIG. 3.

(a) Scheme of time-resolved measurements of the 4-wire voltage drop V(t) across the junction at constant current bias IB vs DC (black dotted line) + AC (square-wave green line) gate voltage. The gate polarization signal is compared with the IC vs VG curve to highlight the dynamic working range of the system. (b) Voltage V vs current I characteristics for selected values of gate voltage VG. The dot pairs show the working points of the system for different values of the bias current IB=18,37,71μA. (c) Evolution of the gate voltage vs time. The signal was obtained by summing a DC voltage VDC = 10 V and an AC square-wave voltage with amplitude VAC = 5 V. (d) Time-resolved voltage drop across the DB for selected values of the current bias. Dashed lines show the corresponding working points in panel (b). These measurements were performed at TB = 3 K.

FIG. 3.

(a) Scheme of time-resolved measurements of the 4-wire voltage drop V(t) across the junction at constant current bias IB vs DC (black dotted line) + AC (square-wave green line) gate voltage. The gate polarization signal is compared with the IC vs VG curve to highlight the dynamic working range of the system. (b) Voltage V vs current I characteristics for selected values of gate voltage VG. The dot pairs show the working points of the system for different values of the bias current IB=18,37,71μA. (c) Evolution of the gate voltage vs time. The signal was obtained by summing a DC voltage VDC = 10 V and an AC square-wave voltage with amplitude VAC = 5 V. (d) Time-resolved voltage drop across the DB for selected values of the current bias. Dashed lines show the corresponding working points in panel (b). These measurements were performed at TB = 3 K.

Close modal

The response of the DB to a sinusoidal gate voltage excitation was measured by applying to the gate an AC sine signal VAC with peak-to-peak amplitudes ranging from 1 to 3.5 V added to a VDC = 11 V DC voltage bias. This results in a periodic sinusoidal oscillation of the critical current IC above and below IB, producing a gate-driven periodic modulation of the device resistance. The voltage drop across the DB as a function of IB and VG is displayed in Fig. 4(a). We note that by choosing IB and VDC, it is possible to set a working point for the rectifier, which determines the intensity and the shape of the resulting output voltage signal through the evolution of the super-to-normal switching point [red dashed line in Fig. 4(a)] and through the dependence of the device resistance on the gate voltage (also see Ref. 2). Indeed, the system lies in the zero voltage state for IC(VG(t))>IB, and a voltage drop across the junction occurs when IC(VG(t))IB. For higher values of VG, i.e., when IC(VG(t))<IB, the voltage drop increases due to the gate-driven evolution of the DB resistance, which eventually will saturate at the asymptotic value of the normal-state resistance. Figure 4(c) shows the voltage drop of the DB as a function of time, in response to the gate signal reported in Fig. 4(b). To clarify the response of the system to such excitation, we highlight in Figs. 4(a) and 4(b) three points corresponding to the minimum of the VG(t) signal (light green dot), the switching point (red dot), and the maximum of VG(t) (dark green dot). The value of the biasing current was set to IB=72μA in order to take advantage of both the sharp normal-to-superconductor transition and the linear relation between R and VG for IB > IC. The latter feature allows us to preserve the sinusoidal shape in the output voltage when the JJ is in the resistive state, as shown in Figs. 4(c)–4(e) for different values of VAC. Indeed, although the output signal resembles that of a conventional half-wave rectifier,43 it is worth emphasizing that, in our systems, it is possible to tune the rectification threshold VG* by changing both IB and VDC.

FIG. 4.

(a). Density plot of the voltage drop across the junction as a function of the gate voltage (x axis) and current bias (y axis). From left to right, the three points represent the minimum zero resistance gate voltage value (light green), the transition between the superconducting and the normal state (red), and the maximum of both gate voltage and voltage drop (dark green). The dashed red curve highlights the critical current IC as a function of the gate voltage VG. (b) Time evolution of the gate voltage. The signal was obtained adding a DC voltage VDC = 11 V and an AC sine wave voltage VAC. The three dots highlight the working points described in panel (a). [(c)–(e)] Time-resolved voltage drop across the DB for VAC=3.5,1.5,1.0 V, respectively. The colormap is the same as panel (a). All these measurements were carried out at TB = 3 K.

FIG. 4.

(a). Density plot of the voltage drop across the junction as a function of the gate voltage (x axis) and current bias (y axis). From left to right, the three points represent the minimum zero resistance gate voltage value (light green), the transition between the superconducting and the normal state (red), and the maximum of both gate voltage and voltage drop (dark green). The dashed red curve highlights the critical current IC as a function of the gate voltage VG. (b) Time evolution of the gate voltage. The signal was obtained adding a DC voltage VDC = 11 V and an AC sine wave voltage VAC. The three dots highlight the working points described in panel (a). [(c)–(e)] Time-resolved voltage drop across the DB for VAC=3.5,1.5,1.0 V, respectively. The colormap is the same as panel (a). All these measurements were carried out at TB = 3 K.

Close modal

In summary, we have demonstrated a half-wave rectifier based on a gated vanadium Dayem nano-bridge Josephson transistor operating at 3 K. We showed the time-resolved response of our device to square-wave and sinusoidal gate signals, demonstrating the ability to convert an AC excitation into a DC one, and with a gain, which is expected to obtain a value close to 7 at a power consumption of 40 nW. The performance of the V DB represents a major improvement in energy efficiency compared to conventional cryogenic semiconductor technology.38,39 V-based gate-controllable Josephson rectifiers could represent a breakthrough and a possible enabling technology to implement innovative all-metallic superconducting photon detectors44–46 based on gating as well as electric diodes suitable for superconducting digital41 and quantum electronics47,48 circuitry.

The authors acknowledge the European Union's Horizon 2020 research and innovation program under Grant No. 777222 ATTRACT (ProjectT-CONVERSE) and under Grant No. 800923-SUPERTED for partial financial support.

The data that support the findings of this study are available from the corresponding author upon reasonable request.

1.
G.
De Simoni
,
F.
Paolucci
,
P.
Solinas
,
E.
Strambini
, and
F.
Giazotto
,
Nat. Nanotechnol.
13
,
802
(
2018
).
2.
F.
Paolucci
,
G.
De Simoni
,
E.
Strambini
,
P.
Solinas
, and
F.
Giazotto
,
Nano Lett.
18
,
4195
(
2018
).
3.
F.
Paolucci
,
G.
De Simoni
,
P.
Solinas
,
E.
Strambini
,
N.
Ligato
,
P.
Virtanen
,
A.
Braggio
, and
F.
Giazotto
,
Phys. Rev. Appl.
11
,
024061
(
2019
).
4.
L.
Bours
,
M. T.
Mercaldo
,
M.
Cuoco
,
E.
Strambini
, and
F.
Giazotto
, arXiv:2003.01655 (
2020
).
5.
G.
De Simoni
,
C.
Puglia
, and
F.
Giazotto
, arXiv:2004.09182 (
2020
).
6.
G.
De Simoni
,
F.
Paolucci
,
C.
Puglia
, and
F.
Giazotto
,
ACS Nano
13
,
7871
(
2019
).
7.
F.
Paolucci
,
G.
De Simoni
,
P.
Solinas
,
E.
Strambini
,
C.
Puglia
,
N.
Ligato
, and
F.
Giazotto
,
AVS Quantum Sci.
1
,
016501
(
2019
).
8.
F.
Paolucci
,
F.
Vischi
,
G.
De Simoni
,
C.
Guarcello
,
P.
Solinas
, and
F.
Giazotto
,
Nano Lett.
19
,
6263
(
2019
).
9.
C.
Puglia
,
G.
De Simoni
, and
F.
Giazotto
,
Phys. Rev. Appl.
13
,
054026
(
2020
).
10.
P.
Virtanen
,
A.
Braggio
, and
F.
Giazotto
,
Phys. Rev. B
100
,
224506
(
2019
).
11.
M. T.
Mercaldo
,
P.
Solinas
,
F.
Giazotto
, and
M.
Cuoco
, arXiv:1907.09227 (
2019
).
12.
L. D.
Alegria
,
C. G.
Bøttcher
,
A. K.
Saydjari
,
A. T.
Pierce
,
S. H.
Lee
,
S. P.
Harvey
,
U.
Vool
, and
A.
Yacoby
, arXiv:2005.00584 (
2020
).
13.
M. F.
Ritter
,
A.
Fuhrer
,
D. Z.
Haxell
,
S.
Hart
,
P.
Gumann
,
H.
Riel
, and
F.
Nichele
, arXiv:2005.00462 (
2020
).
14.
W.
Chen
,
A. V.
Rylyakov
,
V.
Patel
,
J. E.
Lukens
, and
K. K.
Likharev
,
Appl. Phys. Lett.
73
,
2817
(
1998
).
15.
16.
D. N.
Mukherjee
,
S.
Panda
, and
B.
Maji
,
IETE J. Res.
64
,
422
(
2018
).
17.
P.
Valizadeh
,
Field Effect Transistors, a Comprehensive Overview: From Basic Concepts to Novel Technologies
(
John Wiley Sons
,
2016
).
18.
T.
Nishino
,
M.
Hatano
,
H.
Hasegawa
,
F.
Murai
,
T.
Kure
,
A.
Hiraiwa
,
K.
Yagi
, and
U.
Kawabe
,
IEEE Electron Device Lett.
10
,
61
(
1989
).
19.
T. D.
Clark
,
R. J.
Prance
, and
A. D. C.
Grassie
,
J. Appl. Phys.
51
,
2736
(
1980
).
20.
K. E.
Gray
,
Appl. Phys. Lett.
32
,
392
(
1978
).
21.
M.
Okamoto
,
Y.
Ohkura
,
K.
Katayama
, and
T.
Nishino
,
Appl. Phys. Lett.
61
,
2473
(
1992
).
22.
F.
Wen
,
J.
Shabani
, and
E.
Tutuc
,
IEEE Trans. Electron Devices
66
,
5367
(
2019
).
23.
M.
Thomas
,
H.-R.
Blank
,
K. C.
Wong
,
H.
Kroemer
, and
E.
Hu
,
Phys. Rev. B
58
,
11676
(
1998
).
24.
Rectifier Circuits
,” edited by
F. W.
William
(
Macmillan Education
,
London
,
1972
).
25.
CMOS Circuit Design for RF Sensors
,” edited by
G.
Gudnason
and
E.
Bruun
(
Kluwer Academic Publishers
,
Boston
,
2002
).
26.
G.
Kawaguchi
,
A. A.
Bardin
,
M.
Suda
,
M.
Uruichi
, and
H. M.
Yamamoto
,
Adv. Mater.
31
,
1805715
(
2019
).
27.
N. A.
Mortensen
and
G.
Bastian
,
Superlattices Microstruct.
28
,
231
(
2000
).
28.
A. F.
Volkov
and
H.
Takayanagi
,
Phys. Rev. B
53
,
15162
(
1996
).
29.
C. P.
García
and
F.
Giazotto
,
Appl. Phys. Lett.
94
,
132508
(
2009
).
30.
P.
Spathis
,
S.
Biswas
,
S.
Roddaro
,
L.
Sorba
,
F.
Giazotto
, and
F.
Beltram
,
Nanotechnology
22
,
105201
(
2011
).
31.
O.
Quaranta
,
P.
Spathis
,
F.
Beltram
, and
F.
Giazotto
,
Appl. Phys. Lett.
98
,
032501
(
2011
).
32.
F.
Giazotto
,
P.
Spathis
,
S.
Roddaro
,
S.
Biswas
,
F.
Taddei
,
M.
Governale
, and
L.
Sorba
,
Nat. Phys.
7
,
857
(
2011
).
33.
A.
Ronzani
,
M.
Baillergeau
,
C.
Altimiras
, and
F.
Giazotto
,
Appl. Phys. Lett.
103
,
052603
(
2013
).
34.
A.
Ronzani
,
C.
Altimiras
, and
F.
Giazotto
,
Appl. Phys. Lett.
104
,
032601
(
2014
).
35.
N.
Ligato
,
G.
Marchegiani
,
P.
Virtanen
,
E.
Strambini
, and
F.
Giazotto
,
Sci. Rep.
7
,
8810
(
2017
).
36.
W. J.
Skocpol
,
M. R.
Beasley
, and
M.
Tinkham
,
J. Appl. Phys.
45
,
4054
(
1974
).
37.
H.
Courtois
,
M.
Meschke
,
J. T.
Peltonen
, and
J. P.
Pekola
,
Phys. Rev. Lett.
101
,
067002
(
2008
).
38.
B. I.
Ivanov
,
M.
Trgala
,
M.
Grajcar
,
E.
Il'ichev
, and
H.-G.
Meyer
,
Rev. Sci. Instrum.
82
,
104705
(
2011
).
39.
N.
Oukhanski
,
M.
Grajcar
,
E.
Il'ichev
, and
H.-G.
Meyer
,
Rev. Sci. Instrum.
74
,
1145
(
2003
).
40.
F.
Giazotto
,
T. T.
Heikkilä
,
A.
Luukanen
,
A. M.
Savin
, and
J. P.
Pekola
,
Rev. Mod. Phys.
78
,
217
(
2006
).
41.
A.
Barone
and
G.
Paternò
,
Physics and Applications of the Josephson Effect
(
Wiley
,
1982
).
42.
C.-H.
Yeh
,
P.-Y.
Teng
,
Y.-C.
Chiu
,
W.-T.
Hsiao
,
S. S. H.
Hsu
, and
P.-W.
Chiu
,
ACS Appl. Mater. Interfaces
11
,
6336
(
2019
).
43.
P.
Horowitz
,
W.
Hill
, and
R. J.
Rollefson
,
Am. J. Phys.
58
,
702
(
1990
).
44.
F.
Paolucci
,
N.
Ligato
,
V.
Buccheri
,
G.
Germanese
,
P.
Virtanen
, and
F.
Giazotto
, arXiv:2003.05966v2 (
2020
).
45.
K. D.
Irwin
,
Appl. Phys. Lett.
66
,
1998
(
1995
).
47.
J.
Koch
,
T. M.
Yu
,
J.
Gambetta
,
A. A.
Houck
,
D. I.
Schuster
,
J.
Majer
,
A.
Blais
,
M. H.
Devoret
,
S. M.
Girvin
, and
R. J.
Schoelkopf
,
Phys. Rev. A
76
,
042319
(
2007
).
48.
F.
Jazaeri
,
A.
Beckers
,
A.
Tajalli
, and
J.-M.
Sallese
, in
Proceedings of the 26th International Conference “Mixed Design of Integrated Circuits and System” (MIXDES)
(
2019
).