# Full text of "Emission of continuous coherent terahertz waves with tunable frequency by intrinsic Josephson junctions"

## See other formats

Emission of continuous coherent terahertz waves with tunable frequency by intrinsic Josephson junctions Masashi Tachiki National Institute for Materials Science, 1-2-1 Sengen, Tsukuba 305-0047, Japan o o (N ll O o Or c3 I C o o > in o l> o ^t- o -t— > i C O o Mikio Iizuka, Kazuo Minami, Syogo Tejima, and Hisashi Nakamura Research Organization for Information Science and Technology, 2-2-54 Nakameguro, Meguro-ku, Tokyo 153-0061, Japan (Dated: February 2, 2008) We present a mechanism for emission of electromagnetic terahertz waves by simulation. High T c superconductors form naturally stacked Josephson junctions. When an external current and a magnetic field are applied to the sample, fluxon flow induces voltage. The voltage creates oscillating current through the Josephson effect and the current excites the Josephson plasma. The sample works as a cavity, and the input energy is stored in a form of standing wave of the Josephson plasma. A part of the energy is emitted as terahertz waves. PACS numbers: 74.50.+r, 74.25.Gz, 85.25.Cp Continuous coherent terahertz waves have various ap- plications in scientific field such as biology and informa- tion science. One of the hurdles for technological ad- vancements in the terahertz region of electromagnetic wave is the development of sources for intense and con- tinuous coherent terahertz waves. Therefore, we inves- tigate a new mechanism for emitting intense continuous and frequency tunable terahertz waves. In the high tem- perature superconductors, the strongly superconduct- ing Cu02 layers and insulating layers are alternatively stacked along the c-axis of the crystals and form a natu- rally multi-connected Josephson junction called intrinsic Josephson junction (I J J). In the IJJ there appears a new excitation wave called Josephson plasma, the frequency of which is in the range of terahertai*2,. The frequency appears in the region inside the superconducting energy gap and the Landau damping is very weak, and thus the excited plasma decays by emitting a terahertz electro- magnetic wave. For investigating an emission mechanism of terahertz electromagnetic wave from the IJJ, we use the following model shown by Figure 1. In Fig. 1 the IJJ is shown in green and the electrodes of a normal metal (for example gold) are shown in yellow. An external magnetic field B applied in the direction of the y-axis induces fluxons in the direction. The centers of fluxons are in the insu- lating layers. In this system, the superconducting and normal currents almost uniformly flow in the direction indicated by J in Fig.l. The fluxons flow in the direction of the x-axis with a velocity v and induce the flow volt- age in the direction of the z-axis. These voltages creates the oscillating Josephson current along the z-axis by the Josephson effect, when temperature is low enough below T c and the superconducting current is smaller than the superconducting depairing current along the c-axis. This oscillating current interacts strongly with the Josephson plasma due to the nonlinear nature of the system and intensively excites the Josephson plasma wave as shown later. We use Bi 2 Sr 2 CaCu208+5 that is appropriate in DC current Bi2Sr2CaCu208+5 [ GoltJ B y FIG. 1: Schematic diagram of a prototype model for ter- ahertz emission. Bi2Sr2CaCu20g+5 forms the IJJ shown in green, and electrode gold plates form the top and bottom electrodes shown in yellow. A dielectric waveguide shown in orange extends from the left surface of the IJJ. the experiments, and apply a magnetic field and external currents around J c the critical current to the IJJ. Then, the frequency of the plasma waves appears in the tera- hertz frequency range. The plasma wave is converted to an intense terahertz electromagnetic wave in the waveg- uide (dielectric) shown in orange in Fig. 1. In accordance with the mechanism mentioned above, we now derive the equations for the simulation. The superconducting order parameter of the Zth layer is ex- pressed as ipi(r,t) — A;e lw ( r '*). Since superconductivity in the Cu02 layers is strong at low temperatures, we as- sume A/(r, t) to be spatially and time independent. In this system, the current along the z-axis perpendicular to the Cu02 plane is given by a sum of the Josephson, quasi-particle, and displacement currents by Jz,i+i,i (r, t) = J c sin ipi+ij (r, t) + a c E Z, 1 + 1,1 (r,t) 2 +^d t E z , l+hl (r,t), (1) where J c is the critical current in zero magnetic field, a c is the normal conductivity along the c-axis, and ^ z .i+i,i is the electric field between (Z+l)th and Ith layers along the z-axis. The ipi+i t i is the gauge invariant phase difference defined as 2tt r tpi + i : l(r,t) = (p l+1 (r,t) - (fi(r,t) - — / dzA z (r,z,t), 00 J (2) with the vector potential A z (r, z, t) and the flux unit 0o- For the superconducting current densities in the Cu02 plane, we use the generalized London equations, since the Ginzburg-Landau parameter is very large in cuprate su- perconductors. We insert Eqs. (f ) and (2) into Maxwell's equations taking account of the superconducting current in the afo-planes and the charging effect in the CuC>2 plane. Following the calculation procedure^ and after cumbersome calculations, we have (1 - (A^)[{d^ l+ u + + sin *k+i,l) +a S i(d t ,(p' l+1 -p' l )+P(p' l+1 -p' l ))] = + dy,t/ji +lt i, (3) s'(l - aA^)p[ = dt'^i+^i - tPu-i). (4) In the above equation, we use normalized units for length, time, and charge density in the Ith CuC>2 plane, respec- tively defined as i x j i j / AgCJgp; x =— ,t =oj p t,andp l = : , (5) Jc where p\ is the charge density in the Ith Cu02 layer. The plasma gap frequency is given by lu p = c/( v / £c'A c ), e c being the dielectric constant along the z-axis and A c being the magnetic field penetration depth from the be surface plane. The parameters in Eqs. (3) and (4) are defined as , K b £cH 2 , s 4vrcr c A c C = TT ' a = ~Zf' s = t-, a^c? /? = -JTT (6) bll oil A c Y t C C where the s and d are the thickness of the superconduct- ing and insulating layers, respectively. X a b is the London penetration depth of the superconducting layer, and p is the Debye screening length. The operator A^ 2 ) is defined as A( 2 )/ ; = f l+1 - 2f t + fi—\. We use Eqs. (3) and (4) in the IJJ and use Maxwell's equations in the dielectric to simulate the emission of terahertz wave. If Eqs. (3) and (4) are linearized, we obtain the Joseph- son plasma solution. The Josephson plasma is a compos- ite wave of the Josephson current and electromagnetic wave. The amplitudes of the Josephson current and the electric field are always parallel to the c-axis. Depending on the wave propagation directions parallel to the c-axis and the a-axis, there appear respectively the longitudinal and transverse plasma waves. The transverse and longi- tudinal waves have been observed by Tamasaku et al^, and Matsuda et al(&, respectively. Let us show how the plasma wave is converted to the electromagnetic wave at the interface between the IJJ and the dielectric. As mentioned before, the plasma waves have the transverse and longitudinal components. However, the electromagnetic wave has only the trans- verse component. Therefore, only the transverse plasma wave can convert into the electromagnetic wave at the in- terface. When we solve Eqs. (3) and (4), we impose the following boundary condition. To connect the Josephson plasma wave in the IJJ to the electromagnetic wave in the dielectric at the interface, we put the usual electromag- netic boundary condition; the electric and magnetic fields parallel to the interface are continuous at the interface. The electromagnetic wave in the dielectric is assumed to transmit freely to outer space at the end surface of the dielectric. One possible way to achieve this condition is to make a gradation of the dielectric constant near the surface. It is assumed that the opposite surface of the IJJ is exposed to vacuum. Keeping in mind that the IJJ is f^S^CaC^Os+a and the dielectric is MgO, we chose A Q i,=0.4//m, A c =200yum, s=3l, d=12l, n=0.6A, a=0.l, /3=0.01~0.05, and the number of layers=20~1000, and take the dielectric con- stants along the z-axis in the IJJ and that of the dielec- tric wave-guide to be £=10. We apply a magnetic field of lTesla along the y-axis. We supply an external current from the gold electrodes shown in yellow in Fig. 1. We change the normalized external current J/J c from 0.2 to 1.5 in step of 0.0125. The superconductivity in the top and bottom electrodes is assumed to penetrate to 0.075 fim. We note that for the external currents mentioned above, the superconducting current part of J is always less than J c . The length of the IJJ is taken to be 100/im along the x-axis and the length of the dielectric is taken to be 50/mi along the x-axis. For each external current, the time evolution is simulated until the system reaches a stationary state, that is normalized time t'=600, and the emission power is calculated after that time. We assume that the system is uniform along the y- axis and make two-dimensional calculation in the x-z plane. We use the finite difference method to perform the numerical simulation. The simulation uses very large sized nonlinear equations heretofore difficult to compute; for a simulation using 10 6 spatial cells in the x-z two- dimensional model, it would take two years to simulate one case of 10 s time-steps using a personal computer with a 2GHz processor. Therefore, we used the Earth Simula- tor that has a peak performance of 40 teraflops and we carried out this simulation in one day. We simulated a total of sixty cases. Therefore, in sixty days we carried out simulations that would have required one hundred 3 Dielectric Terahertz ^ T ■ t^ t EM wave ■ T O -100 J. > r-f- so "V" mmmw-i o 0.0 2.0 4.0 6.0 X(nm) 8.0 10.0 FIG. 2: A snapshot of moving fluxons and the oscillating electric fields of the standing plasma wave. -40 -30 Dielectric- : , -10 x(|xm) 1 2 3 — MJJ and twenty years using a personal computer. Since the electromagnetic wave is the transverse wave, only the node-less Josephson plasma wave along the z- axis (the transverse plasma) can be converted into the electromagnetic wave with the same frequency. The max- imum stacking number that fulfill the condition is deter- mined by the applied magnetic field, and it increases as the magnetic field and A a & increase. We found that the stacking number that fulfill the condition is 100 layers under the magnetic field of 1 Tesla and X a b of 0.4/zm. Figure 2 shows a snapshot of the moving fluxons and the oscillating electric field. In the figure the static and uniform electric field have been subtracted. The black circles show the centers of the fluxons. The fluxons are almost randomly distributed and move uniformly in the x-axis direction. The speed of the fluxons is about 4% of the speed of light in the current J/ J c =0.6, /3=0.02 and under the applied magnetic field of 1 Tesla. The magnetic field of a fluxon spreads in the range of A c along the ab- plane. In the magnetic field of 1 Tesla, 120 fluxons are contained in the range. Therefore, the magnetic field of the moving fluxons are uniform and thus the fluxon flow create a spatially almost uniform voltage in the direc- tion of the z-axis and the voltage induces the oscillating current due to the Josephson effect and the oscillating current excites the Josephson plasma. The plasma wave is reflected at the interface between the IJJ and the di- electric, and at the surface exposed to vacuum, and it forms a standing wave in the stationary state as shown in Fig. 2. In the figure red and yellow show the electric field directing to the z-axis direction, and blue shows the electric field directing to the -z axis direction. Figure 3 shows a snapshot of the electric and magnetic fields of the standing wave at the 50th insulating layer from the bottom surface in the direction z-axis. The standing wave of the electric field along the z-axis oscil- lates with 2.86 THz around a constant electric field for the normalized current current J/ J c =0.8 and for /3=0.02. FIG. 3: A snapshot of the plasma wave that is a standing wave in the IJJ and the electromagnetic wave that is a propagating wave in the dielectric. The red line shows the interface between the dielectric and the IJJ. The components of the electric and mag- netic fields parallel to the interface fulfill the electromag- netic continuity condition. In Fig. 3 the green curves indicate the magnetic and electric fields at the time of 1.2769ns after the system attains to the stationary state and the brown curves indicate the magnetic and electric fields at the time of 1.2772ns. The belly of an electric field oscillation with the largest amplitude is pinned at the interface, inducing the large amplitude of the electric and magnetic field oscillations in the dielectric. The ratio of the electric field to the magnetic field is equal to that of the electromagnetic wave at 4/im from the interface. The emission mechanism of the coherent terahertz elec- tromagnetic wave is similar to a laser mechanism. The IJJ itself works as a cavity. The energy is stored in the form of a standing Josephson plasma in the IJJ. A few percent of the energy is emitted as a terahertz electro- magnetic wave from the IJJ to the dielectric through the interface. In Fig. 4 we show the frequency spectra calculated by FFT analysis of the electromagnetic wave at a location 4/im from the interface between the IJJ and the dielectric for the case of (3=0.03. When we change the normalized current J/J c from 0.85 to 1.3, the frequencies of the elec- tromagnetic waves change from 2.08 THz to 3.21 THz. The frequency spectra are very sharp. The frequency of the electromagnetic wave in the IJJ varies continuously by changing current J / J c , indicating that emission of the sample is frequency tunable. Some peaks have satellites due to the frequency modulation of the excited plasma wave. The frequencies calculated by the simulation are equal to the AC Josephson frequency estimated by the 4 2.5 3 Frequency (THz) FIG. 4: The frequency spectra analyzed by FFT for 10 nor- malized currents. The intensities are those of electromagnetic waves in the dielectric. The intensity is obtained by calculat- ing the Poynting vector at a location of 4/xm from the inter- face. 1 5 a I" D. 0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 1.4 Normalized Current (J/Jc) along the y-axis is 500/im. The emission occurs in some limited current region as seen in Fig5. When the current is smaller than a critical current, the emission power van- ishes. The reason is that nodes in the Josephson plasma appear along the z-axis, indicating inclusion of the lon- gitudinal component, and thus the emission power van- ishes as mentioned before. When the current is larger than the upper critical value, the fluxons take various dynamical orders 6 . In this state, the emission vanishes since the rigid fluxon lattice flow prohibits occurrence of the standing wave. The sharp dip of power appear- ing at J/J c =0.7 for (3=0.02 comes from the following reason. The fluxons flow with a small amplitude vibra- tion. When the Josephson frequency resonates with the fluxon vibration frequency and this resonance stabilizes the Josephson plasma, it depresses the conversion rate from the Josephson plasma to the electromagnetic wave. When (3 decreases, for example (3=0.01, the oscillations of electric and magnetic fields lose the coherence and the emission power vanishes. When (3 increases, the dumping becomes larger and the emission power decreases. In the present paper we performed the simulation of terahertz wave emission with the parameters correspond- ing to Bi2Sr2CaCu208+<5, and obtained the following re- sults: The sample works as a cavity, and the input energy is stored in a form of standing wave of the Josephson plasma. A part of the energy is emitted as terahertz waves. The emission of electromagnetic waves is in the form of continuous coherent terahertz waves with mW power in the present model. The frequency is tunable by changing the applied current. FIG. 5: The emission power (mW) of the electromagnetic wave as a function of the external current. Acknowledgments fluxon flow voltage within the accuracy of 1.3%. This fact suggests that the plasma wave is excided by the os- cillating current through the AC Josephson effect. Figure 5 shows the emission power (the Poynting vec- tor) measured at the location 4/im from the interface as a function of the external current. The length of the sample We thank Dr. T. Sato to permit us for using the Earth Simulator and Dr. Y. Kanada for the benefit of our use of SR8000. We thank Dr. T. Yamashita, Dr. T. Hatano, Dr. Y. Takano, Dr. K. Hirata, Dr. K. Kadowaki, Dr. H. Wang, Dr. S. Kim, Dr. T. Tachiki, Dr. T. Koyama, Dr. M. Machida, H. Matsumoto, and G. Prabhakar for valuable discussions. 1 M. Tachiki, T. Koyama, and S. Takahashi, Phys. Rev. B 50, 7065 (1994). 2 L. N. Bulaevskii, M. P. Maley, and M. Tachiki, Phys. Rev. Lett 74, 801 (1995). 3 For example, see M. Tachiki, T. Koyama, and Taka- hashi, in Coherence in high temperature superconductor, G. Deutscher, and A. Revcolevschi, Eds, (World Scientific, Singapore, 1996), pp371. 4 K. Tamasaku, Y. Nakamura, and S. Uchida, Phys. Rev. Lett 69, 1455 (1992). 5 Y. Matsuda, et al. Phys. Rev. Lett 75, 4512 (1995). A. E. Koshelev and I. Aranson, Phys. Rev. B 64, 174508 (2001).