Melting of iron at the physical conditions of the Earth's core

Seismological data can yield physical properties of the Earth's core, such as its size and seismic anisotropy. A well-constrained iron phase diagram, however, is essential to determine the temperatures at core boundaries and the crystal structure of the solid inner core. To date, the iron phase diagram at high pressure has been investigated experimentally through both laser-heated diamond-anvil cell and shock-compression techniques, as well as through theoretical calculations. Despite these contributions, a consensus on the melt line or the high-pressure, high-temperature phase of iron is lacking. Here we report new and re-analysed sound velocity measurements of shock-compressed iron at Earth-core conditions. We show that melting starts at 225 +/- 3 GPa (5,100 +/- 500 K) and is complete at 260 +/- 3 GPa (6,100 +/- 500 K), both on the Hugoniot curve-the locus of shock-compressed states. This new melting pressure is lower than previously reported, and we find no evidence for a previously reported solid-solid phase transition on the Hugoniot curve near 200 GPa (ref. 16).     

The plastic deformation of iron at pressures of the Earth's inner core

Soon after the discovery of seismic anisotropy in the Earth's inner core, it was suggested that crystal alignment attained during deformation might be responsible. Since then, several other mechanisms have been proposed to account for the observed anisotropy, but the lack of deformation experiments performed at the extreme pressure conditions corresponding to the solid inner core has limited our ability to determine which deformation mechanism applies to this region of the Earth. Here we determine directly the elastic and plastic deformation mechanism of iron at pressures of the Earth's core, from synchrotron X-ray diffraction measurements of iron, under imposed axial stress, in diamond-anvil cells. The epsilon-iron (hexagonally close packed) crystals display strong preferred orientation, with c-axes parallel to the axis of the diamond-anvil cell. Polycrystal plasticity theory predicts an alignment of c-axes parallel to the compression direction as a result of basal slip, if basal slip is either the primary or a secondary slip system. The experiments provide direct observations of deformation mechanisms that occur in the Earth's inner core, and introduce a method for investigating, within the laboratory, the rheology of materials at extreme pressures.     

Elasticity of iron at the temperature of the Earth's inner core.

Seismological body-wave and free-oscillation studies of the Earth's solid inner core have revealed that compressional waves traverse the inner core faster along near-polar paths than in the equatorial plane. Studies have also documented local deviations from this first-order pattern of anisotropy on length scales ranging from 1 to 1,000 km (refs 3, 4). These observations, together with reports of the differential rotation of the inner core, have generated considerable interest in the physical state and dynamics of the inner core, and in the structure and elasticity of its main constituent, iron, at appropriate conditions of pressure and temperature. Here we report first-principles calculations of the structure and elasticity of dense hexagonal close-packed (h.c.p.) iron at high temperatures. We find that the axial ratio c/a of h.c.p. iron increases substantially with increasing temperature, reaching a value of nearly 1.7 at a temperature of 5,700 K, where aggregate bulk and shear moduli match those of the inner core. As a consequence of the increasing c/a ratio, we have found that the single-crystal longitudinal anisotropy of h.c.p. iron at high temperature has the opposite sense from that at low temperature. By combining our results with a simple model of polycrystalline texture in the inner core, in which basal planes are partially aligned with the rotation axis, we can account for seismological observations of inner-core anisotropy.     

地核物质的状态方程和参数

 解释了3种地核物质状态方程的异同,验证了Vinet et al. EOS(equation of state)比Shanker et al.EOS更接近Stacy EOS的结论,后者与高压极限的热动力学一致、而且与来源于地震波的数据完全匹配;同时与地震波测量数据进行了对照;给出了内核边界压力下物质摩尔体积随温度的变化函数,预言了密度亏损的温度变化.这对于揭示地球深部物质特性、研究地核内的轻物质组分以及深部动力学原理具有重要意义.     

Possible thermal and chemical stabilization of body-centred-cubic iron in the Earth's core

The nature of the stable phase of iron in the Earth's solid inner core is still highly controversial. Laboratory experiments suggest the possibility of an uncharacterized phase transformation in iron at core conditions and seismological observations have indicated the possible presence of complex, inner-core layering. Theoretical studies currently suggest that the hexagonal close packed (h.c.p.) phase of iron is stable at core pressures and that the body centred cubic (b.c.c.) phase of iron becomes elastically unstable at high pressure. In other h.c.p. metals, however, a high-pressure b.c.c. form has been found to become stabilized at high temperature. We report here a quantum mechanical study of b.c.c.-iron able to model its behaviour at core temperatures as well as pressures, using ab initio molecular dynamics free-energy calculations. We find that b.c.c.-iron indeed becomes entropically stabilized at core temperatures, but in its pure state h.c.p.-iron still remains thermodynamically more favourable. The inner core, however, is not pure iron, and our calculations indicate that the b.c.c. phase will be stabilized with respect to the h.c.p. phase by sulphur or silicon impurities in the core. Consequently, a b.c.c.-structured alloy may be a strong candidate for explaining the observed seismic complexity of the inner core.     

Stability of the body-centred-cubic phase of iron in the Earth's inner core

Iron is thought to be the main constituent of the Earth's core, and considerable efforts have therefore been made to understand its properties at high pressure and temperature. While these efforts have expanded our knowledge of the iron phase diagram, there remain some significant inconsistencies, the most notable being the difference between the 'low' and 'high' melting curves. Here we report the results of molecular dynamics simulations of iron based on embedded atom models fitted to the results of two implementations of density functional theory. We tested two model approximations and found that both point to the stability of the body-centred-cubic (b.c.c.) iron phase at high temperature and pressure. Our calculated melting curve is in agreement with the 'high' melting curve, but our calculated phase boundary between the hexagonal close packed (h.c.p.) and b.c.c. iron phases is in good agreement with the 'low' melting curve. We suggest that the h.c.p.-b.c.c. transition was previously misinterpreted as a melting transition, similar to the case of xenon, and that the b.c.c. phase of iron is the stable phase in the Earth's inner core.     

地球内核物态热效应研究

 在考虑Grüneisen系数、Debye温度随体积变化的基础上,详细讨论了各部分热压的贡献,并且结合K-primed冷压方程,得到了地球内核较为全面的物态方程,分析了内核p-V-T的关系,给出了地球内核的温度分布.其结果详实有据,为以后深入研究地球内核甚至液态外核的物态性质提供了坚实的基础.     

Melting of the Earth's inner core

The Earth's magnetic field is generated by a dynamo in the liquid iron core, which convects in response to cooling of the overlying rocky mantle. The core freezes from the innermost surface outward, growing the solid inner core and releasing light elements that drive compositional convection. Mantle convection extracts heat from the core at a rate that has enormous lateral variations. Here we use geodynamo simulations to show that these variations are transferred to the inner-core boundary and can be large enough to cause heat to flow into the inner core. If this were to occur in the Earth, it would cause localized melting. Melting releases heavy liquid that could form the variable-composition layer suggested by an anomaly in seismic velocity in the 150 kilometres immediately above the inner-core boundary. This provides a very simple explanation of the existence of this layer, which otherwise requires additional assumptions such as locking of the inner core to the mantle, translation from its geopotential centre or convection with temperature equal to the solidus but with composition varying from the outer to the inner core. The predominantly narrow downwellings associated with freezing and broad upwellings associated with melting mean that the area of melting could be quite large despite the average dominance of freezing necessary to keep the dynamo going. Localized melting and freezing also provides a strong mechanism for creating seismic anomalies in the inner core itself, much stronger than the effects of variations in heat flow so far considered.     

Hemispherical variations in seismic velocity at the top of the Earth's inner core

Knowledge of the seismic velocity structure at the top of the Earth's inner core is important for deciphering the physical processes responsible for inner-core growth. Previous global seismic studies have focused on structures found 100 km or deeper within the inner core, with results for the uppermost 100 km available for only isolated regions. Here we present constraints on seismic velocity variations just beneath the inner-core boundary, determined from the difference in travel time between waves reflected at the inner-core boundary and those transmitted through the inner core. We found that these travel-time residuals-observed on both global seismograph stations and several regional seismic networks-are systematically larger, by about 0.8 s, for waves that sample the 'eastern hemisphere' of the inner core (40 degrees E to 180 degrees E) compared to those that sample the 'western hemisphere' (180 degrees W to 40 degrees E). These residuals show no correlation with the angle at which the waves traverse the inner core; this indicates that seismic anisotropy is not strong in this region and that the isotropic seismic velocity of the eastern hemisphere is about 0.8% higher than that of the western hemisphere.     

Iron-silica interaction at extreme conditions and the electrically conducting layer at the base of Earth's mantle

The boundary between the Earth's metallic core and its silicate mantle is characterized by strong lateral heterogeneity and sharp changes in density, seismic wave velocities, electrical conductivity and chemical composition. To investigate the composition and properties of the lowermost mantle, an understanding of the chemical reactions that take place between liquid iron and the complex Mg-Fe-Si-Al-oxides of the Earth's lower mantle is first required. Here we present a study of the interaction between iron and silica (SiO2) in electrically and laser-heated diamond anvil cells. In a multianvil apparatus at pressures up to 140 GPa and temperatures over 3,800 K we simulate conditions down to the core-mantle boundary. At high temperature and pressures below 40 GPa, iron and silica react to form iron oxide and an iron-silicon alloy, with up to 5 wt% silicon. At pressures of 85-140 GPa, however, iron and SiO2 do not react and iron-silicon alloys dissociate into almost pure iron and a CsCl-structured (B2) FeSi compound. Our experiments suggest that a metallic silicon-rich B2 phase, produced at the core-mantle boundary (owing to reactions between iron and silicate), could accumulate at the boundary between the mantle and core and explain the anomalously high electrical conductivity of this region.     

Experimental study on electrical conductivity of dunite at high temperature and pressure

The electrical conductivities of the dunite from the Qinghai-Xizang (Tibetan) Plateau were measured with the impedance spectra method at 1.0-4.0 GPa and 643-1093 K. The experimental results indicated that activation enthalpies of the dunite are smaller than 0.9 eV, the conduction mechanism in dunite may be attributed to the mixed electrical conduction involving grain interiors and boundaries. On the basis of the results of this experiment, we can deduce that there exists cold mantle in the area of Gaize-Lugu in the Qinghai-Xizang (Tibetan) Plateau by reverse methods from the magnetotelluric sounding data (conductivity-depths profile) available for western Tibet. The result provides the present cold mantle viewpoint with strong proof on the basis of high temperature and pressure experiments.     

Density of hydrous silicate melt at the conditions of Earth's deep upper mantle

The chemical evolution of the Earth and the terrestrial planets is largely controlled by the density of silicate melts. If melt density is higher than that of the surrounding solid, incompatible elements dissolved in the melt will be sequestered in the deep mantle. Previous studies on dry (water-free) melts showed that the density of silicate melts can be higher than that of surrounding solids under deep mantle conditions. However, melts formed under deep mantle conditions are also likely to contain some water, which will reduce the melt density. Here we present data constraining the density of hydrous silicate melt at the conditions of approximately 410 km depth. We show that the water in the silicate melt is more compressible than the other components, and therefore the effect of water in reducing melt density is markedly diminished under high-pressure conditions. Our study indicates that there is a range of conditions under which a (hydrous) melt could be trapped at the 410-km boundary and hence incompatible elements could be sequestered in the deep mantle, although these conditions are sensitive to melt composition as well as the composition of the surrounding mantle.     

Dislocation creep in MgSiO3 perovskite at conditions of the Earth's uppermost lower mantle

Seismic anisotropy provides an important observational constraint on flow in the Earth's deep interior. The quantitative interpretation of anisotropy, however, requires knowledge of the slip geometry of the constitutive minerals that are responsible for producing rock fabrics. The Earth's lower mantle is mostly composed of (Mg, Fe)SiO3 perovskite, but as MgSiO3 perovskite is not stable at high temperature under ambient pressure, it has not been possible to investigate its mechanical behaviour with conventional laboratory deformation experiments. To overcome this limitation, several attempts were made to infer the mechanical properties of MgSiO3 perovskite on the basis of analogue materials. But perovskites do not constitute an analogue series for plastic deformation, and therefore the direct investigation of MgSiO3 perovskite is necessary. Here we have taken advantage of recent advances in experimental high-pressure rheology to perform deformation experiments on coarse-grained MgSiO3 polycrystals under pressure and temperature conditions of the uppermost lower mantle. We show that X-ray peak broadening measurements developed in metallurgy can be adapted to low-symmetry minerals to identify the elementary deformation mechanisms activated under these conditions. We conclude that, under uppermost lower-mantle conditions, MgSiO3 perovskite deforms by dislocation creep and may therefore contribute to producing seismic anisotropy in rocks at such depths.     

水化早期水泥净浆电导的特征及与凝结时间相关性的研究

通过测定硅酸盐类水泥净浆在水化早期的电导,获知各样品净浆的电导变化均有共同的特点,笔者将水泥净浆电导在水化早期的变化划分为五个阶段,各阶段均能用水化机理或仪器特点进行解释.实验结果表明,水泥净浆电导在水化早期的极值点LK对应的时间TLK与净浆的初、终凝时间存在良好的相关性.     

Texturing of the Earth's inner core by Maxwell stresses.

Elastic anisotropy in the Earth's inner core has been attributed to a preferred lattice orientation, which may be acquired during solidification of the inner core or developed subsequent to solidification as a result of plastic deformation. But solidification texturing alone cannot explain the observed depth dependence of anisotropy, and previous suggestions for possible deformation processes have all relied on radial flow, which is inhibited by thermal and chemical stratification. Here we investigate the development of anisotropy as the inner core deforms plastically under the influence of electromagnetic (Maxwell) shear stresses. We estimate the flow caused by a representative magnetic field using polycrystal plasticity simulations for epsilon-iron, where the imposed deformation is accommodated by basal and prismatic slip. We find that individual grains in an initially random polycrystal become preferentially oriented with their c axes parallel to the equatorial plane. This pattern is accentuated if deformation is accompanied by recrystallization. Using the single-crystal elastic properties of epsilon-iron at core pressure and temperature, we average over the simulated orientation distribution to obtain a pattern of elastic anisotropy which is similar to that observed seismologically.     

外地核物质的状态参数方程

 在整个外地核空间内流体的对流存在狭窄的上升流通道(浮力团块上升至顶部)和宽阔的下降流通道(周围液体缓慢下降),假定这种对流循环相当快,以至于外地核接近于一个混合良好的等熵态,而各种物理量的偏离可以看作等熵态下的微扰.根据这一对流模型给出了外地核内下降部分物质的组分、温度、压力等状态参数的空间变化方程,以及与之密切相关的径向流体速度方程,并定性分析了各参量对外地核下降流体分层稳定性的贡献.     

试论交流铁芯线圈的铁损

论述了引起交流铁芯线圈在工作中铁损的原因，通过周密分析研究推导了计算铁 定量公式，并纠正了《日本实用电工讲座》译著中的错误结论。     

Intense equatorial flux spots on the surface of the Earth's core

A large number of high-accuracy vector measurements of the Earth's magnetic field have recently become available from the satellite Oersted, complementing previous vector data from the satellite Magsat, which operated in 1979/80. These data can be used to infer the morphology of the magnetic field at the surface of the fluid core, approximately 2,900 km below the Earth's surface. Here I apply a new methodology to these data to calculate maps of the magnetic field at the core surface which show intense flux spots in equatorial regions. The intensity of these features is unusually large--some have intensities comparable to high-latitude flux patches near the poles, previously identified as the major component of the dynamo field. The tendency for pairing of some of these spots to the north and south of the geographical equator suggests they might be associated with the tops of equatorially symmetric columnar structures in the fluid, or their antisymmetric equivalents. The drift of the equatorial features may represent material flow or could represent wave motion; discrimination of these two effects based on future data could provide new information on the strength of the hidden toroidal magnetic field of the Earth.