Issue
Natl Sci Open
Volume 2, Number 3, 2023
Special Topic: Glasses—Materials and Physics
Article Number 20220056
Number of page(s) 9
DOI https://doi.org/10.1360/nso/20220056
Published online 27 April 2023

© The Author(s) 2023. Published by Science Press and EDP Sciences

Licence Creative CommonsThis is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

INTRODUCTION

Phase transitions are often accompanied by change of structural orders, with profound impacts on the properties of materials. And changes of structural orders are often sought as indicators for phase transitions. In the cases of melting [1-6] or sublimation [7] of crystalline solids, the material undergoes obvious change of states, when the structures change from ordered crystalline lattice to disordered liquid or gas. In other cases, however, a crystal can be continuously transformed into a disordered solid or glass, by progressively introducing defects into the ordered structure by ion irradiation [8-10], doping [11, 12] or mechanical treatments [13-17] until no more long-range crystalline order left in the system. It is difficult to determine the point where a crystal qualitatively transforms into a glass, as the system remains in the solid state during these processes while the structure evolves continuously. Such methods of creating glasses from crystalline solids are in clear contrast to the conventional way of producing glasses by rapid quenching, when the system continuously transitions from flowing liquids to solid glasses with little changes in the structures. The order-to-disorder transition in solids thus provides a unique opportunity to study the structure of glasses, with the underlying crystalline structure as a well-defined reference.

Colloids are an ideal model system to study structural orders during phase transitions. The melting of colloidal crystals in both 2D and 3D have been intensively studied in-situ, using temperature-sensitive Poly-N-isopropylacrylamide (PNIPAM) particles [1, 18-22]. Yunker et al. [12] show that a small fraction of dopants in a colloidal crystal is sufficient to give rise to strong dynamical heterogeneity, a hallmark of glassy materials. Similar emergence of glassy dynamics is also observed when bond strength disorder is introduced into a crystalline lattice [23]. In both these experiments, the structure of the system remains largely crystalline. Recent computer simulations by Zhang et al. [14] show that sharp transitions between polycrystals and disordered glasses can be achieved by progressively compressing a binary mixture of hard and soft particles. In the simulations, the soft particles can assume one of two diameters depending on the local stress, thus creating two distinctive length scales in the system. It is not clear if similar sharp transition can still be observed when the size of the particles changes continuously.

In this Research Article, we study the transition from a quasi-2D colloidal crystal to a colloidal glass by simultaneously increasing the size mismatch between the constituent particles and the packing fraction of the system, a process analogous to compression-induced crystal-to-glass transition. The order-to-disorder transition is identified by a peak in the fluctuations of hexagonal order parameter and the loss of the long-range correlations of the orientational order in the system. As the system remains in solid state, diffusion contributes little to the transition from ordered crystal to disordered glass, in contrast to the melting of crystals. The transition is facilitated by the generation and growth of disordered clusters with increasing particle size mismatch. The system qualitatively transforms into glass phase when the disordered clusters become connected and percolate the system.

EXPERIMENT

We employ PNIPAM hydrogel particles whose size increases with decreasing temperature. Two types of PNIPAM microspheres, denoted as type A and type B, are synthesized with different crosslinker densities, which result in different thermosensitivities of the particles [24]. A higher crosslinker density leads to a lower thermosensitivity. Figure 1 plots the temperature dependence of diameters of the two types of PNIPAM particles measured by dynamic light scattering. The steeper slope for particle B (red line) indicates higher temperature sensitivity than particle A (black line). The sizes of the two types of particles are close around 28.5℃. In the temperature range between 29℃ and 24℃, the size ratio of the two types of PNIPAM particles varies from 0.94 to 1.33, covering the range needed for crystal-to-glass transition in two dimensions [25-28].

The samples consist of a binary mixture of type A and type B particles, with a number ratio nA/nB= 3:2. 0.5 μL colloidal suspension is placed between two glass coverslips, and spreads under the gravity and capillary forces, forming a monolayer. The samples are then hermetically sealed using optical glue (Norland 63). Before data acquisition, the samples are annealed at 28.5℃ on a microscope stage for 3 h to obtain uniform crystalline structures. The temperature of the sample is controlled by a resistive objective heater (Bioptechs). The sample temperature is gradually lowered from 28.5℃ to 24.0℃ at steps of 0.2℃, and the structure of the sample evolves from ordered triangular lattice to disordered glass. The nominal packing fraction of the sample changes from 0.9 to 1.55 for the same range of temperature. Snapshots of the sample images are shown in the insets of Figure 1. At each temperature, images of the samples are digitally recorded using a CCD camera mounted on the microscope at 60 frames per second for 500 s. There are 4000 particles in the field of view, whose positions are extracted using standard particle tracking technique [29]. During the experiments, the long-time mean-square-displacement (MSD) for all samples are less than 5% of the average particle separation, below the Lindemann melting criterion.

thumbnail Figure 1

Temperature dependence of the diameters of PNIPAM microspheres with different thermosensitivities. The diameters are obtained by dynamic light scattering. Particle B (red squares) is more sensitive to temperature than particle A (black squares), the straightlines are linear fits to the measured data. Snapshots of sample structures at 28.5℃, 27.5℃, and 25.0℃ are shown in the insets from right to left.

RESULTS AND DISCUSSION

Figure 2 shows the measured pair correlation function g(r) and the Fourier transformation of the sample images at different temperatures. At 28.5℃, when the diameters of particles A and B are close, the sample structure exhibits clear features of a triangular crystal. The first three peaks of g(r) correspond to the 1, 3 and 2 times, respectively, of the lattice constant. The bright discrete pattern with six-fold symmetry in the Fourier transformation also confirms that the sample is in crystalline state. As the temperature is lowered, the particles of the two species swell at different rates, which increases the size mismatch between particles and the overall stress in the sample. The g(r) peaks broaden with temperature, and the second and third peaks eventually merge into a round peak at 27.9℃, which is a sign of order-to-disorder transition, as the distortion of the crystalline lattice becomes comparable to the particle separations [30-32]. When the temperature is further lowered, the g(r) peaks at large distances vanish when long-range periodicity is completely lost, and the g(r) curves become indistinguishable from those of typical colloidal glasses. The evolution of the system structure from ordered crystal to glass is also confirmed by the diffraction patterns, when the distinct point patterns blur into halo rings.

The crystalline order in the system is quantified by the bond orientational order parameter ψ6=1Ntot j=1Ntotk=1nn|ei6θjk|/nn. Here, θjk is angle between the bond connecting nearest-neighboring particles j and k, and x-axis, nn is coordination number of reference particle j, and Ntot is the total number of particles [33, 12]. Figure 3A plots the evolution of ψ6 when the sample temperature decreases from 28.5℃ to 24℃. The averaged ψ6 decreases from nearly 1, which corresponds to perfect crystalline order, to about 0.6, the typical value for binary 2D glasses. Figure 3B plots the spatial correlation function of ψ6, g6(r=|rirj|)=Ψ6i*(ri)ψ6j(rj), where ri and rj are positions of particles i and j, and the angle bracket indicates ensemble average. At high temperatures when the particle mismatch is small, the envelope of g6 decays following a power law, indicating long-range correlation of local orientational order. As the size mismatch increases, the g6 starts to decay exponentially, as in the case of glasses. The change from a power law to exponential decay can be roughly determined between 27.9℃ and 27.5℃.

The temporal fluctuation of the bond orientational order is characterized by χ6=Ntot((ψ¯6)2ψ¯62), where ψ¯6 is the averaged ψ6 in a single frame, i.e., ψ¯6=j=1Ntot|ψ6|/Ntot, and bracket 〈〉 indicates time average. To avoid errors from finite size of the sample, we calculated a series of χ6 within different size of sub-boxes, then extrapolated the χ6 to infinite scale limit. The susceptibility χ6 thus obtained is plotted in Figure 3C. A pronounced peak of χ6 is observed at 27.9℃, signaling the order-to-disorder transition, and coincides with the structural crossover observed in g(r) and g6(r).

The structural features shown in Figure 3 are similar to those observed in the melting of 2D colloidal crystals [21], but the microscopic mechanism in the solid order-to-disorder transition is fundamentally different. During the 2D melting, diffusion of the topological defects plays a critical role in creating an isotropic fluid from an anisotropic crystal [3, 4]. The defects are mostly dissociated disclination pairs created by the increased free volumes. During the crystal-to-glass transition, on the other hand, the defects are created by the compression of the lattices that reduces the overall free volume in the system, which then prohibits the diffusion of these defects. During the experiment, all the topological defects are “pinned” in place with no noticeable diffusions.

During the crystal-to-glass transitions, local crystalline order is progressively distorted when the disordered phase grows in size. The glass phase contains not only topological defects with 5 or 7 neighbors, but also distorted hexagonal units. We employ disordered clusters to quantify the structural evolution during the transition [14]. Specifically, a particle is considered ordered if it has three or more nearest neighbors with the product of the local bond order |ψ6i*(ri)ψ6k(rk)|0.65; otherwise, it is considered disordered or a defect. This criterion is based on the observation that in 2D binary glasses, the typical 〈ψ6〉 is 0.6. Figures 4A–4D plot the spatial distribution of the ordered (blue) and disordered (red) particles at different temperatures. The disordered phase emerges from the crystalline background, and the two phases remain separated even after the disordered clusters percolate the whole system. Therefore, the crystal-to-glass transition occurs via one-step transition, instead of a two-step transition as the melting of 2D crystals.

As the temperature is lowered from 28.5℃, both the number and size of the disordered clusters increase, eventually percolating the whole system. Figure 5A plots the typical sizes of disordered clusters, the largest cluster size and the average cluster size, respectively, at different temperatures. The cluster size increases exponentially towards the percolation transition. Figure 5B plots the size distribution of disordered clusters at temperatures above T = 27.7℃. The fraction of large clusters increases as the temperature decreases, and the distributions can be well fitted to a power law truncated by an exponential function. We compare the size distribution of disordered clusters to that of randomly generated clusters. Specifically, after the fraction of disordered particles is determined from experiments, the same fraction of particles in a certain frame is randomly tagged as disordered particles. The particles are then grouped into disordered clusters based on their spatial proximity so that neighboring disordered particles belong to the same disordered cluster. The size distributions of the disordered clusters are qualitatively different from those of randomly generated clusters of the same configurations, as shown in the inset of Figure 5B, with much higher probability at large n, which suggests that as the system is progressively compressed, new disordered particles emerge not randomly, but preferentially next to existing large clusters. The disordered clusters thus act as “nuclei” of the disordered phase until they percolate the whole system, when the system qualitatively transitions from an ordered crystal to a disordered glass. This preferential attachment to existing disordered clusters is probably due to the stress localization near the boundary, which makes it easier to generate a disordered particle near the boundary than inside the crystalline domain. The rapid growth of cluster sizes is driven by the increasing length of boundary between the ordered and disordered phases, and the inter-connection of existing clusters by a few disordered clusters approaching the percolation transition. Figure 6 plots the measured percolation probability of the disordered clusters. For each frame, the percolation status for the disordered clusters is determined. The percolation probability is defined as the fraction of time the disordered clusters are found to percolate the system during the experiment time window, with P=tp/ttot where tp is the number of frames that “percolated disordered cluster" occurs at a certain temperature, and ttot is total frame number of measurement, which is 30,000 in our experiments. The percolation probability jumps sharply from near zero to unity between 27.9℃ and 27.5℃, in the same range where the χ6 indicates a qualitative transition from crystalline phase to glass phase. The range of temperature for the percolation transition is practically the same when the threshold of the criterion for disordered particles changes from 0.6 to 0.7.

thumbnail Figure 2

Evolution of structures during the crystal-to-glass transition. (A) The pair correlation function g(r) at different temperatures. The merge of the second and third peaks (dotted rectangle) indicates the loss of crystalline order. (B) The diffraction patterns of the samples at four representative temperatures.

thumbnail Figure 3

The evolution of the bond orientational order during the crystal-to-glass transition. (A) Bond orientational order ψ6 as a function of sample temperature; (B) the envelope of the spatial correlation function g6 at different temperatures; (C) the susceptibility of the bond orientational order, χ6 as a function of temperature.

thumbnail Figure 4

The distribution of disordered clusters at different temperatures. Particles in disordered clusters are shown in red, while regions with high crystalline order are shown in blue.

thumbnail Figure 5

Sizes of disordered clusters. (A) Typical cluster sizes, including average sizes (black squares) and maximum sizes (red solid circles), as a function of sample temperature; (B) size distribution of disordered clusters at different temperatures. The distributions are fitted to a power law with exponential truncation (solid lines). (Inset) Comparison between the cluster size distributions obtained in experiments (solid triangles) and randomly generated clusters (empty triangles).

thumbnail Figure 6

Percolation probability of the disordered clusters as a function of temperature.

CONCLUSIONS

In summary, we employ PNIPAM colloidal particles with different thermosensitivities to study compression-induced crystal-to-glass transition in solid in a quasi-2D system. Compression creates disordered clusters in the crystalline background, which grow in size and numbers until eventually percolate the system. The percolation of the disordered clusters coincides with the qualitative transition from crystal to glass identified by the structural signatures and corresponding susceptibility. This one-step transition is qualitatively different from the two-step melting of crystals at 2D. It is also in contrast to the conventional glass transition, where the boundary between the liquid and solid phases is vague. In practice, qualitatively altering the structures of a solid by compression usually requires extremely high pressure, and some recent studies have shown evidence of order-to-disorder transition in solid materials [34]. In these experiments, the evolution of the structures on atomic level is difficult to assess. The thermo-sensitive colloids thus provide a convenient model to study the microscopic details of these phenomena. The change of dynamical and thermodynamical properties such as phonon modes during the crystal-glass transition are worthy of further studies in the future.

Funding

This work was supported by the National Natural Science Foundation of China (11874395 and 12174434) and the Strategic Priority Research Program of Chinese Academy of Sciences (XDB33000000).

Author contributions

K.C. and C.W. conceived the concept. C.W. and L.C. carried out the experiments. All the authors participated in the analysis of the experimental data. K.C. and L.C. prepared the manuscript.

Conflict of interest

The authors declare no conflict of interest.

References

All Figures

thumbnail Figure 1

Temperature dependence of the diameters of PNIPAM microspheres with different thermosensitivities. The diameters are obtained by dynamic light scattering. Particle B (red squares) is more sensitive to temperature than particle A (black squares), the straightlines are linear fits to the measured data. Snapshots of sample structures at 28.5℃, 27.5℃, and 25.0℃ are shown in the insets from right to left.

In the text
thumbnail Figure 2

Evolution of structures during the crystal-to-glass transition. (A) The pair correlation function g(r) at different temperatures. The merge of the second and third peaks (dotted rectangle) indicates the loss of crystalline order. (B) The diffraction patterns of the samples at four representative temperatures.

In the text
thumbnail Figure 3

The evolution of the bond orientational order during the crystal-to-glass transition. (A) Bond orientational order ψ6 as a function of sample temperature; (B) the envelope of the spatial correlation function g6 at different temperatures; (C) the susceptibility of the bond orientational order, χ6 as a function of temperature.

In the text
thumbnail Figure 4

The distribution of disordered clusters at different temperatures. Particles in disordered clusters are shown in red, while regions with high crystalline order are shown in blue.

In the text
thumbnail Figure 5

Sizes of disordered clusters. (A) Typical cluster sizes, including average sizes (black squares) and maximum sizes (red solid circles), as a function of sample temperature; (B) size distribution of disordered clusters at different temperatures. The distributions are fitted to a power law with exponential truncation (solid lines). (Inset) Comparison between the cluster size distributions obtained in experiments (solid triangles) and randomly generated clusters (empty triangles).

In the text
thumbnail Figure 6

Percolation probability of the disordered clusters as a function of temperature.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.