Critical current density behaviors across a grain boundary inclined to current with different angles in YBa2Cu3O7−δ bicrystal junctions
Hua Tao1, 2, †, , Xu Wei-Wei2, Ji Zheng-Ming2, Guo Da-Yuan2, Wang Qing-Yun1, Ma Xiang-Rong1, Liang Rui-Yu1
School of Communication Engineering, Nanjing Institute of Technology, Nanjing 211167, China
Research Institute of Superconductor Electronics, School of Electronic Science and Engineering, Nanjing University, Nanjing 210093, China

 

† Corresponding author. E-mail: huataonju@163.com

Project supported by the National Natural Science Foundation of China (Grant Nos. 61501222, 61371036, and 61571219) and the School Scientific Research Fund of Nanjing Institute of Technology, China (Grant Nos. YKJ201418).

Abstract
Abstract

The critical current density behaviors across a bicrystal grain boundary (GB) inclined to the current direction with different angles in YBa2Cu3O7−δ bicrystal junctions in magnetic fields are investigated. There are two main reasons for the difference in critical current density in junctions at different GB inclined angles in the same magnetic field: (i) the GB plane area determines the current carrying cross section; (ii) the vortex motion dynamics at the GB affects the critical current value when the vortex starts to move along the GB by Lorentz force. Furthermore, the vortex motion in a bicrystal GB is studied by investigating transverse (Hall) and longitudinal current–voltage characteristics (IVxx and IVxy). It is found that the IVxx curve diverges from linearity at a high driving current, while the IVxy curve keeps nearly linear, which indicates the vortices inside the GB break out of the GB by Lorentz force.

PACS: 74.25.F–;74.62.Dh;74.25.Sv
1. Introduction

Grain boundaries (GBs) are typical crystalline defects present in high temperature superconductors (HTSs). By growing a film epitaxially on a bicrystalline substrate, a bicrystal GB in the film can be formed. The misorientation and the location of bicrystal GBs can be well controlled. Bicrystal GB Josephson junctions, formed at bicrystal GBs in superconducting films, are widely used in THz wave detection or superconducting quantum interference devices (SQUIDs). The response of YBa2Cu3O7−δ (YBCO) bicrystal Josephson junctions up to 4.25 THz at liquid nitrogen temperature (77 K) has been reported.[1]

In recent years, the fundamental property of vortex motion in HTS with GBs has become a hot topic of vortex dynamics.[29] The GB can be described as a collective line of artificial pinning defects, which suppresses the onset of vortex penetration and increases power dissipation. It has been demonstrated that the critical current density (Jc) across GBs in HTS decreases exponentially with the increase of misorientation angle of the GB.[10,11] However, the complex vortex motion behavior in the GB of HTS is still not well understood. In some work, it was found that GBs can act as easy-flow channels for longitudinal vortex motion,[1214] because the Lorentz force (FL) pushes vortices along the GBs when GBs are perpendicular to the current flow. In that case, Jc should decrease with increasing magnetic field. In contrast, some work shows that GBs act as barriers for vortex motion,[1518] since FL pushes vortices across the GBs transversally when the GBs are parallel to current flow, resulting in the enhancement of Jc. Moreover, lots of theoretical calculations and simulations[1922] suggest that the vortex pinning of the GB depends on the angle between the Lorentz force and the GB orientation. In this paper, we study the dependences of Jc on the magnetic field for YBCO bicrystal Josephson junctions with the GB inclined to the current direction at different angles. We also measure the longitudinal (along current direction) and transversal (vertical to current direction) IV characteristics of the YBCO film bridge with a GB at an inclined angle of 30° to study the vortex motion properties at the GB.

2. Sample and set-up

The YBCO films each with a thickness of 100 nm and a Jc of 3.9 MA/cm2 were deposited on MgO bicrystal substrates (c-axis (001) orientation) using the physical vapor deposition method. The bicrystal substrates were formed by hot pressing two single crystals ((001) orientation) each with a tilt angle of 24°. The following epitaxial growths of YBCO films on these substrates result in the formation of GBs in films each with a misorientation angle (θc) of 24°. The measured superconducting transition temperature (Tc) of YBCO film without a magnetic field is 87 K. Then, three 10-μm-wide bicrystal microbridges are patterned on the same bicrystal film (10 mm × 10 mm) by photolithography and ion beam milling. Here, the GB orientations were inclined to the current direction with angles (θ) of 90°, 60°, and 30° (see Fig. 1), respectively. For comparison, a microbridge of YBCO film onto MgO single-crystal substrate, i.e., no GB existed, was also fabricated. The geometry of this single-crystal microbridge was the same as those of the other three bicrystal microbridges.

Fig. 1. YBCO microbridges with a GB inclined to the current direction at θ = 90° (a), 60° (b), and 30° (c), respectively. (d) A microbridge of YBCO film without GB. The length and width of four microbridges are all 30 μm and 10 μm.

The samples were characterized using a physical property measurement system (PPMS-9). They were located in an adjustable DC mangnet with the direction of magnetic field (H) being perpendicular to the film plane (up to 8 T).

3. Experiments and discussion

We measure Jc values each as a function of H at 70 K for the YBCO microbridges with GBs at θ values of 90°, 60°, and 30° and without GB as illustrated in Fig. 2.

Fig. 2. Plots of critical current density of YBCO microbridges having GBs versus H with different θ values at a temperature of 70 K (left-hand longitudinal axis). The black squares show JcH characteristics of microbridge without GB (right-hand longitudinal axis).

When H = 0, the Jc of YBCO microbridge without a GB is 4.2 MA/cm2, while the Jc values of the other three microbridges with GBs are Jc(θ = 30°/0 T) = 45.0 kA/cm2, Jc(θ = 60°/0 T) = 20.0 kA/cm2, and Jc(θ = 90°/0 T) = 15.5 kA/cm2, respectively.

When H ≠ 0, the Jc values of four microbridges decrease with increasing H, and Jc (no GB) > Jc(θ = 30°) > Jc(θ = 60°) > Jc(θ = 90°). As H increases from 0 T to 1 T, Jc of the microbridge with a GB decreases sharply. When H > 1 T, the reduction rate of Jc is much lower. For the three microbridges with GBs, when H is below 1 T, the DC Josephson effect dominates, leading to a high reduction rate of Jc as H increases. When H is large, the intrinsic properties of YBCO GB play an important role, thus the drop of Jc slows down, and the properties of the microbridge with GBs are similar to those of a strong correlation superconductor. With increasing H, the Jc values of the three microbridges with GBs gradually approach to each other. As H > 7 T, they are nearly the same.

The first and obvious factor contributing to the increase of Jc with the decrease of θ for the three microbridges with GBs is the increase of the GB plane area, which determines the current carrying cross section. The GBs in these three samples are all inclined to the current direction at angles of 30°, 60°, and 90°, respectively, so the corresponding GB plane areas inside bridges are 2tw, 1.15tw, and 1tw, respectively, where t denotes the thickness of film and w is the width of microbridge.

If the GB plane area is the only factor determining critical current (Ic), Ic should be proportional to the GB plane area. In Fig. 3, samples A and B denote YBCO film microbridges with GBs at θ values of 90° and 30° respectively. The widths of the two microbridges are the same. Therefore, the GB plane area of A is tw, while the GB plane area of B is 2tw. Sample C denotes the microbridge with a GB at θ = 90°. The width of microbridge of C is twice that in A, and the GB plane area of C is 2tw. If the relationship of Ic and the GB plane area is linear, the Ic values of three samples should satisfy 2Ic(A) = Ic(B) = Ic(C). Correspondingly, 2Jc(A) = Jc(B) = 2Jc(C). However, the experimental results do not obey the above relationship, and Ic(B) > Ic(C) (the data is not shown here). Based on the data in Fig. 2, we can obtain Jc(θ = 30°/0 T)/Jc(θ = 90°/0 T) = 45.0/15.5 > 2, and Jc(θ = 30°/0.2 T)/Jc(θ = 90°/0.2 T) = 17.0/6.5 > 2. These results indicate that Jc is not proportional to the GB plane area when the width of microbridge is constant. Therefore, the increase of GB plane area is not the only reason for the increase of Jc as θ reduces.

Fig. 3. Diagrams of three different microbridges with GBs at different inclined angles or having different widths. The black line denotes the GB line.

Based on the above analysis, there should be other factors influencing the Jc of a microbridge with a GB. One possible factor is the vortex motion at the GB. Lots of vortices preferentially penetrate into GBs when a magnetic field is applied due to the high energy dissipation of GBs. The GB can be described as an easy-flow channel with a strong pinning effect on vortex motion. Even at 0 T, the self-field of current can create lots of Abrikosov vortices,[13,23] which appear at the GB first. Dissipation in current biased HTS film is caused by the vortex motion. In an HTS film microbridge with a GB, when the applied current is larger than Ic, FL drives vortices to move along the GB, and resistance appears, i.e., the superconducting state disappears. In the GB, the vortex is affected by FL and (the GB pinning force). The has two orthogonal components, i.e., FGB∥ and FGB⊥, which are schematically analyzed in Fig. 4. It is noted that is anisotropic. FGB⊥ is large while FGB∥ is much smaller than FGB⊥.[2] Despite that, the direction of total force F is parallel to the orientation of the GB. As the current increases, FL increases. The vortices do not move until FLcos α > FGB∥, where α denotes the angle between the GB and FL, and α = 90°–θ. Once vortices start to move along the GB, the superconducting state disappears. For the three microbridges having GBs, the FGB∥ values are the same. Therefore, FLcos α keeps constant when the vortices move along the GB. As α increases, or θ reduces, the critical value of the Lorentz force FL increases, resulting in the increase of Jc.

Fig. 4. Force analysis of a vortex in the GB. The vortex is affected by FL and , and the latter can be decomposed into FGB∥ and FGB⊥. Alphabet α denotes the angle between FL and GB.
4. Further experiments

To further investigate the vortex motion in bicrystal GBs, we measure the transverse and longitudinal IV characteristics of the YBCO microbridge with an inclined GB, each as a function of the magnetic field at 75 K. Here, the GB is inclined to the current direction (x) at an angle of 30°. The dimensions of YBCO film microbridge are 30 μm×10 μm×0.1 μm. Six electrodes patterned by photolithography and ion beam milling are used for electrical measurement as shown in Fig. 5(d), and are two leads for longitudinal voltage (Vxx) measurement, whereas and are leads for measuring the transversal voltage, i.e., Hall voltage (Vxy). The current applied along the x direction is slowly increased (10 μA/s) for data acquisition.

Fig. 5. IVxx (a) and IVxy (b) characteristics of a YBCO film bridge with an inclined GB at a angle of 30° in different magnetic fields at 75 K, (c) the close-up of panel (a) at a lower current. (d) A photo of the sample.

Figures 5(a) and 5(b) show that both Vxx and Vxy linearly increase with the increase of current when the current is lower than 13 mA, which indicates that vortex motion appears only inside the GB. When the applied current is higher than 13 mA, FL is larger, IVxx curves in magnetic fields above 0.5 T diverge from a linear relationship, and the transition current for diverging decreases with magnetic field increasing. Meanwhile, IVxy curves almost keep linear. At the transition point where the IVxx curve diverges from linearity, a few vortices inside the GB driven by Lorentz force break out of the GB, and mainly move along the direction of the Lorentz force (negative y-axis direction), which significantly enhances the Vxx value. However, it has little influence on Vxy. As the magnetic field increases, the vortex density increases. As a result, the driving current for the IVxx curve starting to diverge from linearity decreases, and the nonlinearity of the IVxx curve becomes more pronounced.

Figure 5(c) shows a detailed view of Fig. 5(a) at low applied current and 75 K. With increasing magnetic field, the critical driving force for the onset of the vortex motion in the GB becomes smaller. Correspondingly, Ic becomes smaller. Actually, the increase of magnetic field leads to the decrease of . Besides Josephson effects (the negative curvature of the IVxx curve), several step-like jumps are observed in the IVxx curve at 0.01 T, which is most likely due to the change of the vortex lattice in the GB. This phenomenon is often observed in a low magnetic field, and also in good agreement with the simulation results of Ref. [19].

5. Conclusions

In this work, we investigate the dependences of Jc on H for YBCO microbridge with GBs at different θ values at 70 K. The difference in Jc is mainly attributed to two reasons: (i) the different GB plane areas and (ii) the different vortex motion dynamics at the GB. Once the applied current goes beyond a critical value, the vortices pinned in the GB will move along the GB, and the superconducting state of YBCO film will disappear. We also investigate the vortex motions in a YBCO microbridge with an inclined GB at θ = 30° by characterizing its transverse and longitudinal transport property in magnetic fields. It is found that the IVxx curves diverge from linearity at large current and the IVxy curves almost keep linear, indicating that some vortices break out of the GB due to the drive of Lorentz force. Step-like jumps are observed in the IVxx curve at low magnetic field, which is most likely because of the change of the vortex lattice in the GB.

Reference
1Chen JHoriguchi HWang H BNakajima KYamashita TWu P H 2002 Supercond. Sci. Technol. 15 1680
2Hua TXu W WZhang YJiang LAn D YShi J XYu MWang Z HChen JKang LJing B BWang H BWu P H 2015 Supercond. Sci. Technol. 28 025005
3Gurevich ACiovati G 2008 Phys. Rev. 77 104501
4Ding F ZGu H WWang H YZhang H LZhang TQu FDong Z BZhou W W 2015 Chin. Phys. 24 057401
5Naqib S HIslam R S 2015 Chin. Phys. 24 017402
6Puig TGutiérrez JPomar ALlordés AGázquez JRicart SSandiumenge FObradors X 2008 Supercond. Sci. Technol. 21 034008
7Zhu B YDong J MXing D YZhao B R 1998 Phys. Rev. 57 5075
8Horide TMatsumoto KIchinose AMukaida MYoshida YHorii S 2007 Phys. Rev. 75 020504
9Feldmann D MHolesinger T GFeenstra RCantoni CZhang WRupich MLi XDurrell J HGurevich ALarbalestier D C 2007 J. Appl. Phys. 102 083912
10Hilgenkamp HMannhart 2002 Rev. Mod. Phys. 74 485
11Dimos DChaudhari PMannhart J1990Phys. Rev. B414038
12Hogg M JKahlmann FTarte E JBarber Z HEvetts J E 2001 Appl. Phys. Lett. 78 1433
13Jukna ABarboy IJung GBanerjee S SMyasoedov YPlausinaitiene VAbrutis ALi XWang DSobolewski R 2005 Appl. Phys. Lett. 87 192504
14Shlyk LNenkov KKrabbes GFuchs GMickel CRellinghaus B 2006 Supercond. Sci. Technol. 19 S472
15Horide TMatsumoto KYoshida YMukaida MIchinose AHorii S 2008 Phys. Rev. 77 132502
16Durrell J HFeldmann D MCantoni C 2007 Appl. Phys. Lett. 91 182506
17Faleski M CMarchetti M CMiddleton A A 1996 Phys. Rev. 54 12427
18Groth JReichhardt COlson C JField S BNori F1996Phys. Rev. Lett.213625
19Zhu B YXing D YDong J MZhao B R 1999 Physica 311 140
20Brass AJensen H J. 1989 Phys. Rev. 39 9587
21Kwok W KFendrich J AVinokur V MKoshelev A ECrabtree G W 1996 Phys. Rev. Lett. 76 4596
22Wang H YDing F ZGu H WZhang TPeng X Y 2014 Chin. Phys. 23 0107402
23Hilgenkamp HSchneider C WGoetz BSchulz R RSchmehl ABielefeldt H 1999 Supercond. Sci. Technol. 12 1043