We present experimental and numerical results for a flux flow oscillator based on superconducting Josephson
junctions. Our computationally efficient theoretical model takes into account Josephson selfcoupling
of the flux
flow oscillator and is in a good agreement with our experimental results and previous studies

D Gulevich (7013726)

Feodor Kusmartsev (1251207)

P.N. Dmitriev (7172306)

V.P. Koshelets (7171592)

Nanosistemy fizika himiâ matematika

English

Loughborough University Institutional Repository

    This item was submitted to Loughborough’s Institutional Repository (https://dspace.lboro.ac.uk/) by the author and is made available under the following Creative Commons Licence conditions.       For the full text of this licence, please go to: http://creativecommons.org/licenses/by-nc-nd/2.5/  Full details of this paper can be found at: https://dspace.lboro.ac.uk/2134/15463 NANOSYSTEMS: PHYSICS, CHEMISTRY, MATHEMATICS, 2013, 4 (4), P. 507–511EXPERIMENTAL AND THEORETICALCHARACTERISATION OF JOSEPHSONSELF-COUPLING IN SUB-TERAHERTZFLUX-FLOW OSCILLATORD.R. Gulevich1, P. N. Dmitriev2, V. P. Koshelets2, F. V. Kusmartsev11Department of Physics, Loughborough University, United Kingdom2Kotel’nikov Institute of Radio Engineering and Electronics,RAS, Moscow, 125009, RussiaPACS 85.25.Cp, 74.50.+r, 07.57.HmWe present experimental and numerical results for a flux flow oscillator based on superconducting Josephsonjunctions. Our computationally efficient theoretical model takes into account Josephson self-coupling of the fluxflow oscillator and is in a good agreement with our experimental results and previous studies.Keywords: flux flow oscillator, superconducting Josephson junctions, fluxons, Terahertz radiation.1. IntroductionFlux-flow oscillator (FFO) [1] is a long Josephson junction where the flux of Josephsonvortices (fluxons) excite linear modes of electromagnetic waves. Typically, a length of a FFOis made considerably longer than the Josephson penetration length in order to accommodate achain of fluxons. Because of its potential to be used as source of Terahertz and sub-Terahertzwaves [2], it has attracted considerable attention from the scientific community. Investigationof FFO under the influence of different conditions such as noise and inhomogeneous bias hasbeen done in a range of experimental [2–4] and theoretical [5–8] studies. It has been shownthat fluxons in the long Josephson junction exhibit a rich behavior which is not limited to aone-dimensional fluxon motion, but involves a range of essentially two-dimensional effects [9],including excitation of the transverse modes [10]. In order to investigate the influence oftwo-dimensional modes on the dynamics of fluxons in a FFO, a reliable theoretical modelis needed which would take into account the Josephson self-coupling effect [11] and wouldbe efficient enough to be generalized to a model of a two-dimensional FFO which requiressubstantial computational resources. Here, we present our studies of the conventional FFObefore proceeding to the investigation of two-dimensional effects.2. Experimental resultsDesign of the harmonic mixer and all matching circuits was similar to the traditionalone that was successfully used for FFO linewidth measurements (see for details [2, 3] and [4]).Current-voltage characteristics (IVCs) of the FFO measured at different magnetic fields arepresented in Figure 1a. The level of the microwave power delivered to SIS matched to FFOis presented by color palette (blue corresponds to no power, red marks regions where inducesby FFO SIS current exceeds 25 per cent of the quasi-particle SIS jump, see Fig. 1b). Note thatthe frequency of the FFO is determined by its voltage according to the Josephson relation. Onecan see that both SIS and FFO demonstrate perfect tunnel junction behavior with small leakagecurrent and “0” return current of the FFO at large magnetic fields. The FFO provides enough508 D. R. Gulevich, P. N. Dmitriev, V. P. Koshelets, F. V. KusmartsevFIG. 1. (a) Experimental IVCs of the FFO measured at different magnetic fieldscreated by current applied to integrated control line. (b) IVC of the SIS mixer:Blue – autonomous; Red, green and cyan – pumped by the FFO at frequencies500, 600 and 700 GHz correspondingly. Note that the Josephson steps are verysharp and prominent. Positions of the steps exactly correspond to the FFOfrequency. (c) Numerical IVC for FFO according to the model described in thetext. The axes are normalized voltage vdc and bias current γ. Different linescorrespond to different values of the magnetic field h starting from 1.4 (top) to4.0 (bottom) with the step 0.05. (d) Dependence of average damping 〈α〉 versusvoltage in our numerical model. Different lines correspond to different valuesof the magnetic field h starting from 1.4 (top) to 4.0 (bottom) with the step0.05. Steps corresponding to one- two- and three-quanta transitions occurringat vg/3 = 1.9, vg/5 = 1.14 and vg/7 ≈ 0.81 are clearly visible. Four-quantatransition step at vg/9 ≈ 0.63 may also be identified for some values of magneticfield.power to pump SIS in the design frequency range 400-750 GHz (0.8-1.5 mV). Furthermore,narrowband radiation has been measured, as in the previous experimental studies [3].3. Numerical modelTo describe FFO, we use a model of a long Josephson junction [12] with x-dependentdamping parameter α(x),ϕtt + α(x)ϕt − ϕxx − β ϕxxt + sinϕ− γ = 0 (1)Characterisation of Josephson Self-coupling in Sub-Terahertz Flux-flow Oscillator 509with the boundary conditions ϕx(0, t) = ϕx(L, t) = h and homogeneous bias current γ. In ageneral case, the dependence of the bias current on x would be defined by the experimentalimplementation of the system. Our choice of neglecting the x-dependence is justified by ourinterest in universal phenomena, in a FFO rather than those related to a specific implementation.Nevertheless, the generalization to the x-dependent bias current γ = γ(x) is straightforward andwill not affect the computation time of the numerical scheme implemented here. The dampingparameter α(x) in eq. (1) is subjected to Josephson self-coupling (JSC) [11]. JSC arises as aresult of assisted quantum tunneling of quasiparticles in the presence of an AC field. We takethis effect into account by using a damping parameter which is related to the amplitude of theAC field at a each point of the junction, α(x) = γqp/vdc with [13]γqp =n=∞∑n=−∞J2n(vac2ω)γdc (vdc + 2nω) , (2)where vdc and vac are normalized DC and AC voltages: vdc = 2eVdc/h¯ωp, vac = 2eVac/h¯ωp, ωis normalized to the plasma frequency ωp, Jn are Bessel functions and γdc(v) is modelled byγdc(v) = α0 v(1 + b(v/vg)p1 + (v/vg)p).We took values of the parameters b = 35, power index p = 80, and gap voltage vg = 5.7 tobe consistent with the model used in [8]. Also, for the sake of consistency with the previousstudy [8], the length of the FFO was taken to be L = 40 Josephson penetration lengths, dampingparameters α0 = 0.033 and β = 0.035. As it is reasonable to assume that self-coupling of FFOis dominated by a single harmonic the DC and AC components can be found by approximationsvdc = 〈ϕ˙〉 and vac = 12(max ϕ˙−min ϕ˙) ,which save much of the computational time as compared to finding of the amplitudes withthe fast Fourier transform. Equation (1) was solved numerically by the explicit finite differencescheme at fixed function α(x). The self-consistent α(x) has been found by an iterative procedureuntil the desired accuracy in α(x) is reached. The number of discrete points along X was 200,the time step and the integration time were varied in accordance with the specific choice ofparameters. The code has been written in IDL programming language widely used in astronomyapplications [14]. The results of our numerical simulations are presented on Fig. 1c andare in qualitative agreement with our experimental results on Fig. 1a and earlier numericalstudies [6–8].The slopes of the IVC curves differ slightly from the experimental results at highvalues of the bias current. This can possibly be attributed to the fact that we have used ideal,non-radiative boundary conditions, while in the experimental system, a certain amount of theradiation power may escape through the boundary, causing the IVC curves to bend up as seen onthe experimental graph 1a. Fig. 1d shows our numerical results for the dependence on voltageof the average damping parameter defined as:〈α〉 = 1L∫ L0α(x)dxDistribution of the AC amplitude vac and the distribution of the damping parameter α alongthe FFO is shown on Fig. 2a,c and 2b,d for two different sets of values for magnetic fieldand bias current. The spectrum of frequencies at the FFO’s left boundary is shown on Fig. 2eand f. At certain regimes the spectrum becomes sophisticated with more harmonics turning uplike on Fig. 2f. In this case, although the dependence of the spectrum on the choice of the510 D. R. Gulevich, P. N. Dmitriev, V. P. Koshelets, F. V. KusmartsevFIG. 2. (a), (b) - distribution of the amplitude of AC drive vac on the mainharmonics ω = vdc along the FFO; (c), (d) - distribution of damping parameterin the FFO; (e), (f) - frequency spectrum at the FFO’s left boundary; X axis isnormalized frequency and Y axis shows vac. Calculations were performed fortwo set of parameters: h = 3.0, γ = 0.3 (a,c,e) and h = 2.0, γ = 0.2 (b,d,f),correspondingly.Characterisation of Josephson Self-coupling in Sub-Terahertz Flux-flow Oscillator 511calculation parameters, such as time step and maximum integration time, has not been observed,the influence of a numerical artifact may not be completely excluded.4. ConclusionTo summarize, we have presented a theoretical model which is consistent with ourexperimental results and previous theoretical studies. The model is numerically efficient, andtherefore, has the potential to be generalized for the study of the rich two-dimensional dynamicsof magnetic flux in two-dimensional FFOs.AcknowledgementsThe work was supported by the RFBR and the Ministry of Education and Science of theRussian Federation.References[1] T. Nagatsuma et al. J. Appl. Phys., 54, 3302 (1983); ibid 56, 3284 (1984); ibid 58, 441 (1985); ibid 63, 1130(1988).[2] V. P. Koshelets and S. V. Shitov. Supercond. Sci. Technol., 13, R53 (2000).[3] V. P. Koshelets et al. Supercond. Sci. Technol., 17, S127 (2004); V. P. Koshelets et al. IEEE Trans. Appl.Supercond., 17, 336 (2007); M. Yu. Torgashin et al. IEEE Trans. Appl. Supercond., 17, 339 (2007).[4] IRE RAS Lab. of Supercond. Devices web page, http://www.cplire.ru/html/lab234/publications.htm.[5] E. A. Matrozova, A. L. Pankratov and L. S. Revin. J. Appl. Phys., 112, 053905 (2012); L. S. Revin and A.L. Pankratov. Phys. Rev. B, 86, 054501 (2012); E. A. Matrozova et al. J. Appl. Phys., 110, 053922 (2011);A. L. Pankratov. Phys. Rev. B, 78, 024515 (2008); A. L. Pankratov et. al. J. Phys.: Conf. Ser., 97, 012303(2008); A. L. Pankratov, V. L. Vaks and V. P. Koshelets. J. Appl. Phys., 102, 063912 (2007).[6] A. S. Sobolev, A. L. Pankratov, and J. Mygind. Physica C, 435, 112 (2006).[7] M. Jaworski. Phys. Rev. B, 81, 224517 (2010).[8] A. L. Pankratov, A. S. Sobolev, V. P. Koshelets, and J. Mygind. Phys. Rev. B, 75, 184516 (2007).[9] K. Nakajima et al. J. Appl. Phys., 47, 1620 (1976); D. R. Gulevich and F. V. Kusmartsev. Phys. Rev. Lett., 97,017004 (2006); D. R. Gulevich and F. V. Kusmartsev. Supercond. Sci. Tech., 20, S60 (2007); D. R. Gulevichand F. V. Kusmartsev. New J. Phys., 9, 59 (2007); D.R. Gulevich et al. Physica C, 468, 1903 (2008); D. R.Gulevich et al. J. Appl. Phys., 104, 064507 (2008); H. Farhan-Hassan et al. Flux-Flow Oscillator (FFO) Madewith the Fluxon Cloning Circuits. Terahertz and Mid Infrared Radiation, p. 29, Springer (2011).[10] D. R. Gulevich et al. Phys. Rev. Lett., 101, 127002 (2008); D. R. Gulevich et al. Phys. Rev. B, 80, 094509(2009).[11] N.R. Werthamer. Phys. Rev., 147, 255 (1966); L.-E. Hasselberg, M.T. Levinsen, and M.R. Samuelsen. Phys.Rev. B, 9, 3757 (1974); V.P. Koshelets, et al. Phys. Rev. B, 56, 5572 (1997).[12] K. K. Likharev. Dynamics of Josephson Junctions and Circuits. Gordon and Breach, New York, 1986. A.Barone and G. Paterno. Physics and Applications of the Josephson Effect, Wiley, New York, 1982.[13] J. R. Tucker. IEEE Journal of Quantum Electronics, 15, 1234 (1979); J. R. Tucker and M. J. Feldman. Rev.Mod. Phys., 57, 1055 (1985).[14] Interactive Data Language (IDL): http://www.ittvis.com andhttp://idlastro.gsfc.nasa.gov.

Experimental and theoretical characterisation of Josephson self-coupling in sub-terahertz flux-flow oscillator

https://repository.lboro.ac.uk/articles/Experimental_and_theoretical_characterisation_of_Josephson_self-coupling_in_sub-terahertz_flux-flow_oscillator/9411155/files/17028971.pdf

Experimental and theoretical characterisation of Josephson self-coupling in sub-terahertz flux-flow oscillator

Abstract

Similar works

Full text

Available Versions

Loughborough University Institutional Repository