Introduction

From a microscopic point of view, the structure of supercooled liquids and glasses is amorphous. Even though the atomic arrangement shows a local ordering that can extend even beyond first neighbour atoms1, it globally retains both the continuous translational and rotational symmetries that are proper to the liquid state. Crystallization breaks up these symmetries that become finite, and the structure rearranges towards a long-range atomic order and a thermodynamically stable phase. The time evolution of this process depends on the system and on its thermodynamic conditions, and it can be considered as the reverse of the medal of the glass transition and the main limit to glass stability2. In fact, the timescales of crystallization span over a wide interval, ranging from geologically stable systems, e.g. ref. 3, to nanosecond crystallizing materials, e.g. ref. 4.

Crystallization can be described by the combination of two processes: nucleation and growth. Spontaneous fluctuations in the amorphous system lead to the formation of small crystallites. Within the classical nucleation theory (CNT)5,6, these nuclei become stable when a critical size is exceeded, i.e. when the nucleation free energy barrier is overcome, then nuclei grow, leading to the crystallization of the whole system. In CNT, critical nuclei are assumed spherical, and the nucleation barrier is determined by a balance between surface and volume free energy terms. Experiments and simulations in several systems show that this picture seems valid for moderate supercooling, e.g. ref. 7, whereas it starts to fail in the deep supercooled region. Here, nucleation and growth processes seem more complex, involving also collective non-diffusive rearrangements8,9. Finally, recent papers point out the pivotal role of inhomogeneities in the supercooled liquids as triggers for the nucleation10,11,12, since the presence of these heterogeneous regions could ease the nucleation process.

The relevance of crystallization goes well beyond fundamental condensed matter physics and follows from the universality of the glassy state in nature and technology13. Typical examples come from geology, where the crystallization of volcanic magmas strongly affects the eruptive style of volcanoes14, and from material science, where the rapid crystallizing properties of some chalcogenide glasses are considered to develop fast and reliable permanent memories with nanoseconds switching time, e.g. refs 4 and 15.

In this paper we focus on the kinetics of the isothermal crystallization process in the deep supercooled liquid, i.e. for a temperature T ≪ Tm. In this regime, the viscosity is so high that the system is macroscopically solid and structural rearrangements are still so slow that atoms can be thought as frozen. Nevertheless, the structure of an amorphous solid is never really arrested16, and local non diffusive relaxations can lead to crystallization, as observed in metallic glasses17. In general, disordered systems show a hierarchy of excitations down to very small frequency that can contribute to the origin of many of the complex characteristics of the glasses18.

As a benchmark system, we choose vitreous germania v-GeO2. Like v-SiO2, v-GeO2 is a covalent oxide glass and a prototype of the strong network forming systems19. However, with respect to v-SiO2, v-GeO2 has a rather accessible Tm = 1388 K and a Tg ≃ 818 K19. From a structural point of view, v-GeO2 glassy network is based on Ge(O1/2)4 tetrahedra bound together in a corner-sharing network, which is preserved also in the supercooled liquid20. Crystalline GeO2 presents two stable polymorphs at room pressure and temperature: a rutile-like tetragonal structure (P42/mnm)21, and an α-quartz-like structure (P3221)22. The latter is also the stable phase for T ≥ 1281 K.

Starting from the glass, we approached the supercooled liquid by heating the system up to Texp = 1100 K. At this temperature, the viscosity is very high, about 107 Pa · s, while the diffusion coefficient is about 10−18 m2 · s−1, see ref. 23. In this condition, GeO2 is still in a substantially arrested state with dynamical and structural properties very similar to those of the glass24,25. However, with increasing time the system starts to crystallize, and we observed the kinetics of this process by acquiring a set of static structure factors for about 67 h. Results show the emergence of an α-quartz phase in a continuous process that reorganizes the amorphous matrix, eventually leading to a mixed system with a large number of crystalline domains and a small fraction of low-density amorphous regions. The time evolution of both the crystalline and amorphous fraction were interpreted within an empirical model for the crystallization kinetics. This approach provides a very good description of the experimental data and identifies a non linear predator-prey mechanism between crystal and amorphous where the density variation acts as limiting barrier.

Experiments and Results

The static structure factor S(2θ) of v-GeO2 was firstly measured at room temperature to exclude any appreciable crystallization of the original glass. The temperature was then slowly raised up to 975 K, monitoring the structure to detect any trace of crystallization. Finally, the sample was quickly heated to Texp = 1100 K with a slope of 20 K/min, and then the 67 hours long isothermal measurement was initiated. In order to properly trace the time evolution of the crystallization, we chose two different acquisition times. During a first period of 16 h, the acquisition time was set to Δt = 5 minutes. This is the minimum time to get a good statistics, and short enough to provide detailed view of the beginning of the process. Once the crystallized fraction was clearly visible, the acquisition time was lengthened to Δt = 30 min.

The time evolution of the S(2θ) during the isothermal measurement is reported in Fig. 1. On increasing time, the initially amorphous S(2θ) shows the growth of a crystalline phase through the appearance of Bragg peaks, clearly visible in the low 2θ part of the diffraction pattern, below about 60°. The peak intensities increase and eventually saturate but the peaks pattern remains the same, without intermediate phases. At high 2θ the Debye-Waller factor reduces the intensity of the Bragg peaks and the S(2θ) appears substantially unchanged, still keeping glassy-like smooth features. The crystallization process is summarized in Fig. 2(a), that shows a comparison between the first fully amorphous S(2θ) measured at Texp and one acquired after 60 h. The position of the Bragg peaks observed in the latter is compatible with that of the α-quartz22, Fig. 2(b), and no traces of rutile-like structure are visible21, Fig. 2(c). This also implies that the crystallization process preserves the chemical composition without any appreciable phase separation.

Figure 1: Time evolution of the static structure factor of GeO2 at Texp = 1100 K.
figure 1

Time increases from left to right. The color map shows the emergence of the crystalline pattern.

Figure 2
figure 2

(a) Static structure factor S(2θ) measured at t = 0.2 h (black open circles) and at t = 60.0 h (red diamonds). The solid lines are just connections between experimental points. Calculated Bragg peak positions for GeO2 crystalline polymorphs at room temperature: (b) α-quartz-like structure (P3221), ref. 22; (c) rutile-like tetragonal structure (P42/mnm), ref. 21.

Determination of crystalline and amorphous fractions

Assuming that no contributions arise from the crystal-amorphous interfaces, we can write the measured S(2θ) as the sum of an amorphous term and a crystalline one, namely

The first term accounts for the amorphous fraction of the material, and

where Aa is a parameter and Sg is the static structure factor of the fully amorphous system, obtained by considering the first scans at 1100 K, where no trace of crystallization is visible. Since the measurement of the static structure factor implies an integral over the energy, the crystalline term SC(2θ) can be written as the sum of an elastic contribution SB(2θ) accounting for the Bragg peaks, and an inelastic part identified as the thermal diffuse scattering (TDS). Consequently,

where Ac is a parameter. Following ref. 26, we use a very simple approximation for the TDS contribution ST(2θ), i.e.

where is the Debye Waller factor. The parameter B for v-GeO2 at 1100 K was calculated in harmonic approximation using the vibrational density of states from ref. 27. Finally, SB(2θ) is modelled describing each Bragg peak with a Gaussian, whose position 2θi is given by the α-quartz structure using the appropriate lattice parameters a and c. A preliminary analysis of the Bragg diffraction pattern did not show any appreciable t-evolution of the peak width, so we assumed that the peak full width half maximum is given by the instrument resolution, see supplementary info. The Bragg contribution turns out to be:

where Ai is the integrated intensity of the ith reflection which is fitted independently for each peak.

The lattice parameters a and c for the α-quartz GeO2 at Texp were determined by fitting the most crystallized data with Eq. 1. This leads to a(Texp) = 5.053 ± 0.002 Å and c(Texp) = 5.66 ± 0.04 Å, that were then fixed to fit the whole t-evolution. As compared to their room temperature counterparts a(RT) = 4.987 Å and c(RT) = 5.652 Å22, the high-T values are slightly dilated, and the thermal dilatation seems fairly anisotropic, as it affects a more than c.

Equation 1 provides a good fit to data during the whole observed process. This is visible in Fig. 3, where two examples at two different times are reported: Fig. 3(a) shows the early stage of the crystallization, t = 9.5 h, while Fig. 3(b) reports the result after 60 h.

Figure 3
figure 3

Static structure factor of GeO2 at T = 1100 K, measured at t = 9.5 h (a) and t = 60.0 h (b). The red line is the best fit to the data (black circles) according to Eq. 1. The solid blue line is ST(2θ), and the solid green line shows the sum of the amorphous and TDS components, see text.

The time evolution of the integrated intensity of the first two Bragg reflections is reported in Fig. 4. The intensity of the reflections shows a smooth increase as a function of time with an inflection point after about 20 h and a tendency to a long time saturation. However, the (101) reflection displays a step like increase at about t = 36 h, hardly visible in the (100). To trace the origin of this feature, we can analyse the intensity collected on the PSD detector, which is shown at three different times in Fig. 4(b,c,d). A typical amorphous pattern, with broad and regular Debye-Scherrer rings, is visible in the upper panel, corresponding to the beginning of the isotherm. After about 20 h, Bragg peaks emerge and the intensity becomes polycrystal-like, with some high intensity spots located on the Debye-Scherrer rings. This suggests that most part of the crystalline phase is basically a polycrystal, i.e. a spherically averaged assembly of small crystalline domains. However, some domains can grow more than the average and, if conveniently oriented, they produce the observed single crystal diffraction, with Bragg spots on the Debye-Scherrer rings. These bigger grains are then modified by the growth of neighbouring crystals, so the spots can evolve and eventually disappear because of the orientation change of the corresponding crystallite. Consequently, it is clear that the sample cannot be considered as a perfect powder and the fraction of crystallized material is not directly accessible by looking at the intensity of Bragg peaks. However, we can extract the fraction of atoms in the crystalline and amorphous phase by resorting to the coefficients Ac and Aa of Eqs 2 and 3. As a matter of fact, the scattering intensity at high scattering angle is proportional to the number of atoms and the properly normalized S(2θ) is equal to 1. In this limit, Bragg peaks are suppressed by the Debye-Waller factor and smeared out by the instrument resolution, whereas the TDS and the amorphous static structure factor go to unity. This means that Eq. 1 reduces to S(2θ) = Aa + Ac ≃ 1, hence Ac represents the fraction of atoms in the crystalline phase, whereas Aa represents that in the amorphous one. The time evolution of these quantities is reported in Fig. 5 and provides an insight into the kinetics of the crystallization process, as well as into the corresponding decrease of the amorphous matrix. In particular, a qualitative analysis of their shapes confirms that crystallization becomes appreciable after 4 h, and then it rapidly develops by subtracting material from the amorphous phase. After about 30 h, the crystallization rate slows down leading to a final material where a 77% of the atoms is organized in the α-quartz structure, whereas the remaining 23% still shows amorphous features.

Figure 4
figure 4

(a) Time evolution of integrated intensity Ai for the first two reflections of the α-quartz structure, (100) and (101), see legend. (b,c,d) Intensity measured on a portion of the detector at three different times during the isothermal process. The color scale of the PSD maps is the same for all the images.

Figure 5: Time evolution of the crystalline and amorphous fractions, Ac and Aa, open red diamonds and open green circles respectively.
figure 5

The solid line represents the fit with the model, as described in the text.

Discussion

The standard framework to describe the time evolution of the fraction of transformed material during isothermal crystallization is the Johnson-Mehl-Avrami-Kolmogorov (JMAK) model28,29,30,31,32. This approach is based on the nucleation and growth processes, and it assumes that the nucleation occurs randomly with a large number of spherical transforming regions. Growth is the same for all these regions and it stops at points of impingement, continuing elsewhere. This results in a complete crystallization of the starting material, which is in contrast with the results shown in Fig. 5. Consequently, in the present case, the JMAK model fails in describing the observed long-time behaviour. A complete analysis with the JMAK equation is reported in the supplementary info.

Consequently, all the experimental observations need to be combined to develop an empirical model for the kinematic of this crystallization process allowing a non-complete transformation of one phase to the other. This can be done by taking into account the density difference between the amorphous phase and the crystal. We thus consider a system of Ntot atoms at a temperature T ≪ Tm, i.e. where diffusion can be safely neglected. At a given time t the system is composed by Nc atoms in the crystalline phase and Na atoms in the amorphous one. Of course, the relation Nc(t) + Na(t) = Ntot holds at any time. Data suggest the presence of many different crystalline domains evolving in time. According to this, Nc can be written as

where N(t) is the number of crystalline nuclei at a given time t, whereas ni is the number of atoms in the ith domain that originates at ti. Considering that the process is almost continuous and both N(t) ≫ 1 and ni(t) ≫ 1, we can write the sum of Eq. 6 as the time integral

The function n(t) is zero when t ≤ 0, therefore it is convenient to change the integration variable to τ = t − t′, so that Eq. 7 is rewritten as:

where we use the compact notation N′(t) = dN(t)/dt.

Equation 8 describes a process that develops through a nucleation and growth mechanism. The creation of a nucleus is assumed to be a stochastic process that can be triggered by thermal fluctuations in the material and probably eased by even intrinsic heterogeneities10,11. Conversely, since diffusion is practically arrested, the growth of each nucleus proceeds only through structural rearrangements involving the interface between the ordered and disordered regions in a self-limiting process. Indeed, each crystalline domain nucleates and grows at expenses of the surrounding amorphous region. However, since the crystal has a higher density than the supercooled liquid, this mechanism creates high-density fully ordered regions that become surrounded by depleted interfaces. In absence of diffusion, this process slows down and stops the growth of the crystalline nuclei. A similar mechanism applies also to nucleation, that becomes less probable in overcrowded and depleted environments.

According to the previous considerations, the equation for evolution of the number of nuclei N(t) is assumed to be related to a constant rate and can be written as:

where 1/τn is the rate of the nucleation process and Nm is the maximum number of nuclei.

Conversely, for the growth processes, we can write that

where α is a growth parameter and nm is the maximum number of atoms in each nucleus. Of course nm can vary from nucleus to nucleus but, for simplicity, it is assumed to be constant throughout the sample and independent of time. The exponent μ accounts for the fraction of atoms involved in the process, hence in general 0 ≤ μ < 1. Actually, μ conveys information on the geometry of the nucleus and on the dimensionality d of the process, being its minimum value 1 − 1/d. In the case of spherical nuclei, d = 3 and μ turns out to be 2/3. This value seems adequate in the present case.

By changing the variables and considering n* = n/nm and t* = t/τd, with , Eq. 10 turns out to be

This scaled equation does not depend on the actual values of nm and α, and it can be integrated numerically. Typical results for different values of μ are shown in Fig. 6.

Figure 6
figure 6

Solution of Eq. 11 for different values of the exponent μ.

Finally, in a real experimental case, due to the previous thermal history of the sample, we have to consider the possibility that crystallization has already started at t = 0. This can be easily incorporated into the model by integrating Eq. 9 with the condition N(t = 0) = N0. Conversely, for the growth, we can use Eq. 11 by introducing a fictive time t0 such as

Within these assumptions, Eq. 7 can be rewritten as

where is the initial number of crystallized atoms, being n(t = 0) = n0.

Once divided by Ntot, Eq. 13 allows for the fit of the crystalline fraction Ac shown in Fig. 5. Moreover, since Ntot = Na(t) + Nc(t), the amorphous fraction Aa can be also analyzed by considering 1 − Nc(t)/Ntot. The model has thus five fitting parameters: the initial and the final fraction of atoms in the crystalline phase, and , the timescales of nucleation and growth, τn τd, and the fictive starting time t0. The parameters can be fitted to Ac and Aa using an overall procedure. The result is reported in Fig. 5. The model provides an accurate and coherent description of both the crystal growth and the amorphous reduction. At the beginning of the process, the system is almost fully amorphous and the crystalline fraction is about while . However, the presence of a nonzero t0 indicates that the crystallization process is already active, since it starts with a nonzero derivative. The resulting timescales for the nucleation and growth are τn = 16.9 ± 0.2 h and τd = 31.3 ± 0.6 h.

Starting from the density ρm of the material resulting from the experiment, we can evaluate the average density of the residual amorphous medium ρa at room temperature. The density was measured with a pycnometer and resulted ρm = (3.96 ± 0.01) g/cm3. Considering and with density ρc = 4.25 g/cm3, the average density of the amorphous part turns out to be ρa = (3.25 ± 0.03) g/cm3, about 91% of the room temperature glassy value33, ρg = 3.66 g/cm3. Of course ρa is an average value considering both the strongly depleted interfaces and the glassy-like regions. Assuming that the interfaces are fully depleted, we can estimate an upper limit for the true glassy regions as about 89% of the final amorphous material.

Conclusions

We performed a neutron diffraction study of the isothermal crystallization kinetics in the deep supercooled liquid phase of the strong glassformer GeO2. In the experimental time window, the amorphous system evolves towards a material populated by crystalline α-quartz like nuclei plunged into a residual amorphous medium. Due to the difference in density between the two phases, the growth of a crystal causes a rarefaction in the amorphous medium around it. In absence of diffusion, these depleted interfaces act as a barrier for a further growth. The same mechanism involves also the nucleation, that becomes less probable as the nuclei population increases. The density difference between the amorphous and crystalline phases can be considered as a feedback mechanism that controls the growth of the crystal. The so-established predator-prey equilibrium stops further growth at a given crystallite size. In particular, this process introduces a slowing down of the growth function at short time with respect to the simple case, in which diffusion provides material to the nucleating phase. This short time behavior is already evident from the S(2θ) data of Fig. 1, and it is well described by the model, as shown in Fig. 5.

The model can be further extended accounting also for the effects of the diffusion. As a matter of fact, even if its contribution seems negligible in the observed timescale, diffusion could give rise to long-time contributions that are expected to become dominant as Tm is approached.

Materials and Methods

The experiment was carried out at the two axis spectrometer 7C234. This instrument is located on the hot source of the reactor Orphée at the Laboratoire Leon Brillouin (CEA Saclay, France), and it is optimized for structural studies of liquids and amorphous systems. The Cu (111) monochromator of the instrument was set to obtain an incident neutron wavelength λ = 0.724 Å. This value, coupled with the high molecular mass of the sample, is short enough to make the inelastic correction negligible. Scattered neutrons were collected on the recently installed position sensitive detector (PSD) based on an assembly of 256 3He tubes. The whole detector covers a scattering angle of 128°, and each tube has 64 vertical pixels 8 mm height. The calibration of the instrument and the incident wavelength were checked by acquiring the diffraction pattern of a Ni sample. In addition, we measured also a second crystalline standard, a KBr powder, that allowed a thoroughly determination of the resolution function even at low scattering angles, see Supplementary Info.

Vitreous germania samples were prepared by melt-quenching, starting from Aldrich germanium (IV) oxide crystalline powder (purity higher than 99.998%). The powder was melted in an alumina crucible at about 1900 K and then quenched in air. Cylindrical specimens with a diameter of 8 mm were core-drilled from the bulk glass, and piled-up to obtain a 50 mm high sample. The sample was loaded in a vanadium cylindrical cell with an outer diameter of 10 mm and 0.5 mm thick walls. The cell was closed with a steel screw cap, which was carefully shielded with boron nitride masks to minimize its scattering contribution. High temperature measurements were done using a vanadium oven under vacuum, ~10−6 mbar. The temperature was monitored by two thermocouples fixed on the body of the sample holder.

In order to properly evaluate the single scattering intensity of the sample, we collected a set of ancillary measurements including the empty cell, an absorber (a cadmium rod with the same size as the sample container), and the empty beam35,36. The empty cell was also used as vanadium standard to normalize data. The intensity measured on the PSD detector was reduced to I(2θ) using the program ScRiPT provided by LLB. Starting from these data, the properly normalized static structure factor S(2θ) was determined via the procedure outlined in refs 35 and 36. Monte Carlo simulations were exploited to properly estimate the transmission coefficients and the multiple scattering contribution.

Additional Information

How to cite this article: Zanatta, M. et al. Real-time observation of the isothermal crystallization kinetics in a deeply supercooled liquid. Sci. Rep. 7, 43671; doi: 10.1038/srep43671 (2017).

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.