Hyperchaotic behaviours and controlling hyperchaos in an array of RCL-shunted Josephson junctions Ri Ilmyong( ) a)b), Feng Yu-Ling( ) a), Yao Zhi-Hai( ) a), and Fan Jian( ) a) a) Department of Physics, Science College, Changchun University of Science and Technology, Changchun, China b) Department of Physics, University of Science, Unjong Strict, Pyongyang, DPR Korea (Received June ; revised manuscript received July ) This paper deals with dynamical behaviours in an array composed of two resistive-capacitive-inductive-shunted (RCL-shunted) Josephson junctions (RCLSJJs) and a shunted resistor. Numerical simulations show that periodic, chaotic and hyperchaotic states can coexist in this array. Moreover, a scheme for controlling hyperchaos in this array is presented by adjusting the external bias current. Numerical results confirm that this scheme can be effectively used to control hyperchaotic states in this array into stable periodic states, and different stable periodic states with different period numbers can be obtained by appropriately choosing the intensity of the external bias current. Keywords: hyperchaos, controlling hyperchaos, array of RCL-shunted Josephson junctions, Lyapunov exponent PACS:..Pq,..+r DOI:.88/-///. Introduction With the rapid development of fabrication technology and superconducting materials, Josephson junctions have been extensively utilized, and high critical-current junctions are preferred in many applications. [] In order to further study Josephson junctions with higher critical currents, the resistive-capacitive-inductive-shunted Josephson junction (RCLSJJ) model was proposed [] based on previous investigations of the resistive-capacitive shunted Josephson junction (RCSJJ) model. [] Moreover, the RCLSJJ model has avoided the shortcomings of the RCSJJ model and shown the agreement between the experimental and the numerical simulated I V curves when the shunt of the junction contains an inductive component. [,] Following this, there have been some reports on chaos, chaos control and chaotic synchronization in RCLSJJs. [ ] Recently, much attention has been paid to the dynamical characteristics of the Josephson junction array (JJA), since the output power of a single junction is extremely low (typically nw) and higher output power values can be obtained by using JJAs. Chaos and nonlinearity in JJAs have been interesting subjects of research [ ] since JJAs could be used as high precision voltage standards, high power coherent THz sources, a new-type inductor, etc. [,] Basler et al. [] reported on the theory of phase locking in a small Josephson junction cell by using the resistively shunted junction (RSJ) model. Hassel et al. [8] studied self-synchronization in distributed Josephson junction arrays by using the RSJ model. Grib et al. [9] reported on the synchronization of overdamped RCSJJs shunted by a superconducting resonator. Bhagavatula et al. [] numerically investigated spatiotemporal chaos in the JJA composed of RSJs. Chernikov and Schmidt [] studied JJAs composed of RSJs. The dynamical equation set was approximated in their research. They found that when four or more units are globally coupled by a common resistance, the system exhibits adiabatic chaos. Zhou et al. [] studied a JJA composed of RCSJJs. The dynamical equations were not approximated. Direct calculation results showed that phase locking and chaos coexist when there are three junctions in the JJA. Zhang and Shen, [] and Feng et al. [] completed controlling hyperchaos in erbiumdoped fiber laser and degenerate optical parametric oscillators. In the present paper, we investigate chaos, hyperchaos and controlling hyperchaos in a JJA composed of RCLSJJs. The numerical results show that chaos and hyperchaos states can coexist in this array when it contains two RCLSJJs, no matter whether the original states of two RCLSJJs are chaotic. Moreover, we achieve the control of hyperchaos in this array by adjusting the intensity of the external bias current. Corresponding author. E-mail: fyl89@yahoo.com.cn c Chinese Physical Society and IOP Publishing Ltd http://www.iop.org/journals/cpb http://cpb.iphy.ac.cn -
. Array model and equation In this section, we present an array model which consists of two RCLSJJs, a shunted resistance R and a direct current source. The circuit schematic diagram is shown in Fig.. and all currents normalized to I c, with taking I c = I c. The normalized junction voltage: υ n = V n /I c R s. The normalized external dc bias current: i = I dc /I c. The normalized shunt current: i sn = I sn /I c. Other dimensionless parameters are as follows: β Cn = ei cn R snc n /, β Ln = ei cn L n /, g n = R sn /R V, and = = = R V /R = R V /R is used to denote the coupling intensity between two Josephson junctions (JJs) used in the array. Fig.. Schematic diagram of the Josephson junction array. If all the noise terms are neglected, its complete dynamical equations are as follows: C n dv n dt dγ n = V n, () e dt + V n R Vn + I cn (sin γ n ) + I sn + I R = I dc, () L n di sn dt I R = R + I sn R sn = V n, () V n, () n= where subscripts n =, denote the serial numbers of two RCLSJJs used in the array; C n, I cn, and γ n are the junction capacitance, the junction critical current and the phase difference of the superconducting order parameters across the n-th RCLSJJ, respectively; V n is the junction voltage; R V n is the nonlinear resistance; R sn and L n are, respectively, the shunt resistance and shunt inductance; I sn is the shunt current flowing through R sn ; I dc is the external bias dc current; is the Planck constant; e is the electron charge. For the convenience of numerical analysis, Eqs. () () are expressed in dimensionless form as follows: dγ n d = υ n, () dυ n d = [i g n υ n sin(γ n ) i sn β Cn g n n (υ + υ )], () di sn d = (υ n i sn ). β Ln () Here, the normalized time = ω c t with ω c = (ei c R s )/. All voltages are normalized to I c R s. Generation of hyperchaos in the RCLSJJ array To investigate hyperchaotic generation in the RCLSJJ array, we solve Eqs. () () numerically. In the present work, the differential equations are solved using the fourth-order Runge Kutta method with a fixed step = and total steps N step =, the Lyapunov exponents are computed using the algorithms presented by Wolf [] and the voltage bifurcation diagrams are plotted by recording the local maxima in the voltage time series. [8] To illustrate the process of hyperchaotic generation, firstly, we consider case in which two RCLSJJs are in chaotic state before they are coupled and used in the array. The parameter values we used here are obtained from Refs. [] and [] at temperature T =. K and as follows: β C = β C = β C =., β L = β L = β L =.8, g = g = g = 8, and i =.. The initial values are γ =, υ =., i s =, γ =., υ =., and i s =.. The numerical results of Eqs. () (), i.e. the bifurcation diagrams and Lyapunov exponents λ and λ, are shown in Fig.. At the same time, we also compute other Lyapunov exponents λ λ and find that they are negative in the selected range of parameters, thus they are not presented in Fig.. From Figs. and, one can see that for.9 the array and two JJs used in this array are in chaotic states which are occasionally interrupted by periodic windows; for.9 <. two JJs are in the -period (P) state and the array is in -period (P) state. Here, near =.9, with value reducing, the array and two JJs enter into chaotic states through intermittency routes. From Fig. (c), one can find that for. Lyapunov exponents λ and λ are both positive, -
thus this array system is in a hyperchaotic state; for. < Lyapunov exponent λ is positive except 8 regions, i.e., 88 < <, < <, 8 < < 89, and near =.,.8,.8,.9,, 8;,,,, 8,, =,, 8,,, 9,, 9, 8, 9,,, and =, where the λ is near zero and λ is negative, thus this array system is in a chaotic state and occasionally interrupted by these 8 periodic windows, which also can be seen by locally enlarging Figs. and [corresponding local enlarged figures are not presented. vmax (vmax)...... state; for. < <.9 and. <.9, λ is positive and λ is negative, thus this array system is in chaotic states; for.9. and.9 <., λ is near zero and λ is negative, thus this array system and two JJs are in periodic states. Here, we compare the periodic motion behaviours in Figs. and. In the range of.9 <., the relationship of the oscillation characteristics of two junctions and this array will be illustrated. Now =. is chosen, the corresponding time series of phase differences γ and γ, and the time series of junction voltages υ, υ and array output voltage (υ + υ ) are shown in Figs. and, respectively. From Fig., one can see that two junctions are in the anti-phase case[9] in which two junctions have the same normalized oscillation period T but with staggered phases: γ ( ) = γ ( + (T /)). Thus, in Fig., the junction voltages υ and υ undergo the same T periodic oscillation but the time interval of peak values of υ and υ is T /, so the oscillation period of this array output voltage (υ + υ ) is also T /. Thus, for.9 <. both junctions are in P oscillation states and the array is in P oscillation states. T/ T. -.. λ λ v, v, v+v (c). λ, λ γ mod π γ, γ (v+v)max T/ T v v 8 v+v.. -. - -.. 8 Fig.. Time series of phase differences γ and γ ; time series of voltages v, v, and v + v. Fig.. Bifurcation diagram of two junctions versus ; the bifurcation diagram of the JJA versus ; (c) the Lyapunov exponents λ, λ versus. Secondly, we consider case in which two in this paper]; for <., λ is positive and λ is near zero, thus this array system is in chaotic RCLSJJs are in periodic states before being coupled. According to Refs. [] and [], we take βl = βl = -
βl =., and other parametric values are the same. <. and < two JJs are as those in Fig.. The numerical results of Eqs. () in P states, and the array is in P state. In the right (), i.e., the bifurcation diagrams and Lyapunov ex- end of bifurcation diagrams, with value decreasing, ponents λ and λ are shown in Fig.. Similarly, the the array and two JJs enter into chaotic states through other Lyapunov exponents λ λ are negative in the intermittency routes. From Fig. (c), one can find selected range of parameters and are not presented in that for <. Lyapunov exponents λ and Fig.. λ are both positive, thus this array system is in the. hyperchaos state; for. the λ is positive except regions, i.e., < <, vmax(vmax)..99 < < and near =.,.,.,..,.,.8,.,.,,,,. =,, and =, where λ is near. zero or negative and λ is negative or near zero, thus this array system is in chaos state occasionally in- terrupted by these periodic windows, which can..9 also be observed by locally enlarging Figs. and [corresponding local enlarged figures are not pre- (v+v)max. sented in this paper]. For < <., the λ is positive and λ is near zero, thus this array sys- tem is in a chaotic state. For. <.,.9 < <, and <, λ is near zero and λ is negative, thus this array system and two JJs are in periodic states. For..9 and near =, λ is positive and λ is negative,...9 thus this array system and two JJs are in chaos states. Similarly, to illustrate the relationship of periodic (c) λ λ oscillations of two junctions and this array, according to Figs. and, = is chosen. The corre- λ, λ sponding time series of phase differences γ and γ, as well as junction voltages υ, υ, and the array output - voltage υ + υ can be obtained by numerically solving Eqs. () (), and they are similar to Figs. and -.. Thus two junctions are also in anti-phase state -....9 with γ ( ) = γ ( + (T /)) To further show dynamical behaviours of different Fig.. Bifurcation diagrams of two junctions versus ; bifurcation diagram of JJA versus ; (c) Lyapunov exponents λ, λ versus. states in the array, according to Fig. (c), =,,. are chosen, the corresponding attractors, time series, and power spectra are presented in Fig., where From Figs. and, one can see that for figures (a) (a) show that this array is in hyper- <. and., the chaotic state; figures (b) (b) show that this array array and two JJs are in chaotic states, which are is in chaotic state; Figs. (c) (c) show that this occasionally interrupted by periodic windows. array is in P state. For -
(a) γ mod π (b) v+v - γ mod π 8 (c). v+v (c).. γ mod π - Frequency/ Hz (c) - -9 - - - - (b) - Power of (v+v) Frequency/ Hz - - - (b) (a) - 8 Power of (v+v) - - - Power of (v+v) (a) v+v 8 - Frequency/ Hz Fig.. Attractors on plane γ γ, time series of all voltage (v + v ) and power spectra of all voltage (v + v ) in different values. Panels (a) (a) for =, (b) (b) for =, and (c) (c) for =... Controlling hyperchaos in the RCLSJJ array Here we propose a scheme for controlling hyperchaos in the RCLSJJ array by adjusting the external bias dc current i based on the theory of nonfeedback control of chaos.[] To illustrate the process of controlling hyperchaos in this scheme, firstly let this array remain in the hyperchaotic state as shown 8 λ, λ (v+v)max in Figs. (a) (a). Secondly, numerical simulations are carried out by solving Eqs. () () by using i as a control parameter. The bifurcation diagram and Lyapunov exponents are presented in Fig., where one can see that exponent λ is near zero or negative in the selected range of parameters, thus hyperchaotic states disappear; for < i. the maximum Lyapunov exponent λ is near zero, thus this array is in periodic state. λ λ - -8 -. 8. i 8 i Fig.. Bifurcation diagram of the array versus i; Lyapunov exponents λ, λ versus i. -.
To further clarify that the array can be controlled into different period states, according to Fig., i = and i = are chosen. The corresponding attractors, time series and power spectra are shown in Fig., where Figs. (a) (a) and (b) (b) denote that this array is in P state and P state, γ mod π - 8 v+v γ mod π - - - - - (b) - - - - (a) Frequency/ Hz (b) (b) - Power of (v+v) Power of (v+v) (a) (a) v+v respectively. It is remarkable that the array is in hyperchaotic state as shown in Figs. (a) (a) before being controlled, thus the hyperchaotic state in this array can be controlled into different stable period states, such as P and P states as shown in Figs. and, respectively. 8 - Frequency/ Hz Fig.. Attractors on plane γ γ, time series of all voltage (v + v ), and power spectra of all voltage (v + v ) in different i values; panels (a) (a) for i =, and the array in P state; panels (b) (b) for i =, and the array in P state.. Conclusion References [] Benz S P and Buuroughs C 99 Appl. Phys. Lett. 8 In this paper we studied the dynamical characteristics in an array composed of two RCLSJJs and a shunted resistor by using bifurcation diagrams and [] Whan C B and Lobb C J 99 Phys. Rev. E [] Whan C B, Lobb C J and Forrester M G 99 J. Appl. Phys. 8 Lyapunov exponents. Numerical simulations demon- [] Cawthorne A B, Whan C B and Lobb C J 99 IEEE Trans. Appl. Supercond. 9 strated that hyperchaotic, chaotic and periodic states [] Whan C B and Lobb C J 99 Phys. Rev. E can coexist in this array by changing coupling inten- [] Cawthorne A B, Whan C B and Lobb C J 998 J. Appl. Phys. 8 sity, no matter whether two junctions are in chaotic states in advance or not. Moreover a scheme for controlling the hyperchaotic state in this array is presented by adjusting the external bias dc current. Numerical results confirmed that this scheme can be used to effectively control the hyperchaotic state in this array into stable periodic states and the different periodic states with different period numbers can be ob- [] Cawthorne A B, Barbara P, Shitov S V, Lobb C J, Wiesenfeld K and Zangwill A 999 Phys. Rev. B [8] Yang X S and Li Q Chaos, Solitons and Fractals [9] Dana S K, Roy P K, Sethia G C, Sen A and Sengupta D C IEEE Proc. Circ. Dev. Syst. [] Ucar A, Lonngren K E and Bai E W Chaos, Solitons and Fractals [] Feng Y L and Shen K 8 tained by choosing appropriate intensity of the exter- [] Feng Y L, Zhang X H, Jiang Z G and Shen K Int. J. Mod. Phys. B nal dc bias current. Here the hyperchaos state in the [] Feng Y L and Shen K 8 Eur. Phys. J. B JJA is reported for the first time. These results are [] Feng Y L and Shen K 8 useful for fabricating and applying a JJA composed of [] Feng Y L, Zhang X H, Jiang Z G and Shen K 9 Int. J. Mod. Phys. B RCLSJJs. [] Feng Y L and Shen K Chin. J. Phys. 8 -
[] Basler M, Krech W and Platov K Y 99 Phys. Rev. B [8] Hassel J, Gronberg L, Helisto P and Seppa H Appl. Phys. Lett. 89 [9] Grib A N, Seidel P and Scherbel J Phys. Rev. B 98 [] Bhagavatula R, Ebner C and Jayaprakash C 99 Phys. Rev. B 9 [] Chernikov A A and Schmidt G 99 Phys. Rev. E [] Zhou T G, Mao J, Liu T S, Lai Y and Yan S L 9 Chin. Phys. Lett. [] Klushin A M, He M, Yan S L and Klein N Appl. Phys. Lett. 89 [] Zhou T G, Yan S L, Zou X, Li S, Ji L and Zhao X J Chin. Phys. Lett. 9 [] Zhang S H and Shen K Chin. Phys. 9 [] Feng X Q, Yao Z H, Tian Z L and Han X Y Acta Phys. Sin. 9 8 (in Chinese) [] Wolf A, Swift J B, Swinney H L and Vastano J A 98 Physica D 8 [8] Sprott J C 99 Phys. Lett. A 8 [9] Hadley P and Beasley M R 98 Appl. Phys. Lett. [] Hu G, Xiao J H and Zheng Z G Chaos Control (Shanghai: Shanghai Scientific and Technological Education Publishing House) pp. -