Evaluating Superconductors through Current Induced Depairing

TL;DR

This paper reviews the fundamentals of current-induced depairing in superconductors, introduces a pulsed current measurement technique, and discusses its application to topological-insulator/chalcogenide superconducting systems, highlighting its importance for understanding superconducting properties.

Contribution

It presents a detailed review of depairing current measurement techniques and applies them to a specific topological superconductor system, advancing experimental methods in the field.

Findings

Depairing current density Jd is crucial for understanding superconducting boundaries.

Fast-pulsed current technique enables measurement of high Jd values.

Application to topological-insulator/chalcogenide systems reveals new insights.

Abstract

The phenomenon of superconductivity occurs in the phase space of three principal parameters: temperature T, magnetic field B, and current density Jd . The critical temperature Tc is one of the first parameters that is measured and in a certain way defines the superconductor. From the practical applications point of view, of equal importance is the upper critical magnetic field Bc2 and conventional critical current density Jc (above which the system begins to show resistance without entering the normal state). However, a seldom-measured parameter, the depairing current density Jd , holds the same fundamental importance as Tc and Bc2, in that it defines a boundary between the superconducting and normal states. A study of Jd sheds unique light on other important characteristics of the superconducting state such as the superfluid density and the nature of the normal state below Tc, information that can play a key role in better understanding newly-discovered superconducting materials. From a measurement perspective, the extremely high values of Jd make it difficult to measure, which is the reason why it is seldom measured. Here, we will review the fundamentals of current-induced depairing and the fast-pulsed current technique that facilitates its measurement and discuss the results of its application to the topological-insulator/chalcogenide interfacial superconducting system. Keywords: pairbreaking, pair-breaking, vortex, vortices, theory, tutorial, RTS, room-temperature supeconductivity

Equations

$${\rho }_f\sim {\rho }_{n}B/B_{c2}$$

$$\xi =v_F{\tau }_{\Delta }\simeq \hslash v_F/\pi \Delta$$

$$B_{c2}=\frac{{\Phi }_0}{2\pi {\xi }^2}=\frac{{\Phi }_0}{2\pi {\xi }_1{\xi }_2}$$

$${\lambda }_L=\sqrt{\frac{m^{\ast }}{{\mu }_0{n}_s{e}^2}}$$

$${\rho }_s\equiv 1/{\lambda }^2$$

$$E=\left(\frac{m^{\ast }}{n_s{e}^2}\right)\frac{\partial j}{\partial t}$$

$${\nabla }^{2}B=B/{\lambda }_L^2$$

$$B=-{\mu }_0{\lambda }_L^2(\nabla \times j)$$

$$f_c=f_n-f_s=\frac{B_c^2}{2{\mu }_0}$$

$$j_d\le \frac{B_c}{{\mu }_0{\lambda }_L}$$

$$f_s=f_n+\alpha |\psi {|}^2+\frac{\beta }{2}|\psi {|}^4,$$

$$|\psi {|}^2=-\alpha /\beta .$$

$$f_c=\frac{B_c^2}{2{\mu }_0}=\frac{{\alpha }^2}{2\beta }$$

$${\lambda }^2=\frac{m^{\ast }}{2{\mu }_0|\psi {|}^2{e}^2}=\frac{-\beta }{\alpha }\left(\frac{m^{\ast }}{2{\mu }_0{e}^2}\right)$$

$$\alpha =-\frac{2e^2{B}_c^2{\lambda }^2}{m^{\ast }}\text{mbox}and\beta =\frac{4{\mu }_0{e}^4{B}_c^2{\lambda }^4}{m^{\ast 2}}$$

$$f_s=f_n+\alpha |\psi {|}^2+\frac{\beta }{2}|\psi {|}^4+|\psi {|}^2{m}^{\ast }{v}_s^2.$$

$$|{\psi }_{j\ne 0}{|}^2=\frac{-\alpha }{\beta }\left(1-\frac{m^{\ast }{v}_s^2}{|\alpha |}\right)=|{\psi }_{j=0}{|}^2\left(1-\frac{m^{\ast }{v}_s^2}{|\alpha |}\right)$$

$$j=2e|{\psi }_{j\ne 0}{|}^2{v}_s=\frac{-2e\alpha }{\beta }\left(1-\frac{m^{\ast }{v}_s^2}{|\alpha |}\right)v_s$$

$$j_d(T)=\frac{-4e\alpha }{3\beta }{\left(\frac{|\alpha |}{3m^{\ast }}\right)}^{1/2}={\left(\frac{2}{3}\right)}^{3/2}\frac{B_c(T)}{{\mu }_0\lambda (T)}$$

$$j_d(T)\approx j_d(0)[1-(T/T_c{)}^2{]}^{\frac{3}{2}}[1+(T/T_c{)}^2{]}^{\frac{1}{2}}$$

$$j_d(T)\approx \sqrt{2}{j}_d(0)[1-(T/T_c{)}^2{]}^{\frac{3}{2}}.$$

$$B_c=\sqrt{\frac{{\Phi }_0{B}_{c2}}{4\pi {\lambda }^2}}$$

$$j_d(0)=\sqrt{\frac{2{\Phi }_0{B}_{c2}(0)}{27\pi {\mu }_0^2{\lambda }^4(0)}}$$

Full text

vortex, vortices, theory, , tutorial, review, RTS, room-temperature supeconductivity

††preprint: Condens. Matter 4, 54 (2019). Review article in Special Issue ”From Cuprates to Room Temperature Superconductors”

The phenomenon of superconductivity has a long and rich history: from the initial discovery in 1911 of superconductivity in mercury at liquid-helium temperature onnes to the recent discovery of room-temperature superconductivity in lanthanum superhydride soma ; drozdov . Numerous parameters, probed by a variety of techniques, are used to characterize the superconducting state. However, the mixed-state upper critical field $B_{c2}$, reflective of the coherence length $\xi$, and the penetration depth $\lambda_{L}$, reflective of the superfluid density $\rho_{s}$=$1/\lambda^{2}$, are two crucial measurements that are amongst the first to be performed. There are multiple techniques for determining such parameters, each of which has its own advantages and limitations. Our group has developed some uncommon, and in some cases unique, experimental techniques that investigate superconductors at ultra-short time scales, and under unprecedented and extreme conditions of current density $j$, electric fields $E$, and power density $p=\rho j^{2}$ (where $\rho$ is the resistivity). These techniques have led to the discovery or confirmation of several novel phenomena and regimes in superconductors and in addition provide an alternative method to glean information on fundamental superconducting parameters, which in some cases may be hard to obtain by other methods. These methods and approaches are highly relevant in the search for new superconducting materials and in developing an understanding of their fundamental properties. This article discusses the physical meaning of $j_{d}$ and its interrelationships with other basic parameters of the superconducting state, as well as the technical challenges in measuring this important critical parameter. We discuss our results from this approach in the study of the topological-insulator/chalcogenide interfacial superconducting system.

I Introduction

Attractive interactions between charge carriers cause them to condense by pairs into a coherent macroscopic quantum state below some transition temperature $T_{c}$. The formation of this state is governed principally by a competition between four energies: condensation, magnetic-field expulsion, thermal, and kinetic. The order parameter $\Delta$, which describes the extent of condensation and the strength of the superconducting state, is reduced as the temperature $T$, magnetic field $B$, and electric current density $j$ are increased. In type-II superconductors, there is partial flux entry at $B$ values above the lower critical magnetic field $B_{c1}$ and complete destruction of superconductivity above the upper critical field $B_{c2}$ (type-I superconductors can be viewed as a special case where the thermodynamic critical field $B_{c}$ = $B_{c1}$ = $B_{c2}$). The boundary in the $T$-$B$-$j$ phase space that separates the superconducting and normal states is where $\Delta$ vanishes, and the three parameters attain their critical values $T_{c2}(B,j)$, $B_{c2}(T,j)$, and $j_{d}(T,B)$. $j_{d}$ sets the intrinsic upper limiting scale for supercurrent transport in any superconductor, and for $j>j_{d}$, the system attains its normal-state resistivity $\rho_{n}$. $j_{d}$ should be distinguished from the conventional critical value $j_{c}$ (related to extrinsic characteristics such as the depinning of vortices) above which there is partial resistivity $\rho<\rho_{n}$.

The resistivity $\rho$ in the superconducting state is usually less than its normal-state value $\rho_{n}$. The reason for the presence of resistance at all in the superconducting state is because of fluctuations, percolation through junctions (in the case of granular superconductors), and the motion of magnetic flux vortices. For singly-connected superconductors not very close to $T_{c}$, only the last mechanism dominates as the cause of resistance. In the magnetic field region between $B_{c1}$ and $B_{c2}$, a type II superconductor enters a “mixed state” with quantized magnetic flux vortices, each containing an elementary quantum of flux $\Phi_{0}=h/2e$. Under the Lorentz driving force of an applied current, \bm{j}$$\times$$\Phi_{0}, vortices move transverse to $j$, leading to a flux-flow resistivity:

[TABLE]

in the free-flux-flow (large driving force) limit.

Two length scales characterize the superconducting state tinkhamtext . One is the coherence length:

[TABLE]

which is the characteristic length scale for spatial modulations in $\Delta$ (here, $v_{F}$ is the Fermi velocity and $\tau_{\Delta}$ is the order-parameter relaxation time). The normal core of a flux vortex has an approximate effective radius of $\xi$. The destruction of the superconducting state occurs when these normal cores overlap, corresponding to the condition:

[TABLE]

where $\xi$ is the coherence length perpendicular to $B$; the single $\xi$ is replaced by the geometric mean $\sqrt{\xi_{1}\xi_{2}}$ in cases where the plane perpendicular to $B$ is characterized by two anisotropic values.

The other characteristic length scale in a superconductor is the magnetic-field penetration depth $\lambda$, whose London value is given by:

[TABLE]

where $m^{*}$ is the effective electronic mass and $n_{s}$ is the density of superconducting electrons. The theory behind this important quantity and its relationship to $j_{d}$ is described below. Figure 1 shows the profile of the magnetic field as it gets screened from the interior of a superconductor.

I.1 Superfluid Density

For clean metallic superconductors, $n_{s}\rightarrow n$ as $T\rightarrow 0$, where $n$ is the concentration of carriers in the normal state. Non-local effects and other corrections lead to deviations in $\lambda$ from its London value. Hence, the superfluid density $\rho_{s}$ can be conveniently and more completely defined as:

[TABLE]

which includes the effective mass and other corrections to the effective $n_{s}$, rather than the simpler definitions $\rho_{s}=n_{s}$ or $\rho_{s}=n_{s}/m^{*}$ that are sometimes used in the literature. $\rho_{s}$, a quantity of central importance in superconductivity, characterizes the phase stiffness of the condensate emery_kivelson_nature and its effectiveness at screening out magnetic fields and feeds into expressions for the transition temperature (such as the Uemura relation $T_{c}\propto\rho_{s}(0)$ that applies to the underdoped cuprate superconductors Uemura1 ; Uemura2 ; Hardy ).

Traditionally, many common methods for obtaining $\rho_{s}$ do so by directly or indirectly measuring $\lambda$ through its effect on a superconducting sample’s magnetic-field profile and consequent magnetic susceptibility. This category includes methods such as reflection of spin-polarized slow neutrons Felcher , mutual inductance altered by an intervening superconducting film Claassen ; yong-lemberger , changes in the self-inductance of a coil that is part of an LC resonating circuit Boghosian ; Degrift , muon-spin rotation Sonier , magnetic force microscopy Luan , microwave cavity resonance Shibauchi , and measurements of the lower critical field. These measurements are understandably affected if the material’s internal magnetic field is altered, for example by a large paramagnetic background as in the case of the Nd2-xCexCuO4 superconductor because of its Nd3+ magnetic moments.

Another approach to obtaining $\rho_{s}$ is by measuring the inertia of the superfluid (kinetic inductance) during its ballistic acceleration phase lee-lemberger ; ballistic ; Diener ; Jochem . This method requires the sample to be patterned into very high aspect ratio meanders for the highest accuracy.

A measurement of $j_{d}$ provides an alternative to the above approaches for obtaining $\rho_{s}$. It requires a minimal amount of material (typically just a microbridge or nanobridge), does not require the complicated meander patterning needed for a kinetic-inductance measurement, and is unaffected by a material’s normal-state magnetism that affects inductive measurements of $\rho_{s}$ as discussed above. This immunity to material magnetism was used to good advantage for directly obtaining $\rho_{s}$ in the Nd2-xCexCuO4 superconductor for the first time ncco . Furthermore, unlike some of the methods for measuring $\lambda$ that do not provide an accurate absolute value but only provide the temperature variation $\lambda(T)/\lambda(0)$, $j_{d}$ does provide the absolute value of $\lambda$ and $\rho_{s}$ and can hence provide information on the total carrier concentration $n$.

I.2 Normal-State Resistivity

Another valuable byproduct of measuring $j_{d}$ is that it provides a direct measurement of the normal state resistivity $\rho_{n}$ for temperatures below $T_{c}$. One of the starting points in developing an understanding of any newly-discovered superconductor is to understand the underlying normal state: the type of carriers and their concentration, their band-related properties, and the relevant scattering mechanisms and rates. At low applied $B$ and $j$, \rho$$(T) drops precipitously below $T_{c}$, thereby obscuring how $\rho_{n}(T\ll T_{c})$ would have behaved if the superconductivity had not set in. The most common method for measuring $\rho_{n}(T<T_{c})$ utilizes high magnetic fields $B>B_{c2}$ to drive the system normal below $T_{c}$; however, this measurement is subject to magnetoresistance (typically $R(B)\neq$ const) and may require prohibitively high magnetic fields ($B_{c2}>100$ T for some superconductors).

One alternative is to use the core of a magnetic flux vortex as a window to the normal state. From Equation (1), a measurement of $\rho_{f}$ elucidates $\rho_{n}$ metal . However, this extraction of $\rho_{n}$ requires interpretation and modeling, since Equation (1) usually holds only approximately except for very high driving forces, and the exact prefactor depends on the detailed regime of flux flow lo ; blatter ; hote .

Current-induced depairing provides an especially clean method for destroying superconductivity and accessing $\rho_{n}(T\ll T_{0})$. Like the method of applying $B>B_{c2}$ to drive the system normal, applying $j>j_{d}$ is also free of interpretation and modeling, unlike flux-flow dissipation measurements. On the other hand, unlike the potential errors in the $B_{c2}$-based measurement due to normal-state magnetoresistance, the $j_{d}$ method is immune to this issue because the normal-state electroresistance is quite negligible (i.e., $R(E)\simeq$ const) under the electric fields that arise at depairing conditions.

Thus, besides the investigation of interesting phenomena and regimes in superconductivity related directly to current-induced depairing itself, the study of $j_{d}$ provides information on the important parameters of the superconducting state such as \rho_{s}$$(T) and \rho_{n}$$(T).

II Relationship between the Depairing Current and Other Parameters

In a microscopic theory such as the Bardeen–Cooper–Schrieffer (BCS) theory, experimental quantities are calculated from microscopic parameters such as the strength of the effective attractive interaction that leads to Cooper pair formation and the density of states at the Fermi level. Often, these microscopic parameters are not sufficiently well known. In the London and Ginzburg–Landau (GL) phenomenological theories, connections are made between the different observables from constraints based on thermodynamic principles and electrodynamical properties of the superconducting state, leading to an adequate estimation of the depairing current. These phenomenological formulations are described next tinkhamtext ; pbreview .

II.1 London Formulation

The London theory london ; tinkhamtext of superconductivity provides a description of the observed electrodynamical properties by supplementing the basic Maxwell equations by additional equations that constrain the possible behavior to reflect the two hallmarks of the superconducting state: perfect conductivity and the Meissner effect. Note that these properties hold only partially when vortices are present.

An ordinary metal (normal conductor) requires a driving electric field $E$ to maintain a constant current against resistive losses. In the simple Drude picture, this produces Ohm’s law behavior, $j=\sigma E$, with a conductivity given by $\sigma=ne^{2}\tau/m^{*}$. A superconductor can carry a resistanceless current, and so, an electric field is not required for maintaining a persistent current. Instead, $E$ in a perfectly-conducting state causes a ballistic acceleration of charge so that:

[TABLE]

This is the first London equation, which reflects the dissipationless acceleration of the superfluid.

The second property that needs to be accounted for is the expulsion of magnetic flux by a superconductor. The magnetic field is exponentially screened from the interior following a spatial dependence:

[TABLE]

Together with the Maxwell equation $\nabla\times B=\mu_{0}j$, this implies the following condition between $B$ and $j$:

[TABLE]

This is the second London equation, which describes the property of a superconductor to exclude magnetic flux from its interior. Taken together with the Maxwell equation $\nabla\times E=-\partial B/\partial t$, Equations (6) and (8) yield the expression for $\lambda_{L}$ of Equation (4).

Besides the London equations themselves, a third ingredient needed for the estimation of $j_{d}$ in this framework is the thermodynamic critical field $B_{c}$ and its relationship to the Helmholtz free energy density $f$. When flux is expelled, the free energy density is raised by the amount $B^{2}/2\mu_{0}$. The critical flux expulsion energy (for the ideal case of a type-I superconductor with a non-demagnetizing geometry and dimensions large compared to the penetration depth) corresponds to the condition:

[TABLE]

where the L.H.S. of the equation represents the condensation energy density, which is the difference in free energy densities $f_{n}-f_{s}$ between the normal and superconducting states. $j_{d}$ represents the condition when the kinetic energy density equals the condensation energy density: $\frac{1}{2}n_{s}m^{*}v_{s}^{2}=\frac{m^{*}j_{d}^{2}}{2n_{s}e^{2}}=\frac{B_{c}^{2}}{2\mu_{0}}$, where $v_{s}$ is the superfluid speed. Substituting for $\lambda_{L}$ (Equation (4)) gives the London estimate for the depairing current density:

[TABLE]

The inequality reflects the fact that $n_{s}$ does not remain constant, but diminishes as $j$ approaches $j_{d}$.

II.2 Ginzburg–Landau Formulation

There are situations where a system’s quantum wavefunction cannot be solved for by usual means because the Hamiltonian is unknown or not easily approximated. The GL formulation GLpaper is a clever construction that allows useful information and conclusions to be extracted in such a situation where one cannot solve the problem quantum mechanically. For describing macroscopic properties, such as $j_{d}$ that we are about to calculate, the GL theory is in fact more amenable than the microscopic theory tinkhamtext ; bardeen .

The idea is to introduce a complex phenomenological order parameter (pseudo wavefunction) $\psi=|\psi|e^{i\varphi}$ to represent the superconducting state. $|\psi(r)|^{2}$ is assumed to represent the order parameter $\Delta$ introduced earlier, and to approximate the local density of paired superconducting charge carriers (Cooper pairs), which in turn is half the density of superconducting electrons $n_{s}$.

The free energy density $f_{s}$ of the superconducting state is then expressed as a reasonable function of $\psi(r)$ plus other energy terms. A “solution” to $\psi(r)$ is now obtained by the minimization of free energy rather than through quantum mechanics. The unknown parameters of the theory are then solved in terms of measurable physical quantities, thereby providing constraints between the different quantities of the superconducting state.

Close to the phase boundary, $|\psi|$2 is small, and so, $f_{s}$ can be expanded keeping the lowest two orders of $|\psi|$2. First, let us consider the simplest situation where there are no currents, gradients in $|\psi|$, or magnetic fields present. Then, we have:

[TABLE]

where $\alpha$ and $\beta$ are temperature-dependent coefficients whose values are to be determined in terms of measurable parameters. The coefficients can be determined as follows. First of all, for the solution of $|\psi|^{2}$ to be finite at the minimum free energy, $\beta$ must be positive. Second, for the solution of $|\psi|^{2}$ to be non-zero, $\alpha$ must be negative. Since $|\psi|^{2}$ vanishes above $T_{c}$, $\alpha$ must change its sign upon crossing $T_{c}$. The minimum in $f_{s}$ occurs at:

[TABLE]

Substituting this back in Equation (11) and using the definition of $B_{c}$ (Equation (9)), Equation (12) can be written as:

[TABLE]

giving one of the connections between $\alpha$ and $\beta$ and a measurable quantity ($B_{c}$). A second connection can be obtained by noting that $n_{s}$ in Equation (4) can be replaced by 2$|\psi|$2, taking its equilibrium value from Equation (12):

[TABLE]

Solving Equations (13) and (14) simultaneously gives the GL coefficients:

[TABLE]

Note that $e$ and $m^{*}$ refer to single-carrier values and not pair values.

To calculate $j_{d}$, we include the effect of a current in Equation (11) by adding a kinetic energy term $\frac{1}{2}n_{s}m^{*}v_{s}^{2}=|\psi|^{2}m^{*}v_{s}^{2}$ to it:

[TABLE]

For zero $j$ and $v_{s}$, we saw earlier (Equation (12)) that the equilibrium value of $|\psi|^{2}$ that minimizes the free energy is $|\psi_{j=0}|^{2}=-\alpha/\beta$. For a finite $j$ and $v_{s}$, minimization of Equation (16) gives the value of $|\psi|^{2}$ when it is suppressed by a current:

[TABLE]

The corresponding supercurrent density is then:

[TABLE]

The maximum possible value of this expression can now be identified with $j_{d}$:

[TABLE]

where the GL-theory parameters were replaced by their expressions in terms of the physical measurables $B_{c}$ and $\lambda$ through Equation (15). As anticipated at the end of Equation (10) for the London derivation for $j_{d}$, that simpler estimate is indeed larger than this more rigorous GL derivation by the factor $(3/2)^{3/2}$=1.84.

The approximate temperature dependence of $j_{d}$ can be obtained by inserting the generic empirical temperature dependencies $B_{c}(T)\approx B_{c}(0)[1-(T/T_{c})^{2}]$ and $\lambda(T)\approx\lambda(0)/\sqrt{[1-(T/T_{c})^{4}]}$, giving:

[TABLE]

which close to $T_{c}$ reduces to:

[TABLE]

where $j_{d}(0)$ is given by Equation (19) by setting $T=0$ (for high scattering “dirty” superconductors, the $\sqrt{2}$ prefactor can be smaller or absent bardeen ; romijn ).

Since $B_{c}$ is not an easy quantity to measure directly, the relation:

[TABLE]

can be used along with Equation (19) to write the expression for j_{d}$$(0):

[TABLE]

that has the more easily measurable $B_{c2}$. Since both $B_{c2}$ and $j_{d}$ can be obtained from transport measurements, this becomes a convenient way to obtain $\lambda$ and, hence, $\rho_{s}$.

II.3 Microscopic Formulations and Generalizations

Various authors have calculated $j_{d}(T)$ from a microscopic basis bardeen ; maki ; ovchinnikov . For arbitrary temperatures and mean free paths, one must use the Gorkov equations as the starting point. Kupriyanov and Lukichev KL have derived $j_{d}(T)$ from the Eilenberger equations, which are a simplified version of the Gorkov equations. This derivation is beyond the scope of the present review, but a nice shortened version can be found in romijn . The microscopic calculation confirms the overall temperature dependence predicted by GL, and the two normalized curves differ only slightly from each other (e.g., see Figure 4 of romijn ). Thus, the GL theory can be applied over the entire temperature range down to $T\ll T_{c}$. The previous equations relating $j_{d}$ to $B_{c}$ and $\lambda$ are expected to hold in the case of multiple bands and other gap symmetries, as long as one uses the actual empirical temperature dependencies of $B_{c}$ and $\lambda$, which account for modifications in these unconventional cases. This was experimentally demonstrated in the case of MgB2 pbreview , which was recognized as a multi-gap superconductor tuned by strain and doping in the early part of this century; in fact, MgB2 showed superconductivity near a Lifshitz transition as in iron-based superconductors bauer ; agrestini ; kagan .

III Pulsed Measurement Technique

Depairing current densities in superconductors is extremely high: on the order of $j_{d}$($T$=0) = $10^{11}$–$10^{13}$ A/$m^{2}$. If the cross-section of the sample is even as narrow as just 1 mm2, the current required would reach a value of $I=jA\sim 10^{6}$ A. Such a magnitude of current would be exceedingly difficult to produce and control. There are three steps to overcoming this dilemma: (1) Fabricate samples with very narrow cross-sectional areas. This can be achieved by growing nanowires and nanorods or by depositing very thin films and using lithography to pattern narrow bridges (alternatively, the films can be deposited onto nanowires or carbon nanotubes). (2) The next step is by pulsing the current at very low duty cycles so that large values of $I$ can be handled while reducing the time-averaged current and time-averaged power dissipation to manageable levels. (3) The last step is limiting the measurement of $j_{d}$ to the regime close to $T_{c}$. From Equation (21), it would seem that $j_{d}$ can be made arbitrarily small by making $T$ very close to $T_{c}$; however, the $T-T_{c}$ distance needs to be large compared to the transition width for the measurement to be meaningful. Even for this near-$T_{c}$ measurement of $j_{d}$, the current usually will have to be pulsed to avoid significant sample heating. Furthermore, the near-$T_{c}$ measurement will only measure $\rho_{s}$ in that region, and its zero-$T$ value will have to be extrapolated using theory. While this is better than nothing, it will not shed light on any abnormal temperature dependence of $\rho_{s}$ over the entire range, which could be of special interest if the superconductor has some exotic behavior.

Thus, the experimental ingredients needed to conduct a $j_{d}$ measurement are: a superconducting sample with a very narrow cross-section; a means to control the temperature, i.e., a cryostat; and a method for sourcing pulsed signals (current or voltage) and detecting the consequent complementary signal (voltage or current). There are numerous methods for sample fabrication, which vary widely with the different superconducting materials. Some deposition systems for preparing superconducting films can be bought off the shelf. Cryostats also represent standard equipment that can be bought off the shelf. The principal distinguishing the experimental capabilities of our work center on the pulsed electrical measurements. Therefore, the rest of the experimental section will be devoted to describing this unique measurement setup.

Figure 2 shows the overall configuration and functional schematic. The pulsed current/voltage source puts out a time-varying current and voltage. This signal flows through a standard impedance, usually a resistor $R_{std}$ (although an inductor is preferable in some situations) and the superconducting sample of resistance $R$ that are in series. The initial signal can be taken directly from the output of a standard pulse generator (one of the models used was a Wavetek Model 801). These signal generators will typically have an output impedance of $Z_{out}=50\ \Omega$. If a lower $Z_{out}$ is desirable (to allow for constant voltage control), the signal generator’s output can be passed through any standard buffer amplifier (e.g., a transistor-emitter-follower-based circuit, a power-operational-amplifier-based circuit, or an off-the-shelf audio amplifier). If a higher voltage than the signal generator’s output is desirable, its output can be passed through any standard voltage amplifier (fast high quality audio amplifiers can serve this purpose as well). Combining a higher voltage signal with a large series resistor (which can be the $R_{std}$ itself or an additional series resistor) can provide a relatively constant current. In general, the measurement will be in current-controlled or voltage-controlled mode depending on whether the combination of the final $Z_{out}$ (after the amplifier if any) plus $R_{std}$ is greater than or less than $R$. If the current needs to be held constant to a high accuracy (for example, if a series of $R$ vs. $T$ resistive transition curves needs to be traced out at various constant currents, as will be seen later), then it is better to follow the pulse generator with a transconductance amplifier, which converts the generator’s voltage pulse into a constant current pulse. The transconductance amplifier is able to hold the current constant by electronic circuity instead of needing an enormous series resistance. While conducting a pulsed current-voltage ($IV$) curve, which is usually done manually, it is preferable to have the voltage-controlled mode. The reason for this is that as the current and voltage are pushed higher, the sample’s resistance will increase, and at some point, the sample will be driven to normal as $j_{d}$ is exceeded. In constant-current mode, the power dissipation $P=I^{2}R$ rises as $R$ rises, causing an increase in heating and a further rise in $R$. This can lead to a run-away condition, which can destroy the sample. On the other hand, the constant-voltage mode is self-stabilizing since in this case, $P=V^{2}/R$ decreases as $R$ rises, thus reducing heating and controlling the situation.

Once the current pulse flows through the sample and $R_{std}$, the corresponding time-varying voltages, $V(t$) and $V_{std}(t)$, will be developed across them respectively. These must be observed and quantified using an oscilloscope. A digital storage oscilloscope (DSO) allows multiple pulses to be averaged. Since the signal is exactly repetitive, because the DSO is triggered off of the pulse generator’s sync signal, the averaging effectively suppresses random uncorrelated noise. As long as the sample condition ($T$, $B$, etc.) is stable, a very high number of averages can be taken to improve the signal-to-noise ratio (SNR) vastly. Coaxial cables with 50-Ohm characteristic impedance are used between all connection points, including the wiring within the cryostat. Where possible, the originating and/or terminating points at the ends of the cables need to have matching 50-Ohm values to avoid reflections. Multiple ground connections to the circuit must be avoided to prevent ground loops. This means the two channels of the DSO cannot be simultaneously connected to both the sample and $R_{std}$; either a differential instrumentation preamplifier (Princeton Applied Research and Stanford Research Systems are two brands that make instrumentation amplifiers) must be used between the DSO channels and the sample and $R_{std}$, or only one of the two must be measured at a time.

Figure 3 shows the pair of time-varying current $I(t)$=$V_{std}(t)/R_{std}$ and sample-voltage $V(t)$ signals that results. The topmost trace is the scaled calculated resistance $50R(t)=50V(t)/I(t)$. Note that the pulses reach constant plateaus after their initial transients. $R$, $V$, and $I$ are defined by taking the plateau values of the individual quantities. The thermal rise in a sample because of Joule heating involves several processes: thermal diffusion occurs within the sample essentially instantaneously; on the time scale of nanoseconds, phonons transfer heat across the interface between the film and substrate; heat then diffuses within the substrate in a matter of microseconds and finally into the heat sink in milliseconds. For those processes that have time scales comparable to or longer than the pulse duration, there will be a visible rise in $V(t)$, causing the pulse to be distorted. Thus, as long as the $V(t)$ pulse is flat, slow causes of heating that influence the $V(t)$ shape can be assumed to be negligible. The work in pair discusses a method to evaluate a sample’s thermal resistance for pulsed signals quantitatively.

Figure 4 shows an example pbreview of a set of IV curves at various fixed temperatures (in zero magnetic field), where each data point represents a pulsed measurement (plateau values) as described above. As the temperature is increased, $j_{d}$ is reduced, and hence, the “jump” occurs at a lower value of $I$. Notice that the resistance (the V/I slope) jumps from zero (dissipationless superconducting state) to a constant finite value (normal-state) as the current crosses its depairing value. This is one direct way of obtaining $\rho_{n}$ below $T_{c}$. In this particular material, high impurity scattering dominates over electron-phonon scattering at all temperatures, leading to a relatively flat $\rho(T)$. A more interesting application of this technique for elucidating a variable $\rho(T)$ is described in a later section.

Figure 5a represents a set of pulsed constant-current $R(T)$ curves in zero magnetic field slcojd . As the current is increased, the transition is progressively pushed down in temperature. Figure 5b plots these midpoint transition temperatures ($T_{c2}$) versus $I$, and they are seen to follow a $I^{2/3}$ law as per Equation (21). From this measured slope and Equation (21), one can estimate $j_{d}(0)$ without requiring the application of this enormous value of current. This is especially useful for systems (e.g., cuprate high-temperature superconductors) that have a very high $j_{d}(0)$ value.

Figure 6 shows the companion very low DC-current $R(T)$ curves at various constant $B$ values. Here, the shift occurs because of the $B_{c2}(T)$ boundary; the current level is small enough for its depairing to be negligible. Unlike $j_{d}(T)$, $B_{c2}(T)$ has a linear dependence near $T_{c}$ , and that slope can be related to $B_{c2}(0)$ through the WHH (Werthamer, Helfand, and Hohenberg) theory whh and its variations gurevich-Bc2 by relationships such as $B_{c2}(0)\simeq 0.7dB_{c2}/dT$.

The measurements represented by Figure 5 and Figure 6 together with Equation (23) are the key to obtaining $\rho_{s}$ through relatively straightforward transport measurements. We now look at one recent example of a current-induced depairing study of an exotic superconducting system.

IV Investigations in a Topological Insulator/Chalcogenide Interfacial Superconductor

IV.1 Background

The interface between the Bi2Te3 topological insulator and the FeTe chalcogenide provides a fascinating 2D superconducting system, in which neither Bi2Te3, nor FeTe are superconducting by themselves he . While the exact origin of the superconductivity is not known, it has been suggested that the robust topological surfaces states may be doping the FeTe and suppressing the antiferromagnetism in a thin region close to the interface, thus inducing the observed 2D superconductivity. These surface states represent a conducting system with very high normal conductivity because of protection against time-reversal-invariant scattering mechanisms. Therefore, it is of great interest to understand the nature and origin of the charge carriers that underlie this interfacial superconductivity, and in particular, to see if the topologically-protected surface states might be a source of the normal carriers. The relevance of this question is broader than the specific system studied here, since it has been recently proposed that interfacial superconductivity may even play a role in cuprates: for example, in the interface located between charge density wave nanoscale puddles campi2015 and between oxygen-rich grains where the interface is made of a filamentary network with hyperbolic geometry jarlborg2019 ; campi2016 . We describe below how the current-induced depairing approach was used to answer these questions to elucidate the nature of the normal state in the Bi2Te3/FeTe system.

IV.2 Samples and Experimental Information

The Bi2Te3/FeTe samples consisted of a ZnSe buffer layer (50 nm) deposited on a GaAs (001) semi-insulating substrate, followed by a deposition of 220 nm thick FeTe, which was then capped with a 20 nm-thick Bi2Te3 layer. Upper-critical-field measurements he and vortex-explosion measurements FTBTVortexExplosion showed that the superconductivity occurred within an interfacial layer of thickness $d=7$ nm, which was much thinner than both the FeTe and Bi2Te3 layers. Projection photolithography followed by argon-ion milling was used to pattern narrow microbridges optimized for the high current-density pulsed four-probe measurements. Two bridges were studied: Sample A with lateral dimensions of width $w=11.5\ \mu$m and length $l=285\mu$m and Sample B with $w=12\ \mu$m and $l=285\mu$m. The onset $T_{c}$ (defined as the intersection of the extrapolation of the normal-state portion and the extrapolation of the steep transition portion of the $R(T)$ curve) for both bridges, was 11.7 K. Details about sample preparation are provided in he . All measurements were made in zero applied magnetic field. While the very low reference curves at $I\leq 60\ \mu$A were measured using continuous DC signals, the main electrical transport measurements were made with pulsed signals. Contact resistances ($<1\ \Omega$) were much lower than the normal resistance $R_{n}$ of the bridge, and heat generated at contacts did not reach the bridge within the time duration $t$ of each pulse, since the thermal diffusion distance ($\sqrt{Dt}\sim 10\mu$m) was much shorter than the contact-to-bridge distance ($>1$ mm); here, $D$ is the diffusion constant.

IV.3 Results and Discussion

The normal-state resistivity $\rho_{n}(T)$ and depairing current density $j_{d}(T)$ in the Bi2Te3/FeTe samples were extracted over the entire temperature range jdbite , by driving the system normal with high pulsed currents using the methods described earlier and illustrated in Figure 4 and Figure 5. Figure 7 shows the raw depairing current results. The dashed horizontal lines in Panels (a) and (b) provide the values $I_{d}(T\rightarrow 0)\geq 0.131$ A and 0.136 A for Samples A and B, respectively.

In order to obtain more accurate intrinsic $j_{d}$ and $\rho_{n}$ of the 7 nm-thick superconducting interfacial layer itself, we needed to subtract the small parallel current through the normally conductive underlying FeTe layer. For this purpose, a separate measurement of pure FeTe deposited on ZnSe/GaAs, without the Bi2Te3 top layer, was conducted jdbite . With this subtraction, the previous raw $I_{d}(T\rightarrow 0)$ values gave a corrected j_{d}$$(T=0) of $1.5\times 10^{8}$ A/cm2 for both samples (which is a typical value: $j_{d}$ ranges $10^{7}$–$10^{9}$ A/cm2 for most superconductors), and the correction gave the intrinsic $\rho_{n}(T)$ for the two samples, as shown in Figure 8. This absolute value of \rho_{n}(T$$\rightarrow$$0)$$\sim200 n$\Omega$ cm represents an extraordinarily conductive normal state for a superconducting system, as most superconductors are poor conductors in the normal state. This information will be analyzed below within the framework of an anisotropic Ginzburg–Landau (GL) approach jdbite , to obtain information on the superfluid density, carrier concentration, and scattering rate, as well as their implications for the nature of the normal-state.

From the previously-published measurements of He et al. he , we have the following orientation-dependent values of $B_{c2}$: perpendicular-to-interface $B_{c2}^{\perp}(0)\approx 17$ T and parallel-to-interface $B_{c2}^{\parallel}(0)\approx 40$ T. The corresponding coherence lengths from Equation (3) are: in-plane $\xi_{\parallel}(0)\approx 4.4$ nm and perpendicular $\xi_{\perp}(0)\approx 1.9$ nm. Using Equation (23) together with this $B_{c2}^{\perp}(0)$ and our measured in-plane $j_{d}^{\parallel}(0)$ gave $\lambda_{\parallel}(0)=124$ nm and a corresponding $\rho_{s}(0)=1/\lambda^{2}_{\parallel}(0)$. From $\rho_{s}(0)=\mu_{0}n_{s}(0)e^{2}/m^{*}\approx\mu_{0}ne^{2}/m$ applicable in the clean limit at $T=0$, we get $n\approx 1.8\times 10^{21}$ per cm3, approximating $m^{*}\approx m$. This effective single-band value of $n$ evaluated above is similar to $n$ characteristic of high temperature superconductors and about two orders of magnitude lower than $n$ in highly-conductive metals such as copper.

The low value of $n$ together with the very high normal conductivity implies a rather long scattering time $\tau$ and mean-free-path $l$. The Fermi wave number for this $n$ computes to $k_{F}(3D)=m^{*}v_{F}/\hbar=(3\pi^{2}n)^{1/3}=3.8\times 10^{9}$ m*-1* and $k_{F}(2D)=(2\pi nd)^{1/2}=9.0\times 10^{9}$ m*-1* in three and two dimensions, respectively. In both cases, the Fermi wavelength $\lambda_{F}=2\pi/k_{F}\ll d$, validating the continuum approximation for states along the perpendicular direction and justifying the anisotropic 3D treatment of the normal state. Then, from the Drude relationship $\rho\approx m/ne^{2}\tau$, we get $\tau\approx 10$ ps, which agrees well with the scattering rates ($\sim\hbar/0.05$ meV = 13 ps) measured by Pan et al. pan using spin- and angle-resolved photoemission spectroscopy. Combining this value of $\tau$ with the Fermi velocity $v_{F}=\hbar k_{F}/m\approx 440$ km/s, we get $l=v_{F}\tau=4.2\ \mu$m. The very long $l$, which well exceeds the superconducting layer thickness $d$, indicates that scattering from the faces that bound the superconducting layer was of a specular nature. This surprising dramatically low scattering indeed supports the possible role of the topological surface states in the formation of the normal state that underlies this exotic interfacial superconducting system.

V Concluding Remarks

Fast pulsed signals of short duration and low duty cycle make it possible to study transport behavior in superconductors at extreme current densities, power densities, and electric fields. In this article, we focused on the use of this technique in the measurement of one of the fundamental critical parameters of the superconducting state, the depairing current $j_{d}$. It was shown how through $j_{d}$ , one can obtain information on various other key parameters of the superconducting state, in particular the penetration depth and consequent superfluid density, which cast light on the normal state. As an example and illustration of this procedure, we described a recent study of the superconducting system formed at the interface between a topological insulator and a chalcogenide. We hope that the information provided here will encourage other groups to utilize this approach.

VI Acknowledgments

The following are acknowledged for useful discussions and other assistance: Charles L. Dean, Manlai Liang, Gabriel F. Saracila, James M. Knight, Luc Fruchter, Ziang Z. Li, Qing Lin He, Hongchao Liu, Jiannong Wang, Rolf Lortz, Iam Keong Sou, Alex Gurevich, Richard A. Webb, Ken Stephenson, and David K. Christen. This work was supported by the U.S. Department of Energy through Grant Number DE-FG02-99ER45763.