Berezinskii–Kosterlitz–Thouless transition

This article includes a list of general references, but it lacks sufficient corresponding inline citations. Please help to improve this article by introducing more precise citations. (November 2019) (Learn how and when to remove this message)

The Berezinskii–Kosterlitz–Thouless (BKT) transition is a phase transition of the two-dimensional (2-D) XY model in statistical physics. It is a transition from bound vortex-antivortex pairs at low temperatures to unpaired vortices and anti-vortices at some critical temperature. The transition is named for condensed matter physicists Vadim Berezinskii, John M. Kosterlitz and David J. Thouless.^[1] BKT transitions can be found in several 2-D systems in condensed matter physics that are approximated by the XY model, including Josephson junction arrays and thin disordered superconducting granular films.^[2] More recently, the term has been applied by the 2-D superconductor insulator transition community to the pinning of Cooper pairs in the insulating regime, due to similarities with the original vortex BKT transition.

The critical density of the BKT transition in the weakly interacting system reads^[3]

${\displaystyle n_{\text{c}}={\frac {mT}{2\pi }}\ln {\frac {\xi }{mU}}}$

where the dimensionless constant was found to be ${\displaystyle \xi =380\pm 3}$.^[4]

Work on the transition led to the 2016 Nobel Prize in Physics being awarded to Thouless and Kosterlitz; Berezinskii died in 1981.

The XY model is a two-dimensional vector spin model that possesses U(1) or circular symmetry. This system is not expected to possess a normal second-order phase transition. This is because the expected ordered phase of the system is destroyed by transverse fluctuations, i.e. the Nambu-Goldstone modes associated with this broken continuous symmetry, which logarithmically diverge with system size. This is a specific case of what is called the Mermin–Wagner theorem in spin systems.

Rigorously the transition is not completely understood, but the existence of two phases was proved by McBryan & Spencer (1977) and Fröhlich & Spencer (1981).

In the XY model in two dimensions, a second-order phase transition is not seen. However, one finds a low-temperature quasi-ordered phase with a correlation function (see statistical mechanics) that decreases with the distance like a power, which depends on the temperature. The transition from the high-temperature disordered phase with the exponential correlation to this low-temperature quasi-ordered phase is a Kosterlitz–Thouless transition. It is a phase transition of infinite order.

In the 2-D XY model, vortices are topologically stable configurations. It is found that the high-temperature disordered phase with exponential correlation decay is a result of the formation of vortices. Vortex generation becomes thermodynamically favorable at the critical temperature ${\displaystyle T_{\text{c}}}$ of the Kosterlitz–Thouless transition. At temperatures below this, vortex generation has a power law correlation.

Kosterlitz–Thouless transitions is described as a dissociation of bound vortex pairs with opposite circulations, called vortex–antivortex pairs, first described by Vadim Berezinskii. In these systems, thermal generation of vortices produces an even number of vortices of opposite sign. Bound vortex–antivortex pairs have lower energies than free vortices, but have lower entropy as well. In order to minimize free energy, ${\displaystyle F=E-TS}$, the system undergoes a transition at a critical temperature, ${\displaystyle T_{\text{c}}}$. Below ${\displaystyle T_{\text{c}}}$, there are only bound vortex–antivortex pairs. Above ${\displaystyle T_{\text{c}}}$, there are free vortices.

There is an elegant thermodynamic argument for the Kosterlitz–Thouless transition. The energy of a single vortex is ${\displaystyle \kappa \ln(R/a)}$, where ${\displaystyle \kappa }$ is a parameter that depends upon the system in which the vortex is located, ${\displaystyle R}$ is the system size, and ${\displaystyle a}$ is the radius of the vortex core. One assumes ${\displaystyle R\gg a}$. In the 2D system, the number of possible positions of a vortex is approximately ${\displaystyle (R/a)^{2}}$. From Boltzmann's entropy formula, ${\displaystyle S=k_{\rm {B}}\ln W}$ (with W is the number of states), the entropy is ${\displaystyle S=2k_{\rm {B}}\ln(R/a)}$, where ${\displaystyle k_{\rm {B}}}$ is the Boltzmann constant. Thus, the Helmholtz free energy is

${\displaystyle F=E-TS=(\kappa -2k_{\rm {B}}T)\ln(R/a).}$

When ${\displaystyle F>0}$, the system will not have a vortex. On the other hand, when ${\displaystyle F<0}$, entropic considerations favor the formation of a vortex. The critical temperature above which vortices may form can be found by setting ${\displaystyle F=0}$ and is given by

${\displaystyle T_{\text{c}}={\frac {\kappa }{2k_{\rm {B}}}}.}$

The Kosterlitz–Thouless transition can be observed experimentally in systems like 2D Josephson junction arrays by taking current and voltage (I-V) measurements. Above ${\displaystyle T_{\text{c}}}$, the relation will be linear ${\displaystyle V\sim I}$. Just below ${\displaystyle T_{c}}$, the relation will be ${\displaystyle V\sim I^{3}}$, as the number of free vortices will go as ${\displaystyle I^{2}}$. This jump from linear dependence is indicative of a Kosterlitz–Thouless transition and may be used to determine ${\displaystyle T_{\text{c}}}$. This approach was used in Resnick et al.^[5] to confirm the Kosterlitz–Thouless transition in proximity-coupled Josephson junction arrays.

The following discussion uses field theoretic methods. Assume a field φ(x) defined in the plane which takes on values in ${\displaystyle S^{1}}$, so that ${\displaystyle \phi (x)}$ is identified with ${\displaystyle \phi (x)+2\pi }$. That is, the circle is realized as ${\displaystyle S^{1}=\mathbb {R} /2\pi \mathbb {Z} }$.

The energy is given by

${\displaystyle E=\int {\frac {1}{2}}\nabla \phi \cdot \nabla \phi \,d^{2}x}$

and the Boltzmann factor is ${\displaystyle \exp(-\beta E)}$.

Taking a contour integral ${\displaystyle \textstyle \oint _{\gamma }d\phi =\oint _{\gamma }{\frac {d\phi }{dx}}dx}$ over any contractible closed path ${\displaystyle \gamma }$, we would expect it to be zero (for example, by the fundamental theorem of calculus. However, this is not the case due to the singular nature of vortices (which give singularities in ${\displaystyle \phi }$).

To render the theory well-defined, it is only defined up to some energetic cut-off scale ${\displaystyle \Lambda }$, so that we can puncture the plane at the points where the vortices are located, by removing regions with size of order ${\displaystyle 1/\Lambda }$. If ${\displaystyle \gamma }$ winds counter-clockwise once around a puncture, the contour integral ${\displaystyle \textstyle \oint _{\gamma }d\phi }$ is an integer multiple of ${\displaystyle 2\pi }$. The value of this integer is the index of the vector field ${\displaystyle \nabla \phi }$.

Suppose that a given field configuration has ${\displaystyle N}$ punctures located at ${\displaystyle x_{i},i=1,\dots ,N}$ each with index ${\displaystyle n_{i}=\pm 1}$. Then, ${\displaystyle \phi }$ decomposes into the sum of a field configuration with no punctures, ${\displaystyle \phi _{0}}$ and ${\displaystyle \textstyle \sum _{i=1}^{N}n_{i}\arg(z-z_{i})}$, where we have switched to the complex plane coordinates for convenience. The complex argument function has a branch cut, but, because ${\displaystyle \phi }$ is defined modulo ${\displaystyle 2\pi }$, it has no physical consequences.

Now,

${\displaystyle E=\int {\frac {1}{2}}\nabla \phi _{0}\cdot \nabla \phi _{0}\,d^{2}x+\sum _{1\leq i<j\leq N}n_{i}n_{j}\int {\frac {1}{2}}\nabla \ \arg(z-z_{i})\cdot \nabla \arg(z-z_{j})\,d^{2}x}$

If ${\displaystyle \textstyle \sum _{i=1}^{N}n_{i}\neq 0}$, the second term is positive and diverges in the limit ${\displaystyle \Lambda \to \infty }$: configurations with unbalanced numbers of vortices of each orientation are never energetically favoured.

However, if the neutral condition ${\displaystyle \textstyle \sum _{i=1}^{N}n_{i}=0}$ holds, the second term is equal to ${\displaystyle \textstyle -2\pi \sum _{1\leq i<j\leq N}n_{i}n_{j}\ln(|x_{j}-x_{i}|/L)}$, which is the total potential energy of a two-dimensional Coulomb gas. The scale L is an arbitrary scale that renders the argument of the logarithm dimensionless.

Assume the case with only vortices of multiplicity ${\displaystyle \pm 1}$. At low temperatures and large ${\displaystyle \beta }$ the distance between a vortex and antivortex pair tends to be extremely small, essentially of the order ${\displaystyle 1/\Lambda }$. At large temperatures and small ${\displaystyle \beta }$ this distance increases, and the favoured configuration becomes effectively the one of a gas of free vortices and antivortices. The transition between the two different configurations is the Kosterlitz–Thouless phase transition, and the transition point is associated with an unbinding of vortex-antivortex pairs.

• Березинский, В. Л. (1970), "Разрушение дальнего порядка в одномерных и двумерных системах с непрерывной группой симметрии I. Классические системы", ЖЭТФ (in Russian), 59 (3): 907–920. Translation available: Berezinskii, V. L. (1971), "Destruction of long-range order in one-dimensional and two-dimensional systems having a continuous symmetry group I. Classical systems" (PDF), Sov. Phys. JETP, 32 (3): 493–500, Bibcode:1971JETP...32..493B
• Березинский, В. Л. (1971), "Разрушение дальнего порядка в одномерных и двумерных системах с непрерывной группой симметрии II. Квантовые системы", ЖЭТФ (in Russian), 61 (3): 1144–1156. Translation available: Berezinskii, V. L. (1972), "Destruction of long-range order in one-dimensional and two-dimensional systems having a continuous symmetry group II. Quantum systems" (PDF), Sov. Phys. JETP, 34 (3): 610–616, Bibcode:1972JETP...34..610B
• Kosterlitz, J. M.; Thouless, D. J. (1973), "Ordering, metastability and phase transitions in two-dimensional systems", Journal of Physics C: Solid State Physics, 6 (7): 1181–1203, Bibcode:1973JPhC....6.1181K, doi:10.1088/0022-3719/6/7/010
• McBryan, O.; Spencer, T. (1977), "On the decay of correlations in SO(n)-symmetric ferromagnets", Commun. Math. Phys., 53 (3): 299, Bibcode:1977CMaPh..53..299M, doi:10.1007/BF01609854, S2CID 119587247
• B. I. Halperin, D. R. Nelson, Phys. Rev. Lett. 41, 121 (1978)
• A. P. Young, Phys. Rev. B 19, 1855 (1979)
• Resnick, D.J.; Garland, J.C.; Boyd, J.T.; Shoemaker, S.; Newrock, R.S. (1981), "Kosterlitz Thouless Transition in Proximity Coupled Superconducting Arrays", Phys. Rev. Lett., 47 (21): 1542, Bibcode:1981PhRvL..47.1542R, doi:10.1103/PhysRevLett.47.1542
• Fröhlich, Jürg; Spencer, Thomas (1981), "The Kosterlitz–Thouless transition in two-dimensional abelian spin systems and the Coulomb gas", Comm. Math. Phys., 81 (4): 527–602, Bibcode:1981CMaPh..81..527F, doi:10.1007/bf01208273, S2CID 73555642
• Z. Hadzibabic; et al. (2006), "Berezinskii–Kosterlitz–Thouless crossover in a trapped atomic gas", Nature, 41 (7097): 1118–21, arXiv:cond-mat/0605291, Bibcode:2006Natur.441.1118H, doi:10.1038/nature04851, PMID 16810249, S2CID 4314014
• M. Mondal; et al. (2011), "Role of the vortex-core energy on the Beresinkii-Kosterlitz-Thouless transition in thin films of NbN", Phys. Rev. Lett., 107 (21): 217003, arXiv:1108.0912, Bibcode:2011PhRvL.107u7003M, doi:10.1103/PhysRevLett.107.217003, PMID 22181915, S2CID 34729666

• J.V. Jose, 40 Years of Berezinskii–Kosterlitz–Thouless Theory, World Scientific, 2013, ISBN 978-981-4417-65-5
• H. Kleinert, Gauge Fields in Condensed Matter, Vol. I, " SUPERFLOW AND VORTEX LINES", pp. 1–742, World Scientific (Singapore, 1989); Paperback ISBN 9971-5-0210-0 (also available online: Vol. I. Read pp. 618–688);
• H. Kleinert, Multivalued Fields in Condensed Matter, Electrodynamics, and Gravitation, World Scientific (Singapore, 2008) (also available online: here)