Physica C 290, 89 (1997)

Plastic flow, voltage noise and vortex avalanches in superconductors

C. J. Olsona, C. Reichhardta, J. Grotha, Stuart B. Fielda,b, Franco Noria,*
aDepartment of Physics, The University of Michigan, Ann Arbor, MI 48109-1120, USA
bDepartment of Physics, Colorado State University, Fort Collins, CO 80523, USA

Received 16 July 1997; accepted 21 July 1997

Abstract

We use large-scale parallel simulations to compute the motion of superconducting magnetic vortices during avalanches triggered by small field increases. We relate observations of pulsing vortex movement in winding chains to features of statistical distributions and experimentally observable voltage noise. As pin density is increased, the very broad distribution of avalanche sizes begins to develop typical avalanche sizes. We show that this corresponds to a crossover from pin-to-pin motion in broad channels to interstitial motion in narrow easy-flow channels. Our results are consistent with recent experiments. 1997 Elsevier Science B.V.

PACS: 64.60.Ht; 74.60.Ge

I. Introduction and summary of results

Recent studies of a variety of physical systems characterized by avalanche dynamics have revealed many interesting dynamical behaviors. Plastic motion, burst-like behavior and broad distributions of event sizes have been observed in systems such as granular assemblies [1-3], magnetic domains [4], fluid flow [5,6], and flux lines in type-II superconductors [7-13]. Although all of the systems display avalanches, changes in the experimental conditions or variations from sample to sample result in the observation of different forms of avalanche size distributions or different spectral signatures in the same system. There is considerable interest in determining what microscopic properties of the different samples are responsible for the macroscopically observed differences. It has, however, not been feasible experimentally to determine the microscopic motion of the individual grains, domain walls, fluid drops or flux lines, while still collecting macroscopic information over time scales long enough to statistically characterize the system, and so no direct experimental evidence of the microscopic nature of the observed change is available. It is therefore natural to turn to numerical simulations in which it is possible to obtain detailed dynamical information while also precisely controlling the microscopic parameters.

Early avalanche simulations used simple discrete models [14] and the small number of simulations employing more realistic continuous molecular dynamics (MD) models were highly restricted in the number of particles used and avalanche events recorded. None of the MD simulations examined the power spectra of the noise signal, an important quantity that is readily accessible in experiments. The simulation presented in this paper overcomes these limitations. With our highly optimized, parallel MD simulation of superconducting vortex avalanches, we studied samples containing a large number of pins (up to 3700) and vortices (up to 1800). We collected detailed information on both vortex positions and voltage noise information for a much greater number of distinct avalanches (~50 000) than has ever been obtained in an MD avalanche simulation. This provides a unique opportunity to directly compare microscopic information with the longer-time statistical information collected in experiments. By varying the pinning density np and maximum pinning strength fp, we obtain avalanche distributions and voltage noise spectra in good agreement with recent experiments [7,12,13], and quantify how the density and strength of pinning sites affects both the breadth of the distributions and the shape of the noise spectra. Characteristic avalanche sizes and lifetimes, along with a change in the form of the voltage spectra, appear at low pin densities when narrow winding channels of interstitial vortices form. At higher pin densities, pin-to-pin transport through vortex chains occurs, and a crossover in the avalanche size and lifetime distributions occurs to a very broad distribution. We find that there is no universal distribution valid for all pinning parameters.

Although we focus specifically on superconducting avalanches in this paper, our system features many aspects that are relevant to other systems that display avalanche dynamics, including the burst-like behavior, broad-band noise and broad distributions of avalanche sizes. Additionally, the channel structures that we observe in this system have analogs in other non-superconducting systems undergoing avalanches, such as when drops of water roll down a slope, or when a localized portion of sand slides down the edge of a sandpile.

2. Simulation

As an external field is slowly increased, a metastable gradient in vortex density, termed the Bean state [15], forms inside a type-II superconductor as vortices are driven into the sample by their mutual repulsion and are held back by defects in the material [7-12]. The competition between pinning and driving forces resembles that found in many other systems that display intermittent or avalanche-like behavior. To model this system, we simulate an infinite slab with a magnetic field H=H applied parallel to the surface so that there are no demagnetization effects. The rigid vortices and straight columnar pins we consider are all parallel to , so we can obtain all relevant dynamical information from a transverse two-dimensional slice, in the xy plane, of the three-dimensional slab. The flux lines evolve according to a T=0 MD algorithm. The vortex-vortex repulsion, given by the modified Bessel function K1(r/), is cut off beyond r=6, where is the penetration depth, so that each vortex interacts with up to 100 neighbors, and important collective effects, neglected in simulations with shorter interaction ranges, are observed. Each 24 x 26 sample contains up to 3700 attractive parabolic pins of radius =0.15 with pinning densities np=0.96/, np=2.40/, or np=5.93/, and pinning strengths uniformly distributed over the range fpmax/5 to fpmax, where fpmax=0.3f0, 1.0f0, or 3.0f0. All forces are given in units of f0=/8 and lengths in units of .

The total force on vortex i is given by
fi=fivv+fivp=vi,
where the force on vortex i from other vortices is
,
and the force from pinning sites is

Here, is the Heaviside step function, ri (vi) is the location (velocity) of vortex i, rk(p) is the location of pinning site k, there are Np pinning sites and Nv vortices, =(ri-rj)/|ri-rj|, =(ri-rk(p))/|ri-rk(p)|, and we take =1.

We have taken advantage of the cut-off on the vortex interaction range to parallelize our code. Using a one-dimensional domain decomposition, we divide the sample into strips that are multiples of the interaction range in width, place each strip on a separate node, and use message passing techniques at the processor boundaries. Load balancing is simplified by the repulsive nature of the interaction which tends to spread the vortices evenly among the processors. With the flexible domain decomposition, the number of processors can be varied without affecting the results. Using roughly 104 hours on an IBM SP parallel computer, we recorded more than 104 avalanches for each of five combinations of np and fp.

A slowly increasing external field is modeled by adding a single vortex to an unpinned region along the sample edge whenever the system reaches mechanical equilibrium [16-19] (we consider dynamical, not thermal, instabilities). This is analogous to adding a single grain of sand to a pile. Most of these small field increases result in only slight shifts in vortex positions, but occasionally one or more vortices will become depinned, producing an avalanche. We find that avalanche disturbances propagate as an uneven pulse, as seen in fig. 4 of Ref. [16]. Events with longer lifetimes often contain more than one pulse of motion (i.e., multiple oscillations in the total avalanche velocity) [20]. By imaging individual avalanches in our samples, we find that a chain of vortices is displaced in a typical event. Each vortex in the chain is depinned, moves a short distance, and comes to rest in a nearby pinning site. Vortices outside the chain transmit stress by shifting very slightly inside pinning sites, but are not depinned. Chain size varies from event to event: in some cases a chain spans the sample, while in other events the chain contains only three or four vortices. In each case, although the disturbance may cross the sample, an individual vortex does not. Thus, the time span of a typical avalanche is much shorter than the time a single vortex takes to traverse the sample.


Fig. 1. Continuous lines indicate the paths vortices follow over an extended period of time covering many avalanches through two 24 x 26 samples with fpmax=3.0f0 and different pinning densities: (a) np=0.96/ and (b) np=5.93/. Open circles mark pinning sites. The presence or absence of easy-flow channels is clearly dependent on pin density. The channels present in (a) lead to avalanches with characteristic sizes and lifetimes superimposed on a broad distribution. Samples with higher pinning density produce very broad distributions of avalanche sizes.

3. Channels

When the pinning density is high, chains of moving vortices are equally likely to form anywhere in the sample. As the pinning density is lowered, vortices move only through well-defined winding interstitial channels in which vortices are weakly held in place only by the repulsion of other vortices which sit in pinning sites. To identify the cumulative pattern of flow channels for different pinning parameters, in Fig. 1 we plot vortex trajectories with lines over the course of many avalanches. A concentration of trajectory lines indicates a heavily-traveled region of the sample. In Fig. 1(a) we show a sample with low pinning density. We see that all motion occurs through narrow easy-flow interstitial channels. In these channels, mobile interstitial vortices move plastically [21-25] (note, in relation to Ref. [10,11]. As the pin density increases, the amount of interstitial pinning decreases and the number of channels increases, until at high pin densities (Fig. 1b), there is no interstitial flow and the vortices move only from pin to pin. Here, where no well-defined easy-flow channel exists, avalanches are spread evenly throughout the sample.


Fig. 2. Avalanche distributions: (a) lifetimes ; (b) individual vortex displacements di. Inset of (a): number Nf of vortices falling off the edge of the sample. Solid symbols refer to samples with high pin density np=5.93/, and differing pinning strengths: plus signs (heavy solid line), fpmax=3.0f0; filled diamonds (solid line), fpmax=1.0f0; filled circles (dot-dashed line), fpmax=0.3f0. Open symbols refer to samples with fpmax=3.0f0 and varying pinning densities: plus signs (heavy solid line), np=5.93/; open squares (dashed line), np=2.40/; open triangles (heavy dotted line), np=0.96/.

4. Avalanche lifetimes

We determine how the microscopic pinning parameters affect avalanche size by finding the total avalanche lifetime for each event. Since vortices typically hop from pin to pin in samples with high pinning density, a natural unit of time is the interval th a vortex spends hopping between pinning sites, and so we use scaled lifetimes =/th. To find th, we assume that each vortex hops a distance dp=np-1/2, the average distance between pinning sites, and that the vortex speed vc is proportional to the depinning force,
,
where |fi|-fp. This gives

The plot of P() in Fig. 2a shows that the distribution for samples with dense pinning is very broad and can be written as
P()~()-1.4
over a range of . The form of P() changes noticeably when the pinning density np is reduced, and an enhanced probability for avalanches with a characteristic value (arrow in Fig. 2a) arises as a result of the appearance of easy-flow interstitial channels. Many avalanches in these samples consist of a single sample-spanning pulse of motion through one of these channels, seen in Fig. 1a. The estimated pulse lifetime produced by a straight channel is

where Nh is the number of vortices in the channel. Since

where Lx is the sample length, we find
,
which gives a characteristic value of

This value, agrees well with the peak in the distribution of indicated by an arrow in Fig. 2a.

5. Displacements

For all pin densities, only a small fraction of the vortices move significantly while the rest slightly shift in pinning sites, as seen by considering the distance di each vortex is displaced in an avalanche. Vortices that hop from pin to pin create a peak in P(di), marked with arrows in Fig. 2b, at
didp=np-1/2.
For those vortices that remain pinned and accumulate stress in the vortex lattice, di< dp, we can approximate the distribution by

where ~1.4 for all samples except ~1.2 for [np=0.96/, fpmax=3.0f0] and ~0.9 for [np=5.93/, fpmax=0.3f0]. An analytical argument [16], sketched here, predicts a similar for all samples since the distribution is generated only by pinned vortices and is not affected by the presence of easy-flow channels. The addition of a vortex to the sample exerts a small additional force f on an arbitrary vortex located a small distance r away (r), displacing this vortex a distance
u(r)=(f/)t~K1(r/)~1/r.
Since there are
N(r)=2r nvr
vortices a distance r from the added vortex, we find
=-d lnN/d lnu=-(1/r)(-1/r)-1=1,
in general agreement with our computed values of ,
0.9 to 1.4.

6. Number of vortices leaving the sample

Altering the pinning parameters affects the number of vortices Nf that exit the sample during an avalanche, as shown in the plot of P(Nf) in the inset of Fig. 2a. If we approximate P(Nf) for small Nf by the form

we find that all samples with high pinning density, np=5.93/, have ~2.4. As np decreases, increases: ~3.4 for np=2.40/ and ~4.4 for np=0.96/. When all vortex motion occurs in an easy-flow interstitial channel, increases since the channel does not build up enough stress to allow events with large Nf to occur. Smaller events continually relieve the accumulated stress instead. For example, with the small number of flux paths samples with np=0.96/, as in Fig. 1a, events with large Nf are rare, and ~3.4. In samples with high pin density, even after the stress in one vortex path has been depleted by a large avalanche, other areas still contain enough stress to remain active in large and small events while the depleted regions build up stress again. This leads to a greater likelihood of large events and correspondingly smaller values of , 2.4.

7. Relation to experiment

Distributions similar to the high pin density case considered here have been obtained experimentally in [7], where values of ranging from

are observed. The pinning density from grain boundaries in the experimental sample is very high, np~100/, so it is reasonable that the values observed in [7] are similar to the values produced by our most densely pinned samples. Broad distributions with

were also observed in [7], in good agreement with both [7] and our results. In addition, [12] finds a regime where avalanches of a characteristic size occur, offering an experimental example in which samples with a lower density of weaker pinning sites produce narrower distributions, as we also observe.

8. Voltage noise power

Voltage noise power can be used to probe the interesting peak effect that appears as an increase in the critical current with increasing applied magnetic field. The resulting increase in noise power has been postulated to result from plastic motion of the vortex lattice [13]. We compare our simulation to experiment by using the analog of an experimentally measured voltage signal, the average vortex velocity

Our raw signal vav, seen in fig. 2 of Ref. [16], strongly resembles experimentally observed voltages, as in fig. 1 of Ref. [7]. To find the power spectra, we measure time in units of the average time of flight tf of a single vortex across the sample. Such units make it possible to compare our work to experiments in which the time of flight is known. For example, in Ref. [26] vortices travelled 0.02 mm in about 7.2 ms. Taking for YBCO to be =140 nm, the vortices in our simulation travel the equivalent of 0.0036 mm, and if they are considered to move at the same speed as the vortices of Ref. [26], the time of flight becomes tf1.3 ms. In these units, the spectra shown here range from 300 Hz to 3 MHz in frequency, so the lower part of the frequency range may be compared with experiments.


Fig. 3. Voltage noise spectrum from each type of sample, with frequencies measured in units of the inverse of the time of flight 1/tf. Heavy solid line, fpmax=3.0f0, np=5.93/; solid line, fpmax=1.0f0, np=5.93/; dashed line, fpmax=3.0f0, np=2.40/; dot-dashed line, fpmax=0.3f0, np=5.93/; heavy dotted line, fpmax=3.0f0, np=0.96/. Each power spectrum was obtained from time series totaling 106 MD steps per sample. Inset: the spectrum of the voltage signal for a sample with high pinning strength fpmax=3.0f0 and low pinning density np=0.96/ has more than one slope.

From the power spectra of vav for each of our samples, shown in Fig. 3, it is clear that our simulated plastically moving vortex lattice produces broad-band noise that changes in form in a manner consistent with experiment [7,13]. The overall noise power in the simulations is reduced both as the pinning density is reduced and as the pinning strength is lowered. For frequencies 10-3, the spectrum is flat, indicating that avalanches separated by times 103 MD steps are uncorrelated. For 10-3, samples with a high density of pinning sites produce spectra of the form

where increases as the pinning strength decreases: =1.54 for fpmax=3.0f0, =1.66 for fpmax=1.0f0, and =1.93 for fpmax=0.3f0. Since samples with stronger pinning effectively have a lower driving rate than samples with weaker pinning, the change in slope appears consistent with Ref. [7], where the slope increased from ~1.5 to 2.0 as the driving rate increased. The exponents observed here are also similar to those found by Marley et al. [13], who obtained exponents ~1.5 to 2.0 in the peak effect region near depinning where plastic flow should be occurring.

The linear form of the spectra breaks down in samples with a lower density of pinning sites, which have more correlated vortex motion due to the presence of interstitial channels with typical avalanche lifetimes (102103). This tends to produce a steeper slope in the power spectrum at low frequencies (10-310-2). This is shown in the inset of Fig. 3, where the spectra of these samples have a region of relatively steep slope at lower frequencies and a region with a more gentle slope at higher frequencies.

9. Conclusions

We have quantitatively shown how pinning determines the nature of vortex avalanches. By using large-scale MD simulations, we image pulses of motion in chain-like disturbances through the sample. The presence or absence of distinct channels for flow leads to a crossover from broad distributions of avalanche size to characteristic sizes. Pinning strength causes a transition between strongly plastic flow to mildly plastic, 'semi-elastic' flow. Lowering the pinning density causes the very broad avalanche distributions to develop characteristic sizes, as well as changing the nature of the voltage noise spectrum.

Acknowledgments

Computer services were provided by: the Maui High Performance Computing Center, sponsored in part by the Phillips Laboratory, Air Force Materiel Command, USAF, under cooperative agreement number F29601-93-2-0001; and by the University of Michigan Center for Parallel Computing, partially funded by NSF grant CDA-92-14296. CO was supported by the NASA Graduate Student Researchers Program.

References

  1. G.A. Held et al., Phys. Rev. Lett. 65 (1990) 1120.
  2. M. Bretz et al., Phys. Rev. Lett. 69 (1992) 2431.
  3. H.M. Jaeger, S.R. Nagel, R.P. Behringer, Rev. Mod. Phys. 68 (1996) 1259, and references therein.
  4. K.L. Babcock, R.M. Westervelt, Phys. Rev. Lett. 64 (1990) 2168.
  5. B. Plourde, F. Nori, M. Bretz, Phys. Rev. Lett. 71 (1993) 2749.
  6. O. Narayan, D.S. Fisher, Phys. Rev. B 49 (1994) 9469.
  7. S. Field et al., Phys. Rev. Lett. 74 (1995) 1206.
  8. R.J. Zieve et al., Phys. Rev. B 53 (11)(1996) 849.
  9. O. Pla, F. Nori, Phys. Rev. Lett. 67 (1991) 919.
  10. T. Matsuda et al., Science 271 (1996) 1393.
  11. F. Nori, Science 271 (1996) 1373.
  12. E.R. Nowak et al., Phys. Rev. B 55 (11)(1997) 702.
  13. A.C. Marley et al., Phys. Rev. Lett. 74 (1995) 3029.
  14. P. Bak, C. Tang, K. Wiesenfeld, Phys. Rev. A 38 (1988) 364.
  15. C.P. Bean, Rev. Mod. Phys. 36 (1964) 31.
  16. C.J. Olson, et al., Phys. Rev. B 56, in press, available in Ref. [20].
  17. C. Reichhardt et al., Phys. Rev. B 52 (10) (1995) 441.
  18. C. Reichhardt et al., Phys. Rev. B 53 (1996) R8898.
  19. C. Reichhardt et al., Phys. Rev. B 54 (16) (1996) 108.
  20. Brief vortex avalanche videos are available at http://www-personal.engin.umich.edu/~nori.
  21. R.D. Merithew et al., Phys. Rev. Lett. 77 (1996) 3197.
  22. M.C. Faleski, M.C. Marchetti, A.A. Middleton, Phys. Rev. B 54 (12) (1996) 427.
  23. H.J. Jensen et al., Phys. Rev. B 38 (1988) 9235.
  24. A.E. Koshelev, V.M. Vinokur, Phys. Rev. Lett. 73 (1994) 3580.
  25. N. Gronbech-Jensen, A.R. Bishop, D. Dominguez, Phys. Rev. Lett. 76 (1996) 2985.
  26. G. D'Anna et al., Phys. Rev. Lett. 75 (1995) 3521.
Physica C 290, 89 (1997).

Back to Home CJOR
Back to Home CR
Last Modified: 1/1/02