|
|
|
|||||||||||||||||
| JOURNAL HOME | HELP | CONTACT PUBLISHER | SUBSCRIBE | ARCHIVE | SEARCH | TABLE OF CONTENTS |

1 High Pressure Science and Engineering Center, Department of Physics, University of Nevada, Las Vegas, Nevada 89154, U.S.A.
2 Division of Geology and Planetary Sciences, California Institute of Technology, Pasadena, California 91125, U.S.A.
3 Department of Physics, New Mexico State University, Las Cruces, New Mexico 88003, U.S.A.
4 HPCAT, Advanced Photon Source, Argonne National Laboratory, Argonne, Illinois 60439, U.S.A.
5 Center for Condensed Matter Science and Technology, Research Academy of Science and Technology, Harbin 150080, China
6 Geophysical Laboratory, Carnegie Institution of Washington, Washington, D.C. 20015, U.S.A.
7 Physics Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, U.S.A.
Correspondence: * E-mail: olivert{at}physics.unlv.edu
| ABSTRACT |
|---|
|
|
|---|
Key Words: Crystal structure post-perovskite XRD data quantum mechanical calculation high-pressure studies Raman spectroscopy
| INTRODUCTION |
|---|
|
|
|---|
| THEORY |
|---|
|
|
|---|
,0 – CaIrO3-type) were relaxed to explore their relative energies and crystal chemistry and as initial structure models for the interpretation of the experimental diffraction data.
|
|
|
| EXPERIMENTAL METHODS |
|---|
|
|
|---|
In the first set of experiments (Sample 1), we changed the pressure several times between 100 and 120 GPa and heated at each pressure point to homogenize the stress conditions in the sample. In detail the loading and heating cycles during the first experiment were as follows: Initially, the sample was pressurized to 80 GPa and heated to ~ 2000 K for 300 s to convert the orthopyroxene sample to the perovskite structure. Then the sample was pressurized to between 110 and 120 GPa and heated to 2200 (±100) K for 600 s. The pressure after heating was 105 ± 5 GPa. Because ruby was incidentally heated by the laser beam upon alignment and underwent a structural transition upon heating (Lin et al. 2004), we used the volume of NaCl to calculate the pressure (Fei et al. 2007; Sata et al. 2002) and the given uncertainties are from the combined error in volume determination and in calibration of the NaCl-scale relative to the Pt-scale (Sata et al. 2002) and to the gold scale (Fei et al. 2007). The volume was obtained upon Rietveld refinement of the NaCl structure using the experimental diffraction pattern (RF2 = 0.04). While there is no absolute calibration of pressure at 100 GPa better than ±10 GPa (Holzapfel 2003), the volume of NaCl will allow for more accurate determination of pressure once a better pressure scale is established. Thereafter, we reduced and raised pressure between 98 and 115 ± 3 GPa two further times and heated the sample each time to 2000–2200 K.
In the second experiment (Sample 2) we raised the pressure initially to 80 GPa based on ruby fluorescence (Mao et al. 1990) and heated the sample for 300 s to 1800 K. At this pressure and temperature perovskite is the stable structure, and it is known that this phase has very limited solubility of water (Bolfan-Casanova et al. 2003; Ross et al. 2003). We assume that the water of the gedrite starting material has been mostly removed during this heating. Afterward, we raised the pressure to around 100 GPa and heated the sample to 2200 ± 50 K for 380 s. Diffraction data were collected after reducing temperature to 300 K, and the pressure was found to be 93 ± 2 GPa from the platinum pressure marker (Fei et al. 2007) and 95.0 ± 1 GPa from ruby (Mao et al. 1990).
X-ray diffraction patterns were collected at the 16ID-B undulator beamline at the High Pressure Collaborative Access Team (HPCAT), section 16 of the APS-ANL synchrotron, using a monochromatic beam of 32 or 36 keV energy focused with Kirkpatrick-Baez mirrors to 10 x 14 µm2 and a Mar345 image plate detector. We integrated the diffraction images and corrected for geometric distortions using the Fit2D software package (Hammersley et al. 1996).
Raman spectra were collected immediately after the diffraction experiment using the micro-Raman spectrometer of GSECARS at Sector 13 of the APS-ANL synchrotron operated with the 514.5 nm excitation line of an Ar-ion laser at 60 mW power focused to a spot of about 10 µm diameter in 180° backscattering geometry, a Spex 0.5 m monochromator with a single 1200 grooves/mm grating, 20 µm slit-width, and a liquid nitrogen cooled CCD camera detector from Oxford Instruments. We collected spectra at different locations of the silicate sample and in the pressure medium. There is no variation of the spectral features in the sample while the medium does not yield any Raman signal even at a distance of one diameter of the laser focal spot from the silicate. The recorded Raman spectra therefore belong to the silicate samples.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
|
2 of 14.5. More importantly, this fit gives a RF2 of 0.09 for the silicate-perovskite (Fig. 3a
angles between 3.5 and 5.5° indicates the existence of diffraction planes in the sample assembly that can neither be explained by perovskite nor by CaIrO3 type—post-perovskite or alumina. In contrast, kinked structures (with larger unit cells) allow for matching of these two low angle peaks. We note that marked diffraction at such low angles have been observed by Tateno et al. (2005) in the patterns of mixed perovskite and post-perovskite assemblies from laser heating experiments. Marked peaks at 12.8° and between 13 and 13.9 °2
do not correspond to reflections of perovskite but are matched by predicted peaks of the 3 x 1 type structure and even peaks of low intensity seem to match weaker features in the observed pattern. Subsequent Rietveld refinement (Fig. 3aBecause of the rather large instrumental contribution and the complex background in our data Rwp is not a good measure for judging the validity of a structure model. On the other hand, the modeled silicate structures are of low symmetry and strict overlap of reflections is not significant. Therefore, RF2 is a more trustworthy measure of the proximity between the modeled structures and the observation than Rwp in the present case.
In summary, the X-ray diffraction pattern in Figure 3a
shows features in ferrous aluminous MgSiO3 at 91–95 GPa that cannot be explained by perovskite, CaIrO3 type post-perovskite, or alumina. In contrast an intermediate 3 x 1 kinked post-perovskite structure provides a more consistent interpretation of the unexplained diffraction features especially at low 2
angles. The remaining mismatch between observed and predicted structure factors may be due to several reasons: The total fraction of the new phase is only a third of the fraction of perovskite. Thus, refinement based on whole profile fitting is more sensitive to slight mismatches in background and profiles of coexisting phases and this affects extraction of |F|2 values. Furthermore, the diffraction peaks show variable width and this, again, is difficult to fit if the signal to noise ratio is not very high. Finally it is possible that the observed pattern actually belongs to a mixture of similar kinked post-perovskite structures although we can preclude presence of alumina as free phase (see above).
The diffraction image of Sample 1 indicates noticeable texture of a fine-grained polycrystalline aggregate (Fig. 2b
). After integration we ran an automatic indexing of the observed sample Bragg peaks using the Jade 7.5+ package. To our surprise, we found a small number of orthorhombic and monoclinic cells of mostly equal cell dimensions, volumes, and angle. This suggests that the observed peaks belong to a single phase. We examined these candidate cells with respect to CaIrO3-type post-perovskite, our calculated post-perovskite polymorphs, and perovskite. The diffraction data from the present experiment are clearly different from perovskite and CaIrO3-type post-perovskite (Table 1
). However, all observed reflections including 12 distinct reflections below 12 °2
, can be indexed by the cell of the 2 x 1 type post-perovskite structure with dimensions 6.76 x 6.32 x 4.65, β = 99.2°, which is similar to the cell of the 2 x 1 structure as obtained from the ab initio calculations for pure MgSiO3 (Table 1
).
We calculated the powder diffraction pattern of 2 x 1 kinked post-perovskite based on the atomic positions of the ab-initio calculated structure and compared it to the observed pattern. Subsequent Rietveld refinement focuses on matching the apparent texture by spherical harmonics up to second order (Fig. 3b
) using the method by von Dreele (1997). Along with the profile weighted refinement parameters (see caption Fig. 3b
) this procedure shows a match of calculated and observed structure factor moduli of an RF2 of 0.25. As for Sample 1, we find that the diffraction pattern of this sample is incompatible with perovskite, CaIrO3-type post-perovskite, and alumina, while the observed pattern in this case is consistent with a 2 x 1 kinked post-perovskite. Both Samples 1 and 2 show well-resolved low angle diffraction peaks in the range of 3.5 to 6 in °2
. This observation is incompatible with perovskite and CaIrO3-type post-perovskite while it is consistent with kinked post-perovskite.
Upon subsequent extended heating cycles of Sample 1 (see experimental section) the size of crystallites increased, and the diffraction pattern changes from that of a textured powder to that of a few larger crystallites that give rise to several individual Bragg reflections at different locations of the image plate detector. The textural pattern shows that during the high P-T experiment the silicate sample had passed from a regime of formation of randomly oriented crystalline seeds into a regime of collective recrystallization and re-orientation. This state is the result of repetitive heating under conditions of small temperature gradients, due to the use of a comparatively thick envelope of pressure medium of apparently lower thermal conductivity than the sample itself (Kiefer and Duffy 2005) and due to direct heating of the silicate sample rather than using mixed-in metallic absorbers. We observe 52 reflections that belong to various crystallites. Out of this set we select seven reflections which seem to belong to one crystallite since they establish two intersecting lattice trajectories which can be mapped onto each other (Fig. 1 Appendix
). Because our diffraction patterns were collected at a fixed angular position, we sample only a two-dimensional slice of the reciprocal space. It is the narrow backside aperture of the Mao-Bell cell that does not permit sampling of further angular settings. Because of the two-dimensionality of the observable portion of reciprocal space we are constrained to compare the 2
angles of these reflections to silicate-perovskite, CaIrO3-post-perovskite, and the kinked post-perovskite structures (Table 1 Appendix
). This comparison is presented in table 1 of the Appendix. It shows consistence of the observed set of 2
values with various n x 1 type kinked post-perovskites while it is inconsistent with either perovskite or CaIrO3-type post-perovskite. This is discussed in more detail in the Appendix.
|
|
| RAMAN SPECTRA |
|---|
|
|
|---|
|
|
Figure 4
, lower part, shows the Raman spectrum of Sample 1 (textured) at 115 ± 5 GPa. The spectrum also shows a strong broad, structured peak between 700 and 1000 cm–1 while there are sharp and intense peaks between 250 and 500 cm–1. For this textured sample intensities are likely controlled by crystallite orientation (see Figs. 2b
and 3b
caption), while the spectrum of Sample 2 gives an average over a powder.
Sample 1 shows at least 22 Raman peaks (Table 3
) none of which adheres to perovskite, whereas the CaIrO3-structure exhibits only 12 Raman active modes according to the Bilbao Crystallographic Server (Kroumova et al. 2003). The broad peaks between 700 and 1000 cm–1 are again consistent with the concept of Raman-active modes involving silicate-octahedral breathing motions. The more complex structure of the Raman signal in this range (Table 3
) showing at least seven peaks rather than three (Caracas and Cohen 2006) indicates a symmetry reduction that may be consistent with kinked post-perovskite structures.
The larger than expected number of Raman peaks suggests a lower symmetry and larger unit cell than that of the CaIrO3-type phase, not inconsistent with the 2 x 1 and the 3 x 1 structure, which have 42 and 54 Raman active modes, respectively (Kroumova et al. 2003).
The spectra from Samples 1 and 2 are similar but not identical. However, their composition and, according to our analysis of the X-ray diffraction data, their structures are different also while sharing a common construction principle (Figs. 1b and 1c
). Both spectra share the occurrence of Raman peaks above 800 cm–1, which are not expected for perovskite (Serghiou et al. 1998; Karki et al. 2000; Caracas and Cohen 2006), but may be a signature of post-perovskite structures. The larger half width of these high-energy Raman peaks in comparison to the lower energy peaks in the same spectra may be due to disorder and to Al substitution on the Si sites, but this remains speculative as long as we are lacking definitive structure solutions for the observed phases.
In summary, our ab-initio calculations have reproduced the n x m type kinked post-perovskite structures found in metadynamics calculations by Oganov et al. (2005), and we extended our calculations to further polymorphs of this type. We have shown that kinked post-perovskites correspond to local minima in the energy hypersurface at 0 K. Above 100 GPa, the calculated enthalpy of the kinked structures is lower than that of perovskite but remains higher than that of CaIrO3-type post-perovskite. Thus kinked type structures are not stable at 0 K, but they may be stabilized at high temperature due to thermal contributions to the Gibbs free energy that are neglected in our calculations. The energetic preference for the CaIrO3-type at 0 K holds for pure MgSiO3 as well as for charge-coupled aluminum substitutions in MgSiO3. Both X-ray diffraction data and Raman spectra on aluminous magnesium metasilicate show features that are not expected for perovskite or CaIrO3-type post-perovskite. Our Raman spectra show similarities with the predicted spectrum for CaIrO3-type post-perovskite. However, the large number of Raman peaks suggests lower symmetry and larger unit cell size than that of CaIrO3-type post-perovskite. These observations from Raman and X-ray diffraction data and the comparison to the diffraction patterns derived from our ab-initio calculations suggest that the structures in our experiments are more consistent with kinked structures of the type n x 1 than with perovskite or CaIrO3-type post-perovskite.
Our calculations show that at 0 K and for pure and aluminous MgSiO3, kinked structures become more stable than perovskite at Mbar pressures but remain less stable than CaIrO3-type post-perovskite, consistent with previous predictions from metadynamics simulations (Oganov et al. 2005). In the absence of detailed knowledge of the activation barriers between the local minima that belong to different n x m structures, it is not possible to assess the thermal stability of the kinked structures from theory alone. Thus, it is not clear if kinked structures are stable or metastable for sufficiently long time to be relevant for the evolution of the Earths lowermost mantle or if they are remnants of an incomplete phase transition. However, independent of the possibility of a natural occurrence of kinked post-perovskite structures in the deep mantle our results support the previous conjecture that the perovskite
post-perovskite structural transition proceeds via intermediate kinked structures rather than by homogeneously straining the perovskite structure.
| APPENDIX |
|---|
|
|
|---|
values of these reflections are compared to calculated values of CaIrO3-type post-perovskite and various kinked post-perovskites. Some n x 1 type kinked phases match the observed 2x values but CaIrO3-type post-perovskite does not (Table 1
values) does not improve the match to a satisfactory degree (Table 1
| ACKNOWLEDGMENTS |
|---|
|
|
|---|
| Footnotes |
|---|
Open Access: Thanks to the authors generous funding this article is available to all online (http://ammin.geoscienceworld.org) and the MSA web site, which also has information about the MSA Open Access policy at http://www.minsocam.org/MSA/ammin/e-pub_policy.htm (scroll down). MANUSCRIPT HANDLED BY PRZEMYSLAW DERA
MANUSCRIPT RECEIVED June 30, 2006; MANUSCRIPT ACCEPTED November 12, 2007
| REFERENCES CITED |
|---|
|
|
|---|
Blöchl, P.E. (1994) Projector augmented-wave method. Physical Review B, 50, 17953–17979.[Medline]
Bolfan-Casanova, N., Keppler, H., and Rubie, D.C. (2003) Water partitioning at 660 km depth and evidence for very low water solubility in magnesium silicate perovskite. Geophysical Research Letters, 30, Art. No. 1905.
Caracas, R. and Cohen, R.E. (2006) Theoretical determination of the Raman spectra of MgSiO3 perovskite and post-perovskite at high pressure. Geophysical Research Letters, 33, Art. No. L12S05.
Fei, Y., Ricolleau, A., Frank, M., Mibe, K., Shen, G., and Prakapenka, V. (2007) Toward an internally consistent pressure scale. Proceedings of the National Academy of Sciences of the United States of America, 114, 9182–9184.
Hammersley, A.P., Svensson, S.O., Hanfland, M., Fitch, A.N., and Hausermann, D. (1996) Two-dimensional detector software: From real detector to idealised image or two-theta scan. High Pressure Research, 14, 235–248.[CrossRef][ISI]
Helmberger, D., Lay, T., Ni, S., and Gurnis, M. (2005) Deep mantle structure and the postperovskite phase transition. Proceedings of the National Academy of Sciences of the United States of America, 102, 17257–17263.
Hohenberg, P. and Kohn, W. (1964) Inhomogeneous electron gas. Physical Review B, 136, 864–871.
Holzapfel, W.B. (2003) Refinement of the ruby luminescence pressure scale. Journal of Applied Physics, 93, 1813–1818.[CrossRef][ISI]
Kameyama, M. and Yuen, D.A. (2006) 3-D convection studies on the thermal state in the lower mantle with post-perovskite phase transition. Geophysical Research Letters, 33, Art. No. L12S10.
Karki, B.B., Wentzcovitch, R.M., de Gironcoli, S., and Baroni, S. (2000) Ab initio lattice dynamics of MgSiO3 perovskite at high pressure. Physical Review B, 62, 14750–14756.
Kiefer, B. and Duffy, T.S. (2005) Finite element simulations of the laser-heated diamond-anvil cell. Journal of Applied Physics, 97, 114902.[CrossRef]
Kresse, G. and Furthmüller, J. (1996) Efficient iterative schemes for ab initio total-energy calculations using a plane wave basis set. Physical Review B, 54, 11169–11186.[Medline]
Kresse, G. and Hafner, J. (1993) Ab-initio molecular-dynamics for open shell transition metals. Physical Review B, 48, 13115–13118.[Medline]
Kresse, G. and Joubert, D. (1999) From ultrasoft pseudopotentials to the projector augmented-wave method. Physical Review B, 59, 1758–1775.
Kroumova, E., Aroyo, M.I., Perez-Mato, J.M., Kirov, A., Capillas, C., Ivantchev, S., and Wondratschek, H. (2003) Bilbao crystallographic server: Useful databases and tools for phase-transition studies. Phase Transitions, 76, 155–170.[CrossRef][ISI]
Lin, J.F., Degtyareva, O., Prewitt, C.T., Dera, P., Sata, N., Gregoryanz, E., Mao, H.K., and Hemley, R.J. (2004) Crystal structure of a high-pressure/high-temperature phase of alumina by in situ X-ray diffraction. Nature Materials, 3, 389–393.[CrossRef][ISI][Medline]
Mao, H.K., Wu, Y., Chen, L.C., Shu, J.F., and Jephcoat, A.P. (1990) Static compression of iron to 300 GPa and Fe0.8Ni0.2 alloy to 260 GPa: Implications for composition of the core. Journal of Geophysical Research, 95, 21737–21742.[CrossRef]
Monkhorst, M.J. and Pack, J.D. (1976) Special points for Brillouin-zone integrations. Physics Reviews B, 13, 5188–5192.[CrossRef]
Murakami, M., Hirose, K., Kawamura, K., Sata, N., and Ohishi, Y. (2004) Post-perovskite phase transition in MgSiO3. Science, 304, 855–858.
Nakagawa, T. and Tackley, P.J. (2004) Effects of a perovskite-post perovskite phase change near core-mantle boundary in compressible mantle convection. Geophysical Research Letters, 31, 4–5.
Oganov, A.R. and Ono, S. (2004) Theoretical and experimental evidence for a post-perovskite phase of MgSiO3 in Earths D layer. Nature, 430, 445–448.[CrossRef][Medline]
Oganov, A.R., Martonak, R., Laio, A., Raiteri, P., and Parinello, M. (2005) Anisotropy of Earths D layer and stacking faults in the MgSiO3 post-perovskite phase. Nature, 438, 1142–1144.[CrossRef][Medline]
Perdew, J.P., Burke, K., and Enzerhof, M. (1996) Generalized gradient approximation made simple. Physical Review Letters, 77, 3865–3868.[CrossRef][ISI][Medline]
Ross, N.L., Gibbs, G.V., and Rosso, K.M. (2003) Potential docking sites and positions of hydrogen in high-pressure silicates. American Mineralogist, 88, 1452–1459.
Sata, N., Shen, G.Y., Rivers, M.L., and Sutton, S.R. (2002) Pressure-volume equation of state of the high-pressure B2 phase of NaCl. Physical Review B, 65, no. 104114.
Serghiou, G., Zerr, A., and Boehler, R. (1998) (Mg,Fe)SiO3-perovskite stability under lower mantle conditions. Science, 280, 2093–2095.
Tateno, S., Hirose, K., Sata, N., and Ohishi, Y. (2005) Phase relations Mg3Al2Si3O12 to 180 GPa: Effect of Al on post-perovskite phase transition. Geophysical Research Letters, 32, Art. No. L15306.
Tsuchiya, T., Tsuchiya, J., Umemoto, K., and Wentzcovitch, R.M. (2004a) Phase transition in MgSiO3 perovskite in the Earths lower mantle. Earth and Planetary Science Letters, 224, 241–248.[CrossRef][ISI][GeoRef]
von Dreele, R.B. (1997) Quantitative texture analysis by Rietveld refinement. Journal of Applied Crystallography, 30, 517–525.[CrossRef][ISI]
| ||||||||||||||||||||||||||||||||||||||||||||||||||||
| JOURNAL HOME | HELP | CONTACT PUBLISHER | SUBSCRIBE | ARCHIVE | SEARCH | TABLE OF CONTENTS |