arxiv: v1 [cond-mat.supr-con] 2 Jan 2019

Transcription

1 Scaling Behavior of Quasi-One-Dimensional Vortex Avalanches in Superconducting Films A. J. Qviller nsolution AS, Maries gt. 6, 0368 Oslo, Norway T. Qureishy and P. Mikheenko Department of Physics, University of Oslo, P. O. Box 1048 Blindern, 0316 Oslo, Norway arxiv: v1 [cond-mat.supr-con] 2 Jan 2019 Y. Xu The Key Laboratory of Advanced Functional Materials, Ministry of Education, Beijing University of Technology, Beijing , China and Department of Energy Conversion and Storage, Technical University of Denmark, Roskilde 4000, Denmark H. Suo The Key Laboratory of Advanced Functional Materials, Ministry of Education, Beijing University of Technology, Beijing , China P. B. Mozhaev Institute of Physics and Technology of the Russian Academy of Sciences, Moscow, , Russia J. B. Hansen Department of Physics, Technical University of Denmark, Kongens Lyngby, DK-2800, Denmark T. H. Johansen Department of Physics, University of Oslo, P. O. Box 1048 Blindern, 0316 Oslo, Norway and Institute for Superconducting and Electronic Materials, University of Wollongong, Northfields Avenue, Wollongong, NSW 2522, Australia J. I. Vestgården Department of Physics, University of Oslo, P. O. Box 1048 Blindern, 0316 Oslo, Norway and Norwegian Defence Research Establishment (FFI), P. O. Box 25, 2027 Kjeller, Norway (Dated: January 3, 2019) Scaling behaviour of dynamically driven vortex avalanches in superconducting YBa 2Cu 3O 7 δ films deposited on tilted crystalline substrates has been observed using quantitative magneto-optical imaging. Two films with different tilt angles are characterized by the probability distributions of avalanche length, area, and sizes in terms of numbers of vortices. It was found in both samples that these distributions follow power laws over up to three decades, and have exponents ranging between 1.0 and 1.4. The distributions also show clear finite-size scaling, when the system size is defined by the depth of the flux penetration front a signature of self-organized criticality. A scaling relation between the avalanche size exponent and the corresponding fractal dimension, previously derived in the literature and shown in numerics, is also found to be satisfied. Qualitative measurements on an additional sample covered with a gold layer give evidence that dynamically driven vortex avalanches are not nucleated by thermomagnetic instabilities. I. INTRODUCTION Avalanche behaviour is found in a wide range of natural systems, and is commonly observed, e.g., as abrupt displacement in granular media, earthquakes, Barkhausen noise caused by sudden motion of magnetic domain walls in ferromagnetic materials, and abrupt displacement of vortices in superconductors triggering magnetic flux avalanches [1 3]. Such avalanche activity is in general highly unwanted in practical applications of superconductors. However, its origin, which lies in the competition between intervortex interactions and quenched disorder, makes vortex matter an interesting system for studies of non-equilibrium dynamics. Power laws in avalanche sizes, avalanche durations and temporal power spectra, together with finite-size scaling, are often interpreted as manifestations of self-organized criticality (SOC). This framework was originally developed from studies of a cellular automaton model of a sandpile [4]. Although SOC seems to capture basic aspects of sandpile physics, both rolling and inertial effects tend to create system-spanning avalanches, which introduce cutoffs in the power law probability distributions already after approximately one decade [1]. The same effect is seen also in results obtained by numerical simulations [5]. In contrast, studies of piles of rice grains have demonstrated power law probability behaviour over nearly four decades [6 8]. This increase was explained to bearesultofricegrainshavingelongatedshapesandalso

2 2 being lighter, reducing both rolling and inertial effects. In type-ii superconductors an applied magnetic field is penetrating the material in the form of quantized vortices. These lines, each carrying one quantum of magnetic flux, Φ 0 = h/2e, with h being Planck s constant and e the elementary charge, have zero interia. Moreover, also rolling effects are here nonexistent. Thus, the vortex matter represents a highly favourable system for studies of avalanche dynamics and SOC. The vortices are affected by the Lorentz force from electrical currents flowing inside the superconductor. However, this force can be counterbalanced by forces from microscopic material defects, which then serve as pinning sitesforthevortexmotion. Thus, inaslabgeometry,vortices build up a metastable state [9] with a density profile quite analogous to the shape of a sandpile [10]. Just as a sandpile with constant slope is uniquely characterized by its angle of repose, the flux pile in a slab geometry has a constant flux density gradient, which by a Maxwell equation implies a current flow of constant magnitude, the so-called critical current density, j c. For thin films in perpendicular magnetic fields, the analogy is less obvious [11, 12] as the flux density gradient here is not constant, although the critical state is still characterized by j c. Flux avalanches have been reported to occur in a wide variety of superconductors including Nb, Pb, Nb 3 Sn, NbN, MgB 2, YNi 2 B 2 C, MoGe, MoSi and YBa 2 Cu 3 O 7 δ (YBCO) [13 24]. The avalanches may be of thermomagnetic origin or be caused by dynamically driven vortex rearrangements [3, 23]. Both kinds of avalanches can occur in the same sample [25, 26]. There is extensive experimental and numerical evidence that under some circumstances, quantities measuring the size of dynamically driven vortex avalanches are distributed according to power laws and may also exibit finite-size scaling [3, 22, 23, 27, 28]. For a recent review of SOC experiments in granular media and superconductors, see ref. [29]. In the present work, quantitative magneto-optical imaging (MOI) was used to measure the probability distributions of avalanche length, area and size in two YBCO films deposited on tilted substrates with different tilt angles. It is shown that the probability distributions of avalanche length and areal extent, as well as size in term of number of vortices, all follow power laws. Additionally, these distributions obey finite-size scaling for both samples. The avalanche size exponent and its corresponding fractal dimension is found to satisfy a scaling relation. Finally, dynamically driven vortex avalanches are observed in a film covered with a gold layer. This is in contrast to thermomagnetic avalanches, whose nucleation in thin films is generally suppressed by metal coatings [30]. II. EXPERIMENTAL When YBCO is deposited epitaxially on slightly tilted (vicinal) substrates, planar defects in the form of antiphase boundaries are introduced with a period of 2-5 nm [31]. For tilt angles of θ 10, grain alignment and current-carrying abilities of the films are improved, thus making them interesting for power applications [31, 32]. Such films also have potential applications in Josephson junction circuits [33 35]. YBCO films on tilted substrates also exhibit other extended planar defects due to lattice mismatches. This results in anisotropic flux penetration and different critical currents in the direction parallel and normal to the defects. The extended defects also facilitate easy vortex motion [31, 32, 36], and at low temperatures, the flux penetration becomes strongly quasi-one-dimensional, forming straight channels. For this work, one sample with the shape of a strip of 0.9 mm width was deposited by laser ablation on a NdGaO 3 substrate with a tilt angle θ = 14. Another sample shaped as a 4x4 mm 2 square was deposited by spin coating on a LaAlO 3 substrate with a tilt angle of θ = 20. Both films are 200 nm thick and have critical temperatures, T c s, of 88 K and 90 K, respectively. For sample preparation details, see [35] and [37], respectively (the θ = 20 sample is prepared like sample B in ref. [37]). A third YBCO sample with θ = 10 and the same substrate, thickness and T c as the 14 -sample was covered with gold by sputtering and used for qualitative measurements. The flux penetration in all three samples is strongly intermittent at low temperatures in the form of quasi-1d avalanches along the channels. The samples described above have different j c s, and thus, two slightly different experimental procedures were used to investigate the intermittent flux penetration into them. Both samples were first zero-field-cooled (zfc) to T = 4 K. The 14 -sample was then subjected to a field ramptob a =17.0mTin400steps, whilethe20 -sample was subjected to a field ramp to B a = 8.5 mt in 200 steps. Thus, the same field step B a = 42.5 µt was used in both cases. Five magneto-optical images were collected at each field, after waiting 5 seconds for the vortex medium to relax. After reaching the maximum field in the ramp, B a was set to zero and the temperature raised above T c in order to anneal the samples, and ensure they were in a virgin state before the next field ramp. In total, image series from 10 ramps were collected for the 14 -sample, and from 20 ramps for the 20 -sample. These image series were then averaged over the five images taken at each field to reduce the noise. Qualitativemeasurementsonthegold-coated10 -sample were made after zfc to T = 4 K and subjecting it to a field ramp to B a = 25.5 mt in 600 steps. MOI was performed using an in-plane magnetized ferrite-garnet film as Faraday rotating sensor [38, 39]. This technique allows realtime visualization of flux patterns in superconductors, and allows high spatial (µm) and temporal(picosecond) resolution. Moreover, by MOI the entire flux pile can be seen, and thus, the statistics of internal avalanches and their morphology be precisely measured. The Faraday-rotating indicator film has a non-linear response to magnetic field and the lamp illuminationinthe setupisalsouneven. In ordertocalibrate

3 3 FIG. 1. (a) Magneto-optical image of flux penetration in the 14 -sample at T = 4 K and B a = 15.0 mt. (b) Differential image at the same temperature and fieldas in(a) with B a = 42.5 µt. An avalanche of length 170 µm, area 3700 µm 2 and size 1280 Φ 0 is marked with the red box. (c) Differential image of the gold-coated 10 -sample at T = 4 K, B a = 23.7 mt, and B a = 42.5 µt. All scale bars are 150 µm long. the image series, an additional ramp from zero up B a = 17.0 mt was performed a few Kelvin above T c for both samples. The light response I(x, y) of the indicator film was fitted by a polynomial in B a (x,y) for every pixel of the CCD chip of the camera. III. RESULTS Presented in Fig. 1 (a) is a calibrated magneto-optical image of the flux density distribution near the long edge of the 14 -sample at T = 4 K. The edge itself is seen here as the bright horizontal line, caused by the external magnetic field piling up along the rim of the diamagnetic sample. One can see a finger-like pattern of easy flux penetration, corresponding to the extended defects induced by the substrate tilt. Flux penetration in the form of quasi-one-dimensional avalanches inside these channels is visible as bright lines in Fig. 1 (b), where a differential image of the same sample as in Fig. 1 (a) is shown. This differential image was obtained by taking the next image in the series after the ones shown in Fig. 1 (a) and subtracting the latter from the former, thus it shows the change in the flux penetration pattern after an increase of the applied magnetic field by B a. In Fig. 1 (c), a differential image of the gold-coated 10 -sample recorded at T = 4 K, B a = 23.7 mt, and B a = 42.5 µt is shown. Vortex avalanches are clearly visible as stripes lighter than the background, which shows that the gold layer is not effective at suppressing the dynamic vortex avalanches. The image is less sharp than the differential image in Fig. 1 (b) as the gold layer introduces a greater distance between the magneto-optical indicator film and the sample. A computer program was used to extract the probability distributions of the flux avalanches from the differential images throughout all the quantitative series, segmented into B a -intervals corresponding to 40 differential images in the case of the 14 -sample and 20 differential images in the case of the 20 -sample. Overall 400 differentialimagesfor both samplesper B a -intervalwereavailable for statistics. In order to ensure that a proper critical state was built up before measurements, the first 120 images were discarded for the 14 -sample and the first 80 for the 20 -sample in every ramp. A threshold value of brightness was used to separate vortex avalanches from background in the differential images, and a mean filter was used to remove noise. The same threshold brightness and mean filter were used for both samples, resulting in directly comparable probability distributions for the measured parameters: avalanche length, avalanche area and size of the observed avalanches. Avalanche length is defined as the length of the avalanches in the longitudinal direction of the channels (vertical in Figs. 1 (a) and (b)), avalanche area as the total area of pixels above the threshold brightness, and avalanche size as the number of vortices in a given avalanche, calculated by dividing the amount of flux in an avalanche by the magnitude of one flux quantum Φ 0 = Wb. The absolute frequencies of these parameters were binned in histograms and each bin was subsequently divided by the total number of identified avalanches to yield probability distributions. Shown in Figs. 2 (a), (b) and (c) are the probability distributions of length, area and size for the 14 -sample, respectively, while the corresponding probability distributions for the 20 -sample are shown in Figs. 3 (a), (b) and (c). The different colour-coded graphs correspond todifferent B a -intervalsas indicated in the figures. The length distributions show power-law behaviour over about one and a half decade, while the size and area dis-

4 4 10 (a) 0 10 (a) 0 P(l) (b) l ( m) P(l) (b) l ( m) P(a) (c) a ( m 2 ) P(a) (c) a ( m 2 ) P(s) s ( ) 0 P(s) s ( ) FIG. 2. Probability distributions of (a) avalanche lengths l, (b) avalanche areas a and (c) avalanche sizes s in the 14 - sample. FIG. 3. Probability distributions of (a) avalanche lengths l, (b) avalanche areas a and (c) avalanche sizes s in the 20 - sample. tributions show such behaviour over about three decades for both samples. The smallest avalanches consist of about 10 flux quanta, whereas the largest are rearrange- ments of vortices or more. The distributions in Figs. 2 and 3 look like approximate power laws with cutoffs. We will thus investigate

5 5 if the functions P(l,L), P(a,L) and P(s,L) are able to capture the functional forms of the distributions of length l, area a and size s, respectively. The functions are described as ( ) l P(l,L) l ρ f L D, (1) ( a ) P(a,L) a σ g L D, (2) ( s ) P(s,L) s τ h L D. (3) Here, ρ, σ and τ are the avalanche exponents for the length, area and size distributions. The scaling functions f, g and h provide cutoffs corresponding to the system size L raised to the power of the fractal dimension D of the avalanches. In this system, L is identified, and manually measured, as depth of the flux front in the middle of each investigated B a -interval. In order to check for finite-size scaling of the distributions measured at different B a -intervals and thus corresponding to different system sizes, one starts by rearranging Eqs. 1, 2 and 3, such that scaling functions f, g and h are isolated on the right-hand side. For example, in the case of Eq. 3, one plots s τ P(s,L) against s/l D to investigate the functional form of h. If the distributions collapse onto one curve when this is done for all the B a -intervals, finite-size scaling is present and the better the curve collapse, the more exact is the scaling. Such curve collapses are shown in Figs. 4-5 (a), (b), and (c) for the distributions in Figs. 2-3 (a), (b), and (c). The figures show expected deviations from power law behaviour and scaling near the lower cutoff. Furthermore, the avalanche exponents and fractal dimensions may be determined from careful adjustment and inspection of the curve collapses. When the avalanche exponent is adjusted to its correct value, a flat plateau will appear in the collapsed curves corresponding to the power law region between the lower and upper cutoffs of the distribution. The fractal dimension is, on the other hand, determined by adjustment until the upper cutoffs are aligned as well as possible. This was done for all the probability distributions measured in this work and the results are stated in Table I. For the length distributions, the linear depth of the flux front envelope is the obvious cutoff parameter and D was therefore fixed to 1. The validity of the data analysis is supported by considering the flux conservation scaling relation D(2 τ) = 1 (4) This equation was derived from flux conservation and shown to be satisfied numerically in [27]. See [40] for a statement on its validity in the experimental context of this work. For the 14 -sample, good curve collapses are obtained in the range from D = 1.35 to D = 1.5. TABLE I. Avalanche exponents ρ, σ, τ and fractal dimensions D. Sample Parameter ρ, σ, τ 14 Length 1.35 (±0.03) 1 14 Area 1.1 (±0.1) 1.4 (±0.1) 14 Size 1.30 (±0.03) 1.4 (+0.1, 0.05) 20 Length 1.1 (±0.1) 1 20 Area 1.2 (±0.1) 1.4 (±0.1) 20 Size 1.06 (±0.03) 1.4 (±0.1) Insertingthe values D = 1.35, τ = 1.33and D = 1.5, τ = 1.27, correspondingtothe uncertaintiesofd and τ in the left-handsideofeq. 4, givesd(2 τ)=0.9045and1.095, respectively. Thus, this relation seems to be obeyed with an accuracy of approximately 10 %. In the case of the 20 -sample, the agreement is not so good, D = 1.3, τ = 1.09 and D = 1.5, τ = 1.03 gives D(2 τ) = and 1.455, respectively. This is most likely because its avalanche size distribution shows some deviations from a pure power law. Eq. 4 is not expected to be valid for avalanche length or area distributions, as it is derived from integration of the size distribution. D IV. DISCUSSION It is clear that the probability distributions of avalanche lengths, areas and sizes demonstrate power law behaviour and finite-size scaling consistent with SOC in both samples, but it is expressed more clearly in the 14 - sample than in the 20 -sample. This can be explained by the more regular and ordered channel structure of the former sample in comparison with the latter one and thus by the introduction of multiple length scales in the latter sample, which renders the scale invariance only approximate. Molecular dynamics simulations have also shown a similar breakdown to dirty power laws when length scales, in addition to the system size, are introduced by low density of pinning sites and the resulting formation of vortex rivers [41]. For the length distributions in this work, power law behaviour and scaling are observed over one decade of value, while the various area and size distributions show this behaviour and scaling over two to almost three decades. The area distributions, however, show the greatest deviation from the pure power laws. Other MOI experiments such as ref. [22], with τ = 1.29 ± 0.02 and D = 1.89 ± 0.03, and ref. [28] with τ = 1.07 ± 0.02 and D = 2.25 ± 0.05, give similar avalanche size exponents, but considerably different fractal dimensions in comparison with this work. The famous experiment of Field et al. [42] utilizing a hollow-cylinder geometry, gave a range of exponents τ = , but the exponents found in off-edge avalanches are not necessarily directly comparable with internal avalanches [29]. Values differing considerably from this work were found in [26] with τ = 2.05 and [43] with τ = 3.0, but these experiments had single-vortex resolution and probed a

6 (a) 10 1 (a) 6 l P(l,L) (b) l/l l P(l,L) (b) l/l a P(a,L) (c) a/l D a P(a,L) (c) a/l D s P(s,L) 10 1 s/l D FIG. 4. Finite-size scaling of (a) avalanche length probability distributions P(l, L), (b) avalanche area probability distributions P(a, L) and (c) avalanche size probability distributions P(s,L) in the 14 -sample. s P(s,L) 10 1 s/l D FIG. 5. Finite-size scaling of (a) avalanche length probability distributions P(l, L), (b) avalanche area probability distributions P(a, L) and (c) avalanche size probability distributions P(s,L) in the 20 -sample. different range of the probability distribution functions than that covered by MOI experiments (from tens to tens of thousands of vortices). Molecular dynamics simulations give exponents in the rangeofτ = for the numberofdisplaced vortices (avalanche size) in systems with high pinning density, where the bigger exponent corresponds to higher pinning strength [41]. The avalanche size exponents found in this

7 7 work are also not very different from that in the Bassler- Paczuski model, a cellular automaton for 2D flux penetration [27] with τ = 1.63±0.02 and D = 2.75±0.1 and even closer to a variation of it with strong, periodic and densely spaced pinning sites [44], in which τ = 1.45±0.02 and D = 2.2 ± 0.1. The Bassler-Paczuski model has been shown describing robust SOC behaviour over four decades for a variety of different parameters in the pinning landscape [27, 44]. However, both the molecular dynamics simulations and the cellular automata models assume short-ranged vortex-vortex interactions for the bulk state, whereas in thin films the vortex-vortex interactions fall off as 1/r and probably result in a different universality class. The intermittent flux penetration seen in the investigated YBCO samples differs from thermomagnetic instabilities in the sense that it neither has threshold fields nor threshold temperatures [45, 46], but rather appears gradually when the sample is cooled down and is seen for all accessible fields (0-85 mt). Additionally, thermomagnetic microscopic avalanches were analyzed in [17] from both a theoretical and experimental viewpoint and were found to have a field-dependent characteristic size and the same threshold temperature as the macroscopic, dendritic flux avalanches. The avalanche lengths of the intermittent flux penetration are also limited by the envelope of the flux front [23], which is not the case for thermomagnetic instabilities. Furthermore, nucleation of thermomagnetic instabilities is suppressed by metal coatings, and thus the measurements on an additional sample covered with a gold layer give evidence that dynamically driven vortex avalanches are not nucleated by thermomagnetic instabilities. Fast flux avalanches induce voltage pulses when they propagate under metal coatings, which allows their speed to be measured from a combination of MOI and voltage pulse measurements [47]. This interaction with metal coatings also causes fast flux avalanches to be refracted [48]. No signs of avalanche suppression, induced voltage pulses or refraction were seen for the intermittent flux penetration in the gold-coated 10 -sample, which was investigated by the same method as in ref. [47]. Unpublished work on the NbN sample investigated in ref. [48] also shows intermittent flux penetration with no interaction with the metal coating, easily visible above the threshold temperature of thermomagnetic instability present in the sample. This provides evidence that the intermittent flux penetration is relatively slow. V. CONCLUSIONS In summary, the probability distributions of flux avalanches were measured by quantitative MOI for two YBCO samples deposited on substrates cut with tilt angles of 14 and 20. Probability distributions of avalanche lengths, areas and sizes in terms of numbers of vortices were extracted from the data. These parameters followed approximate power laws in distributions over up to three decades and demonstrated finite-size scaling. Avalanche exponents and fractal dimensions were obtained by careful inspection of the collapse of finite-size scaling curves. The obtained exponents were between 1.0 and 1.4, and the avalanche size exponents determined in this work are similar to those that has been found in other MOI experiments on superconductors, molecular dynamics simulations and cellular automata models. For bothsamples, thefractaldimensiond wassetto1forthe avalanche length distribution, while the best collapse of the avalanche area and size distributions were obtained with D = 1.4. The scaling relation D(2 τ) = 1 between the avalanche size exponent and its corresponding fractal dimension, previously derived theoretically from flux conservation and shown to be satisfied in numerics, was also experimentally validated with an accuracy of 10% for the 14 sample. Qualitative work on an additional gold-coated YBCO sample with tilt angle of 10 showed no interaction between the intermittent flux penetration and the gold coating, providing further evidence that these avalanches are not nucleated by thermomagnetic instabilities and move relatively slow. VI. ACKNOWLEDGMENTS The Authors express their gratitude to Berit Løken Berg at the Department of Geosciences, University of Oslo for assistance with gold deposition. The authors also want to thank Julia E. Mozhaeva for her kind assistance during sample preparation. This work was financially supported by the National Natural Science Foundation of China ( ) and Beijing Natural Science Foundation ( ). REFERENCES [1] H. J. Jensen, Self-Organized Criticality Emergent Complex Behavior in Physical and Biological Systems, Cambridge University Press, pp (1998). [2] G. Pruessner, Self-Organised Criticality: Theory, Models and Characterisation, Cambridge University Press (2012). [3] E. Altshuler and T. H. Johansen, Colloquium: Experiments in vortex avalanches, Rev. Mod. Phys. 76, 471 (2004). [4] P. Bak, C. Tang and K. Wiesenfeld, Self-organized criticality, Phys. Rev. A 38, 364, (1988). [5] C. Prado and Z. Olami, Inertia and break of selforganized criticality in sandpile cellular-automata models, Phys. Rev. A 45, 665 (1992). [6] V. Frette, K. Christensen, A. Malthe-Sørensen, J. Feder, T. Jøssang and P. Meakin, Avalanche dynamics in a pile

8 8 of rice Nature 379, 49 (1996). [7] C. M. Aegerter, R. Günther and R. J. Wijngaarden, Avalanche dynamics, surface roughening, and self-organized criticality: Experiments on a threedimensional pile of rice, Phys. Rev. E 67, (2003). [8] K. A. Lörincz and R. J. Wijngaarden, Edge effect on the power law distribution of granular avalanches, Phys. Rev. E 76, (R) (2007). [9] C. P. Bean, Magnetization of High-Field Superconductors, Rev. Mod. Phys. 36, 31 (1964). [10] de Gennes, P. G., Superconductivity of Metals and Alloys, W. A. Benjamin, New York (1966). [11] E. H. Brandt and M. Indenbom, Type-IIsuperconductor strip with current in a perpendicular magnetic field, Phys. Rev. B 48, (1993). [12] E. Zeldov, J. R. Clem, M. McElfresh and M. Darwin, Magnetization and transport currents in thin superconducting films, Phys. Rev. B 49, 9802 (1994). [13] C. A. Duran, P. L. Gammel, R. E. Miller and D. J. Bishop, Observation of magnetic-field penetration via dendritic growth in superconducting niobium films, Phys. Rev. B 52, 75 (1995). [14] M. Menghini, R. J. Wijngaarden, A. V. Silhanek, S. Raedts and V. V. Moshchalkov, Dendritic flux penetration in Pb films with a periodic array of antidots, Phys. Rev. B 71, (2005). [15] I. A. Rudnev, S. V. Antonenko, D. V. Shantsev, T. H. Johansen and A. E. Primenko, Dendritic flux avalanches in superconducting Nb 3Sn films, Cryogenics 43, 663 (2003). [16] I. A. Rudnev, D. V. Shantsev, T. H. Johansen and A. E. Primenko, Avalanche-driven fractal flux distributions in NbN superconducting films, Appl. Phys. Lett. 87, (2005). [17] D. V. Shantsev, A. V. Bobyl, Y. M. Galperin, T. H. Johansen, and S. I. Lee, Size of flux jumps in superconducting films, Phys. Rev. B 72, (2005). [18] S. C. Wimbush, B. Holzapfel and Ch. Jooss, Magnetic instability in YNi 2B 2C thin films: Observation of dendritic flux instabilities in borocarbides, J. Appl. Phys. 96, 3589 (2004). [19] M. Motta, F. Colauto, J. I. Vestgården, J. Fritzche, M. Timmermans, J. Cuppens, C. Attanasio, C. Cirillo, V. V. Moschalkov, J. Van de Vondel, T. H. Johansen, W. A. Ortiz, A. V. Silhanek, Controllable morphology of flux avalanches in microstructured superconductors, Phys. Rev. B 89, (2014). [20] F. Colauto, M. Motta, A. Palau, M. G. Blamire, T. H. Johansen and W. A. Ortiz, First Observation of Flux Avalanches in a-mosi Superconducting Thin Films, IEEE Trans. on Appl. Supercond., 25, , (2015). [21] P. Leiderer, J. Boneberg, P. Brüll, V. Bujok, and S. Herminghaus, Nucleation and growth of a flux instability in superconducting YBa 2Cu 3O 7 x films, Phys. Rev. Lett. 71, 2646 (1993). [22] C. M. Aegerter, M. S. Welling and R. J. Wijngaarden, Self-organized criticality in the Bean state in YBa 2Cu 3O 7 x thin films, Europhys. Lett., 65 (6), pp (2004). [23] A. J. Qviller, V. V. Yurchenko, Y. M. Galperin, J. I.Vestgården, P. B. Mozhaev, J. B. Hansen and T.H. Johansen, Quasi-One-Dimensional Intermittent Flux Behavior in Superconducting Films, Phys. Rev. X 2, (2012). [24] M. Baziljevich, E. Baruch-El, T. H. Johansen and Y. Yeshurun, Dendritic instability in YBa 2Cu 3O 7 δ films triggered by transient magnetic fields, Appl. Phys. Lett. 105, (2014). [25] V. K. Vlasko-Vlasov, U. Welp, V. Metlushko and G. W. Crabtree, Experimental test of the self-organized criticality of vortices in superconductors, Phys. Rev. B 69, (R) (2004). [26] Behnia, K., C. Capan, D. Mailly, and B. Etienne, Internal avalanches in a pile of superconducting vortices, Phys. Rev. B 61, 3815, (2000). [27] K. E. Bassler and M. Paczuski, Simple Model of Superconducting Vortex Avalanches, Phys. Rev. Lett. 81, 17, 3761 (1996). [28] M. S. Welling, C. M. Aegerter and R. J. Wijngaarden, Self-organized criticality induced by quenched disorder: Experiments on flux avalanches in NbH x films Phys. Rev. B 71, (2005). [29] G. Pruessner, Self-Organised Criticality: Theory, Models and Characterisation, Cambridge University Press, pp (2012). [30] M. Baziljevich, A. V. Bobyl, D. V. Shantsev, E. Altshuler, T. H. Johansen and S. I. Lee, Origin of dendritic flux patterns in MgB 2 films, Physica C 369, pp (2002). [31] Ch. Jooss, R. Warthmann and H. Kronmüller, Pinning mechanism of vortices at antiphase boundaries in YBaCuO, Phys. Rev. B 61, (2000). [32] A. Polyanskii, R. L. S. Emergo, J. Z. Wu, T. Aytug, D. K. Christen, G. K. Perkins, and D. Larbalestier, Magneto-optical imaging and electromagnetic study of YBa 2Cu 3O 7 vicinal films of variable thickness, Phys. Rev. B 72, (2005). [33] J. H. Kim and D. Youm, Growth of textured c-axis normal Y 1Ba 2Cu 3O 7 δ thin films on yttria-stabilized zirconia substrates with crystalline axes tilted with respect to the surfaces, Physica C 275, 273 (1997). [34] I. K. Bdikin, P. B. Mozhaev, G. A. Ovsyannikov, P. V. Komissinski, and I. M. Kotelyanskii, Growth and domain structure of YBa 2Cu 3O x thin films and YBa 2Cu 3O x/ceo 2 heterostructures on tilted NdGaO 3 substrates, Physica C 377, 26 (2002). [35] P. B. Mozhaev, J. E. Mozhaeva, I. K. Bdikin, I. M. Kotelyanskii, V.A. Lusanov, J. B. Hansen, C. S. Jacobsen and A. L. Kholkin, Out-of-substrate plane orientation control of thin YBa 2Cu 3O x films on NdGaO 3 Tilted-Axes Substrates, Physica C (Amsterdam) 434, 105 (2006). [36] M. Djupmyr, G. Cristiani, H.-U. Habermeier, and J. Albrecht, Anisotropic temperature-dependent current densities in vicinal YBCO, Phys. Rev. B 72, (R) (2005). [37] T. Qureishy, Y. Zhao, Y. Xu, J. I. Vestgården, T. H. Johansen, J-C. Grivel, H. Suo and P. Mikheenko, Magnetic fluxchannelling in YBa 2Cu 3O 7 δ films grown byachemical solution deposition technique on vicinal and nonvicinal substrates, Advanced Materials Letters 8(12), (2017). [38] L. E. Helseth, R. W. Hansen, E. I. Il yashenko, M. Baziljevich and T. H. Johansen, Faraday rotation spectra of bismuth-substituted ferrite garnet films with in-plane magnetization, Phys. Rev B 64, (2001). [39] P. E. Goa, H. Hauglin, Å. A. F. Olsen, M. Baziljevich and T. H. Johansen, Magneto-optical imaging setup for single vortex observation, Rev. Sci. Instrum. 74, 141-

9 9 146 (2003). [40] G. Pruessner, Self-Organised Criticality: Theory, Models and Characterisation, Cambridge University Press, note p. 61 (2012). [41] C. J. Olson, C. Reichhardt, and F. Nori, Superconducting vortex avalanches, voltage bursts, and vortex plastic flow: Effect of the microscopic pinning landscape on the macroscopic properties, Phys. Rev. B 56, 6175 (1997). [42] S. Field, J. Witt, F.Nori and X. Ling, Superconducting Vortex Avalanches, Phys. Rev. Lett (1995). [43] E. Altshuler, T. H. Johansen, Y. Paltiel, P. Jin, K. E. Bassler, O.Ramos, Q. Y. Chen, G. F. Reiter, E. Zeldov and C. W. Chu, Vortex avalanches with robust statistics observed in superconducting niobium, Phys. Rev. B 70, (R) (2004). [44] R. Cruz, R. Mulet and E. Altshuler, Universality of vortex avalanches in a type II superconductor with periodic pinning, Physica A 275, (2000). [45] D. V. Denisov, D. V. Shantsev, Y. M. Galperin, Eun- Mi Choi, Hyun-Sook Lee, Sung-Ik Lee, A. V. Bobyl, P. E. Goa, A. A. F. Olsen, and T. H. Johansen, Onset of dendritic flux avalanches in superconducting films, Phys. Rev. Lett. 97, (2006). [46] V. V. Yurchenko, D. V. Shantsev, T. H. Johansen, M. R. Nevala, I. J. Maasilta, K. Senapati and R. C. Budhani, Reentrant stability of superconducting films and the vanishing of dendritic flux instability, Phys. Rev. B 76, (2007). [47] P. Mikheenko, A. J. Qviller, J. I. Vestgården, S. Chaudhuri, I. J. Maasilta, Y. M. Galperin, and T. H. Johansen, Nanosecond voltage pulses from dendritic flux avalanches in superconducting NbN films, Appl. Phys. Lett. 102, (2013). [48] P. Mikheenko, T. H. Johansen, S. Chaudhuri, I. J. Maasilta and Y. M. Galperin, Ray optics behaviour of flux avalanche propagation in superconducting films, Phys. Rev. B 91, (R) (2015).