Fossilized high pressure from the Earth's deep interior: the coesite-in-diamond barometer

Natural diamond formation by self-redox of ferromagnesian carbonate

The presence of extra reducer was thought to be essential for producing natural diamonds from reduction of carbonates. The present study of the Xiuyan meteoritic crater, however, finds natural diamond formation via a subsolidus self-redox of a ferromagnesian carbonate during shock compression to 25–45 GPa and 800–900 °C without melting, fluid, and another reductant. The ability of carbonate to produce diamond by itself implies that diamond would be a very common mineral in the lower mantle where the carbonates are abundant and pressures and temperatures are sufficiently high.

Formation of natural diamonds requires the reduction of carbon to its bare elemental form, and pressures (P) greater than 5 GPa to cross the graphite–diamond transition boundary. In a study of shocked ferromagnesian carbonate at the Xiuyan impact crater, we found that the impact pressure–temperature (P-T) of 25–45 GPa and 800–900 °C were sufficient to decompose ankerite Ca(Fe²⁺,Mg)(CO₃)₂ to form diamond in the absence of another reductant. The carbonate self-reduced to diamond by concurrent oxidation of Fe²⁺ to Fe³⁺ to form a high-P polymorph of magnesioferrite, MgFe³⁺₂O₄. Discovery of the subsolidus carbonate self-reduction mechanism indicates that diamonds could be ubiquitously present as a dominant host for carbon in the Earth’s lower mantle.

Tiny timekeepers witnessing high-rate exhumation processes

Tectonic forces and surface erosion lead to the exhumation of rocks from the Earth's interior. Those rocks can be characterized by many variables including peak pressure and temperature, composition and exhumation duration. Among them, the duration of exhumation in different geological settings can vary by more than ten orders of magnitude (from hours to billion years). Constraining the duration is critical and often challenging in geological studies particularly for rapid magma ascent. Here, we show that the time information can be reconstructed using a simple combination of laser Raman spectroscopic data from mineral inclusions with mechanical solutions for viscous relaxation of the host. The application of our model to several representative geological settings yields best results for short events such as kimberlite magma ascent (less than ~4,500 hours) and a decompression lasting up to ~17 million years for high-pressure metamorphic rocks. This is the first precise time information obtained from direct microstructural observations applying a purely mechanical perspective. We show an unprecedented geological value of tiny mineral inclusions as timekeepers that contributes to a better understanding on the large-scale tectonic history and thus has significant implications for a new generation of geodynamic models.

Multiple pathways in pressure-induced phase transition of coesite

High-pressure single-crystal X-ray diffraction method with precise control of hydrostatic conditions, typically with helium or neon as the pressure-transmitting medium, has significantly changed our view on what happens with low-density silica phases under pressure. Coesite is a prototype material for pressure-induced amorphization. However, it was found to transform into a high-pressure octahedral (HPO) phase, or coesite-II and coesite-III. Given that the pressure is believed to be hydrostatic in two recent experiments, the different transformation pathways are striking. Based on molecular dynamic simulations with an ab initio parameterized potential, we reproduced all of the above experiments in three transformation pathways, including the one leading to an HPO phase. This octahedral phase has an oxygen hcp sublattice featuring 2 × 2 zigzag octahedral edge-sharing chains, however with some broken points (i.e., point defects). It transforms into α-PbO₂ phase when it is relaxed under further compression. We show that the HPO phase forms through a continuous rearrangement of the oxygen sublattice toward hcp arrangement. The high-pressure amorphous phases can be described by an fcc and hcp sublattice mixture.

Mineral inclusions in diamonds may be synchronous but not syngenetic

It is widely assumed that mineral inclusions and their host diamonds are ‘syngenetic' in origin, which means that they formed simultaneously and from the same chemical processes. Mineral inclusions that, instead, were formed earlier with respect to diamonds are termed protogenetic. However, minerals can have the same age as the diamonds in that they become enclosed in and isolated from any further isotopic exchange. But this is termed ‘synchronous' not ‘syngenetic'. Here we demonstrate conclusively the protogenesis of inclusions in diamonds, based upon data from an exceptional fragment of a diamond-bearing peridotite, its clinopyroxene and a gem-quality diamond. Clinopyroxenes in the xenolith had the same chemistry and crystallographic orientation as those for inclusions in the diamond. With our results with garnets, olivines and sulfides, we can state that a major portion of the mineral inclusions in non-coated, monocrystalline-lithospheric diamonds are protogenetic. Our discovery here presented has implications for all genetic aspects of diamond growth, including their ages.

Trapped mineral inclusions in diamonds give information on diamond crystallization and ages, under the assumption that they are syngenetic (formed simultaneously). Here, the authors show evidence that many mineral inclusions are protogenetic (formed at different times) thus undermining previous diamond ages.

Lattice dynamics of coesite

The lattice dynamics of coesite has been studied by a combination of diffuse x-ray scattering, inelastic x-ray scattering and ab initio lattice dynamics calculations. The combined technique gives access to the full lattice dynamics in the harmonic description and thus eventually provides detailed information on the elastic properties, the stability and metastability of crystalline systems. The experimentally validated calculation was used for the investigation of the eigenvectors, mode character and their contribution to the density of vibrational states. High-symmetry sections of the reciprocal space distribution of diffuse scattering and inelastic x-ray scattering spectra as well as the density of vibrational states and the dispersion relation are reported and compared to the calculation. A critical point at the zone boundary is found to contribute strongly to the main peak of the low-energy part in the density of vibrational states. Comparison with the most abundant SiO2 polymorph--α-quartz--reveals similarities and distinct differences in the low-energy vibrational properties.

Oriented graphite single-crystal inclusions in diamond

A crystallographic study of 13 specimens of diamonds with euhedral single-crystalline graphite inclusions in their centres is presented. All inclusions belong to the hexagonal graphite modification (space group P63/mmc; a 0 ≈ 2.46 Å, c 0 ≈ 6.70 Å) and are up to 300 μm in diameter. Comparison of the measured c-lattice parameters of the graphite crystals with lattice parameters of graphite at elevated isostatic pressures indicates remnant pressures of up to 2.6 GPa acting on the inclusions. All samples exhibit distinct orientation relations between graphite and diamond. In 12 samples the direction [001] of graphite (G[001]) approximately parallels one of the directions 〈111〉 of cubic diamond (D〈111〉). The largest deviation is about 4°. A further division of these 12 samples is: (a) The three G〈100〉 directions have angles of about 4° with three D〈110〉 directions. This orientation relation is observed in five samples. (b) The angle between the same directions is about 34° in six samples. (c) In one sample this angle is about 16°. In a 13th specimen G[001] approximately parallels one of the three D〈100〉. The deviation is of about 6°, and one of the G〈100〉 directions nearly parallels one of the D〈110〉. These orientation relations are analysed with a simplified application of the “coincidence site lattice (CSL)” concept. A 14th specimen is different to all others, as it exhibits a large (about 300 μm edge length) region with a sharp hexagonal borderline filled with a “patchwork” of tiny graphite “islands”. Several parallel lamellae of up to 5 μm thickness result in a hexagonal pyramidal form. This graphite inclusion is suggested to be protogenetic with respect to the diamond.

High precision in Raman frequency achieved using real-time calibration with a neon emission line: application to three-dimensional stress mapping observations

A three-dimensional (3D) Raman mapping system with a real-time calibration function was developed for detecting stress distributions in solid materials from subtle frequency shifts in Raman spectra. An atomic emission line of neon at 918.3 cm(-1) when excited at 514.5 nm was used as a wavenumber standard. An emission spectrum of neon and a Raman spectrum from a sample were introduced into a single polychromator using a bifurcated optical fiber. These two spectra were recorded simultaneously on a charge-coupled device (CCD) detector using double-track mode. Energy deviation induced by the fluctuation of laboratory temperature, etc., was removed effectively using the neon emission line. High stability during long measurements was achieved. By applying curve fitting, positions of the Raman line were determined with precision of about 0.05 cm(-1). The present system was applied to measurements of residual pressure around mineral inclusions in a natural diamond: 3D stress mapping was achieved.

Raman mapping of coesite inclusions in garnet from the Kokchetav Massif (Northern Kazakhstan)

Coesite inclusions occur in a wide range of lithologies and coesite is therefore a powerful ultrahigh-pressure (UHP) indicator. The transformation of coesite to quartz is evidenced by three optically well identifiable characteristics (e.g. palisade textures, radial crack patterns, polycrystalline quartz pseudomorphs). Under overpressure monomineralic coesite (on an optical basis), lacking the above transformation characteristics may survive. Raman micro-spectroscopy was applied on monomineralic coesite inclusions in garnet porphyroblasts from diamond-bearing garnet-clinozoisite-biotite gneisses of the Barchi-Kol area (Kokchetav Massif, Northern Kazakhstan). These coesite inclusions are euhedral and display a characteristic anisotropic hallo. However, Raman maps and separate spectra of these inclusions display shifted bands for coesite and quartz. Microscopically undetectable, quartz shows on the Raman map as a thin shell around coesite inclusion. Shift of the main coesite band allows to estimate their overpressure: coesite inclusions record 0-2.4 GPa in garnet and zircon. The quartz shell remains under lower pressure 0-1.6 GPa. The possible application of coesite and quartz Raman geobarometers for UHP metamorphic rocks is discussed.

Factors determining the stability, resolution, and precision of a conventional Raman spectrometer

We verified the performance of a conventional Raman spectrometer, which is composed of a 30 cm single polychromator, a Si based charge-coupled device (CCD) camera, and a holographic supernotch filter. For that purpose, the time change of the peak positions of Raman spectra of naphthalene and fluorescence spectra of ruby (Cr-doped Al(2)O(3)) were monitored continually. A time-dependent deviation composed of two components was observed: a monotonous drift up to 0.4 cm(-1) and a periodic oscillation with a range of 0.15 cm(-1). The former component was stabilized at approximately 2000 s after the CCD detector was cooled, indicating that incomplete refrigeration of the CCD detector induced the drift. The latter component synchronized with the periodic oscillation of the room temperature, indicating that thermal expansion or contraction of the whole apparatus induced this oscillation. The implemental deviation is reduced when measurements are conducted using a sufficiently cooled CCD detector at a constant room temperature. Moreover, the effect of the room temperature oscillation is lowered in a spectrum acquired over a duration that is longer than one cycle of this oscillation. Applying the least squares fitting method to carefully measured spectra enhanced the precision of the determination of the peak position to 0.05 cm(-1) using the spectrometer with pixel resolution of 1.5 cm(-1).

Near-infrared spectroscopic determination of salinity and internal pressure of fluid inclusions in minerals

A near-infrared (NIR) spectroscopic method is proposed to achieve the simultaneous determination of salinity and internal pressure of fluid inclusions in natural minerals. A combination band between the anti-symmetric stretching and bending vibrations of molecular water at approximately 5180 cm-1 was observed for standard salt solutions and natural minerals containing fluid inclusions with known salinities. A curve-fitting procedure was used to analyze the change in the band shape of the combination. Justification of the calibration was confirmed by observation of fluid inclusions in natural minerals whose salinities had already been determined using microthermometry. The detection limit of the present method is 1 NaCl-eq wt. %. The minimum size of fluid inclusions that produced well-resolved spectra was approximately 30 microm. This method was applied to assess micro fluid inclusions in a natural diamond with cubic growth habit (cuboid). The salinity and residual pressure of those fluid inclusions were estimated respectively as 4.4 wt. % NaCl-eq and 0.6-0.8 GPa. The present method is complementary to Raman microscopy and microthermometry for the determination of salinity in fluid inclusions of geological samples.

Crystal structure of a high-pressure/high-temperature phase of alumina by in situ X-ray diffraction

Alumina (alpha-Al(2)O(3)) has been widely used as a pressure calibrant in static high-pressure experiments and as a window material in dynamic shock-wave experiments; it is also a model material in ceramic science. So understanding its high-pressure stability and physical properties is crucial for interpreting such experimental data, and for testing theoretical calculations. Here we report an in situ X-ray diffraction study of alumina (doped with Cr(3+)) up to 136 GPa and 2,350 K. We observe a phase transformation that occurs above 96 GPa and at high temperatures. Rietveld full-profile refinements show that the high-pressure phase has the Rh(2)O(3) (II) (Pbcn) structure, consistent with theoretical predictions. This phase is structurally related to corundum, but the AlO(6) polyhedra are highly distorted, with the interatomic bond lengths ranging from 1.690 to 1.847 A at 113 GPa. Ruby luminescence spectra from Cr(3+) impurities within the quenched samples under ambient conditions show significant red shifts and broadening, consistent with the different local environments of chromium atoms in the high-pressure structure inferred from diffraction. Our results suggest that the ruby pressure scale needs to be re-examined in the high-pressure phase, and that shock-wave experiments using sapphire windows need to be re-evaluated.

Micro-Raman densimeter for CO2 inclusions in mantle-derived minerals

We investigated the applicability of Raman microprobe spectroscopy for determining the density of CO2 in fluid inclusions in minerals of mantle-derived xenolith samples. A separation (delta) between two Raman bands of CO2 due to Fermi resonance can be a reliable densimeter for CO2 fluid. The relationship between the density of CO2 (g/cm3) and delta (cm-1) can be expressed as: d = -0.03238697 delta 3 + 10.08428 delta 2 - 1046.189 delta + 36163.67. This equation was obtained from the Raman data on CO2 fluid with densities from 0.1 to 1.21 g/cm3, including super critical fluids at 58-59 degrees C. The delta value was constant with increasing temperature from room temperature to 200 degrees C. This indicates that the Raman densimeter is not affected by a possible rise in temperature, an artifact induced by the high flux of the incident laser. The minimum size of measurable inclusions is 1 micron, and the precision in the determination of delta is 0.1 cm-1, corresponding to 0.02 g/cm3 for inclusions of 1 micron in size. The precision can be better for larger inclusions. The micro-Raman densimeter can determine the density of CO2 fluid inclusions over a wide range. In particular, densities of gas and mixtures of gas and liquid phases, which cannot be measured by microthermometry, can be determined.