Phase equilibria in the MgO–SiO<Subscript>2</Subscript>–ZrO<Subscript>2</Subscript> system

The results of study of phase equilibria in the MgO–SiO₂–ZrO₂ system at 1450–1550°C are reported. The studied system contains two eutectic points and six fields: (I) MgSiO₃ + SiO₂; (II) MgSiO₃ + ZrO₂; (III) ZrSiO₄ + SiO₂; (IV) MgSiO₃ + Mg₂SiO₄; (V) ZrO₂ + MgO; (VI) ZrSiO₄ + ZrO₂. The presence of fields (II) and (III) on the diagram shows that zircon in equilibrium with olivine and pyroxene crystallizes at very low concentrations of ZrO₂ in the system. This provides a solution for one of the most important problems in zirconology of dunites: the probability of the formation and preservation of zircon in the course of the formation and evolution of dunite.     

The high-pressure stability of chlorite and other hydrates in subduction mélanges: experiments in the system Cr<Subscript>2</Subscript>O<Subscript>3</Subscript>–MgO–Al<Subscript>2</Subscript>O<Subscript>3</Subscript>–SiO<Subscript>2</Subscript>–H<Subscript>2</Subscript>O

The solubility of chromium in chlorite as a function of pressure, temperature, and bulk composition was investigated in the system Cr₂O₃–MgO–Al₂O₃–SiO₂–H₂O, and its effect on phase relations evaluated. Three different compositions with X _Cr = Cr/(Cr + Al) = 0.075, 0.25, and 0.5 respectively, were investigated at 1.5–6.5 GPa, 650–900 °C. Cr-chlorite only occurs in the bulk composition with X _Cr = 0.075; otherwise, spinel and garnet are the major aluminous phases. In the experiments, Cr-chlorite coexists with enstatite up to 3.5 GPa, 800–850 °C, and with forsterite, pyrope, and spinel at higher pressure. At P > 5 GPa other hydrates occur: a Cr-bearing phase-HAPY (Mg_2.2Al_1.5Cr_0.1Si_1.1O₆(OH)₂) is stable in assemblage with pyrope, forsterite, and spinel; Mg-sursassite coexists at 6.0 GPa, 650 °C with forsterite and spinel and a new Cr-bearing phase, named 11.5 Å phase (Mg:Al:Si = 6.3:1.2:2.4) after the first diffraction peak observed in high-resolution X-ray diffraction pattern. Cr affects the stability of chlorite by shifting its breakdown reactions toward higher temperature, but Cr solubility at high pressure is reduced compared with the solubility observed in low-pressure occurrences in hydrothermal environments. Chromium partitions generally according to \(X_{\text{Cr}}^{\text{spinel}}\) ? \(X_{\text{Cr}}^{\text{opx}}\) > \(X_{\text{Cr}}^{\text{chlorite}}\) ≥ \(X_{\text{Cr}}^{\text{HAPY}}\) > \(X_{\text{Cr}}^{\text{garnet}}\). At 5 GPa, 750 °C (bulk with X _Cr = 0.075) equilibrium values are \(X_{\text{Cr}}^{\text{spinel}}\) = 0.27, \(X_{\text{Cr}}^{\text{chlorite}}\) = 0.08, \(X_{\text{Cr}}^{\text{garnet}}\) = 0.05; at 5.4 GPa, 720 °C \(X_{\text{Cr}}^{\text{spinel}}\) = 0.33, \(X_{\text{Cr}}^{\text{HAPY}}\) = 0.06, and \(X_{\text{Cr}}^{\text{garnet}}\) = 0.04; and at 3.5 GPa, 850 °C \(X_{\text{Cr}}^{\text{opx}}\) = 0.12 and \(X_{\text{Cr}}^{\text{chlorite}}\) = 0.07. Results on Cr–Al partitioning between spinel and garnet suggest that at low temperature the spinel- to garnet-peridotite transition has a negative slope of 0.5 GPa/100 °C. The formation of phase-HAPY, in assemblage with garnet and spinel, at pressures above chlorite breakdown, provides a viable mechanism to promote H₂O transport in metasomatized ultramafic mélanges of subduction channels.     

Effect of KCl on melting in the Mg<Subscript>2</Subscript>SiO<Subscript>4</Subscript>–MgSiO<Subscript>3</Subscript>–H<Subscript>2</Subscript>O system at 5 GPa

To examine the effect of KCl-bearing fluids on the melting behavior of the Earth’s mantle, we conducted experiments in the Mg₂SiO₄–MgSiO₃–H₂O and Mg₂SiO₄–MgSiO₃–KCl–H₂O systems at 5 GPa. In the Mg₂SiO₄–MgSiO₃–H₂O system, the temperature of the fluid-saturated solidus is bracketed between 1,200–1,250°C, and both forsterite and enstatite coexist with the liquid under supersolidus conditions. In the Mg₂SiO₄–MgSiO₃–KCl–H₂O systems with molar Cl/(Cl + H₂O) ratios of 0.2, 0.4, and 0.6, the temperatures of the fluid-saturated solidus are bracketed between 1,400–1,450°C, 1,550–1,600°C, and 1,600–1,650°C, respectively, and only forsterite coexists with liquid under supersolidus conditions. This increase in the temperature of the solidus demonstrates the significant effect of KCl on reducing the activity of H₂O in the fluid in the Mg₂SiO₄–MgSiO₃–H₂O system. The change in the melting residues indicates that the incongruent melting of enstatite (enstatite = forsterite + silica-rich melt) could extend to pressures above 5 GPa in KCl-bearing systems, in contrast to the behavior in the KCl-free system.     

Mechanisms and rates of quartz dissolution in melts in the CMAS (CaO–MgO–Al<Subscript>2</Subscript>O<Subscript>3</Subscript>–SiO<Subscript>2</Subscript>) system

The dissolution rate of minerals in silicate melts is generally assumed to be a function of the rate of mass transport of the released cations in the solvent. While this appears to be the case in moderately to highly viscous solvents, there is some evidence that the rate-controlling step may be different in very fluid, highly silica undersaturated melts such as basanites. In this study, convection-free experiments using solvent melts with silica activity from 0.185–0.56 and viscosity from 0.03–4.6 Pa s show that the dissolution rate is strongly dependent on the degree of superheating, silica activity and the viscosity of the solvent. Dissolution rates increase with increasing melt temperature and decreasing silica activity and viscosity. Quartz dissolution in melts with viscosity <0.59–1.9 Pa s and silica activity <0.47 is controlled by the rate of interface reaction as shown by the absence of steady state composition and silica saturation in the interface melts. Only in the most viscous melt with the highest silica activity is quartz dissolution controlled by the rate of diffusion in the melt and only after a long initiation time. The results of this study indicate that although a diffusion-based model may be applicable to dissolution in viscous magmas, a different approach that combines the interplay between the degree of undersaturation of the melt and its viscosity is required in very fluid melts.This revised version was published online September 2004 with a correction to Figure 8.     

Hydrosilicate liquids in the system rare-metal granite–Na<Subscript>2</Subscript>O–SiO<Subscript>2</Subscript>–H<Subscript>2</Subscript>O as accumulators of ore components at high pressure and temperature

Experimental investigations in the system rare-metal granite–Na₂O–SiO₂–H₂O with the addition of aqueous solutions containing Rb, Cs, Sn, W, Mo, and Zn at 600°C and 1.5 kbar showed that the typical elements of rare-metal granites (Li, Rb, Cs, Be, Nb, and Ta) are preferentially concentrated in hydrosilicate liquids coexisting with aqueous fluid. The same behavior is characteristic of Zn and Sn, the minerals of which are usually formed under hydrothermal conditions. In contrast, Mo and W are weakly extracted by hydrosilicate liquids and almost equally distributed between them and aqueous fluids. Liquids similar to those described in this paper are formed during the final stages of magmatic crystallization in granite and granitepegmatite systems. The formation of hydrosilicate liquids in late magmatic and postmagmatic processes will be an important factor controlling the redistribution of metal components between residual magmatic melts, minerals, and aqueous fluids and, consequently, the mobility of these components in fluid-saturated magmatic systems enriched in rare metals.     

Mineral-fluid equilibria in the system CaO–MgO–SiO<Subscript>2</Subscript>–H<Subscript>2</Subscript>O–CO<Subscript>2</Subscript>–NaCl and the record of reactive fluid flow in contact metamorphic aureoles

Phase equilibria in the system CaO–MgO–SiO₂–CO₂–H₂O–NaCl are calculated to illustrate phase relations in metacarbonates over a wide-range of P–T–X[H₂O–CO₂–NaCl] conditions. Calculations are performed using the equation of state of Duan et al. (Geochim Cosmochim Acta 59:2869–2882, 1995) for H₂O–CO₂–NaCl fluids and the internally consistent data set of Gottschalk (Eur J Mineral 9:175–223, 1997) for thermodynamic properties of solids. Results are presented in isothermal-isobarical plots showing stable mineral assemblages as a function of fluid composition. It is shown that in contact-metamorphic P–T regimes the presence of very small concentrations of NaCl in the fluid causes almost all decarbonation reactions to proceed within the two fluid solvus of the H₂O–CO₂–NaCl system. Substantial flow of magma-derived fluids into marbles has been documented for many contact aureoles by shifts in stable isotope geochemistry of the host rocks and by the progress of volatile-producing mineral reactions controlled by fluid compositions. Time-integrated fluid fluxes have been estimated by combining fluid advection/dispersion models with the spatial arrangement of mineral reactions and isotopic resetting. All existing models assume that minerals react in the presence of a single phase H₂O–CO₂ fluid and do not allow for the effect that fluid immiscibility has on the flow patterns. It is shown that fluids emanating from calc-alkaline melts that crystallize at shallow depths are brines. Their salinity may vary depending mainly on pressure and fraction of crystallized melt. Infiltration-driven decarbonation reactions in the host rocks inevitably proceed at the boundaries of the two fluid solvus where the produced CO₂ is immiscible and may separate from the brine as a low salinity, low density H₂O–CO₂ fluid. Most parameters of fluid–rock interaction in contact aureoles that are derived from progress of mineral reactions and stable isotope resetting are probably incorrect because fluid phase separation is disregarded.     

Measurements and calculations of solid-liquid equilibria in the ternary systems NaBr–Na<Subscript>2</Subscript>SO<Subscript>4</Subscript>–H<Subscript>2</Subscript>O and KBr–K<Subscript>2</Subscript>SO<Subscript>4</Subscript>–H<Subscript>2</Subscript>O at 348 K

According to the compositions of the underground brine resources in the west of Sichuan Basin, solubilities of the ternary systems NaBr–Na₂SO₄–H₂O and KBr–K₂SO₄–H₂O were investigated by isothermal method at 348 K. The equilibrium solid phases, solubilities of salts, and densities of the solutions were determined. On the basis of the experimental data, the phase diagrams and the density-composition diagrams were plotted. In the two ternary systems, the phase diagrams consist of two univariant curves, one invariant point and two crystallization fields. Neither solid solution nor double salts were found. The equilibrium solid phases in the ternary system NaBr–Na₂SO₄–H₂O are NaBr and Na₂SO₄, and those in the ternary system KBr–K₂SO₄–H₂O are KBr and K₂SO₄. Using the solubilities data of the two ternary subsystems at 348 K, mixing ion-interaction parameters of Pitzer’s equation θxxx, Ψxxx and Ψxxx were fitted by multiple linear regression method. Based on the chemical model of Pitzer’s electrolyte solution theory, the solubilities of phase equilibria in the two ternary systems NaBr–Na₂SO₄–H₂O and KBr–K₂SO₄–H₂O were calculated with corresponding parameters. The calculation diagrams were plotted. The results showed that the calculated values have a good agreement with experimental data.     

Attikaite,Ca<Subscript>3</Subscript>Cu<Subscript>2</Subscript>Al<Subscript>2</Subscript>(AsO<Subscript>4</Subscript>)<Subscript>4</Subscript>(OH)<Subscript>4</Subscript> · 2H<Subscript>2</Subscript>O,a new mineral species

Attikaite, a new mineral species, has been found together with arsenocrandalite, arsenogoyazite, conichalcite, olivenite, philipsbornite, azurite, malachite, carminite, beudantite, goethite, quartz, and allophane at the Christina Mine No. 132, Kamareza, Lavrion District, Attiki Prefecture (Attika), Greece. The mineral is named after the type locality. It forms spheroidal segregations (up to 0.3 mm in diameter) consisting of thin flexible crystals up to 3 × 20 × 80 μm in size. Its color is light blue to greenish blue, with a pale blue streak. The Mohs’ hardness is 2 to 2.5. The cleavage is eminent mica-like parallel to {001}. The density is 3.2(2) g/cm³ (measured in heavy liquids) and 3.356 g/cm³ (calculated). The wave numbers of the absorption bands in the infrared spectrum of attikaite are (cm^?1; sh is shoulder; w is a weak band): 3525sh, 3425, 3180, 1642, 1120w, 1070w, 1035w, 900sh, 874, 833, 820, 690w, 645w, 600sh, 555, 486, 458, and 397. Attikaite is optically biaxial, negative, α = 1.642(2), β = γ = 1.644(2) (X = c) 2V _means = 10(8)°, and 2V _calc = 0°. The new mineral is microscopically colorless and nonpleochroic. The chemical composition (electron microprobe, average over 4 point analyses, wt %) is: 0.17 MgO, 17.48 CaO, 0.12 FeO, 16.28 CuO, 10.61 Al₂O₃, 0.89 P₂O₅, 45.45 As₂O₅, 1.39 SO₃, and H₂O (by difference) 7.61, where the total is 100.00. The empirical formula calculated on the basis of (O,OH,H₂O)₂₂ is: Ca_2.94Cu _1.93 ²⁺ Al_1.97Mg_0.04Fe _0.02 ²⁺ [(As_3.74S_0.16P_0.12)_Σ4.02O_16.08](OH)_3.87 · 2.05H₂ O. The simplified formula is Ca₃Cu₂Al₂(AsO₄)₄(OH)₄ · 2H₂O. Attikaite is orthorhombic, space group Pban, Pbam or Pba2; the unit-cell dimensions are a = 10.01(1), b = 8.199(5), c = 22.78(1) Å, V = 1870(3) ų, and Z = 4. In the result of the ignition of attikaite for 30 to 35 min at 128–140°, the H₂O bands in the IR spectrum disappear, while the OH-group band is not modified; the weight loss is 4.3%, which approximately corresponds to two H₂O molecules per formula; and parameter c decreases from 22.78 to 18.77 Å. The strongest reflections in the X-ray powder diffraction pattern [d, Å (I, %)((hkl)] are: 22.8(100)(001), 11.36(60)(002), 5.01(90)(200), 3.38(5)(123, 205), 2.780(70)(026), 2.682(30)(126), 2.503(50)(400), 2.292(20)(404). The type material of attikaite is deposited in the Fersman Mineralogical Museum, Russian Academy of Sciences, Moscow. The registration number is 3435/1.     

Behaviour of Fe<Subscript>4</Subscript>O<Subscript>5</Subscript>–Mg<Subscript>2</Subscript>Fe<Subscript>2</Subscript>O<Subscript>5</Subscript> solid solutions and their relation to coexisting Mg–Fe silicates and oxide phases

Experiments at high pressures and temperatures were carried out (1) to investigate the crystal-chemical behaviour of Fe₄O₅–Mg₂Fe₂O₅ solid solutions and (2) to explore the phase relations involving (Mg,Fe)₂Fe₂O₅ (denoted as O₅-phase) and Mg–Fe silicates. Multi-anvil experiments were performed at 11–20 GPa and 1100–1600 °C using different starting compositions including two that were Si-bearing. In Si-free experiments the O₅-phase coexists with Fe₂O₃, hp-(Mg,Fe)Fe₂O₄, (Mg,Fe)₃Fe₄O₉ or an unquenchable phase of different stoichiometry. Si-bearing experiments yielded phase assemblages consisting of the O₅-phase together with olivine, wadsleyite or ringwoodite, majoritic garnet or Fe³⁺-bearing phase B. However, (Mg,Fe)₂Fe₂O₅ does not incorporate Si. Electron microprobe analyses revealed that phase B incorporates significant amounts of Fe²⁺ and Fe³⁺ (at least ~?1.0 cations Fe per formula unit). Fe-L_2,3-edge energy-loss near-edge structure spectra confirm the presence of ferric iron [Fe³⁺/Fe_tot?=?~?0.41(4)] and indicate substitution according to the following charge-balanced exchange: [4]Si⁴⁺?+?[6]Mg²⁺?=?2Fe³⁺. The ability to accommodate Fe²⁺ and Fe³⁺ makes this potential “water-storing” mineral interesting since such substitutions should enlarge its stability field. The thermodynamic properties of Mg₂Fe₂O₅ have been refined, yielding H°_1bar,298?=???1981.5 kJ mol^??1. Solid solution is complete across the Fe₄O₅–Mg₂Fe₂O₅ binary. Molar volume decreases essentially linearly with increasing Mg content, consistent with ideal mixing behaviour. The partitioning of Mg and Fe²⁺ with silicates indicates that (Mg,Fe)₂Fe₂O₅ has a strong preference for Fe²⁺. Modelling of partitioning with olivine is consistent with the O₅-phase exhibiting ideal mixing behaviour. Mg–Fe²⁺ partitioning between (Mg,Fe)₂Fe₂O₅ and ringwoodite or wadsleyite is influenced by the presence of Fe³⁺ and OH incorporation in the silicate phases.     

Elastic and thermoelastic properties of synthetic Ca<Subscript>2</Subscript>MgSi<Subscript>2</Subscript>O<Subscript>7</Subscript> (åkermanite) and Ca<Subscript>2</Subscript>ZnSi<Subscript>2</Subscript>O<Subscript>7</Subscript> (hardystonite)

Elastic and thermoelastic constants of large single crystals of Ca₂MgSi₂O₇ and Ca₂ZnSi₂O₇ have been derived from ultrasonic resonance frequencies of plane-parallel plates and their shift upon variation of temperature, respectively. In addition, coefficients of thermal expansion and dielectric constants were determined. Both species possess quite similar properties. As observed in other isotypic magnesium and zinc compounds, the mean elastic stiffness and the deviation from the Cauchy relations are significantly larger in the zinc compound, due to a covalent contribution of the Zn–O bond. Positive thermoelastic constants T₄₄ and T₆₆ in Ca₂MgSi₂O₇ allow temperature-independent ultrasonic generators and oscillators to be manufactured.     

Structural changes and valence states in the MgCr<Subscript>2</Subscript>O<Subscript>4</Subscript>–FeCr<Subscript>2</Subscript>O<Subscript>4</Subscript> solid solution series

The influence on the structure of Fe²⁺ Mg substitution was studied in synthetic single crystals belonging to the MgCr₂O₄–FeCr₂O₄ series produced by flux growth at 900–1200 °C in controlled atmosphere. Samples were analyzed by single-crystal X-ray diffraction, electron microprobe analyses, optical absorption-, infrared- and Mössbauer spectroscopy. The Mössbauer data show that iron occurs almost exclusively as ^IVFe²⁺. Only minor Fe³⁺ (<0.005 apfu) was observed in samples with very low total Fe. Optical absorption spectra show that chromium with few exceptions is present as a trivalent cation at the octahedral site. Additional absorption bands attributable to Cr²⁺ and Cr³⁺ at the tetrahedral site are evident in spectra of end-member magnesiochromite and solid-solution crystals with low ferrous contents. Structural parameters a₀, u and T–O increase with chromite content, while the M–O bond distance remains nearly constant, with an average value equal to 1.995(1) Å corresponding to the Cr³⁺ octahedral bond distance. The ideal trend between cell parameter, T–O bond length and Fe²⁺ content (apfu) is described by the following linear relations: a₀=8.3325(5) + 0.0443(8)Fe²⁺ (Å) and T–O=1.9645(6) + 0.033(1)Fe²⁺ (Å) Consequently, Fe²⁺ and Mg tetrahedral bond lengths are equal to 1.998(1) Å and 1.965(1) Å, respectively.     

Calcioolivine, γ-Ca<Subscript>2</Subscript>SiO<Subscript>4</Subscript>, an old and New Mineral species

Calcioolivine has been included into the MDI mineral database in the list of grandfathered minerals. Its history, together with related artificial compounds, is extremely complex: various minerals and compounds received this name, including natural orthorhombic Ca orthosilicate. In this paper, the crystal structure and properties of natural calcioolivine are described for the first time. The new mineral has been found at Mt. Lakargi, Upper Chegem Plateau, the northern Caucasus, Kabarda-Balkaria Republic, Russia. It has been identified in skarnified, primary carbonate xenoliths entrained by middle to late Pliocene silicic ignimbrites of the Upper Chegem caldera. These xenoliths of a few centimeters to a few meters in size are located close to the volcanic vent. Calcioolivine rims relics of larnite and occurs as relict grains among crystals of spurrite, rondorfite, wadalite or secondary hillebrandite, afwillite, thaumasite, and ettringite. Hillebrandite is the major product of alteration of calcioolivine; larnite is relatively more resistant to low-temperature alteration. Spurrite, larnite, tilleyite, kilchoanite, cuspidine, wadalite, rondorfite, reinhardbraunsite, lakargiite (CaZrO₃), members of ellestadite series, afwillite, ettringite, katoite, and thaumasite are associated minerals. It is inferred that calcioolivine has been produced as a result of interaction of host carbonate rocks in xenoliths with volcanic lava and gases during eruption. The name calcioolivine was approved by the Commission on New Minerals and Mineral Names, International Mineralogical Association, September 6, 2007 (no. 07-B).     

Titanium solubility in olivine in the system TiO<Subscript>2</Subscript>–MgO–SiO<Subscript>2</Subscript>: no evidence for an ultra-deep origin of Ti-bearing olivine

The finding of ilmenite rods in olivine from orogenic peridotites has sparked a discussion about the processes of incorporation and exsolution of titanium in olivine. We have experimentally investigated the solubility of Ti in olivine as a function of composition, temperature and pressure in the synthetic TiO₂–MgO–SiO₂ system. Experiments at atmospheric pressure in the temperature range 1,200–1,500°C showed that the highest concentration of TiO₂ is obtained when olivine coexists with spinel (Mg₂TiO₄). The amount of TiO₂ in olivine in the assemblages olivine + spinel + periclase and olivine + spinel + ilmenite at 1,500°C was 1.25 wt.%. Changes in the coexisting phases and decreasing temperature result in a significant reduction of the Ti solubility. Olivine coexisting with pseudobrookite (MgTi₂O₅) and a Ti–Si-rich melt at 1,500°C displays a fourfold lower TiO₂ content than when buffered with spinel. A similar decrease in solubility is obtained by a decrease in temperature to 1,200°C. There is a negative correlation between Ti and Si and no correlation between Ti and Mg in Ti-bearing olivine. Together with the established phase relations this suggests that there is a direct substitution of Ti for Si at these temperatures, such that the substituting component has the stoichiometry Mg₂TiO₄. The unit cell volume of olivine increases systematically with increasing TiO₂ content demonstrating that the measured TiO₂ contents in olivine are not caused by micro-inclusions but by incorporation of Ti in the olivine structure. Least squares fitting of 20 olivine unit cell volumes against the Ti content yield the relation: V (ų)=290.12(1) + 23.67(85) N_Ti. The partial molar volume of end-member Mg₂TiO₄ olivine (N_Ti=1) is thus 47.24±0.13 cm³. The change of the Ti solubilty in olivine coexistent with rutile and orthopyroxene with pressure was investigated by piston cylinder experiments at 1,400°C from 15 to 55 kbar. There is no increase in TiO₂ contents with pressure and in all the experiments olivine contains ~0.2 wt.% TiO₂. Moreover, a thermodynamic analysis indicates that Ti contents of olivine coexisting with rutile and orthopyroxene should decrease rather than increase with increasing pressure. These data indicate that the ilmenite exsolution observed in some natural olivine does not signify an ultra-deep origin of peridotite massifs.     

Thermophysical properties of the α–β–γ polymorphs of Mg<Subscript>2</Subscript>SiO<Subscript>4</Subscript>: a computational study

Thermophysical properties of the various polymorphs (i.e. α-, β- and γ) of Mg₂SiO₄ were computed with the CRYSTAL06 code within the framework of CO-LCAO-GTF approach by using the hybrid B3LYP density functional method. Potential wells were calculated through a symmetry preserving, variable cell-shape structure relaxation procedure. Vibrational frequencies were computed at the long-wavelength limit corresponding to the center of the Brillouin zone (k → 0). Thermodynamic properties were estimated through a semiclassical approach that combines B3LYP vibrational frequencies for optic modes and the Kieffer’s model for the dispersion relation of acoustic modes. All computed values except volume (i.e. electronic energy, zero point energy, optical vibrational modes, thermal corrections to internal energy, standard state enthalpy and Gibbs free energy of reaction, bulk modulus and its P and T derivatives, entropy, C _V, C _P) are consistent with available experimental data and/or reasonable estimates. Volumes are slightly overestimated relative to those determined directly by X-ray diffraction. A set of optimized volumetric properties that are consistent with the other semiclassical properties of the phases α, β and γ have been derived by optimization procedure such that the calculated boundaries for the α/β and β/γ equilibria have the best overall agreement with the experimental data for these transitions. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.

Optical spectroscopic study of synthetic NaScSi<Subscript>2</Subscript>O<Subscript>6</Subscript>–CaNiSi<Subscript>2</Subscript>O<Subscript>6</Subscript> pyroxenes at normal and high pressures

Six synthetic NaScSi₂O₆–CaNiSi₂O₆ pyroxenes were studied by optical absorption spectroscopy. Five of them of intermediate (Na_1−x , Ca _x )(Sc_1−x , Ni _x )Si₂O₆ compositions show spectra typical of Ni²⁺ in octahedral coordination, more precise Ni²⁺ at the M1 site of the pyroxene structure. The common feature of all spectra is three broad absorption bands with maxima around 8,000, 13,000 and 24,000 cm⁻¹ assigned to ³ A _2g → ³ T _2g, ³ A _2g → ³ T _1g and →³ T _1g (³ P) electronic spin-allowed transitions of ^VINi²⁺. A weak narrow peak at ∼14,400 cm⁻¹ is assigned to the spin-forbidden ³ A _2g → ¹ T _2g (¹ D) transition of Ni²⁺. Under pressure the spin-allowed bands shift to higher energies and change in intensity. The octahedral compression modulus, calculated from the shift of the ³ A _2g → ³ T _2g band in the (Na_0.7Ca_0.3)(Sc_0.7Ni_0.3)Si₂O₆ pyroxene is evaluated as 85±20 GPa. The Racah parameter B of Ni²⁺(M1) is found gradually changing from ∼919 cm⁻¹ at ambient pressure to ∼890 cm⁻¹ at 6.18 GPa. The Ni end-member pyroxene [(Ca_0.93 Ni_0.07)NiSi₂O₆] has a spectrum different from all others. In addition to the above mentioned bands of Ni²⁺(M1) it displays several new relatively intense and broad extra bands, which were attributed to electronic transitions of Ni²⁺ at the M2 site. In difference to CaO₈ polyhedron geometry of an eightfold coordination, Ni²⁺(M2)O₈ polyhedra are assumed to be relatively large distorted octahedra. Due to different distortions and different compressibilities of the M1 and M2 sites the Ni²⁺(M1)- and Ni²⁺(M2)-bands display rather different pressure-induced behaviors, becoming more resolved in the high-pressure spectra than in that measured at atmospheric pressure. The octahedral compression modulus of Ni²⁺(M1) in this end-member pyroxene is evaluated as 150 ± 25 GPa, which is noticeably larger than in Ni_0.3 pyroxene. This is due to a smaller size and, thus, a stiffer character of Ni²⁺(M1)O₆ octahedron in the (Ca_0.93Ni_0.07)NiSi₂O₆ pyroxene compared to (Na_0.7Ca_0.3)(Sc_0.7Ni_0.3)Si₂O₆.

Single-crystal X-ray diffraction study of cation distribution in MgAl<Subscript>2</Subscript>O<Subscript>4</Subscript>–MgFe<Subscript>2</Subscript>O<Subscript>4</Subscript> spinel solid solution

Synthesis experiments in the system MgAl₂O₄–MgFe₂O₄ [MgAl_2–xFe_xO₄ (0 x 2)] were carried out using a PbF₂ flux. The crystalline products synthesized in the compositional range of 0.6 <x 1.2 consisted of two spinel phases, whereas those synthesized in the compositional ranges of 0.0 x 0.6 and 1.2 < x 2.0 crystallized as single spinel phases. Structure refinements of the spinel single crystals, which grew in the ranges of 0.0 x 0.6 and 1.2 < x 2.0, show that the degree of randomness of cation distribution between A and B sites increases as x approaches the two-phase region. This means that the degree of the size mismatch among Mg²⁺, Fe³⁺ and Al³⁺occupying each equivalent mixing site increases as x approaches the two-phase region. Consequently, if the coexistence of two spinels observed in the intermediate compositions reveals the existence of a miscibility gap at low temperatures, this increase in the degree of the size mismatch among the three cations is suggested as a factor of energetic destabilization to form the miscibility gap.     

Structural characterization of high-pressure C–Na<Subscript>2</Subscript>Si<Subscript>2</Subscript>O<Subscript>5</Subscript> by single-crystal diffraction and <Superscript>29</Superscript>Si MAS NMR

Single crystals of C–Na₂Si₂O₅ have been synthesized from the hydrothermal recrystallization of a glass. The title compound is monoclinic, space group P2₁/c with Z= 8 and unit-cell parameters a= 4.8521 (4)Å, b=23.9793(16)Å, c=8.1410(6)Å, β=90.15(1)° and V=947.2(2)ų. The structure has been determined by direct methods and belongs to the group of phyllosilicates. It is based on layers of tetrahedra with elliptically six-membered rings in chair conformation. The sequence of directedness within a single ring is UDUDUD. The sheets are parallel to (010) with linking sodium cations in five- and sixfold coordination. Concerning the shape and the conformation of the rings, C–Na₂Si₂O₅ is closely related to β-Na₂Si₂O₅. However, both structures differ in the stacking sequences of the layers. A possible explanation for the frequently observed polysynthetic twinning of phase C is presented. In the ²⁹Si MAS-NMR spectrum of C–Na₂Si₂O₅ four well-resolved lines of equal intensity are observed at ?86.0, ?86.3, ?87.4, and ?88.2?ppm. The narrow range of isotropic chemical shifts reflects the great similarity of the environments of the different Si sites. This lack of pronounced differences in geometry renders a reliable assignment of the resonance lines to the individual sites on the basis of known empiric correlations and geometrical features impossible.     

Fluid–mineral reactions and melting of orthopyroxene–cordierite–biotite gneiss in the presence of H<Subscript>2</Subscript>O-CO<Subscript>2</Subscript>-NaCl and H<Subscript>2</Subscript>O-CO<Subscript>2</Subscript>-KCl fluids under parameters of granulite-facies metamorphism

Reactions and partial melting of peraluminous rocks in the presence of H₂O-CO₂–salt fluids under parameters of granulite-facies metamorphism were modeled in experiments on interaction between orthopyroxene–cordierite–biotite–plagioclase–quartz metapelite with H₂O, H₂O-CO₂, H₂O-CO₂-NaCl, and H₂O-CO₂-KCl fluids at 600 MPa and 850°C. Rock melting in the presence of H₂O and equimolar H₂O-CO₂ fluids generates peraluminous (A/CNK¹ > 1.1) melts whose composition corresponds to magnesian calcic or calc–alkaline S-type granitoids. The melts are associated with peritectic phases: magnesian spinel and orthopyroxene containing up to 9 wt % Al₂O₃. In the presence of H₂O-CO₂-NaCl fluid, cordierite and orthopyroxene are replaced by the association of K-Na biotite, Na-bearing gedrite, spinel, and albite. The Na₂O concentrations in the biotite and gedrite are functions of the NaCl concentrations in the starting fluid. Fluids of the composition H₂O-CO₂-KCl induce cordierite replacement by biotite with corundum and spinel and by these phases in association with potassium feldspar at X _KCl = 0.02 in the fluid. When replaced by these phases, cordierite is excluded from the melting reactions, and the overall melting of the metapelite is controlled by peritectic reactions of biotite and orthopyroxene with plagioclase and quartz. These reactions produce such minerals atypical of metapelites as Ca-Na amphibole and clinopyroxene. The compositions of melts derived in the presence of salt-bearing fluids are shifted toward the region with A/CNK < 1.1, as is typical of so-called peraluminous granites of type I. An increase in the concentrations of salts in the fluids leads to depletion of the melts in Al₂O₃ and CaO and enrichment in alkalis. These relations suggest that the protoliths of I-type peraluminous granites might have been metapelites that were melted when interacting with H₂O-CO₂-salt fluids. The compositions of the melts can evolve from those with A/CNK > 1.1 (typical of S-type granites) toward those with A/CNK = 1.0–1.1 in response to an increase in the concentrations of alkali salts in the fluids within a few mole percent. Our experiments demonstrate that the origin of new mineral assemblages in metapelite in equilibrium with H₂O-CO₂-salt fluids is controlled by the activities of alkaline components, while the H₂O and CO₂ activities play subordinate roles. This conclusion is consistent with the results obtained by simulating metapelite mineral assemblages by Gibbs free energy minimization (using the PERPE_X software), as shown in log(\({a_{{H_2}O}}\))–log(\({a_{N{a_2}O}}\)) and log(\({a_{{H_2}O}}\))–log(\({a_{{K_2}O}}\)) diagrams.     

Trace element partitioning between orthopyroxene and anhydrous silicate melt on the lherzolite solidus from 1.1 to 3.2 GPa and 1,230 to 1,535°C in the model system Na<Subscript>2</Subscript>O–CaO–MgO–Al<Subscript>2</Subscript>O<Subscript>3</Subscript>–SiO<Subscript>2</Subscript>

We determined experimentally the Nernst distribution coefficient between orthopyroxene and anhydrous silicate melt for trace elements i in the system Na₂O–CaO–MgO–Al₂O₃–SiO₂ (NCMAS) along the dry model lherzolite solidus from 1.1 GPa/1,230°C up to 3.2 GPa/1,535°C in a piston cylinder apparatus. Major and trace element composition of melt and orthopyroxene were determined with a combination of electron microprobe and ion probe analyses. We provide partitioning data for trace elements Li, Be, B, K, Sc, Ti, V, Cr, Co, Ni, Rb, Sr, Y, Zr, Nb, Cs, Ba, La, Ce, Sm, Nd, Yb, Lu, Hf, Ta, Pb, U, and Th. The melts were chosen to be boninitic at 1.1 and 2.0 GPa, picritic at 2.3 GPa and komatiitic at 2.7 and 3.2 GPa. Orthopyroxene is Tschermakitic with 8 mol% Mg-Tschermaks MgAl[AlSiO₆] at 1.1 GPa while at higher pressure it has 18–20 mol%. The rare earth elements show a continuous, significant increase in compatibility with decreasing ionic radius from D _La^opx−melt ∼ 0.0008 to D _Lu^opx−melt ∼ 0.15. For the high-field-strength elements compatibility increases from D _Th^opx−melt ∼ 0.001 through D _Nb^opx−melt ∼ 0.0015, D _U^opx−melt ∼ 0.002, D _Ta^opx−melt ∼ 0.005, D _Zr^opx−melt ∼ 0.02 and D _Hf^opx−melt ∼ 0.04 to D _Ti^opx−melt ∼ 0.14. From mathematical and graphical fits we determined best-fit values for D ₀^M1, D ₀^M2, r ₀^M1, r ₀^M2, E ₀^M1, and E ₀^M2 for the two different M sites in orthopyroxene according to the lattice strain model and calculated the intracrystalline distribution between M1 and M2. Our data indicate extreme intracrystalline fractionation for most elements in orthopyroxene; for the divalent cations D _i ^M2−M1 varies by three orders of magnitude between D _Co^M2−M1 = 0.00098–0.00919 and D _Ba^M2−M1 = 2.3–28. Trivalent cations Al and Cr almost exclusively substitute on M1 while the other trivalent cations substitute on M2; D _La^M2−M1 reaches extreme values between 6.5 × 10⁷ and 1.4 × 10¹⁶. Tetravalent cations Ti, Hf, and Zr almost exclusively substitute on M1 while U and Th exclusively substitute on M2. Our new comprehensive data set can be used for polybaric-polythermal melting models along the Earth’s mantle solidus. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.