Ginzburg–Landau theory

From testwiki
Jump to navigation Jump to search

Template:Short description Template:Distinguish Template:For

In physics, Ginzburg–Landau theory, often called Landau–Ginzburg theory, named after Vitaly Ginzburg and Lev Landau, is a mathematical physical theory used to describe superconductivity. In its initial form, it was postulated as a phenomenological model which could describe type-I superconductors without examining their microscopic properties. Alexei Abrikosov later used the theory to describe type-II superconductors as well.[1] One GL-type superconductor is YBCO, and generally all cuprates.[2]

Later, a version of Ginzburg–Landau theory was derived from the Bardeen–Cooper–Schrieffer microscopic theory by Lev Gor'kov,[3] thus showing that it also appears in a limit of the microscopic theory and giving it a microscopic interpretation of all its parameters.Template:Not in body The theory can also be given a general geometric setting, placing it in the context of Riemannian geometry, where in many cases exact solutions can be given. This general setting then extends to quantum field theory and string theory, again owing to its solvability, and its close relation to other, similar systems.

Ginzburg was awarded a third of the 2003 Nobel prize in physics for his development of the theory. Abrikosov was awarded another third of this prize for his description of type-II superconductors using the theory.[4]

Introduction

Free energy

Based on Landau's previously established theory of second-order phase transitions Ginzburg and Landau argued that the free energy density fs of a superconductor near the superconducting transition can be expressed in terms of a complex order parameter field ψ(r)=|ψ(r)|eiϕ(r), where the quantity |ψ(r)|2 is a measure of the local density of superconducting electrons ns(r) analogous to a quantum mechanical wave function.[3] In the convention used in this article, the relation is |ψ(r)|2=ns(r)/2, which is the density of Cooper pairs.[5]Template:Rp While ψ(r) is nonzero below a phase transition into a superconducting state, no direct interpretation of this parameter was given in the original paper.Template:Cn Assuming smallness of |ψ| and smallness of its gradients, the free energy density has the form of a field theory and exhibits U(1) gauge symmetry. In Gaussian units, it is written as

fs=fn+α(T)|ψ|2+12β(T)|ψ|4+14m|(−iℏ∇−2ec𝐀)ψ|2+𝐇28π,

where

  • fn is the free energy density of the normal phase,
  • α(T) and β(T) are phenomenological parameters that are functions of T (and often written just α and β).
  • 2m the mass of a Cooper pair,
  • 2e is the charge of a Cooper pair,
  • 𝐀 is the magnetic vector potential, and
  • 𝐇=∇×𝐀 is the magnetic field.[5]Template:Rp

The total free energy is given by F=∫fsd3r.

Ginzburg–Landau equations

By minimizing F with respect to variations in the order parameter ψ and the vector potential 𝐀, one arrives at the Ginzburg–Landau equations:[5]Template:Rp

αψ+β|ψ|2ψ+14m(−iℏ∇−2ec𝐀)2ψ=0

∇×𝐇=4πc𝐉,𝐉=−ieℏ2m(ψ∗∇ψ−ψ∇ψ∗)−2e2mc|ψ|2𝐀,

where 𝐉 denotes the dissipation-free electric current density, or supercurrent density. The first equation, which bears some similarities to the time-independent Schrödinger equation but is principally different due to a nonlinear term, determines the order parameter, ψ. The second equation is instead Maxwell-like and provides the superconducting current and the magnetic field.

Additionally, the boundary condition at the surface of the superconductor can be determined to be[5]Template:Rp

𝐧(−iℏ∇−2ec𝐀)ψ|b=0,

where 𝐧 is the vector normal to the surface and b denotes the location of the boundary.

Close to the superconducting transition, it is assumed that β(T) is a positive constant and that α(T)=aT−TCTC, where a is a positive constant and TC is the critical temperature.[5]Template:Rp

Simple interpretation

Visualization of the free energy density of the Ginzburg–Landau theory as a function of a real order parameter above the critical temperature (left) and below the critical temperature (right).

Consider a homogeneous superconductor where ψ(𝐫) is constant and there is no external magnetic field. The equation for the free energy then simplifies to[5]Template:Rp fs=fn+α(T)|ψ|2+12β|ψ|4. The equilibrium order parameter is found by minimizing this equation with respect to ψ.

Above the critical temperature, α(T) is positive and the only minimum is given by ψ=0. This corresponds to the absence of superconductivity.[5]Template:Rp

Below the superconducting transition, α(T) is negative, and the minima are given by[5]Template:Rp |ψ0|2=−αβ. The density of electrons participating in the supercurrent ns is given by |ψ|2=12ns. Thus, this corresponds to the superconducting phase.[5]Template:Rp

Characteristic quantities

Coherence length and penetration depth

The Ginzburg–Landau equations predicts two characteristic lengths in a superconductor. The first characteristic length is termed coherence length. In one dimension at zero magnetic field, the first GL equation becomes[6]Template:Rp αψ+β|ψ|2ψ−ℏ24md2ψdx2=0. At an interface between a superconductor and a normal conductor, this equation is solved by ψ(x)=αβtanh⁡(xξ2), where ξ≡ℏ24m|α| is the coherence length.[6]Template:Rp This length scale also sets the exponential law according to which small perturbations of density of superconducting electrons recover their equilibrium value ψ0.Template:Cn

The magnetic field strength as a function of position at the boundary between a normal conductor and a superconductor given an external magnetic field H0.

The second length scale is the penetration depth. Consider an external magnetic field H0 applied to a semi-infinite superconducting slab. Assuming that ψ=ψ0 is constant on the length scale we are considering,Template:Efn the second GL equation becomes d2Aydx2=8πe2mc2|ψ|2Ay. This is equation is solved for the magnetic field by H=H0e−x/λ, where λ≡mc28πe2|ψ|2=mc28πe2βα is the penetration depth. By identifying the order parameter with the density of superconducting electrons, it becomes apparent that this is equivalent to the penetration depth in the London theory.[5]Template:Rp

The ratio between these two length scales, κ≡λξ=βc2m22πe2ℏ2 is called the Ginzburg–Landau parameter.[5]Template:Rp The importance of this property is that it distinguishes between type I and type II superconductors; superconductors for which κ<1/2 are of type I, and those for which κ>1/2 are of type II.[6]Template:Rp

Critical magnetic field

To calculate the magnetic field Hcm above which superconductivity is destroyed in the bulk, ψ is considered to be constant and in equilibrium. Using these assumptions, the free energy is fs−fn=α|ψ0|2+12β|ψ0|4+𝐇28π=−α22β+𝐇28π. At the phase transition, the free energies of the normal and superconducting phases should be equal. This gives[5]Template:RpHcm2=4πα2β.

Dimensionless form of the Ginzburg–Landau equations

In Gaussian units, the GL equations read αψ+β|ψ|2ψ+14m(−iℏ∇−2ec𝐀)2ψ=0

∇×𝐇=4πc𝐉,𝐉=−ieℏ2m(ψ∗∇ψ−ψ∇ψ∗)−2e2mc|ψ|2𝐀, with the boundary condition 𝐧(−iℏ∇−2ec𝐀)ψ|b=0.

By changing to a unit system which is defined by the characteristic quantities, the equations can be made dimensionless. For this, the units become[5]Template:Rp ψ′=ψψ0,𝐫′=𝐫λ,𝐇′=𝐇Hcm2.

Then, the GL equations become dimensionless:[5]Template:Rp (−i∇κ−𝐀)2ψ−ψ+|ψ|2ψ=0, ∇×(∇×𝐀)=−i2κ(ψ∗∇ψ−ψ∇ψ∗)−|ψ|2𝐀, with the boundary condition 𝐧(−i∇κ−𝐀)ψ|b=0.

In this unit system, the equations then depend only on the Ginzburg–Landau parameter κ.[5]Template:Rp

Applications

Thin films

One of the systems that can be analysed with the GL theory is a thin film. Consider a type-I (κ≪1) superconducting film extending infinitely far in the y- and z-directions, with a thickness d≪λ. In addition, an external magnetic field H0 is applied in the z-direction.[5]Template:Rp

By symmetry, only derivatives to x can be nonzero. This means that the gauge of the vector potential can be chosen such that[5]Template:Rp 𝐀=A(x)ŷ,𝐇(x)=dA(x)dxẑ. Using this gauge, and choosing ψ to be real, the first dimensionless GL equations becomes[5]Template:Rp d2ψdx2=κ2ψ(A2−1+ψ2).

Because of the assumption that the superconductor is of type-I, it is reasonable to assume that ψ is approximately constant over the thickness of the thin film. The second GL equation then becomes[5]Template:Rp d2Adx2=ψ2A.

By using that the magnetic field must be continuous at the boundary, i.e. H(±d/2)=H0, the solution to for the magnetic field is[5]Template:Rp H=H0cosh⁡(ψx)cosh⁡(ψd/2).

In the second-order approximation of the first GL equation, the relation between ψ and H0 can be found to be[5]Template:Rp H02=2ψ2(1−ψ2)cosh2(12dψ)[sinh⁡(dψ)/(dψ)]−1≈12(1−ψ2)d2.

In this system, the difference in the free energies of the superconducting and normal state is[5]Template:Rp ∫−12d12dfs−fndx=Hcm24π∫−12d12d[−ψ2+12ψ4+(H−H0)2+A2ψ2]dx.

Using the expression for H0 and that at the critical magnetic field, this has to be zero. Therefore,[5]Template:Rp Hc2=ψ2(2−ψ2)2(1−tanh⁡(12dψ)12dψ)≈6(2−ψ2)d2.

At the transition point, |ψ|2=0, so[5]Template:Rp Hc=26d=26λdHcm.

Therefore, the critical magnetic field is higher for a thin film than for a bulk superconductor. Moreover, it can be derived that for a thin film, the phase transition between the normal and superconducting phase is of second order, while it is of first order in a bulk material.[5]Template:Rp This change in the order of the phase transition happens for d=5λ.[5]Template:Rp

Flux quantization

The order parameter can be written in terms of the modulus and the phase as ψ=|ψ|eiϕ. The supercurrent is then[5]Template:Rp 𝐣=ℏem|ψ|2(∇ϕ−2eℏc𝐀).

The setup for the derivation of magnetic flux quantization, showing the external magnetic field inside the hollow superconductor and the integration contour 𝒞.

Now consider a hollow superconducting cylinder that encloses a non-superconducting material. In the non-superconducting region, and external magnetic field 𝐇 is applied. Integrating the supercurrent over a contour inside the hollow cylinder that encloses the non-superconducting region, one gets[5]Template:Rp ∮𝒞𝐣d𝐥=|ψ|2ℏem(∮𝒞∇ϕd𝐥−2eℏc∮𝒞𝐀d𝐥)

Inside the wall, there is no current and the first integral vanishes. From electromagnetism, it is known that ∮𝒞𝐀d𝐥=∫𝒞𝐇d𝐒=Φ, where Φ is the magnetic flux. Finally, the first integral on the left side is constrained by the fact that ψ must be a single-valued function. Because eix=eix+2πn for any integer n,[5]Template:Rp ∮𝒞∇ϕd𝐥=2πn.

Putting it together, these constraints imply that 2eℏcΦ≡2πΦΦ0=2πn.

Thus, the flux inside a hollow cylinder can only be a multiple of the magnetic flux quantum Φ0=πℏc/e. This phenomenon is known as flux quantization.[5]Template:Rp

Fluctuations

The phase transition from the normal state is of second order for Type II superconductors, taking into account fluctuations, as demonstrated by Dasgupta and Halperin, while for Type I superconductors it is of first order, as demonstrated by Halperin, Lubensky and Ma.[7]

Classification of superconductors

In the original paper Ginzburg and Landau observed the existence of two types of superconductors depending on the energy of the interface between the normal and superconducting states. The Meissner state breaks down when the applied magnetic field is too large. Superconductors can be divided into two classes according to how this breakdown occurs. In Type I superconductors, superconductivity is abruptly destroyed when the strength of the applied field rises above a critical value Hc. Depending on the geometry of the sample, one may obtain an intermediate state[8] consisting of a pattern of regions of normal material carrying a magnetic field mixed with regions of superconducting material containing no field. In Type II superconductors, raising the applied field past a critical value Hc1 leads to a mixed state (also known as the vortex state) in which an increasing amount of magnetic flux penetrates the material, but there remains no resistance to the flow of electric current as long as the current is not too large. At a second critical field strength Hc2, superconductivity is destroyed. The mixed state is actually caused by vortices in the electronic superfluid, sometimes called fluxons because the flux carried by these vortices is quantized. Most pure elemental superconductors, except niobium and carbon nanotubes, are Type I, while almost all impure and compound superconductors are Type II.

The most important finding from Ginzburg–Landau theory was made by Alexei Abrikosov in 1957. He used Ginzburg–Landau theory to explain experiments on superconducting alloys and thin films. He found that in a type-II superconductor in a high magnetic field, the field penetrates in a triangular lattice of quantized tubes of flux vortices.[9]

Geometric formulation

The Ginzburg–Landau functional can be formulated in the general setting of a complex vector bundle over a compact Riemannian manifold.[10] This is the same functional as given above, transposed to the notation commonly used in Riemannian geometry. In multiple interesting cases, it can be shown to exhibit the same phenomena as the above, including Abrikosov vortices (see discussion below).

For a complex vector bundle E over a Riemannian manifold M with fiber ℂn, the order parameter ψ is understood as a section of the vector bundle E. The Ginzburg–Landau functional is then a Lagrangian for that section:

ℒ(ψ,A)=∫M|g|dx1∧⋯∧dxm[|F|2+|Dψ|2+14(σ−|ψ|2)2]

The notation used here is as follows. The fibers ℂn are assumed to be equipped with a Hermitian inner product ⟨⋅,⋅⟩ so that the square of the norm is written as |ψ|2=⟨ψ,ψ⟩. The phenomenological parameters α and β have been absorbed so that the potential energy term is a quartic mexican hat potential; i.e., exhibiting spontaneous symmetry breaking, with a minimum at some real value σ∈ℝ. The integral is explicitly over the volume form

∗(1)=|g|dx1∧⋯∧dxm

for an m-dimensional manifold M with determinant |g| of the metric tensor g.

The D=d+A is the connection one-form and F is the corresponding curvature 2-form (this is not the same as the free energy F given up top; here, F corresponds to the electromagnetic field strength tensor). The A corresponds to the vector potential, but is in general non-Abelian when n>1, and is normalized differently. In physics, one conventionally writes the connection as d−ieA for the electric charge e and vector potential A; in Riemannian geometry, it is more convenient to drop the e (and all other physical units) and take A=Aμdxμ to be a one-form taking values in the Lie algebra corresponding to the symmetry group of the fiber. Here, the symmetry group is SU(n), as that leaves the inner product ⟨⋅,⋅⟩ invariant; so here, A is a form taking values in the algebra 𝔰𝔲(n).

The curvature F generalizes the electromagnetic field strength to the non-Abelian setting, as the curvature form of an affine connection on a vector bundle . It is conventionally written as

F=D∘D=dA+A∧A=(∂Aν∂xμ+AμAν)dxμ∧dxν=12(∂Aν∂xμ−∂Aμ∂xν+[Aμ,Aν])dxμ∧dxν

That is, each Aμ is an n×n skew-symmetric matrix. (See the article on the metric connection for additional articulation of this specific notation.) To emphasize this, note that the first term of the Ginzburg–Landau functional, involving the field-strength only, is

ℒ(A)=YM(A)=∫M∗(1)|F|2

which is just the Yang–Mills action on a compact Riemannian manifold.

The Euler–Lagrange equations for the Ginzburg–Landau functional are the Yang–Mills equations [11]

D∗Dψ=12(σ−|ψ|2)ψ

and

D∗F=−Re⁡⟨Dψ,ψ⟩

where D∗ is the adjoint of D, analogous to the codifferential δ=d∗. Note that these are closely related to the Yang–Mills–Higgs equations.

Specific results

In string theory, it is conventional to study the Ginzburg–Landau functional for the manifold M being a Riemann surface, and taking n=1; i.e., a line bundle.[12] The phenomenon of Abrikosov vortices persists in these general cases, including M=ℝ2, where one can specify any finite set of points where ψ vanishes, including multiplicity.[13] The proof generalizes to arbitrary Riemann surfaces and to Kähler manifolds.[14][15][16][17] In the limit of weak coupling, it can be shown that |ψ| converges uniformly to 1, while Dψ and dA converge uniformly to zero, and the curvature becomes a sum over delta-function distributions at the vortices.[18] The sum over vortices, with multiplicity, just equals the degree of the line bundle; as a result, one may write a line bundle on a Riemann surface as a flat bundle, with N singular points and a covariantly constant section.

When the manifold is four-dimensional, possessing a spinc structure, then one may write a very similar functional, the Seiberg–Witten functional, which may be analyzed in a similar fashion, and which possesses many similar properties, including self-duality. When such systems are integrable, they are studied as Hitchin systems.

Self-duality

When the manifold M is a Riemann surface M=Σ, the functional can be re-written so as to explicitly show self-duality. One achieves this by writing the exterior derivative as a sum of Dolbeault operators d=∂+∂‾. Likewise, the space Ω1 of one-forms over a Riemann surface decomposes into a space that is holomorphic, and one that is anti-holomorphic: Ω1=Ω1,0⊕Ω0,1, so that forms in Ω1,0 are holomorphic in z and have no dependence on z‾; and vice-versa for Ω0,1. This allows the vector potential to be written as A=A1,0+A0,1 and likewise D=∂A+∂‾A with ∂A=∂+A1,0 and ∂‾A=∂‾+A0,1.

For the case of n=1, where the fiber is ℂ so that the bundle is a line bundle, the field strength can similarly be written as

F=−(∂A∂‾A+∂‾A∂A)

Note that in the sign-convention being used here, both A1,0,A0,1 and F are purely imaginary (viz U(1) is generated by eiθ so derivatives are purely imaginary). The functional then becomes

ℒ(ψ,A)=2πσdeg⁡L+∫Σi2dz∧dz‾[2|∂‾Aψ|2+(∗(−iF)−12(σ−|ψ|2)2]

The integral is understood to be over the volume form

∗(1)=i2dz∧dz‾,

so that

Area⁡Σ=∫Σ∗(1)

is the total area of the surface Σ. The ∗ is the Hodge star, as before. The degree deg⁡L of the line bundle L over the surface Σ is

deg⁡L=c1(L)=12π∫ΣiF

where c1(L)=c1(L)[Σ]∈H2(Σ) is the first Chern class.

The Lagrangian is minimized (stationary) when ψ,A solve the Ginzberg–Landau equations

∂‾Aψ=0∗(iF)=12(σ−|ψ|2)

Note that these are both first-order differential equations, manifestly self-dual. Integrating the second of these, one quickly finds that a non-trivial solution must obey

4πdeg⁡L≤σArea⁡Σ.

Roughly speaking, this can be interpreted as an upper limit to the density of the Abrikosov vortices. One can also show that the solutions are bounded; one must have |ψ|≤σ.

In string theory

In particle physics, any quantum field theory with a unique classical vacuum state and a potential energy with a degenerate critical point is called a Landau–Ginzburg theory. The generalization to N = (2,2) supersymmetric theories in 2 spacetime dimensions was proposed by Cumrun Vafa and Nicholas Warner in November 1988;[19] in this generalization one imposes that the superpotential possess a degenerate critical point. The same month, together with Brian Greene they argued that these theories are related by a renormalization group flow to sigma models on Calabi–Yau manifolds.[20] In his 1993 paper "Phases of N = 2 theories in two-dimensions", Edward Witten argued that Landau–Ginzburg theories and sigma models on Calabi–Yau manifolds are different phases of the same theory.[21] A construction of such a duality was given by relating the Gromov–Witten theory of Calabi–Yau orbifolds to FJRW theory an analogous Landau–Ginzburg "FJRW" theory.[22] Witten's sigma models were later used to describe the low energy dynamics of 4-dimensional gauge theories with monopoles as well as brane constructions.[23]

See also

Template:Div col

Template:Div col end

Notes

Template:Notelist

References

Template:Reflist

Papers

  • V.L. Ginzburg and L.D. Landau, Zh. Eksp. Teor. Fiz. 20, 1064 (1950). English translation in: L. D. Landau, Collected papers (Oxford: Pergamon Press, 1965) p. 546
  • A.A. Abrikosov, Zh. Eksp. Teor. Fiz. 32, 1442 (1957) (English translation: Sov. Phys. JETP 5 1174 (1957)].) Abrikosov's original paper on vortex structure of Type-II superconductors derived as a solution of G–L equations for κ > 1/√2
  • L.P. Gor'kov, Sov. Phys. JETP 36, 1364 (1959)
  • A.A. Abrikosov's 2003 Nobel lecture: pdf file or video
  • V.L. Ginzburg's 2003 Nobel Lecture: pdf file or video

Template:Quantum field theories Template:Superconductivity Template:Authority control