Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Dresner, Stability of superconductors.2002

.pdf
Скачиваний:
24
Добавлен:
27.08.2013
Размер:
1.68 Mб
Скачать

20

CHAPTER 2

low temperatures than was observed. This defect was remedied by improvements to Einstein’s theory made in 1912 by Debye. Einstein’s and Debye’s theories are a well-known part of modern physics and will not be described in detail; only Debye’s result will be quoted here. The interested reader can find the details of the Einstein and Debye theories in many places; two good references are Born (1946) and Richtmyer (1947).

According to Debye, the molar specific heat Cv of solids is given by

Cv /rR = 12[(3/ z³) {t³/(et – 1)} dt] – 9z /(ez – 1); z = θ/T

(2.1.1)

where r is the number of atoms per molecule, R is the universal gas constant (8.317 J mol-1 K-1), T is the absolute temperature, and q is the so-called Debye temperature, an empirical parameter related to the phonon spectrum of the solid. (Phonons are the quantized vibrations of the solid lattice.) The integral in brackets has been tabulated by Abramowitz and Stegun (1965), among others.

Figure 2.1. A log-log plot of Cv/r versus thedimensionless temperature T/θ according to Debye’s theory. See Eq. (2.1.2) and the text after Eq. (2.1.4) for the definition of b. (Redrawn from an original appearing in Dresner (1993, “Stability”) with permission of Butterworth-Heineman, Oxford, England.)

Material Properties

21

Figure 2.1 shows a log-log plot of Cv /r calculated using Eq.

(2.1.1) versus

x = 1/ z = T/θ. Two straight-line asymptotes are evident. When

T << q, Cv is

asymptotic to

 

Cv1 = (12π4/5)rRx3 = 1944rx3 J mol¯1 1

(2.1.2)

When T << q, Cv is asymptotic to

 

Cv2 = 3rR = 24.95r J mol¯1 1

(2.1.3)

A family of curves that spans these two asymptotic limits is

 

Cv=

(2.1.4)

 

The value n = 0.85 gives an excellent fit to the exact curve and is a useful approximation for numerical calculation. Eq. (2.1.2), known as the Debye T3-law, is often written Cv1 = βT3, where b = (12π4/5)rR/θ3 = 1944r/θ3 J mol¯1 4.

The only physical property of the solid that enters the Debye theory is the Debye temperature q. Table 2.1 shows its value for several materials of interest in applied superconductivity.

Debye’s theory accounts for the contribution to the specific heat of the vibrations of the atomic lattice of which the solid is composed. (These vibrations are known conventionally as phonons, and Debye’s theory is said to give the phonon contribution to the specific heat.) If the solid is a metal, as are all the entries in Table 2.1, there is an additional contribution to the specific heat from the electrons. This additive contribution, which plays a significant role only at low temperatures, has

the form

 

Cve = γ T

(2.1.5)

where g is the so-called Sommerfeld constant. Its value is given in Table 2.2 for the same materials that appear in Table 2.1.

Table 2.1. The Debye Temperature qfor

 

Various Materials

Material

Debye temperature(K)

 

 

Ag

226

Cu

344

Al

430

Fe

472

 

 

22

 

 

CHAPTER 2

 

Table2.2. The Sommerfeld Constant g for

 

 

Various Materials

 

 

 

 

 

Material

Sommerfeld constant (J mol¯ 1 2)

 

Ag

6.09 x 10¯4

 

Cu

7.44 x 10¯4

 

Al

1.36 x 10¯3

 

Fe

5.02x 10¯3

 

2.2. SPECIFIC HEATS OF TYPE-II SUPERCONDUCTORS

The specific heat of a type-II superconductor in the normal state is also described by the Debye-Sommerfeld theory given in the previous section. In the superconducting state, it is given approximately by

Cs(T,H) = Cn(T) + 3γ T3 /Tco2 γ T+ γ TH/Hc2(0)

(2.2.1)

= (β + 3γ /Tco2 )T 3+ γ TH/Hc2(0)

(2.2.2)

where Cs(T,H) is the specific heat in the superconducting state at temperature T and applied magnetic field H, Cn(T) is the specific heat in the normal state at temperature T, Tco is the critical temperature at zero applied field, and Hc2(0) is the (upper) critical field at zero temperature. The second line applies at low temperatures, when

the Debye T 3-law, Cphonon = βT 3, is a satisfactory approximation. At the critical temperature Tc (at applied field H), where the phase change occurs, the specific heat

undergoes a jump

Cs(Tc,H) – Cn(Tc) = γ Tc[3Tc2 /T co2 – 1 + H/Hc 2(0)]

(2.2.3)

When there is no applied field,

 

Cs(Tco,H) – Cn(Tco) = 2γ Tco

(2.2.4)

Such jumps in the specific heat (so-called l-anomalies) are characteristic of second-order phase changes in which there is no latent heat.

Figure 2.2 shows a plot of measured values of C/T versus T 2 for the low-tem- perature superconductor Nb 44.6-wt-% Ti reported by Elrod, Miller, and Dresner (1982). The lines are best fits to Eq. (2.2.2) using the values γ = 0.145 J kg¯1 2 , β = 2.3 x 10¯3 J kg¯1 4 , and µ oHc2(0) = 14 T. The fit is excellent and reinforces our belief that Eq. (2.2.2) is a description of superconductor specific heats at low temperatures adequate for use in the theory of superconductor stability. It is perhaps worth noting here that our goal in this chapter is to find, where possible, simple, widely applicable formulas that allow approximate computation of the physical properties we use to study the stability of superconductors.

The Debye theory is based on an approximate phonon spectrum in which the density of phonon states varies as the square of the phonon frequency. At low

Material Properties

23

Figure 2.2. Measured values of C/T versus T2 for Nb 44.6-wt-% Ti (Elrod et al., 1982). (Redrawn from an original appearing in Elrod et al. (1982) with permission of Plenum Publishing Corp., New York.)

temperatures, only low-frequency phonons are excited, and at this end of the phonon spectrum, the quadratic approximation is accurate. At higher temperature, such as might typify the new high-temperature superconductors, the fit to Eq. (2.2.1) might not be as good. Indeed, deficiencies in the Debye theory due to its simplified phonon spectrum have been known for some time—long ago, Born and von Karman (Reed and Clark, 1983) sought to improve the Debye theory by invoking a more realistic phonon spectrum, but only at the cost of increased computational complexity. Occasionally, the difficulties with the Debye theory have been circumvented in an ad hoc way by letting the Debye temperature q vary with T, but such a “cure” is only applicable when the specific heat is already known experimentally over a wide range of temperatures.

This is exactly the situation with recent measurements of the specific heat of YBa2Cu3O7 recently published by Pavese and Malishev (1994): see Fig. 2.3. Shown in Fig. 2.4 are plots of the specific heat calculated using Eqs. (2.1.1) and (2.2.1) for five values of the Debye temperature together with seven points taken from Fig. 2.3. At low temperatures q ~ 375 K, whereas at high temperatures q ~ 525 K. Choosing the smaller value leads to a maximum error in the neighborhood of 100 K of about 30%.

We shall not abandon Eqs. (2.1.1) and (2.2.1) for high-temperature superconductors, but rather accept the facts of the last paragraph as a caution on the limitation of their accuracy.

The next four optional sections present a discussion of the thermodynamics of superconductors that culminates in the derivation of Eq. (2.2.1).

24

CHAPTER 2

Figure 2.3. Measured values of C/T versus T for YBa2Cu3O7 (Pavese and Malishev, 1994). (Redrawn from an original appearing in Pavese and Malishev (1994) with permission of Plenum Publishing Corp., New York.)

Figure 2.4. Comparison of the measurements of (Pavese and Malishev, 1994) for YBa2Cu3O7 with the Debye theory. g = 0.04 J/kg·K2, r = 13, MCLWT = 0.6662 kg/mol, Tco = 93 K.

Material Properties

25

2.3. FIRST LAW OF THERMODYNAMICS FOR A MAGNETIZABLE BODY

Presented below is a phenomenological theory of the specific heat of the superconducting state. In this theory, the specific heat is deduced by thermodynamic arguments from the measured value of another, more accessible quantity, the magnetization M of the superconductor in the presence of an applied field

H.

The state of a magnetizable body is described by Maxwell’s equations

 

(2.3.1)

x H = J

(2.3.2)

B = µ o (H+ M)

(2.3.3)

where E is the electric field, B is the magnetic induction, H is the applied magnetic field, created by J, the current density produced by the external power supply, and

M is the magnetization of the body. Because measurements of

the magnetization

are carried out slowly, the displacement current density

has been omitted

from the right-hand side of Eq. (2.3.2).

In MKS units, M has the dimensions of A/m = Am2/m3 or the dimensions of dipole moment per unit volume. In a conceptual sense, M can be thought of as an increase in the total field inside the body brought about when the applied field orients the body’s atomic dipoles.

To write the first law of thermodynamics for a magnetizable body, we need to know the increment of work the body does when its magnetic state changes. Change in the magnetic state, i.e., change in H and M, causes change in the induction B, which in turn creates the electric field E. This field applies a force Ee to the charge carriers produced by the external power supply that create the external field H. This force (the back emf) is opposed by an equal and opposite force –Ee produced by the external power supply. If the velocity of these charge carriers is v, eE·v is the rate at which the external power supply does work on them. If their density is n, neE·v = –E·J is the rate at which the external power supply does work on the charge carriers in a unit volume. The quantity d(vol) is thus the rate at which the external power supply does work when the state (H,M) is changing.

To relate the quantity E·J to the state variables H and M, we scalar multiply

Eq. (2.3.1) by H and Eq. (2.3.2) by E and subtract:

 

H.xEE. xH+H.B/t =-E.J

(2.3.4)

26

CHAPTER 2

Now the first two terms in Eq. (2.3.4) together equal

(E x H), so that if we

integrate Eq. (2.3.4) over a large volume V whose surface S surrounds the body and the external power supply, we find

(2.3.5)

As the volume V grows large, the magnitude of the product E x H falls faster than the surface S increases, and the surface integral vanishes in the limit. Thus

(2.3.6)

(2.3.7)

The passage from Eq. (2.3.6) to Eq. (2.3.7) has been achieved by substituting the right-hand side of Eq. (2.3.3) for B. The second volume integral only extends over the volume of the body, for only there is M different from zero.

The left-hand side of Eq. (2.3.6) is the rate at which the external power source is supplying power. The first term on the right-hand side of Eq. (2.3.7) is the rate of increase of the energy of the magnetic field and the second term is the rate of increase of the energy of the body, i.e., the rate as which the external power supply does work on the body. For thermodynamic purposes, we need the rate at which work is done by the body during a change of state; it is

(2.3.8)

In the simplest experiment to measureM, the external field H is applied parallel to the axis of a long wire so that both H and M are parallel and uniform over the body. Then the rate of work done by the body is

oH d(MV )/dt

(2.3.9)

where V is the volume of the body. Thus, finally, the work done by the body in a change of state is

(2.3.10)

Material Properties

27

2.4. GlBBS FREE ENERGY OF A MAGNETIZABLE BODY

The two laws of thermodynamics may be written together for the magnetizable body as

(2.4.1)

where T is the temperature, s the entropy of the body per unit mass, u its internal energy per unit mass, and δ its density. The most convenient independent variables to use in the thermodynamic discussion are the applied field H and the temperature T since they are under our control. Accordingly, we introduce the function

(2.4.2)

which is the Gibbs free energy per unit mass of the body. Then

(2.4.3)

so that

(2.4.4)

Finally, the specific heat at constant applied field CH is given by

(2.4.5)

According to Eq. (2.4.4),

(2.4.6)

where the subscript s on g indicates the superconducting state. When the pair of values (T,H) lies on the phase boundary between the mixed and the normal states, the Gibbs free energy of the normal and superconducting states are equal (remember, the transition is second-order). Thus

(2.4.7)

where Hc2 is the upper critical field at temperature T (see Fig. 1.2). If we subtract Eq. (2.4.7) from Eq. (2.4.6), we find

28

CHAPTER 2

(2.4.8)

It is important to note here that because superconductors are diamagnetic (i.e., exclude the magnetic field), M is oppositely directed to H. Thus gs(T,H) < gn(T) when H < Hc2 as required.

2.5. SPEClFlC HEAT AT ZERO FlELD

At zero applied field, Eq. (2.4.8) takes the form

(2.5.1)

Now the measured curves of M versus H look like the sketch in Fig. 2.5. Each curve corresponds to a particular fixed temperature. As H increases from zero, M starts out equal to –H : this is the Meissner state, which persists until H = Hc1. Then M begins slowly to rise until it reaches zero at H = Hc2. The curves at different temperatures resemble one another; we make the assumption now that they are strictly similar geometrically to one another. Then since both the height Hc1 and

the base Hc2 of the M-H curve scale as (1 – T 2/T co2 ), where Tco is the critical temperature at zero field,

(2.5.2)

Figure 2.5. A sketch of the magnetization M versus the applied magnetic field H. Each curve corresponds to a particular fixed temperature. In the theory it is assumed that these curves are strictly similar geometrically to one another, and the sketch is drawn consistently with this assumption.

Material Properties

29

For convenience, we denote the integral on the right-hand side of Eq. (2.5.2) as –a, where a > 0. Then

gs(T,0) = gn(T) – (µ0 /δ)a(1 – T 2 /T 2co )2

(2.5.3)

so that the specific heats in the normal and superconducting states are related by

(2.5.4)

Now it is experimentally known that Cs(T,0) approaches zero with decreasing temperature faster than linearly, so the linear term -4µo /d)aT/T co2 must cancel the

Sommerfeld term γ T in Cn(T). Therefore, 4(µo /δ)a/T 2co = γ, and Eq. (2.5.4) becomes

Cs(T,0) = Cn(T) + 3γ T 3/T 2co γ T= (β + 3γ /T 2co )T 3

(2.5.5)

the last equality applying when Cn(T) can be represented as γ T+βT 3 (i.e., at low temperatures).

2.6.CONTRIBUTION OF THE MAGNETIC FIELD TO THE SPEClFlC HEAT

To determine the effect of the magnetic field on the specific heat we must calculate

(2.6.1)

Now we make the further assumption that the shallow curves in the M-H diagram of Fig. 2.5 are in fact parallel straight lines. The slope of these parallel lines is Hc1(0)/[Hc2(0)–Hc1(0)] ~ Hc1(0) /Hc2(0). Thus,

M =

[Hc1(0)/Hc2(0)][H Hc2(T)]

 

=

[Hc1(0)/Hc2(0)][HHc2(0)(1–T 2/T co2 )]

(2.6.2)

Furthermore,

(2.6.3)

so that