Experimental determination of liquidus H<Subscript>2</Subscript>O contents of haplogranite at deep-crustal conditions

The liquidus water content of a haplogranite melt at high pressure (P) and temperature (T) is important, because it is a key parameter for constraining the volume of granite that could be produced by melting of the deep crust. Previous estimates based on melting experiments at low P (≤0.5 GPa) show substantial scatter when extrapolated to deep crustal P and T (700–1000 °C, 0.6–1.5 GPa). To improve the high-P constraints on H₂O concentration at the granite liquidus, we performed experiments in a piston–cylinder apparatus at 1.0 GPa using a range of haplogranite compositions in the albite (Ab: NaAlSi₃O₈)—orthoclase (Or: KAlSi₃O₈)—quartz (Qz: SiO₂)—H₂O system. We used equal weight fractions of the feldspar components and varied the Qz between 20 and 30 wt%. In each experiment, synthetic granitic composition glass + H₂O was homogenized well above the liquidus T, and T was lowered by increments until quartz and alkali feldspar crystalized from the liquid. To establish reversed equilibrium, we crystallized the homogenized melt at the lower T and then raised T until we found that the crystalline phases were completely resorbed into the liquid. The reversed liquidus minimum temperatures at 3.0, 4.1, 5.8, 8.0, and 12.0 wt% H₂O are 935–985, 875–900, 775–800, 725–775, and 650–675 °C, respectively. Quenched charges were analyzed by petrographic microscope, scanning electron microscope (SEM), X-ray diffraction (XRD), and electron microprobe analysis (EMPA). The equation for the reversed haplogranite liquidus minimum curve for Ab_36.25Or_36.25Qz_27.5 (wt% basis) at 1.0 GPa is \(T = - 0.0995 w_{{{\text{H}}_{ 2} {\text{O}}}}^{ 3} + 5.0242w_{{{\text{H}}_{ 2} {\text{O}}}}^{ 2} - 88.183 w_{{{\text{H}}_{ 2} {\text{O}}}} + 1171.0\) for \(0 \le w_{{{\text{H}}_{ 2} {\text{O}}}} \le 17\) wt% and \(T\) is in °C. We present a revised \(P - T\) diagram of liquidus minimum H₂O isopleths which integrates data from previous determinations of vapor-saturated melting and the lower pressure vapor-undersaturated melting studies conducted by other workers on the haplogranite system. For lower H₂O (<5.8 wt%) and higher temperature, our results plot on the high end of the extrapolated water contents at liquidus minima when compared to the previous estimates. As a consequence, amounts of metaluminous granites that can be produced from lower crustal biotite–amphibole gneisses by dehydration melting are more restricted than previously thought.     

Thermal equation of state of CaIrO<Subscript>3</Subscript> post-perovskite

The pressure–volume–temperature (P–V–T) relation of CaIrO₃ post-perovskite (ppv) was measured at pressures and temperatures up to 8.6 GPa and 1,273 K, respectively, with energy-dispersive synchrotron X-ray diffraction using a DIA-type, cubic-anvil apparatus (SAM85). Unit-cell dimensions were derived from the Le Bail full profile refinement technique, and the results were fitted using the third-order Birth-Murnaghan equation of state. The derived bulk modulus \( K_{T0} \) at ambient pressure and temperature is 168.3 ± 7.1 GPa with a pressure derivative \( K_{T0}^{\prime } \) = 5.4 ± 0.7. All of the high temperature data, combined with previous experimental data, are fitted using the high-temperature Birch-Murnaghan equation of state, the thermal pressure approach, and the Mie-Grüneisen-Debye formalism. The refined thermoelastic parameters for CaIrO₃ ppv are: temperature derivative of bulk modulus \( (\partial K_{T} /\partial T)_{P} \) = ?0.038 ± 0.011 GPa K^?1, \( \alpha K_{T} \) = 0.0039 ± 0.0001 GPa K^?1, \( \left( {\partial K_{T} /\partial T} \right)_{V} \) = ?0.012 ± 0.002 GPa K^?1, and \( \left( {\partial^{2} P/\partial T^{2} } \right)_{V} \) = 1.9 ± 0.3 × 10^?6 GPa² K^?2. Using the Mie-Grüneisen-Debye formalism, we obtain Grüneisen parameter \( \gamma_{0} \) = 0.92 ± 0.01 and its volume dependence q = 3.4 ± 0.6. The systematic variation of bulk moduli for several oxide post-perovskites can be described approximately by the relationship K _T0  = 5406.0/V(molar) + 5.9 GPa.     

The thermal expansion and crystal structure of mirabilite (Na<Subscript>2</Subscript>SO<Subscript>4</Subscript>·10D<Subscript>2</Subscript>O) from 4.2 to 300 K,determined by time-of-flight neutron powder diffraction

We have collected high resolution neutron powder diffraction patterns from Na₂SO₄·10D₂O over the temperature range 4.2–300 K following rapid quenching in liquid nitrogen, and over a series of slow warming and cooling cycles. The crystal is monoclinic, space-group P2₁/c (Z = 4) with a = 11.44214(4) Å, b = 10.34276(4) Å, c = 12.75486(6) Å, β = 107.847(1)°, and V = 1436.794(8) Å³ at 4.2 K (slowly cooled), and a = 11.51472(6) Å, b = 10.36495(6) Å, c = 12.84651(7) Å, β = 107.7543(1)°, V = 1460.20(1) Å³ at 300 K. Structures were refined to R _P (Rietveld powder residual, \( R_{P} = {{\sum {\left| {I_{\text{obs}} - I_{\text{calc}} } \right|} } \mathord{\left/ {\vphantom {{\sum {\left| {I_{\text{obs}} - I_{\text{calc}} } \right|} } {\sum {I_{\text{obs}} } }}} \right. \kern-\nulldelimiterspace} {\sum {I_{\text{obs}} } }} \)) better than 2.5% at 4.2 K (quenched and slow cooled), 150 and 300 K. The sulfate disorder observed previously by Levy and Lisensky (Acta Cryst B34:3502–3510, 1978) was not present in our specimen, but we did observe changes with temperature in deuteron occupancies of the orientationally disordered water molecules coordinated to Na. The temperature dependence of the unit-cell volume from 4.2 to 300 K is well represented by a simple polynomial of the form V = ? 4.143(1) × 10^?7 T ³ + 0.00047(2) T² ? 0.027(2) T + 1437.0(1) Å³ (R ² = 99.98%). The coefficient of volume thermal expansion, α _V , is positive above 40 K, and displays a similar magnitude and temperature dependence to α _V in deuterated epsomite and meridianiite. The relationship between the magnitude and orientation of the principal axes of the thermal expansion tensor and the main structural elements are discussed; freezing in of deuteron disorder in the quenched specimen affects the thermal expansion, manifested most obviously as a change in the behaviour of the unit-cell parameter β.     

Thermoelastic properties of chromium oxide Cr<Subscript>2</Subscript>O<Subscript>3</Subscript> (eskolaite) at high pressures and temperatures

A new synchrotron X-ray diffraction study of chromium oxide Cr₂O₃ (eskolaite) with the corundum-type structure has been carried out in a Kawai-type multi-anvil apparatus to pressure of 15 GPa and temperatures of 1873 K. Fitting the Birch–Murnaghan equation of state (EoS) with the present data up to 15 GPa yielded: bulk modulus (K _0,T0), 206 ± 4 GPa; its pressure derivative K′_0,T , 4.4 ± 0.8; (?K _0,T /?T) = ?0.037 ± 0.006 GPa K^?1; a = 2.98 ± 0.14 × 10^?5 K^?1 and b = 0.47 ± 0.28 × 10^?8 K^?2, where α _0,T  = a + bT is the volumetric thermal expansion coefficient. The thermal expansion of Cr₂O₃ was additionally measured at the high-temperature powder diffraction experiment at ambient pressure and α _0,T0 was determined to be 2.95 × 10^?5 K^?1. The results indicate that coefficient of the thermal expansion calculated from the EoS appeared to be high-precision because it is consistent with the data obtained at 1 atm. However, our results contradict α ₀ value suggested by Rigby et al. (Brit Ceram Trans J 45:137–148, 1946) widely used in many physical and geological databases. Fitting the Mie–Grüneisen–Debye EoS with the present ambient and high-pressure data yielded the following parameters: K _0,T0 = 205 ± 3 GPa, K′_0,T  = 4.0, Grüneisen parameter (γ ₀) = 1.42 ± 0.80, q = 1.82 ± 0.56. The thermoelastic parameters indicate that Cr₂O₃ undergoes near isotropic compression at room and high temperatures up to 15 GPa. Cr₂O₃ is shown to be stable in this pressure range and adopts the corundum-type structure. Using obtained thermoelastic parameters, we calculated the reaction boundary of knorringite formation from enstatite and eskolaite. The Clapeyron slope (with \({\text{d}}P/{\text{d}}T = - 0.014\) GPa/K) was found to be consistent with experimental data.     

Synthesis,TEM characterization and thermal behaviour of LiNiSi<Subscript>2</Subscript>O<Subscript>6</Subscript> pyroxene

A pyroxene with composition LiNiSi₂O₆ was synthesized at T = 1,473 K and P = 2.0 GPa; the cell parameters at T = 298 K are a = 9.4169(6) Å, b = 8.4465(7) Å, c = 5.2464(3) Å, β = 110.534(6)°, V = 390.78(3) Å³. TEM examination of the LiNiSi₂O₆ pyroxene showed the presence of h + k odd reflections indicative of a primitive lattice, and of antiphase domains obtained by dark field imaging of the h + k odd reflections. A HT in situ investigation was performed by examining TEM selected area diffraction patterns collected at high temperature and synchrotron radiation powder diffraction. In HTTEM the LiNiSi₂O₆ was examined together with LiCrSi₂O₆ pyroxene. In LiCrSi₂O₆ the h + k odd critical reflections disappear at about 340 K; they are sharp up to the transition temperature and do not change their shape until they disappear. In LiNiSi₂O₆ the h + k odd reflections are present up to sample deterioration at 650 K. A high temperature synchrotron radiation powder diffraction investigation was performed on LiNiSi₂O₆ between 298 and 773 K. The analysis of critical reflections and of changes in cell parameters shows that the space group is P-centred up to the highest temperature. The comparative analysis of the thermal and spontaneous strain contributions in P2₁/c and C2/c pyroxenes indicates that the high temperature strain in P-LiNiSi₂O₆ is very similar to that due to thermal strain only in C2/c spodumene and that a spontaneous strain contribution related to pre-transition features is not apparent in LiNiSi₂O₆. A different high-temperature behaviour in LiNiSi₂O₆ with respect to other pyroxenes is suggested, possibly in relation with the presence of Jahn–Teller distortion of the M1 polyhedron centred by low-spin Ni³⁺.     

Equation of state of synthetic qandilite Mg<Subscript>2</Subscript>TiO<Subscript>4</Subscript> at ambient temperature

Using a diamond-anvil cell and synchrotron X-ray diffraction, the compressional behavior of a synthetic qandilite Mg_2.00(1)Ti_1.00(1)O₄ has been investigated up to about 14.9 GPa at 300 K. The pressure–volume data fitted to the third-order Birch–Murnaghan equation of state yield an isothermal bulk modulus (K _T0) of 175(5) GPa, with its first derivative \(K_{T0}^{{\prime }}\) attaining 3.5(7). If \(K_{T0}^{{\prime }}\) is fixed as 4, the K _T0 value is 172(1) GPa. This value is substantially larger than the value of the adiabatic bulk modulus (K _S0) previously determined by an ultrasonic pulse echo method (152(7) GPa; Liebermann et al. in Geophys J Int 50:553–586, 1977), but in general agreement with the K _T0 empirically estimated on the basis of crystal chemical systematics (169 GPa; Hazen and Yang in Am Miner 84:1956–1960, 1999). Compared to the K _T0 values of the ulvöspinel (Fe₂TiO₄; ~148(4) GPa with \(K_{T0}^{{\prime }} = 4\)) and the ringwoodite solid solutions along the Mg₂SiO₄–Fe₂SiO₄ join, our finding suggests that the substitution of Mg²⁺ for Fe²⁺ on the T sites of the 4–2 spinels can have more significant effect on the K _T0 than that on the M sites.     

New insights on pressure,temperature, and chemical stability of CsAlSi<Subscript>5</Subscript>O<Subscript>12</Subscript>, a potential host for nuclear waste

A Cs-bearing polyphase aggregate with composition (in wt%): 76(1)_CsAlSi5O12 + 7(1)_CsAlSi2O6 + 17(1)_amorphous, was obtained from a clinoptilolite-rich epiclastic rock after a beneficiation process of the starting material (aimed to increase the fraction of zeolite to 90 wt%), cation exchange and then thermal treatment. CsAlSi₅O₁₂ is an open-framework compound with CAS topology; CsAlSi₂O₆ is a pollucite-like material with ANA topology. The thermal stability of this polyphase material was investigated by in situ high-T X-ray powder diffraction, the combined P–T effects by a series of runs with a single-stage piston cylinder apparatus, and its chemical stability following the “availability test” (“AVA test”) protocol. A series of additional investigations were performed by WDS–electron microprobe analysis in order to describe the P–T-induced modification of the material texture, and to chemically characterize the starting material and the run products. The “AVA tests” of the polyphase aggregate show an extremely modest release of Cs⁺: 0.05 mg/g. In response to applied temperature and at room P, CsAlSi₅O₁₂ experiences an unquenchable and displacive Ama2-to-Amam phase transition at about 770 K, and the Amam polymorph is stable in its crystalline form up to 1600 K; a crystalline-to-amorphous phase transition occurs between 1600 and 1650 K. In response to the applied P = 0.5 GPa, the crystalline-to-amorphous transition of CsAlSi₅O₁₂ occurs between 1670 and 1770 K. This leads to a positive Clapeyron slope (i.e., dP/dT > 0) of the crystalline-to-amorphous transition. When the polyphase aggregate is subjected at P = 0.5 GPa and T > 1770 K, CsAlSi₅O₁₂ melts and only CsAlSi₂O₆ (pollucite-like; dominant) and Cs-rich glass (subordinate) are observed in the quenched sample. Based on its thermo-elastic behavior, P–T phase stability fields, and Cs⁺ retention capacity, CsAlSi₅O₁₂ is a possible candidate for use in the immobilization of radioactive isotopes of Cs, or as potential solid hosts for ¹³⁷Cs γ-radiation source in sterilization applications. More in general, even the CsAlSi₅O₁₂-rich aggregate obtained by a clinoptilolite-rich epiclastic rock appears to be suitable for this type of utilizations.     

Compressibility of strontium orthophosphate Sr<Subscript>3</Subscript>(PO<Subscript>4</Subscript>)<Subscript>2</Subscript> at high pressure

High pressure in situ synchrotron X-ray diffraction experiment of strontium orthophosphate Sr₃(PO₄)₂ has been carried out to 20.0 GPa at room temperature using multianvil apparatus. Fitting a third-order Birch–Murnaghan equation of state to the P–V data yields a volume of V ₀ = 498.0 ± 0.1 Å³, an isothermal bulk modulus of K _T  = 89.5 ± 1.7 GPa, and first pressure derivative of K _T ′ = 6.57 ± 0.34. If K _T ′ is fixed at 4, K _T is obtained as 104.4 ± 1.2 GPa. Analysis of axial compressible modulus shows that the a-axis (K _a  = 79.6 ± 3.2 GPa) is more compressible than the c-axis (K _c  = 116.4 ± 4.3 GPa). Based on the high pressure Raman spectroscopic results, the mode Grüneisen parameters are determined and the average mode Grüneisen parameter of PO₄ vibrations of Sr₃(PO₄)₂ is calculated to be 0.30(2).     

Water diffusion in Mount Changbai peralkaline rhyolitic melt

Diffusion couple experiments with wet half (up to 4.6 wt%) and dry half were carried out at 789–1,516 K and 0.47–1.42 GPa to investigate water diffusion in a peralkaline rhyolitic melt with major oxide concentrations matching Mount Changbai rhyolite. Combining data from this work and a related study, total water diffusivity in peralkaline rhyolitic melt can be expressed as:

$ D_{{{\text{H}}_{ 2} {\text{O}}_{\text{t}} }} = D_{{{\text{H}}_{ 2} {\text{O}}_{\text{m}} }} \left( {1 - \frac{0.5 - X}{{\sqrt {[4\exp (3110/T - 1.876) - 1](X - X^{2} ) + 0.25} }}} \right), $

$ {\text{with}}\;D_{{{\text{H}}_{ 2} {\text{O}}_{\text{m}} }} = \exp \left[ { - 1 2. 7 8 9- \frac{13939}{T} - 1229.6\frac{P}{T} + ( - 27.867 + \frac{60559}{T})X} \right], $

where D is in m² s^?1, T is the temperature in K, P is the pressure in GPa, and X is the mole fraction of water and calculated as X = (C/18.015)/(C/18.015 + (100 ? C)/33.14), where C is water content in wt%. We recommend this equation in modeling bubble growth and volcanic eruption dynamics in peralkaline rhyolitic eruptions, such as the ~1,000-ad eruption of Mount Changbai in North East China. Water diffusivities in peralkaline and metaluminous rhyolitic melts are comparable within a factor of 2, in contrast with the 1.0–2.6 orders of magnitude difference in viscosities. The decoupling of diffusivity of neutral molecular species from melt viscosity, i.e., the deviation from the inversely proportional relationship predicted by the Stokes–Einstein equation, might be attributed to the small size of H₂O molecules. With distinct viscosities but similar diffusivity, bubble growth controlled by diffusion in peralkaline and metaluminous rhyolitic melts follows similar parabolic curves. However, at low confining pressure or low water content, viscosity plays a larger role and bubble growth rate in peralkaline rhyolitic melt is much faster than that in metaluminous rhyolite.

     

An experimental study of Ti and Zr partitioning among zircon, rutile, and granitic melt

In order to evaluate the effect of trace and minor elements (e.g., P, Y, and the REEs) on the high-temperature solubility of Ti in zircon (zrc), we conducted 31 experiments on a series of synthetic and natural granitic compositions [enriched in TiO₂ and ZrO₂; Al/(Na + K) molar ~1.2] at a pressure of 10 kbar and temperatures of ~1,400 to 1,200 °C. Thirty of the experiments produced zircon-saturated glasses, of which 22 are also saturated in rutile (rt). In seven experiments, quenched glasses coexist with quartz (qtz). SiO₂ contents of the quenched liquids range from 68.5 to 82.3 wt% (volatile free), and water concentrations are 0.4–7.0 wt%. TiO₂ contents of the rutile-saturated quenched melts are positively correlated with run temperature. Glass ZrO₂ concentrations (0.2–1.2 wt%; volatile free) also show a broad positive correlation with run temperature and, at a given T, are strongly correlated with the parameter (Na + K + 2Ca)/(Si·Al) (all in cation fractions). Mole fraction of ZrO₂ in rutile $ \left( {\mathop X\nolimits_{{{\text{ZrO}}_{ 2} }}^{\text{rt}} } \right) $ in the quartz-saturated runs coupled with other 10-kbar qtz-saturated experimental data from the literature (total temperature range of ~1,400 to 675 °C) yields the following temperature-dependent expression: $ {\text{ln}}\left( {\mathop X\nolimits_{{{\text{ZrO}}_{ 2} }}^{\text{rt}} } \right) + {\text{ln}}\left( {a_{{{\text{SiO}}_{2} }} } \right) = 2.638(149) - 9969(190)/T({\text{K}}) $ , where silica activity $ a_{{{\text{SiO}}_{2} }} $ in either the coexisting silica polymorph or a silica-undersaturated melt is referenced to α-quartz at the P and T of each experiment and the best-fit coefficients and their uncertainties (values in parentheses) reflect uncertainties in T and $ \mathop X\nolimits_{{{\text{ZrO}}_{2} }}^{\text{rt}} $ . NanoSIMS measurements of Ti in zircon overgrowths in the experiments yield values of ~100 to 800 ppm; Ti concentrations in zircon are positively correlated with temperature. Coupled with values for $ a_{{{\text{SiO}}_{2} }} $ and $ a_{{{\text{TiO}}_{2} }} $ for each experiment, zircon Ti concentrations (ppm) can be related to temperature over the range of ~1,400 to 1,200 °C by the expression: $ \ln \left( {\text{Ti ppm}} \right)^{\text{zrc}} + \ln \left( {a_{{{\text{SiO}}_{2} }} } \right) - \ln \left( {a_{{{\text{TiO}}_{2} }} } \right) = 13.84\left( {71} \right) - 12590\left( {1124} \right)/T\left( {\text{K}} \right) $ . After accounting for differences in $ a_{{{\text{SiO}}_{2} }} $ and $ a_{{{\text{TiO}}_{2} }} $ , Ti contents of zircon from experiments run with bulk compositions based on the natural granite overlap with the concentrations measured on zircon from experiments using the synthetic bulk compositions. Coupled with data from the literature, this suggests that at T ≥ 1,100 °C, natural levels of minor and trace elements in “granitic” melts do not appear to influence the solubility of Ti in zircon. Whether this is true at magmatic temperatures of crustal hydrous silica-rich liquids (e.g., 800–700 °C) remains to be demonstrated. Finally, measured $ D_{\text{Ti}}^{{{\text{zrc}}/{\text{melt}}}} $ values (calculated on a weight basis) from the experiments presented here are 0.007–0.01, relatively independent of temperature, and broadly consistent with values determined from natural zircon and silica-rich glass pairs.     

Experimental calibration of Forsterite–Anorthite–Ca-Tschermak–Enstatite (FACE) geobarometer for mantle peridotites

The crystallization of plagioclase-bearing assemblages in mantle rocks is witness of mantle exhumation at shallow depth. Previous experimental works on peridotites have found systematic compositional variations in coexisting minerals at decreasing pressure within the plagioclase stability field. In this experimental study we present new constraints on the stability of plagioclase as a function of different Na₂O/CaO bulk ratios, and we present a new geobarometer for mantle rocks. Experiments have been performed in a single-stage piston cylinder at 5–10 kbar, 1050–1150?°C at nominally anhydrous conditions using seeded gels of peridotite compositions (Na₂O/CaO?=?0.08–0.13; X _Cr = Cr/(Cr?+?Al)?=?0.07–0.10) as starting materials. As expected, the increase of the bulk Na₂O/CaO ratio extends the plagioclase stability to higher pressure; in the studied high-Na fertile lherzolite (HNa-FLZ), the plagioclase-spinel transition occurs at 1100?°C between 9 and 10 kbar; in a fertile lherzolite (FLZ) with Na₂O/CaO?=?0.08, it occurs between 8 and 9 kbar at 1100?°C. This study provides, together with previous experimental results, a consistent database, covering a wide range of P–T conditions (3–9 kbar, 1000–1150?°C) and variable bulk compositions to be used to define and calibrate a geobarometer for plagioclase-bearing mantle rocks. The pressure sensitive equilibrium:

$$\mathop {{\text{M}}{{\text{g}}_{\text{2}}}{\text{Si}}{{\text{O}}_{\text{4}}}^{{\text{Ol}}}}\limits_{{\text{Forsterite}}} +\mathop {{\text{CaA}}{{\text{l}}_{\text{2}}}{\text{S}}{{\text{i}}_{\text{2}}}{{\text{O}}_{\text{8}}}^{{\text{Pl}}}}\limits_{{\text{Anorthite}}~} =\mathop {{\text{CaA}}{{\text{l}}_{\text{2}}}{\text{Si}}{{\text{O}}_{\text{6}}}^{{\text{Cpx}}}}\limits_{{\text{Ca-Tschermak}}} +{\text{ }}\mathop {{\text{M}}{{\text{g}}_{\text{2}}}{\text{S}}{{\text{i}}_{\text{2}}}{{\text{O}}_{\text{6}}}^{{\text{Opx}}}}\limits_{{\text{Enstatite}}} ,$$

has been empirically calibrated by least squares regression analysis of experimental data combined with Monte Carlo simulation. The result of the fit gives the following equation:

$$P=7.2( \pm 2.9)+0.0078( \pm 0.0021)T{\text{ }}+0.0022( \pm 0.0001)T{\text{ }}\ln K,$$

$${R^2}=0.93,$$

where P is expressed in kbar and T in kelvin. K is the equilibrium constant K?=?a _CaTs × a _en/a _an × a _fo, where a _CaTs, a _en, a _an and a _fo are the activities of Ca-Tschermak in clinopyroxene, enstatite in orthopyroxene, anorthite in plagioclase and forsterite in olivine. The proposed geobarometer for plagioclase peridotites, coupled to detailed microstructural and mineral chemistry investigations, represents a valuable tool to track the exhumation of the lithospheric mantle at extensional environments.

     

Laboratory parallel-beam transmission X-ray powder diffraction investigation of the thermal behavior of nitratine NaNO<Subscript>3</Subscript>: spontaneous strain and structure evolution

Present work provides in-situ structural data at a fine temperature scale from RT to the melting point of nitratine, NaNO₃. From the analysis of log e ₃₃ versus log t plots, it is possible to prove that an univocal indication on the R \( \overline{3} \) c (low temperature, LT) → R \( \overline{3} \) m (high temperature, HT) transition mechanism cannot be obtained because of the relevant role played by the arbitrary assumptions required for defining the c ₀ dependence from temperature of the HT phase. This is due to the occurrence of excess thermal expansion for the HT phase. A significantly better fit for an Ising-spin structural model over a non-Ising rigid-body one has been obtained for the LT phase. Moreover, the Ising model led to a smooth variation of the oxygen site x fractional coordinate throughout the transition. The structure of the HT polymorph has been successfully refined considering an oxygen site at x, 0, ½, with 50% occupancy. Such model was the only acceptable one from the crystal chemical point of view as the alternative model (oxygen site at x, y, z with 25% occupancy) led to unrealistically aplanar \( {\text{NO}}_{3}^{ - } \) groups.     

H<Subscript>2</Subscript>O diffusion in peralkaline to peraluminous rhyolitic melts

The diffusion of water in a peralkaline and a peraluminous rhyolitic melt was investigated at temperatures of 714–1,493 K and pressures of 100 and 500 MPa. At temperatures below 923 K dehydration experiments were performed on glasses containing about 2 wt% H₂O _t in cold seal pressure vessels. At high temperatures diffusion couples of water-poor (<0.5 wt% H₂O _t ) and water-rich (~2 wt% H₂O _t ) melts were run in an internally heated gas pressure vessel. Argon was the pressure medium in both cases. Concentration profiles of hydrous species (OH groups and H₂O molecules) were measured along the diffusion direction using near-infrared (NIR) microspectroscopy. The bulk water diffusivity () was derived from profiles of total water () using a modified Boltzmann-Matano method as well as using fittings assuming a functional relationship between and Both methods consistently indicate that is proportional to in this range of water contents for both bulk compositions, in agreement with previous work on metaluminous rhyolite. The water diffusivity in the peraluminous melts agrees very well with data for metaluminous rhyolites implying that an excess of Al₂O₃ with respect to alkalis does not affect water diffusion. On the other hand, water diffusion is faster by roughly a factor of two in the peralkaline melt compared to the metaluminous melt. The following expression for the water diffusivity in the peralkaline rhyolite as a function of temperature and pressure was obtained by least-squares fitting:

Cu diffusivity in granitic melts with application to the formation of porphyry Cu deposits

We report new experimental data of Cu diffusivity in granite porphyry melts with 0.01 and 3.9 wt% H₂O at 0.15–1.0 GPa and 973–1523 K. A diffusion couple method was used for the nominally anhydrous granitic melt, whereas a Cu diffusion-in method using Pt₉₅Cu₅ as the source of Cu was applied to the hydrous granitic melt. The diffusion couple experiments also generate Cu diffusion-out profiles due to Cu loss to Pt capsule walls. Cu diffusivities were extracted from error function fits of the Cu concentration profiles measured by LA-ICP-MS. At 1 GPa, we obtain \({D_{{\text{Cu, dry, 1 GPa}}}}=\exp \left[ {( - {\text{13.89}} \pm {\text{0.42}}) - \frac{{{\text{12878}} \pm {\text{540}}}}{T}} \right],\) and \({D_{{\text{Cu, 3}}{\text{.9 wt\% }}{{\text{H}}_{\text{2}}}{\text{O}},{\text{ 1 GPa}}}}=\exp \left[ {( - 16.31 \pm 1.30) - \frac{{{\text{8148}} \pm {\text{1670}}}}{T}} \right],\) where D is Cu diffusivity in m²/s and T is temperature in K. The above expressions are in good agreement with a recent study on Cu diffusion in rhyolitic melt using the approach of Cu₂S dissolution. The observed pressure effect over 0.15–1.0 GPa can be described by an activation volume of 5.9 cm³/mol for Cu diffusion. Comparison of Cu diffusivity to alkali diffusivity and its variation with melt composition implies fourfold-coordinated Cu⁺ in silicate melts. Our experimental results indicate that in the formation of porphyry Cu deposits, the diffusive transport of magmatic Cu to sulfide liquids or fluid bubbles is highly efficient. The obtained Cu diffusivity data can also be used to assess whether equilibrium Cu partitioning can be reached within certain experimental durations.     

Elastic behavior of vanadinite,Pb<Subscript>10</Subscript>(VO<Subscript>4</Subscript>)<Subscript>6</Subscript>Cl<Subscript>2</Subscript>, a microporous non-zeolitic mineral

The high-pressure behavior of a vanadinite (Pb₁₀(VO₄)₆Cl₂, a = b = 10.3254(5), c = 7.3450(4) Å, space group P6₃/m), a natural microporous mineral, has been investigated using in-situ HP-synchrotron X-ray powder diffraction up to 7.67 GPa with a diamond anvil cell under hydrostatic conditions. No phase transition has been observed within the pressure range investigated. Axial and volume isothermal Equations of State (EoS) of vanadinite were determined. Fitting the P–V data with a third-order Birch-Murnaghan (BM) EoS, using the data weighted by the uncertainties in P and V, we obtained: V ₀ = 681(1) Å³, K ₀ = 41(5) GPa, and K′ = 12.5(2.5). The evolution of the lattice constants with P shows a strong anisotropic compression pattern. The axial bulk moduli were calculated with a third-order “linearized” BM-EoS. The EoS parameters are: a ₀ = 10.3302(2) Å, K ₀(a) = 35(2) GPa and K′(a) = 10(1) for the a-axis; c ₀ = 7.3520(3) Å, K ₀(c) = 98(4) GPa, and K′(c) = 9(2) for the c-axis (K ₀(a):K ₀(c) = 1:2.80). Axial and volume Eulerian-finite strain (fe) at different normalized stress (Fe) were calculated. The weighted linear regression through the data points yields the following intercept values: Fe _a (0) = 35(2) GPa for the a-axis, Fe _c (0) = 98(4) GPa for the c-axis and Fe _V (0) = 45(2) GPa for the unit-cell volume. The slope of the regression lines gives rise to K′ values of 10(1) for the a-axis, 9(2) for the c-axis and 11(1) for the unit cell-volume. A comparison between the HP-elastic response of vanadinite and the iso-structural apatite is carried out. The possible reasons of the elastic anisotropy are discussed.     

Fe<Superscript>3+</Superscript> reduction during biotite melting in graphitic metapelites: another origin of CO<Subscript>2</Subscript> in granulites

The Fe³⁺/Fe_tot of all Fe-bearing minerals has been analysed by Mössbauer spectroscopy in a suite of biotite-rich to biotite-free graphitic metapelite xenoliths, proxies of an amphibolite-granulite transition through progressive biotite melting. Biotite contains 9 to 16% Fe³⁺/Fe_tot, whereas garnet, cordierite and ilmenite are virtually Fe³⁺ -free (0–1% Fe³⁺/Fe_tot) in all samples, regardless of biotite presence. Under relatively reducing conditions (graphite-bearing assemblages), biotite is the only carrier of Fe³⁺ during high-temperature metamorphism; therefore, its disappearance by melting represents an important event of iron reduction during granulite formation, because haplogranitic melts usually incorporate small amounts of ferric iron. Iron reduction is accompanied by the oxidation of carbon and the production of CO₂, according to the redox reaction:

Pressure induced phase transition in Pb<Subscript>6</Subscript>Bi<Subscript>2</Subscript>S<Subscript>9</Subscript>

The crystal structure of Pb₆Bi₂S₉ is investigated at pressures between 0 and 5.6 GPa with X-ray diffraction on single-crystals. The pressure is applied using diamond anvil cells. Heyrovskyite (Bbmm, a = 13.719(4) Å, b = 31.393(9) Å, c = 4.1319(10) Å, Z = 4) is the stable phase of Pb₆Bi₂S₉ at ambient conditions and is built from distorted moduli of PbS-archetype structure with a low stereochemical activity of the Pb²⁺ and Bi³⁺ lone electron pairs. Heyrovskyite is stable until at least 3.9 GPa and a first-order phase transition occurs between 3.9 and 4.8 GPa. A single-crystal is retained after the reversible phase transition despite an anisotropic contraction of the unit cell and a volume decrease of 4.2%. The crystal structure of the high pressure phase, β-Pb₆Bi₂S₉, is solved in Pna2 ₁ (a = 25.302(7) Å, b = 30.819(9) Å, c = 4.0640(13) Å, Z = 8) from synchrotron data at 5.06 GPa. This structure consists of two types of moduli with SnS/TlI-archetype structure in which the Pb and Bi lone pairs are strongly expressed. The mechanism of the phase transition is described in detail and the results are compared to the closely related phase transition in Pb₃Bi₂S₆ (lillianite).     

High-pressure phase transitions of CaRhO<Subscript>3</Subscript> perovskite

High-pressure phase transitions of CaRhO₃ perovskite were examined at pressures of 6–27 GPa and temperatures of 1,000–1,930°C, using a multi-anvil apparatus. The results indicate that CaRhO₃ perovskite successively transforms to two new high-pressure phases with increasing pressure. Rietveld analysis of powder X-ray diffraction data indicated that, in the two new phases, the phase stable at higher pressure possesses the CaIrO₃-type post-perovskite structure (space group Cmcm) with lattice parameters: a = 3.1013(1) Å, b = 9.8555(2) Å, c = 7.2643(1) Å, V _m  = 33.43(1) cm³/mol. The Rietveld analysis also indicated that CaRhO₃ perovskite has the GdFeO₃-type structure (space group Pnma) with lattice parameters: a = 5.5631(1) Å, b = 7.6308(1) Å, c = 5.3267(1) Å, V _m  = 34.04(1) cm³/mol. The third phase stable in the intermediate P, T conditions between perovskite and post-perovskite has monoclinic symmetry with the cell parameters: a = 12.490(3) Å, b = 3.1233(3) Å, c = 8.8630(7) Å, β = 103.96(1)°, V _m  = 33.66(1) cm³/mol (Z = 6). Molar volume changes from perovskite to the intermediate phase and from the intermediate phase to post-perovskite are –1.1 and –0.7%, respectively. The equilibrium phase relations determined indicate that the boundary slopes are large positive values: 29 ± 2 MPa/K for the perovskite—intermediate phase transition and 62 ± 6 MPa/K for the intermediate phase—post-perovskite transition. The structural features of the CaRhO₃ intermediate phase suggest that the phase has edge-sharing RhO₆ octahedra and may have an intermediate structure between perovskite and post-perovskite.