首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 458 毫秒
1.
  相似文献   

2.
We perform first principles molecular dynamics simulations of Mg2SiO4 liquid and crystalline forsterite. On compression by a factor of two, we find that the Grüneisen parameter of the liquid increases linearly from 0.6 to 1.2. Comparison of liquid and forsterite equations of state reveals a temperature-dependent density crossover at pressures of ∼12-17 GPa. Along the melting curve, which we calculate by integration of the Clapeyron equation, the density crossover occurs within the forsterite stability field at P = 13 GPa and T = 2550 K. The melting curve obtained from the root mean-square atomic displacement in forsterite using the Lindemann law fails to match experimental or calculated melting curves. We attribute this failure to the liquid structure that differs significantly from that of forsterite, and which changes markedly upon compression, with increases in the degree of polymerization and coordination. The mean Si coordination increases from 4 in the uncompressed system to 6 upon twofold compression. The self-diffusion coefficients increase with temperature and decrease monotonically with pressure, and are well described by the Arrhenian relation. We compare our equation of state to the available highpressure shock wave data for forsterite and wadsleyite. Our theoretical liquid Hugoniot is consistent with partial melting along the forsterite Hugoniot at pressures 150-170 GPa, and complete melting at 170 GPa. The wadsleyite Hugoniot is likely sub-liquidus at the highest experimental pressure to date (200 GPa).  相似文献   

3.
We investigate two key transport properties, self-diffusion and viscosity, of Mg2SiO4 liquid as a function of temperature and pressure using density functional theory-based molecular dynamics method. Liquid dynamics in a 224-atom supercell was captured in equilibrium simulations of relatively long durations (50-300 ps) to obtain an acceptable convergence. Our results show that Mg and Si are, respectively, the most and least mobile species at most conditions studied and all diffusivities become similar at high pressure. With increasing temperature from 2200 to 6000 K at ambient pressure, the self-diffusivities increase by factors of 25 (Mg), 80 (Si) and 65 (O), and the viscosity decreases by a factor of 30. The predicted temperature variations of all transport coefficients closely follow the Arrhenian law. However, their pressure variations show a significant non-Arrhenian behavior and also are sensitive to temperature. At 3000 K, the diffusivity (viscosity) decreases (increases) by more than one order of magnitude between 0 and 50 GPa with their activation volumes increasing on compression. Over the entire mantle pressure range, the variations at 4000 K are of two orders of magnitude with nearly constant activation volumes whereas the variations at 6000 K are within one order of magnitude with decreasing activation volumes. The predicted complex dynamical behavior of Mg2SiO4 liquid can be associated with the structural changes occurring on compression. We also estimate the diffusivity and viscosity profiles along a magma ocean isentrope, which suggest that the melt transport properties vary modestly over the relevant magma ocean depth ranges.  相似文献   

4.
A thermodynamic model and equation of state (EOS) is developed from the molecular dynamics simulation experiments of Spera et al. (2009) for CaAl2Si2O8 liquid over the temperature range 3500-6000 K and pressure interval 0-125 GPa. The model is constructed utilizing the isothermal Universal EOS of Vinet et al. (1986) combined with an expression for the temperature-dependence of the internal energy derived from density functional theory (Rosenfeld and Tarazona, 1998). It is demonstrated that this model is more successful at reproducing the data than the temperature-dependent Universal EOS (Vinet et al., 1987) or the volume-explicit EOS of Ghiorso (2004a). Distinct parameterizations are required to model low (<20 GPa) and high (>20 GPa) pressure regimes. This result is ascribed to the affect of liquid structure on macroscopic thermodynamic properties, specifically the interdependence of average cation-oxygen coordination number on the bulk modulus. The thermodynamic transition between the high- and low-pressure parameterizations is modeled as second order, although the nature of the transition is open to question and may well be first order or lambda-like in character.Analysis of the thermodynamic model reveals a predicted region of liquid-liquid un-mixing at low-temperatures (<1624 K) and pressures (<1.257 GPa). These pressure-temperature conditions are above the glass transition temperature but within the metastable liquid region. They represent the highest temperatures yet suggested for liquid-liquid un-mixing in a silicate bulk composition. A shock wave Hugoniot curve is calculated for comparison with the experimental data of Rigden et al. (1989) and of Asimow and Ahrens (2008). The comparison suggests that the model developed in this paper underestimates the density of the liquid by roughly 10% at pressures greater than ∼20 GPa.  相似文献   

5.
Liquid MgSiO3 is a model for the Earth’s magma ocean and of remnant melt present near the core-mantle boundary. Here, models for molten MgSiO3 are computed employing empirical potential molecular dynamics (EPMD) and results are compared to published results including two EPMD studies and three first-principles molecular dynamics (FPMD) models and to laboratory data. The EPMD results derived from the Oganov (OG) potential come closest to the density of MgSiO3 liquid at the 1-bar melting point inferred from the melting curve. At higher P, EPMD densities calculated from the OG potential and FPMD broadly match shock wave studies, with the OG potential yielding the better comparison. Matsui (M) potential results deviate from other studies above ∼50 GPa. Overall, results based on the OG potential compare best to experimental densities over the P-T range of the mantle. Isothermally, upon increasing P the mean coordination numbers () of oxygen around Si and Mg monotonically increase with pressure. Tetrahedral Si and octahedral Si monotonically increase and decrease, respectively, whereas pentahedral Si maximizes at 10-20 GPa. Tetrahedral Mg decreases monotonically as P increases whereas pentahedral, octahedral and higher coordination polyhedra each show similar behavior first increasing and then decreasing after attaining a maximum; the P of the maximum for each polyhedra type migrates to higher P as the CN increases. Free oxygen and oxygen with one nearest neighbor of either Si or Mg decreases whereas Si or Mg with two or three nearest oxygens (i.e., tricluster oxygen) increases with increasing P isothermally. The increase of tricluster oxygen is consistent with spectroscopy on MgSiO3 glass quenched from 2000 K and 0-40 GPa and high-energy X-ray studies constraining the coordination of O around Mg and around Si at 2300 K and 1 bar. Coordination statistics from FPMD studies for O around Si and Si around O are in agreement with the EPMD results based on the M and OG potentials. Mg self-diffusivity is greater than O and Si self-diffusivities for both the M and OG potentials. All D values monotonically decrease with increasing pressure isothermally and all atoms are more diffusive in the M liquid compared to the OG liquid except at T > ∼5000 K and P > 100 GPa. Previously published EPMD diffusivities fall between values given by the M and OG potentials, at least up to 45 GPa. The M liquid is generally less viscous than the OG liquid except at P > ∼80 GPa. Activation energy and volume are around 96 kJ/mol and 1.5 cm3/mol, respectively. The FPMD viscosity results at 120 GPa and 4000 and 4500 K are essentially identical to the values from the M and OG potentials. FPMD viscosity results are similar to the OG results for P < 60 GPa; at higher P, the FPMD viscosities are higher. At 4000 K and 100 GPa the shear viscosity of liquid MgSiO3 is ∼0.1 Pa s. More extensive laboratory results are required to better define the thermodynamic, transport and structural properties of MgSiO3 liquids and for comparison with computational studies.  相似文献   

6.
Armenite, ideal formula BaCa2Al6Si9O30·2H2O, and its dehydrated analog BaCa2Al6Si9O30 and epididymite, ideal formula Na2Be2Si6O15·H2O, and its dehydrated analog Na2Be2Si6O15 were studied by low-temperature relaxation calorimetry between 5 and 300 K to determine the heat capacity, Cp, behavior of their confined H2O. Differential thermal analysis and thermogravimetry measurements, FTIR spectroscopy, electron microprobe analysis and powder Rietveld refinements were undertaken to characterize the phases and the local environment around the H2O molecule.The determined structural formula for armenite is Ba0.88(0.01)Ca1.99(0.02)Na0.04(0.01)Al5.89(0.03)Si9.12(0.02)O30·2H2O and for epididymite Na1.88(0.03)K0.05(0.004)Na0.01(0.004)Be2.02(0.008)Si6.00(0.01)O15·H2O. The infrared (IR) spectra give information on the nature of the H2O molecules in the natural phases via their H2O stretching and bending vibrations, which in the case of epididymite only could be assigned. The powder X-ray diffraction data show that armenite and its dehydrated analog have similar structures, whereas in the case of epididymite there are structural differences between the natural and dehydrated phases. This is also reflected in the lattice IR mode behavior, as observed for the natural phases and the H2O-free phases. The standard entropy at 298 K for armenite is S° = 795.7 ± 6.2 J/mol K and its dehydrated analog is S° = 737.0 ± 6.2 J/mol K. For epididymite S° = 425.7 ± 4.1 J/mol K was obtained and its dehydrated analog has S° = 372.5 ± 5.0 J/mol K. The heat capacity and entropy of dehydration at 298 K are Δ = 3.4 J/mol K and ΔSrxn = 319.1 J/mol K and Δ = −14.3 J/mol K and ΔSrxn = 135.7 J/mol K for armenite and epididymite, respectively. The H2O molecules in both phases appear to be ordered. They are held in place via an ion-dipole interaction between the H2O molecule and a Ca cation in the case of armenite and a Na cation in epididymite and through hydrogen-bonding between the H2O molecule and oxygen atoms of the respective silicate frameworks. Of the three different H2O phases ice, liquid water and steam, the Cp behavior of confined H2O in both armenite and epididymite is most similar to that of ice, but there are differences between the two silicates and from the Cp behavior of ice. Hydrogen-bonding behavior and its relation to the entropy of confined H2O at 298 K is analyzed for various microporous silicates.The entropy of confined H2O at 298 K in various silicates increases approximately linearly with increasing average wavenumber of the OH-stretching vibrations. The interpretation is that decreased hydrogen-bonding strength between a H2O molecule and the silicate framework, as well as weak ion-dipole interactions, results in increased entropy of H2O. This results in increased amplitudes of external H2O vibrations, especially translations of the molecule, and they contribute strongly to the entropy of confined H2O at T < 298 K.  相似文献   

7.
We use a lattice vibrational technique to derive thermophysical and thermochemical properties and phase equilibria in the system MgO-SiO2 at pressures and temperatures relevant to Earth’s mantle. The technique is based on an extension of Kieffer’s model to incorporate details of the phonon spectrum, and it includes treatment of intrinsic anharmonicity. We use a least squares inversion technique applied to available experimental data, and show that it results in an accurate representation of thermodynamic properties and sound wave velocities of high-pressure phases in the system MgSiO3. The vibrational method has been validated against laboratory experimental data in the temperature range between 0 and 2500 K and at pressures between 1 bar and 30 GPa. The technique results in a phase diagram consistent with the majority of thermophysical and thermochemical data. It is shown that intrinsic anharmonicity affects significantly slopes and positions of the phase boundaries. Our analysis indicates inconsistencies in a number of data sets of thermophysical properties for stishovite, majorite and ortho-enstatite necessitating new measurements. For akimotoite elasticity data at high-pressure and high-temperature conditions and 1 bar heat capacity measurements are needed. For stishovite elasticity measurements are necessary to reconcile elasticity data with V-P-T measurements. Additionally V-P-T measurements at pressures higher than 10 GPa are needed, which should be reconciled with V-P-T data at lower pressures. Raman and infrared spectroscopic data are necessary for both clino-enstatite and majorite. Additionally structural data are needed to resolve the discrepancy between values for the degree of disorder in majorite. Volume expansion data for ortho-enstatite are needed and effects causing differences in measured volume expansion should be elucidated.  相似文献   

8.
Petrological and geochemical study of volatile bearing phases (fluid inclusions, amphibole, and nominally anhydrous minerals) in a spinel lherzolite xenolith suite from Quaternary lavas at Injibara (Lake Tana region, Ethiopian plateau) shows compelling evidence for metasomatism in the lithospheric mantle in a region of mantle upwelling and continental flood basalts. The xenolith suite consists of deformed (i.e., protogranular to porphyroclastic texture) Cl-rich pargasite lherzolites, metasomatized (LILE and Pb enrichment in clinopyroxene and amphibole) at T ? 1000 °C. Lherzolites contain chlorine-rich H2O-CO2 fluid inclusions, but no melt inclusions. Fluid inclusions are preserved only in orthopyroxene, while in olivine, they underwent extensive interaction with the host mineral. The metasomatic fluid composition is estimated: XCO2 = 0.64, XH2O = 0.33, XNa = 0.006, XMg = 0.006, XCl = 0.018, (salinity = 14-10 NaCl eq. wt.%, aH2O = 0.2, Cl = 4-5 mol.%). Fluid isochores correspond to trapping pressures of 1.4-1.5 GPa or 50-54 km depth (at T = 950 °C). Synchrotron sourced micro-infrared mapping (ELECTRA, Trieste) shows gradients for H2O-distribution in nominally anhydrous minerals, with considerable enrichment at grain boundaries, along intragranular microfractures, and around fluid inclusions. Total water amounts in lherzolites are variable from about 150 up to 400 ppm. Calculated trace-element pattern of metasomatic fluid phases, combined with distribution and amount of H2O in nominally anhydrous minerals, delineate a metasomatic Cl- and LILE-rich fluid phase heterogeneously distributed in the continental lithosphere. Present data suggest that Cl-rich aqueous fluids were important metasomatic agents beneath the Ethiopian plateau, locally forming a source of high water content in the peridotite, which may be easily melted. High Cl, LILE, and Pb in metasomatic fluid phases suggest the contribution of recycled altered oceanic lithosphere component in their source.  相似文献   

9.
Based on our previous development of the molecular interaction potential for pure H2O and CO2 [Zhang, Z.G., Duan, Z.H. 2005a. Isothermal-isobaric molecular dynamics simulations of the PVT properties of water over wide range of temperatures and pressures. Phys. Earth Planet Interiors149, 335-354; Zhang, Z.G., Duan, Z.H. 2005b. An optimized molecular potential for carbon dioxide. J. Chem. Phys.122, 214507] and the ab initio potential surface across CO2-H2O molecules constructed in this study, we carried out more than one thousand molecular dynamics simulations of the PVTx properties of the CO2-H2O mixtures in the temperature-pressure range from 673.15 to 2573.15 K up to 10.0 GPa. Comparison with extensive experimental PVTx data indicates that the simulated results generally agree with experimental data within 2% in density, equivalent to experimental uncertainty. Even the data under the highest experimental temperature-pressure conditions (up to 1673 K and 1.94 GPa) are well predicted with the agreement within 1.0% in density, indicating that the high accuracy of the simulation is well retained as the temperature and pressure increase. The consistent and stable predictability of the simulation from low to high temperature-pressure and the fact that the molecular dynamics simulation resort to no experimental data but to ab initio molecular potential makes us convinced that the simulation results should be reliable up to at least 2573 K and 10 GPa with errors less than 2% in density. In order to integrate all the simulation results of this study and previous studies [Zhang and Duan, 2005a, 2005b] and the experimental data for the calculation of volumetric properties (volume, density, and excess volume), heat properties, and chemical properties (fugacity, activity, and possibly supercritical phase separation), an equation of state (EOS) is laboriously developed for the CO2, H2O, and CO2-H2O systems. This EOS reproduces all the experimental and simulated data covering a wide temperature and pressure range from 673.15 to 2573.15 K and from 0 to 10.0 GPa within experimental or simulation uncertainty.  相似文献   

10.
The solubility and stability of synthetic grossular were determined at 800 °C and 10 kbar in NaCl-H2O solutions over a large range of salinity. The measurements were made by evaluating the weight losses of grossular, corundum, and wollastonite crystals equilibrated with fluid for up to one week in Pt capsules and a piston-cylinder apparatus. Grossular dissolves congruently over the entire salinity range and displays a large solubility increase of 0.0053 to 0.132 molal Ca3Al2Si3O12 with increasing NaCl mole fraction (XNaCl) from 0 to 0.4. There is thus a solubility enhancement 25 times the pure H2O value over the investigated range, indicating strong solute interaction with NaCl. The Ca3Al2Si3O12 mole fraction versus NaCl mole fraction curve has a broad plateau between XNaCl = 0.2 and 0.4, indicating that the solute products are hydrous; the enhancement effect of NaCl interaction is eventually overtaken by the destabilizing effect of lowering H2O activity. In this respect, the solubility behavior of grossular in NaCl solutions is similar to that of corundum and wollastonite. There is a substantial field of stability of grossular at 800 °C and 10 kbar in the system CaSiO3-Al2O3-H2O-NaCl. At high Al2O3/CaSiO3 bulk compositions the grossular + fluid field is limited by the appearance of corundum. Zoisite appears metastably with corundum in initially pure H2O, but disappears once grossular is nucleated. At XNaCl = 0.3, however, zoisite is stable with corundum and fluid; this is the only departure from the quaternary system encountered in this study. Corundum solubility is very high in solutions containing both NaCl and CaSiO3: Al2O3 molality increases from 0.0013 in initially pure H2O to near 0.15 at XNaCl = 0.4 in CaSiO3-saturated solutions, a >100-fold enhancement. In contrast, addition of Al2O3 to wollastonite-saturated NaCl solutions increases CaSiO3 molality by only 12%. This suggests that at high pH (quench pH is 11-12), the stability of solute Ca chloride and Na-Al ± Si complexes account for high Al2O3 solubility, and that Ca-Al ± Si complexes are minor. The high solubility and basic dissolution reaction of grossular suggest that Al may be a very mobile component in calcareous rocks in the deep crust and upper mantle when migrating saline solutions are present.  相似文献   

11.
High-pressure liquids in the MgO-SiO2-H2O (MSH) system have been investigated at 11 and 13.5 GPa and between 1000 and 1350 °C. A bulk composition more magnesian than the tie-line forsterite-H2O was employed for the study. Rocking multi-anvil experiments were combined with a diamond trap set-up. After termination of the experiments, the liquid trapped in the diamond layer was analysed by laser ablation ICP-MS using the ‘freezing’ technique. At 11 GPa, liquids coexist with one or two of phase A, clinohumite, chondrodite, and forsterite. A marked discontinuity in the evolution of liquid compositions near 1100 °C is observed at 11 GPa. A step of ∼13 wt% H2O and 13 wt% MgO is interpreted to result from overstepping the fluid-saturated solidus reaction mass balanced to 1.00(18) phase A + 1.07(4) fluid = 0.63(15) chondrodite + 1.44(2) melt. At 13.5 GPa liquids coexist with one or two of hydrous wadsleyite, clinohumite, superhydrous B, phase B, and forsterite. The discontinuity in liquid composition is no longer present, indicating that the second critical endpoint of the solidus has been overstepped. Thus, hydrous melts in the Mg-rich part of the MSH system (molar bulk Mg/Si > 2) are chemically distinct from aqueous fluids at pressure up to 11 GPa. Convergence of fluid and melt compositions along the solidus resulting in a supercritical liquid occurs between 11 and 13.5 GPa, at which pressure the entire MSH system becomes supercritical.  相似文献   

12.
This work reports the application of thermodynamic models, including equations of state, to binary (salt-free) CH4-H2O fluid inclusions. A general method is presented to calculate the compositions of CH4-H2O inclusions using the phase volume fractions and dissolution temperatures of CH4 hydrate. To calculate the homogenization pressures and isolines of the CH4-H2O inclusions, an improved activity-fugacity model is developed to predict the vapor-liquid phase equilibrium. The phase equilibrium model can predict methane solubility in the liquid phase and water content in the vapor phase from 273 to 623 K and from 1 to 1000 bar (up to 2000 bar for the liquid phase), within or close to experimental uncertainties. Compared to reliable experimental phase equilibrium data, the average deviation of the water content in the vapor phase and methane solubility in the liquid phase is 4.29% and 3.63%, respectively. In the near-critical region, the predicted composition deviations increase to over 10%. The vapor-liquid phase equilibrium model together with the updated volumetric model of homogenous (single-phase) CH4-H2O fluid mixtures (Mao S., Duan Z., Hu J. and Zhang D. (2010) A model for single-phase PVTx properties of CO2-CH4-C2H6-N2-H2O-NaCl fluid mixtures from 273 to 1273 K and from 1 to 5000 bar. Chem. Geol.275, 148-160), is applied to calculate the isolines, homogenization pressures, homogenization volumes, and isochores at specified homogenization temperatures and compositions. Online calculation is on the website: http://www.geochem-model.org/.  相似文献   

13.
Based on our previous study of the intermolecular potential for pure H2O and the strict evaluation of the competitive potential models for pure CH4 and the ab initio fitting potential surface across CH4-H2O molecules in this study, we carried out more than two thousand molecular dynamics simulations for the PVTx properties of pure CH4 and the CH4-H2O mixtures up to 2573 K and 10 GPa. Comparison of 1941 simulations with experimental PVT data for pure CH4 shows an average deviation of 0.96% and a maximum deviation of 2.82%. The comparison of the results of 519 simulations of the mixtures with the experimental measurements reveals that the PVTx properties of the CH4-H2O mixtures generally agree with the extensive experimental data with an average deviation of 0.83% and 4% in maximum, which is equivalent to the experimental uncertainty. Moreover, the maximum deviation between the experimental data and the simulation results decreases to about 2% as temperature and pressure increase, indicating that the high accuracy of the simulation is well retained in the high temperature and pressure region.After the validation of the simulation method and the intermolecular potential models, we systematically simulated the PVTx properties of this binary system from 673 K and 0.05 GPa to 2573 K and 10 GPa. In order to integrate all the simulation results and the experimental data for the calculation of thermodynamic properties, an equation of state (EOS) is developed for the CH4-H2O system covering 673-2573 K and 0.01-10 GPa. Isochores for compositions <4 mol% CH4 up to 773 K and 600 MPa are also determined in this paper. The program for the EOS can be downloaded from www.geochem-model.org/programs.htm.  相似文献   

14.
Longitudinal acoustic velocities were measured at 1 bar in 10 Na2O-TiO2-SiO2 (NTS) liquids for which previous density and thermal expansion data are reported in the literature. Data were collected with a frequency-sweep acoustic interferometer at centered frequencies of 4.5, 5, and 6 MHz between 1233 and 1896 K; in all cases, the sound speeds decrease with increasing temperature. Six of the liquids have a similar TiO2 concentration (∼25 mol %), so that the effect of varying Na/Si ratio on the partial molar compressibility of the TiO2 component can be evaluated. Theoretically based models for βT and (∂V/∂P)T as a function of composition and temperature are presented. As found previously for the partial molar volume of TiO2 in sodium silicate melts, values of (13.7-18.8 × 10−2/GPa) vary systematically with the Na/Si and Na/(Si + Ti) ratio in the liquid. In contrast values of for the SiO2 and Na2O components (6.6 and 8.0 × 10−2/GPa, respectively, at 1573 K) are independent of composition. Na2O is the only component that contributes to the temperature dependence of the compressibility of NTS liquids (1.13 ± 0.04 × 10−4/GPa K). The results further indicate that the TiO2 component is twice as compressible as the Na2O and SiO2 components. The enhanced compressibility of TiO2 appears to be related to the abundance of five-coordinated Ti ([5]Ti) in these liquids, but not with a change in Ti coordination. Instead, it is proposed that the asymmetric geometry of [5]Ti in a square pyramidal site promotes different topological rearrangements in alkali titanosilicate liquids, which lead to the enhanced compressibility of TiO2.  相似文献   

15.
Superliquidus metal-silicate partitioning was investigated for a number of moderately siderophile (Mo, As, Ge, W, P, Ni, Co), slightly siderophile (Zn, Ga, Mn, V, Cr) and refractory lithophile (Nb, Ta) elements. To provide independent constrains on the effects of temperature, oxygen fugacity and silicate melt composition, isobaric (3 GPa) experiments were conducted in piston cylinder apparatus at temperature between 1600 and 2600 °C, relative oxygen fugacities of IW−1.5 to IW−3.5, and for silicate melt compositions ranging from basalt to peridotite. The effect of pressure was investigated through a combination of piston cylinder and multi-anvil isothermal experiments between 0.5 and 18 GPa at 1900 °C. Oxidation states of siderophile elements in the silicate melt as well as effect of carbon saturation on partitioning are also derived from these results. For some elements (e.g. Ga, Ge, W, V, Zn) the observed temperature dependence does not define trends parallel to those modeled using metal-metal oxide free energy data. We correct partitioning data for solute interactions in the metallic liquid and provide a parameterization utilized in extrapolating these results to the P-T-X conditions proposed by various core formation models. A single-stage core formation model reproduces the mantle abundances of several siderophile elements (Ni, Co, Cr, Mn, Mo, W, Zn) for core-mantle equilibration at pressures from 32 to 42 GPa along the solidus of a deep peridotitic magma ocean (∼3000 K for this pressure range) and oxygen fugacities relevant to the FeO content of the present-day mantle. However, these P-T-fO2 conditions cannot produce the observed concentrations of Ga, Ge, V, Nb, As and P. For more reducing conditions, the P-T solution domain for single stage core formation occurs at subsolidus conditions and still cannot account for the abundances of Ge, Nb and P. Continuous core formation at the base of a magma ocean at P-T conditions constrained by the peridotite liquidus and fixed fO2 yields concentrations matching observed values for Ni, Co, Cr, Zn, Mn and W but underestimates the core/mantle partitioning observed for other elements, notably V, which can be reconciled if accretion began under reducing conditions with progressive oxidation to fO2 conditions consistent with the current concentration of FeO in the mantle as proposed by Wade and Wood (2005). However, neither oxygen fugacity path is capable of accounting for the depletions of Ga and Ge in the Earth’s mantle. To better understand core formation, we need further tests integrating the currently poorly-known effects of light elements and more complex conditions of accretion and differentiation such as giant impacts and incomplete equilibration.  相似文献   

16.
Precise determination of the partitioning of Mg and Fe2+ between olivine and ultramafic melt has been made at pressures from 5 to 13 GPa using a MA-8 type multi-anvil high-pressure apparatus (PREM) installed at Earthquake Research Institute, University of Tokyo. A very short rhenium capsule (<100 μm sample thickness) was adopted to minimize temperature variation within the sample container. Synthetic gels with the composition of the upper mantle peridotite were used as starting materials to promote the homogeneity. Analyses of quenched melts and coexisting olivines were made with an electron probe microanalyzer. The obtained partition coefficient, KD [=(FeO/MgO)ol/(FeO/MgO)melt], decreases from 0.35 to 0.25 with increasing pressure from 5 to 13 GPa, suggesting a negative correlation between pressure and KD above 5 GPa. Our result is consistent with a parabolic relationship between KD and degree of polymerization (NBO/T) of melts reported by previous studies at lower pressures. The negative correlation between pressure and KD suggests that olivine crystallizing in a magma ocean becomes more Mg-rich with depth and that primary magmas generated in the upper mantle become more Fe-rich with depth than previously estimated.  相似文献   

17.
We present new high-pressure temperature experiments on melting phase relations of Fe-C-S systems with applications to metallic core formation in planetary interiors. Experiments were performed on Fe-5 wt% C-5 wt% S and Fe-5 wt% C-15 wt% S at 2-6 GPa and 1050-2000 °C in MgO capsules and on Fe-13 wt% S, Fe-5 wt% S, and Fe-1.4 wt% S at 2 GPa and 1600 °C in graphite capsules. Our experiments show that: (a) At a given P-T, the solubility of carbon in iron-rich metallic melt decreases modestly with increasing sulfur content and at sufficiently high concentration, the interaction between carbon and sulfur can cause formation of two immiscible melts, one rich in Fe-carbide and the other rich in Fe-sulfide. (b) The mutual solubility of carbon and sulfur increases with increasing pressure and no super-liquidus immiscibility in Fe-rich compositions is likely expected at pressures greater than 5-6 GPa even for bulk compositions that are volatile-rich. (c) The liquidus temperature in the Fe-C-S ternary is significantly different compared to the binary liquidus in the Fe-C and Fe-S systems. At 6 GPa, the liquidus of Fe-5 wt% C-5 wt% S is 150-200 °C lower than the Fe-5 wt% S. (d) For Fe-C-S bulk compositions with modest concentration of carbon, the sole liquidus phase is iron carbide, Fe3C at 2 GPa and Fe7C3 at 6 GPa and metallic iron crystallizes only with further cooling as sulfur is concentrated in the late crystallizing liquid. Our results suggest that for carbon and sulfur-rich core compositions, immiscibility induced core stratification can be expected for planets with core pressure less than ∼6 GPa. Thus planetary bodies in the outer solar system such as Ganymede, Europa, and Io with present day core-mantle boundary (CMB) pressures of ∼8, ∼5, and 7 GPa, respectively, if sufficiently volatile-rich, may either have a stratified core or may have experienced core stratification owing to liquid immiscibility at some stage of their accretion. A similar argument can be made for terrestrial planetary bodies such as Mercury and Earth’s Moon, but no such stratification is predicted for cores of terrestrial planets such as Earth, Venus, and Mars with the present day core pressure in the order ?136 GPa, ?100 GPa, and ?23 GPa. (e) Owing to different expected densities of Fe-rich (and carbon-bearing) and sulfur-rich metallic melts, their settling velocities are likely different; thus core formation in terrestrial planets may involve rain of more than one metallic melt through silicate magma ocean. (f) For small planetary bodies that have core pressures <6 GPa and have a molten core or outer core, settling of denser carbide-rich liquid or flotation of lighter, sulfide-rich melt may contribute to an early, short-lived geodynamo.  相似文献   

18.
Molecular diffusion in natural volcanic liquids discriminates between isotopes of major ions (e.g., Fe, Mg, Ca, and Li). Although isotope separation by diffusion is expected on theoretical grounds, the dependence on mass is highly variable for different elements and in different media. Silicate liquid diffusion experiments using simple liquid compositions were carried out to further probe the compositional dependence of diffusive isotopic discrimination and its relationship to liquid structure. Two diffusion couples consisting of the mineral constituents anorthite (CaAl2Si2O8; denoted AN), albite (NaAlSi3O8; denoted AB), and diopside (CaMgSi2O6; denoted DI) were held at 1450 °C for 2 h and then quenched to ambient pressure and temperature. Major-element as well as Ca and Mg isotope profiles were measured on the recovered quenched glasses. In both experiments, Ca diffuses rapidly with respect to Si. In the AB-AN experiment, DCa/DSi ≈ 20 and the efficiency of isotope separation for Ca is much greater than in natural liquid experiments where DCa/DSi ≈ 1. In the AB-DI experiment, DCa/DSi ≈ 6 and the efficiency of isotope separation is between that of the natural liquid experiments and the AB-AN experiment. In the AB-DI experiment, DMg/DSi ≈ 1 and the efficiency of isotope separation for Mg is smaller than it is for Ca yet similar to that observed for Mg in natural liquids.The results from the experiments reported here, in combination with results from natural volcanic liquids, show clearly that the efficiency of diffusive separation of Ca isotopes is systematically related to the solvent-normalized diffusivity - the ratio of the diffusivity of the cation (DCa) to the diffusivity of silicon (DSi). The results on Ca isotopes are consistent with available data on Fe, Li, and Mg isotopes in silicate liquids, when considered in terms of the parameter Dcation/DSi. Cations diffusing in aqueous solutions display a similar relationship between isotopic separation efficiency and Dcation/DH2O, although the efficiencies are smaller than in silicate liquids. Our empirical relationship provides a tool for predicting the magnitude of diffusive isotopic effects in many geologic environments and a basis for a more comprehensive theory of isotope separation in liquid solutions. We present a conceptual model for the relationship between diffusivity and liquid structure that is consistent with available data.  相似文献   

19.
The global rise in atmospheric greenhouse gas concentrations calls for practicable solutions to capture CO2. In this study, a mineral carbonation process was applied in which CO2 reacts with alkaline lignite ash and forms stable carbonate solids. In comparison to previous studies, the assays were conducted at low temperatures and pressures and under semi-dry reaction conditions in an 8 L laboratory mixing device. In order to find optimum process conditions the pCO2 (10-20%), stirring rate (500-3000 rpm) and the liquid to solid ratio (L/S = 0.03-0.36 L kg−1) were varied. In all experiments a considerable CO2 uptake from the gas phase was observed. Concurrently the solid phase contents of Ca and Mg (hydr)oxides decreased and CaCO3 and MgCO3 fractions increased throughout the experiments, showing that CO2 was stabilized as a solid carbonate. The carbonation reaction depends on three factors: Dissolution of CO2 in the liquid phase, mobilization of Ca and Mg from the mineral surface and precipitation of the carbonate solids. Those limitations were found to depend strongly on the variation of the process parameters. Optimum reaction conditions could be found for L/S ratios between 0.12 and 0.18, medium stirring velocities and pCO2 between 10% and 20%.Maximum CO2 uptake by the solid phase was 4.8 mmol g−1 after 120 min, corresponding to a carbonation efficiency for the alkaline material of 53% of the theoretical CO2 binding capacity. In comparison to previous studies both CO2 uptake and carbonation efficiencies were in a similar range, but the reaction times in the semi-dry process were considerably shorter. The proposed method additionally allows for a more simple carbonation setup due to low T and P, and produces an easier to handle product with low water content.  相似文献   

20.
The sodium solubility in silicate melts in the CaO-MgO-SiO2 (CMS) system at 1400 °C has been measured by using a closed thermochemical reactor designed to control alkali metal activity. In this reactor, Na(g) evaporation from a Na2O-xSiO2 melt imposes an alkali metal vapor pressure in equilibrium with the molten silicate samples. Because of equilibrium conditions in the reactor, the activity of sodium-metal oxide in the molten samples is the same as that of the source, i.e., aNa2O(sample) = aNa2O(source). This design also allows to determine the sodium oxide activity coefficient in the samples. Thirty-three different CMS compositions were studied. The results show that the amount of sodium entering from the gas phase (i.e., Na2O solubility) is strongly sensitive to silica content of the melt and, to a lesser extent, the relative amounts of CaO and MgO. Despite the large range of tested melt compositions (0 < CaO and MgO < 40; 40 < SiO2 < 100; in wt%), we found that Na2O solubility is conveniently modeled as a linear function of the optical basicity (Λ) calculated on a Na-free basis melt composition. In our experiments, γNa2O(sample) ranges from 7 × 10−7 to 5 × 10−6, indicating a strongly non-ideal behavior of Na2O solubility in the studied CMS melts (γNa2O(sample) ? 1). In addition to showing the effect of sodium on phase relationships in the CMS system, this Na2O solubility study brings valuable new constraints on how melt structure controls the solubility of Na in the CMS silicate melts. Our results suggest that Na2O addition causes depolymerization of the melt by preferential breaking of Si-O-Si bonds of the most polymerized tetrahedral sites, mainly Q4.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号