To save content items to your account,
please confirm that you agree to abide by our usage policies.
If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your account.
Find out more about saving content to .
To save content items to your Kindle, first ensure no-reply@cambridge.org
is added to your Approved Personal Document E-mail List under your Personal Document Settings
on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part
of your Kindle email address below.
Find out more about saving to your Kindle.
Note you can select to save to either the @free.kindle.com or @kindle.com variations.
‘@free.kindle.com’ emails are free but can only be saved to your device when it is connected to wi-fi.
‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.
A statistical method that uses a generalized linear model is presented to provide best-fit phase boundaries to experimentally determined phase-diagram data. The experimental determination of the exact locus of a phase boundary is inherently uncertain because of increasing difficulty in determining the presence of one phase and/or the absence of another as the phase boundary is approached from either side. We present a logistic model, which states that a phase transformation α → β can be expressed in terms of the probability of observing β at a given set of measured state variables X and Y such that the phase boundary is defined as P(β | X, Y) = 0.5. As an example, high-pressure and high-temperature data for the melting curve of platinum are analyzed with this method.
Introduction
Phase diagrams are determined by recording the state of a material at a series of thermodynamic conditions, for example, whether the sample is molten or crystalline as a function of temperature, composition, and pressure. Although phase boundaries are thermodynamically well defined, the experimental observations that constrain phase diagrams often display significant scatter about the equilibrium boundaries. The scatter may be evident either within a single set of measurements or between distinct experiments from different laboratories.
We extend a recent comparison of the abilities of different pair-potential models for SiO2 to reproduce the IR spectrum of the glassy state (a short-time dynamic property) to other properties. The comparison is extended to include static (thermodynamic and structural) and long-time dynamic (diffusivity) properties. As with the IR spectrum, none of the existing potentials succeeds in correctly locating the temperature of the density maximum Tp(max)- Although the modified Matsui potential, which best represented the IR spectrum, also performs well in describing Tp(max) and the intertetrahedral bond angle, it seriously overestimates the diffusivity and hence incorrectly represents the shape of the potential far from the minimum. The best performance across the range of properties examined is obtained with the van Beest–Kramer–van Santen potential. Some evidence for a higher-order phase transition below the density maximum is given.
Introduction
It is now two decades since it was found (unexpectedly to many) that key properties of liquid and glassy SiO2 could be reproduced semiquantitatively by molecular-dynamics (MD) computer simulations in which the Si and the O components are treated as simple ions with full formal charges [1]. The internal energy at 300 K was obtained to within 3%, and the qualitatively unique features of (1) high liquid compressibility contrasting with low liquid expansivity, (2) anomalous increase of diffusivity with pressure, and (3) large pressure-induced glass densification were also demonstrated.
The career and achievements of Yoshito Matsui have emphasized the attainment of accuracy in experimental studies. To honor him, this contribution presents an analysis of several sources of systematic errors in powder-diffraction experiments and evaluates their effects on the results. In particular, the error arising from uncorrected axial divergence in an energy-dispersive experiment may be as large as 1%. Smaller errors of 0.3% may arise from sample positioning errors, and errors as large as 2% can arise from the effects of nonhydrostatic stress.
Introduction
Although an x-ray diffraction experiment may have very high precision, there are many factors that need to be controlled if an accurate result is to be obtained. This paper discusses the potential systematic errors arising from peak asymmetry that are due to axial divergence, sample positioning errors, and nonhydrostatic stress distributions.
Axial Divergence
Axial divergence, the effect that causes peak asymmetry for diffraction maxima near 0° or 180° 2θ in powder-diffraction patterns, arises from the elliptical shape of the intersection of diffraction cones with the cylinder that describes the opening of the detector as it is swept through 2θ. The description of this effect and a method for calculating the resulting profile have been presented by Finger et al. [1]. This correction can be applied directly for data measured with monochromatic radiation and a diffractometer. Axial divergence should not affect monochromatic powder patterns measured with flat-plate detectors, such as imaging plates or CCD detectors. Such patterns are usually processed to convert the circular rings on the flat two-dimensional detector into a pseudo-one-dimensional pattern by dividing the pattern into a number of wedges and integrating each piece [2].
If you look at the cosmic abundance of elements, oxygen, carbon, neon, nitrogen, magnesium, and silicon are the top six (apart from hydrogen and helium); most of them are indeed the main constituents of rock-forming minerals. If the goal of modern mineralogy is to understand classes of materials on an atomic scale, this is also the aim of modern condensed-matter physics. Therefore it is natural that condensed-matter physics should be applied to mineralogy and earth sciences, just as particle physics is a key tool in understanding cosmology.
This book is intended to give a state-of-the-art description of intensive interactions between geophysics and condensed-matter physics that recent years have witnessed. Although you might assume that, given the maturity of the solid-state physics, the crystalline structures of materials must have been readily understood in atomistic, nonempirical terms, this rather naive expectation has not been satisfied until quite recently.
Although traditional mineralogical principles had remained largely empirical, a few individuals pursued the idea that modern theoretical solid-state physics must become the foundation of mineralogy and geophysical sciences. One of the pioneers in this respect is Professor Yoshito Matsui in Japan. Following the tradition of Goldschmidt, he realized the importance of identifying and characterizing the underlying physics that controls geochemical, geophysical, and geological phenomena over the entire range of pressure and temperature relevant to the Earth. In so doing, he, in collaboration with Dr. Eiji Ito, directed the development of new experimental high-pressure facilities that would be required for understanding the rich and sometimes unexpected behaviour of minerals under conditions found deep within the Earth.
The diversity in crystal structures of materials of interest in condensed-matter physics and mineral sciences is reviewed from the viewpoint that they are an outcome of multiple ways in which building blocks (polyhedra, clusters, chains, layers, etc.) can be arranged in ambient or high pressures. Opportunities for interesting electronic properties expected from unconventional as well as conventional crystal structures are outlined, which may even include electron correlation engineering when those properties arise from the electron correlation.
Introduction
Although quantum mechanics laid a solid foundation for condensed-matter physics earlier this century, we are now witnessing a new era in that branch of physics in which we begin to understand crystal structures in nonempirical ways. As stressed in the Preface of this volume, this has been particularly fruitful in understanding rock-forming minerals. It is indeed fascinating if we can predict crystal structures of minerals or other materials in general from a knowledge of the chemical composition alone. This also opens up a way to search for novel crystal structures, either in denser or open-structured forms.
The purpose of this chapter is twofold. One is to take silica as an example for elaborating on various crystal structures, obtained from a nonempirical computer-simulation study. Silica comprises the two most abundant elements on the Earth, i.e., silicon and oxygen. Silica, along with silicates, is an important ingredient in rock-forming minerals. Crystal-to-crystal phase transformations, as well as pressure-induced amorphisation, are described.
Periclase, or MgO, is a simple ionic mineral, but one in which many-body interactions are important. The thermal equation of state, elasticity, melting, thermal conductivity, and diffusivity have been studied by use of first-principles methods, i.e., from fundamental physics. There is generally good agreement with experiment. Studies of MgO provide fundamental information on the high-pressure behavior of minerals in the deep Earth.
MgO, periclase, is the simplest oxide and has been a subject of intense experimental and theoretical study. Oxides and silicates make up the bulk of the Earth's mantle and crust, and thus it is important to understand and predict their behavior. The behavior of MgO, as the prototypical oxide, is the key to understanding mineral and rock behavior in the bulk of the Earth. An important feature of MgO is the nonrigid behavior of the oxygen O2– ion, which makes the interactions not describable by pairwise interactions. All other more complex oxides share this feature and add additional complications as well. MgO is simple in that it is an ionic material with no solid-state phase transitions until over 500 GPa according to the best computations. If we understand MgO we will not immediately understand all other oxides and silicates, but if we cannot understand MgO, we cannot understand any other oxide or silicate.
There is an enormous literature on MgO, and a comprehensive review is impossible here. Instead the main thrust is on the underlying fundamental physics of MgO. First-principles methods, as opposed to empirical methods, are emphasized. Properties that are addressed range from the thermal equation of state and elasticity to properties of defects, surfaces, and impurities.
Diffusivity measurements in silicate liquids at high pressures have now been obtained up to 15 GPa and 2800 K by use of a multianvil apparatus. Most striking from these investigations is that oxygen and silicon self-diffusivities in Na2Si4O9 liquid increase continuously as a function of pressure from 2.5 to 15 GPa. According to the Eyring relation, this would suggest that the viscosity of the liquid decreases over the same pressure range. However, as the composition of the liquid becomes more polymerized, approaching metaluminous albite, diffusivities reach a maximum at pressures below 10 GPa and decrease with further increases in pressure. These results demonstrate the importance of considering compositional parameters, such as degree of polymerization and Al/Al + Si ratio, as well as pressure, when we are modeling magmatic processes in the Earth's deep interior.
Introduction
Magmatic processes within the Earth are controlled to a large extent by transport properties of silicate liquids. For example, both rates of magma ascent and rates of crystal settling during fractionation processes are controlled by the viscosity of silicate melts. In addition, ionic diffusion in silicate liquids controls rates of magma mixing, homogenization, and chemical equilibration and is therefore an important controlling parameter at all stages of magma evolution, from initial partial melting at depth to eruption or emplacement at or near the Earth's surface. Because a knowledge of the transport properties of silicate melts is essential for quantifying magmatic processes in the Earth, the viscosities of such melts have been studied extensively as functions of temperature and chemical composition at 1 bar with a range of experimental techniques [1].
Cristobalite was compressed at room temperature up to 30 GPa under quasi-hydrostatic and nonhydrostatic conditions. The structure of the high-pressure phase differs completely, depending on the hydrostaticity of applied pressure. Under quasi-hydrostatic conditions, a structure very similar to that of stishovite was formed above ∼20 GPa, which can be quenched to ambient conditions. When this phase was heated to several hundred degrees centigrade under pressure, normal stishovite was formed. Under nonhydrostatic conditions, on the other hand, cristobalite transformed into the unidentified X-I phase at ∼15 GPa, in accordance with previous study. This X-I phase transforms into stishovite on heating to above 1000°C. The apparent compressibility of cristobalite also differs considerably, depending on the nature of applied pressure. The nature of the high-pressure phases formed by the room-temperature compression of cristobalite is discussed.
Introduction
Silica is one of the most intensively studied materials under a wide range of pressure and temperature conditions because of its importance in earth science as well as in material science. In spite of its simple chemical formula, many different polymorphs are formed, depending on the pressure and the temperature conditions. Moreover, recent room-temperature compression studies of various polymorphs clarified the formation of variety of metastable phases [1–5], and the situation is complicated. It is likely that the formation of some of the metastable phases is related to the nonhydrostatic nature of applied pressure. Here we report a new example of the formation of completely different high-pressure structures that depend on the hydrostaticity of applied pressure.
Observations of order–disorder phenomena at high pressure in hydrous phases are reinterpreted with the results of Rietveld analysis and neutron-diffraction data. The reported partial amorphization of the hydrogen sublattice in β-Co(OD)2 at 11.2 GPa was not confirmed in powder-diffraction data collected with the Paris-Edinburgh cell to 15.5 GPa. The diffraction data, and perhaps the spectroscopic data on which the observations of amorphization are based, are consistent with an increase in the H…H repulsion with pressure. The structural consequences of competition between H…H repulsion and H-bond (O—H…O) formation is observed in the M(OH)2 compounds in general. It is also observed in the dense high-pressure phases recovered from high-pressure synthetic experiments.
Introduction
The hydrogen bond (X—H… Y) is one of the most studied bond geometries in the mineralogical, biological, and solid-state organic chemical communities [1, 2]. For nonmineral and mineral structures alike, the published literature, consisting mainly of crystal-structure determinations at ambient pressure, provides a means to study the bond as donor (X) and acceptor (Y) vary over a variety of structures and chemistries [3, 4]. The secondary environment, however, is important in considering the effects of structure on H-bond geometry [5]; in many cases gross changes in this environment from one structure type to the next make it difficult to separate the effects of the relatively weak H bonding from the steric effects because of the framework making up the remainder of the structure [5].
The interplay among geosciences, mineralogy, and condensed-matter physics is briefly overviewed in the context of the chapters contained in this volume. The chapters report current developments along many lines of inquiry introduced by one of the leading figures, Yoshito Matsui, for which we start with a brief history. The diverse body of work summarized in the book includes ranges from recent advances in theoretical condensed-matter physics, crystallography and crystal chemistry, high-pressure transformations, and the behaviour of high-temperature melts.
Introduction
It has long been recognized that atomic interactions are key in understanding large-scale geological phenomena, which is an approach dating at least as far back as the days of Goldschmidt [1]. Conversely, a study of the materials that comprise the planets can tell us much about fundamental physics and chemistry. This line of approach is subsequently exemplified by Pauling's development of the theory of chemical bonding from the structural studies of minerals [2].
Among many facets that revolutionized modern earth science are plate tectonic theory, discoveries in planetary astronomy, and extraordinary technical advances that provide windows on our planet's deep interior. Modern condensed-matter physics, on the other hand, has acquired much impetus from the discovery of new phenomena, the development of new experimental methods, and vast improvements in computational methods. Although physics and mineralogy have continued largely on distinct and disparate trajectories, many recent advances in both earth science and physics stem from science done at the crossroads of the two fields.
By
Z. Fang, Angstrom Technology Partnership,
K. Terakura, National Institute for Advanced Interdisciplinary Research,
H. Sawada, Angstrom Technology Partnership,
I. Solovyev, Angstrom Technology Partnership,
T. Miyazaki, National Research Institute for Metals
The phase stability of transition-metal monoxides MnO and FeO under ultrahigh pressure, which reaches the range in the Earth's lower mantle, was studied with the first-principles calculations based on density function theory. The plane-wave basis pseudopotential method was used to perform the structure optimization efficiently, and the electron–electron interaction was treated by the generalized gradient approximation (GGA) supplemented by the LDA + U method (LDA is local-density approximation). Two related structures, normal NiAs (nB8) and inverse NiAs (iB8) types, are emphasized. Our results predict that the high-pressure phase of MnO should take the nB8 structure rather than the CsCl (B2) structure and that a metastable nonmagnetic Bl structure can be realized for MnO in the intermediate pressure range. A very unique iB8 structure rather than the nB8 structure is predicted as the high-pressure phase of FeO, although no materials have ever been known to take the iB8 structure. The novel feature of the iB8 FeO is that the system should be a band insulator in the antiferromagnetic state and that the existence of a bandgap leads to special stability of the phase. The larger c/a ratios for both nB8 MnO and iB8 FeO were explained based on our analysis of the cation radius-anion radius ratios versus c/a for series of similar materials.
Introduction
The high-pressure phases of metal monoxides [including alkaline-earth-metal monoxides, transition-metal monoxides (TMMOs), etc.], with the rock-salt (Bl) structure at normal pressure and room temperature, are important for both condensed-matter physics and earth science because their crystal structure is simple and moreover MgO and FeO are considered to be important constituents of the Earth's deep mantle.
To understand the anisotropic nature of the shock-induced phase transition of TiO2 (rutile), the transition to postrutile phases under isotropic and several stress-field conditions was studied by molecular-dynamics calculations. The rutile structure was shown to transform to the fluorite structure by a displacive mechanism under the isotropic compression within the time scale of the shock transition. Calculations with variable stress fields showed anisotropic behavior. The rutile–fluorite transition occurred smoothly with [100] compression of rutile. In the case of the [110] compression, the rutile structure transformed to a twinned fluorite structure. However, the rutile structure transformed to the CaCl2-type structure instead of the fluorite structure under the [001] compression. These results were in good agreement with actual shock experiments.
Introduction
Shock Compression Method for Solids
The dynamic compression method using shock wave has been utilized for highpressure research on solids. The method can easily generate high pressures to more than 100 GPa. In principle, the nature of the dynamic compression is understood by hydrodynamical considerations [1]: When a shock wave travels in solids with a supersonic speed (U) and accelerates particles to a particle velocity (w), the shock front is generated as a discontinuity boundary of pressure and density. In real solids, the discontinuity is observed as a very steep change of pressure in solids and the transition interval is only several nanoseconds.