# Simulation of the vortex dynamics in a real pinning landscape

of YBaCuO coated conductors

###### Abstract

The ability of high-temperature superconductors (HTSs) to carry very large currents with almost no dissipation makes them irreplaceable for high-power applications. The development and further improvement of HTS-based cables requires an in-depth understanding of the superconducting vortex dynamics in presence of complex pinning landscapes. We present a critical current analysis of a real HTS sample in a magnetic field by combining state-of-the-art large-scale Ginzburg-Landau simulations with reconstructive three-dimensional scanning transmission electron microscopy tomography of the pinning landscape in Dy-doped YBaCuO. This methodology provides a unique look at the vortex dynamics in the presence of a complex pinning landscape, responsible for the high current-carrying capacity characteristic of commercial HTS wires. Our method demonstrates very good functional and quantitative agreement of the critical current between simulation and experiment, providing a new predictive tool for HTS wires design.

###### pacs:

74.20.De, 74.25.Sv, 74.25.Wx, 05.10.a## I Introduction

Commercial high-temperature superconducting (HTS) wires are being successfully applied in a variety of electric power equipment for the power grid. Their high current-carrying capacity and low dissipation provide multiple advantages over conventional conductors.Malozemoff (2012); Shiohara et al. (2012) However, even commercial second-generation HTS wires have large headroom for improvement.Senatore et al. (2014); Obradors and Puig (2014) In particular, the critical current in these YBaCuO (YBCO) HTS wires decreases rapidly in magnetic fields, prohibiting them from use in broader applications such as in superconducting electrical power generators that could enable lightweight and efficient compact systems, e.g., for wind turbines. In this work we present the first vortex pinning simulation in a reconstructed mixed-pinning landscape obtained by three-dimensional (3D) STEM tomography of an actual section of a Dy-doped YBCO wire, enabling a new strategy for optimizing the critical current in HTS wires.

Energy dissipation in superconductors in the presence of an applied field arises from the motion of vortices driven by the current-induced Lorentz force,Bardeen and Stephen (1965) thus restricting their mobility through pinning by admixed inclusions is the main route to minimize dissipation and increase the critical current.Foltyn et al. (2007); Holesinger et al. (2008); Maiorov et al. (2009); Jia et al. (2013) At present, the quest for higher critical current in HTS is carried out mostly via the laborious process of empirical trial and error. (Only few systematic studies of the critical current dependence on sizes and densities of defects have been published.Miura et al. (2013a); Selvamanickam et al. (2015); Sadovskyy et al. (2015a)) However, at the fundamental level, the basic principles of vortex pinning have been established, at least for simple idealized situations, see, e.g., Refs. Blatter et al., 1994; Brandt, 1995; Blatter and Geshkenbein, 2003. A major impediment to rapid progress in improving the performance of superconducting wires for applications is an insufficient understanding of vortex dynamics in the complex pinning landscape of real materials.

Vortex pinning is a complicated collective phenomenon controlled by the interaction of vortices with pinning centers as well as the flexibility of vortex lines and inter-vortex interactions. The analytical treatment of this problem has been limited to qualitative estimates of the critical current in simple defect environments such as weak pinning by large densities of atomic impuritiesLarkin and Ovchinnikov (1979) or strong pinning by low densities of strong inclusions.Ovchinnikov and Ivlev (1991); Blatter et al. (2004) On the other hand, numerical simulations of vortex dynamics have been used to provide a better insight into the pinning mechanisms and to improve the quantitative characterization. Here, the choice among several models is determined by a trade-off between complexity and fidelity. For example, the minimal approach is to consider only the vortex degrees of freedom and treat vortices as particles (in thin films) or elastic strings (in the bulk). In this case, the dynamics is described by an overdamped equation of motion, which takes into account thermal Langevin forces, see, e.g., Refs. Ertaş and Kardar, 1996; Bustingorry et al., 2007; Luo and Hu, 2007; Koshelev and Kolton, 2011. This approach provides a reasonable description of vortex dynamics at small magnetic fields and small densities of pinning sites. However, vortex-vortex and vortex-defect interactions are treated only approximately in the Langevin dynamics approach and vortex cuttings and reconnections cannot be described within this model.

In contrast, the time-dependent Ginzburg-Landau (TDGL) approachSchmid (1966); Aranson and Kramer (2002) provides an ideal compromise between approximate phenomenological and exact microscopic descriptions of vortex matter. The TDGL equations automatically account for the flexibility of the vortex lines, the long-range mutual vortex repulsion, vortex cutting and reconnection. Pinning defects of arbitrary shape and size can be incorporated via spatial modulation of the critical temperature. Furthermore, the interruption of current paths by replacement of a current-carrying superconductor with non-superconducting pinning defects is automatically accounted for by the model. TDGL-based numerical simulations have been used several times in the past to study various properties of vortex matter.Doria et al. (1990); Machida and Kaburaki (1993); Crabtree et al. (1996); Aranson et al. (1996); Crabtree et al. (2000); Winiecki and Adams (2002); Vodolazov (2013); Berdiyorov et al. (2014)

Recent advancements in computational capabilities in combination with efficient parallel solvers for the TDGL equation on modern graphics processing units (GPUs)Sadovskyy et al. (2015b) enabled the simulation of rather large 3D samples with different pinning landscapes making meaningful predictions for the behavior of critical currents possible. Practically all theoretical and numerical studies of vortex pinning dealt with idealized models where only one type of pinning center is typically considered. However, commercial HTS wires have been engineered with a complex variety of pinning defects of different sizes and shapes, which have empirically been found beneficial for high critical current densities. A direct simulation of the vortex behavior in such mixed landscape is usually not possible because information about internal structures remains qualitative and scarce. A recent scanning transmission electron microscopy (STEM) tomography study,Ortalan et al. (2009) has determined the location and size of inclusions within a superconducting Dy-doped YBCO compound used for second-generation HTS wires, see Fig. LABEL:sub@fig:setup. This technique allows for unprecedented exact mapping of inclusion sizes and locations within the sample.

Here, we combine the advances in large-scale TDGL simulations and STEM tomography information to produce the first comparison of the computed critical currents for a realistic sample with experimental measurements at different magnitudes and orientations of the external magnetic field. We simulate the complete volume and nanoparticles pinning structure of the experimental system down to the resolution limit of the tomography.

## Ii Model

We use TDGL equations as the main tool for the numerical analysis of vortex dynamics. In the infinite- limit, these equations provide a quantitatively adequate description of strong type-II superconductors at high magnetic fields. In this limit, the Maxwell-type equation for the vector potential is eliminated and we concentrate on the remaining equation for the superconducting order parameter ,

(1) |

where is the chemical potential, is the vector potential associated with the external magnetic field as , and is the temperature-dependent -correlated Langevin term.Sadovskyy et al. (2015b) In Eq. (1), written in the dimensionless form, the unit of length is given by the superconducting coherence length , the unit of time is , where is the London penetration depth, the normal state conductance, and the unit of magnetic field is given by the -axis upper critical field . Here is the electron’s charge and the speed of light. In our notations of the TDGL equations, the -plane of the HTS is in the -plane of the simulation and the -axis is along the -direction.

To account for the anisotropy and layered structure of the material, we replace the -component of the Laplacian term in Eq. (1), , with the discrete second derivative , where is the anisotropy factor and is the grid point spacing in -direction, . While similar numerical discretizations are used for the in-plane directions ( and ), the discreteness in direction is an essential part of the model. This description is known as the Lawrence-Doniach model.Lawrence and Doniach (1970) In this model the effective interlayer coupling is controlled by the parameter . In particular, the layered structure leads to intrinsic pinning of vortices for magnetic fields applied in-plane.

The dimensionless function vanishes at the local critical temperature . The explicit dependence of the critical temperature on spatial coordinate is used to model large-scale inhomogeneities in the superconductor and is a convenient way to introduce pinning effects. We use inside the superconductor. This corresponds to setting the coherence length at a given temperature as the unit of length. Inside non-superconducting inclusions we choose .

The total (normal and superconducting) in-plane reduced current density in units of is given by the expression

(2) |

The maximum theoretical depairing current density is then . In addition to Eq. (1), we solve the Poisson equation, , for the scalar potential .

To determine the critical current value, , we apply an external current, , in the -plane (in the -direction) and ramp it down from the resistive to the superconducting state of the sample. Each time step we calculate the electric field in direction of the applied current averaged across the sample, i.e., . The - curve is then obtained by averaging the electric field over the steady states reached after the transient regime following each current ramping event. We use a finite electric field criterion to determine the critical current, where is the electric field unit. In other words, we compare electrical field induced by the external current with a certain electric field level . The defined level is sufficiently low for a reasonable definition of the critical current. However, it is much higher than the level of dissipation corresponding to the value of V/cm, routinely used as a practical criterion for . Therefore, the simulated critical currents are expected to be somewhat higher than the experimental ones.

To simulate the anisotropy of the critical current we apply and rotate the magnetic field from the -plane to the -axis of the HTS keeping it perpendicular to the external current. In the coordinate system of the simulation, the applied field can be described by the gauge , where is the absolute value of the field and the angle with respect to the -axis.

In the simulated volume, pinning centers are positioned in strict accordance with their positions and sizes as measured by 3D STEM tomographyOrtalan et al. (2009) in the experimental sample. The original reconstruction of the superconductor volume together with the particles obtained by the segmentation of the STEM tomogram using the IMOD software is shown in Fig. LABEL:sub@fig:setup. The defect inclusions are confined within a rectangular box of size nm. This box contains 71 almost spherical particles with sizes ranging from 12.2 to 100 nm. For the numerical analysis we use a coherence length nm in the -plane, which is close to the experimental value at 77 K and an anisotropy factor of suitable for YBCO. We simulate a volume of size , which corresponds to nm with our choice of . The non-superconducting particles with diameters ranging from 2.90 to occupy about 8.1% of the simulated volume. A typical vortex configuration for magnetic field T is shown in Fig. LABEL:sub@fig:vortices, where T is the -axis upper critical field at 77 K. Isosurfaces of the order parameter, shown in red in Fig. LABEL:sub@fig:vortices, reveal the vortex positions and the contours of the pinning landscape. On the back and sides of the simulated volume, a color-code indicates the amplitude of .

## Iii Results

In Fig. 2, we present dependence of the critical current, , on the magnetic field applied parallel to the -axis of the HTS (). Experimentally,Ortalan et al. (2009) at K and in a magnetic field ranging from 0.02 to T the critical current follows a power-law dependence with an exponent decreasing with increasing Dy doping (green curve) and (yellow curve) for Dy doping and respectively.^{1}^{1}1Note, that we refitted the experimental data in the same field range as for the simulation data. Therefore the exponents are slightly different from the ones reported in Ref. Long et al., 2005. At lower magnetic field, T, the critical current is more or less independent of the external field, since self-field effects dominate. In this work we concentrate on higher-field region, T, where this effect is negligible. The red curve in Fig. 2 presents the simulated dependence from the 0.5 Dy doped sample shown in Fig. LABEL:sub@fig:setup, using only the STEM-resolved nanoparticles. For ranging from (0.1 T) to (1.5 T), corresponding to the range used in the experimental measurements, the simulated is also described by a power law but with the slightly larger exponent, . To show the effects of defects smaller than the resolution of the tomography, we also added 1500 spherical particles of size (8.4 nm) to the reconstructed defects. These background inclusions reduce the exponent of the dependence to , see blue curve in Fig. 2.

In Fig. 3, we present the anisotropy of the critical current in a tilted magnetic field, always perpendicular to the applied current. The dependencies on the tilt angle of the magnetic field with respect to the -axis for three different magnetic fields are shown. The red curves correspond to the simulation with only the defects detected by STEM tomography, while the blue lines correspond to the simulation with the additional background inclusions. The anisotropy of increases with magnetic field and at the highest fields, , is nearly flat at all angles except for a narrow peak around caused by the simulated intrinsic pinning due to the layered structure of the material. The ratio varies depending on the field: for [Fig. LABEL:sub@fig:Jc_angle_B=0p005], for [Fig. LABEL:sub@fig:Jc_angle_B=0p05], and for [Fig. LABEL:sub@fig:Jc_angle_B=0p1]. This ratio is mainly determined by the Lawrence-Doniach parameter .

## Iv Discussions

The almost quantitative agreement between the exponents in the simulation (0.80) and the measurements (0.74) is quite striking for this first of its kind numerical simulation of a real 3D sample. The small discrepancy in the exponents is most likely related to the presence of background inclusions that are smaller than the resolution of the STEM tomogram. In fact, we have shown that those additional small defects reduce the field exponent in the simulations. Such small inclusions, that several groups are trying to control through chemical processes,Gutierrez et al. (2007); Miura et al. (2013a, b) become especially relevant at high fields when all strong pinning sites are already occupied. In this regime it was also found in experiments that their effect is to reduce the exponent of the field dependence of the critical current.Jia et al. (2013); Leroux et al. (2015) Specifically, the addition of spherical nanoparticles with diameter , occupying only 1.2% of the sample volume, reduces the exponent from (red curve in Fig. 2) to (blue curve). The latter essentially coincides with the experimental value for 0.75 Dy doping. The increased Dy content is expected to create more nanoparticles and thus to reduce the exponent . Our simulation therefore appears to correctly reproduce this trend.

Next, we compare the absolute values of the critical current. In the experiment,Ortalan et al. (2009) the critical current density varies from MA/cm at magnetic field T to MA/cm at T, for applied along the -axis. From our numerical simulations for these two fields we obtain MA/cm and MA/cm, respectively. Here we estimate MA/cm using nm and nm at 77 K. The absolute values are also in a reasonable agreement (within a factor of two) with the experiment. This semi-quantitative agreement is quite remarkable considering that we included only the Dy nanoparticles as pinning centers but neglected other types of defects present in YBCO films such as twin boundaries, atomic point defects, stacking faults, and dislocations. In addition, at least some of the larger simulated values of may be attributed to the higher critical electric field criterion used for their determination, as explained in the model section. This suggests that the nanoparticle defects play the dominant role in the field dependence of the critical current in the measured commercial tape.

As shown in Fig. LABEL:sub@fig:Jc_angle_B=0p05, the overall shape of the angular dependence appears consistent with the experiment,Ortalan et al. (2009) in particular the flat plateau around () is clearly reproduced. The overall shape of angular dependences is expected. As the nanoparticles have an isotropic shape and no -axis correlated pinning centers (e.g. dislocations or twin boundaries) were introduced, there is no sharp peak in for . The only correlated pinning is the intrinsic pinning in the -plane which does yield a peak in . In addition, in an anisotropic material is expected to have smooth angular dependence with the typical angle , as we observe at small field, Fig. LABEL:sub@fig:Jc_angle_B=0p005. With increasing magnetic field, the simulated anisotropy is found to increase, a trend also in agreement with experiments.Holesinger et al. (2008) In absolute terms, the simulated anisotropy is for T, whereas the experimental anisotropy is somewhat smaller for 0.5 Dy doping and for the undoped sample. This discrepancy is easily explained since our choice of the Lawrence-Doniach parameter yields qualitatively reasonable anisotropy but would need to be fine-tuned against experimental measurements of the intrinsic pinning alone in clean single crystal for really quantitative agreement. More importantly, beyond intrinsic pinning we did not take into account other types of -plane pinning centers, such as flat precipitates, as they do not appear in the tomogram. Nonetheless, our model of -plane pinning is reasonable on a qualitative level, and the simulated angular dependence does exhibit the proper trends for the combined effects of nanoparticles and intrinsic pinning.

Finally, from a practical point of view, the most relevant question is probably how “optimal” the observed pinning configuration is in terms of absolute critical current values. The question can be answered in part by looking at a recent study,Koshelev et al. (2015) by some of the present authors, on monodisperse spherical defects using the same approach. It was found that the optimal critical current for such defects, in magnetic fields ranging from to , was achieved for defect diameters ranging between to and occupying a volume fraction of about . In the sample measured by tomography, the nanoparticles occupy a total volume fraction of %. The distribution of diameters of these particles peaks in the range to (72% of all particle are in this range) [see histogram Fig. 5(d) in Ref. Ortalan et al., 2009]. However, the volume fraction occupied by the above mentioned “optimal” particles is only . Even counting all particles from to and taking into account the added background inclusions, barely increases the volume fraction of these pinning centers to 2%, which is ten times less than the optimal density found in Ref. Koshelev et al., 2015. Most of the defect volume fraction is actually due to a few large nanoparticles, which may capture several vortices at a time [see Fig. LABEL:sub@fig:vortices]. According to our simulations, these large defects contribute little to pinning and also reduce the effective cross-section of the sample, thus negatively affecting the critical current.

In fact, the critical currents found in this study are about three times smaller than for the optimal pinning configuration for monodisperse spherical defects mentioned above. This therefore suggests that the present configuration is far from optimal and could be improved by skewing the distribution of nanoparticles toward smaller sizes, for instance 10–20 nm in diameter, and simultaneously by raising their density. This is in line with recent experimental results that find extremely effective pinning with nanoparticles of size 15–30 nmGutierrez et al. (2007) or 20–80 nmMiura et al. (2013a) in YBCO, and 8 nm in pnictides.Miura et al. (2013b) Also, similar results were observed with smaller irradiation-induced defects, 5 nm in size.Jia et al. (2013); Leroux et al. (2015)

## V Conclusions

In summary, we numerically simulated the superconducting vortex dynamics in a real pinning landscape of nanoparticles inside a commercial Dy-doped YBCO tape, using the same sample size as in the experiment. The positions and sizes of these nanoparticles were obtained directly through STEM tomography reconstruction of the sample nanostructure. We obtained good qualitative and almost quantitative agreement in the functional dependencies and absolute values of the critical current between simulation and experiment. One can expect that more detailed STEM tomography studies, and in particular more quantitative characterizations of pinning centers, should improve the quantitative agreement. In addition, there is still room for improving the simulation in terms of exact treatment of intrinsic pinning as well as lowering the voltage criterion towards the one used in the industry. Our results, however, show a promising pathway for the quantitative analysis and optimization of vortex dynamics in various realistic defects environment and physical conditions, using a combination of 3D tomography analysis and large-scale time-dependent Ginzburg-Landau simulations.

### Acknowledgements

We are delighted to thank A. P. Malozemoff, W.-K. Kwok, and U. Welp for careful reading of the manuscript and numerous useful comments. This work was supported by the Scientific Discovery through Advanced Computing (SciDAC) program funded by U.S. Department of Energy Office of Science, Advanced Scientific Computing Research. Simulations were performed at Argonne LCF supported by DOE under contract DE-AC02-06CH11357. A. E. K. and M. L. acknowledge support by the Center for Emergent Superconductivity, an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences Award No. DEAC0298CH1088.

## References

- Malozemoff (2012) A. Malozemoff, Second-generation high-temperature superconductor wires for the electric power grid, Annu. Rev. Mater. Res. 42, 373 (2012).
- Shiohara et al. (2012) Y. Shiohara, T. Taneda, and M. Yoshizumi, Overview of materials and power applications of coated conductors project, Jpn. J. Appl. Phys. 51, 010007 (2012).
- Senatore et al. (2014) C. Senatore, M. Alessandrini, A. Lucarelli, R. Tediosi, D. Uglietti, and Y. Iwasa, Progresses and challenges in the development of high-field solenoidal magnets based on re123 coated conductors, Supercond. Sci. Technol. 27, 103001 (2014).
- Obradors and Puig (2014) X. Obradors and T. Puig, Coated conductors for power applications: materials challenges, Supercond. Sci. Tech. 27, 044003 (2014).
- Bardeen and Stephen (1965) J. Bardeen and M. J. Stephen, Theory of the motion of vortices in superconductors, Phys. Rev. 140, A1197 (1965).
- Foltyn et al. (2007) S. Foltyn, L. Civale, J. MacManus-Driscoll, Q. Jia, B. Maiorov, H. Wang, and M. Maley, Materials science challenges for high-temperature superconducting wire, Nat. Mat. 6, 631 (2007).
- Holesinger et al. (2008) T. G. Holesinger, L. Civale, B. Maiorov, D. M. Feldmann, J. Y. Coulter, D. J. Miller, V. A. Maroni, Z. Chen, D. C. Larbalestier, R. Feenstra, et al., Progress in nanoengineered microstructures for tunable high-current, high-temperature superconducting wires, Adv. Mat. 20, 391 (2008).
- Maiorov et al. (2009) B. Maiorov, S. A. Baily, H. Zhou, O. Ugurlu, J. A. Kennison, P. C. Dowden, T. G. Holesinger, S. R. Foltyn, and L. Civale, Synergetic combination of different types of defect to optimize pinning landscape using BaZrO-doped YBaCuO, Nat. Mater. 8, 398 (2009).
- Jia et al. (2013) Y. Jia, M. LeRoux, D. J. Miller, J. G. Wen, W. K. Kwok, U. Welp, M. W. Rupich, X. Li, S. Sathyamurthy, S. Fleshler, A. P. Malozemoff, A. Kayani, O. Ayala-Valenzuela, and L. Civale, Doubling the critical current density of high temperature superconducting coated conductors through proton irradiation, Appl. Phys. Lett. 103, 122601 (2013).
- Miura et al. (2013a) M. Miura, B. Maiorov, J. O. Willis, T. Kato, M. Sato, T. Izumi, Y. Shiohara, and L. Civale, The effects of density and size of BaMO (MZr, Nb, Sn) nanoparticles on the vortex glassy and liquid phase in (Y,Gd)BaCuO coated conductors, Supercond. Sci. Technol. 26, 035008 (2013a).
- Selvamanickam et al. (2015) V. Selvamanickam, M. H. Gharahcheshmeh, A. Xu, E. Galstyan, L. Delgado, and C. Cantoni, High critical currents in heavily doped (Gd,Y)BaCuO superconductor tapes, Appl. Phys. Lett. 106, 032601 (2015).
- Sadovskyy et al. (2015a) I. A. Sadovskyy, Y. Jia, M. Leroux, J. Kwon, H. Hu, L. Fang, C. Chaparro, S. Zhu, U. Welp, J.-M. Zuo, V. Selvamanickam, G. W. Crabtree, A. E. Koshelev, A. Glatz, and W.-K. Kwok, Towards superconducting critical current by design (2015a), arXiv:1509.06446 .
- Blatter et al. (1994) G. Blatter, M. V. Feigel’man, V. B. Geshkenbein, A. I. Larkin, and V. M. Vinokur, Vortices in high-temperature superconductors, Rev. Mod. Phys. 66, 1125 (1994).
- Brandt (1995) E. H. Brandt, The flux-line-lattice in superconductors, Rep. Prog. Phys. 58, 1465 (1995).
- Blatter and Geshkenbein (2003) G. Blatter and V. Geshkenbein, Vortex matter, in The physics of superconductors, edited by K. Bennemann and J. Ketterson (Springer Berlin Heidelberg, 2003) pp. 725–936.
- Larkin and Ovchinnikov (1979) A. I. Larkin and Y. N. Ovchinnikov, Pinning in type II superconductors, J. Low Temp. Phys. 34, 409 (1979).
- Ovchinnikov and Ivlev (1991) Y. N. Ovchinnikov and B. I. Ivlev, Pinning in layered inhomogeneous superconductors, Phys. Rev. B 43, 8024 (1991).
- Blatter et al. (2004) G. Blatter, V. B. Geshkenbein, and J. A. G. Koopmann, Weak to strong pinning crossover, Phys. Rev. Lett. 92, 067009 (2004).
- Ertaş and Kardar (1996) D. Ertaş and M. Kardar, Anisotropic scaling in threshold critical dynamics of driven directed lines, Phys. Rev. B 53, 3520 (1996).
- Bustingorry et al. (2007) S. Bustingorry, L. F. Cugliandolo, and D. Domínguez, Langevin simulations of the out-of-equilibrium dynamics of vortex glasses in high-temperature superconductors, Phys. Rev. B 75, 024506 (2007).
- Luo and Hu (2007) M.-B. Luo and X. Hu, Depinning and creep motion in glass states of flux lines, Phys. Rev. Lett. 98, 267002 (2007).
- Koshelev and Kolton (2011) A. E. Koshelev and A. B. Kolton, Theory and simulations on strong pinning of vortex lines by nanoparticles, Phys. Rev. B 84, 104528 (2011).
- Ortalan et al. (2009) V. Ortalan, M. Herrera, M. W. Rupich, and N. D. Browning, Three dimensional analyses of flux pinning centers in Dy-doped YBaCuO coated superconductors by STEM tomography, Phys. C 469, 2052 (2009).
- Schmid (1966) A. Schmid, A time dependent Ginzburg-Landau equation and its applications to a problem of resistivity in the mixed state, Phys. kondens. Materie 5, 302 (1966).
- Aranson and Kramer (2002) I. S. Aranson and L. Kramer, The world of the complex Ginzburg-Landau equation, Rev. Mod. Phys. 74, 99 (2002).
- Doria et al. (1990) M. M. Doria, J. E. Gubernatis, and D. Rainer, Solving the Ginzburg-Landau equations by simulated annealing, Phys. Rev. B 41, 6335 (1990).
- Machida and Kaburaki (1993) M. Machida and H. Kaburaki, Direct simulation of the time-dependent Ginzburg-Landau equation for type-II superconducting thin film: Vortex dynamics and V-I characteristics, Phys. Rev. Lett. 71, 3206 (1993).
- Crabtree et al. (1996) G. Crabtree, G. Leaf, H. Kaper, V. Vinokur, A. Koshelev, D. Braun, D. Levine, W. Kwok, and J. Fendrich, Time-dependent Ginzburg-Landau simulations of vortex guidance by twin boundaries, Phys. C 263, 401 (1996).
- Aranson et al. (1996) I. Aranson, B. Y. Shapiro, and V. Vinokur, Nucleation and growth of the normal phase in thin superconducting strips, Phys. Rev. Lett. 76, 142 (1996).
- Crabtree et al. (2000) G. W. Crabtree, D. O. Gunter, H. G. Kaper, A. E. Koshelev, G. K. Leaf, and V. M. Vinokur, Numerical simulations of driven vortex systems, Phys. Rev. B 61, 1446 (2000).
- Winiecki and Adams (2002) T. Winiecki and C. S. Adams, Time-dependent Ginzburg-Landau simulations of the voltage-current characteristic of type-II superconductors with pinning, Phys. Rev. B 65, 104517 (2002).
- Vodolazov (2013) D. Y. Vodolazov, Vortex-induced negative magnetoresistance and peak effect in narrow superconducting films, Phys. Rev. B 88, 014525 (2013).
- Berdiyorov et al. (2014) G. Berdiyorov, K. Harrabi, F. Oktasendra, K. Gasmi, A. I. Mansour, J. P. Maneval, and F. M. Peeters, Dynamics of current-driven phase-slip centers in superconducting strips, Phys. Rev. B 90, 054506 (2014).
- Sadovskyy et al. (2015b) I. A. Sadovskyy, A. E. Koshelev, C. L. Phillips, D. A. Karpeyev, and A. Glatz, Stable large-scale solver for ginzburg-landau equations for superconductors, J Comp. Phys. 294, 639 (2015b).
- Lawrence and Doniach (1970) W. E. Lawrence and S. Doniach, in Proceedings of Twelfth Conference on Low Temperature Physics, edited by E. Kanda (Keigaku, Tokyo, 1970) p. 361.
- Long et al. (2005) N. Long, N. Strickland, B. Chapman, N. Ross, J. Xia, X. Li, W. Zhang, T. Kodenkandath, Y. Huang, and M. Rupich, Enhanced in-field critical currents of YBCO second-generation (2G) wire by Dy additions, Supercond. Sci. Technol. 18, S405 (2005).
- Gutierrez et al. (2007) J. Gutierrez, A. Llordes, J. Gazquez, M. Gibert, N. Roma, A. Pomar, F. Sandiumenge, N. Mestres, T. Puig, and X. Obradors, Strong isotropic flux pinning in solution-derived YBaCuO nanocomposite superconductor films, Nat. Mater. 6, 367 (2007).
- Miura et al. (2013b) M. Miura, B. Maiorov, T. Kato, T. Shimode, K. Wada, S. Adachi, and K. Tanabe, Strongly enhanced flux pinning in one-step deposition of BaFe(AsP) superconductor films with uniformly dispersed BaZrO nanoparticles, Nat. Comm. 4, 2499 (2013b).
- Leroux et al. (2015) M. Leroux, K. J. Kihlstrom, S. Holleis, M. W. Rupich, S. Sathyamurthy, S. Fleshler, H. Sheng, D. J. Miller, S. Eley, L. Civale, A. Kayani, P. M. Niraula, U. Welp, and W.-K. Kwok, Rapid doubling of the critical current of YBaCuO coated conductors for viable high-speed industrial processing, Appl. Phys. Lett. 107, 192601 (2015).
- Koshelev et al. (2015) A. E. Koshelev, I. A. Sadovskyy, C. L. Phillips, and A. Glatz, Optimization of vortex pinning by metallic nanoparticles using large-scale numerical simulations (2015), arXiv:1509.04212 .