Heat capacity of γ-Fe2SiO4 between 5 and 303 K and derived thermodynamic properties

A multi-anvil device was used to synthesize 24 mg of pure γ-Fe₂SiO₄ crystals at 8.5 GPa and 1,273 K. The low-temperature heat capacity (C _p) of γ-Fe₂SiO₄ was measured between 5 and 303 K using the heat capacity option of a physical properties measurement system. The measured heat capacity data show a broad λ-transition at 11.8 K. The difference in the C _p between fayalite and γ-Fe₂SiO₄ is reduced as the temperature increases in the range of 50–300 K. The gap in C _p data between 300 and 350 K of γ-Fe₂SiO₄ is an impediment to calculation of a precise C _p equation above 298 K that can be used for phase equilibrium calculations at high temperatures and high pressures. The C _p and entropy of γ-Fe₂SiO₄ at standard temperature and pressure (S°₂₉₈) are 131.1 ± 0.6 and 140.2 ± 0.4 J mol⁻¹ K⁻¹, respectively. The Gibbs free energy at standard pressure and temperature (ΔG° _f,298) is calculated to be −1,369.3 ± 2.7 J mol⁻¹ based on the new entropy data. The phase boundary for the fayalite–γ-Fe₂SiO₄ transition at 298 K based on current thermodynamic data is located at 2.4 ± 0.6 GPa with a slope of 25.4 bars/K, consistent with extrapolated results of previous experimental studies.     

A computer simulation study of (OH) defects in olivine

Recent experiments (Miller et al. 1987; Bell and Rossman 1992; Bai and Kohlstedt 1992) have shown that olivine, the dominant mineral in the Earth's upper mantle, can contain substantial amounts of water in the form of OH. There is uncertainty, however, as to the mechanisms by which water dissolves into the mineral structure and as to the site it occupies in the host lattice. We have therefore used atomistic computer simulation techniques based on the Born model of solids, to investigate the structures and energies of OH defects in olivine. Our calculations suggest that the most favourable route for incorporation of OH into olivine involves reactions with water accompanied by the reduction of ferric iron for which we obtain a solution energy of 0.46 eV. We propose therefore that the OH content will be largely controlled by the concentration of Fe³⁺.     

Experimental Investigation and Optimization of Thermodynamic Properties and Phase Diagrams in the Systems CaO-SiO2, MgO-SiO2, CaMgSi2O6-SiO2 and CaMgSi2O6-Mg2SiO4 to 1{middle dot}0 GPa

Using experimental results at 1·0 GPa for the systemsCaO–SiO₂, MgO–SiO₂, CaMgSi₂O₆–SiO₂ and CaMgSi₂O₆–Mg₂SiO₄,and all the currently available phase equilibria and thermodynamicdata at 1 bar, we have optimized the thermodynamic propertiesof the liquid phase at 1·0 GPa. The new optimized thermodynamicparameters indicate that pressure has little effect on the topologyof the CaO–SiO₂, CaMgSi₂O₆–SiO₂, and CaMgSi₂O₆–Mg₂SiO₄systems but a pronounced one on the MgO–SiO₂ binary. Themost striking change concerns passage of the MgSiO₃ phase fromperitectic melting at 1 bar to eutectic melting at 1·0GPa. This transition is estimated to occur at 0·41 GPa.For the CaMgSi₂O₆–SiO₂ and CaMgSi₂O₆–Mg₂SiO₄ pseudo-binaries,the size of the field clinopyroxene + liquid increases withincreasing pressure. This change is related to the shift ofthe piercing points clinopyroxene + silica + liquid (from 0·375mol fraction SiO₂ at 1 bar to 0·414 at 1·0 GPa)and clinopyroxene + olivine + liquid (from 0·191 molfraction SiO₂ at 1 bar to 0·331 at 1·0 GPa) thatbound the clinopyroxene + liquid field in the CaMgSi₂O₆·SiO₂and CaMgSi₂O₆·Mg₂SiO₄ pseudo-binaries, respectively. KEY WORDS: CaO–SiO₂; CaMgSi₂O₆–Mg₂SiO₄; CaMgSi₂O₆–SiO₂; experiments; MgO–SiO₂     

高温高压下Mg2SiO4-MgAlO4体系相变的实验研究

采用多顶砧高压实验装置研究了Mg2SiO4-MgAl2O4体系在压力为22 GPa,温度为1550～1750℃条件下的相变,并考查了Al2O3在γ相中的固溶度.结果表明,随着体系中MgAl2O4组分含量的增加,相组合发生了变化,依次为γ相+镁铝硅酸盐固溶体+方镁石→镁铝硅酸盐固溶体+方镁石→镁铝硅酸盐固溶体+方镁石+刚玉固溶体;镁铝硅酸盐固溶体具有石榴子石结构,其化学成分随着体系中共存相的改变而有所变化;Al2O3在γ相中的固溶度很低(其重量百分比＜0.8%),因此,在Mg2SiO4-MgAl2O4体系中Al2O3可能对γ相超尖晶石分解转变的压力不会有很大的影响.     

The phase boundary between wadsleyite and ringwoodite in Mg2SiO4 determined by in situ X-ray diffraction

The phase boundary between wadsleyite and ringwoodite in Mg₂SiO₄ has been determined in situ using a multi-anvil apparatus and synchrotron X-rays radiation at SPring-8. In spite of the similar X-ray diffraction profiles of these high-pressure phases with closely related structures, we were able to identify the occurrence of the mutual phase transformations based on the change in the difference profile by utilizing a newly introduced press-oscillation system. The boundary was located at ~18.9 GPa and 1,400°C when we used Shim’s gold pressure scale (Shim et al. in Earth Planet Sci Lett 203:729–739, 2002), which was slightly (~0.8 GPa) lower than the pressure as determined from the quench experiments of Katsura and Ito (J Geophys Res 94:15663–15670, 1989). Although it was difficult to constrain the Clapeyron slope based solely on the present data due to the kinetic problem, the phase boundary [P (GPa)=13.1+4.11×10⁻³×T (K)] calculated by a combination of a P–T position well constrained by the present experiment and the calorimetric data of Akaogi et al. (J Geophys Res 94:15671–15685, 1989) reasonably explains all the present data within the experimental error. When we used Anderson’s gold pressure scale (Anderson et al. in J Appl Phys 65:1535–1543, 1989), our phase boundary was located in ~18.1 GPa and 1,400°C, and the extrapolation boundary was consistent with that of Kuroda et al. (Phys Chem Miner 27:523–532, 2000), which was determined at high temperature (1,800–2,000°C) using a calibration based on the same pressure scale. Our new phase boundary is marginally consistent with that of Suzuki et al. (Geophys Res Lett 27:803–806, 2000) based on in situ X-ray experiments at lower temperatures (<1,000°C) using Brown’s and Decker’s NaCl pressure scales.     

High-T phase transition of synthetic ANaB(LiMg)CMg5Si8O22(OH)2 amphibole: an X-ray synchrotron powder diffraction and FTIR spectroscopic study

The synthetic amphibole Na_0.95(Li_0.95Mg_1.05)Mg₅Si₈O₂₂(OH)₂ was studied in situ at high-T, using IR OH-stretching spectroscopy and synchrotron X-ray powder diffraction. At room-T the sample has P2₁ /m symmetry, as shown by the FTIR spectrum. It shows in the OH region two well-defined and intense absorptions at 3,748 and 3,712 cm⁻¹, respectively, and two minor bands at 3,667 and 3,687 cm⁻¹. The main bands are assigned to the two independent O–H groups in the primitive structure. The two minor bands evidencing the presence of small amount of vacant A-site (^A□_0.05). With increasing T, these bands shift continuously and merge into a unique absorption at high temperature. A change as a function of increasing T is revealed by the evolution of the refined unit-cell parameters, whose trend shows a transition to C2/m at about 320–330°C. The spontaneous scalar strain, fitted with a tricritical 2–6 Landau potential, gives a T _c of 325(10)°C (β parameter = 0.27). Comparison with the second-order P2₁ /m ⇔ C2/m phase transition at 255°C for synthetic amphibole ^ANa_0.8^B(Na_0.8Mg_1.2)^CMg₅Si₈O₂₂(OH)₂ indicates that the substitution of Na with Li at the B-sites strongly affects the thermodynamic character and the T _c of the phase transition. The comparison of LNMSH amphiboles with cummingtonitic ones shows that the high-T thermodynamic behaviour is affected by A-site occupancy.     

Low-temperature heat capacities of MgAl<Subscript>2</Subscript>O<Subscript>4</Subscript> and spinels of the MgCr<Subscript>2</Subscript>O<Subscript>4</Subscript>–MgAl<Subscript>2</Subscript>O<Subscript>4</Subscript> solid solution

We present results from low-temperature heat capacity measurements of spinels along the solid solution between MgAl₂O₄ and MgCr₂O₄. The data also include new low-temperature heat capacity measurements for MgAl₂O₄ spinel. Heat capacities were measured between 1.5 and 300 K, and thermochemical functions were derived from the results. No heat capacity anomaly was observed for MgAl₂O₄ spinel; however, we observe a low-temperature heat capacity anomaly for Cr-bearing spinels at temperatures below 15 K. From our data we calculate standard entropies (298.15 K) for Mg(Cr,Al)₂O₄ spinels. We suggest a standard entropy for MgAl₂O₄ of 80.9 ± 0.6 J mol⁻¹ K⁻¹. For the solid solution between MgAl₂O₄ and MgCr₂O₄, we observe a linear increase of the standard entropies from 80.9 J mol⁻¹ K⁻¹ for MgAl₂O₄ to 118.3 J mol⁻¹ K⁻¹ for MgCr₂O₄.     

Low-temperature heat capacities,entropies and enthalpies of Mg<Subscript>2</Subscript>SiO<Subscript>4</Subscript> polymorphs,and α−β−γ and post-spinel phase relations at high pressure

The low-temperature isobaric heat capacities (C _p) of β- and γ-Mg₂SiO₄ were measured at the range of 1.8–304.7 K with a thermal relaxation method using the Physical Property Measurement System. The obtained standard entropies (S°₂₉₈) of β- and γ-Mg₂SiO₄ are 86.4 ± 0.4 and 82.7 ± 0.5 J/mol K, respectively. Enthalpies of transitions among α-, β- and γ-Mg₂SiO₄ were measured by high-temperature drop-solution calorimetry with gas-bubbling technique. The enthalpies of the α−β and β−γ transitions at 298 K (ΔH°₂₉₈) in Mg₂SiO₄ are 27.2 ± 3.6 and 12.9 ± 3.3 kJ/mol, respectively. Calculated α−β and β−γ transition boundaries were generally consistent with those determined by high-pressure experiments within the errors. Combining the measured ΔH°₂₉₈ and ΔS°₂₉₈ with selected data of in situ X-ray diffraction experiments at high pressure, the ΔH°₂₉₈ and ΔS°₂₉₈ of the α−β and β−γ transitions were optimized. Calculation using the optimized data tightly constrained the α−β and β−γ transition boundaries in the P, T space. The slope of α−β transition boundary is 3.1 MPa/K at 13.4 GPa and 1,400 K, and that of β−γ boundary 5.2 MPa/K at 18.7 GPa and 1,600 K. The post-spinel transition boundary of γ-Mg₂SiO₄ to MgSiO₃ perovskite plus MgO was also calculated, using the optimized data on γ-Mg₂SiO₄ and available enthalpy and entropy data on MgSiO₃ perovskite and MgO. The calculated post-spinel boundary with a Clapeyron slope of −2.6 ± 0.2 MPa/K is located at pressure consistent with the 660 km discontinuity, considering the error of the thermodynamic data.     

Breakdown of hydrous ringwoodite to pyroxene and spinelloid at high P and T and oxidizing conditions

To get deeper insight into the phase relations in the end-member system Fe₂SiO₄ and in the system (Fe, Mg)₂SiO₄ experiments were performed in a multi-anvil apparatus at 7 and 13 GPa and 1,000–1,200°C as a function of oxygen fugacity. The oxygen fugacity was varied using the solid oxygen buffer systems Fe/FeO, quartz–fayalite–magnetite, MtW and Ni/NiO. The run products were characterized by electron microprobe, Raman- and FTIR-spectroscopy, X-ray powder diffraction and transmission electron microscopy. At fO₂ corresponding to Ni/NiO Fe-ringwoodite transforms to ferrosilite and spinelloid according to the reaction: 9 Fe₂SiO₄ + O₂ = 6 FeSiO₃ + 5 Fe_2.40Si_0.60O₄. Refinement of site occupancies in combination with stoichiometric Fe³⁺ calculations show that 32% of the total Fe is incorporated as Fe³⁺ according to From the Rietveld refinement we identified spl as spinelloid III (isostructural with wadsleyite) and/or spinelloid V. As we used water in excess in the experiments the run products were also analyzed for structural water incorporation. Adding Mg to the system increases the stability field of ringwoodite to higher oxygen fugacity and the spinel structure seems to accept higher Fe³⁺ but also water concentrations that may be linked. At oxygen fugacity corresponding to MtW conditions similar phase relations in respect to the breakdown reaction in the Fe-end-member system were observed but with a strong fractionation of Fe into spl and Mg into coexisting cpx. Thus, through this strong fractionation it is possible to stabilize very Fe-rich wadsleyite with considerable Fe³⁺ concentrations even at an intermediate Fe–Mg bulk composition: assuming constant K _D independent on composition and a bulk composition of x _Fe = 0.44 this fractionation would stabilize spl with x _Fe = 0.72. Thus, spl could be a potential Fe³⁺ bearing phase at P–T conditions of the transition zone but because of the oxidizing conditions and the Fe-rich bulk composition needed one would expect it more in subduction zone environments than in the transition zone in senso stricto.

First-principles simulation of high-pressure polymorphs in MgAl<Subscript>2</Subscript>O<Subscript>4</Subscript>

We have used density functional theory to investigate the stability of MgAl₂O₄ polymorphs under pressure. Our results can reasonably explain the transition sequence of MgAl₂O₄ polymorphs observed in previous experiments. The spinel phase (stable at ambient conditions) dissociates into periclase and corundum at 14 GPa. With increasing pressure, a phase change from the two oxides to a calcium-ferrite phase occurs, and finally transforms to a calcium-titanate phase at 68 GPa. The calcium-titanate phase is stable up to at least 150 GPa, and we did not observe a stability field for a hexagonal phase or periclase + Rh₂O₃(II)-type Al₂O₃. The bulk moduli of the phases calculated in this study are in good agreement with those measured in high-pressure experiments. Our results differ from those of a previous study using similar methods. We attribute this inconsistency to an incomplete optimization of a cell shape and ionic positions at high pressures in the previous calculations.     

A MEG study of the olivine and spinel forms of Mg2SiO4

In this paper we present a theoretical investigation of the structures and relative stability of the olivine and spinel phases of Mg₂SiO₄. We use both a purely ionic model, based on the Modified Electron Gas (MEG) model of intermolecular forces, and a bond polarization model, developed for low pressure silica phases, to investigate the role of covalency in these compounds. The standard MEG ionic model gives adequate structural results for the two phases but incorrectly predicts the spinel phase to be more stable at zero pressure. This is mainly because the ionic modeling of Mg₂SiO₄ only accounts for 95 percent of the lattice energy. The remainder can be attributed to covalency and many-body effects. An extension of the MEG ionic model using “many-body” pair potentials corrects the phase stability error, but predicts structures which are in poorer agreement with experiment than the standard ionic approach. In addition, calculations using these many-body pair potentials can only account for 10 percent of the missing lattice energy. This model predicts an olivine-spinel phase transition of 8 GPa, below the experimental value of 20 GPa. Therefore, in order to understand more fully the stability of these structures we must consider polarization. A two-shell bond polarization model enhances the stability of both structures, with the olivine structure being stabilized more. This model predicts a phase transition at about 80 GPa, well above the observed value. Also, the olivine and spinel structures calculated with this approach are in poorer agreement with experiment than the ionic model. Therefore, based on our investigations, to properly model covalency in Mg₂SiO₄, a treatment more sophisticated than the two-shell model is needed.     

Reaction between liquid iron and (Mg,Fe)SiO3-perovskite and solubilities of Si and O in molten iron at 27 GPa

The Earth’s core contains light elements and their identification is essential for our understanding of the thermal structure and convection in the core that drives the geodynamo and heat flow from the core to the mantle. Solubilities of Si and O in liquid iron coexisting with (Mg,Fe)SiO₃-perovskite, a major constituent of the lower mantle, were investigated at temperatures between 2,320 and 3,040 K at 27 GPa. It was observed that Si dissolved in the liquid iron up to 1.70 wt% at 3,040 K and O dissolved in the liquid iron up to 7.5 wt% at 2,800 K. It was also clearly seen that liquid iron reacts with (Mg,Fe)SiO₃-perovskite to form magnesiowüstite and it contains Si and O at 27 GPa and at 2,640 and 3,040 K. The amounts of Si and O in the liquid iron are 1.70 and 2.25 wt% at 3,040 K, respectively. The solubilities of Si and O in liquid iron coexisting with (Mg,Fe)SiO₃-perovskite have strong positive temperature dependency. Hence, they can be plausible candidates for the light elements in the core.     

The stability and compressibility of MgAl2O4 high-pressure polymorphs

High-pressure and high-temperature experiments using a laser-heated diamond anvil cell (LHDAC) and synchrotron X-ray diffraction have revealed a phase transition in MgAl₂O₄. CaTi₂O₄-type MgAl₂O₄ was found to be stable at pressures between 45 and at least 117 GPa. The transition pressure of CaTi₂O₄-type phase in MgAl₂O₄ is much lower than that in the natural N-type mid-oceanic ridge basalt composition. The Birch–Murnaghan equation of state for CaTi₂O₄-type MgAl₂O₄ was determined from the experimental unit cell parameters with K ₀=219(±6) GPa, K ₀′=4(constrained value), and V ₀=238.9(±9) ų. The observed compressibility was in agreement with the theoretical compressibility calculated in a previous study. ε-MgAl₂O₄ was observed at pressures between 40 and 45 GPa, which has not been reported in natural rock compositions. The gradient (dP/dT slope) of the transition from the ε-type to CaTi₂O₄-type MgAl₂O₄ had a positive value. These results should resolve the dispute regarding the stable high-pressure phase of MgAl₂O₄, which has been reported in earlier studies using both the multi-anvil press and the diamond anvil cell.     

Heat capacity and phase equilibria of wadeite-type K2Si4O9

The low-temperature heat capacity (C _p) of Si-wadeite (K₂Si₄O₉) synthesized with a piston cylinder device was measured over the range of 5–303 K using the heat capacity option of a physical properties measurement system. The entropy of Si-wadeite at standard temperature and pressure calculated from the measured heat capacity data is 253.8 ± 0.6 J mol⁻¹ K⁻¹, which is considerably larger than some of the previous estimated values. The calculated phase transition boundaries in the system K₂O–Al₂O₃–SiO₂ are generally consistent with previous experimental results. Together with our calculated phase boundaries, seven multi-anvil experiments at 1,400 K and 6.0–7.7 GPa suggest that no equilibrium stability field of kalsilite + coesite intervenes between the stability field of sanidine and that of coesite + kyanite + Si-wadeite, in contrast to previous predictions. First-order approximations were undertaken to calculate the phase diagram in the system K₂Si₄O₉ at lower pressure and temperature. Large discrepancies were shown between the calculated diagram compared with previously published versions, suggesting that further experimental or/and calorimetric work is needed to better constrain the low-pressure phase relations of the K₂Si₄O₉ polymorphs. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

直至核幔边界的Mg2SiO4成分熔体剪切粘滞度的平衡分子动力学研究

熔体的粘滞度极大地影响熔体的各处动力学行为。虽然有一些有关酸碱生熔体粘滞度的数据但涉及地球更深部过程的基性和超基性熔体的粘滞度的数据很少。本项研究利用平衡分子动力学方法Ｍｇ２ＳｉＯ４成分熔体从２１６９Ｋ到５５９５Ｋ。压力从１０＾５Ｐａ到１３５ＧＰａ的粘滞度,所获得的粘滞度数据介于Ｕｒｂａｉｎｅｔ ａｌ,（１９８２）和Ｉｖａｎｏｖ ａｎｄ Ｓｔｅｎｇｅｌｍｅｙｅｒ（１９８２）的实验数据之间,比Ｓｈａ     

Variation of the oxygen content and point defects in olivines, (Fe x Mg 1−x )2SiO4, 0.2 ≤ x ≤ 1.0

 The variation of the oxygen content in olivines, (Fe _x Mg_{1− x} )₂SiO₄, with 0.2 ≤ x ≤ 1.0, was investigated by thermogravimetric measurements. Mass changes occurring upon oxygen activity changes were measured as a function of oxygen activity and cationic composition at 1130 and 1200 °C. During the measurements the samples were in direct contact with gases containing CO, CO₂ and N₂ and, at a few spots at the bottom of the sample stack, also with SiO₂. By fitting experimental data of mass changes to equations derived using point defect thermodynamics, it was shown for olivines with 0.2 ≤ x ≤ 1.0 at 1130 °C and 0.2 ≤ x ≤ 0.7 at 1200 °C within the oxygen activity ranges investigated that the observed variations in the oxygen contents are compatible with cation vacancies and Fe³⁺ ions on M sites and Fe³⁺ ions on silicon sites as majority defects if it is assumed that only three types of point defects occur as majority defects. Different cases were considered, closed systems, taking into account that ξ=[Si]/([Si]+[Fe]+[Mg]) is not necessarily equal to 1/3, and olivines in equilibrium with SiO₂ or pyroxenes. The oxygen content variations observed in this study are significantly smaller than those reported previously in the literature. It is proposed that these differences are related to the dissolution of Fe into noble metal containers used as sample holders in earlier studies and/or to the presence of secondary phases. Received: 1 November 1995 / Accepted: 15 September 2002 Acknowledgements This work was supported by the Cornell Center for Materials Research (CCMR), a Materials Research Science and Engineering Center of the National Science Foundation (DMR-0079992). The authors thank Mr. Daniel M. DiPasquo and Mr. Jason A. Schick for helping in experimental work.     

Ab initio study of the high-pressure behavior of CaSiO3 perovskite

Using density functional simulations, within the generalized gradient approximation and projector-augmented wave method, we study structures and energetics of CaSiO₃ perovskite in the pressure range of the Earths lower mantle (0–150 GPa). At zero Kelvin temperature the cubic CaSiO₃ perovskite structure is unstable in the whole pressure range, at low pressures the orthorhombic (Pnam) structure is preferred. At 14.2 GPa there is a phase transition to the tetragonal (I4/mcm) phase. The CaIrO₃-type structure is not stable for CaSiO₃. Our results also rule out the possibility of decomposition into oxides.

Radiation-induced defects in quartz. II. Single-crystal W-band EPR study of a natural citrine quartz

Single-crystal electron paramagnetic resonance (EPR) spectra of a natural citrine quartz without any artificial irradiation, measured at W-band frequencies (∼94 GHz) and temperatures of 77, 110 and 298 K, allow better characterization of three previously-reported Centers (#6, #7 and B) and discovery of three new defects (B′, C′ and G′). The W-band EPR spectra reveal that Centers #6 and #7 do not reside on twofold symmetry axes, contrary to results from a previous X-band EPR study. The W-band spectra also show that the previously reported Center B is a mixture of two defects (B and B′) with similar g matrices but different-sized ²⁷Al hyperfine structures. Center C′ has similar principal g values to the previously reported Center C but is distinct from the latter by a larger ²⁷Al hyperfine structure with splittings from 0.10 to 0.22 mT. Also, Center G′ has a similar g matrix to the previously reported Center G but a different ²⁷Al hyperfine structure with splittings from 0.41 to 0.53 mT. These spin-Hamiltonian parameters, together with observed thermal properties and microwave-power dependence, suggest that Centers #6 and #7 probably represent O₂³⁻ type defects. Centers B and B′ are probably superoxide radicals (O₂⁻) with the unpaired spin localized on the same pair of oxygen atoms around a missing Si atom but linked to a substitutional Al³⁺ ion each at different neighboring tetrahedral sites. Similarly, Centers G and G′ are most likely superoxide radicals with the unpaired spin localized on another pair of oxygen atoms around a missing Si atom and linked to a substitutional Al³⁺ ion each at different neighboring tetrahedral sites. Center C′ is probably an ozonide radical associated with a missing Si atom and linked to a substitutional Al³⁺ ion at the neighboring tetrahedral site. This study exemplifies the value of  high-frequency EPR for discrimination of  similar defect centers and determination of  small local structural distortions that are often difficult to resolve in conventional  X- and Q-band EPR studies.     

溶出珊瑚砂中Ca~(2+)、Mg~(2+)的动力学实验研究及其影响因素

通过改变固液比、摇床转速、珊瑚砂粒径、温度、溶液pH值及溶液含盐量等参数,对珊瑚砂在水溶液中溶出Ca~(2+)、Mg~(2+)进行了实验,以探讨溶出过程中的动力学规律和影响因素。实验结果表明,珊瑚砂中Ca~(2+)、Mg~(2+)溶出量随反应时间逐渐增大;摇床转速越快、固液比越大、温度越高、溶液pH值越低,Ca~(2+)、Mg~(2+)溶出量越大;当珊瑚砂粒径为2.36~4.75 mm、溶液含盐量为100 mg/L时,Ca~(2+)、Mg~(2+)溶出量最大。统计分析表明,摇床转速、温度及溶液p H值均对珊瑚砂溶出有显著影响,但溶液p H值影响最大。珊瑚砂在水中的溶出过程符合收缩核内扩散模型,表明控制整个溶出过程反应速率的决定因素是内扩散速率;在15~40℃时,珊瑚砂在纯水中溶出Ca~(2+)、Mg~(2+)的活化能分别为78.07和74.91 k J/mol。