# Signature of a polyamorphic transition in the THz spectrum of vitreous GeO2

Alessandro Cunsolo, Yan Li, Chaminda N. Kodituwakku, Shibing Wang, Daniele Antonangeli, Filippo Bencivenga, Andrea Battistoni, Roberto Verbeni, Satoshi Tsutsui, Alfred Q. R. Baron, et al.

► To cite this version:

Alessandro Cunsolo, Yan Li, Chaminda N. Kodituwakku, Shibing Wang, Daniele Antonangeli, et al. Signature of a polyamorphic transition in the THz spectrum of vitreous GeO2. Scientific Reports, 2015, 5, pp.2015. 10.1038/srep14996 . hal-01216568

HAL Id: hal-01216568

https://hal.science/hal-01216568v1

Submitted on 6 Jan 2025

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

<figure><img src="image_2.png" /></figure>

OPEN

# Signature of a polyamorphic transition in the THz spectrum of vitreous GeO₂

Received: 09 July 2015

Accepted: 15 September 2015

Published: 13 October 2015

Alessandro Cunsolo<fnref n="1" />, Yan Li², Chaminda N. Kodituwakku<fnref n="1" />, Shibing Wang³, Daniele Antonangeli⁴, Filippo Bencivenga⁵, Andrea Battistoni⁵,⁶, Roberto Verbeni⁷, Satoshi Tsutsui⁸, Alfred Q. R. Baron⁸,⁹, Ho-Kwang Mao¹⁰,¹¹, Dima Bolmatov<fnref n="1" /> & Yong Q. Cai<fnref n="1" />

The THz spectrum of density fluctuations, S(Q, ω), of vitreous GeO₂ at ambient temperature was measured by inelastic x-ray scattering from ambient pressure up to pressures well beyond that of the known α-quartz to rutile polyamorphic (PA) transition. We observe significant differences in the spectral shape measured below and above the PA transition, in particular, in the 30–80 meV range. Guided by first-principle lattice dynamics calculations, we interpret the changes in the phonon dispersion as the evolution from a quartz-like to a rutile-like coordination. Notably, such a crossover is accompanied by a cusp-like behavior in the pressure dependence of the elastic response of the system. Overall, the presented results highlight the complex fingerprint of PA phenomena on the high-frequency phonon dispersion.

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 water<fnref n="1" />, liquid silicon², germanium³, and phosphorus⁴, as well as in amorphous SiO₂⁵ and GeO₂⁶ (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.

<footnote n="1">¹National Synchrotron Light Source II, Brookhaven National Laboratory, Upton, NY 11973, USA. ²American Physical Society, 1 Research Road, Ridge, New York 11961, USA. ³Department of Geological and Environmental Sciences, Stanford University, Stanford, CA 94305, USA. ⁴Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, UMR CNRS 7590, Sorbonne Universités - UPMC, Muséum National d'Historie Naturelle, IRD Unité 206, 75252 Paris, France. ⁵Sincrotrone Trieste, S.S. 14 km 163,5 in AREA Science Park 34012 Basovizza, Trieste, Italy. ⁶Dipartimento di Fisica, Università degli Studi di Trieste, I-34127 Trieste, Italy. ⁷European Synchrotron Radiation Facility (ESRF), 71 avenue des Martyrs, 38043 Grenoble, France. ⁸Japan Synchrotron Radiation Research Institute, SPring-8, Sayo, Hyogo 679-5198, Japan. ⁹SPring-8/RIKEN, Hyogo 679-5148, Japan. ¹⁰Geophysical Laboratory Carnegie Institution of Washington, 5251 Broad Branch Road, NW, Washington, DC 20015, USA. ¹¹Center for High Pressure Science & Technology Advanced Research, Pudong, Shanghai, 201203, China. Correspondence and requests for materials should be addressed to A.C. (email: acunsolo@bnl.gov)</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 1.** IXS spectra below and above the PA crossover. Representative IXS spectra of v-GeO₂ 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.</figcaption></figure>

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 water⁸ and silica⁹ 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 plane¹⁰.

Compared to SiO₂, its structural analogous GeO₂ 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)¹¹. Accordingly, we have studied the pressure-dependent spectrum of density fluctuation, *S*(*Q*, *ω*), of vitreous (v-)GeO₂ 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 structure⁶,¹², or, in other terms, from an α-quartz-like to a rutile-like local lattice organization. In a more recent x-ray diffraction work¹³ 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 measurements¹⁴. Furthermore, previous studies based on classical molecular dynamics¹⁵ 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 measurements¹⁶, signatures of a PA transition in v-GeO₂ have been sought for by investigating the vibrational behavior by Raman and infrared techniques¹⁶,¹⁷, 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-GeO₂.

## 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 GeO₂ is here identified with the large transformation of the average Ge-O bond
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Dispersion curves below and above the PA crossover. Best-fit dispersion of Ω<sub>LF</sub> (blue symbols) and Ω<sub>HF</sub> (red symbols) as measured at the indicated pressures below (left panel) and above (right panel) the PA transition.</figcaption></figure>

distance, $d_{O-Ge}$, in the $5 \lesssim P \lesssim 9$ GPa range(see refs 6,12). Although, as mentioned, such a crossover was also either observed as a sharp, yet continuous, P-increase<sup>13,14</sup> 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 ($\Omega_{HF}$) and low frequency ($\Omega_{LF}$) excitations (see Methods for further details).

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

Even more striking is the transformation of the P-evolution of the generalized sound velocities $c_s$ (Fig. 3) extracted from the low Q ($\le 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 = \rho c_s^2$ - with $\rho$ being the mass density, here derived from Eq. 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 $GeO_2$ and rutile $GeO_2$ 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, \omega)$ (see Methods for further details). The green lines in Figs 4 and 5 denote the ($\omega, 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$ GPa<sup>18</sup> and $32$ GPa<sup>19</sup> 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.

## 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: $\Omega_{LF}$ and $\Omega_{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-$\omega$ phonon
<figure><img src="image_5.png" /><figcaption>**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.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 4.** Comparison between experimental and theoretical results below the PA crossover. Left: Dispersion curves along high-symmetry directions of α-quartz GeO₂ computed at *P*=0 GPa (black lines) and *P*=0 GPa values of Ω<sub>LF</sub> (blue symbols) and Ω<sub>HF</sub> (red symbols) measured in v-GeO₂. 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 GeO₂ at *P*=0 GPa, compared with Raman scattering measurements in v-GeO₂ at *P*=0 GPa¹⁸ (dashed line), after re-scaling for an arbitrary factor. The inset displays the calculated v-DOS in an extended energy range.</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Comparison between experimental and theoretical results above the PA crossover. Same as Fig. 4 for the rutile $GeO_2$ at $P=26$ GPa. The DFT computed dispersions are compared to the $P=26$ GPa values of $\Omega_{LF}$ (blue symbols) and $\Omega_{HF}$ (red symbols) measured in v-GeO$_2$ 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.</figcaption></figure>

branch and by a mildly Q-dependent low-$\omega$ 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 $c_s$.

The drastic change is emphasized by the different slopes of the straight lines best fitting all reported $c_s$ 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 works$^{20,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 $c_s$ data seems more pronounced. This is the likely consequence of the mentioned rough matching between probed distances and first neighboring molecules$^1$.

As the crystalline α-quartz and rutile structures are known to approximate first neighbor Ge-O bonding arrangements of v-GeO$_2$ for P below and above the PA crossover$^6$, 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 rutile lattice arrangement, the main contribution to $S(Q, \omega)$ comes from highly dispersive high-$\omega$ 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 $\omega$'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-GeO$_2$$^6$ 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-GeO₂ 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 Ω<sub>LF</sub> 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) coupling<sup>22</sup> 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 water<sup>23-26</sup>, GeSe<sub>2</sub><sup>27</sup> and v-GeO<sub>2</sub> as well<sup>28</sup>. 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 Ω<sub>LF</sub> 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-GeO<sub>2</sub>/E<sup>2</sup>, 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 hypothesis<sup>29,30</sup>. We observe that the intensity ratio between TA and LA modes decreases upon increasing the pressure consistently with the observed pressure trend of the BP<sup>31</sup>.

## Conclusion

In summary, we investigated the evolution of the THz spectrum of density fluctuations in v-GeO₂ 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-GeO₂ 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 BL35XU<sup>32</sup> at SPring-8 (Hyogo, Japan) and beamline ID28<sup>33</sup> 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<sup>-1</sup>). The spectrometer was rotated to probe within a single energy scan the 4.5-7.8 nm<sup>-1</sup> Q-range (BL35XU measurements) or the 5.4-13.8 nm<sup>-1</sup> Q range (ID28 measurements). The focal spot of the beam at sample position was 35 × 15 µm<sup>2</sup> and 30 × 60 µm<sup>2</sup> (horizontal × vertical FWHM) for measurements at BL35XU and ID28, respectively.

The v-GeO₂ sample was synthesized as described in a previous work<sup>34</sup> (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 gauge³⁵. 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 exper-
iment we collected data at room pressure and 26 GPa, in both cases the sample was at ambient tem-
perature. 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

$$
\frac{S(Q, \omega)}{S(Q)} = f(\omega) \left[ A \delta(\omega) + I_{LF} \frac{(2\Omega_{LF} \Gamma_{LF})^2}{(\omega^2 - \Omega_{LF}^2)^2 + 4\Gamma_{LF}^2 \omega^2} + I_{HF} \frac{(2\Omega_{HF} \Gamma_{HF})^2}{(\omega^2 - \Omega_{HF}^2)^2 + 4\Gamma_{HF}^2 \omega^2} \right], \quad (1)
$$

where $f(\omega) = \hbar\omega/k_B T [1 - \exp(-\hbar\omega/k_B T)]$ is the detailed balance factor accounting for the frequency-dependent statistical population of the states, with $k_B$ being the Boltzmann constant and $\hbar = h/2\pi$ with $h$ being the Planck constant. The parameters $\Omega_{LF}(\Omega_{HF})$ and $\Gamma_{LF}(\Gamma_{HF})$ represent the shift and the width the low (high) frequency inelastic mode, respectively, while $A$, $I_{LF}$ and $I_{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 $I_{LF}$ 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
$\chi^2$ variable defined as the normalized distance between the experimental line shape and the following
model profile:

$$
I(Q, \omega) = KS(Q, \omega) \otimes R(\omega) + B, \quad (2)
$$

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 com-
plemented our experimental investigation with a numerical study. DFT calculations within the local den-
sity approximations were carried out using the ABINIT package³⁶ and norm-conserving pseudopotentials
in the Troullier-Martins scheme. The crystal structures of bulk GeO₂ 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 GeO₂, 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 values³⁸,³⁹ and previous DFT results⁴⁰. The weight of the individual contributions to
S(Q, ω) at a discrete Qᵢ on the j-th phonon branch (represented in Figs 3 and 4 as the thickness of the
corresponding line after multiplying with ω(Qᵢ)) is denoted as Sⱼ(Qᵢ):

$$
S(Q, \omega) = \sum_{i,j} S_j(Q_i) \delta(\omega - \omega_j(Q_i)), \quad (3)
$$

with:

$$
S_j(\mathbf{Q}_i) = \frac{1 + n_j(\mathbf{Q}_i)}{2\omega_j(\mathbf{Q}_i)} \left| \sum_d \frac{f_d(\mathbf{Q}_i)}{\sqrt{M_d}} \mathbf{Q}_i \cdot \sigma_d^j(\mathbf{Q}_i) \exp^{i\mathbf{Q}_i \cdot \mathbf{d}} \right|^2 \quad (4)
$$
where $f_d(Q)$ is the x-ray form factor and $\sigma_d^j(Q)$ 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.

## References

1. Mishima, O. & Stanley, H. E. The relationship between liquid, supercooled and glassy water. *Nature* **396**, 329–335 (1998).

2. Beye, M., Sorgenfrei, F., Schlotter, W. F., Wurth, W. & Föhlisch, A. The liquid-liquid phase transition in silicon revealed by snapshots of valence electrons. *Proc. Natl. Acad. Sci. USA* **107**, 16772–16776 (2010).

3. Köga, J., Okumura, H., Nishio, K., Yamaguchi, T. & Yonezawa, F. Simulational analysis of the local structure in liquid germanium under pressure. *Phys. Rev. B* **66**, 064211 (2002).

4. Katayama, Y. et al. A first-order liquid-liquid phase transition in phosphorus. *Nature* **403**, 170–173 (2000).

5. Lacks, D. J. First-Order Amorphous-Amorphous Transformation in Silica. *Phys. Rev. Lett.* **84**, 4629–4632 (2000).

6. Itie, J. P. et al. Pressure-induced coordination changes in crystalline and vitreous GeO₂. *Phys. Rev. Lett.* **63**, 398–401 (1989).

7. Brazhkin, V. V. & Lyapin, A. G. High-pressure phase transformations in liquids and amorphous solids. *J. Phys. Condens. Matter* **15**, 6059–6084 (2003).

8. Cunsolo, A. The THz Spectrum of Density Fluctuations of Water: The Viscoelastic Regime. *Advances in Condensed Matter Physics* **2015**, 137435 (2015).

9. Baldi, G., Giordano, V. M., Monaco, G. & Ruta, B. Sound attenuation at terahertz frequencies and the boson peak of vitreous silica. *Phys. Rev. Lett.* **104**, 195501 (2010).

10. Poole, P. H., Grande, T., Angell, C. A. & McMillan, P. F. Polymorphic phase transitions in liquids and glasses. *Science* **275**, 322–323 (1997).

11. Bencivenga, F. & Antonangeli, D. Positive sound dispersion in vitreous GeO₂ at high pressure. *Phys. Rev. B* **90**, 134310 (2014).

12. Smith, K. H., Shero, E., Chizmeshya, A. & Wolf, G. H. The equation of state of polyamorphic germania glass: A two-domain description of the viscoelastic response. *J. Chem. Phys.* **102**, 6851–6851 (1995).

13. Drewitt, J. W. E. et al. Structure of GeO₂ glass at pressures up to 8.6 GPa. *Phys. Rev. B* **81**, 014202 (2010).

14. Lelong, G. et al. Evidence of fivefold-coordinated Ge atoms in amorphous GeO₂ under pressure using inelastic x-ray scattering. *Phys Rev B* **85**, 134202 (2012).

15. Peralta, J. & Gutiérrez, G. Pressure-induced structural transition in amorphous GeO₂: a molecular dynamics simulation. *Eur. Phys. J. B* **87**, 1–9 (2014).

16. Micoulaut, M., Cormier, L. & Henderson, G. S. The structure of amorphous, crystalline and liquid GeO₂. *J. Phys: Condens. Matter* **18**, R753–R784 (2006).

17. Galeener, F. L., Geissberger, A. E., Ogar Jr., G. W. & Loehman, R. E. Vibrational dynamics in isotopically substituted vitreous GeO₂. *Phys. Rev. B* **28**, 4768–4773 (1983).

18. Deschamps, T., Martinet, C., de Ligny, D., Bruneel, J. L. & Champagnon, B. Correlation between boson peak and anomalous elastic behavior in GeO₂ glass: An *in situ* Raman scattering study under high-pressure. *J. Chem. Phys.* **134**, 234503 (2011).

19. Durben, D. J. & Wolf, G. H. Raman spectroscopic study of the pressure-induced coordination change in GeO₂ glass. *Phys. Rev. B* **43**, 2355–2363 (1991).

20. Murakami, M. & Bass, J. D. Spectroscopic Evidence for Ultrahigh-Pressure Polymorphism in SiO₂ Glass. *Phys. Rev. Lett.* **104**, 025504 (2010).

21. Huang, L., Nicholas, J., Kieffer, J. & Bass, J. D. Polyamorphic transitions in vitreous B₂O₃ under pressure. *J. Phys. Condens. Matter* **20**, 075107 (2008).

22. Sampoli, M., Ruocco, G. & Sette, F. Mixing of longitudinal and transverse dynamics in liquid water. *Phys. Rev. Lett.* **79**, 1678–1681 (1997).

23. Ruocco, G. et al. Equivalence of the sound velocity in water and ice at mesoscopic wavelengths. *Nature (London)* **379**, 521–523 (1996).

24. Pontecorvo, E. et al. G. High-frequency longitudinal and transverse dynamics in water. *Phys. Rev. E* **71**, 011501 (2005).

25. Cimatoribus, A. et al. C. The mixed longitudinal-transverse nature of collective modes in water. *New J. Phys.* **12**, 053008 (2010).

26. Cunsolo, A. et al. Transverse dynamics of water across the melting point: A parallel neutron and x-ray inelastic scattering study. *Phys. Rev. B* **85**, 174305 (2012).

27. Orsinger, L. et al. F. High-frequency dynamics of vitreous GeSe₂. *Phys. Rev. B* **82**, 115201 (2010).

28. Bove, L. E. et al. Brillouin neutron scattering of v-GeO₂. *Europhys. Lett.* **71**, 563 (2005).

29. Schirmacher, W., Diezemann, G. & Ganter, C. Harmonic vibrational excitations in disordered solids and the 'boson peak'. *Phys. Rev. Lett.* **81**, 136–139 (1998).

30. Chumakov, A. I. et al. Equivalence of the boson peak in glasses to the transverse acoustic van hove singularity in crystals. *Phys. Rev. Lett.* **106**, 225501 (2011).

31. Niss, K. et al. Influence of Pressure on the Boson Peak: Stronger than Elastic Medium Transformation. *Phys. Rev. Lett.* **99**, 055502 (2007).

32. Baron, A. Q. R. et al. An X-ray scattering beamline for studying dynamics. *J. Phys. Chem. Solids* **61**, 461–465 (2000).

33. Krisch, M. Status of phonon studies at high pressure by inelastic x-ray scattering. *J. Raman Spectrosc.* **34**, 628–632 (2003).

34. Mei, Q. et al. High-pressure x-ray diffraction measurements on vitreous GeO₂ under hydrostatic conditions. *Phys. Rev. B* **81**, 174113 (2010).

35. Mao, H.-K., Xu, J. & Bell, P. M. Calibration of the ruby pressure gauge to 800 kbar under quasi-hydrostatic conditions. *J. Geophys. Res.* **91** 4673–4676 (1986).

36. Gonze, X. et al. ABINIT: First-principles approach to material and nanosystem properties. *Comput. Phys. Commun.* **180**, 2582 (2009).

37. Gonze, X. & Lee, C. Dynamical matrices, Born effective charges, dielectric permittivity tensors, and interatomic force constants from density-functional perturbation theory. *Phys. Rev. B* **55**, 10355 (1997).

38. Madon, M., Gillet, Ph., Julien, Ch. & Price, G. D. A vibrational study of phase transitions among the GeO₂ polymorphs. *Phys. Chem. Miner.* **18**, 7–18 (1991).

39. Sharma, S. K. Applications of advanced Raman spectroscopic techniques in Earth Sciences, in Raman Spectroscopy: sixty years on vibrational spectra and structure ed. by Bist, H. D., Durig, J. R. & Sullivan, J. F. Elsevier Science, New York, vol. 17B, pp. 513–568 (1989).

40. Hermet, P., Lignie, A., Fraysse, G., Armand, P. & Papet, Ph. Thermodynamic properties of the quartz-type and rutile-type GeO₂ from first-principles calculations. *Phys. Chem. Chem. Phys.* **15**, 15943–15948 (2013).
## Acknowledgements

This work was supported by U. S. Department of Energy, Office of Science, Office of Basic Energy Sciences, under Contract No. DE-SC0012704. The authors acknowledge the technical and scientific support from both Spring-8 and ESRF staff. A.C., C.N.K., D.B. and Y.Q.C. acknowledge NSLS-II project for funding the travels for the measurements. D.A. acknowledge financial contribution from the French National Research Agency (ANR) through Grant 2010-JCJC-604-1. Calculations were performed at the National Energy Research Scientific Computing Center. Measurements at the beamline at BL35XU of SPring-8 were carried out using beamtime granted to the 2012A1122 research proposal.

## Author Contributions

A.C., D.B., D.A., F.B., Y.L. and Y.Q.C. designed the research. A.B., A.C., A.Q.B., C.N.K., D.A., F.B., H.-K.M., R.V., S.W., S.T. and Y.Q.C. prepared the high pressure cell and performed IXS measurements. Y.L. performed the numerical simulations. F.B. carried out the data analysis and A.C. wrote the manuscript. All authors discussed the results and commented on the manuscript.

## Additional Information

**Competing financial interests:** The authors declare no competing financial interests.

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

<figure><img src="image_8.png" /></figure>

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/