Abstract
We analyze experimentally and theoretically mutual phase locking and electromagnetic interaction between two linear arrays with a large number of Josephson junctions. Arrays with different separation, either on the same chip or on two separate substrates are studied. We observe a large coherent gain, up to a factor of three, of emitted power from two simultaneously biased arrays, compared to the sum of powers from two individually biased arrays. The phenomenon is attributed to the phase locking of junctions in different arrays via a common electromagnetic field. Remarkably, the gain can exceed the factor of two expected for a simple constructive interference of two oscillators. The larger gain is explained by an additional consequence of mutual interaction between two large arrays. Mutual phase locking of large arrays does not only result in constructive interference outside the arrays, but also improved synchronization of junctions inside each array. Our conclusion is supported by numerical modelling.
Introduction
A Josephson junction (JJ) has the unique ability to transform an applied constant voltage V into electromagnetic (EM) oscillations. The fundamental Josephson frequency, f_{J}, is connected to V via the acJosephson relation, hf_{J} = 2eV, where h is the Planck constant and e is the elementary charge. Josephson generation occurs up to the superconducting gap voltage. Therefore, f_{J} can be up to about 1 THz for lowT_{c} JJs [1] and can reach tens of terahertz for highT_{c} JJs [2,3]. Thus, a JJ has the potential to be the basis of compact, continuouswave and tunable terahertz generators, which would facilitate solving the problem of socalled “THz gap” [4].
A single JJ emits only a very small offchip power, typically in the picowatt range. To enhance it to a practical level of about 1 mW, it is necessary to combine many equivalent junctions in an array [5]. However, synchronization of a large number of radiation sources is a serious electrodynamical problem. This problem can be resolved for a group of JJs that extends in subwavelength dimensions. Such a configuration is realized for intrinsic JJs formed in a BiSCCO crystal where almost 700 JJs are localized within 1 μm [6]. The generation in a wide frequency range of 1–11 THz has been demonstrated from BiSCCO mesas containing up to 250 JJs [2]. A maximal emission of about 1 μW corresponds to inphase cavity modes in the mesas, indicating the coherent superradiant nature of the emission.
The arrays based on intrinsic JJs suffer from overheating, which impedes a raise of radiation power. This problem manifests considerably less in discrete JJ arrays. Modern lithographic technologies allow for fabricating discrete JJs with dimensions down to the sub100 nm scale [7]. However, the size of very large arrays with thousands of JJs may exceed the radiation wavelength. For such superwavelength systems, delay effects become dominant for synchronization. To reach the radiation power maximum, JJs should be synchronized with the EM mode excited within the resonator. Such large JJ arrays become similar to a laser where the junctions play the role of atoms in an active medium. The advantage of large JJ arrays working similar to lasers is discussed in more detail in [8]. The resonator can be a cavity of the JJs itself [2], an electrode with embedded JJs [9], or the dielectric substrate on which the JJ array is arranged [10].
Coherent superradiant amplification of emitted power is caused by a constructive interference of EM fields from phaselocked oscillators [11]. For two oscillators, the EM field in the farfield maximally doubles. Hence, the superradiant power, proportional to the square of the EM field, is at most four times larger than that from a single oscillator. For incoherent emission from two unlocked oscillators the power just adds up and is twice the power from a single oscillator. Therefore, the total superradiant power gain for two oscillators, defined as the ratio of coherenttoincoherent emission, is at most two. For N oscillators, the supperradiant power increase is at most N^{2} times the power from one oscillator and the superradiant gain factor is at most N. For large N, this could greatly enhance the emitted power. This is the main motivation for the development of Josephson oscillators based on arrays with many JJs [2,3,6,9,10,1215].
Resonant modes formed along five straight electrodes with niobium JJs have been directly visualized recently using lowtemperature scanning laser microscopy [13]. The scans revealed that the standing waves can provide the global coupling of all junctions in the array, that is, extended parts of the array can interact with each other. This generates the twodimensional resonant mode that should lead to the increased output power. Therefore, along with the interaction between individual JJs, there is also a mutual coupling between different arrays (including JJs and electrodes) [14]. In the ideal case of two perfectly phaselocked arrays with N JJs each, the total superradiant power is proportional to 4N^{2}, which is two times larger than the sum of powers from two incoherent arrays. Therefore, such arrays can be considered as individual oscillators, and the gain factor for two phaselocked arrays equals two. Yet, the physics of interarray coupling is much more complicated. In reality, synchronization between junctions in each array is not perfect due to the insufficient amplitude of the resonant mode. In this case, resonant coupling of two arrays may improve the state within each array. As we will show, this could increase the gain factor well above a factor of two. The goal of our work is to study interarray coupling and its manifestations.
In this work we study the interaction between two linear arrays of Nb/NbSi/Nb JJs. The arrays have a singleline geometry with 332 or 380 JJs embedded in a straight electrode. We analyze the mutual interaction between two independently biased arrays oriented parallel to each other. First, we study arrays on the same chip for different distances of 4 and 238 μm between them. Then, we consider two arrays on different chips, stacked on top of each other. We perform simultaneously the measurement of current–voltage characteristic (IVC) and bolometric analysis of the emitted radiation. In all cases, we observe clear signatures of interarray interaction. They occur when both arrays are biased at the same voltage and oscillate at the same frequency, coinciding with one of the cavity modes in the array electrodes. This leads to a profound enhancement of resonant step amplitudes in the IVCs of the arrays, indicating that the state of one array is strongly affected by oscillations in another array. The interarray coupling is manifested by a significant amplification of emitted power with a gain factor of up to three. It is well above the factor of two expected in the simple case of bare coherent superposition of oscillations. This result points out that phase locking of oscillations in the two arrays not only leads to coherent amplification of radiation. It also can improve the synchronization inside each array. The latter effect removes the limit of two for the gain factor. Finally, for better understanding, we performed numerical simulations of the inner dynamics for two interacting arrays. Our simulations confirm that two arrays can be phaselocked by a common EM field. They also provide estimation of the resulting superradiant gain. The performed experimental investigations and numerical calculations can give new ideas about the design of discrete JJ arrays that would provide more effective synchronization of JJs in order to get an output power sufficient for practical applications.
Experimental
Samples
We study samples containing one or several straight strips with embedded Nb/NbSi/Nb overlap JJs connected in series. The samples were fabricated by Oliver Kieler (Braunschweig, Germany) and were measured in AlbaNova University Center (Stockholm, Sweden). The fabrication is a selfaligning process using ebeam lithography and reactive ion etching [16,17]. Similar arrays were studied earlier in [9,12,13], where additional information about sample characterization can be found.
Figure 1a,b shows the layout of “sample1”. It has been fabricated on a 1 × 1 cm^{2} silicon substrate with the thickness 0.38 mm. It contains three closely located straight strips with a separation of only 4 μm. Each strip has the length L = 5 mm and the width w = 14 μm and contains 332 JJs distributed uniformly along the strip. The junction area is 8 × 8 μm^{2}. Contact electrodes are connected to each strip, allowing for independent biasing of each of these three arrays. Below, we will analyze the interaction between the leftmost “arraya” biased with a variable dc current and “arrayb” in the middle biased with a fixed current (Figure 1b).
Figure 1c shows the layout of “sample2”. It has a significantly larger separation of 238 μm between the adjacent linear arrays. In total, it contains 17 similar lines with 380 JJs and a total length of L = 5.7 mm. The junction area is 6 × 6 μm^{2}. Below, we will show data for the case when the rightmost “arraya” is biased with a variable dc current and the nearby “arrayb” is biased with a fixed current.
Sample1 and sample2 were used for onchip analysis where two linear arrays are placed on the same substrate. We also present data for offchip synchronization. To this end, two linear arrays were stacked on top of each other.
Radiation detection
An InSb bolometer is used for the detection of Josephson radiation. The detector and measurement procedure are the same as described in [9,12], where additional information can be found. The bolometer is based on a highpurity ndoped InSb crystal with dimension of 2–3 mm, which is placed approximately 0.5 cm above the array. The absorbed radiation causes an increase of charge carriers in the conduction band and leads to a decrease of the dc voltage at a fixed bias current. Therefore, we take the negative change of the dc voltage of the bolometer ΔU as a measure of the absorbed power, that is, ΔU = U_{0} − U > 0 where U_{0} and U are the voltages on the crystal in absence and in presence of the radiation, respectively. All measurements were performed in a liquid helium dewar at a temperatures T 4.2 K both for the samples and the detector. The calibrated responsivity of the detector at this temperature was estimated as ≃300 V/W.
Results
Figure 2a shows the individually measured IVCs of arraya and arrayb of sample1 (Figure 1a,b). The critical current in both arrays is I_{c} = 2.0–2.1 mA and the characteristic frequency, estimated within the resistively shunted junction (RSJ) model, is in the range of f_{c} ∼ 100–120 GHz. Figure 2b shows similar data for sample2. Here, for both arrays, I_{c} ≈ 2.9 mA and f_{c} ∼ 80–100 GHz. Although the area of junctions in sample2 is smaller, the critical current is slightly larger than in sample1 due to higher doping of the NbSi interlayers. Note that the abrupt transition from a superconductive to a resistive state observed in all IVCs is typical for niobium junctions with medium doping Si interlayer 11% [18].
Resonant steps at similar voltages are observed for both pairs of arrays. As shown earlier [9,13], they are caused by standing wave (cavity mode) resonances in the whole length of strips of JJs. However, for both samples the steps are more pronounced in the inner arrayb than in the outer arraya. Presumably, this is due to a more beneficial EM environment for the inner strip, which has two adjacent strips on both sides operating as additional singlestrip line resonators (see Figure 1).
In Figure 2a, we marked bias points 1–8 in arrayb, at which detailed measurements are reported below for sample1. Using the value of EM wave speed along a singlestrip line obtained in [13], we can estimate the corresponding numbers of cavity modes m = 11–18. Similarly, bias points 1–10 in Figure 2b are used for analysis of sample2. According to our estimation, they correspond to numbers of cavity modes m = 10–19.
Synchronization of two arrays on the same chip
In this work, we aim to study the EM interaction between two independently biased arrays. For the onchip measurement, the bias current in the outer arraya I_{a} was varied while the middle/inner arrayb was biased at fixed values of current I_{b} = const. All other arrays remain unbiased and, therefore, inactive. Thus, we measure the dc voltage of arrayb V_{b} and the full IVC of arraya V_{a}(I_{a}). Simultaneously, the radiation signal ΔU_{ab}(V_{a}, V_{b}) is detected. The dc bias point of arrayb (I_{b}, V_{b}) is chosen at one of the resonant steps. These points are marked by numbers in the IVCs of Figure 2.
Synchronization of closely spaced arrays (sample1)
We start with sample1, which has the shorter separation between the strips. Figure 3a shows two IVCs of arraya. The green IVC is measured with a passive arrayb, V_{b} = 0, and the red is measured with a fixed bias I_{b} = 2.98 mA and V_{b} = 123 mV, corresponding to the bias point 5 in Figure 2a. The inset shows a closeup of the voltage range V_{a} ∼ V_{b}, which demonstrates that oscillations in arrayb lead to a pronounced enhancement of the resonant step in arraya at V_{a} ≈ V_{b} while other steps are practically unaffected. The differential resistance in the center of this step R_{d} decreases by a factor of four, from R_{d} = 16 Ω to R_{d} = 4 Ω. Since the step amplitude reflects (approximately proportionally) the amplitude of the EM field in the cavity mode, this clearly demonstrates that the active arrayb amplifies the EM oscillation amplitude in arraya under the condition V_{a} ≃ V_{b}. This means almost exact equality of the mean Josephson frequencies averaged over all JJs ≃ ≈ 179 GHz, which is the necessary condition for phase locking.
In Figure 3b, we present results of the radiation detection measured simultaneously with the IVCs in Figure 3a. The lower green curve shows the detector signal ΔU_{a}(V_{a}, V_{b} = 0) as a function of the voltage V_{a} in arraya, for the unbiased arrayb, V_{b} = 0. It represents the emission power solely from arraya. The upper red curve shows similar data, ΔU_{ab}(V_{a}, V_{b}), when arrayb is biased to the point 5 in Figure 2a. It is seen that here ΔU_{ab}(V_{a} = 0) = 0.07 mV. This signal offset represents the emission ΔU_{b}(V_{b}) from arrayb alone. It can be seen that the shapes of the two curves in Figure 3b are quite similar. At almost all V_{a}, they simply differ by a constant offset, ΔU_{ab}(V_{a}, V_{b} = const) ≃ ΔU_{a}(V_{a}, V_{b} = 0) + ΔU_{b}(V_{a} = 0, V_{b} = const), as indicated by the dashed horizontal line. This implies that, usually, the powers from the two arrays simply add up, which is typical for the incoherent state. However, a remarkable peak is observed when the voltages of the two arrays practically coincide, V_{a} ≃ V_{b} = 123 mV. At this point, ΔU_{ab}(V_{a} ≃ V_{b}) = 0.34 mV, which is 2.4 times larger than the sum of individual arrays ΔU_{a}(V_{b}) + ΔU_{b}(V_{b}) ≃ 2ΔU_{b}(V_{b}) = 0.14 mV. To quantify this effect, we consider the gain factor
which describes the coherent superradiant amplification of the radiation power.
In Figure 3c,d, we show a similar analysis for all explored bias points in arrayb. Figure 3c shows the IVC of arraya I_{a}(V_{a}) without bias in arrayb (green line) and portions of the IVCs at bias points 1–8 of arrayb indicated in Figure 2a (red and blue curves). Here, we show eight fragments of separately measured IVCs close to the condition of phase locking V_{a} ∼ V_{b}. It can be seen that at all bias points, V_{b}, a strong enhancement of the resonant step in arraya occurs compared to the case without bias, V_{b} = 0. This is particularly clear for higher bias points 7 and 8, for which the steps without bias in arrayb are barely visible, but with bias they are well developed.
Figure 3d shows the detector response ΔU_{a}(V_{a}, V_{b}) measured simultaneously with the IVCs from Figure 3c. The lower green curve is measured at an unbiased arrayb. The upper red/blue curves correspond to bias points 1–8 in arrayb (Figure 2a). At low bias, bias points 1–3, we observe a multimode excitation, that is, some gain occurs even at modes adjacent to V_{b}. At higher bias, single mode amplification takes place, as for the bias point 5 discussed above (Figure 3b). The gain factors for all eight bias points are given in Table 1. The highest gain is observed at point 7 with g = 2.9.
Synchronization of more distant arrays (sample2)
Next, we consider sample2 with significantly larger separation between the arrays, namely 238 μm. Figure 4a shows IVCs of the outer arraya. The green curve shows the result without bias in the inner arrayb. The blue and red curves show results with bias at points 1–10 indicated in Figure 2b. Comparable to sample1, Figure 3c, we also can see a significant enhancement of steps in the array at V_{a} ≈ V_{b}. Figure 4b shows the detector response measured simultaneously with the IVCs from Figure 4a. A significant enhancement of emission occurs practically under the condition of phase locking, V_{a} ≃ V_{b}. The corresponding gain factors are listed in Table 2. They are only slightly smaller than those for sample1, indicating that the interarray coupling is not a shortrange phenomenon.
Offchip synchronization
The third series of measurements was performed using arrays at two different substrates. The substrates were stacked on top of each other, as sketched in Figure 5a. To facilitate access to the contact pads of the bottom array, the substrate of the top array was trimmed to a width of about 3 mm. The arrays are similar to those in sample1, but with a different shape of connecting electrodes (cf. Figure 1a,b and Figure 5a), which does not influence the measurements. The distance between two arrays is approximately equal to the thickness of substrate plus the glue layer and, possibly, a slight misalignment in the lateral direction. Overall, it is about 0.4 mm through the silicon substrate. During measurements, the bottom arraya is biased with a variable dc current and the top arrayb is biased with a fixed current. To obtain the most prominent effect, we slightly increased the temperature of the stack to T = 4.4 K by placing it above the surface of liquid helium.
Figure 5b shows the IVC of the individually biased top arrayb together with the simultaneously measured detector signal. The IVC has the characteristic parameters I_{c} = 2.1 mA and f_{c} ≈ 120 GHz. The maximum radiation signal ΔU_{b} = 0.11 mV is observed at I_{b} = 2.34 mA, V_{b} = 84 mV, which corresponds to the mean Josephson frequency = 121 GHz. At this bias point, a weak current step is observed (Figure 5b, inset), which actually is the only one in the whole IVC curve. This bias point is chosen for the subsequent measurements.
In Figure 5c, the green curve represents the individual IVC of arraya without bias in arrayb. The red curve shows the part of IVC with bias in arrayb at the maximum emission point 1 indicated in the inset in Figure 5b. This IVC has I_{c} = 1.9 mA and f_{c} ≈ 100 GHz, which are close to the values for arrayb. Figure 5d shows the simultaneously measured detector signal. The lower green curve represents the emission signal from the individual arraya, ΔU_{a}(V_{a}, V_{b} = 0). Note that it has a clearly different shape compared to ΔU_{b}(V_{a} = 0, V_{b}) shown in Figure 5b. The upper red curve represents the combined emission, ΔU_{ab}(V_{a}, V_{b}). The general behavior is quite similar to that found in the onchip experiments. A sharp peak is revealed under the condition of phase locking, V_{a} = V_{b}, with ΔU_{ab} = 0.28 mV. It is considerably larger than the sum ΔU_{a} + ΔU_{b} = 0.17 mV. The gain factor is g = 1.7. Note that, in contrast to onchip measurements, ΔU_{a} ≠ ΔU_{b}, which can depend on the difference in position of the two arrays with respect to the detector, different geometries of substrates and connecting electrodes, and the stacking arrangement of the samples.
The obtained results show that the coupling between JJ arrays can be realized due to EM waves propagating inside the substrate. Note that the distance between the arrays is close to the half wavelength in the substrate, λ_{Si}/2 = (ε_{Si} = 11.9 is the dielectric permittivity of silicon). Under this condition, the fundamental resonant mode can be excited in the substrate between the arrays. This condition is beneficial for interarray coupling.
Numerical Calculations
The experimental results presented above show that phase locking of two large JJ arrays is a complex phenomenon, which cannot be reduced to a simple constructive interference of two independent sources. For a better understanding of the phase locking dynamics, we perform numerical modelling. Figure 6a demonstrates the general view of the considered model for onchip synchronization. It contains two identical JJ arrays arranged on a common substrate with a dielectric permittivity of ε = 12, close to that of silicon. The lateral dimensions of the substrate are 2 × 0.6 mm while the thickness is 0.3 mm. We chose such a narrow substrate to avoid excitation of transverse resonant modes inside the substrate. The substrate is surrounded by vacuum, which is terminated by a perfectly matched layer (PML) to cancel back reflection. The PML conditions, in fact, simulate the walls in an anechoic chamber [19].
The circuit of each JJ array has the form of a rectangle with L = 1.8 mm in length and 0.2 mm in width (Figure 6a). The long side close to the center of substrate contains N = 100 JJs and two identical inductances at the ends. They are needed to electromagnetically decouple the lines with junctions from other parts of the circuits. A power supply of each circuit is located in the middle of the opposite long side. All lumped elements are connected by ideal conductors located at the edges of mesh. The gap between two arrays is 0.1 mm.
The JJs are described by the RSJ model [20]. The corresponding equations of junction dynamics are solved selfconsistently with Maxwell equations, which are calculated by the finitedifference timedomain (FDTD) method [21], as described in [8,12,22]. We used the following junction parameters: I_{c} = 2.5 mA, normal resistance R_{n} = 0.1 Ω, and McCumber parameter β = 2. These parameters are consistent with experimental data for Nb/NbSi/Nb junctions (Figure 2, Figure 5b,c). The inductances are equal to 100 pH while the internal resistance of the power supplies has the value of 90 Ω. The latter allows for measurements of IVCs close to the regime of constant bias current. However, the electromotive force of the power supply is, in fact, the primary source of biasing. The algorithm of the numerical calculations allows one to obtain both transport and radiation characteristics of the lines with JJs. A more detailed description of the model and the calculation procedure can be found in [8,22].
Figure 6b shows a part of the IVC of arrayb and the corresponding dependence of the radiation power P_{b}(V_{b}) for an inactive arraya. The power P is calculated by integration of the radiation pattern at the specific frequencies where the maximum in the spectrum of the ac current averaged over the JJs is observed [22]. This and subsequent simulations are performed for an upward bias sweep in a range of V that is 1.5–1.7 larger than the characteristic voltage V_{c} = I_{c}R_{n}N = 25 mV. As seen from Figure 6b, in this range the form of the IVC is close to a straight line excluding the range V ≈ 40–41 mV, where a current step is observed. The step amplitude is ΔI = 0.08 mA, and the lowest differential resistance is R_{d} = 2.8 Ω. The radiation power increases abruptly at the step and reaches the maximal value P_{b} = 0.32 μW at V_{b} = 40.06 mV, corresponding to an averaged Josephson frequency of = 193.7 GHz. The indicated value V_{b} corresponds to the electromotive force of the power supply = 398 mV. Arrayb was then constantly biased at for the subsequent analysis of interarray coupling.
The results of the simulation with two biased arrays are represented in Figure 6c. They show the calculated IVC of arraya at constantly biased arrayb as well as the calculated radiation power. It can be seen that, similar to the experimental observations (Figure 3a,c, Figure 4a, and Figure 5c), the step in the IVC becomes more pronounced compared to the previous simulations in Figure 6b. The amplitude has doubled, ΔI = 0.17 mA, and the differential resistance decreased by nearly 5.5 times. The total emitted power P_{ab}(V_{a}, V_{b}) has a nonzero offset P_{ab}(V_{a} = 0, V_{b}) ≈ 0.32 μW, corresponding to the power of the individually biased arrayb. The maximum total power of P_{ab} = 1.32 μW is observed at V_{a} = 40.19 mV (point 1 in Figure 6c). The gain factor is g ≃ 2.1. This value is consistent with the experimental values reported in Table 1 and Table 2.
Figure 6d shows the distribution of the average dc voltage on the JJs for both simulations. The green line is are for the individually biased arrayb and the red/blue line is for the collectively biased arraya and arrayb. Here, we can clearly see a signature of standing waves along the arrays. We can see eight flat regions with almost equal junction voltages and frequencies = 193.3 GHz. These junctions are in the antinodes of the cavity mode and are synchronized by the EM field of the standing wave oscillating at the frequency .
The junctions located at the nodes of the resonant mode are asynchronous. This means that their Josephson frequencies differ from the radiation frequency , actually > . A similar pattern was obtained in [8], but that inequality was opposite, that is, < . As follows from [8], the latter relation occurs at β ≪ 1. The small value of β also establishes the inverted shape of current steps in the IVC compared to that shown in Figure 6b,c. It can be shown in the same manner as in [8] that the relation between f_{J} of synchronous and asynchronous junctions as well as the shape of current steps changes to the opposite at β > 1. Also note that the difference of asynchronous regions for arraya and arrayb, which is clearly seen in Figure 6d, is caused by the different biasing sequences and the corresponding historydependent dynamics.
Discussion
Our experimental data and numerical simulations demonstrate that large JJ arrays can be effectively coupled to each other, resulting in a coherent superradiant enhancement of the emission power. The amplification is observed in the frequency range of 100–200 GHz, both for arrays on a common substrate and for arrays on different substrates formed in a stack. We explain this effect by the interaction between JJ arrays via an EM field. This field is excited along the surface of the substrate as well as inside the substrate. The amplification tends to grow with an increase of the frequency although the overall radiation power decreases. Similar results were obtained for Bi_{2}Sr_{2}CaCu_{2}O_{8+}_{x} mesa structures in [14,15]. In [14], three simultaneously biased mesas emit a high power of 610 μW while each mesa alone emits a maximum of 120 μW. Following our terminology, this corresponds to a gain factor of g = 1.7. In [15], the interaction between two mesas has been revealed via the study of polarization of EM emission. Similar to the present work, the obtained data allows one to conclude that such intrinsic JJ arrays have a mutual coupling through the common substrate.
As described in the Introduction, the simple constructive interference of two oscillators, be it single JJs [23] or arrays, should result in a gain of g 2. However, in this work, we observe also significantly larger gains. Although this clearly indicates that coherent emission from both arrays takes place, it also indicates that additional more complex phenomena are involved. The clue to understanding is provided by the inset in Figure 3a, which demonstrates that phase locking of the two arrays leads to enhancement of the oscillation amplitude in arraya. From Figure 2, we observe that resonant steps in individually biased arrays are more pronounced in the inner arrayb than in the outer arraya. At higher bias, steps in arraya are almost invisible. Let us suppose that, initially, only arrayb is synchronized at the cavity mode and emits radiation while arraya is not synchronized and, therefore, practically not emitting. In this case, if interarray coupling totally synchronizes arraya, then the gain factor would become four. This explanation is consistent with the observation that g > 2 is observed for higher steps (Table 1 and Table 2), which are less pronounced in the individually biased arraya. This is also confirmed by numerical simulations where we also observed g > 2. Therefore, the gain is caused both by the coherent superradiant effect and by the enhancement of the oscillating EM field in each array. The latter is an additional factor that explains why/how the gain factor could be larger than two. The key is that, in our case, the interaction takes place between large arrays with many oscillating junctions. This is a much more complex phenomenon than locking of two oscillators. Here, a mutual synchronization of the two arrays assists also in better internal synchronization within each of the arrays.
The simulated voltage profile in Figure 6d clearly indicates that the cavity mode is playing a decisive role for synchronization of the array. Junctions in the antinodal regions are phaselocked by the driving EM field of the cavity mode. In the nodal regions, the driving force is very small, and, therefore, JJs are unsynchronized there. These asynchronous nodal regions make a vanishingly small contribution to the radiation power (see Appendix). Therefore, the overall emission spectrum remains very sharp and is practically not influenced by the voltage/frequency deviation at nodal JJs. In fact, it is the cavity mode in the electrode, rather than individual junctions, that is causing the emission. The role of the JJs is just to excite and pump energy into the mode.
We have observed similar coherent gains for the onchip case with small, 4 μm (sample1), and significantly larger, 238 μm (sample2), separation between the arrays as well as for offchip measurements on different chips with even larger separation ∼400 μm. This demonstrates that arrays can effectively interact at a fairly long range. It would be interesting to study in more detail how the amplification depends on the separation between arrays. We have already access to a suitable Nb array consisting of several subarrays with different distances between each. Hence, the corresponding measurements may be carried out in close future.
In the simulations, we see the same effect of amplification for the total power radiated by JJ arrays at frequencies near 190 GHz. The amplification is comparable to that from the measurements. We present in the Appendix the distribution of power generated by each junction and of the phase shift between ac voltage and ac current. This analysis gives a visual pattern of the considered effect of amplification as well as of the synchronization of JJs.
Conclusion
We explored experimentally and numerically the EM interaction between large JJ arrays containing more than 300 JJs and having superwavelength dimensions. The studied Nb/NbSi/Nb JJ arrays exhibit strong cavity mode resonances, caused by the formation of standing waves along the whole length of the array [13]. We observed mutual coupling between the arrays both on the same chip and on different chips and at different separations between the arrays. We reported significant coherent amplification of radiation emission when both arrays are brought to the same cavity mode resonance. A coherent gain factor, that is, the ratio of the joint emission power from the two arrays to the sum of powers from individually biased arrays, as large as 2.9 was observed. This is well beyond the limit of two, characteristic for the bare constructive interference of two oscillators. The large gain factor indicates that additional effects are taking place. The key is that, in our case, the interaction takes place between large arrays with many oscillating junctions. This is a much more complex phenomenon than locking of two oscillators. Here, a mutual synchronization of the two arrays assists also in better internal synchronization within each of the arrays. This facilitates gains larger than two. This conclusion has been supported by the performed numerical simulations.
Finally, we note that the effect of coherent amplification of radiation from several coupled JJ arrays depends on a number of geometrical factors and material parameters. Those should be accounted for in the design and fabrication of large JJ arrays. Proper design, in which this effect is maximally manifested, allows for increasing the output radiation power, which will facilitate the implementation of JJ arrays in practical applications.
Appendix: Additional Information about Numerical Simulations
The numerical algorithm allows for calculating also the energy parameters for all JJs that can facilitate the study of synchronization of junctions in the arrays (Figure 6a). We can define the work of the nth JJ under an EM field per unit time, or the generated power, as P_{n} = −Re()/2 = , where φ_{n} is the phase shift between and , that is, φ_{n} = arg() − arg(). As well as for the radiation power P, these amplitudes are taken at the frequency corresponding to the maximum of
As seen from Figure 7a, if arraya is inactive and arrayb is biased to the point 2 in Figure 6b, out of the current step, then P_{n} is distributed along arrayb rather chaotically taking values of both signs. Almost half of the junctions have a negative sign of P_{n}. This means that the field does a positive work under these junctions. Hence, these junctions operate as consumers, not as generators. The total generated power in the array is only P = = 0.3 nW. As seen from Figure 6b the radiation power in point 2 is also practically zero on the scale of microwatts. The phase shift φ_{n} also has a chaotic character taking values in a wide range (Figure 7b). For 77 junctions, φ_{n} ranges from 180° to 360°. Thus, the differential impedance Z of most of the junctions acquires a capacitive character.
When arrayb is biased within the current step, and arraya is inactive (point 1 in Figure 6b), the distributions P_{n} and φ_{n} along arrayb become structured. As seen from Figure 7c,d, they have an oscillation form similar to the distribution of the dc voltage in JJs (Figure 6d). Therefore, P_{n} and φ_{n} exhibit, here, the resonant mode excited in the array. As seen from Figure 7c, the junctions in the nodes do a significantly smaller work under the field compared to the junctions in the antinodes, which generate a power of P_{n} = 4–6 nW. Regarding the phase shift, the junctions in the antinodes have φ_{n} 180° (Figure 7d), which is characteristic for a generator with small capacitance. However, when moving towards the nodes, φ_{n} decreases down to 110–120°, that is, the differential impedance Z of the JJs acquires an inductive character.
Similar patterns of P_{n} and φ_{n} are observed for both arrays when arrayb is biased within the current step (point 1 in Figure 6c). However, the junctions in the antinodes of the resonant mode now generate a larger power: P_{n} = 7–9 nW for arrayb and P_{n} = 9–11 nW for arraya (Figure 7c). Moreover, the range of junctions that do a large work under the field slightly widens in each antinode. This is in agreement with the conclusion made from the radiative analysis about the amplification of Josephson radiation when two arrays are biased (Figure 6b,c). Each array does a larger work under the field when two arrays are simultaneously biased compared to the case in which only one array is active. Calculating the total generated power P for all three cases in Figure 7c, we obtain 0.4 μW from arrayb with inactive arraya and 1.35 μW from the arrays when they are both biased. These values are slightly less than those of the radiation power P_{b} and P_{ab} (Figure 6b,c) calculated by the integration of radiation pattern. This decrease can be caused by dissipations of the radiation in the active loads of the power supplies (Figure 6a).
The phase shift φ_{n} in the arrays undergoes slight changes when two arrays are biased instead of one (Figure 7d). In each antinode of the mode of arrayb, φ_{n} has two local maxima instead of one with the widening of range where φ_{n} 180°. For arraya, the phase shift in the antinodes decreases to φ_{n} ≈ 160–170°, that is, Z of the effectively generating JJs acquires an inductive character.
Keep in mind that I_{n} is determined by the resonant mode. Thus, it has an equal oscillation frequency throughout each array [8]. Moreover, when the active arrays are coupled, the common resonant mode is formed. Therefore, the data presented in Figure 6d and Figure 7c,d allow one to conclude that the JJs in antinodes of the resonant mode become phaselocked via the common EM field. Such phase locking is fully constructive, that is, all junctions in the antinodes do a positive work under the field. This is what essentially provides the amplification of the Josephson radiation. The present conclusion is in accordance with the data of direct visualisation of the modes presented in [13], in which the global synchronization of the whole JJ array was indicated.
References

Villegier, J.C. Refractory Niobium Nitride NbN Josephson Junctions and Applications. In Josephson Junctions: History, Devices, and Applications; Wolf, E. L.; Arnold, G. B.; Gurvitch, M. A.; Zasadzinski, J. F., Eds.; Pan Stanford Publishing: Singapore, 2017; pp 147–183.
Return to citation in text: [1] 
Borodianskyi, E. A.; Krasnov, V. M. Nat. Commun. 2017, 8, 1742. doi:10.1038/s41467017018884
Return to citation in text: [1] [2] [3] [4] 
Kleiner, R.; Wang, H. J. Appl. Phys. 2019, 126, 171101. doi:10.1063/1.5116660
Return to citation in text: [1] [2] 
Tonouchi, M. Nat. Photonics 2007, 1, 97–105. doi:10.1038/nphoton.2007.3
Return to citation in text: [1] 
Darula, M.; Doderer, T.; Beuven, S. Supercond. Sci. Technol. 1999, 12, R1–R25. doi:10.1088/09532048/12/1/001
Return to citation in text: [1] 
Kakeya, I.; Wang, H. Supercond. Sci. Technol. 2016, 29, 073001. doi:10.1088/09532048/29/7/073001
Return to citation in text: [1] [2] 
Morosh, V.; Linek, J.; Müller, B.; MartínezPérez, M. J.; Wolter, S.; Weimann, T.; Beyer, J.; Schurig, T.; Kieler, O.; Zorin, A. B.; Kleiner, R.; Koelle, D. Phys. Rev. Appl. 2020, 14, 054072. doi:10.1103/physrevapplied.14.054072
Return to citation in text: [1] 
Galin, M. A.; Shereshevsky, I. A.; Vdovicheva, N. K.; Kurin, V. V. Supercond. Sci. Technol. 2021, 34, 075005. doi:10.1088/13616668/abfd0b
Return to citation in text: [1] [2] [3] [4] [5] [6] [7] 
Galin, M. A.; Klushin, A. M.; Kurin, V. V.; Seliverstov, S. V.; Finkel, M. I.; Goltsman, G. N.; Müller, F.; Scheller, T.; Semenov, A. D. Supercond. Sci. Technol. 2015, 28, 055002. doi:10.1088/09532048/28/5/055002
Return to citation in text: [1] [2] [3] [4] [5] 
Song, F.; Müller, F.; Scheller, T.; Semenov, A.; He, M.; Fang, L.; Hübers, H.W.; Klushin, A. M. Appl. Phys. Lett. 2011, 98, 142506. doi:10.1063/1.3576910
Return to citation in text: [1] [2] 
Rabinovich, M. I.; Trubetskov, D. I. SelfExcited Oscillations in Multifrequency Systems. In Oscillations and Waves in Linear and Nonlinear Systems; Hazewinkel, M., Ed.; Kluwer Academic Publishers: Dordrecht, Netherlands, 1989; pp 329–352. doi:10.1007/9789400910331_16
Return to citation in text: [1] 
Galin, M. A.; Borodianskyi, E. A.; Kurin, V. V.; Shereshevskiy, I. A.; Vdovicheva, N. K.; Krasnov, V. M.; Klushin, A. M. Phys. Rev. Appl. 2018, 9, 054032. doi:10.1103/physrevapplied.9.054032
Return to citation in text: [1] [2] [3] [4] 
Galin, M. A.; Rudau, F.; Borodianskyi, E. A.; Kurin, V. V.; Koelle, D.; Kleiner, R.; Krasnov, V. M.; Klushin, A. M. Phys. Rev. Appl. 2020, 14, 024051. doi:10.1103/physrevapplied.14.024051
Return to citation in text: [1] [2] [3] [4] [5] [6] [7] 
Benseman, T. M.; Gray, K. E.; Koshelev, A. E.; Kwok, W.K.; Welp, U.; Minami, H.; Kadowaki, K.; Yamamoto, T. Appl. Phys. Lett. 2013, 103, 022602. doi:10.1063/1.4813536
Return to citation in text: [1] [2] [3] [4] 
Kobayashi, R.; Hayama, K.; Fujita, S.; Tsujimoto, M.; Kakeya, I. Phys. Rev. Appl. 2022, 17, 054043. doi:10.1103/physrevapplied.17.054043
Return to citation in text: [1] [2] [3] 
Mueller, F.; Behr, R.; Weimann, T.; Palafox, L.; Olaya, D.; Dresselhaus, P. D.; Benz, S. P. IEEE Trans. Appl. Supercond. 2009, 19, 981–986. doi:10.1109/tasc.2009.2017911
Return to citation in text: [1] 
Kieler, O.; Wendisch, R.; Gerdau, R.W.; Weimann, T.; Kohlmann, J.; Behr, R. IEEE Trans. Appl. Supercond. 2021, 31, 1100705. doi:10.1109/tasc.2021.3060678
Return to citation in text: [1] 
Gudkov, A. L.; Kupriyanov, M. Y.; Samus’, A. N. J. Exp. Theor. Phys. 2012, 114, 818–829. doi:10.1134/s1063776112030144
Return to citation in text: [1] 
Taflove, A.; Hagness, S. C. Perfectly Matched Layer Absorbing Boundary Conditions. Computational Electrodynamics: The FiniteDifference TimeDomain Method; Artech House Inc.: Boston/London, 2005; pp 273–328.
Return to citation in text: [1] 
Likharev, K. K. The Josephson effect. Dynamics of Josephson Junctions and Circuits; Gordon and Breach Science Publishers: Paris, France, 1986; pp 1–26.
Return to citation in text: [1] [2] 
Yee, K. IEEE Trans. Antennas Propag. 1966, 14, 302–307. doi:10.1109/tap.1966.1138693
Return to citation in text: [1] 
Kurin, V. V.; Vdovicheva, N. K.; Shereshevskii, I. A. Radiophys. Quantum Electron. (Engl. Transl.) 2017, 59, 922–936. doi:10.1007/s1114101797626
Return to citation in text: [1] [2] [3] 
Finnegan, T. F.; Wahlsten, S. Appl. Phys. Lett. 1972, 21, 541–544. doi:10.1063/1.1654250
Return to citation in text: [1]
22.  Kurin, V. V.; Vdovicheva, N. K.; Shereshevskii, I. A. Radiophys. Quantum Electron. (Engl. Transl.) 2017, 59, 922–936. doi:10.1007/s1114101797626 
8.  Galin, M. A.; Shereshevsky, I. A.; Vdovicheva, N. K.; Kurin, V. V. Supercond. Sci. Technol. 2021, 34, 075005. doi:10.1088/13616668/abfd0b 
8.  Galin, M. A.; Shereshevsky, I. A.; Vdovicheva, N. K.; Kurin, V. V. Supercond. Sci. Technol. 2021, 34, 075005. doi:10.1088/13616668/abfd0b 
1.  Villegier, J.C. Refractory Niobium Nitride NbN Josephson Junctions and Applications. In Josephson Junctions: History, Devices, and Applications; Wolf, E. L.; Arnold, G. B.; Gurvitch, M. A.; Zasadzinski, J. F., Eds.; Pan Stanford Publishing: Singapore, 2017; pp 147–183. 
6.  Kakeya, I.; Wang, H. Supercond. Sci. Technol. 2016, 29, 073001. doi:10.1088/09532048/29/7/073001 
14.  Benseman, T. M.; Gray, K. E.; Koshelev, A. E.; Kwok, W.K.; Welp, U.; Minami, H.; Kadowaki, K.; Yamamoto, T. Appl. Phys. Lett. 2013, 103, 022602. doi:10.1063/1.4813536 
8.  Galin, M. A.; Shereshevsky, I. A.; Vdovicheva, N. K.; Kurin, V. V. Supercond. Sci. Technol. 2021, 34, 075005. doi:10.1088/13616668/abfd0b 
5.  Darula, M.; Doderer, T.; Beuven, S. Supercond. Sci. Technol. 1999, 12, R1–R25. doi:10.1088/09532048/12/1/001 
16.  Mueller, F.; Behr, R.; Weimann, T.; Palafox, L.; Olaya, D.; Dresselhaus, P. D.; Benz, S. P. IEEE Trans. Appl. Supercond. 2009, 19, 981–986. doi:10.1109/tasc.2009.2017911 
17.  Kieler, O.; Wendisch, R.; Gerdau, R.W.; Weimann, T.; Kohlmann, J.; Behr, R. IEEE Trans. Appl. Supercond. 2021, 31, 1100705. doi:10.1109/tasc.2021.3060678 
13.  Galin, M. A.; Rudau, F.; Borodianskyi, E. A.; Kurin, V. V.; Koelle, D.; Kleiner, R.; Krasnov, V. M.; Klushin, A. M. Phys. Rev. Appl. 2020, 14, 024051. doi:10.1103/physrevapplied.14.024051 
2.  Borodianskyi, E. A.; Krasnov, V. M. Nat. Commun. 2017, 8, 1742. doi:10.1038/s41467017018884 
3.  Kleiner, R.; Wang, H. J. Appl. Phys. 2019, 126, 171101. doi:10.1063/1.5116660 
6.  Kakeya, I.; Wang, H. Supercond. Sci. Technol. 2016, 29, 073001. doi:10.1088/09532048/29/7/073001 
9.  Galin, M. A.; Klushin, A. M.; Kurin, V. V.; Seliverstov, S. V.; Finkel, M. I.; Goltsman, G. N.; Müller, F.; Scheller, T.; Semenov, A. D. Supercond. Sci. Technol. 2015, 28, 055002. doi:10.1088/09532048/28/5/055002 
10.  Song, F.; Müller, F.; Scheller, T.; Semenov, A.; He, M.; Fang, L.; Hübers, H.W.; Klushin, A. M. Appl. Phys. Lett. 2011, 98, 142506. doi:10.1063/1.3576910 
12.  Galin, M. A.; Borodianskyi, E. A.; Kurin, V. V.; Shereshevskiy, I. A.; Vdovicheva, N. K.; Krasnov, V. M.; Klushin, A. M. Phys. Rev. Appl. 2018, 9, 054032. doi:10.1103/physrevapplied.9.054032 
13.  Galin, M. A.; Rudau, F.; Borodianskyi, E. A.; Kurin, V. V.; Koelle, D.; Kleiner, R.; Krasnov, V. M.; Klushin, A. M. Phys. Rev. Appl. 2020, 14, 024051. doi:10.1103/physrevapplied.14.024051 
14.  Benseman, T. M.; Gray, K. E.; Koshelev, A. E.; Kwok, W.K.; Welp, U.; Minami, H.; Kadowaki, K.; Yamamoto, T. Appl. Phys. Lett. 2013, 103, 022602. doi:10.1063/1.4813536 
15.  Kobayashi, R.; Hayama, K.; Fujita, S.; Tsujimoto, M.; Kakeya, I. Phys. Rev. Appl. 2022, 17, 054043. doi:10.1103/physrevapplied.17.054043 
23.  Finnegan, T. F.; Wahlsten, S. Appl. Phys. Lett. 1972, 21, 541–544. doi:10.1063/1.1654250 
2.  Borodianskyi, E. A.; Krasnov, V. M. Nat. Commun. 2017, 8, 1742. doi:10.1038/s41467017018884 
3.  Kleiner, R.; Wang, H. J. Appl. Phys. 2019, 126, 171101. doi:10.1063/1.5116660 
13.  Galin, M. A.; Rudau, F.; Borodianskyi, E. A.; Kurin, V. V.; Koelle, D.; Kleiner, R.; Krasnov, V. M.; Klushin, A. M. Phys. Rev. Appl. 2020, 14, 024051. doi:10.1103/physrevapplied.14.024051 
13.  Galin, M. A.; Rudau, F.; Borodianskyi, E. A.; Kurin, V. V.; Koelle, D.; Kleiner, R.; Krasnov, V. M.; Klushin, A. M. Phys. Rev. Appl. 2020, 14, 024051. doi:10.1103/physrevapplied.14.024051 
2.  Borodianskyi, E. A.; Krasnov, V. M. Nat. Commun. 2017, 8, 1742. doi:10.1038/s41467017018884 
10.  Song, F.; Müller, F.; Scheller, T.; Semenov, A.; He, M.; Fang, L.; Hübers, H.W.; Klushin, A. M. Appl. Phys. Lett. 2011, 98, 142506. doi:10.1063/1.3576910 
14.  Benseman, T. M.; Gray, K. E.; Koshelev, A. E.; Kwok, W.K.; Welp, U.; Minami, H.; Kadowaki, K.; Yamamoto, T. Appl. Phys. Lett. 2013, 103, 022602. doi:10.1063/1.4813536 
8.  Galin, M. A.; Shereshevsky, I. A.; Vdovicheva, N. K.; Kurin, V. V. Supercond. Sci. Technol. 2021, 34, 075005. doi:10.1088/13616668/abfd0b 
11.  Rabinovich, M. I.; Trubetskov, D. I. SelfExcited Oscillations in Multifrequency Systems. In Oscillations and Waves in Linear and Nonlinear Systems; Hazewinkel, M., Ed.; Kluwer Academic Publishers: Dordrecht, Netherlands, 1989; pp 329–352. doi:10.1007/9789400910331_16 
15.  Kobayashi, R.; Hayama, K.; Fujita, S.; Tsujimoto, M.; Kakeya, I. Phys. Rev. Appl. 2022, 17, 054043. doi:10.1103/physrevapplied.17.054043 
7.  Morosh, V.; Linek, J.; Müller, B.; MartínezPérez, M. J.; Wolter, S.; Weimann, T.; Beyer, J.; Schurig, T.; Kieler, O.; Zorin, A. B.; Kleiner, R.; Koelle, D. Phys. Rev. Appl. 2020, 14, 054072. doi:10.1103/physrevapplied.14.054072 
8.  Galin, M. A.; Shereshevsky, I. A.; Vdovicheva, N. K.; Kurin, V. V. Supercond. Sci. Technol. 2021, 34, 075005. doi:10.1088/13616668/abfd0b 
2.  Borodianskyi, E. A.; Krasnov, V. M. Nat. Commun. 2017, 8, 1742. doi:10.1038/s41467017018884 
9.  Galin, M. A.; Klushin, A. M.; Kurin, V. V.; Seliverstov, S. V.; Finkel, M. I.; Goltsman, G. N.; Müller, F.; Scheller, T.; Semenov, A. D. Supercond. Sci. Technol. 2015, 28, 055002. doi:10.1088/09532048/28/5/055002 
14.  Benseman, T. M.; Gray, K. E.; Koshelev, A. E.; Kwok, W.K.; Welp, U.; Minami, H.; Kadowaki, K.; Yamamoto, T. Appl. Phys. Lett. 2013, 103, 022602. doi:10.1063/1.4813536 
15.  Kobayashi, R.; Hayama, K.; Fujita, S.; Tsujimoto, M.; Kakeya, I. Phys. Rev. Appl. 2022, 17, 054043. doi:10.1103/physrevapplied.17.054043 
18.  Gudkov, A. L.; Kupriyanov, M. Y.; Samus’, A. N. J. Exp. Theor. Phys. 2012, 114, 818–829. doi:10.1134/s1063776112030144 
9.  Galin, M. A.; Klushin, A. M.; Kurin, V. V.; Seliverstov, S. V.; Finkel, M. I.; Goltsman, G. N.; Müller, F.; Scheller, T.; Semenov, A. D. Supercond. Sci. Technol. 2015, 28, 055002. doi:10.1088/09532048/28/5/055002 
12.  Galin, M. A.; Borodianskyi, E. A.; Kurin, V. V.; Shereshevskiy, I. A.; Vdovicheva, N. K.; Krasnov, V. M.; Klushin, A. M. Phys. Rev. Appl. 2018, 9, 054032. doi:10.1103/physrevapplied.9.054032 
13.  Galin, M. A.; Rudau, F.; Borodianskyi, E. A.; Kurin, V. V.; Koelle, D.; Kleiner, R.; Krasnov, V. M.; Klushin, A. M. Phys. Rev. Appl. 2020, 14, 024051. doi:10.1103/physrevapplied.14.024051 
9.  Galin, M. A.; Klushin, A. M.; Kurin, V. V.; Seliverstov, S. V.; Finkel, M. I.; Goltsman, G. N.; Müller, F.; Scheller, T.; Semenov, A. D. Supercond. Sci. Technol. 2015, 28, 055002. doi:10.1088/09532048/28/5/055002 
12.  Galin, M. A.; Borodianskyi, E. A.; Kurin, V. V.; Shereshevskiy, I. A.; Vdovicheva, N. K.; Krasnov, V. M.; Klushin, A. M. Phys. Rev. Appl. 2018, 9, 054032. doi:10.1103/physrevapplied.9.054032 
8.  Galin, M. A.; Shereshevsky, I. A.; Vdovicheva, N. K.; Kurin, V. V. Supercond. Sci. Technol. 2021, 34, 075005. doi:10.1088/13616668/abfd0b 
12.  Galin, M. A.; Borodianskyi, E. A.; Kurin, V. V.; Shereshevskiy, I. A.; Vdovicheva, N. K.; Krasnov, V. M.; Klushin, A. M. Phys. Rev. Appl. 2018, 9, 054032. doi:10.1103/physrevapplied.9.054032 
22.  Kurin, V. V.; Vdovicheva, N. K.; Shereshevskii, I. A. Radiophys. Quantum Electron. (Engl. Transl.) 2017, 59, 922–936. doi:10.1007/s1114101797626 
8.  Galin, M. A.; Shereshevsky, I. A.; Vdovicheva, N. K.; Kurin, V. V. Supercond. Sci. Technol. 2021, 34, 075005. doi:10.1088/13616668/abfd0b 
22.  Kurin, V. V.; Vdovicheva, N. K.; Shereshevskii, I. A. Radiophys. Quantum Electron. (Engl. Transl.) 2017, 59, 922–936. doi:10.1007/s1114101797626 
20.  Likharev, K. K. The Josephson effect. Dynamics of Josephson Junctions and Circuits; Gordon and Breach Science Publishers: Paris, France, 1986; pp 1–26. 
21.  Yee, K. IEEE Trans. Antennas Propag. 1966, 14, 302–307. doi:10.1109/tap.1966.1138693 
19.  Taflove, A.; Hagness, S. C. Perfectly Matched Layer Absorbing Boundary Conditions. Computational Electrodynamics: The FiniteDifference TimeDomain Method; Artech House Inc.: Boston/London, 2005; pp 273–328. 
20.  Likharev, K. K. The Josephson effect. Dynamics of Josephson Junctions and Circuits; Gordon and Breach Science Publishers: Paris, France, 1986; pp 1–26. 
9.  Galin, M. A.; Klushin, A. M.; Kurin, V. V.; Seliverstov, S. V.; Finkel, M. I.; Goltsman, G. N.; Müller, F.; Scheller, T.; Semenov, A. D. Supercond. Sci. Technol. 2015, 28, 055002. doi:10.1088/09532048/28/5/055002 
13.  Galin, M. A.; Rudau, F.; Borodianskyi, E. A.; Kurin, V. V.; Koelle, D.; Kleiner, R.; Krasnov, V. M.; Klushin, A. M. Phys. Rev. Appl. 2020, 14, 024051. doi:10.1103/physrevapplied.14.024051 
13.  Galin, M. A.; Rudau, F.; Borodianskyi, E. A.; Kurin, V. V.; Koelle, D.; Kleiner, R.; Krasnov, V. M.; Klushin, A. M. Phys. Rev. Appl. 2020, 14, 024051. doi:10.1103/physrevapplied.14.024051 
© 2022 Galin et al.; licensee BeilsteinInstitut.
This is an open access article licensed under the terms of the BeilsteinInstitut Open Access License Agreement (https://www.beilsteinjournals.org/bjnano/terms), which is identical to the Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0). The reuse of material under this license requires that the author(s), source and license are credited. Thirdparty material in this article could be subject to other licenses (typically indicated in the credit line), and in this case, users are required to obtain permission from the license holder to reuse the material.