BCS Theory and Superconductivity
Total Page:16
File Type:pdf, Size:1020Kb
BCS Theory and Superconductivity 1. Introduction of the superconductor (7) (11). This information based on experiments on naturally occurring mercury and its Superconductivity discovered in 1911 by Onnes (9), is the isotope, observed a decrease in Tc with an increase in quantum phenomena that certain materials exhibit un- isotopic mass. This allows for the assumption that the der particular magnetic and temperature regimes. There basis of superconductivity relies on electron-phonon in- existed no consistent microscopic theory that described teractions, an assumption that would later lead to the why superconductivity arose, from the time it was dis- formation of the BCS theory. covered until the 1950's, only macroscopic theories that allowed you to calculate certain thermodynamic and elec- 2.2 In search of a microscopic theory trodynamic quantities. In 1957, two papers released by Bardeen, Cooper, and Schrieffer (1) (2) described the Because superconductivity was found in materials before conceptual and mathematical foundation for conventional the physics community predicted the phenomena, theo- superconductivity, the Bardeen-Cooper-Schrieffer (BCS) ries were formed to attempt to explain, match and pre- theory, for which they later received the Nobel Prize for dict the characteristics of these materials that undergo in 1972. the phase transition. London and London were inter- We will look at features of superconductors before the ested in superconductivity but did not attempt to de- discovery of the BCS theory, and examine the assump- scribe the reason for it. Instead, they derived an equa- tions and methods used to develop the theory. We will tion for the penetration depth, λ, of the superconductor then calculate and study interesting quantities of the su- (4), but their results consistently overestimated the ex- perconducting system, and finally describe how the re- perimentally found values, and so their assumptions were sults predicted by the BCS theory fare against experi- discarded. mental evidence obtained about superconductors. The Ginzburg-Landau theory in 1950 (6), was a phe- nomenological theory using physical intuition and the 2. Before BCS Theory variational principle of quantum mechanics. It allowed the calculation of macroscopic quantities of the material 2.1 Aspects of Superconductivity in the superconducting state if one assumed the phase transition to be of second order. His results were able Onness discovery of superconductivity came when he wit- to accurately match the experimental results of the time, nessed a sudden drop in the resistance of solid mercury and were later shown to be a specific form of the BCS at 4.2 K. All superconductors show this drop of resis- theory. While useful and accurate for macroscopic quan- tance, either gradually or suddenly, at a particular tran- tities, like the London-London attempt, it did not explain sition temperature, Tc. Infinite conductivity implies that the foundation for superconductivity in these materials. if a current were passed through the material during its superconducting phase, the current would follow forever 3. Foundations of BCS Theory without any dissipation. Another characteristic, later found by Meissner and In this section, we lay out the theoretical grounds for BCS Ochsenfeld (8), is that all superconductors are diamag- theory. All derivations may be referenced in \Theory of nets. Diamagnetism occurs when an external magnetic Superconductivity" by Bardeen, Cooper, and Schrieffer field penetrates only a finite, amount of the material, and or Tinkham's Introduction to Superconductivity. does not hinder the remaining inner parts of the material. The penetration depth is usually small compared to the 3.1 Cooper Pairs width of the material. Also known as the Meissner effect, this event also indicates that a particular magnetic field The BCS theory relies on the assumption that supercon- would destroy the superconductivity of a material. ductivity arises when the attractive Cooper pair interac- The isotope effect discovered by Maxwell and Reynolds tion dominates over the repulsive Coulomb force (2). A describes the relation between Tc and the isotopic mass Cooper pair is a weak electron-electron bound pair medi- 1 ated by a phonon interaction. Although somewhat am- for N0V 1. Thus, the energy of the pair satisfies biguous, one can visualize this pairing by the following E = 2 − 2 ! e−2=N0V < 2 : (1) explanation. Imagine an electron moving within a mate- F ~ c F rial. The Coulomb attraction between the electron and 3.2 The model the positively charged cores of ions in the material will leave a net positive charge in the vicinity. A \paired" We now proceed to write down the model Hamiltonian electron is one with opposite momentum and spin that is for the theory. This is most easily done in the language y attracted to this force. of second quantization. Let ckσ and ckσ to be electron This heuristic explanation is somewhat incomplete, annihilation and creation operators of momentum k and because at the heart of the phonon-mediated interaction spin σ =" or #. The usual commutation relations are: is a long range attraction and thus, requires quantum y (3) 0 fc ; c g = δ (k − k )δ 0 mechanics for a full explanation. Cooper's monumental kσ k0σ σσ 1956 work showed that due to the fermi statistics of the and − − y y electron, this paired e − e state can have energy less fckσ; ck0σg = 0 = fckσ; ck0σg: than that of the Fermi-energy of the material(3). Thus, The proposed Hamiltonian is taken to be at adequately low temperatures, when thermal energy is − − 0 X y X y y not a factor, bound e − e states can form. H = kckσckσ + Vklck"c−k#c−l#cl": (2) We give a short, simplified argument of for this fact. k,σ k;l Suppose we have two electrons interacting with this at- The first term is the usual kinetic energy of the electrons. tractive Cooper force with a background Fermi sea at The second term is the translation of the phonon me- T = 0 by which these electrons only interact via Pauli- diated electron-electron interaction into this framework. exclusion. We look for a zero-momentum wave-function The matrix element Vkk0 may be taken to be general, but of the form: we shall simplify it by using the mean field approximation X mentioned earlier in later calculations. Ψ(r ; r ) = g eik·(r1−r2) (j"#i − j#"i) : 0 1 k Now, in a normal state we would expect no formation k of Cooper pairs, hence the operator ck"c−k# should aver- Antisymmetry demands gk = g−k. Placing this in the age out to zero. It is then natural to define the quantity Schr¨odingerequation HΨ = EΨ yields the following con- b = hc c i : dition: k k" −k# X (E − 2k)gk = Vkk0 gk0 The so-called gap energy is then defined to be 0 k >kF X ∆k = − Vkk0 bk0 : (3) 1 0 R 3 i(k−k )·r 0 where Vkk0 = vol vol d r V (r)e . Now the follow- k ing mean field approximation is made: To allow the exchange of particles it makes sense to con- sider the new Hamiltonian H = H0 + µN where µ is the −V for F < k < F + ~!c V 0 = chemical potential. kk 0 else The Hamiltonian can be diagonalized following the method of Bogoliubov. We defined the linearly trans- where F is the Fermi energy and !c is a cutoff frequency. This indicates that we only consider interactions that are formed states γk0 and γk1 by allowed by the metal's frequency range, similar to the ∗ y ck" = ukγk0 + vkγ assumptions made in the Debeye model. We then have k1 y ∗ y c = −v γk0 + ukγ 1 1 k# k k1 X 2 2 = such that jukj + jvkj = 1. In practice, we can let one of V 2k − E k>k F uk or vk be real. To complete the diagonaization we look Z F +~!c d at the the Hamiltonian in this basis: ! N0 2 − E 0 X n 2 2 y y F H = ξk (jukj − jvkj )(γ γk1 + γ γ−k0) k1 −k0 N0 2F − E + 2~!c k = ln 2 2F − E 2 ∗ ∗ y y +2jvkj + 2ukvkγ−k0γk1 + 2ukvkγk1γ−k0 Simplifying: ∗ ∗ ∗ ∗ y y (∆kukvk + ∆kukvk)(γk1γk1 + γ−k0γ−k0 − 1) 1 1 ∗ ∗ ∗ (2 − E) = ≈ e−2=N0V +(∆k(vk) − ∆k(uk) )γ−k0γk1 F 2=N V 2~!c e 0 − 1 ∗ 2 2 y y ∗ +(∆kvk − ∆kuk)γk1γ−k0 +∆kbk)g : 2 This is analogous to the raising and lowering operators acting on the ground state of the one-dimensional simple harmonic oscillator. 4. Thermodynamic Calculations 4.1 Evaluating Tc By definition of ∆k and γkσ, we have X ∆k = − Vkk0 bk0 k0 X ∗ D y y E = − Vkk0 ukvk0 1 − γk00γk00 − γk01γk01 : k0 Figure 1: The energy gap seen between the superconduct- y 0 y Now, γ γ = δ(k − k ) = γ γ 0 thus, ing state (top) and the normal state (bottom) produces k0 k0 k01 k 1 the order parameter of the system. X ∗ D y y E ∆k = − Vkk0 ukvk0 1 − γk00γk00 − γk01γk01 k0 X ∗ We determine the u and v terms by demanding the coef- = − Vkk0 uk0 vk0 (1 − 2(f(Ek0 )) ficients of γ γ and γy γy . This is equivalent to k0 −k0 k1 k1 −k0 X ∆k0 βEk ∆∗ v q = − Vkk0 tanh k k 2 2 2Ek0 2 = ξk + j∆kj − ξk ≡ Ek − ξk k0 uk V X ∆k0 βEk = tanh : where 2 E 0 2 q 0 k 2 2 k Ek ≡ ξ + j∆kj : k In the last step, we have made the usual mean field ap- From this we observe that the gap energy ∆k is the proximation Vkk0 = −V . In this approximation we know order parameter for this interacting theory as seen in Fig.