Evidence of a two-dimensional glass transition in graphene: Insights from molecular simulations

Liquids exhibit a sudden increase in viscosity when cooled fast enough, avoiding thermodynamically predicted route of crystallization. This phenomenon, known as glass transition, leads to the formation of non-periodic structures known as glasses. Extensive studies have been conducted on model materials to understand glass transition in two dimensions. However, despite the synthesis of disordered/amorphous single-atom thick structures of carbon, little attention has been given to glass transition in realistic two-dimensional materials such as graphene. Herein, using molecular dynamics simulation, we demonstrate the existence of glass transition in graphene leading to a realistic two-dimensional glassy structure, namely glassy graphene. We show that the resulting glassy structure exhibits excellent agreement with experimentally realized disordered graphene. Interestingly, this glassy graphene exhibits a wrinkled but stable structure, with reduced thermal vibration in comparison to its crystalline counterpart. We suggest that the topological disorder induced by glass transition governs the unique properties of this structure.

Relationship between bond-breakage correlations and four-point correlations in heterogeneous glassy dynamics: configuration changes and vibration modes

We investigate the dynamic heterogeneities of glassy particle systems in the theoretical schemes of bond breakage and four-point correlation functions. In the bond-breakage scheme, we introduce the structure factor S(b)(q,t) and the susceptibility χ(b)(t) to detect the spatial correlations of configuration changes. Here χ(b)(t) attains a maximum at t=t(b)(max) as a function of time t, where the fraction of the particles with broken bonds φ(b)(t) is about 1/2. In the four-point scheme, treating the structure factor S(4)(q,t) and the susceptibility χ(4)(t), we detect superpositions of the heterogeneity of bond breakage and that of thermal low-frequency vibration modes. While the former grows slowly, the latter emerges quickly to exhibit complex space-time behavior. In two dimensions, the vibration modes extending over the system yield significant contributions to the four-point correlations, which depend on the system size logarithmically. A maximum of χ(4)(t) is attained at t=t(4)(max), where these two contributions become of the same order. As a result, t(4)(max) is considerably shorter than t(b)(max).

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.

Theory of specific heat in glass-forming systems

Experimental measurements of the specific heat in glass-forming systems are obtained from the linear response to either slow cooling (or heating) or to oscillatory perturbations with a given frequency about a constant temperature. The latter method gives rise to a complex specific heat with the constraint that the zero frequency limit of the real part should be identified with thermodynamic measurements. Such measurements reveal anomalies in the temperature dependence of the specific heat, including the so called "specific heat peak" in the vicinity of the glass transition. The aim of this paper is to provide theoretical explanations of these anomalies in general and a quantitative theory in the case of a simple model of glass formation. We first present interesting simulation results for the specific heat in a classical model of a binary mixture glass former. We show that in addition to the formerly observed specific heat peak there is a second peak at lower temperatures which was not observable in earlier simulations. Second, we present a general relation between the specific heat, a caloric quantity, and the bulk modulus of the material, a mechanical quantity, and thus offer a smooth connection between the liquid and amorphous solid states. The central result of this paper is a connection between the micromelting of clusters in the system and the appearance of specific heat peaks; we explain the appearance of two peaks by the micromelting of two types of clusters. We relate the two peaks to changes in the bulk and shear moduli. We propose that the phenomenon of glass formation is accompanied by a fast change in the bulk and the shear moduli, but these fast changes occur in different ranges of the temperature. Last, we demonstrate how to construct a theory of the frequency dependent complex specific heat, expected from heterogeneous clustering in the liquid state of glass formers. A specific example is provided in the context of our model for the dynamics of glycerol. We show that the frequency dependence is determined by the same alpha -relaxation mechanism that operates when measuring the viscosity or the dielectric relaxation spectrum. The theoretical frequency dependent specific heat agrees well with experimental measurements on glycerol. We conclude the paper by stating that there is nothing universal about the temperature dependence of the specific heat in glass formers-unfortunately, one needs to understand each case by itself.

Quantitative theory of a relaxation function in a glass-forming system

We present a quantitative theory for a relaxation function in a simple glass-forming model (binary mixture of particles with different interaction parameters). It is shown that the slowing down is caused by the competition between locally favored regions (clusters) that are long-lived but each of which relaxes as a simple function of time. Without the clusters, the relaxation of the background is simply determined by one typical length, which we deduce from an elementary statistical mechanical argument. The total relaxation function (which depends on time in a nontrivial manner) is quantitatively determined as a weighted sum over the clusters and the background. The "fragility" in this system can be understood quantitatively since it is determined by the temperature dependence of the number fractions of the locally favored regions.

Ergodicity and slowing down in glass-forming systems with soft potentials: no finite-temperature singularities

The aim of this paper is to discuss some basic notions regarding generic glass-forming systems composed of particles interacting via soft potentials. Excluding explicitly hard-core interaction, we discuss the so-called glass transition in which a supercooled amorphous state is formed, accompanied by a spectacular slowing down of relaxation to equilibrium, when the temperature is changed over a relatively small interval. Using the classical example of a 50-50 binary liquid of N particles with different interaction length scales, we show the following. (i) The system remains ergodic at all temperatures. (ii) The number of topologically distinct configurations can be computed, is temperature independent, and is exponential in N. (iii) Any two configurations in phase space can be connected using elementary moves whose number is polynomially bounded in N, showing that the graph of configurations has the small world property. (iv) The entropy of the system can be estimated at any temperature (or energy), and there is no Kauzmann crisis at any positive temperature. (v) The mechanism for the super-Arrhenius temperature dependence of the relaxation time is explained, connecting it to an entropic squeeze at the glass transition. (vi) There is no Vogel-Fulcher crisis at any finite temperature T>0 .

Aging and relaxation in glass-forming systems

We propose that there exists a generic class of glass-forming systems that have competing states (of crystalline order or not) which are locally close in energy to the ground state (which is typically unique). Upon cooling, such systems exhibit patches (or clusters) of these competing states which become locally stable in the sense of having a relatively high local shear modulus. It is in between these clusters where aging, relaxation, and plasticity under strain can take place. We demonstrate explicitly that relaxation events that lead to aging occur where the local shear modulus is low (even negative) and result in an increase in the size of local patches of relative order. We examine the aging events closely from two points of view. On the one hand we show that they are very localized in real space, taking place outside the patches of relative order, and from the other point of view we show that they represent transitions from one local minimum in the potential surface to another. This picture offers a direct relation between structure and dynamics, ascribing the slowing down in glass-forming systems to the reduction in relative volume of the amorphous material which is liquidlike. While we agree with the well-known Adam-Gibbs proposition that the slowing down is due to an entropic squeeze (a dramatic decrease in the number of available configurations), we do not agree with the Adam-Gibbs (or the Volger-Fulcher) formulas that predict an infinite relaxation time at a finite temperature. Rather, we propose that generically there should be no singular crisis at any finite temperature: the relaxation time and the associated correlation length (average cluster size) increase at most superexponentially when the temperature is lowered.

Correlation between dynamic heterogeneity and medium-range order in two-dimensional glass-forming liquids

A glassy state of matter results if crystallization is avoided upon cooling or increasing density. However, the physical factors controlling the ease of vitrification and nature of the glass transition remain elusive. Using numerical simulations of polydisperse hard disks, we find a direct relation between medium-range crystalline ordering and the slow dynamics which characterizes the glass transition. This suggests an intriguing scenario that the strength of frustration controls both the ease of vitrification and nature of the glass transition. Vitrification may be a process of hidden crystalline ordering under frustration, at least in our system.

Statistical mechanics of the glass transition as revealed by a Voronoi tesselation

The statistical mechanics of simple glass forming systems in two dimensions is worked out. The glass disorder is encoded via a Voronoi tesselation, and the statistical mechanics is performed directly in this encoding. The theory provides, without free parameters, an explanation of the glass transition phenomenology, including the identification of two different temperatures, T(g) and T(c) , the first associated with jamming and the second associated with crystallization at very low temperatures.

High-density liquidlike component facilitates plastic flow in a model amorphous silicon system

Molecular dynamics simulations show that plastic deformation of amorphous Si modeled by the Stillinger-Weber potential is very sensitive to the density of the initial unstressed state. Low-density systems exhibit a pronounced yield phenomenon, strain softening, and a dramatic drop in pressure during deformation at a constant volume. This behavior is explained by the presence in every system of a certain concentration of solidlike and liquidlike components, the latter being denser and more amenable to plastic flow.

Crystal phases of a glass-forming Lennard-Jones mixture

We compare the potential energy at zero temperature of a range of crystal structures for a glass-forming binary mixture of Lennard-Jones particles. The lowest-energy ordered state consists of coexisting phases of a single component face centered cubic structure and an equimolar cesium chloride structure. An infinite number of layered crystal structures are identified with energies close to this ground state. We demonstrate that the finite size increase of the energy of the coexisting crystal with incoherent interfaces is sufficient to destabilize this ordered phase in simulations of typical size. Two specific local coordination structures are identified as of possible structural significance in the amorphous state. We observe rapid crystal growth in the equimolar mixture.

Stability and structure of a supercooled liquid mixture in two dimensions

The structural and thermodynamic properties of a two-dimensional binary mixture of soft discs are reported over a range of temperatures down to large supercoolings using constant NPT molecular dynamics simulations. It is shown that the four orders of magnitude increase in the structural relaxation time is not accompanied by any significant increase in translational or orientational order. The phase diagram, calculated in the temperature/composition plane using thermodynamic integration, exhibits a deep eutectic point that is responsible for stabilizing the amorphous state. Voronoi analysis of the low-temperature ground state reveals a structure characterized by a network of linear arrays of fivefold and sevenfold sites. The heat capacity C(P) exhibits an asymmetric peak with a maximum at T(*)=0.55. It is argued that the initial rapid drop in C(P) for T*<0.55 is an equilibrium result and, hence, the peak in the heat capacity corresponds to the existence of an "enthalpy gap" with a characteristic temperature of T* approximately 0.35. This gap results from a minimum volume change associated with an anharmonic fluctuation.

Topological features of structural relaxations in a two-dimensional model atomic glass II

The topological features of atom motions in a high-temperature melt, a sub-cooled melt above T g , and a glass below T g , were analysed in detail by means of a two-dimensional molecular dynamics simulation. A striking analogy was observed between the structure and properties of the liquid-like material separating quasi-ordered domains of atom clusters, and high-angle grain boundaries. The main feature of the structural relaxation below the melting point, both above and below T g was the gradual dissolution and disappearance of the liquid-like material, permitting increasing order in the previously quasi-ordered domains and a growth in their sizes. In these processes, many sequences reminiscent of cancellation of dislocation pairs, or mutual reactions to give more stable sets, were observed.

A molecular dynamics model of melting and glass transition in an idealized two-dimensional material I

In a preparatory study of structural relaxations and plastic flow in a two-dimensional idealized atomic glass, the process of melting and quenching through a glass transition has been studied by computer simulation using a molecular dynamics model. In this model, the transition from a solid to a melt was observed to take place when liquid-like structural elements composed of dipoles of five- and seven-sided Voronoi polygons percolate through the two-dimensional structure of distorted hexagons in the form of strings. Such dipoles constitute discrete elements of excess free volume within which liquid like behaviour is established in the sense of reduced cohesion or local elastic moduli. Upon quenching the melt, the percolation condition of liquid-like regions is retained for under-cooled melts between the melting point and a glass transition temperature below which the percolation condition is broken and the thermal expansion is sharply reduced. The simulation that has used empirical pair potentials characteristic of Cu and Zr has substantially underpredicted the melting and glass transition temperatures and overpredicted the thermal expansion of C u x Zr 1-x type glasses. These defects of the model can be partly attributed to the two-dimensional nature of the material, which stores larger concentrations of free volume than a corresponding three-dimensional material. In spite of these quantitative shortcomings, the model gives valuable insight into the topological features of the local atomic configurations at melting and upon vitrification.