KTHNY theory

The KTHNY-theory describes melting of crystals in two dimensions (2D). The name is derived from the initials of the surnames of John Michael Kosterlitz, David J. Thouless, Bertrand Halperin, David R. Nelson, and A. Peter Young, who developed the theory in the 1970s. It is, beside the Ising model in 2D and the XY model in 2D, one of the few theories, which can be solved analytically and which predicts a phase transition at a temperature formula_1.
Main idea.
Melting of 2D crystals is mediated by the dissociation of topological defects, which destroy the order of the crystal. In 2016, Michael Kosterlitz and David Thouless were awarded with the Nobel prize in physics for their idea, how thermally excited pairs of ``virtual´´ dislocations induce a softening (described by renormalization group theory) of the crystal during heating. The shear elasticity disappears simultaneously with the dissociation of the dislocations, indicating a fluid phase. Based on this work, David Nelson and Bertrand Halperin showed, that the resulting hexatic phase is not yet an isotropic fluid. Starting from a hexagonal crystal (which is the densest packed structure in 2D), the hexatic phase has a six-folded director field, similar to liquid crystals. Orientational order only disappears due to the dissociations of a second class of topological defects, named disclinations. Peter Young calculated the critical exponent of the diverging correlations length at the transition between crystalline and hexatic.
KTHNY theory predicts two continuous phase transitions, thus latent heat and phase coexistence is ruled out. The thermodynamic phases can be distinguished based on discrete versus continuous translational and orientational order. One of the transitions separates a solid phase with quasi-long range translational order and perfect long ranged orientational order from the hexatic phase. The hexatic phase shows short ranged translational order and quasi-long ranged orientational order. The second phase transition separates the hexatic phase from the isotropic fluid, where both, translational and orientational order is short ranged. The system is dominated by critical fluctuations, since for continuous transitions, the difference of energy between the thermodynamic phases disappears in the vicinity of the transition. This implies, that ordered and disordered regions fluctuate strongly in space and time. The size of those regions grows strongly near the transitions and diverges at the transition itself. At this point, the pattern of symmetry broken versus symmetric domains is fractal. Fractals are characterized by a scaling invariance – they appear similar on an arbitrary scale or by arbitrarily zooming in (this is true on any scale larger than the atomic distance). The scale invariance is the basis to use the renormalization group theory to describe the phase transitions. Both transitions are accompanied by spontaneous symmetry breaking. Unlike for melting in three dimensions, translational and orientational symmetry breaking does not need to appear simultaneously in 2D, since two different types of topological defects destroy the different types of order.
Background.
Michael Kosterlitz and David Thouless tried to resolve a contradiction about 2D crystals: on one hand side, the Mermin-Wagner theorem claims that symmetry breaking of a continuous order-parameter cannot exist in two dimensions. This implies, that perfect long range positional order is ruled out in 2D crystals. On the other side, very early computer simulations of Berni Alder and Thomas E. Wainwright indicated crystallization in 2D. The KTHNY theory shows implicitly that periodicity is not a sufficient criterion for a solid (this is already indicated by the existence of amorphous solids like glasses). Following M. Kosterlitz, a finite shear elasticity defines a 2D solid, including quasicrystals in this description.
Structure factor in 2D.
All three thermodynamic phases and their corresponding symmetries can be visualized using the structure factor:formula_2. The double sum runs over all positions of particle pairs I and j and the brackets denote an average about various configurations. The isotropic phase is characterized by concentric rings at formula_3, if formula_4 is the average particle distance calculated by the 2D particle density formula_5. The (closed packed) crystalline phase is characterized by six-fold symmetry based on the orientational order. Unlike in 3D, where the peaks are arbitrarily sharp (formula_6-peaks), the 2D peaks have a finite width described with a Lorenz-curve. This is due to the fact, that the translational order is only quasi-long ranged as predicted by the Mermin-Wagner theorem. The hexatic phase is characterized by six segments, which reflect the quasi-long ranged orientational order. The structure factor of Figure 1 is calculated from the positions of a colloidal monolayer (crosses at high intensity are artefacts from the Fourier transformation due to the finite (rectangular) field of view of the ensemble).
Interaction between dislocations.
To analyse melting due to the dissociation of dislocations, one starts with the energy formula_7 as function of distance between two dislocations. An isolated dislocation in 2D is a local distortions of the six-folded lattice, where neighbouring particles have five- and seven nearest neighbours, instead of six. It is important to note, that dislocations can only be created in pairs, due to topological reasons. A bound pair of dislocations is a local configuration with 5-7-7-5 neighbourhood.
The double sum runs over all positions of defect pairs formula_9 and formula_10, formula_11 measures the distance between the dislocations. formula_12 is the Burgers vector and denotes the orientation of the dislocation at position Orte formula_13. The second term in the brackets brings dislocations to arrange preferentially antiparallel due to energetic reasons. Its contribution is small and can be neglected for large distance between defects. The main contribution stems from the logarithmic term (the first one in the brackets) which describes, how the energy of a dislocation pair diverges with increasing distance. Since the shortest distance between two dislocations is given approximatively by the average particle distance formula_14, the scaling of distances with formula_14 prevents the logarithm formula_16 to become negative. The strength of the interaction is proportional to Young's modulus formula_17 given by the stiffness of the crystal lattice. To create a dislocation from an undisturbed lattice, a small displacement on a scale smaller than the average particle distance formula_14 is needed. The discrete energy associated with this displacement is usually called core energy Energie formula_19 and has to be counted for each of the formula_20 dislocations individually (last term).
An easy argument for the dominating logarithmic term is, that the magnitude of the strain induced by an isolated dislocation decays according to formula_21 with distance. Assuming Hooke's approximation, the associated stress is linear with the strain. Integrating the strain ~1/r gives the energy proportional to the logarithm. The logarithmic distance dependence of the energy is the reason, why KTHNY-theory is one of the few theories of phase transitions which can be solved analytically: in statistical physics one has to calculate partition functions, e.g. the probability distribution for ``all´´ possible configurations of dislocation pairs given by the Boltzmann distribution formula_22. Here, formula_23 is the thermal energy with Boltzmann constant formula_24. For the majority of problems in statistical physics one can hardly solve the partition function due to the enormous amount of particles and degrees of freedoms. This is different in KTHNY theory due to the logarithmic energy functions of dislocations formula_7 and the e-function from the Boltzmann factor as inverse which can be solved easily.
Example.
We want to calculate the mean squared distance between two dislocations considering only the dominant logarithmic term for simplicity:
This mean distance formula_27 tends to zero for low temparatures – dislocations will annihilate and the crystal is free of defects. The expression diverges formula_28, if the denominator tends to zero. This happens, when
formula_29. A diverging distance of dislocations implies, that they are dissociated and do not form a bound pair. The crystal is molten, if several isolated dislocations are thermally excited and the melting temperature formula_30 is given by Young’s modulus:
The dimensionless quantity formula_32 is a universal constant for melting in 2D and is independent of details of the system under investigation. 
This example investigated only an isolated pair of dislocations. In general, a multiplicity of dislocations will appear during melting. The strain field of an isolated dislocation will be shielded and the crystal will get softer in the vicinity of the phase transition; Young’s modulus will decrease due to dislocations. In KTHNY theory, this feedback of dislocations on elasticity, and especially on Young’s modulus acting as coupling constant in the energy function, is described within the framework of renormalization group theory.
Renormalization of elasticity.
If a 2D crystal is heated, ``virtual´´ dislocation pairs will be excited due to thermal fluctuations in the vicinity of the phase transition. Virtual means, that the average thermal energy is not large enough to overcome (two times) the core-energy and to dissociate (unbind) dislocation pairs. Nonetheless, dislocation pairs can appear locally on very short time scales due to thermal fluctuations, before they annihilate again. Although they annihilate, they have a detectable impact on elasticity: they soften the crystal. The principle is completely analogue to calculating the bare charge of the electron in quantum electrodynamics (QED). In QED, the charge of the electron is shielded due to virtual electron-positron pairs due to quantum fluctuations of the vacuum. Roughly spoken one can summarize: If the crystal is softened due to the presence of virtual pairs of dislocation, the probability (fugacity) formula_33 for creating additional virtual dislocations is enhanced, proportional to the Boltzmann factor of the core-energy of a dislocation formula_34. If additional (virtual) dislocations are present, the crystal will get additionally softer. If the crystal is additionally softer, the fugacity will increase further... and so on and so forth.
David Nelson, Bertrand Halperin and independently Peter Young formulated this in a mathematically precise way, using renormalization group theory for the fugacity and the elasticity: In the vicinity of the continuous phase transition, the system becomes critical – this means that it becomes self-similar on all length scales formula_35. Executing a transformation of all length scales by an factor of formula_36, the energy formula_37 and fugacity formula_38 will depend on this factor, but the system has to appear identically, simultaneously due to the self similarity. Especially the energy function (Hamiltonian) of the dislocations have to be invariant in structure. The softening of the system after a length scale transformation (zooming out to visualize a larger area implies to count more dislocations) is now covered in a renormalized (reduced) elasticity. The recursion relation for elasticity and fugacity are:
Similar recursion relations can be derived for the shear modulus and the bulk modulus. formula_41 and formula_42 are Bessel functions, respectively. Depending on the starting point, the recursion relation can run into two directions. formula_43 implies no defects, the ensemble is crystalline. formula_44, implies arbitrary many defects, the ensemble is fluid. The recursion relation have a fix-point at formula_45 with formula_46. Now, formula_47 is the renormalized value instead of the bare one. Figure 2 shows Youngs’modulus as function of the dimensionless control parameter formula_48. It measures the ratio of the repelling energy between two particles and the thermal energy (which was constant in this experiment). It can be interpreted as pressure or inverse temperature. The black curve is a thermodynamic calculation of a perfect hexagonal crystal at formula_49. The blue curve is from computer simulations and shows a reduced elasticity due to lattice vibrations at formula_50. The red curve is the renormalization following the recursion relations, Young's modulus disappears discontinuously to zero at formula_32. Turquoise symbols are from measurements of elasticity in a colloidal monolayer, and confirm the melting point at formula_52.
Interaction between disclinations.
The system enters the hexatic phase after the dissociation of dislocations. To reach the isotropic fluid, dislocations (5-7-pairs) have to dissociate into disclinations, consisting of isolated 5-folded and isolated 7-folded particles. Similar arguments for the interaction of disclinations compared to dislocations can be used. Again, disclinations can only be created as pairs due to topological reasons. Starting with the energy formula_53 as function of distance between two disclinations one finds:
The logarithmic term is again dominating. The sign of the interaction gives attraction or repulsion for the winding numbers formula_55 and formula_56 of the five- and seven-folded disclinations in a way that ``charges´´ with opposite sign have attraction. The overall strength is given by the stiffness against twist. The coupling constant formula_57 is called Frank's constant, following the theory of liquid crystals. 
formula_58 is the discrete energy of a dislocation to dissociate into two disclinations. The squared distance of two disclinations can be calculated the same way, as for dislocations, only the prefactor, denoting the coupling constant, has to be changed accordingly. It diverges for formula_59. The system is molten from the hexatic phase into the isotropic liquid, if unbound disclinations are present. This transition temperature formula_60 is given by Frank's constant:
formula_62 is again a universal constant. Figure 3 shows measurements of the orientational stiffness of a colloidal monolayer; Frank's constant drops below this universal constant at formula_60.
Critical exponents.
Typically, Kosterlitz–Thouless transitions have a continuum of critical points which can be characterised by self-similar grains of disordered and ordered regions. In second order phase transitions, the correlation length measuring the size of those regions diverges algebraically:
Here, formula_65 is the transition temperature and formula_66 is a critical exponent. Another special feature of Kosterlitz–Thouless transitions is, that translational and orientational correlation length in 2D diverge exponentially (see also hexatic phase for the definition of those correlation functions):
The critical exponent becomes formula_68 for the diverging translational correlation length at the hexatic – crystalline transition. D. Nelson and B. Halperin predicted, that Frank's constant diverges exponentially with formula_69 at formula_60, too. The red curve shows a fit of experimental data covering the critical behaviour; the critical exponent is measured to be formula_71. This value is compatible with the prediction of KTHNY theory within the error bars. The orientational correlation length at the hexatic – isotropic transition is predicted to diverge with an exponent formula_72. This rational value is compatible with mean-field-theories and implies that a renormalization of Frank's constant is not necessary. The increasing shielding of orientational stiffness due to disclinations has not to be taken into account – this is already done by dislocations which are frequently present at formula_60. Experiments measured a critical exponent of formula_74.
KTHNY-theory has been tested in experiment and in computer simulations. For short range particle interaction (hard discs), simulations found a weakly first order transition for the hexatic – isotropic transition, slightly beyond KTHNY-theory.