Dendritic flux avalanches in superconducting films

Thermomagnetic instability in general, and dendritic flux avalanches in particular, have attracted considerable attention of both scientists and engineers working on superconductor applications. Though being harmful for the performance of many superconducting devices, the avalanches provide a fruitf...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Datum:2009
Hauptverfasser: Yurchenko, V.V., Johansen, T.H., Galperin, Y.M.
Format: Artikel
Sprache:English
Veröffentlicht: Фізико-технічний інститут низьких температур ім. Б.І. Вєркіна НАН України 2009
Schriftenreihe:Физика низких температур
Schlagworte:
Online Zugang:http://dspace.nbuv.gov.ua/handle/123456789/117305
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Назва журналу:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Zitieren:Dendritic flux avalanches in superconducting films / V.V. Yurchenko, T.H. Johansen, Y.M. Galperin // Физика низких температур. — 2009. — Т. 35, № 8-9. — С. 789-797. — Бібліогр.: 36 назв. — англ.

Institution

Digital Library of Periodicals of National Academy of Sciences of Ukraine
id irk-123456789-117305
record_format dspace
spelling irk-123456789-1173052017-05-22T03:02:48Z Dendritic flux avalanches in superconducting films Yurchenko, V.V. Johansen, T.H. Galperin, Y.M. Электронные свойства проводящих систем Thermomagnetic instability in general, and dendritic flux avalanches in particular, have attracted considerable attention of both scientists and engineers working on superconductor applications. Though being harmful for the performance of many superconducting devices, the avalanches provide a fruitful playground for experimental and theoretical studies of complex dynamics of the vortex matter. In this paper, we report on the progress in understanding the mechanisms responsible for the development of the giant magnetic avalanches. We review recent results on magneto-optical imaging of the fingering instability in superconducting films and analyze them basing on the recent theoretical model that establishes criteria for onset of the dendritic avalanches. 2009 Article Dendritic flux avalanches in superconducting films / V.V. Yurchenko, T.H. Johansen, Y.M. Galperin // Физика низких температур. — 2009. — Т. 35, № 8-9. — С. 789-797. — Бібліогр.: 36 назв. — англ. 0132-6414 PACS: 74.25.Qt, 74.25.Ha, 68.60.Dv http://dspace.nbuv.gov.ua/handle/123456789/117305 en Физика низких температур Фізико-технічний інститут низьких температур ім. Б.І. Вєркіна НАН України
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection DSpace DC
language English
topic Электронные свойства проводящих систем
Электронные свойства проводящих систем
spellingShingle Электронные свойства проводящих систем
Электронные свойства проводящих систем
Yurchenko, V.V.
Johansen, T.H.
Galperin, Y.M.
Dendritic flux avalanches in superconducting films
Физика низких температур
description Thermomagnetic instability in general, and dendritic flux avalanches in particular, have attracted considerable attention of both scientists and engineers working on superconductor applications. Though being harmful for the performance of many superconducting devices, the avalanches provide a fruitful playground for experimental and theoretical studies of complex dynamics of the vortex matter. In this paper, we report on the progress in understanding the mechanisms responsible for the development of the giant magnetic avalanches. We review recent results on magneto-optical imaging of the fingering instability in superconducting films and analyze them basing on the recent theoretical model that establishes criteria for onset of the dendritic avalanches.
format Article
author Yurchenko, V.V.
Johansen, T.H.
Galperin, Y.M.
author_facet Yurchenko, V.V.
Johansen, T.H.
Galperin, Y.M.
author_sort Yurchenko, V.V.
title Dendritic flux avalanches in superconducting films
title_short Dendritic flux avalanches in superconducting films
title_full Dendritic flux avalanches in superconducting films
title_fullStr Dendritic flux avalanches in superconducting films
title_full_unstemmed Dendritic flux avalanches in superconducting films
title_sort dendritic flux avalanches in superconducting films
publisher Фізико-технічний інститут низьких температур ім. Б.І. Вєркіна НАН України
publishDate 2009
topic_facet Электронные свойства проводящих систем
url http://dspace.nbuv.gov.ua/handle/123456789/117305
citation_txt Dendritic flux avalanches in superconducting films / V.V. Yurchenko, T.H. Johansen, Y.M. Galperin // Физика низких температур. — 2009. — Т. 35, № 8-9. — С. 789-797. — Бібліогр.: 36 назв. — англ.
series Физика низких температур
work_keys_str_mv AT yurchenkovv dendriticfluxavalanchesinsuperconductingfilms
AT johansenth dendriticfluxavalanchesinsuperconductingfilms
AT galperinym dendriticfluxavalanchesinsuperconductingfilms
first_indexed 2025-07-08T12:00:35Z
last_indexed 2025-07-08T12:00:35Z
_version_ 1837080025028886528
fulltext Fizika Nizkikh Temperatur, 2009, v. 35, Nos. 8/9, p. 789–797 Dendritic flux avalanches in superconducting films V.V. Yurchenko1, T.H. Johansen1,2, and Y.M. Galperin1,3 1 Department of Physics and Centre for Materials Science and Nanotechnology, University of Oslo P.O. Box 1048, Blindern, 0316 Oslo, Norway E-mail: Vitaliy.Yurchenko@smn.uio.no 2 Institute for Superconducting and Electronic Materials, University of Wollongong, NSW, Australia 3 A.F. Ioffe Physico-Technical Institute, St Petersburg 194021, Russia Received February 13, 2009 Thermomagnetic instability in general, and dendritic flux avalanches in particular, have attracted con- siderable attention of both scientists and engineers working on superconductor applications. Though being harmful for the performance of many superconducting devices, the avalanches provide a fruitful playground for experimental and theoretical studies of complex dynamics of the vortex matter. In this paper, we report on the progress in understanding the mechanisms responsible for the development of the giant magnetic ava- lanches. We review recent results on magneto-optical imaging of the fingering instability in superconduct- ing films and analyze them basing on the recent theoretical model that establishes criteria for onset of the dendritic avalanches. PACS: 74.25.Qt Vortex lattices, flux pinning, flux creep; 74.25.Ha Magnetic properties; 68.60.Dv Thermal stability; thermal effects. Keywords: superconducting films, flux avalanche, vortices in superconductor. 1. Introduction 1.1. Vortex pinning and motion One of the characteristic features of superconductors is perfect diamagnetism or expulsion of a weak external magnetic field, H, from their interior, the so-called Meiss- ner and Ochsenfeld effect. In type-II superconductors, the perfect diamagnetism exists for applied fields below a lower critical field, Hc1, and there is a broad domain of magnetic fields, H H Hc c1 2� � , where the field pene- trates into the sample in a form of quantized flux lines. In a perfect sample these lines-vortices form an Abrikosov lattice [1]. Above the upper critical field, Hc2 , the bulk superconductivity seizes to exist. An isolated vortex con- sists of a core where the superconducting order parameter is suppressed, while the magnetic field reaches a local maximum. The radius of the core is of the order of the co- herence length, � . The magnetic field outside the core de- cays exponentially over a distance of a magnetic penetra- tion depth � , where also electrical current circulates around the core. Each vortex carries one flux quantum �0 152 2 07 10� � � �h e/ . Wb. The magnetic field is sup- ported by circular currents around the core. The flux lines repel each other, an interaction that can be understood ap- plying the Ampere’s law to the circular currents. This re- pulsion leads to a formation of a flux line lattice, most typically a hexagonal and in some cases square. A number of phases and dynamic effects in the flux line lattice was reviewed in Refs. 2, 3. The vortices interact with a transport current via the Lorentz force per unit length f j n� �0 [ ] , where j is the current density and n is the unit vector along the flux line. Since the vortex motion implies mov- ing of normal cores it is accompanied with dissipation, which is proportional to the velocity. This dissipation can be described by an effective viscosity, and the velocity is determined by the balance between the Lorentz force and the viscous force. Therefore a free vortex lattice would move as a whole with a constant velocity resulting in a fi- nite resistance of the sample. Fortunately, in all real superconductors the flux lines interact with material defects that act as pinning sites and © V.V. Yurchenko, T.H. Johansen, and Y.M. Galperin, 2009 hamper their motion. Besides the pinning barriers that arise from rather inevitable structural irregularities such as vacancies, dislocations, grain boundaries, etc., there exists a rich zoo of artificially introduced pinning sites like magnetic inclusions, phases of weaker (stronger) su- perconductivity, lithographically patterned «antidots», magnetic dots, etc. According to the particular nature and dimensionality of the defects the pinning potential has different spatial extent, magnetic field and temperature dependencies, see Ref. 4 for a review. When a supercon- ductor is exposed to an increasing magnetic field (or self field of a transport current), vortices are formed at the edges and then propagate inwards. The presence of pin- ning leads to formation of an inhomogeneous distribution of the magnetic flux. According to the critical state model [5] the stationary distribution can be found from Am- pere’s law with the condition that the current density at each point is equal to its local critical value, j Tc ( , )B , i.e., �B j� 0 , | | ( , )j B� j Tc . The case where jc is independent on magnetic field is called the Bean model [5]. The critical current density is in practice a key measure of the applicability of supercon- ductors since jc is the maximum current density in a ma- terial for transport without energy loss. Microscopic eva- luation of the critical current requires direct summation of vortex–vortex interactions and all elementary pinning forces, which all together is a rather unfeasible task. Therefore the critical state model has become a major pa- radigm in the studies of electromagnetic properties of type-II superconductors. The critical state emerging from the balance between the Lorentz and the pinning forces is metastable. Even a slow increase in the external magnetic field may lead to total collapse of the critical state and sudden large-scale redistribution of the flux. Such dramatic events observed experimentally as abrupt drops in the magnetization are ascribed to a thermomagnetic instability where the local heat dissipation associated with vortex motion reduces the pinning, which in turn facilitates further vortex mo- tion. With this positive feedback, a small perturbation can quickly evolve into a flux avalanche [6–19]. 1.2. Experimental methods Experimental techniques used to investigate the ava- lanches in the vortex matter are conventionally subdi- vided in two groups: integral and spatially resolved. The former one includes all types of magnetometry: inductive coils, vibrating sample magnetometry, SQUID mag- netometry [20]. These methods are sensitive to global re- distributions of the flux, current flow, and the resulting change in the total magnetic moment caused by an ava- lanche in the sample. A drawback of the integral methods is a lack of detailed information about the jumps, their lo- cation in the sample, their morphology, etc. Moreover, the relatively low sampling rate in these methods prevents identification of jumps occurring with short time inter- vals. Obviously there is no chance to separate simul- taneous avalanches. Thus, measurement of the volume averaged magnetization excludes a possibility to study detailed statistics of flux jumps. Yet another disadvantage is that it is impossible also to discriminate between small jumps and instrument noise. These problems are partly solved in locally resolved magnetometry, an overview of available methods can be found elsewhere [21]. One of the most powerful tools is magneto-optical imaging (MOI), which combines high magneto-spatial resolution and short acquisition times (restrained only by the exposure time of a camera). This makes MOI the most suitable real-time technique for investigation of flux avalanches. In our experiments we use MOI based on the Faraday effect where the sensing element is an in-plane magne- tized ferrite garnet film with strong Faraday rotation. The film, which is grown separately on a transparent gadolin- ium gallium garnet substrate, is placed directly on the top of the sample under investigation. Stray fields from the magnetic flux pattern of the sample locally remagnetize the indicator, and when shining polarized light through the film, the out-of-plane component of magnetization in- duces Faraday rotation. The effect is doubled after the re- flected light (either from the mirror layer or, in more rare cases, from the sample surface) passes through the garnet film again. Figure 1 illustrates the principal experimental scheme. An analyzer set at 90° relative to the polarizer fil- ters out the nonrotated light. The rotation angle, �, is a monotonous function of the local magnetic field, B. The transmitted light is projected on the matrix of a CCD camera, its intensity, I, being related to the rotation angle according to Malus’ law I I Bo� sin [ ( ) ]2 � � . Here � is an offset angle (equal to 0 when the polarizer and the analyzer are crossed). The dependence �( )B can be determined by a proper calibration of the indicator film. Using this dependence together with Malus’ law one can recalculate the distribution of the magnetic field in the sample from the intensity map of the projected light. Since MOI was invented in the 1950s several materials have been used as indicator films [22]. The MOI results reported in the present paper were obtained utilizing in-plane magnetized garnet films (Lu,Bi)3(Fe,Ga)5O12, grown as a few micron thick epitaxial layer on gadolin- ium gallium garnet (transparent) substrates. One limita- tion of these indicator films is a relatively small magnetic 790 Fizika Nizkikh Temperatur, 2009, v. 35, Nos. 8/9 V.V. Yurchenko, T.H. Johansen, and Y.M. Galperin field range: above � 100 mT the films reach saturation, i.e., they become magnetized totally out-of-plane. 2. Dendritic flux avalanches Thermomagnetic instabilities have been in focus of many experimental and theoretical studies for decades [11,23]. However, only with the use of MOI it has been discovered that in most superconducting films, such as MgB2, Nb, Pb, Nb3Sn, NbN, YBa Cu O2 3 x , and YNi2B2C, the avalanches create dendritic flux patterns (see Fig. 2) [6,11,12,17]. From the massive experimental data collected on the subject one can identify some com- mon features of the dendritic instability: (i) It occurs below a certain temperature Tth where T Tcth � . (ii) It occurs in a limited range of applied fields: H H H 1 2 th th� � , where H 1 th and H 2 th are the so-called lower and upper threshold fields, respectively. Shown in Fig. 3 is a set of images recorded during a slow ramping up of the applied magnetic field in a thin film of NbN after the sample was initially zero-field cooled to T � 4 K. In the image taken at H � 8 Oe, there is only a shallow penetration, with a flux front that gradu- ally advanced during the field increase. This smooth mode of penetration ended abruptly when reaching H � = 12 Oe, where a flux avalanche suddenly occurred. As the applied field increased further, more and more ava- Dendritic flux avalanches in superconducting films Fizika Nizkikh Temperatur, 2009, v. 35, Nos. 8/9 791 CCD camera Analyzer Beam splitter Coil Light source MO film Mirror Sample Cold finger Polarizer Fig. 1. Schematic draft of a MOI setup. A sample is mounted on a massive cold finger of a He-flow cryostat. Resistive coils are used as a source of an external magnetic field. The light from a mercury lamp shines through a polarizer and is guided onto an indicator film, where it experiences the Faraday rota- tion. The light becomes reflected by a thin mirror and pro- jected through an analyzer on a CCD matrix of a computer-op- erated camera. Fig. 2. Dendritic flux avalanche observed in MgB2 film by MOI. From Ref. 6. 8 Oe 42.5 Oe 400 Oe 12 Oe 357 Oe 425 Oe H Th 1 H Th 2 Fig. 3. Magneto-optical images of a NbN film at 4 K in in- creasing applied field illustrating development of the dendritic flux avalanches. From Ref. 24. lanches took place, creating a complex pattern of flux dendrites covering most of the sample area, Fig. 3 (H � = 42.5 Oe). Then, as the field reached H � 357Oe, the av- alanche activity stopped entirely although the field con- tinued to increase. The advancing flux front now erased the previously formed dendritic structures, see images taken at H � 400 and 425 Oe, and the critical state mode of flux penetration was restored. The bright fan-like features in the lower part of the images are due to tiny film defects, and are not related to the flux instability. (iii) The formation of the dendritic instability is a sto- chastic process. Usually indentations on the sample edges serve as the most probable origins of the avalanches, as they give a significant local magnetic field amplification and in- crease in the electric field [26]. Nevertheless, the exact nucleation place of the next dendrite, field interval be- tween two consecutive events, and the final shape of the dendritic structure are nonreproducible. To illustrate this, Fig. 4 shows three magneto-optical images of flux pene- tration in a MgB2 film taken after three repeated expe- rimental runs with initial zero-field-cooling the sample 9.2 K and applying a field of 20 mT. Image (b) was ob- tained by superimposing the three original images co- lored red, green, and blue, respectively. In those parts of the final image where they overlap, i.e., the flux front’s behavior is reproducible, one observes a gradation of white. Different colors show the places where there is no overlap, or the overlap is only partial. Strong irrep- roducibility is seen in the dendrite shapes, while the pene- tration near the edge and along static defects is repeat- able. Randomness in the behavior of the vortex matter also explains a large statistical variation of the measured H 1 2, th values, which will be discussed later. (iv) The degree of branching of the dendrites, some- times represented by their fractal dimensionality, and size vary with temperature and the field. Figure 5 illustrates typical behaviors of the dendrites in a NbN film at different temperatures. At T � 4 K the number of the dendritic avalanches per interval of the field was higher compared to T � 6 K. The size of the den- drites shows opposite trend — it increases when the tem- perature approached Tth . Regarding the dependences on magnetic field, at lower (slowly increasing) fields the avalanches are smaller in size, but appear in larger num- bers per field interval. At larger fields, the dendrites are less frequent, but larger in size. Images in Fig. 3 taken at H �12 and 357 Oe illustrate this property. 3. Theory of the dendritic flux jumps There exist numerous experimental indications that the instability is a rather unavoidable generic feature of the vortex matter in superconductor thin films. The con- ventional theory of the thermomagnetic instability con- sidered only «uniform» flux jumps, where the flux front is smooth and essentially straight [23,27]. The problem of the flux pattern formation was examined in the slab ge- ometry [28] and more recently in thin films placed in per- pendicular magnetic field [29,30]. The criteria for the in- stability onset were obtained from the linear analysis of small coordinate-dependent perturbations in the electric field, E, and temperature, T [25,31]. The problem was solved for a thin (thickness d is much less than the 792 Fizika Nizkikh Temperatur, 2009, v. 35, Nos. 8/9 V.V. Yurchenko, T.H. Johansen, and Y.M. Galperin a b 500 m� Fig. 4. MOI demonstrating irreproducibility of the flux distri- bution that results from the dendritic avalanches. (a) Three im- ages were taken under identical experimental conditions. (b) Image obtained by superimposing the three images above. Repeatable parts of the flux front appear as gradations of white. From Ref. 25. a b Fig. 5. Magneto-optical images of dendritic flux avalanches in a NbN film taken at a) T � 4 K and b) T = 6 K. width 2w) superconducting strip infinite in y direction, placed in a transverse magnetic field, H, and thermally coupled to the substrate (Fig. 6). 3.1. Model and basic equations The thermomagnetic instability can be described by the set including the Maxwell equations and the equation for thermal conduction: curl B j�� 0 , (1) curl E B � � � �t , (2) C T t T � � � �� 2 jE . (3) In the above equations, j is the current density, E is the electric field, B is the magnetic induction, B H�� 0 ; C is the specific heat, � is the thermal conductivity. The heat release into substrate is taken into account by applying the boundary condition ��nT h T T� � �0 0( ) (4) to both film surfaces. Here T0 is temperature of the sub- strate, while h0 characterizes heat transfer through the in- terface between the sample and the substrate. The above set should be supplemented by a relation between the cur- rent density, and the electric field, j j E B T� ( , , ). This rela- tion is chosen as j E� j T g E Ec ( ) ( ) ( / ) . (5) Here g E( ) represents a very steep E j( ) curve having large logarithmic derivative, n E E j j Ec( ) ln / ln /� � � � ��� 1, (6) where � is the differential electrical conductivity, �( ) /E j E� � � . Parameter n stands for the exponent in the frequently used power-law relation E jn� . A strong nonlinearity of the function g E( ) leads to a formation of a quasistatic critical state with j j Tc� ( ), where jc is the critical current density. Note that in this model B-depen- dence of jc is neglected, as in the Bean’s model [5]. We will come back to this assumption later. The key dimensionless parameter of the model is the ratio of the thermal and magnetic diffusion coefficients [23], � � ��� 0 / C . (7) The smaller is �, the slower heat diffuses from the pertur- bation region into the surrounding areas. Hence, one can expect that for smaller � the superconductor is more un- stable, and the formation of instability-induced nonuni- form structures is more likely. Solutions of Eqs. (1)–(5) were sought within linear perturbation analysis in the form T T x y z t � ( , , , ), E E � ( , , , )x y z t , j j � ( , , , )x y z t with perturbations � �T Ex y, , , and �jx y, determined via corresponding dimensionless Fourier amplitudes �, �, and i: � � � � �T T t t ik ikx y� * exp ( / )0 , � � � � �E E t t ik ikx y x y x y, , exp ( / )� 0 , � � � �j j i t t ik ikx y c x y x y, , exp ( / )� 0 . (8) These components can be regarded as z-independent for a thick slab. For the case of a thin film they should be prop- erly averaged over z, see Ref. 31 for details. 3.2. Criteria for dendritic instability in thin films Here we will focus on the analysis of the main results of Ref. 31. Provided there exists a set of parameters that satisfies Re ( , )� k kx y � 0, the superconductor becomes unstable. The wave vectors kx and k y characterize the extent of the flux front perturbations in corresponding directions. Instability would develop to a fingering pat- tern if k kx y�� � 0. The critical k y * and kx * for the finger- ing instability can be found from the requirement max{Re ( )}� k y � 0 for k y � 0 in the limit d l/ 2 1�� (nota- tions follow from Fig. 6): k n nh n x * � �1 � � , k nh n nh n y * / /( ) ( ) � � � � � 1 1 11 4 1 2 , h h C j d j Tc c � � � 2 0 0 2� � ( ln / ) . (9) Dendritic flux avalanches in superconducting films Fizika Nizkikh Temperatur, 2009, v. 35, Nos. 8/9 793 z l no flux j, E d/2 –d/2 H 0 x y flux substrate w m id dl e of th e st rip Fig. 6. One half of a superconductor strip on a substrate. The dark gray area is the flux-penetrated region (Ref. 31). Since the front penetration depth, l, is related to the ap- plied field according to Ref. 32, the quantity k lx * /�1 de- termines the applied magnetic field at which the instabi- lity first takes place. The quantity k y * determines spatial scale of the fingers. Dispersion relations Re ( , )� k kx y for arbitrary k y are presented in Fig. 7. For small � and large kx , i.e., small magnetic field, Re � is negative for all k y (Fig. 7,a), meaning the superconductor is stable. At small kx , Re � becomes positive in some finite range of k y , which implies that perturbations with finger patterns of elevated T and E will start growing perpendicular to the flux front, giving rise to the dendritic instability. For large �, i.e., large heat diffusion, at small kx the maximum of Re � corresponds to k y � 0, and the uniform perturba- tion will be dominant (Fig. 7,b). One of the most important results that can be directly tested in experiments is the threshold flux penetration depth l* , when the superconductor becomes unstable [25]. Putting in Eqs. (9) l k y * */�1 one obtains: l j E h nd j Ec c * | | | | � � � � � � � � � ! ! � " � 2 1 2 0 1 . (10) Here �jc is the temperature derivative of the critical current density. By combining Eq. (10) with the Bean model expression for the flux penetration depth of a thin superconducting strip in transverse field [32] one imme- diately arrives at the expression for the threshold field: H j d w w l c th arccosh� � � � �� � !!" * . (11) Equations (10), (11) set the most practically important relations between the external parameters, such as the field and temperature, on one hand, and the sample prop- erties (�, h0 , jc , w) on the other. In order to check the va- lidity of the model, the threshold field was measured first in a series of MgB2 samples of different width. A set of eight thin film samples was deposited in one run to allow 794 Fizika Nizkikh Temperatur, 2009, v. 35, Nos. 8/9 V.V. Yurchenko, T.H. Johansen, and Y.M. Galperin � = 0.01 � = 7 H 0.50.5 0 –0.5 –1.0 k <x k* x k >x k* x 0 –0.2 –0.4 k* y ky ky kx k* x a b R e � R e � Fig. 7. Dispersion relations Re ( , )� k kx y for small and large �, for d l/ .2 0 001� and n � 20. From Ref. 31. 1.6 mm 0.2 mm Fig. 8. MOI of flux distribution in MgB2 samples of different width at 4 K and magnetic field of 15 mT. From Ref. 25. for a direct comparative analysis and shaped later by photolithography into 3 mm long rectangles having dif- ferent widths ranging from 0.2 to 1.6 mm. In Fig. 8 pre- senting the results of MOI experiments, the dendrites pro- nouncedly invade the wider strips, while there are fewer of them in the samples with smaller widths, and no den- drites at all in the one with w � 02. mm. The quantitative results for H wth ( ) are summarized in Fig. 9. The fitted data reproduce the experimentally observed tendency of the threshold field to diverge as the strip becomes nar- rower. Rephrasing the latter statement, the narrower is the superconductor the more thermomagnetically stable it is. A good agreement between the theory and experiment was also obtained for the threshold field as a function of temperature. To fit the H Tth ( ), observed in Nb and MgB2 thin films (Fig. 10), some assumptions for temperature dependent model parameters had to be made. The as- sumed cubic dependencies for the thermal conductivity, � �� ~( / )T Tc 3 , and the heat transfer coefficient, h0 � � ~ ( / )h T Tc0 3 , were based on the low-temperature data for MgB2 [33] and the acoustic mismatch model [34], respec- tively. The critical current density was assumed to depend linearly on temperature: j j T Tc c c� �0 1( / ). Given a pin- ning potential, U T Tc� �( / )1 , the exponent n U kT~ / be- comes T-dependent, n n� ~ ( / )T Tc �1 . The model clearly reproduces the property (i) mentioned above — existence of a threshold temperature Tth such that for T T� th the in- stability is absent, and a steep increase of the threshold field H th when T approaches Tth . 4. Reentrant stability of superconducting films Although the existence of the lower threshold mag- netic field, H 1 th , was explained in Sec. 3, the origin of the upper threshold field, H 2 th , above which the flux distribu- tions are stable, was not discussed. Recently, we managed to put the existence of an upper threshold into the frame- work of the model described in the previous section [24]. As it follows from Eqs. (10) and (11), the threshold field depends nonmonotonously on the critical current density, as shown in Fig. 11. For intermediate jc the threshold field is nearly constant, and increases slowly as jc becomes larger giving eventually a sublinear asymp- totic dependence, H jcth ~ /3 4 . More importantly, when jc decreases the threshold field diverges at some finite jc , Dendritic flux avalanches in superconducting films Fizika Nizkikh Temperatur, 2009, v. 35, Nos. 8/9 795 0 5 10 15 20 25 w, mm 82 m� 0.2 0.4 0.6 0.8 Dendritic instability No dendrities � 0 th H , m T Fig. 9. Threshold magnetic field for onset of the dendritic in- stability in MgB2 strips of different width. Experimental data (symbols) are in a good agreement with a fitted theoretical curve (full line), for jc � �0 9 1011. A/m 2 . From Ref. 25. 0 2 4 6 8 10 5 10 15 20 25 � 0 th H , m T Nb MgB2 T, K Fig. 10. Temperature dependence of the threshold magnetic field (from Ref. 25). Experimental data obtained for the 5 mm wide MgB2 sample (�) and for a 1.8 mm wide Nb film (�) [13]. The full lines are theoretical fits. The dashed lines show the limiting temperature above which the instability vanishes. j (B)c sm la ler j c(B ) – al a w ys stab el H2 th H1 th Instability region j , arb. unitsc H , ar b . u n it s H (j )th c Fig. 11. Schematic plot illustrating the existence of two thresh- old fields for the dendritic instability adapted from Ref. 24. The two main curves represent dependencies of the threshold filed on the critical current density H jcth ( ) and a typical mo- notonous field dependence of jc. The intersection of the two curves defines the two thresholds H 1 th and H 2 th for the onset and vanishing of the instability. corresponding to l jc * ( ) approaching w. Taking into ac- count that the critical current density most commonly de- creases with the field, it follows that the field range with unstable behavior can indeed have both a lower and upper limit, H 1 th and H 2 th , as indicated in the figure. Evidently, it is essential how fast jc is decreasing with the field. In par- ticular, along the Bean model (constant jc), or with a jc having only weak field dependence, the conditions for having an upper threshold will never be met. Note also that if jc is sufficiently small, e.g., because of a larger temperature, the thermomagnetic avalanches will not occur at any magnetic field. This idea was supported by MOI experiments on a NbN superconductor thin film. To measure both jc and the upper and lower instability threshold fields the sample was field-cooled to 4 K in various constant magnetic fields, H fc . The field-dependent critical current density j Hc ( )fc was determined from the depth of the penetra- tion front in the middle section of the rectangular sample [32,35]. A nearly exponential decay was found, with jc decreasing from 76 1010. � A/m2 at zero field to almost one half at H fc � 300 Oe, which agrees well with the results obtained earlier from ac susceptibility measurements [36] on the same type of films. The threshold fields were deter- mined by slowly ramping up the additional field after an initial field-cooling in extern field H fc . Figures 12 and 13 present the experimentally obtained functions j Hc ( )fc and H H 1 2, ( )th fc , which were used to make the parametric plot of the instability onset field #H � � �H H 1 th fc shown in Fig. 14. Note that under field-cool- ing the vortices are distributed homogeneously inside this superconductor where pinning is fairly strong. Application of additional field leads to a build-up of the critical state, which collapses due to the thermomag- netic instability as soon as the flux front advances over the threshold length l* corresponding to applying the ad- ditional field #H. Figure 14 demonstrates an excellent quantitative agreement between the experimental data for NbN (symbols) and theoretical curve obtained from Eqs. (10) and (11), using only two fitting parameters, �T Eth / � 500 mA and h T nE0 5700th / � A/m (these pa- rameter values can be achieved, e.g., with the combi- nation � � 25 mW/(K·m), h0 1424� W/(K·m2), n � 5, E � 055. V/m, Tth �11 K). As mentioned above, the den- drite formation process is stochastic, which is reflected in a noticeable dispersion of the measured threshold fields and the error bars in the figures, calculated as a standard deviation for repeated experiments. 796 Fizika Nizkikh Temperatur, 2009, v. 35, Nos. 8/9 V.V. Yurchenko, T.H. Johansen, and Y.M. Galperin 7.5 7.0 6.5 6.0 5.0 5.5 4.5 4.0 0 50 100 150 200 250 300 j , 1 0 A /m c 1 0 2 H , Oefc Fig. 12. Field dependence of the critical current density at T � 4 K measured using MOI in NbN thin film (Ref. 24). H2 th H1 th #H Instability 400 300 200 100 0 Hfc Hfc, Oe 0 100 200 300 H , O e th Fig. 13. Thermomagnetic instability diagram of NbN at T � 4 K (Ref. 24): the lower and upper threshold fields H 1 th and H 2 th , measured at different H fc, limit the region where the dendritic avalanches were observed. 4 5 6 7 8 0 10 20 30 40 50 60 70 80 90 # H , O e j , 10 A/mc 10 2 Fig. 14. Instability onset field #H as a function of the critical current at 4 K in NbN thin film. From Ref. 24. 5. Summary and conclusions In this work, we review both well established facts and the results of very recent experiments on dendritic insta- bilities in superconductor films. The most common prop- erties of the dendritic instability are summarized in Sec. 2. In Sec. 3 we introduced the recent theoretical model allowing to understand key features of the instabil- ity and link them to the samples’ parameters. These links are provided by Eqs. (10) and (11) for the threshold field H th as a function of temperature, the critical current den- sity jc , the sample width w, the thermal conductivity �. It was demonstrated that the model explains the observed reentrant behavior of the dendritic instability and allows to quantify its upper and lower critical fields. The validity of the model is illustrated by several experiments confirming its predictions. Acknowledgments This work was supported by the Norwegian Research Council through the grant 176740/I30. One of the authors — YMG — is thankful to Ben-Gurion University of Ne- gev, Israel, and to University of Lancaster, UK, for hospi- tality, and THJ acknowledges support from the Australian Research Council. We are thankful to E.I. Il’yashenko and A. Solovyev for their help in producing magneto-op- tical indicators. 1. A.A. Abrikosov, Sov. Phys. JETP 5, 1174 (1957). 2. G. Blatter, M.V. Feigel’man, V.B. Geshkenbein, A.I. Lar- kin, and V.M. Vinokur, Rev. Mod. Phys. 66, 1125 (1994). 3. E.H. Brandt, Rep. Prog. Phys. 58, 1465 (1995). 4. R. W�rdenweber, Rep. Prog. Phys. 62 187 (1999). 5. C.P. Bean, Rev. Mod. Phys. 36, 31 (1964). 6. T.H. Johansen, M. Baziljevich, D.V. Shantsev, P.E. Goa, Y.M. Galperin, W.N. Kang, H.J. Kim, E.M. Choi, M.-S. Kim, and S.I. Lee, Europhys. Lett. 59, 599 (2002). 7. F.L. Barkov, D.V. Shantsev, T.H. Johansen, P.E. Goa, W.N. Kang, H.J. Kim, E.M. Choi, and S.I. Lee, Phys. Rev. B67, 064513 (2003). 8. Z.X. Ye, Q. Li, Y.F. Hu, A.V. Pogrebnyakov, Y. Cui, X.X. Xi, J.M. Redwing, and Q. Li, Appl. Phys. Lett. 85, 5284 (2004). 9. J. Albrecht, A.T. Matveev, M. Djupmyr, G. Schütz, B. Stuhlhofer, and H.-U. Habermeier, Appl. Phys. Lett. 87, 182501 (2005). 10. F. Laviano, D. Botta, C. Ferdeghini, V. Ferrando, L. Goz- zelino, and E. Mezzetti, in: Magneto-Optical Imaging, T.H. Johansen and D. Shantsev (eds.), Kluwer Academic (2004), p. 237. 11. M.R. Wertheimer and J. de G. Gilchrist, J. Phys. Chem. Solids 28, 2509 (1967). 12. C.A. Duran, P.L. Gammel, R.E. Miller, and D.J. Bishop, Phys. Rev. B52, 75 (1995). 13. M.S. Welling, R.J. Westerwaal, W. Lohstroh, and R.J. Wijngaarden, Physica C411, 11 (2004). 14. M. Menghini, R.J. Wijngaarden, A.V. Silhanek, S. Raedts, and V.V. Moshchalkov, Phys. Rev. B71, 104506 (2005). 15. I.A. Rudnev, S.V. Antonenko, D.V. Shantsev, T.H. Johansen, and A.E. Primenko, Cryogenics 43, 663 (2003). 16. I.A. Rudnev, D.V. Shantsev, T.H. Johansen, and A.E. Pri- menko, Appl. Phys. Lett. 87, 042502 (2005). 17. P. Leiderer, J. Boneberg, P. Brull, V. Bujok, and S. Her- minghaus, Phys. Rev. Lett. 71, 2646 (1993). 18. B. Biehler, B.-U. Runge, P. Leiderer, and R.G. Mints, Phys. Rev. B72, 024532 (2005). 19. S.C. Wimbush, B. Holzapfel, and Ch. Jooss, J. Appl. Phys. 96, 3589 (2004). 20. E. Altshuler and T.H. Johansen, Rev. Mod. Phys. 76, 471 (2004) and references therein. 21. S.J. Bending, Adv. Phys. 48, 449 (1999). 22. M.R. Koblischka and R.W. Wijngaarden, Supercond. Sci. Technol. 8, 199 (1995). 23. R.G. Mints and A.L. Rakhmanov, Rev. Mod. Phys. 53, 551 (1981). 24. V.V. Yurchenko, D.V. Shantsev, T.H. Johansen, M.R. Ne- vala, I.J. Maasilta, K. Senapati, and R.C. Budhani, Phys. Rev. B76, 092504 (2007). 25. D.V. Denisov, D.V. Shantsev, Y.M. Galperin, E.-M. Choi, H.-S. Lee, S.-Ik Lee, A.V. Bobyl, P.E. Goa, A.A.F. Olsen, and T.H. Johansen, Phys. Rev. Lett. 97, 077002 (2006). 26. J.I. Vestgaarden, D.V. Shantsev, Y.M. Galperin, and T.H. Johansen, Phys. Rev. B76, 174509 (2007). 27. A.Vl. Gurevich, R.G. Mints, and A.L. Rakhmanov, The Physics of Composite Superconductors, Begell House Inc., NY (1997). In Russian: A.Vl. Gurevich, R.G. Mints, and A.L. Rakhmanov, The Physics of Composite Superconduc- tors, Science, Moscow (1987). 28. A.L. Rakhmanov, D.V. Shantsev, Y.M. Galperin, and T.H. Johansen, Phys. Rev. B70, 224502 (2004). 29. I.S. Aranson, A. Gurevich, M.S. Welling, R.J. Wijngaar- den, V.K. Vlasko-Vlasov, V.M. Vinokur, and U. Welp, Phys. Rev. Lett. 94, 037002 (2005). 30. C. Baggio, R.E. Goldstein, A.I. Pesci, and W. van Saar- loos, Phys. Rev. B72, 060503 (2005). 31. D.V. Denisov, A.L. Rakhmanov, D.V. Shantsev, Y.M. Gal- perin, and T.H. Johansen, Phys. Rev. B73, 014512 (2006). 32. E.H. Brandt and M. Indenbom, Phys. Rev. B48, 12893 (1993). 33. M. Schneider, D. Lipp, A. Gladun, P. Zahn, A. Handstein, G. Fuchs, S. Drechsler, M. Ritcher, K. Muller, and H. Ros- ner, Physica C363, 3 (2001). 34. E.T. Swartz and R.O. Pohl, Rev. Mod. Phys. 61, 605 (1989). 35. A more precise way to determine j Bc ( ) should be based on the self-consistent relation between the depth of the flux front and the j Bc ( ) law derived in J. McDonald and J.R. Clem, Phys. Rev. B53, 8643 (1996). 36. K. Senapati, N.K. Pandey, R. Nagar, and R.C. Budhani, Phys. Rev. B74, 104514 (2006). Dendritic flux avalanches in superconducting films Fizika Nizkikh Temperatur, 2009, v. 35, Nos. 8/9 797