Introduction

Pressure (P) or temperature (T) induced modifications in crystal structures and associated effects on the lattice dynamics are commonly observed and reasonably well understood. On the contrary, transformations in amorphous systems between distinct aggregates having different local structure and density are more elusive. These polyamorphic (PA) transitions are often difficult to observe, since hampered by several concomitant factors. For instance, when the density is the order parameter, extreme thermodynamic conditions are required to significantly alter this variable due to the low compressibility of amorphous, non-gaseous, systems. Furthermore, PA phenomena often happen in metastable thermodynamic regions, where they are overshadowed by competing effects, such as glass transition or crystal nucleation.

On a general ground, the best candidates to observe PA transitions are systems with an intrinsically open, often tetrahedral, local structure. In fact, the large free volume available in tetrahedral arrangements can in principle allow structural modifications even at moderate thermodynamic conditions. This was demonstrated to be the case in water1, liquid silicon2, germanium3 and phosphorus4, as well as in amorphous SiO25 and GeO26 (see ref. 7 for a review on the topic). In spite of a thorough experimental scrutiny, some general aspects of PA transitions are still obscure, including the possible influence on the propagation of collective excitations. This can be particularly relevant at mesoscopic (~nm) length-scales, where the dynamics is known to be strongly coupled with local atomic arrangements.

Inelastic neutron (INS) and x-ray (IXS) scattering are two classic experimental techniques commonly used to probe atomic and lattice motions; however severe technical difficulties hinder the observation of PA transitions. A major one relates to the fact that modifications in the local order involving, e.g. the coordination number, usually disappear when the sample is recovered to ambient conditions (this doesn’t apply to long-range “densification” effects, which in glasses are usually permanent). This imposes in situ high-pressure experiments on very small samples for direct observation of PA transitions. This, in most practical situations, rules out the possibility of using INS and often causes problems of spectral background in IXS measurements. Although water8 and silica9 are the two polyamorphic materials that have been most extensively investigated by inelastic spectroscopies, no signature of PA transitions has been reported in the THz spectrum in either case, mainly owing to two different reasons: in silica PA phenomena happen at pressures still prohibitively high for scattering measurements, while in water the PA transition is expected to take place in a deeply supercooled region, representing a sort of no man’s land in the thermodynamic plane10.

Compared to SiO2, its structural analogous GeO2 has proven to be a better candidate for IXS investigations of PA phenomena due to both the larger tetrahedral cell, which shifts the onset of PA transitions to lower P’s and the higher electronic number—and consequently shorter x-ray absorption length—which substantially enhances the IXS signal from the small-sized sample suited for the use of Diamond Anvil Cells (DAC)11. Accordingly, we have studied the pressure-dependent spectrum of density fluctuation, S(Q, ω), of vitreous (v-)GeO2 at ambient temperature by in situ IXS measurements from ambient P up to 26 GPa (see Methods for further details). This P range has been chosen to well track the P-dependence below and above 9 GPa, pressure around which a sudden jump of the bond distance is reported and commonly ascribed to a transition from a tetrahedral to an octahedral local structure6,12, or, in other terms, from an α-quartz-like to a rutile-like local lattice organization. In a more recent x-ray diffraction work13 important structural changes have been observed to continuously occur for pressure spanning the 5–8.6 GPa range, as later confirmed by oxygen K-edge IXS measurements14. Furthermore, previous studies based on classical molecular dynamics15 predicted a main structural change from tetrahedral to octahedral arrangement at 3–7 GPa, slightly lower than observed here and reported in previous works.

In addition to diffraction and absorption measurements16, signatures of a PA transition in v-GeO2 have been sought for by investigating the vibrational behavior by Raman and infrared techniques16,17, which provided evidences of possible dynamic counterparts of the aforementioned PA transition. Combining experimental results with first-principle density functional theory (DFT) calculations, we aim at detecting and explaining the signature of the PA crossover in the THZ spectrum of v-GeO2.

Results

Typical IXS spectra measured at ambient pressure and at P = 26 GPa, i.e. respectively well below and above the PA transition, are shown in Fig. 1 along with the experimental resolution function. We recall here that PA transition of GeO2 is here identified with the large transformation of the average Ge-O bond distance, dO-Ge, in the range(see refs 6,12). Although, as mentioned, such a crossover was also either observed as a sharp, yet continuous, P-increase13,14 or to occur at slightly lower P values, it is unanimously found that it takes place for P < 10 GPa. A clear excess of scattering intensity with respect to the resolution profile can be observed in all spectra, indicating the presence of an appreciable contribution of collective excitations to the experimental spectra in all the probed (Q, P)-values. Such inelastic spectral wings turn into a “double shoulder” line-shape at higher Q’s for both P values, hence suggesting the presence of two distinct collective modes. In order to determine the Q-dispersion relations of such modes and gain insights into their correlations to structural changes within the probed P-range, we performed a best-fit line-shape analysis to determine the characteristic frequencies of the high (ΩHF) and low frequency (ΩLF) excitations (see Methods for further details).

Figure 1
figure 1

IXS spectra below and above the PA crossover.

Representative IXS spectra of v-GeO2 measured at low and high pressures for selected Q values (open circles). The thick black line corresponds to the best-fit model line-shape in Eq. 1, with its low frequency (blue line) and high frequency (red line) DHO components. The dashed and cyan lines represent respectively the resolution function and a DHO profile accounting for the transverse mode of the diamond anvils.

The dispersion curves obtained as the best-fit values of the ΩLF and ΩHF parameters in Eq. 1 are displayed in Fig. 2. It may be noticed that at the lowest Q’s some ΩLF are not reported in the plot; this owes to the corresponding vanishing intensity of the low frequency peak in the spectrum.

Figure 2
figure 2

Dispersion curves below and above the PA crossover.

Best-fit dispersion of ΩLF (blue symbols) and ΩHF (red symbols) as measured at the indicated pressures below (left panel) and above (right panel) the PA transition.

Even more striking is the transformation of the P-evolution of the generalized sound velocities cs (Fig. 3) extracted from the low Q (≤10 nm−1) slope of the dispersions curves in Fig. 2. The fact that IXS probes sound propagation at THz frequencies and nm distances makes IXS measurements more sensitive to local molecular arrangements and interactions than ultrasound or Brillouin light scattering (BLS) techniques. Derived velocities can thus differ from ultrasonic and hypersonic determinations (see for instance the discussion in ref. 11), bringing complementary information. To further stress the distinctive elastic behavior of the low-pressure and high-pressure polyamorph, the inset of Fig. 3 shows the value of the longitudinal modulus M = ρcs2 - with ρ being the mass density, here derived from from ref. 12.

Figure 3
figure 3

P-dependence of the sound velocity across the PA crossover.

The sound velocities extracted from the low Q slope of the dispersion curves in Fig. 2 (see text) are reported (blue squares) as a function of pressure. The dashed lines are the outcome of a linear fit to data below and above the PA crossover. The pressure value at which the two straight lines intersect (~8.7 GPa) is highlighted by a vertical dotted line. The red dots are the velocities extracted from ref. 11. The inset shows the corresponding longitudinal modulus M, as derived using density data in ref. 12.

Figures 4 and 5 compare, respectively, the lowest P (P = 0 GPa) and highest P (P = 26 GPa) dispersions shown in Fig. 2 with the phonon dispersions computed for crystalline α-quartz GeO2 and rutile GeO2 at the corresponding Ps along high symmetry directions; here the thickness of the dispersion curves is proportional to the weight of the corresponding contributions to the S(Q, ω) (see Methods for further details). The green lines in Figs 4 and 5 denote the (ω, Q)-range dominated by the phonon modes of diamond. Within this region the very intense scattering contribution from the diamond anvils hampers the detection of the signal from the sample. Right panels in Figs 4 and 5 show the computed total vibrational density of state (v-DOS). The v-DOS measured by Raman scattering at P = 0 GPa18 and 32 GPa19 are also included for comparison in Figs 4 and 5; such profiles were derived from Raman measurements as prescribed by Eq. 2 of ref. 18, yet without correction for the light-vibration coupling factor.

Figure 4
figure 4

Comparison between experimental and theoretical results below the PA crossover.

Left: Dispersion curves along high-symmetry directions of α-quartz GeO2 computed at P = 0 GPa (black lines) and P = 0 GPa values of ΩLF (blue symbols) and ΩHF (red symbols) measured in v-GeO2. The thickness of computed dispersions is proportional to the corresponding S(Q, ω) contribution (see text). In the three plots the lowest frequency TA branch is also reported for comparison as dashed line. The green curves represent the boundary of the zone dominated by the diamond’s phonons. Right: computed v-DOS (red line) of α-quartz GeO2 at P = 0 GPa, compared with Raman scattering measurements in v-GeO2 at P = 0 GPa18 (dashed line), after re-scaling for an arbitrary factor. The inset displays the calculated v-DOS in an extended energy range.

Figure 5
figure 5

Comparison between experimental and theoretical results above the PA crossover.

Same as Fig. 4 for the rutile GeO2 at P = 26 GPa. The DFT computed dispersions are compared to the P = 26 GPa values of ΩLF (blue symbols) and ΩHF (red symbols) measured in v-GeO2 at P = 26 GPa. The corresponding v-DOS in the right plot is compared to the Raman scattering profile taken from ref. 19 at P = 32 GPa.

Discussion

The quality of the data analysis can be readily judged from the overall good agreement between measured and best-fit line-shapes in Fig. 1, where the spectral contributions of the low and high frequency modes to the total scattering intensity are also shown. By comparing inelastic modes in the spectra measured at P = 0 GPa and P = 26 GPa one immediately observes that both inelastic shift and relative intensity of the high frequency mode increase dramatically with pressure.

A substantial difference between low- and high-P sound dispersions readily emerges from Fig. 2: ΩLF and ΩHF at low P (left panel) span over lower energy values and exhibit a moderate Q-dependence. In contrast, the high-P dispersions (right panel) are featured by a strongly Q-dependent high-ω phonon branch and by a mildly Q-dependent low-ω branch. Such features become even more striking when considering that P effects are relatively moderate within each subset of data. This result is in good agreement/compatible with the expected transition value of 9 GPa.

Further transformations of the acoustic properties of the sample across the PA crossover clearly emerge in Fig. 3 through the cusp-like P dependence of cs.

The drastic change is emphasized by the different slopes of the straight lines best fitting all reported cs data either below or above the crossover (dashed lines). These lines intercept at a P value (≈8.7 GPa) close to the known PA crossover pressure. The lower high P slope clearly indicates a higher resistance of the sample to P-induced modifications of elastic properties. In a similar manner, in recent BLS works20,21 changes in the slope of the P-evolution of sound velocity and/or elastic moduli were interpreted as signatures of PA phenomena; more specifically, they have been correlated to changes in the coordination number. Here the cusp-like behavior exhibited by cs data seems more pronounced. This is the likely consequence of the mentioned rough matching between probed distances and first neighboring molecules’.

As the crystalline α-quartz and rutile structures are known to approximate first neighbor Ge-O bonding arrangements of v-GeO2 for P below and above the PA crossover6, we use DFT results as an interpretative guidance for the observed P-evolution of the collective dynamics measured on glasses. Bearing in mind the different sample nature (crystalline vs glassy), experimental and computational results shown in Figs 4 and 5 are overall consistent. In particular, it is worth noticing that both computed dispersions and v-DOS exhibit distinctive features, markedly different between the α-quartz and rutile phases. Indeed, dispersion curves of α-quartz clearly span over lower energies and optical phonon modes exhibit a classical relatively mild Q dependence. Accordingly, the v-DOS of this structure is characterized by the expected parabolic low-frequency trend, followed by van Hove singularities of comparable intensity and, finally, by a clear phonon gap in the 45–55 meV interval (see inset in Fig. 3). In contrast, for the the rutile lattice arrangement, the main contribution to S(Q, ω) comes from highly dispersive high-ω phonon branches. The v-DOS is dominated by a large van Hove singularity at about 45 meV and presents no evidence of phonon gaps below 60 meV. As a consequence, within the reported 0–60 meV range, the v-DOS of rutile spans comparatively higher ω’s than that of α-quartz. In view of these qualitative differences, the overall agreement of experimental data collected at pressures below 9 GPa with calculations for the α-quartz structure (Fig. 4) and of the experimental data collected at pressure above 9 GPa with the calculations for the rutile crystal structure (Fig. 5) provide strong indications that the PA quartz-like to rutile-like transition in the local structure of v-GeO26 has a well-defined counterpart in the phonon dispersion behavior.

A closer inspection of Fig. 5 reveals a few additional interesting aspects. The IXS measurements unmistakably show that, at high P, the S(Q, ω) of v-GeO2 is dominated by two modes with very different energy, one highly Q-dispersive and a second weakly Q-dispersive. The presence of the latter was particularly evident in the Spring-8 data (P = 13 GPa) most likely because of the better contrast in the low-frequency spectral range (see Methods). The comparison with the computed dispersion curves of the rutile crystal allows us to tentatively associate these two branches to longitudinal (LA) and transverse acoustic (TA) modes, respectively.

In order to facilitate the comparison between ΩLF values and the TA modes in the crystal, in Fig. 4 the lowest frequency TA branches of α-quartz are also indicated as dashed lines in the three crystalline paths considered. However, the spectral contribution of these branches was found negligible in the crystal. Whereas for the crystal phase TA contributions to S(Q, ω) are forbidden within the first Brillouin zone, they may become sizable in the glassy phase. This directly relates to the absence of long-range translational symmetries, or equivalently well-defined Brillouin zones. “Pure” symmetric TA modes can be “contaminated” and acquire a mixed longitudinal and transverse character. Such mode-mixing is often referred to as the longitudinal-transverse (L-T) coupling22 and is the physical rationale behind the appearance of a shear mode in the S(Q, ω), which primarily couples with longitudinal movements only. It is worth noticing that a similar L-T coupling has been reported in various systems sharing a tetrahedral molecular arrangements, such as water23,24,25,26, GeSe227 and v-GeO2 as well28. It can thus be envisaged that the open and highly directional nature of such arrangement fosters the onset of an L-T coupling.

Furthermore, the comparison between left and right panels in Fig. 5 indicates that the value of ΩLF well corresponds to the low frequency excess of the v-DOS in the glass as compared to the crystal. This intensity excess relates to almost universal feature of glasses essentially amounting in a peak in the reduced density of states v-GeO2/E2, customarily referred to as the Boson peak (BP).

The equivalence between the BP of a glass and the Van Hove singularity of the TA branch of the corresponding crystal has been recently demonstrated to be a sound hypothesis29,30. We observe that the intensity ratio between TA and LA modes decreases upon increasing the pressure consistently with the observed pressure trend of the BP31.

Conclusion

In summary, we investigated the evolution of the THz spectrum of density fluctuations in v-GeO2 as a function of pressure. We observed a clear transformation in the Q-dispersion behavior upon crossing the known polyamorphic transition occurring in the glassy phase at ≈9 GPa. Supported by DFT calculations, we interpreted this as the abrupt evolution from a quartz-like to a rutile-like behavior, concluding that the collective dynamics of v-GeO2 in the THz range is strongly sensitive to the undergoing changes in the local structure. This is clearly indicated by a cusp-like behavior of the pressure dependence of generalized sound velocity and longitudinal modulus across the PA transition. Both these trends suggest that the high-P polyamorphic phase is characterized by a higher resistance to pressure induced modification of elastic properties, likely due to the more packed first neighbor arrangement. Furthermore, presented data indicates that the inherent disorder characteristics of the glassy phase seems to foster the visibility of a low frequency transverse modes, which, especially at low pressures, is evident at high Q’s but it is still appreciable within the first pseudo-Brillouin zone. This mode in the glass is here found to be possibly related to the low frequency excess intensity in the vibrational density of state.

We finally remark that, when compared to more traditional structural measurements, investigations of the THz dynamics as the one presented in this work provide a complementary insight onto the PA transition linking it to transformations of elastic properties.

More in general, the results presented in this work pave the way toward future investigations of the dynamic aspects of polyamorphism phenomena by using THz probes of phonon-like modes.

Methods

IXS measurements and data analysis

Two independent IXS experiments were carried out on two different IXS beamlines: beamline BL35XU32 at SPring-8 (Hyogo, Japan) and beamline ID2833 at the European Synchrotron Radiation Facility (ESRF, Grenoble, France). These two triple-axis spectrometers have the same working principle, based upon high order reflections from nearly perfect Si crystals. The instruments were operated at the Si(9, 9, 9) configuration for both monochromator and analyzer crystals, corresponding to an incident energy of 17.947 KeV and an overall resolution bandwidth of 3.0 meV (FWHM). The profile of the instrumental resolution function, R(ω) was estimated through the measurement of the spectral line-shape of an essentially elastic scatterer, specifically a cryogenically cooled sample of perspex at the Q position of its first sharp diffraction peak (10 nm−1). The spectrometer was rotated to probe within a single energy scan the 4.5–7.8 nm−1 Q-range (BL35XU measurements) or the 5.4–13.8 nm−1 Q range (ID28 measurements). The focal spot of the beam at sample position was 35 × 15 μm2 and 30 × 60 μm2 (horizontal × vertical FWHM) for measurements at BL35XU and ID28, respectively.

The v-GeO2 sample was synthesized as described in a previous work34 (BL35XU experiments) or purchased (Sigma-Aldrich; ID28 experiments). In both cases, the samples were loaded in DAC without any pressure transmitting medium to maximize the scattering volume for high-pressure measurements. Samples used for BL35XU experiments were conditioned in Mao type DAC, equipped with tungsten gaskets, while samples used for ESRF experiments were conditioned in membrane driven Le Toullec type DAC, equipped with rhenium gasket. We stress that the metallic gasket, beside radially containing the sample, also acts as cleaning pinhole in the x-ray path. To further minimize spurious quasi-elastic scattering contributions, DACs were placed in a specifically designed vacuum chamber equipped with motorized exit slits. We indifferently used 350 or 300 μm flat culets diamond to cover the pressure range of interest. A small chip of ruby placed in the sample chamber served as pressure gauge35. In order to measure the spectral line-shape with the needed statistical accuracy, we performed energy transfer scans typically lasting 18–24 hours. Prior to each spectral acquisition we measured the Q-dependent elastic scattering intensity I(Q, ω = 0), systematically confirming the vitreous nature of the investigated phase.

The pressures probed in the SPing-8 experiment were 3.7, 13 and 20 GPa, while in the ESRF experiment we collected data at room pressure and 26 GPa, in both cases the sample was at ambient temperature. SPring-8 measurements were performed with a better spectral contrast in the low-frequency spectral region, likely due to the lower quasi-elastic scattering from the diamonds DAC windows, while in spectra collected at ESRF had substantially higher count-rate and, consequently, a better statistical accuracy.

The best-fit line-shape modeling of IXS spectra is based on the sum of two Damped Harmonic Oscillator (DHO) functions plus an essentially elastic central peak. Overall the used profile reads as

where is the detailed balance factor accounting for the frequency-dependent statistical population of the states, with kB being the Boltzmann constant and with h being the Planck constant. The parameters ΩLFHF) and ΓLFHF) represent the shift and the width the low (high) frequency inelastic mode, respectively, while A, ΙLF and ΙHF are frequency independent scaling factors. A Q-dependence is assumed implicitly for all parameters. The double DHO profile of Eq. 1 was used to fit the measured spectral line-shape, even though at the lowest Q’s best-fit values of ILF were negligible. Additional DHO terms were added to account for the inelastic scattering from the phonon modes of diamond. Due to the very high sound speed of diamond, such well resolved spectral features locate in the high-frequency side of the spectra and disperse out of the probed window at high Q’s (see Fig. 1).

Ultimately, the best fit of measured spectra (see Fig. 1) was achieved through the minimization a χ2 variable defined as the normalized distance between the experimental line shape and the following model profile:

The symbol “” represents the convolution operator, while K and B are two frequency-independent constants representing, respectively, an overall intensity factor and a flat background, which also includes the electronic noise of the detectors (<1 mHz).

Numerical computations

To get further insights into the dynamic response of the sample, we complemented our experimental investigation with a numerical study. DFT calculations within the local density approximations were carried out using the ABINIT package36 and norm-conserving pseudopotentials in the Troullier-Martins scheme. The crystal structures of bulk GeO2 were optimized using a kinetic energy cutoff of 58 Hatree and the first Brillouin zone were sampled using a 4 × 4 × 4 and 6 × 6 × 6 Q-grid for the α-quartz and rutile structures, respectively. Phonon dispersion curves were interpolated over a 4 × 4 × 4 and 3 × 3 × 3 Q-grid for α-quartz and rutile structures of GeO2, respectively, using the density-functional perturbation theory scheme described in ref. 37. Nonanalytic corrections due to the long-ranged anisotropic dipole-dipole interactions were applied to the dynamic matrix. The computed zone-center phonon frequencies at P = 0 GPa for α-quartz and rutile structures are within 0.6 meV of the experimental values38,39 and previous DFT results40. The weight of the individual contributions to S(Q, ω) at a discrete Qi on the j-th phonon branch (represented in Figs 3 and 4 as the thickness of the corresponding line after multiplying with ω(Qi)) is denoted as Sj(Qi):

with:

where fd(Q) is the x-ray form factor and is the projection on to atom d of eigenvector in the j-th branch. Effects due to the Debye-Waller factor were ignored since it induces a negligible Q dependence within the probed momentum transfer interval.

Additional Information

How to cite this article: Cunsolo, A. et al. Signature of a polyamorphic transition in the THz spectrum of vitreous GeO2. Sci. Rep. 5, 14996; doi: 10.1038/srep14996 (2015).