Ginzburg–Landau theory

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. One GL-type superconductor is the famous YBCO, and generally all cuprates.^[1]

Later, a version of Ginzburg–Landau theory was derived from the Bardeen–Cooper–Schrieffer microscopic theory by Lev Gor'kov,^[2] thus showing that it also appears in some limit of microscopic theory and giving microscopic interpretation of all its parameters. 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.

Introduction[edit]

Based on Landau's previously established theory of second-order phase transitions, Ginzburg and Landau argued that the free energy density ${\displaystyle f_{s}}$ of a superconductor near the superconducting transition can be expressed in terms of a complex order parameter field ${\displaystyle \psi (r)=|\psi (r)|e^{i\phi (r)}}$, where the quantity ${\displaystyle |\psi (r)|^{2}}$ is a measure of the local density of superconducting electrons ${\displaystyle n_{s}(r)}$ analogous to a quantum mechanical wave function.^[2] While ${\displaystyle \psi (r)}$ is nonzero below a phase transition into a superconducting state, no direct interpretation of this parameter was given in the original paper. Assuming smallness of ${\displaystyle |\psi |}$ and smallness of its gradients, the free energy density has the form of a field theory and exhibits U(1) gauge symmetry:

$${\displaystyle f_{s}=f_{n}+\alpha (T)|\psi |^{2}+{\frac {1}{2}}\beta (T)|\psi |^{4}+{\frac {1}{2m^{*}}}\left|\left(-i\hbar \nabla -{\frac {e^{*}}{c}}\mathbf {A} \right)\psi \right|^{2}+{\frac {\mathbf {B} ^{2}}{8\pi }},}$$

where

• ${\displaystyle f_{n}}$ is the free energy density of the normal phase,
• ${\displaystyle \alpha (T)}$ and ${\displaystyle \beta (T)}$ are phenomenological parameters that are functions of T (and often written just ${\displaystyle \alpha }$ and ${\displaystyle \beta }$).
• ${\displaystyle m^{*}}$ is an effective mass,
• ${\displaystyle e^{*}}$ is an effective charge (usually ${\displaystyle 2e}$, where e is the charge of an electron),
• ${\displaystyle \mathbf {A} }$ is the magnetic vector potential, and
• ${\displaystyle \mathbf {B} =\nabla \times \mathbf {A} }$ is the magnetic field.

The total free energy is given by ${\displaystyle F=\int f_{s}d^{3}r}$. By minimizing ${\displaystyle F}$ with respect to variations in the order parameter ${\displaystyle \psi }$ and the vector potential ${\displaystyle \mathbf {A} }$, one arrives at the Ginzburg–Landau equations

$${\displaystyle \alpha \psi +\beta |\psi |^{2}\psi +{\frac {1}{2m^{*}}}\left(-i\hbar \nabla -{\frac {e^{*}}{c}}\mathbf {A} \right)^{2}\psi =0}$$

$${\displaystyle \nabla \times \mathbf {B} ={\frac {4\pi }{c}}\mathbf {J} \;\;;\;\;\mathbf {J} ={\frac {e^{*}}{m^{*}}}\operatorname {Re} \left\{\psi ^{*}\left(-i\hbar \nabla -{\frac {e^{*}}{c}}\mathbf {A} \right)\psi \right\},}$$

where ${\displaystyle J}$ denotes the dissipation-free electric current density and Re the real part. 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, ${\displaystyle \psi }$. The second equation then provides the superconducting current.

Simple interpretation[edit]

Consider a homogeneous superconductor where there is no superconducting current and the equation for ψ simplifies to:

$${\displaystyle \alpha \psi +\beta |\psi |^{2}\psi =0.}$$

This equation has a trivial solution: ψ = 0. This corresponds to the normal conducting state, that is for temperatures above the superconducting transition temperature, T > T_c.

Below the superconducting transition temperature, the above equation is expected to have a non-trivial solution (that is ${\displaystyle \psi \neq 0}$). Under this assumption the equation above can be rearranged into:

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

When the right hand side of this equation is positive, there is a nonzero solution for ψ (remember that the magnitude of a complex number can be positive or zero). This can be achieved by assuming the following temperature dependence of ${\displaystyle \alpha :\alpha (T)=\alpha _{0}(T-T_{\rm {c}})}$ with ${\displaystyle \alpha _{0}/\beta >0}$:

• Above the superconducting transition temperature, T > T_c, the expression α(T) / β is positive and the right hand side of the equation above is negative. The magnitude of a complex number must be a non-negative number, so only ψ = 0 solves the Ginzburg–Landau equation.
• Below the superconducting transition temperature, T < T_c, the right hand side of the equation above is positive and there is a non-trivial solution for ψ. Furthermore,
  $${\displaystyle |\psi |^{2}=-{\frac {\alpha _{0}(T-T_{c})}{\beta }},}$$

  
  that is ψ approaches zero as T gets closer to T_c from below. Such a behavior is typical for a second order phase transition.

In Ginzburg–Landau theory the electrons that contribute to superconductivity were proposed to form a superfluid.^[3] In this interpretation, |ψ|² indicates the fraction of electrons that have condensed into a superfluid.^[3]

Coherence length and penetration depth[edit]

The Ginzburg–Landau equations predicted two new characteristic lengths in a superconductor. The first characteristic length was termed coherence length, ξ. For T > T_c (normal phase), it is given by

${\displaystyle \xi ={\sqrt {\frac {\hbar ^{2}}{2m^{*}|\alpha |}}}.}$

while for T < T_c (superconducting phase), where it is more relevant, it is given by

${\displaystyle \xi ={\sqrt {\frac {\hbar ^{2}}{4m^{*}|\alpha |}}}.}$

It sets the exponential law according to which small perturbations of density of superconducting electrons recover their equilibrium value ψ₀. Thus this theory characterized all superconductors by two length scales. The second one is the penetration depth, λ. It was previously introduced by the London brothers in their London theory. Expressed in terms of the parameters of Ginzburg–Landau model it is

${\displaystyle \lambda ={\sqrt {\frac {m^{*}}{\mu _{0}e^{*2}\psi _{0}^{2}}}}={\sqrt {\frac {m^{*}\beta }{\mu _{0}e^{*2}|\alpha |}}},}$

where ψ₀ is the equilibrium value of the order parameter in the absence of an electromagnetic field. The penetration depth sets the exponential law according to which an external magnetic field decays inside the superconductor.

The original idea on the parameter κ belongs to Landau. The ratio κ = λ/ξ is presently known as the Ginzburg–Landau parameter. It has been proposed by Landau that Type I superconductors are those with 0 < κ < 1/√2, and Type II superconductors those with κ > 1/√2.

Fluctuations[edit]

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.^[4]

Classification of superconductors[edit]

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 H_c. Depending on the geometry of the sample, one may obtain an intermediate state^[5] consisting of a baroque pattern^[6] 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 H_c1 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 H_c2, 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.^[7]

Geometric formulation[edit]

The Ginzburg–Landau functional can be formulated in the general setting of a complex vector bundle over a compact Riemannian manifold.^[8] 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 ${\displaystyle E}$ over a Riemannian manifold ${\displaystyle M}$ with fiber ${\displaystyle \mathbb {C} ^{n}}$, the order parameter ${\displaystyle \psi }$ is understood as a section of the vector bundle ${\displaystyle E}$. The Ginzburg–Landau functional is then a Lagrangian for that section:

${\displaystyle {\mathcal {L}}(\psi ,A)=\int _{M}{\sqrt {|g|}}dx^{1}\wedge \dotsm \wedge dx^{m}\left[\vert F\vert ^{2}+\vert D\psi \vert ^{2}+{\frac {1}{4}}\left(\sigma -\vert \psi \vert ^{2}\right)^{2}\right]}$

The notation used here is as follows. The fibers ${\displaystyle \mathbb {C} ^{n}}$ are assumed to be equipped with a Hermitian inner product ${\displaystyle \langle \cdot ,\cdot \rangle }$ so that the square of the norm is written as ${\displaystyle \vert \psi \vert ^{2}=\langle \psi ,\psi \rangle }$. The phenomenological parameters ${\displaystyle \alpha }$ and ${\displaystyle \beta }$ 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 ${\displaystyle \sigma \in \mathbb {R} }$. The integral is explicitly over the volume form

${\displaystyle *(1)={\sqrt {|g|}}dx^{1}\wedge \dotsm \wedge dx^{m}}$

for an ${\displaystyle m}$-dimensional manifold ${\displaystyle M}$ with determinant ${\displaystyle |g|}$ of the metric tensor ${\displaystyle g}$.

The ${\displaystyle D=d+A}$ is the connection one-form and ${\displaystyle F}$ is the corresponding curvature 2-form (this is not the same as the free energy ${\displaystyle F}$ given up top; here, ${\displaystyle F}$ corresponds to the electromagnetic field strength tensor). The ${\displaystyle A}$ corresponds to the vector potential, but is in general non-Abelian when ${\displaystyle n>1}$, and is normalized differently. In physics, one conventionally writes the connection as ${\displaystyle d-ieA}$ for the electric charge ${\displaystyle e}$ and vector potential ${\displaystyle A}$; in Riemannian geometry, it is more convenient to drop the ${\displaystyle e}$ (and all other physical units) and take ${\displaystyle A=A_{\mu }dx^{\mu }}$ 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 ${\displaystyle \langle \cdot ,\cdot \rangle }$ invariant; so here, ${\displaystyle A}$ is a form taking values in the algebra ${\displaystyle {\mathfrak {su}}(n)}$.

The curvature ${\displaystyle 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

${\displaystyle {\begin{aligned}F&=D\circ D\\&=dA+A\wedge A\\&=\left({\frac {\partial A_{\nu }}{\partial x^{\mu }}}+A_{\mu }A_{\nu }\right)dx^{\mu }\wedge dx^{\nu }\\&={\frac {1}{2}}\left({\frac {\partial A_{\nu }}{\partial x^{\mu }}}-{\frac {\partial A_{\mu }}{\partial x^{\nu }}}+[A_{\mu },A_{\nu }]\right)dx^{\mu }\wedge dx^{\nu }\\\end{aligned}}}$

That is, each ${\displaystyle A_{\mu }}$ is an ${\displaystyle n\times 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

${\displaystyle {\mathcal {L}}(A)=YM(A)=\int _{M}*(1)\vert F\vert ^{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 ^[9]

${\displaystyle D^{*}D\psi ={\frac {1}{2}}\left(\sigma -\vert \psi \vert ^{2}\right)\psi }$

and

${\displaystyle D^{*}F=-\operatorname {Re} \langle D\psi ,\psi \rangle }$

where ${\displaystyle D^{*}}$ is the adjoint of ${\displaystyle D}$, analogous to the codifferential ${\displaystyle \delta =d^{*}}$. Note that these are closely related to the Yang–Mills–Higgs equations.

Specific results[edit]

In string theory, it is conventional to study the Ginzburg–Landau functional for the manifold ${\displaystyle M}$ being a Riemann surface, and taking ${\displaystyle n=1}$; i.e., a line bundle.^[10] The phenomenon of Abrikosov vortices persists in these general cases, including ${\displaystyle M=\mathbb {R} ^{2}}$, where one can specify any finite set of points where ${\displaystyle \psi }$ vanishes, including multiplicity.^[11] The proof generalizes to arbitrary Riemann surfaces and to Kähler manifolds.^{[12][13][14][15]} In the limit of weak coupling, it can be shown that ${\displaystyle \vert \psi \vert }$ converges uniformly to 1, while ${\displaystyle D\psi }$ and ${\displaystyle dA}$ converge uniformly to zero, and the curvature becomes a sum over delta-function distributions at the vortices.^[16] 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 spin^c 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[edit]

When the manifold ${\displaystyle M}$ is a Riemann surface ${\displaystyle M=\Sigma }$, 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 ${\displaystyle d=\partial +{\overline {\partial }}}$. Likewise, the space ${\displaystyle \Omega ^{1}}$ of one-forms over a Riemann surface decomposes into a space that is holomorphic, and one that is anti-holomorphic: ${\displaystyle \Omega ^{1}=\Omega ^{1,0}\oplus \Omega ^{0,1}}$, so that forms in ${\displaystyle \Omega ^{1,0}}$ are holomorphic in ${\displaystyle z}$ and have no dependence on ${\displaystyle {\overline {z}}}$; and vice-versa for ${\displaystyle \Omega ^{0,1}}$. This allows the vector potential to be written as ${\displaystyle A=A^{1,0}+A^{0,1}}$ and likewise ${\displaystyle D=\partial _{A}+{\overline {\partial }}_{A}}$ with ${\displaystyle \partial _{A}=\partial +A^{1,0}}$ and ${\displaystyle {\overline {\partial }}_{A}={\overline {\partial }}+A^{0,1}}$.

For the case of ${\displaystyle n=1}$, where the fiber is ${\displaystyle \mathbb {C} }$ so that the bundle is a line bundle, the field strength can similarly be written as

${\displaystyle F=-\left(\partial _{A}{\overline {\partial }}_{A}+{\overline {\partial }}_{A}\partial _{A}\right)}$

Note that in the sign-convention being used here, both ${\displaystyle A^{1,0},A^{0,1}}$ and ${\displaystyle F}$ are purely imaginary (viz U(1) is generated by ${\displaystyle e^{i\theta }}$ so derivatives are purely imaginary). The functional then becomes

${\displaystyle {\mathcal {L}}\left(\psi ,A\right)=2\pi \sigma \operatorname {deg} L+\int _{\Sigma }{\frac {i}{2}}dz\wedge d{\overline {z}}\left[2\vert {\overline {\partial }}_{A}\psi \vert ^{2}+\left(*(-iF)-{\frac {1}{2}}(\sigma -\vert \psi \vert ^{2}\right)^{2}\right]}$

The integral is understood to be over the volume form

${\displaystyle *(1)={\frac {i}{2}}dz\wedge d{\overline {z}}}$,

so that

${\displaystyle \operatorname {Area} \Sigma =\int _{\Sigma }*(1)}$

is the total area of the surface ${\displaystyle \Sigma }$. The ${\displaystyle *}$ is the Hodge star, as before. The degree ${\displaystyle \operatorname {deg} L}$ of the line bundle ${\displaystyle L}$ over the surface ${\displaystyle \Sigma }$ is

${\displaystyle \operatorname {deg} L=c_{1}(L)={\frac {1}{2\pi }}\int _{\Sigma }iF}$

where ${\displaystyle c_{1}(L)=c_{1}(L)[\Sigma ]\in H^{2}(\Sigma )}$ is the first Chern class.

The Lagrangian is minimized (stationary) when ${\displaystyle \psi ,A}$ solve the Ginzberg–Landau equations

${\displaystyle {\begin{aligned}{\overline {\partial }}_{A}\psi &=0\\*(iF)&={\frac {1}{2}}\left(\sigma -\vert \psi \vert ^{2}\right)\\\end{aligned}}}$

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

${\displaystyle 4\pi \operatorname {deg} L\leq \sigma \operatorname {Area} \Sigma }$.

Roughly speaking, this can be interpreted as an upper limit to the density of the Abrikosov vortecies. One can also show that the solutions are bounded; one must have ${\displaystyle |\psi |\leq \sigma }$.

In string theory[edit]

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;^[17] 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.^[18] 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.^[19] 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.^[20] 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.^[21]

See also[edit]

References[edit]

Papers[edit]

• 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