# Assessing the spatial and field dependence of the critical current density in YBCO bulk superconductors by scanning Hall probes

## Abstract

Although the flux density map of a bulk superconductor provides in principle sufficient information for calculating the magnitude and the direction of the supercurrent flow, the inversion of the Biot-Savart law is ill conditioned for thick samples, thus rendering this method unsuitable for state of the art bulk superconductors. If a thin () slab is cut from the bulk, the inversion is reasonably well conditioned and the variation of the critical current density in the sample can be calculated with adequate spatial resolution. Therefore a novel procedure is employed, which exploits the symmetry of the problem and solves the equations non-iteratively, assuming a planar -independent current density. The calculated current density at a certain position is found to depend on the magnetic induction. In this way the average field dependence of the critical current density is obtained also at low fields, which is not accessible to magnetisation measurements due to the self-field of the sample. It is further shown that an evaluation of magnetisation loops, taking the self-field into account, results in a similar dependence in the field range accessible to this experiment.

PACS: 74.25.Qt, 74.25.Sv, 74.72.Bk, 74.81.Bd

## 1 Introduction

At present, bulk superconductors with several centimeters in diameter and about one centimeter in thickness, trapping remanent magnetic fields exceeding 1 T at 77 K, can be reproducibly grown [1]. Texturing of the monolith is achieved by top seeded melt texture growth (TSMG), where crystallisation evolves from a seed crystal placed on top of the bulk. Five growth sectors are formed, propagating from the facets of the crystal through the entire material [2]. The critical currents achieved in each growth sector, especially as a function of the seed distance, are therefore of particular interest.

A straightforward approach is to cut the sample and to characterise small pieces by magnetometry. However, this procedure is destructive and takes extensive measurement time. On the other hand, scanning Hall probe techniques can be employed to analyse the local properties of the bulk [3]. Among them the magnetoscan technique [5] proved to provide detailed information on the critical current flow on a local scale. However, even the strongest permanent magnets used in the magnetoscan device activate currents only in the uppermost layer of the bulk. The information refers to depths of less than 1 mm [6] and the inside of the bulk cannot be probed.

Currents flowing in the entire volume of the bulk can be activated by performing a (zero) field cooled hysteresis loop in a magnet. In principle sufficient data to calculate the current density on the sub-mm length scale, where substantial changes in are expected, can be obtained from trapped field maps.

Although elaborate and numerically stable techniques exist [7, 8, 9], it was shown [10] that due to the large thickness of those samples the matrix inverted in such procedures is notoriously ill conditioned and can even be numerically singular. Thus, a reduction in thickness, either by grinding or preferably by cutting the bulk into disks, is mandatory to assess the local critical current distribution in this way. If the disk is sufficiently thin, a scanning mesh can be found which allows both an analysis of the current density on the sub-mm scale and a comparison to the magnetoscan signal at certain positions.

## 2 Numerics

### 2.1 General

In the coordinate system used in the following (cf. figure 4) the top sample surface lies in the -plane and the perpendicular component of the magnetic induction is assessed. Similar to [7] the current density is expressed as the curl of a magnetisation density pointing in -direction . This implicitly satisfies current conservation for the -independent planar current distribution . Unlike in [7] the basic equation of the problem is derived by splitting the magnetic induction into two components and using a scalar magnetic potential (see the Appendix for a detailed derivation). Discretisation of the integral equation using cubic volume elements with constant results in the 2D matrix equation

(1) | |||||

(2) | |||||

(3) |

Here denotes the sample thickness, the distance between the active area of the Hall probe and the top sample surface (gap), the step width, and , the indices on the mesh; the antiderivative is evaluated at the eight corners of the cubes in (2).

Equation 1 can be mapped one-to-one to 1D by substituting . The equation now reads

(4) |

and the problem can be tackled by matrix inversion algorithms.

As pointed out earlier, the matrix is in fact a Toeplitz block Toeplitz matrix [7]. This is a consequence of the translation invariance of the Biot-Savart law () and results in a highly symmetric matrix (cf. inset of figure 1). Therefore, both efficient storage and a fast algorithm for solving the system can be expected. However, [7] exploits the symmetry only for the storage of the matrix elements and the method of conjugated gradients using the fast Fourier transform (FFT) is employed to solve the linear equation. Although the procedure is fast, employing FFT [7, 8] implies the unnecessary assumption of periodicity in outside the measurement area, which may create artefacts, if currents are flowing close to the edge of the scanning area.

Block Toeplitz matrices occur in a number of problems, such as image reconstruction or system identification. Fortunately, an efficient and fast algorithm, which exploits the symmetry of the structured matrix, is provided in [11]. The computation time is approximately 10 minutes for a system on a desktop PC. This is presumably longer than FFT based algorithms, but still much less than the actual measurement time.

### 2.2 Condition number

The residuum of the solution is of the order of the machine error. However, it is shown in standard numerical algebra textbooks, that matrix inversions can amplify a relative (measurement) error in the right-hand side of (4), leading to an unknown error in the calculated magnetisation density . This behaviour is described by the condition number of a linear system

(5) |

where the errors are measured in Euclidean norm. It was shown in [10] that the condition number of the inverted matrix is in notoriously high for bulk superconductors. Therefore, special attention has to be paid to the choice of the step width , defining at a given gap . With a realistic distance of about 0.15 mm (see section 3) and a sample thickness of 0.55 mm, the condition number was estimated using the numerical computation package Scilab [12] (cf. Fig. 1).

As a rule of thumb one aims at a condition number such that

(6) |

holds. It was found that the matrix inversion is reasonably well conditioned as long as the step width is larger than the gap. The actual choice of mm is adequate to resolve the spatial variations of the critical current on a sub-mm scale. The condition number combined with a rough estimate of the relative measurement error results in , satisfying (6).

Calculating the current density involves computing a (numerical) derivative, which implies that the relative error in will be high, if the change in (the current at this position) vanishes, as for example outside of the bulk or in defects. This is a peculiarity of measuring the relative error of a quantity close to zero. Note that the absolute error in is bounded by (5) and expected to be acceptably low for most of the current distribution inside the sample as long as (6) holds.

## 3 Experimental

A thin disk was cut from an undoped YBCO bulk superconductor with a diameter of , which was grown by the top seeded melt growth technique [4], using a diamond saw. The cut was made near the upper surface of the bulk and the disk was polished with abrasive paper to a thickness of .

Commercial Hall probes from Arepoc with active areas of m (trapped field map) and m (magnetoscan) were used in the experiments. All scans were carried out with a step width of 0.2 mm, which allows to apply the inversion (see section 2.2).

For the trapped field maps a scanning area of mm was used. The sample was field cooled in a split coil magnet in a field of 1.4 T. In order to minimise relaxation effects during the measurement, the scan was started 10 min after sweeping the magnetic field to zero. The temperature of the liquid nitrogen bath was recorded prior to and after the scans and was found to be 77.2 K, increasing due to oxygen uptake by about 0.1 K during the measurement. The resistive offset of the Hall probe was determined at the end of each measurement, with the Hall probe still immersed in liquid nitrogen, but at a large distance from the bulk.

A somewhat wider scanning area of mm was used for the magnetoscans. This precaution was taken to assure that the magnet would not produce artefacts by stopping over the bulk after finishing a single line of the scan. The magnetoscan was carried out using a small SmCo permanent magnet of 2 mm in diameter, applying an induction of around 150 mT at the top sample surface of the bulk. Further details of the technique can be found in [5].

## 4 Results

### 4.1 Inversion of the Trapped Flux Density Profile

Assuming that the variation of the critical current over the sample thickness can be neglected immediately implies that the trapped field profiles recorded on the bottom and top surface are identical, which was confirmed by the experiment. The maxima of the trapped field were equal within one percent on both sides of the bulk and found to be . Strong negative fields of up to were detected close to the edge of the sample, which result from the high diameter to thickness aspect ratio of the disk (cf. Fig. 2a).

The high reproducibility between several measurements shows that the gap between the Hall probe and the sample surface () remains unchanged in subsequent runs.

The inversion of the trapped flux density profile of the bottom surface of the disk is depicted in figure 2b. Most of the bulk carries a current density of around in the remanent state. The defects (blue arrows) close to the edge appear as regions, where the critical current is drastically reduced, most likely due to cracks or large scale inclusions. A remarkable result of the inversion is the detection of the -growth sector at the center of the bulk. It is clearly displayed in the current density map as a rectangular area with low critical current density (green region).

In addition, a clear negative correlation between the critical current density and the magnitude of the perpendicular magnetic induction is found, for example close to one of the – growth sector boundary, where both the magnetic field and the current density change simultaneously (small white arrows). The correlation is most prominent at low fields, i.e. high current densities of up to flow close to the sample edge, where the magnetic induction changes sign and therefore . Moreover, a point inside the bulk (lowest white arrow) with reduced magnetic induction is reproducibly detected, where the current density significantly exceeds the nearby current densities. This demonstrates that a strong field dependence of the critical current is present, especially at low fields.

### 4.2 Comparison to Magnetoscan

Due to the strong field dependence of the critical current, it is difficult to compare different regions of the bulk, as the self-field of the bulk is position-dependent in the remanent state. Contrary, in the magnetoscan the background field of the permanent magnet is constant and the self-field is smaller, since the currents are activated only in an area of about the magnet’s diameter [6].

All except one of the prominent defects are found in the magnetoscan (cf. Fig. 2c). A possible explanation would be a large defect or inclusion at a depth, which exceeds the penetration depth of the permanent magnet. Similar to the inversion, the reduction of the critical current is evident from the low magnetic response in the -axis growth sector. Moreover, the observable granularity on a sub-mm scale of the bulk’s response indicates strong local variations in the critical current density.

### 4.3 Comparison to Magnetometry

The calculated current densities can be correlated with the -component of the magnetic induction at their position. For this purpose the magnetic induction in the central plane inside the bulk was calculated after the inversion (see the Appendix), resulting in a plot (cf. figure 3). Only points with positive induction were considered, which effectively cuts off the noise due to the extended defects close to the sample edge.

For comparison four cubes () were cut from different positions of the bulk (cf. Fig. 2c) and analysed in a SQUID magnetometer at the temperature of the liquid nitrogen bath during the previous scans. The magnetisation loops were evaluated assuming the Bean model. A large sample to sample variation was found, at this point the lowest critical current density being located at the -axis growth sector in the center of the bulk (cubes #3 and #4), in agreement with the above results. Although numerical errors in the inversion procedure (mainly due to the numerical derivation) cannot completely be excluded, the variation in the magnetisation data elucidates the scatter in the calculated current densities, if one takes into account that the averaged volume is much larger in the SQUID measurements. This suggests that the scatter is primarily due to strong local variations of the critical current in the bulk, which is further supported by the granularity observed in the magnetoscan signal (cf. Fig. 2c).

The average field dependence of the critical current density from the magnetisation loops shows good agreement with obtained from the trapped flux density maps, except at low fields, where the curves clearly differ (cf. Fig. 3a). This can be explained in terms of the sample’s self-field, which is not considered in the evaluation of the magnetisation experiment, since it takes only the externally applied field into account and neglects the field generated by the currents flowing in the sample.

To estimate this contribution the magnetoscan signal serves as a first approximation, because the currents are induced in an area determined by the magnet’s diameter and the volume of current flow mainly contributing to the magnetoscan signal is therefore similar to the cubes used in magnetometry. Indeed the average induction of the magnetoscan () is close to the field, where the and start to differ.

For a more detailed analysis an evaluation method was applied (again based on the Bean model), which accounts for both the externally applied field and the mean self-field in the sample, thus providing an approximate dependence [13]. Especially for low applied fields, where self-field effects become dominant, the evaluation shifts all data points to higher fields (cf. Fig. 3b). This effect is particularly clear for the remanent state at zero applied field (arrow), where the magnetic induction is solely due to the trapped self-field. Consequently, a sample cannot be probed at zero magnetic induction by SQUID magnetisation experiments.

When accounting for the self-field contribution, good agreement between the curve from the SQUID loops and the average curve from the inversion is found in the range accessible to the magnetisation experiment (depicted in figure 3b for one of the cubes). Moreover, the field, where the correction starts to become effective and the self-field becomes important, is equal to the simple estimate made above using the magnetoscan data ().

In contrast to SQUID magnetometry, the current calculation allows one to analyse the dependence in fields ranging from zero induction to the maximum trapped field. A clear dependence of the critical current density on the magnetic induction is revealed in this way. There is no indication for a flattening in the average at low fields.

## 5 Summary

The inversion of the Biot-Savart law represents an ill conditioned problem for bulk superconductors, but parameters, e.g. the step width of scan and the thickness of the sample, can be found, which allow its application to thin disks cut from the sample. The matrix equation was solved without any additional assumptions by a fast algorithm, which exploits the symmetry of the problem. In this way, for example the -axis growth sector was clearly identified as a distinct region of low critical current density, a result, which is also obtained by the magnetoscan and confirmed by magnetisation measurements.

From a technological point of view both methods, the magnetoscan and the inversion of the flux density map, are complementary. The magnetoscan provides information on the local critical current density at a certain constant background field, which is of interest for superconducting bearings. On the other hand, information about the remanent current flow at the bulk’s self-field is important for magnet applications and can be obtained by the inversion of the trapped flux density map of thin disks cut from the sample.

In addition, the current calculation was found to provide the important average field dependence of the critical current also at low fields (below the self-field), a region, which is not accessible to magnetic measurements, even if the self-field is explicitly accounted for in the evaluation procedure. There is no indication of a plateau in and the critical current density is found to depend on the field in a continuous way also at the lowest magnetic inductions.

### Acknowledgement

The authors wish to thank Vasile Sima for valuable discussions. This work was supported in part by the Austrian Science Fund under contract #17443.

## Appendix

The magnetostatic Maxwell equations ( and Ampere’s law ) are treated by splitting the magnetic induction into a sum of two fields

(A-1) |

where is chosen to satisfy the differential equation

(A-2) |

in entire space, which satisfies current conservation . Equation (A-2) defines apart from a gradient field and an arbitrary constant, which are both chosen to be zero. Consequently vanishes outside the sample and can therefore be interpreted as a magnetisation, which is confined to the sample volume and establishes the current density by spatial variations. Ampere’s law results now in a homogeneous equation for the field

(A-3) |

which can thus be derived from a scalar potential

(A-4) |

Taking the divergence of this expression and substituting the second Maxwell equation () results in

(A-5) |

and the problem can be solved by the method of Green’s functions

(A-6) |

The planar -independent current flow can be represented by and consequently the first term vanishes since in the sample volume. Further, only the surfaces perpendicular to contribute to the second term. Using (A-6) the measured induction is expressed solely as a function of

(A-7) |

Here and the notation for evaluating antiderivatives are used to indicate the positive contribution from the top () and the negative from the bottom () surface.

Note, that vanishes for the case of an infinitely long slab (zero demagnetisation) as in (A-7). In this case the induction is simply determined by the magnitude of at a certain position and inside and outside the slab, where . However, for finite geometries there will always be a contribution from in (A-7) and the relation between and is non-local. Therefore the integral equation (A-7) must be solved in order to attain .

The corresponding set of linear equations for the discrete measured data is formulated by approximating the sample as an array of cubes with constant . Summing over all elements results in the matrix equation

(A-8) |

Here, the matrix entries are calculated by evaluating the antiderivative

(A-9) | |||||

(A-10) |

at the eight corners of the cubes

(A-11) |

Once the components are obtained by inverting (A-8), the current density can be calculated by employing (A-2). Further, the induction at any distance outside the sample can be obtained from (A-8). Note, that the magnetic induction in the central plane () of the sample, which is needed to obtain , is given by

(A-12) |

as does not vanish inside the sample.

### References

- Cardwell D A, Hari Babu N, 2006 Physica C 445 1–7
- Cardwell D A, 1998 Mater. Sci. Eng. 53 1–10
- Haindl S, Eisterer M, Hörhager N, Weber H W, Walter H, Shlyk L, Krabbes G, Hari Babu N and Cardwell D A 2005 Physica C 426 625–631
- Krabbes G, Schätzle P, Bieger W, Wiesner U, Stöver G, Wu M, Strasser T, Köhler A, Litzkendorf D, Fischer K and Gr̈onert P 1995 Physica C 244 145–152
- Eisterer M, Haindl S, Wojcik T and Weber H W 2003 Supercond. Sci. Technol. 16 1282–1285
- Zehetmayer M, Eisterer M and Weber H W 2006 Supercond. Sci. Technol. 19 429–437
- Wijngaarden R J, Spoelder H J W, Surdenau R and Griessen R 1996 Phys. Rev. B 17 6742–6749
- Jooss C, Warthmann R, Forkl A and Kronmüller H 1998 Physica C 299 215–230
- Carrera M, Amorós J, Obradors X and Fontcuberta J 2003 Supercond. Sci. Technol. 16 1187–1194
- Eisterer M, 2005 Supercond. Sci. Technol. 18 58–62
- http://www.icm.tu-bs.de/NICONET 2008 According to the online manual subroutine MB02ED uses “Householder transformations, modified hyperbolic rotations and block Gaussian eliminations in the Schur algorithm to solve the linear equation”. All matrices were found to meet the prerequisite of being symmetric and numerically positive definite.
- http://www.scilab.org 2008 The condition number is obtained by computing the singular value decomposition (SVD) of the matrix.
- Wiesinger H P, Sauerzopf F M and Weber H W 1992 Physica C 203 121–128