Hexatic phase

The hexatic phase is a state of matter that is between the solid and the isotropic liquid phases in two dimensional systems of particles. It is characterized by two order parameters: a short-range positional and a quasi-long-range orientational (sixfold) order. More generally, a hexatic is any phase that contains sixfold orientational order, in analogy with the nematic phase (with twofold orientational order).

It is a fluid phase, since the shear modulus and the Young's modulus vanish due to the dissociation of dislocations. It is an anisotropic phase, since there exists a director field with sixfold symmetry. The existence of the director field implies that an elastic modulus against drilling or torsion exists within the plane, that is usually called Frank's constant after Frederick C. Frank in analogy to liquid crystals. The ensemble becomes an isotropic liquid (and Frank's constant becomes zero) after the dissociation of disclinations at a higher temperature (or lower density). Therefore, the hexatic phase contains dislocations but no disclinations.

The KTHNY theory of two-step melting by i) destroying positional order and ii) destroying orientational order was developed by John Michael Kosterlitz, David J. Thouless, Bertrand Halperin, David Robert Nelson and A. P. Young in theoretical studies about topological defect unbinding two dimensions. In 2016, Kosterlitz and Thouless were awarded with the Nobel Prize in Physics (together with Duncan Haldane) for the idea that melting in 2D is mediated by topological defects. The hexatic phase was predicted by D. Nelson and B. Halperin; it does not have a strict analogue in three dimensions.

Order parameter

The hexatic phase can be described by two order parameters, where the translational order is short ranged (exponential decay) and the orientational order is quasi-long ranged (algebraic decay).

phase	translational order	orientational order	defects
crystalline	quasi-long range: $G_{\vec {G}}({\vec {R}})\propto R^{-\eta _{\vec {G}}}$	long range: $\lim _{r\to \infty }G_{6}({\vec {r}})\propto const.$	no defects
hexatic (anisotropic fluid)	short range: $G_{\vec {G}}({\vec {R}})\propto e^{-R/\xi _{\vec {G}}}$	quasi-long range: $G_{6}({\vec {r}})\propto r^{-\eta _{6}}$	dislocations
isotropic fluid	short range: $G_{\vec {G}}({\vec {R}})\propto e^{-R/\xi _{\vec {G}}}$	short range: $G_{6}({\vec {r}})\propto e^{-r/\xi _{6}}$	dislocations and disclinations

Translational order

Dislocations destroy translational order (shearing along the red arrows), but orientational order is still visible, indicated by black lines in one direction (upper part). Disclinations additionally destroy orientational order (lower part).

If the position of atoms or particles is known, then the translational order can be determined with the translational correlation function $G_{\vec {G}}({\vec {R}})$ as function of the distance between lattice site at place ${\vec {R}}$ and the place ${\vec {0}}$ , based on the two-dimensional density function $\rho _{\vec {G}}({\vec {R}})=e^{i{\vec {G}}\cdot [{\vec {R}}+{\vec {u}}({\vec {R}})]}$ in reciprocal space:

G_{\vec {G}}({\vec {R}})=\langle \rho _{\vec {G}}({\vec {R}})\cdot \rho _{\vec {G}}^{\ast }({\vec {0}})\rangle

The vector ${\vec {R}}$ points to a lattice site within the crystal, where the atom is allowed to fluctuate with an amplitude ${\vec {u}}({\vec {R}})$ by thermal motion. ${\vec {G}}$ is a reciprocal vector in Fourier space. The brackets denote a statistical average about all pairs of atoms with distance R.

The translational correlation function decays fast, i. e. exponential, in the hexatic phase. In a 2D crystal, the translational order is quasi-long range and the correlation function decays rather slow, i. e. algebraic; It is not perfect long range, as in three dimensions, since the displacements ${\vec {u}}({\vec {R}})$ diverge logarithmically with systems size at temperatures above T=0 due to the Mermin-Wagner theorem.

A disadvantage of the translational correlation function is, that it is strictly spoken only well defined within the crystal. In the isotropic fluid, at the latest, disclinations are present and the reciprocal lattice vector is not defined any more.

Orientational order

The j nearest neighbours of the i-th particle define the sixfold director field.

The orientational order can be determined by the local director field of a particle at place ${\vec {r}}_{i}$ , if the angles $\theta _{ij}$ are taken, given by the bond to the $N_{i}$ nearest neighbours in sixfolded space, normalized with the number of nearest neighbours:

\Psi ({\vec {r}}_{i})={\frac {1}{N_{i}}}\sum _{j=1}^{N_{i}}e^{i6\theta _{ij}}

$\Psi$ is a complex number of magnitude $|\Psi ({\vec {r}})|\leq 1$ and the orientation of the six-folded director is given by the phase. In a hexagonal crystal, this is nothing else but the crystal-axes. The local director field disappears for a particle with five or seven nearest neighbours, as given by dislocations and disclinations $\Psi \sim 0$ , except a small contribution due to thermal motion. The orientational correlation function between two particles i and k at distance ${\vec {r}}={\vec {r}}_{i}-{\vec {r}}_{k}$ is now defined using the local director field:

G_{6}({\vec {r}})=\langle \Psi ({\vec {r}}_{i})\cdot \Psi ^{\ast }({\vec {r}}_{k})\rangle

Again, the brackets denote the statistical average about all pairs of particles with distance $|{\vec {r}}|=r$ . All three thermodynamic phases can be identified with this orientational correlation function: it does not decay in the 2D crystal but takes a constant value (shown in blue in the figure). The stiffness against local torsion is arbitrarily large, Franks's constant is infinity. In the hexatic phase, the correlation decays with a power law (algebraic). This gives straight lines in a log-log-plot, shown in green in the Figure. In the isotropic phase, the correlations decay exponentially fast, this are the red curved lines in the log-log-plot (in a lin-log-plot, it would be straight lines). The discrete structure of the atoms or particles superimposes the correlation function, given by the minima at half integral distances $a$ . Particles which are poorly correlated in position, are also poorly correlated in their director.

Orientational correlation of the local director field as function of distance, plotted in blue for the crystal, green for the hexatic phase and red for the isotropic fluid.

External links

Nobel-talk 2016 from M. Kosterlitz

References

Kosterlitz, J M; Thouless, D J (12 June 1972). "Long range order and metastability in two dimensional solids and superfluids. (Application of dislocation theory)". Journal of Physics C: Solid State Physics. 5 (11): L124–L126. Bibcode:1972JPhC....5L.124K. doi:10.1088/0022-3719/5/11/002. ISSN 0022-3719.
Kosterlitz, J M; Thouless, D J (12 April 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. ISSN 0022-3719.
Kosterlitz, J M (21 March 1974). "The critical properties of the two-dimensional xy model". Journal of Physics C: Solid State Physics. 7 (6): 1046–1060. Bibcode:1974JPhC....7.1046K. doi:10.1088/0022-3719/7/6/005. ISSN 0022-3719.
Nelson, David R.; Kosterlitz, J. M. (7 November 1977). "Universal Jump in the Superfluid Density of Two-Dimensional Superfluids". Physical Review Letters. 39 (19): 1201–1205. Bibcode:1977PhRvL..39.1201N. doi:10.1103/physrevlett.39.1201. ISSN 0031-9007.
Halperin, B. I.; Nelson, David R. (10 July 1978). "Theory of Two-Dimensional Melting". Physical Review Letters. 41 (2): 121–124. Bibcode:1978PhRvL..41..121H. doi:10.1103/physrevlett.41.121. ISSN 0031-9007.
Nelson, David R.; Halperin, B. I. (1 February 1979). "Dislocation-mediated melting in two dimensions". Physical Review B. 19 (5): 2457–2484. Bibcode:1979PhRvB..19.2457N. doi:10.1103/physrevb.19.2457. ISSN 0163-1829.
Young, A. P. (15 February 1979). "Melting and the vector Coulomb gas in two dimensions". Physical Review B. 19 (4): 1855–1866. Bibcode:1979PhRvB..19.1855Y. doi:10.1103/physrevb.19.1855. ISSN 0163-1829.
Jaster, A. (2004). "The hexatic phase of the two-dimensional hard disk system". Physics Letters A. 330 (1–2): 120–125. arXiv:cond-mat/0305239. Bibcode:2004PhLA..330..120J. doi:10.1016/j.physleta.2004.07.055. ISSN 0375-9601. S2CID 119522893.
Keim, P.; Maret, G.; Grünberg, H.H.v. (2007). "Frank's constant in the hexatic phase". Physical Review E. 75 (3): 031402. arXiv:cond-mat/0610332. Bibcode:2007PhRvE..75c1402K. doi:10.1103/PhysRevE.75.031402. PMID 17500696. S2CID 5886990.
Gasser, U.; Eisenmann, C.; Maret, G.; Keim, P. (2010). "Melting of crystals in two dimensions". ChemPhysChem. 11 (5): 963–970. doi:10.1002/cphc.200900755. PMID 20099292.
Kosterlitz, M. (2016). "Commentary on Ordering, metastability and phase transitions in two-dimensional systems". Journal of Physics C. 28 (48): 481001. doi:10.1088/0953-8984/28/48/481001. PMID 27665689. S2CID 46754095.