# Abstract

Elastic geobarometry for host-inclusion systems can provide new constraints to assess the pressure and temperature conditions attained during metamorphism. Current experimental approaches and theory are developed only for crystals immersed in a hydrostatic stress field, whereas inclusions experience deviatoric stress. We have developed a method to determine the strains in quartz inclusions from Raman spectroscopy using the concept of the phonon-mode Grüneisen tensor. We used ab initio Hartree-Fock/Density Functional Theory to calculate the wavenumbers of the Raman-active modes as a function of different strain conditions. Least-squares fits of the phonon-wavenumber shifts against strains have been used to obtain the components of the mode Grüneisen tensor of quartz (

## Introduction

Mineral inclusions entrapped in ultrahigh-pressure metamorphic rocks can provide fundamental information about geological processes such as subduction and continental collision. When a host-inclusion pair is exhumed from depth to the Earth’s surface non-lithostatic stresses are developed in the inclusion because of the contrast in their elastic properties (Angel et al. 2015 and references therein). The inclusion is under compressive stress if it is more compressible than the host. If correctly interpreted, these stresses on the inclusion allow the stress conditions at entrapment to be estimated. However, the theory for elastic geobarometry for host-inclusion pairs is only developed for crystals immersed in a hydrostatic stress field. This is valid for an isotropic and spherical inclusion entrapped in an isotropic host as it is subject to isotropic strains imposed by the host and will, therefore, exhibit isotropic stresses; that means the inclusion will be under hydrostatic pressure. However, if inclusions of elastically anisotropic minerals such as quartz, coesite, or zircon are entrapped in a cubic host such as garnet the same isotropic strain imposed by the garnet will result in anisotropic (i.e., deviatoric) stresses in the inclusion (see Fig. 1). Thus, the stress state of the inclusion will be different from the hydrostatic case, and it cannot be characterized by a single “pressure” value (Anzolini et al. 2018). Furthermore, the shifts of Raman mode frequencies under deviatoric stress are unknown in general. The limited experimental evidence (e.g., Briggs and Ramdas 1977) is that the Raman modes change differently from those measured under hydrostatic pressure. Therefore, the key question remains, posed by Korsakov et al. (2010), can shifts of Raman modes measured in hydrostatic experiments be used to interpret the Raman shifts from an inclusion under non-hydrostatic stress? And what errors does this introduce into estimates of inclusion stress and entrapment conditions?

### Figure 1

Since it is challenging to perform experiments under controlled deviatoric stress conditions, we have performed ab initio Hartree-Fock/Density Functional Theory (HF/DFT) calculations to determine how the Raman modes of quartz change under both hydrostatic pressure and deviatoric stress conditions. We show how the HF/DFT simulations can be used to determine the strains (and by inference the stresses) within a crystal inclusion via mode Grüneisen tensors and measurement of its Raman shifts.

## Grüneisen tensor and strains

It is well known that when stresses are applied to a crystal, for example by changing the pressure, the phonon wavenumbers are shifted. This is most easily seen in the change in Raman peak positions measured under hydrostatic pressure. The shifts of certain modes of materials (e.g., the 464 cm^{–1} mode of quartz) have frequently been used as secondary pressure standards. However, there is a common misconception that the observed wavenumber shifts of the Raman-active modes are directly related to the magnitude of the applied pressure or stress. If that were true, then Raman modes would not show a change in wavenumber when the temperature of a free crystal is changed at ambient pressure. Instead, Raman modes generally exhibit a decrease in wavenumber as the temperature is increased. As Figure 2 shows, the shift of the 464 cm^{–1} mode of quartz with both pressure and temperature shows the same dependence upon volume. These observations point to the correct interpretation that the shift Δω of the wavenumber ω of a vibrational mode is primarily due to the strains of the crystal induced by the applied temperature or pressure. Thus, the Raman peak position of a crystal under a strain ɛ (i.e., the full strain tensor) is determined by the second-rank symmetric tensor: the mode Grüneisen tensor γ^{m} (Ziman 1960; Key 1967; Cantrell 1980) characteristic of each phonon mode *m*, which can be written in Voigt (1910) notation as:

### Figure 2

This equation means that the changes in the Raman peak positions, in general, depend on all the components of the strain tensor (see Supplemental material^{[1]}), not just on the relative change in the volume (i.e., ɛ_{1} + ɛ_{2} + ɛ_{3}). Because the mode Grüneisen tensor is a symmetric second-rank property tensor, it is subject to the same symmetry constraints on its component values as other second-rank property tensors, such as the thermal expansion and compressibility tensors. For the trigonal symmetry of quartz _{1} = ɛ_{2} for which the shift in the phonon wavenumbers should be given by:

To determine the values of ^{13}it is not sufficient to measure the Raman shifts under hydrostatic pressure, because only a single series of values of ɛ_{1} and ɛ_{3} are measured. We, therefore, use HF/DFT simulations to calculate the Raman modes at different imposed strains on the crystal to determine its structure and properties at those strains.

## Raman shifts and strains from ab initio calculations

Ab initio HF/DFT simulations have been performed by means of CRYSTAL 14 code (Dovesi et al. 2014) employing the hybrid Hamiltonian WC1LYP (Wu and Cohen 2006) that has been demonstrated to be particularly suitable for the accurate reproduction of the elastic and vibrational properties (Zicovich-Wilson et al. 2004; Prencipe et al. 2014). Further computational details are reported in the Supplemental material^{[1]}.

We will only discuss in detail the Raman-active modes near 464 cm^{–1} (non-degenerate A_{1} mode) and 696 cm^{–1} (doubly degenerate E mode) because they give peaks in the Raman spectra that are easily resolved from the peaks of the host garnets. Note that all E modes in quartz are polar, and therefore have longitudinal optical (LO) and transverse optical (TO) components that generate two separate Raman peaks, whose intensity ratio varies depending on the scattering geometry. Because of the polarization mixing (Loudon 1964), the wavenumber of the LO component depends on the angle between the triad axis of quartz and the phonon-propagation direction that in the case of backscattering geometry coincides with the direction of the laser beam. For the E mode of interest in this study, the LO-TO splitting is small. At ambient conditions w_{ETO} = 696 cm^{–1}, while the maximum w_{ELO} ~ 697.5 cm^{–1} and it corresponds to the case when the **c** axis is perpendicular to the laser beam. The two components are thus very close to each other and in specific experimental geometries may generate peaks of similar intensity. To avoid possible contribution from the LO component that may lead to a subtle artificial shift of the corresponding Raman peak toward higher wavenumbers, one should either rotate the sample about the direction of the laser beam to verify that the peak position does not change, or find the orientation at which the wavenumber of the Raman peak is lowest (which should be when the **c** axis is perpendicular to the polarization of the incident light).

Figures 3a and 3b are contour maps that display the HF/DFT wavenumber shifts of the A1 mode near 464 cm^{–1} and ETO mode near 696 cm^{–1} as a function of the two independent strain components. For small strains (close to the origin) the iso-shift lines are parallel to one another and equally spaced, confirming that the values of γ_{1}^{m} and γ_{3}^{m} are constants over these strain ranges. The contour lines for the two modes have different slopes (e.g., Figs. 3a, 3b, and 3c) indicating that the values of their Grüneisen components are different. Furthermore, Figure 3a shows that the iso-shift lines are not parallel to isochors, which are represented by lines of equal volume strain [ɛ*V* = 2ɛ_{1} + ɛ_{3}]. If the wavenumber shifts are plotted against the stress components σ_{1} and σ_{3} (Fig. S4, Supplemental^{[1]} material) then the iso-shift lines are not parallel to isobars [lines of constant pressure (2σ_{1} + σ_{3})/3 = –*P*]. Therefore, in general, Raman shifts do not measure either volume or pressures; for uniaxial crystals like quartz, they indicate the principal normal strain components ɛ_{1} and ɛ_{3}. The Grüneisen tensor components for all modes were then determined by fitting Equation 2 by least-squares to the wavenumber shifts at different strain states simulated by HF/DFT calculations. For the 464 mode, the maximum misfit was 1.76 cm^{–1} (see Supplemental material^{[1]} for residual plots and full details).

### Figure 3

HF/DFT calculations are performed under static stress at absolute-zero temperature. Inclusions are measured at room temperature. To use the mode Grüneisen tensors to determine strains in inclusions, one has to also demonstrate that their values are independent of *P* and *T*. This is equivalent to the solid behaving according to the quasiharmonic approximation (QHA). We achieve this by calculating the strains of a free quartz crystal relative to room conditions from the known unit-cell parameter variation of quartz (Angel et al. 2017) with *P* and *T*. We then compare the wavenumber shifts from room-condition values using the Grüneisen components (^{–1} (i.e., 0.60 and 1.19; 0.50 and 0.36) determined by HF/DFT and Equation 2 with experimental data. The line in Figure 2 shows that for the mode near 464 cm^{–1} the experimentally measured wavenumber shifts under high pressure at room *T* and those up to ~400 °C at room *P* are reproduced by the Grüneisen components determined by HF/DFT. At higher temperatures, the predicted shifts differ from the experimental data because of pre-transition effects associated with the α-β quartz transition that cannot be accounted using QHA. Furthermore, the fundamental soft mode near 206 cm^{–1}, which is heavily involved in the temperature α- to β-quartz phase transition (Scott 1968), clearly violates the QHA, and HF/DFT simulations cannot be used to determine the mode Grüneisen components. To use this band to determine strains in inclusions, the variation of its position with strain must be determined experimentally.

## Validation from Raman scattering and X‑ray diffraction

To validate this approach, we performed micro-Raman spectroscopy and X‑ray diffraction measurements on quartz inclusion in garnet from a diamond-grade eclogite xenolith (TM90-1) from the Mir kimberlite pipe (Yakutiya). We selected this example as it has the highest Raman shifts reported for quartz inclusions in garnet (Korsakov et al. 2009; Zhukov and Korsakov 2015). Parallel polarized Raman spectra were collected in backscattering geometry with a Horiba Jobin-Yvon T64000 triple-monochromator spectrometer (spectral resolution of ~2 cm^{–1}, instrumental accuracy in peak positions of ~0.35 cm^{–1} and 2 μm spot size) following the same protocol reported in Campomenosi et al. (in review). Figure 4 shows that the shift of the 464 cm^{–1} Raman line changes significantly across the crystal as a consequence of edge and corner effects (Campomenosi et al. in review; Mazzucchelli et al. 2018), and is lowest at the center of the inclusion. We determined the unit-cell parameters of the inclusion by single-crystal X‑ray diffraction measurements using the 8-position centering method (SINGLE, Angel and Finger 2011) using a newly developed Huber 4-circle Eulerian cradle diffractometer equipped with point detector and microfocus source (120 μm spot size). The difference between the unit-cell parameters of the inclusion *a* = 4.86669(44) Å, *c* = 5.35408(14) Å, and those of a free crystal measured on the same instrument [*a* = 4.91160(7) Å, *c* = 5.40325(9) Å] shows that the inclusion is under strains ɛ_{1} = –0.00914(9), ɛ_{3} = –0.00910(26), and ɛ_{V} = –0.02714(30). Relative to room conditions the quartz inclusion is currently under isotropic strain within the experimental uncertainties. From these strains and the Grüneisen-tensor components for the mode near 464 cm^{–1} (i.e., ^{–1}. This is between the minimum and maximum shifts actually measured on the inclusion (Fig. 4) because the unit-cell parameters measured by X‑ray diffraction are an average over the entire volume of the inclusion. Conversely, if we take the shifts of the Raman modes near 464 and 696 cm^{–1} measured at the center of the inclusion and the calculated mode Grüneisen components we predict the following strains: ɛ_{1} = –0.0093(5) and ɛ_{3} = –0.0070(5). In this case, the total strain is lower than that measured by X‑ray diffraction, again because this method provides an average over the inclusion.

### Figure 4

## Implications

We have demonstrated by a combination of HF/DFT simulations and comparison to experimental Raman data that quartz has some Raman-active modes whose wavenumbers are only a function of strain (for small strains) and are not directly dependent on *P* and *T*. The HF/DFT simulations show also that in general the shift of a Raman line does not indicate either the volume change (or volume strain) or the mean normal stress but is a more complex function of the linear strains (Fig. 3). The linear strains of a quartz inclusion can be determined from the wavenumber shifts of at least two Raman peaks by using the mode Grüneisen components determined by HF/DFT. Conversely, a wavenumber shift for a specific Raman-active mode can be determined from the linear strains measured by X‑ray diffraction and by using the mode Grüneisen components determined by HF/DFT. We proved that strains and Raman shifts vary significantly across an inclusion (Fig. 4) as a result of the influence of shape combined with elastic anisotropy (Mazzucchelli et al. 2018; Campomenosi et al. in review), and should never be averaged. The Raman shifts measured in the center of inclusions are least affected by the presence of edges and should be the only ones used to infer, after shape corrections (Mazzucchelli et al. 2018), the entrapment conditions. On the other hand, X‑ray diffraction measurements give strain components that are averaged over the entire volume of the investigated crystal. For the example discussed, the interpretation of wavenumber shift of the mode near 464 cm^{–1} measured at the center of the quartz inclusion with the hydrostatic calculation (Schmidt and Ziemann 2000) leads to only a small error in the estimation of the mean normal stress (i.e., pressure) of ca. 0.03 GPa. This would lead to an error in the calculated entrapment pressure of 0.06 GPa. However, this is not always the case. When the deviatoric stress (Fig. 1) in the inclusion is higher, the error is larger. For example, a quartz inclusion in garnet reset at 1.5 GPa and 1100 °C during exhumation will have a wavenumber shift of the Raman mode near 464 cm^{–1} that would yield an entrapment pressure of 0.6 GPa too low when the hydrostatic calibration is used. In anisotropic hosts, the discrepancy from the hydrostatic calibration also depends on the relative orientation of the host and inclusion. For example, for quartz in zircon the error in the calculated entrapment pressure can vary from 0.2 to 0.6 GPa for the same entrapment condition (i.e., 0.3 GPa and 800 °C) depending on orientation.

# Acknowledgments

This project received funding from the European Research Council under the European Union’s Horizon 2020 research and innovation program grant agreements 714936 and from the MIUR-SIR (Ministry of Education, University and Research–Scientific Independence of Young Researchers, Italy) grant MILE DEEp (Mineral Inclusion Elasticity for a New Deep Subduction Geobarometer; RBSI140351) to Matteo Alvaro. Andrey V. Korsakov was supported by the Russian Federation state assignment project no. 0330-2016-0006.

### References cited

Angel, R.J., and Finger, L.W. (2011) SINGLE: a program to control single-crystal diffractometers. Journal of Applied Crystallography, 44(1), 247–251.10.1107/S0021889810042305Search in Google Scholar

Angel, R.J., Nimis, P., Mazzucchelli, M.L., Alvaro, M., and Nestola, F. (2015) How large are departures from lithostatic pressure? Constraints from host–inclusion elasticity. Journal of Metamorphic Geology, 33(8), 801–813.10.1111/jmg.12138Search in Google Scholar

Angel, R.J., Alvaro, M., Miletich, R., and Nestola, F. (2017) A simple and generalised P–T–V EoS for continuous phase transitions, implemented in EosFit and applied to quartz. Contributions to Mineralogy and Petrology, 172(5), 29.10.1007/s00410-017-1349-xSearch in Google Scholar

Anzolini, C., Prencipe, M., Alvaro, M., Romano, C., Vona, A., Lorenzon, S., Smith, E.M., Brenker, F.E., and Nestola, F. (2018) Depth of formation of super-deep diamonds: Raman barometry of CaSiO_{3}-walstromite inclusions. American Mineralogist, 103(1), 69–74.10.2138/am-2018-6184Search in Google Scholar

Briggs, R.J., and Ramdas, A.K. (1977) Piezospectroscopy of the Raman spectrum of α-quartz. Physical Review B, 16(8), 3815.10.1103/PhysRevB.16.3815Search in Google Scholar

Cantrell, J.H. Jr. (1980) Generalized Grüneisen tensor from solid nonlinearity parameters. Physical Review B, 21(10), 4191.10.1103/PhysRevB.21.4191Search in Google Scholar

Civalleri, B., D’Arco, P., Orlando, R., Saunders, V., and Dovesi, R. (2001) Hartree–Fock geometry optimisation of periodic systems with the CRYSTAL code. Chemical Physics Letters, 348, 131–138.10.1016/S0009-2614(01)01081-8Search in Google Scholar

Dovesi, R., Orlando, R., Erba, A., Zicovich-Wilson, C.M., Civalleri, B., Casassa, S., Maschio, L., Ferrabone, M., De La Pierre, M., D’Arco, P., and Noël, Y. (2014) CRYSTAL14: A program for the ab initio investigation of crystalline solids. International Journal of Quantum Chemistry, 114, 1287–1317.10.1002/qua.24658Search in Google Scholar

Key, S.W. (1967) Grüneisen tensor for anisotropic materials. Journal of Applied Physics, 38(7), 2923–2928.10.1063/1.1710025Search in Google Scholar

Korsakov, A.V., Perraki, M., Zhukov, V.P., De Gussem, K., Vandenabeele, P., and Tomilenko, A.A. (2009) Is quartz a potential indicator of ultrahigh-pressure metamorphism? Laser Raman spectroscopy of quartz inclusions in ultrahigh-pressure garnets. European Journal of Mineralogy, 21(6), 1313–1323.10.1127/0935-1221/2009/0021-2006Search in Google Scholar

Korsakov, A.V., Zhukov, V.P., and Vandenabeele, P. (2010) Raman-based geobarometry of ultrahigh-pressure metamorphic rocks: applications, problems, and perspectives. Analytical and Bioanalytical Chemistry, 397(7), 2739–2752.10.1007/s00216-010-3831-4Search in Google Scholar PubMed

Kuzmany, H. (2009) Solid-State Spectroscopy, An Introduction, second edition, cap-9, p. 185–200. Springer. doi: 10.1007/978-3-642-01479-6.10.1007/978-3-642-01479-6Search in Google Scholar

Lakshtanov, D.L., Sinogeikin, S.V., and Bass, J.D. (2007) High-temperature phase transitions and elasticity of silica polymorphs. Physics and Chemistry of Minerals, 34(1), 11–22.10.1007/s00269-006-0113-ySearch in Google Scholar

Loudon, R. (1964) The Raman effect in crystals. Advances in Physics, 13(52), 423–482.10.1080/00018736400101051Search in Google Scholar

Mazzucchelli, M.L., Burnley, P., Angel, R.J., Morganti, S., Domeneghetti, M.C., Nestola, F., and Alvaro, M. (2018) Elastic geothermobarometry: Corrections for the geometry of the host-inclusion system. Geology, 46(3), 231–234.10.1130/G39807.1Search in Google Scholar

Milani, S., Nestola, F., Alvaro, M., Pasqual, D., Mazzucchelli, M.L., Domeneghetti, M.C., and Geiger, C.A. (2015) Diamond–garnet geobarometry: The role of garnet compressibility and expansivity. Lithos, 227, 140–147.10.1016/j.lithos.2015.03.017Search in Google Scholar

Prencipe, M., Maschio, L., Kirtman, B., Salustro, S., Erba, A., and Dovesi, R. (2014) Raman spectrum of NaAlSi_{2}O_{6} jadeite. A quantum mechanical simulation. Journal of Raman Spectroscopy, 45(8), 703–709.10.1002/jrs.4519Search in Google Scholar

Schmidt, C., and Ziemann, M.A. (2000) In-situ Raman spectroscopy of quartz: A pressure sensor for hydrothermal diamond-anvil cell experiments at elevated temperatures. American Mineralogist, 85, 1725–1734.10.2138/am-2000-11-1216Search in Google Scholar

Scott, J.F. (1968) Evidence of coupling between one- and two-phonon excitations in quartz. Physical Review Letters, 21(13), 907–910.10.1103/PhysRevLett.21.907Search in Google Scholar

Sinogeikin, S.V., and Bass, J.D. (2002) Elasticity of majorite and a majorite-pyrope solid solution to high pressure: Implications for the transition zone. Geophysical Research Letters, 29(2).10.1029/2001GL013937Search in Google Scholar

Voigt, W. (1910). Lehrbuch der Kristallphysik, p.754. B.G. Teubner, Berlin.Search in Google Scholar

Wu, Z., and Cohen, R.E. (2006) More accurate generalized gradient approximation for solids. Physical Review B, 73, 235116.10.1103/PhysRevB.73.235116Search in Google Scholar

Zhukov, V.P., and Korsakov, A.V. (2015) Evolution of host-inclusion systems: a viscoelastic model. Journal of Metamorphic Geology, 33(8), 815–828.10.1111/jmg.12149Search in Google Scholar

Zicovich-Wilson, C.M., Pascale, F., Roetti, C., Saunders, V.R., Orlando, R., and Dovesi, R. (2004) Calculation of the vibration frequencies of α-quartz: The effect of Hamiltonian and basis set. Journal of Computational Chemistry, 25(15), 1873–1881.10.1002/jcc.20120Search in Google Scholar PubMed

Ziman, J.M. (1960) Electrons and Phonons: The theory of transport phenomena in solids. Oxford University Press, p. 469.Search in Google Scholar

**Received:**2018-05-10

**Accepted:**2018-05-31

**Published Online:**2018-10-30

**Published in Print:**2018-11-27

© 2018 Walter de Gruyter GmbH, Berlin/Boston

This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 3.0 License.