Melting and Liquid Structure in two Dimensions

Melting, freezing, and dynamics of two-dimensional dipole systems in screening bulk media

This paper reports on the molecular dynamics simulations of classical two-dimensional (2D) electric dipole systems. The properties of 2D systems with bare (nonscreened) and screened dipole-dipole interactions have been investigated. Based on the polygon construction method, we present simulation results on the phase transition, and we locate the melting and freezing points of 2D dipole systems in terms of a polygon disorder parameter, with the polygon disorder parameter being the sum of nontriangular polygon order parameters. It was found that the phase transition of the system occurs when the polygon disorder parameter has a value 0.165. This result was cross-checked by using both local and overall orientational order parameters. We also identified that the value of the average local orientational order parameter at the phase transition point is 0.67. These results are valid for the ordinary (bare) dipole-dipole interaction as well as the screened dipole-dipole interaction, and they are expected to be general for other 2D systems with repulsive pair interaction. We observed that both melting and freezing points shift to lower values of temperature due to screening. In the liquid state, the radial distribution function and polygon construction method show the loss of order in a structure as screening becomes more severe. Furthermore, the impact of screening on the system's collective excitation spectra and diffusive characteristics at liquid and solid states has been studied. Results show the decrease in the values of both longitudinal and transverse sound speeds and the emergence of anomalous superdiffusive motion in the liquid state due to screening.

Two-Dimensional Nanoparticle Supracrystals: A Model System for Two-Dimensional Melting

In a Langmuir trough, successive compression cycles can drive a two-dimensional (2D) nanoparticle supracrystal (NPSC) closer to its equilibrium structure. Here, we show a series of equilibrated 2D NPSCs consisting of gold NPs of uniform size, varying solely in the length of their alkanethiol ligands. The ordering of the NPSC is governed by the ligand length, thus providing a model system to investigate the nature of 2D melting in a system of NPs. As the ligand length increases the supracrystal transitions from a crystalline to a liquid-like phase with evidence of a hexatic phase at an intermediate ligand length. The phase change is interpreted as an entropy-driven phenomenon associated with steric constraints between ligand shells. The density of topological defects scales with ligand length, suggesting an equivalence between ligand length and temperature in terms of melting behavior. On the basis of this equivalence, the experimental evidence indicates a two-stage 2D melting of NPSCs.

Fluctuations of orientational order and clustering in a two-dimensional colloidal system under quenched disorder

Using both video microscopy of superparamagnetic colloidal particles confined in two dimensions and corresponding computer simulations of repulsive parallel dipoles, we study the formation of fluctuating orientational clusters and topological defects in the context of the KTHNY-like melting scenario under quenched disorder. We analyze cluster densities, average cluster sizes, and the population of noncluster particles, as well as the development of defects, as a function of the system temperature and disorder strength. In addition, the probability distribution of clustering and orientational order is presented. We find that the well-known disorder-induced widening of the hexatic phase can be traced back to the distinct development characteristics of clusters and defects along the melting transitions from the solid phase to the hexatic phase to the isotropic fluid.

Comparison of 2D melting criteria in a colloidal system

We use super-paramagnetic spherical particles which are arranged in a two-dimensional monolayer at a water/air interface to investigate the crystal to liquid phase transition. According to the KTHNY theory a crystal melts in thermal equilibrium by two continuous phase transitions into the isotropic liquid state with an intermediate phase, commonly known as the hexatic phase. We verify the significance of several criteria based on dynamical and structural properties to identify the crystal-hexatic and hexatic-isotropic liquid phase transitions for the same experimental data of the given setup. The criteria are the bond orientational correlation function, the Larson-Grier criterion, the 2D dynamic Lindemann parameter, the bond orientational susceptibility, the 2D Hansen-Verlet rule, the Löwen-Palberg-Simon criterion as well as a criterion based on the shape factor of Voronoi cells and Minkowski functionals. For our system with long-range repulsion, the bond order correlation function and bond order susceptibility work best to identify the hexatic-isotropic liquid transition and the 2D dynamic Lindemann parameter identifies unambiguously the hexatic-crystalline transition.

Geometrical defects in two-dimensional melting of many-particle Yukawa systems

We present a theoretical polygon construction analysis of two-dimensional melting and freezing transitions in many-particle Yukawa systems. Two-dimensional melting transitions can be characterized as proliferation of geometrical defects--nontriangular polygons, obtained by removing unusually long bonds in the triangulation of particle positions. A liquid state is characterized by the temperature-independent number of quadrilaterals and linearly increasing number of pentagons. We analyze specific types of vertices, classified by the type and distribution of polygons surrounding them, and determine temperature dependencies of their concentrations. Solid-liquid phase transitions are followed by the peaks in the abundances of certain types of vertices.

Morphology of monolayer films on quasicrystalline surfaces from the phase field crystal model

We present a computational study of the morphology of adsorbed monolayers on quasicrystalline surfaces with five- and seven-fold symmetry. The phase field crystal model is employed to first simulate the growth of the quasicrystal surfaces and then to elastically couple a two-dimensional film to the substrate. We find several distinct pseudomorphic phases that depend on the position of adsorption sites as well as the strength of the monolayer/substrate interaction, and quantify them by computing local order parameters. In qualitative agreement with recent experiments using colloids in quasiperiodic light fields, we find that the formation of quasicrystalline order is greatly inhibited on seven-fold surfaces.

Roles of bond orientational ordering in glass transition and crystallization

It is widely believed that crystallization in three dimensions is primarily controlled by positional ordering, and not by bond orientational ordering. In other words, bond orientational ordering is usually considered to be merely a consequence of positional ordering and thus has often been ignored. This one-order-parameter (density) description may be reasonable when we consider an equilibrium liquid-solid transition, but may not be enough to describe a metastable state and the kinetics of the transition. Here we propose that bond orientational ordering can play a key role in (i) crystallization, (ii) the ordering to quasi-crystal and (iii) vitrification, which occurs under rather weak frustration against crystallization. In a metastable supercooled state before crystallization, a system generally tends to have bond orientational order at least locally as a result of a constraint of dense packing. For a system interacting with hard-core repulsions, the constraint is intrinsically of geometrical origin and thus the basic physics is the same as nematic ordering of rod-like particles upon densification. Furthermore, positional ordering is easily destroyed even by weak frustration such as polydispersity and anisotropic interactions which favour a symmetry not consistent with that of the equilibrium crystal. Thus we may say that vitrification can be achieved by disturbing and prohibiting long-range positional ordering. Even in such a situation, bond orientational ordering still survives, accompanying its critical-like fluctuations, which are the origin of dynamic heterogeneity for this case. This scenario naturally explains both the absence of positional order and the development of bond orientational order upon cooling in a supercooled state. Although our argument is speculative in nature, we emphasize that this physical picture can coherently explain crystallization, vitrification, quasi-crystallization and their relationship in a natural manner. For a strongly frustrated system, even bond orientational order can be destroyed. Even in such a case there may still appear a structural signature of dense packing, which is linked to slow dynamics.

Polygon construction to investigate melting in two-dimensional strongly coupled dusty plasma

The polygon construction method of Glaser and Clark is used to characterize melting and crystallization in a two-dimensional (2D) strongly coupled dusty plasma. Using particle positions measured by video microscopy, bonds are identified by triangulation, and unusually long bonds are deleted. The resulting polygons have three or more sides. Geometrical defects, which are polygons with more than three sides, are found to proliferate during melting. Pentagons are found in liquids, where they tend to cluster with other pentagons. Quadrilaterals are a less severe defect, so that disorder can be characterized by the ratio of quadrilaterals to pentagons. This ratio is found to be less in a liquid than in a solid or a superheated solid. Another measure of disorder is the abundance of different kinds of vertices, according to the type of polygons that adjoin there. Unexpectedly, spikes are observed in the abundance of certain vertex types during rapid temperature changes. Hysteresis, revealed by a plot of a disorder parameter vs temperature, is examined to study sudden heating. The hysteresis diagram also reveals features suggesting a possibility of latent heat in the melting and rapid cooling processes.

Structural signature of slow dynamics and dynamic heterogeneity in two-dimensional colloidal liquids: glassy structural order

Glassy states are formed if crystallization is avoided upon cooling or increasing density. However, the physical factors controlling the ease of vitrification and the nature of glass transition remain elusive. Among various glass-forming systems, colloidal liquids are one of the most ideal glass-forming systems because of the simplicity and controllability of the interactions. We use numerical simulations of two-dimensional polydisperse and binary hard discs to tackle both of these longstanding questions. For polydisperse systems, we systematically control the polydispersity, which can be regarded as the strength of frustration effects on crystallization. We reveal that crystal-like hexatic order grows in size and lifetime with an increase in the colloid volume fraction or with a decrease in polydispersity (or frustration). We stress that hexatic ordering in hard disc systems is a direct consequence of dense packing and a manifestation of low configurational entropy. Our study suggests an intriguing scenario that the strength of frustration controls both the ease of vitrification and the nature of the glass transition. Vitrification may be a process of hidden crystal-like ordering under frustration for this system. This may provide not only a physical basis for glass formation, but also an answer to another longstanding question on the structure of amorphous materials: 'order in disorder' may be an intrinsic feature of a glassy state of many materials. For binary mixtures, on the other hand, the relevant structural feature linked to slow dynamics is not hexatic order, but an amorphous structure of low structural entropy. These results suggest that slow dynamics is associated with bond orientational order linked to the crystal for a weakly frustrated system, whereas to amorphous structures of low configurational entropy for a strongly frustrated system. This suggests an intrinsic link between structure and dynamics in glass-forming materials: slow dynamics is linked to structuring ('glassy ordering') towards low configurational entropy. We discuss the nature of 'glassy order' responsible for slow dynamics.

Melting of multilayer colloidal crystals confined between two walls

Video microscopy is employed to study the melting behaviors of multilayer colloidal crystals composed of diameter-tunable microgel spheres confined between two walls. We systematically explore film thickness effects on the melting process and on the phase behaviors of single crystal and polycrystalline films. Thick films (>4 layers) are observed to melt heterogeneously, while thin films (≤4 layers) melt homogeneously, even for polycrystalline films. Grain-boundary melting dominates other types of melting processes in polycrystalline films thicker than 12 layers. The heterogeneous melting from dislocations is found to coexist with grain-boundary melting in films between 5- and 12-layers. In dislocation melting, liquid nucleates at dislocations and forms lakelike domains embedded in the larger crystalline matrix; the "lakes" are observed to diffuse, interact, merge with each other, and eventually merge with large strips of liquid melted from grain boundaries. Thin film melting is qualitatively different: thin films homogeneously melt by generating many small defects which need not nucleate at grain boundaries or dislocations. For three- and four-layer thin films, different layers are observed to have the same melting point, but surface layers melt faster than bulk layers. Within our resolution, two- to four-layer films appear to melt in one step, while monolayers melt in two steps with an intermediate hexatic phase.

Structural origin of dynamic heterogeneity in three-dimensional colloidal glass formers and its link to crystal nucleation

The physical understanding of glass transition remains a major challenge of physics and materials science. Among various glass-forming liquids, a colloidal liquid interacting with hard-core repulsion is now regarded as one of the most ideal model systems. Here we study the structure and dynamics of three-dimensional polydisperse colloidal liquids by Brownian dynamics simulations. We reveal that medium-range crystalline bond orientational order of the hexagonal close packed structure grows in size and lifetime with increasing packing fraction. We show that dynamic heterogeneity may be a direct consequence of this transient structural ordering, which suggests its origin is thermodynamic rather than kinetic. We also reveal that nucleation of crystals preferentially occurs in regions of high medium-range order, reflecting the low crystal-liquid interfacial energy there. These findings may shed new light not only on the fundamental nature of the glass transition, but also the mechanism of crystal nucleation.

Crystallization in three- and two-dimensional colloidal suspensions

Despite progress in the understanding of crystal nucleation and crystal growth since the first theories for nucleation were developed, an exact quantitative prediction of the nucleation rates in most systems has remained an unsolved problem. Colloidal suspensions show a phase behavior that is analogous to atomic or molecular systems and serve accordingly as ideal model systems for studying crystal nucleation with an accuracy and depth on a microscopic scale that is hard to reach for atomic or molecular systems. Due to the mesoscopic size of colloidal particles they can be studied in detail on the single-particle level and their dynamics is strongly slowed down in comparison with atomic or molecular systems, such that the formation of a crystal nucleus can be followed in detail. In this review, recent progress in the study of homogeneous and heterogeneous crystal nucleation in colloids and the controlled growth of crystalline colloidal structures is reviewed. All this work has resulted in unprecedented insights into the early stage of nucleation and it is also relevant for a deeper understanding of soft matter materials in general as well as for possible applications based on colloidal suspensions.

Melting of two-dimensional tunable-diameter colloidal crystals

Melting of two-dimensional colloidal crystals is studied by video microscopy. The samples were composed of microgel spheres whose diameters could be temperature tuned, and whose pair potentials were measured to be short ranged and repulsive. We observed two-step melting from the crystal to a hexatic phase and from the hexatic to the liquid phase as a function of the temperature-tunable volume fraction. The translational and orientational susceptibilities enabled us to definitively determine the phase transition points, avoiding ambiguities inherent in other analyses and resolving a "dislocation precursor stage" in the solid phase that some of the traditional analyses may incorrectly associate with the hexatic phase. A prefreezing stage of the liquid with ordered patches was also found.

Crystallization and melting of bacteria colonies and Brownian bugs

Motivated by the existence of remarkably ordered cluster arrays of bacteria colonies growing in Petri dishes and related problems, we study the spontaneous emergence of clustering and patterns in a simple nonequilibrium system: the individual-based interacting Brownian bug model. We map this discrete model into a continuous Langevin equation which is the starting point for our extensive numerical analyses. For the two-dimensional case we report on the spontaneous generation of localized clusters of activity as well as a melting-freezing transition from a disordered or isotropic phase to an ordered one characterized by hexagonal patterns. We study in detail the analogies and differences with the well-established Kosterlitz-Thouless-Halperin-Nelson-Young theory of equilibrium melting, as well as with another competing theory. For that, we study translational and orientational correlations and perform a careful defect analysis. We find a nonstandard one-stage, defect-mediated transition whose nature is only partially elucidated.

Frank's constant in the hexatic phase

Using videomicroscopy data of a two-dimensional colloidal system the bond-order correlation function G{6} is calculated and used to determine both the orientational correlation length xi{6} in the liquid phase and the modulus of orientational stiffness, Frank's constant F{A}, in the hexatic phase. The latter is an anisotropic fluid phase between the crystalline and the isotropic liquid phase. F{A} is found to be finite within the hexatic phase, takes the value 72/pi at the hexatic<-->isotropic liquid phase transition, and diverges at the hexatic<-->crystal transition as predicted by the Kosterlitz-Thouless-Halperin-Nelson-Young theory. This is a quantitative test of the mechanism of breaking the orientational symmetry by disclination unbinding.

Sixfold bond orientational properties of a model liquid crystal in the dimensional crossover of B phases: A computer simulation study

A wide range of NPT simulations of a bead necklace liquid crystal model in the crystal B, smectic B, smectic A, and nematic phases have been performed. Systems with up to 21 600 molecules have been studied to observe the behavior of slowly decaying spatial correlation functions. The pair correlation function and its in-plane restriction are consistent with a crystalline phase made of independent two-dimensional crystalline layers. Smectic B phase is studied by the bond orientational pair correlation functions g(6) and its extension g(6ext). The first reaches a constant value, which seems to rule out a classical hexatic phase. The latter shows a power-law decay within the layers: its typical decay exponent (eta(6ext)) is evaluated. Relationships between multiple harmonics of the C(6n) order parameter have been evaluated through the whole range of existence of B phases (crystalline and smectic): the extension to the crystalline phase holds and provides an excellent fit of the simulation data.

Computer simulations of two-dimensional melting with dipole-dipole interactions

We perform molecular dynamics and Monte Carlo simulations of two-dimensional melting with dipole-dipole interactions. Both static and dynamic behaviors are examined. In the isotropic liquid phase, the bond orientational correlation length xi 6 and susceptibility chi 6 are measured, and the data are fitted to the theoretical ansatz. An algebraic decay is detected for both spatial and temporal bond orientational correlation functions in an intermediate temperature regime, and it provides an explicit evidence for the existence of the hexatic phase. From the finite-size scaling analysis of the global bond orientational order parameter, the disclination unbinding temperature Ti is estimated. In addition, from dynamic Monte Carlo simulations of the positional order parameter, we extract the critical exponents at the dislocation unbinding temperature Tm. All the results are in agreement with those from experiments and support the Kosterlitz-Thouless-Halperin-Nelson-Young (KTHNY) theory.

Discrete elastic model for two-dimensional melting

We present a network model for the study of melting and liquid structure in two dimensions, the first in which the presence and energy of topological defects (dislocations and disclinations) and of geometrical defects (elemental voids) can be independently controlled. Interparticle interaction is via harmonic springs and control is achieved by Monte Carlo moves which springs can either be orientationally "flipped" between particles to generate topological defects, or can be "popped" in force-free shape, to generate geometrical defects. With the geometrical defects suppressed the transition to the liquid phase occurs via disclination unbinding, as described by the Kosterlitz-Thouless-Halperin-Nelson-Young model and found in soft potential two-dimensional (2D) systems, such as the dipole-dipole potential [H. H. von Grünberg, Phys. Rev. Lett. 93, 255703 (2004)]. By contrast, with topological defects suppressed, a disordering transition, the Glaser-Clark condensation of geometrical defects [M. A. Glaser and N. A. Clark, Adv. Chem. Phys. 83, 543 (1993); M. A. Glaser, (Springer-Verlag, Berlin, 1990), Vol. 52, p. 141], produces a state that accurately characterizes the local liquid structure and first-order melting observed in hard-potential 2D systems, such as hard disk and the Weeks-Chandler-Andersen (WCA) potentials (M. A. Glaser and co-workers, see above). Thus both the geometrical and topological defect systems play a role in melting. The present work introduces a system in which the relative roles of topological and geometrical defects and their interactions can be explored. We perform Monte Carlo simulations of this model in the isobaric-isothermal ensemble, and present the phase diagram as well as various thermodynamic, statistical, and structural quantities as a function of the relative populations of geometrical and topological defects. The model exhibits a rich phase behavior including hexagonal and square crystals, expanded crystal, dodecagonal quasicrystal, and isotropic liquid phases. In this system the geometrical defects effectively control the melting, reducing the solid-liquid transition temperature by a factor of relative to the topological-only case. The local structure of the dense liquid has been investigated and the results are compared to that from simulations of WCA systems.

Melting of quasi-two-dimensional crystalline Pb supported on liquid Ga

Experimental studies have shown that the Pb monolayer that segregates in the PbGa alloy liquid-vapor interface forms a two-dimensional hexagonal crystal that melts at 341 K, and it has been speculated that the disordered phase formed is hexatic. This paper reports the results of simulation studies of the in-plane structure of the outermost stratum of the liquid-vapor interface of a dilute Pb in Ga alloy. These simulations are based on four major improvements to a previous study. First, the simulation studies involve considerably more atoms and considerably longer equilibration runs than considered in the previous work of Chekmarev, Oxtoby, and Rice. Second, a more accurate nonlocal pseudopotential representation of the interactions in the system is used. Third, the amplitude of the out-of-plane motion of the Pb atoms is constructed to have the observed value. Fourth, an approximation to the role of the liquid Ga substrate is provided by adding a layer of Ga atoms to the layer of Pb atoms. The results of our simulation studies show that the Ga layer adjacent to the Pb layer has a profound influence on that layer's properties. In particular, it is shown that in the two-layer PbGa system the Pb layer forms, at low temperature, a stable two-dimensional crystal on top of liquid Ga. This two-dimensional crystal melts at a temperature close to that found experimentally. It is found that the crystalline Pb layer is transformed to the liquid state via two intermediate hexatic phases that differ in the magnitude of the bond orientation order. Each of the phase transitions along this melting pathway is first order. The temperature range over which each hexatic phase is stable is small. The profound influence of out-of-plane motion is demonstrated by a comparison of the results of simulations of a quasi-two-dimensional (Q2D) and of a strictly two-dimensional monolayer of Pb. The melting transition in the Q2D one-layer system is first order, directly to the liquid, with no intervention of a hexatic phase. The melting transition in the strictly 2D system involves two stages: a first-order transition to an intermediate hexatic phase followed by melting of the hexatic to a liquid phase. The latter transition is continuous over a small temperature range. An examination of the role of defects in the melting process reveals a picture rather different from that postulated in the Kosterlitz-Thouless-Halperin-Nelson-Young theory of 2D melting.

Voronoi cell volume distribution and configurational entropy of hard-spheres

The Voronoi cell volume distributions for hard-disk and hard-sphere fluids have been studied. The distribution of the Voronoi free volume vf, which is the difference between the actual cell volume and the minimal cell volume at close packing, is well described by a two-parameter (2gamma) or a three-parameter (3gamma) gamma distribution. The free parameter m in both the 2gamma and 3gamma models is identified as the "regularity factor." The regularity factor is the ratio of the square of the mean and the variance of the free volume distribution, and it increases as the cell volume distribution becomes narrower. For the thermodynamic structures, the regularity factor increases with increasing density and it increases sharply across the freezing transition, in response to the onset of order. The regularity factor also distinguishes between the dense thermodynamic structures and the dense random or quenched structures. The maximum information entropy (max-ent) formalism, when applied to the gamma distributions, shows that structures of maximum information entropy have an exponential distribution of vf. Simulations carried out using a swelling algorithm indicate that the dense random-packed states approach the distribution predicted by the max-ent formalism, though the limiting case could not be realized in simulations due to the structural inhomogeneities introduced by the dense random-packing algorithm. Using the gamma representations of the cell volume distribution, we check the numerical validity of the Cohen-Grest expression [M. H. Cohen and G. S. Grest, Phys. Rev. B 20, 1077 (1979)] for the cellular (free volume) entropy, which is a part of the configurational entropy. The expression is exact for the hard-rod system, and a correction factor equal to the dimension of the system, D, is found necessary for the hard-disk and hard-sphere systems. Thus, for the hard-disk and hard-sphere systems, the present analysis establishes a relationship between the precisely defined Voronoi free volume (information) entropy and the thermodynamic entropy. This analysis also shows that the max-ent formalism, when applied to the free volume entropy, predicts an exponential distribution which is approached by disordered states generated by a swelling algorithm in the dense random-packing limit.

Structural transitions in two-dimensional hard-sphere systems

We spread randomly noncharged steel particles (diameter, 1.59 mm) on a silicon wafer to form a two-dimensional hard-sphere system. The particle structure versus the particle coverage was monitored. We observed the particle structural transition from liquidlike to triangular-lattice crystal-like with increasing particle coverage by analyzing the particle structure factor. The particle coverage at which the structural transition occurs was quantified by the curves of S(max) (A) and G6 (A); S(max) is the amplitude of the first peak of the structure factor (depicting the particle positional order), and G6 is the bond orientation order parameter. We also conducted a Monte Carlo simulation study. The Monte Carlo simulation results show good agreement with the experimental results at low particle area fractions. However, at high area fractions, the experimentally observed particle structure is less organized than that generated by simulations.

Detection and characterization of structural changes in the hard-disk fluid under freezing and melting conditions

The fluid of two-dimensional hard disks is investigated over a range of densities by Monte Carlo simulations in order to detect and characterize structural changes which take place when the condition of freezing and melting is approached. A novel method is proposed based on the use of the Voronoi tessellation and a certain shape factor which turns out to be a clear indicator of the presence of different underlying substructures (domains). Close to the freezing condition the probability distribution of the shape factor develops a second distinct maximum corresponding to a predominant presence of near-regular hexagons, whereas the original peak, having its origin primarily in pentagons and distorted hexagons, diminishes and disappears at melting density.

Elastic behavior of a two-dimensional crystal near melting

Using positional data from video microscopy, we determine the elastic moduli of two-dimensional colloidal crystals as a function of temperature. The moduli are extracted from the wave-vector-dependent normal-mode spring constants in the limit q-->0 and are compared to the renormalized Young's modulus of the Kosterlitz-Thouless-Halperin-Nelson-Young theory. An essential element of this theory is the universal prediction that Young's modulus must approach 16 pi at the melting temperature. This is indeed observed in our experiment.