High-pressure densification of amorphous silica

Structure, Properties, and Applications of Silica Nanoparticles: Recent Theoretical Modeling Advances, Challenges, and Future Directions

Silica nanoparticles (SNPs), one of the most widely researched materials in modern science, are now commonly exploited in surface coatings, biomedicine, catalysis, and engineering of novel self-assembling materials. Theoretical approaches are invaluable to enhancing fundamental understanding of SNP properties and behavior. Tremendous research attention is dedicated to modeling silica structure, the silica-water interface, and functionalization of silica surfaces for tailored applications. In this review, the range of theoretical methodologies are discussed that have been employed to model bare silica and functionalized silica. The evolution of silica modeling approaches is detailed, including classical, quantum mechanical, and hybrid methods and highlight in particular the last decade of theoretical simulation advances. It is started with discussing investigations of bare silica systems, focusing on the fundamental interactions at the silica-water interface, following with a comprehensively review of the modeling studies that examine the interaction of silica with functional ligands, peptides, ions, surfactants, polymers, and carbonaceous species. The review is concluded with the perspective on existing challenges in the field and promising future directions that will further enhance the utility and importance of the theoretical approaches in guiding the rational design of SNPs for applications in engineering and biomedicine.

Thermobaric history as a tool to govern properties of glasses: case of dipropylene glycol

The elastic properties of high- and low-pressure glasses of dipropylene glycol were determined for the first time under conditions of isothermal compression up to 1 GPa at 77 K and isobaric heating of 77-300 K at 0.05 GPa and 1 GPa. A strong dependence of the elastic properties of glasses on their thermobaric history has been revealed: glasses obtained at high pressure have not only higher densities (3.9%), but also noticeably higher elastic moduli. This effect is especially pronounced in the shear modulus: high-pressure glass has a 30% higher shear modulus than low-pressure glass. The behavior of elastic moduli during the glass-to-liquid transition also depends on the thermobaric history. Glass produced at low pressure but heated at high pressure has anomalous temperature dependences of the elastic moduli. Heating dipropylene glycol glasses at different pressures allowed us to refine the Tg(P) dependence.

Atomistic insights into the structure and elasticity of densified 45S5 bioactive glasses

Glasses have applications in regenerative medicine due to their bioactivity, enabling interactions with hard and soft tissues. Soda-lime phosphosilicate glasses, such as 45S5, represent a model system of bioactive glasses. Regardless of their importance as bioactive materials, the relationship between the structure, density, and cooling process has not been studied in detail. This hinders further development of glasses as biomaterials. We used molecular dynamics simulations to study the elastic and structural properties of densified 45S5 bioactive glass and liquids over a wide range of densities. We performed a systematic analysis of the glass structure to density relationship to correlate the change in the properties with the structural change to enhance the mechanical properties of bioactive glasses while preserving their bioactive nature. The results show that the glass structure tends to be repolymerized, as indicated by increased network connectivity and a tetrahedral to octahedral polyhedral transition. We were able to tailor the elastic properties while keeping the bioactivity of the glass. The results presented here will provide some guidance to develop bioactive glasses with enhanced mechanical properties.

Negative density-dependence of the structural relaxation time of liquid silica: insights from a comparative molecular dynamics study

In many tetrahedral network-forming liquids, structural relaxation is anomalously accelerated by compression over relatively low pressure ranges. Here, for silica, we study this problem through comparative molecular dynamics simulations using two different models. Under compression, the network structures are compacted by slight tuning of the intertetrahedral bond angles while nearly preserving the unit tetrahedral structure. The consequent structural changes are remarkable for length scales larger than the nearest neighbor ion-pair distances. Accompanying with such structural changes, the interactions of the nearest Si-O pairs remain almost unchanged, whereas those of other ion pairs are, on average, strengthened by the degree of compression. In particular, the enhancement of the net Si-O interactions at the next nearest neighbor distance, which assist an ion in escaping from the potential well, reduces the activation energy, leading to a significant acceleration of structural relaxation. The results of our comparative molecular dynamics simulations are compatible with the scenario proposed by Angell, and further indicate that the structural relaxation dynamics cannot be uniquely determined by the configurations but strongly depends on the details of the coupling between the structure and the interaction.

In situ high pressure neutron diffraction and Raman spectroscopy of 20BaO-80TeO2 glass

The short-range structure of 20BaO-80TeO2 glass was studied in situ by high pressure neutron diffraction and high pressure Raman spectroscopy. Neutron diffraction measurements were performed at the PEARL instrument of the ISIS spallation neutron source up to a maximum pressure of 9.0 ± 0.5 GPa. The diffraction data was analysed via reverse Monte Carlo simulations and the changes in the glass short-range structural properties, Ba-O, Te-O and O-O bond lengths and speciation were studied as a function of pressure. Te-O co-ordination increases from 3.51 ± 0.05 to 3.73 ± 0.05, Ba-O coordination from 6.24 ± 0.19 to 6.99 ± 0.34 and O-O coordination from 6.00 ± 0.05 to 6.69 ± 0.06 with an increase in pressure from ambient to 9.0 GPa. In situ high pressure Raman studies found that the ratio of intensities of the two bands at 668 cm-1 and 724 cm-1 increases from 0.99 to 1.18 on applying pressure up to 19.28 ± 0.01 GPa, and that these changes are due to the conversion of TeO3 into TeO4 structural units in the tellurite network. It is found that pressure causes densification of the tellurite network by the enhancement of co-ordination of cations, and an increase in distribution of Te-O and Ba-O bond lengths. The original glass structure is restored upon the release of pressure.

Degree of Permanent Densification in Oxide Glasses upon Extreme Compression up to 24 GPa at Room Temperature

During the decompression of plastically deformed glasses at room temperature, some aspects of irreversible densification may be preserved. This densification has been primarily attributed to topological changes in glass networks. The changes in short-range structures like cation coordination numbers are often assumed to be relaxed upon decompression. Here the NMR results for aluminosilicate glass upon permanent densification up to 24 GPa reveal noticeable changes in the Al coordination number under pressure conditions as low as ∼6 GPa. A drastic increase in the highly coordinated Al fraction is evident over only a relatively narrow pressure range of up to ∼12 GPa, above which the coordination change becomes negligible up to 24 GPa. In contrast, Si coordination environments do not change, highlighting preferential coordination transformation during deformation. The observed trend in the coordination environment shows a remarkable similarity to the pressure-induced changes in the residual glass density, yielding a predictive relationship between the irreversible densification and the detailed structures under extreme compression. The results open a way to access the nature of plastic deformation in complex glasses at room temperature.

Oxygen Quadclusters in SiO_{2} Glass above Megabar Pressures up to 160 GPa Revealed by X-Ray Raman Scattering

As oxygen may occupy a major volume of oxides, a densification of amorphous oxides under extreme compression is dominated by reorganization of oxygen during compression. X-ray Raman scattering (XRS) spectra for SiO_{2} glass up to 1.6 Mbar reveal the evolution of heavily contracted oxygen environments characterized by a decrease in average O-O distance and the potential emergence of quadruply coordinated oxygen (oxygen quadcluster). Our results also reveal that the edge energies at the centers of gravity of the XRS features increase linearly with bulk density, yielding the first predictive relationship between the density and partial density of state of oxides above megabar pressures. The extreme densification paths with densified oxygen in amorphous oxides shed light upon the possible existence of stable melts in the planetary interiors.

Amorphous boron oxide at megabar pressures via inelastic X-ray scattering

Structural transition in amorphous oxides, including glasses, under extreme compression above megabar pressures (>1 million atmospheric pressure, 100 GPa) results in unique densification paths that differ from those in crystals. Experimentally verifying the atomistic origins of such densifications beyond 100 GPa remains unknown. Progress in inelastic X-ray scattering (IXS) provided insights into the pressure-induced bonding changes in oxide glasses; however, IXS has a signal intensity several orders of magnitude smaller than that of elastic X-rays, posing challenges for probing glass structures above 100 GPa near the Earth's core-mantle boundary. Here, we report megabar IXS spectra for prototypical B₂O₃ glasses at high pressure up to ∼120 GPa, where it is found that only four-coordinated boron (^[4]B) is prevalent. The reduction in the ^[4]B-O length up to 120 GPa is minor, indicating the extended stability of sp³-bonded ^[4]B. In contrast, a substantial decrease in the average O-O distance upon compression is revealed, suggesting that the densification in B₂O₃ glasses is primarily due to O-O distance reduction without the formation of ^[5]B. Together with earlier results with other archetypal oxide glasses, such as SiO₂ and GeO₂, the current results confirm that the transition pressure of the formation of highly coordinated framework cations systematically increases with the decreasing atomic radius of the cations. These observations highlight a new opportunity to study the structure of oxide glass above megabar pressures, yielding the atomistic origins of densification in melts at the Earth's core-mantle boundary.

Nanosecond homogeneous nucleation and crystal growth in shock-compressed SiO2

Understanding the kinetics of shock-compressed SiO2 is of great importance for mitigating optical damage for high-intensity lasers and for understanding meteoroid impacts. Experimental work has placed some thermodynamic bounds on the formation of high-pressure phases of this material, but the formation kinetics and underlying microscopic mechanisms are yet to be elucidated. Here, by employing multiscale molecular dynamics studies of shock-compressed fused silica and quartz, we find that silica transforms into a poor glass former that subsequently exhibits ultrafast crystallization within a few nanoseconds. We also find that, as a result of the formation of such an intermediate disordered phase, the transition between silica polymorphs obeys a homogeneous reconstructive nucleation and grain growth model. Moreover, we construct a quantitative model of nucleation and grain growth, and compare its predictions with stishovite grain sizes observed in laser-induced damage and meteoroid impact events.

Networks under pressure: the development of in situ high-pressure neutron diffraction for glassy and liquid materials

The pressure-driven collapse in the structure of network-forming materials will be considered in the gigapascal (GPa) regime, where the development of in situ high-pressure neutron diffraction has enabled this technique to obtain new structural information. The improvements to the neutron diffraction methodology are discussed, and the complementary nature of the results is illustrated by considering the pressure-driven structural transformations for several key network-forming materials that have also been investigated by using other experimental techniques such as x-ray diffraction, inelastic x-ray scattering, x-ray absorption spectroscopy and Raman spectroscopy. A starting point is provided by the pressure-driven network collapse of the prototypical network-forming oxide glasses B2O3, SiO2 and GeO2. Here, the combined results help to show that the coordination number of network-forming structural motifs in a wide range of glassy and liquid oxide materials can be rationalised in terms of the oxygen-packing fraction over an extensive pressure and temperature range. The pressure-driven network collapse of the prototypical chalcogenide glass GeSe2 is also considered where, as for the case of glassy GeO2, site-specific structural information is now available from the method of in situ high-pressure neutron diffraction with isotope substitution. The application of in situ high-pressure neutron diffraction to other structurally disordered network-forming materials is also summarised. In all of this work a key theme concerns the rich diversity in the mechanisms of network collapse, which drive the changes in physico-chemical properties of these materials. A more complete picture of the mechanisms is provided by molecular dynamics simulations using theoretical schemes that give a good account of the experimental results.

High-pressure transformation of SiO₂ glass from a tetrahedral to an octahedral network: a joint approach using neutron diffraction and molecular dynamics

A combination of in situ high-pressure neutron diffraction at pressures up to 17.5(5) GPa and molecular dynamics simulations employing a many-body interatomic potential model is used to investigate the structure of cold-compressed silica glass. The simulations give a good account of the neutron diffraction results and of existing x-ray diffraction results at pressures up to ~60  GPa. On the basis of the molecular dynamics results, an atomistic model for densification is proposed in which rings are "zipped" by a pairing of five- and/or sixfold coordinated Si sites. The model gives an accurate description for the dependence of the mean primitive ring size ⟨n⟩ on the mean Si-O coordination number, thereby linking a parameter that is sensitive to ordering on multiple length scales to a readily measurable parameter that describes the local coordination environment.

Progressive transformations of silica glass upon densification

The elastic and plastic behaviors of silica glasses densified at various maximum pressure reached (12 GPa, 15 GPa, 19 GPa, and 22 GPa), were analyzed using in situ Raman and Brillouin spectroscopies. The elastic anomaly was observed to progressively vanish up to a maximum pressure reached of 12 GPa, beyond which it is completely suppressed. Above the elastic anomaly the mechanical behavior of silica glass, as derived from Brillouin measurements, is interpreted in terms of pressure induced transformation of low density amorphous silica into high density amorphous silica.

Effect of helium on structure and compression behavior of SiO2 glass

The behavior of volatiles is crucial for understanding the evolution of the Earth's interior, hydrosphere, and atmosphere. Noble gases as neutral species can serve as probes and be used for examining gas solubility in silicate melts and structural responses to any gas inclusion. Here, we report experimental results that reveal a strong effect of helium on the intermediate range structural order of SiO(2) glass and an unusually rigid behavior of the glass. The structure factor data show that the first sharp diffraction peak position of SiO(2) glass in helium medium remains essentially the same under pressures up to 18.6 GPa, suggesting that helium may have entered in the voids in SiO(2) glass under pressure. The dissolved helium makes the SiO(2) glass much less compressible at high pressures. GeO(2) glass and SiO(2) glass with H(2) as pressure medium do not display this effect. These observations suggest that the effect of helium on the structure and compression of SiO(2) glass is unique.

In situ high pressure and high temperature Raman studies of (1-x)SiO(2)xGeO(2) glasses

The structure of glasses in the binary system SiO(2)-GeO(2) has been studied by Raman spectroscopy. Our results are consistent with mixing of SiO(2) and GeO(2) tetrahedra. The changes induced by temperature and by pressure on the structure are monitored by in situ measurements on the same mixed glasses. Anomalous temperature dependences are observed not only for SiO(2) glass and GeO(2) glass but also for mixed glasses. Particular attention is focused on the pressure densification mechanism in mixed glasses. Via the pressure dependence of the width of the main Raman band, we show that the compression mechanism in mixed glasses is intermediate between that of the end members.

Sixfold-coordinated amorphous polymorph of SiO2 under high pressure

We have developed synchrotron x-ray absorption and diffraction techniques for measuring the density and structure of noncrystalline materials at high pressures and have applied them to studying the behavior of SiO2 glass. The density, coordination number, and Si-O bond length at a pressure of 50 GPa were measured to be 4.63 g/cm;{3}, 6.3, and 1.71 A, respectively. Based on the density data measured in this study and the sound velocity data available in the literature, the bulk modulus at 50 GPa was estimated to be 390 GPa, which is consistent with the pressure dependence of the density in the vicinity of 50 GPa. These results, together with the knowledge from our exploratory study, suggest that SiO2 glass behaves as a single amorphous polymorph having a sixfold-coordinated structure at pressures above 40-45 GPa up to at least 100 GPa.

Poisson's ratio and the densification of glass under high pressure

Because of a relatively low atomic packing density, (Cg) glasses experience significant densification under high hydrostatic pressure. Poisson's ratio (nu) is correlated to Cg and typically varies from 0.15 for glasses with low Cg such as amorphous silica to 0.38 for close-packed atomic networks such as in bulk metallic glasses. Pressure experiments were conducted up to 25 GPa at 293 K on silica, soda-lime-silica, chalcogenide, and bulk metallic glasses. We show from these high-pressure data that there is a direct correlation between nu and the maximum post-decompression density change.

Comparison of model potentials for molecular-dynamics simulations of silica

Structural, thermomechanical, and dynamic properties of pure silica SiO2 are calculated with three different model potentials, namely, the potential suggested by van Beest, Kramer, and van Santen (BKS) [Phys. Rev. Lett. 64, 1955 (1990)], the fluctuating-charge potential with a Morse stretch term for the short-range interactions proposed by Demiralp, Cagin, and Goddard (DCG)[Phys. Rev. Lett. 82, 1708 (1999)], and a polarizable force field proposed by Tangney and Scandolo (TS) [J. Chem. Phys. 117, 8898 (2002)]. The DCG potential had to be modified due to flaws in the original treatment. While BKS reproduces many thermomechanical properties of different polymorphs rather accurately, it also shows qualitatively wrong trends concerning the phononic density of states, an absence of the experimentally observed anomaly in the c/a ratio at the quartz alpha-beta transition, pathological instabilities in the beta-cristobalite phase, and a vastly overestimated transition pressure for the stishovite I --> II transition. These shortcomings are only partially remedied by the modified DCG potential but greatly improved by the TS potential. DCG and TS both reproduce a pressure-induced transition from alpha-quartz to quartz II, predicted theoretically based on the BKS potential.

Densification of silica glass at ambient pressure

We show that densification of silica glass at ambient pressure as observed in irradiation experiments can be attributed to defect generation and subsequent structure relaxation. In our molecular dynamics simulations, defects are created by randomly removing atoms, by displacing atoms from their nominal positions in an otherwise intact glass, and by assigning certain atom excess kinetic energy (simulated ion implantation). The former forms vacancies; displacing atoms and ion implantation produce both vacancies and "interstitials." Appreciable densification is induced by these defects after equilibration of the defective glasses. The structural and vibrational properties of the densified glasses are characterized, displaying resembling features regardless of the means of densification. These results indicate that relaxation of high free-energy defects into metastable amorphous structures enriched in atomic coordination serves as a common mechanism for densification of silica glass at ambient pressure.

Molecular dynamics study of CaSiO(3)-MgSiO(3) glasses under high pressure

The pressure-induced structural evolutions of CaSiO(3)-MgSiO(3) glasses have been examined by means of molecular dynamics simulation. Our calculations revealed that Si coordination remained unchanged up to 15 GPa, while modifier cations caused significant changes in the short-range order structure. In the present study, we conclude that the main compression mechanisms for CaSiO(3)-MgSiO(3) glasses are: (1) the Si-O-Si angle reduction, (2) the coordination increase of Ca and Mg cations, and (3) the compaction in the medium-range scale. Furthermore, small changes in the Q(n) distribution suggest pressure-induced disproportionation reactions. Similar pressure responses between CaSiO(3)-MgSiO(3) glasses may imply that the structural changes of SiO(4) framework units are more significant than those of interstitial cations, Ca and Mg.

A QUICKSTEP-based quantum mechanics/molecular mechanics approach for silica

Quantum mechanics/molecular mechanics (QM/MM) approaches are currently used to describe several properties of silica-based systems, which are local in nature and require a quantum description of only a small number of atoms around the site of interest, e.g., local chemical reactivity or spectroscopic properties of point defects. We present a QM/MM scheme for silica suitable to be implemented in the general QM/MM framework recently developed for large scale molecular dynamics simulations, within the QUICKSTEP approach to the description of the quantum region. Our scheme has been validated by computing the structural and dynamical properties of an oxygen vacancy in alpha-quartz, a prototypical defect in silica. We have found that good convergence in the Si-Si bond length and formation energy is achieved by using a quantum cluster of only eight atoms in size. We check the suitability of the method for molecular dynamics and evaluate the Si-Si bond frequency from the velocity-velocity correlation function.

Ab initio simulations of photoinduced interconversions of oxygen deficient centers in amorphous silica

We have studied by ab initio molecular dynamics the interconversion between oxygen deficient centers (Si-Si bond, dicoordinated silicon =Si:, and E' centers) induced by UV irradiation in a-SiO2. By dynamical simulations in the excited state of a periodic model of a-SiO2 we have identified the reaction path and activation barrier for the Si-Si-->=Si: interconversion. A new competitive transformation of the excited, neutral Si-Si bond into two E' centers has been identified. Our results provide strong theoretical support to the viability of these processes, proposed experimentally on the basis of optical data only.

Relaxation dynamics of a viscous silica melt: the intermediate scattering functions

We use molecular dynamics computer simulations to study the relaxation dynamics of a viscous melt of silica. The coherent and incoherent intermediate scattering functions, F(q,t) and F(s)(q,t), show a crossover from a nearly exponential decay at high temperatures to a two-step relaxation at low temperatures. Close to the critical temperature of mode-coupling theory (MCT) the correlators obey in the alpha regime the time temperature superposition principle (TTSP) and show a weak stretching. We determine the wave-vector dependence of the stretching parameter and find that for F(q,t) it shows oscillations that are in phase with the static structure factor. The temperature dependence of the alpha-relaxation times tau shows a crossover from an Arrhenius law at low temperatures to a weaker T dependence at intermediate and high temperatures. At the latter temperatures the T dependence is described well by the power law proposed by MCT with the same critical temperature that has previously been found for the diffusion constant D and the viscosity. We find that the exponent gamma of the power law for tau are significantly larger than the one for D. The wave-vector dependence of the alpha-relaxation times for F(q,t) oscillates around tau(q) for F(s)(q,t) and is in phase with the structure factor. Due to the strong vibrational component of the dynamics at short times the TTSP is not valid in the beta-relaxation regime. We show, however, that in this time window the shape of the curves is independent of the correlator and is given by a functional form proposed by MCT. We find that the value of the von Schweidler exponent and the value of gamma for finite q are compatible with the expression proposed by MCT. Finally we discuss the q dependence of the critical amplitude and the correction term and find that they are qualitatively similar to the ones for simple liquids and the prediction of MCT. We conclude that, in the temperature regime where the relaxation times are mesoscopic, many aspects of the dynamics of this strong glass former can be rationalized very well by MCT.

Low-temperature behavior of core-softened models: water and silica behavior

A core-softened model of a glass forming fluid is numerically studied in the limit of very low temperatures. The model shows two qualitatively different behaviors depending on the strength of the attraction between particles. For no or low attraction, the changes of density as a function of pressure are smooth, although hysteretic due to mechanical metastabilities. For larger attraction, sudden changes of density upon compressing and decompressing occur. This global mechanical instability is correlated to the existence of a thermodynamic first-order amorphous-amorphous transition. The two different behaviors obtained correspond qualitatively to the different phenomenology observed in silica and water.

Computer simulations of liquid silica: equation of state and liquid-liquid phase transition

We conduct extensive molecular dynamics computer simulations of two models for liquid silica [the model of Woodcock, Angell and Cheeseman, J. Phys. Chem. 65, 1565 (1976); and that of van Beest, Kramer, and van Santen, Phys. Rev. Lett. 64, 1955 (1990)] to determine their thermodynamic properties at low temperature T across a wide density range. We find for both models a wide range of states in which isochores of the potential energy U are a linear function of T(3/5), as recently proposed for simple liquids [Rosenfeld and P. Tarazona, Mol. Phys. 95, 141 (1998)]. We exploit this behavior to fit an accurate equation of state to our thermodynamic data. Extrapolation of this equation of state to low T predicts the occurrence of a liquid-liquid phase transition for both models. We conduct simulations in the region of the predicted phase transition, and confirm its existence by direct observation of phase separating droplets of atoms with distinct local density and coordination environments.

First-order amorphous-amorphous transformation in silica

Molecular simulations predict that a first-order amorphous-amorphous transformation occurs in SiO2 under pressure, analogous to the first-order amorphous-amorphous transformation known to occur in H2O. At low temperatures the first-order transformation is kinetically hindered, and an amorphous-amorphous transformation occurs instead by gradual spinodal decomposition at higher pressures. We suggest that previous experiments have observed the spinodal decomposition pathway in SiO2 and that the predicted first-order transformation will be observed in experiments carried out at higher temperatures.

Two-membered silicon rings on the dehydroxylated surface of silica

We present extensive modeling of the amorphous silica surface, aimed at connecting its structural and chemical features. beta-cristobalite surfaces are initially studied to model the hydroxylated surfaces. A model reconstruction of the (111) surface is used to define a path leading to the formation of two-membered silicon rings upon dehydroxylation. Subsequently, a realistic model of the amorphous dehydroxylated (dry) surface is produced, by full ab initio annealing of an initial model generated by classical simulation. The presence of surface two-membered silicon rings emerges naturally. A calculation of IR activity yields an associated peak doublet in agreement with experimental data.

Molecular dynamics simulation of structural transformation in silicon carbide under pressure

Pressure-induced structural transformation in cubic silicon carbide is studied with the isothermal-isobaric molecular-dynamics method using a new interatomic potential scheme. The reversible transformation between the fourfold coordinated zinc-blende structure and the sixfold coordinated rocksalt structure is successfully reproduced by the interatomic potentials. The calculated volume change at the transition and hysteresis are in good agreement with experimental data. The atomistic mechanisms of the structural transformation involve a cubic-to-monoclinic unit-cell transformation and a relative shift of Si and C sublattices in the 100 direction.

High-pressure infrared and FT-Raman investigation of a dental composite

Composite resins are often used as filling materials on load-bearing surfaces of teeth. As masticatory stresses can be high, here, we study the effect of pressure on the behaviour of a dental composite. Using a polymerized wafer, the IR and FT-Raman spectra of a zirconia-containing proprietary composite (Z100, 3M, Minneapolis, MN, USA) were recorded. The high-pressure IR spectra were also recorded. Band assignments were made for the main peaks of both organic and inorganic components. Breaks in the pressure dependences (dv/dP) of the organic components were found at 22 kbar. Different pressure dependences for different vibrational modes of inorganic components were also observed. These data suggest that the network structure of the composite is compacted under high pressure and that both the atomic distance and bonding angles in the network are altered.