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.