H-bonding scheme and cation partitioning in axinite: a single-crystal neutron diffraction and Mössbauer spectroscopic study

The crystal chemistry of a ferroaxinite from Colebrook Hill, Rosebery district, Tasmania, Australia, was investigated by electron microprobe analysis in wavelength-dispersive mode, inductively coupled plasma–atomic emission spectroscopy (ICP–AES), ⁵⁷Fe Mössbauer spectroscopy and single-crystal neutron diffraction at 293 K. The chemical formula obtained on the basis of the ICP–AES data is the following: \( ^{X1,X2} {\text{Ca}}_{4.03} \,^{Y} \left( {{\text{Mn}}_{0.42} {\text{Mg}}_{0.23} {\text{Fe}}^{2 + }_{1.39} } \right)_{\varSigma 2.04} \,^{Z1,Z2} \left( {{\text{Fe}}^{3 + }_{0.15} {\text{Al}}_{3.55} {\text{Ti}}_{0.12} } \right)_{\varSigma 3.82} \,^{T1,T2,T3,T4} \left( {{\text{Ti}}_{0.03} {\text{Si}}_{7.97} } \right)_{\varSigma 8} \,^{T5} {\text{B}}_{1.96} {\text{O}}_{30} \left( {\text{OH}} \right)_{2.18} \). The ⁵⁷Fe Mössbauer spectrum shows unambiguously the occurrence of Fe²⁺ and Fe³⁺ in octahedral coordination only, with Fe²⁺/Fe³⁺ = 9:1. The neutron structure refinement provides a structure model in general agreement with the previous experimental findings: the tetrahedral T1, T2, T3 and T4 sites are fully occupied by Si, whereas the T5 site is fully occupied by B, with no evidence of Si at the T5, or Al or Fe³⁺ at the T1–T5 sites. The structural and chemical data of this study suggest that the amount of B in ferroaxinite is that expected from the ideal stoichiometry: 2 a.p.f.u. (for 32 O). The atomic distribution among the X1, X2, Y, Z1 and Z2 sites obtained by neutron structure refinement is in good agreement with that based on the ICP–AES data. For the first time, an unambiguous localization of the H site is obtained, which forms a hydroxyl group with the oxygen atom at the O16 site as donor. The H-bonding scheme in axinite structure is now fully described: the O16–H distance (corrected for riding motion effect) is 0.991(1) Å and an asymmetric bifurcated bonding configuration occurs, with O5 and O13 as acceptors [i.e. with O16···O5 = 3.096(1) Å, H···O5 = 2.450(1) Å and O16–H···O5 = 123.9(1)°; O16···O13 = 2.777(1) Å, H···O13 = 1.914(1) Å and O16–H···O13 = 146.9(1)°].     

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.

     

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.

     

Earthquake hazard and risk assessment based on unified scaling law for earthquakes: Altai–Sayan Region

We apply the general concept of seismic risk analysis based on morphostructural analysis of the territory, pattern recognition of earthquake-prone nodes, and the Unified Scaling Law for Earthquakes, USLE, in another seismic region of Russia to the west from Lake Baikal, i.e., Altai–Sayan Region. The USLE generalizes the empirical Gutenberg–Richter relationship making use of apparently fractal distribution of earthquake sources of different size: \( \log_{10} N\left( {M,L} \right)\, = \,A\, + \,B \cdot \left( {5\, - \,M} \right)\, + \,C \cdot \log_{10} L, \) where N (M, L) is the expected annual number of earthquakes of a certain magnitude M within an seismically prone area of linear dimension L. The local estimates of A, B, and C allow determination of the expected maximum credible magnitude in a given time interval and the associated spread around ground shaking parameters (e.g., peak ground acceleration, PGA, or macroseismic intensity, I₀). Compilation of the corresponding seismic hazard map of Altai–Sayan Region and its rigorous testing against the available seismic evidences in the past is used to model regional maps of specific earthquake risks for population, cities, and infrastructures.     

The influence of materials on the breaching process of natural dams

In this study, a series of natural dam overtopping laboratory tests are reported. In these tests, the effect of seven different sediment mixtures on the breaching process was investigated. According to the test results, three stages of the breaching process of natural dams made of different materials were observed. Backward erosion was the primary cause for the incising slopes. The effects of backward erosion became stronger with the larger fines contents of the materials. With an increase in the median diameter (d ₅₀) of particles, the breaching time became longer. However, the peak discharge became smaller. With an increase in the fines contents (p), the median diameter of the particles and the void ratio were changed, which resulted in a decrease in the breaching time and an increase in the peak discharge. The breaching time and peak discharge were more sensitive to the median diameter than to the fines contents. The relation between breach width and depth was found to follow a logistic function \( W\kern0.5em =\kern0.5em \frac{\zeta }{1\kern0.5em +\kern0.5em {e}^{\left(-k\left(D\kern0.5em -\kern0.5em {D}_0\right)\right)}} \). The parameters ζ, k, and D ₀ are defined by a linear relationship with the median diameter and fines content. A breach of the side slope occurred as a tensile failure when the fines contents of the materials were large; otherwise, shear failure occurred. Furthermore, when the materials had fewer fines contents, the volume of the collapsed breach side slope became larger.     

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.     

Trend analysis of evapotranspiration over India: Observed from long-term satellite measurements

Owing to the lack of consistent spatial time series data on actual evapotranspiration (ET), very few studies have been conducted on the long-term trend and variability in ET at a national scale over the Indian subcontinent. The present study uses biome specific ET data derived from NOAA satellite’s advanced very high resolution radiometer to investigate the trends and variability in ET over India from 1983 to 2006. Trend analysis using the non-parametric Mann–Kendall test showed that the domain average ET decreased during the period at a rate of \(0.22\,\hbox {mm year}^{-1}\). A strong decreasing trend (\(m = -1.75\, \hbox {mm year}^{-1}\), \(F = 17.41\), \(P\) 0.01) was observed in forest regions. Seasonal analyses indicated a decreasing trend during southwest summer monsoon (\(m= -0.320\, \hbox {mm season}^{-1}\,\hbox {year}^{-1})\) and post-monsoon period (\(m= -0.188\, \hbox {mm season}^{-1 }\,\hbox {year}^{-1})\). In contrast, an increasing trend was observed during northeast winter monsoon (\(m = 0.156 \,\hbox {mm season}^{-1 }\,\hbox {year}^{-1})\) and pre-monsoon (\(m = 0.068\, \hbox {mm season}^{-1 }\,\hbox {year}^{-1})\) periods. Despite an overall net decline in the country, a considerable increase ( \(4 \,\hbox {mm year}^{-1}\)) was observed over arid and semi-arid regions. Grid level correlation with various climatic parameters exhibited a strong positive correlation (\(r \!>\!0.5\)) of ET with soil moisture and precipitation over semi-arid and arid regions, whereas a negative correlation (\(r\) \(-0.5\)) occurred with temperature and insolation in dry regions of western India. The results of this analysis are useful for understanding regional ET dynamics and its relationship with various climatic parameters over India. Future studies on the effects of ET changes on the hydrological cycle, carbon cycle, and energy partitioning are needed to account for the feedbacks to the climate.     

Fabric evolution of granular materials along imposed stress paths

The stress–strain behavior of a granular material is dominated by its internal structure, which is related to the spatial connectivity of particles, and the force chain network. In this study, a series of discrete element simulations were carried out to investigate the evolution of internal structure and force chain networks in initially isotropic granular materials along various imposed stress paths. The fabric tensor of the strong sub-network, which is the bearing network toward loading, can be related to the applied stresses uniquely. The principal directions of fabric tensor of the strong sub-network coincide with those of stress tensor during the loading process in the Lode coordinate system. The fabric of the whole contact network in the pre- and post-peak deformation stages can be related to the applied stresses as \(q_{\phi } = B\left( {q/p} \right)^{z}\) (B and z are constants depending on loading condition, such as the stress paths and mean stress level) and \(\phi_{1} :\phi_{2} :\phi_{3} \approx \left( {\sigma_{1} } \right)^{0.4} :\left( {\sigma_{2} } \right)^{0.4} :\left( {\sigma_{3} } \right)^{0.4}\), respectively. At the critical stress state, the deviator of fabric tensor of the strong sub-network is much larger than that of the whole contact network. When plotted on the π-plane, the fabric state of the strong sub-network can be expressed as a Lade’s surface, while the fabric state of the whole network corresponds to an inverted Lade’s surface.     

The potential use of OH-defects in enstatite as geobarometer

In this study, single crystals of pure enstatite (Mg₂Si₂O₆) were synthesised under water-saturated conditions at 4 and 8 GPa and 1,150°C with variable silica activity, leading to phase assemblages enstatite + forsterite, enstatite or enstatite + coesite. Run products were investigated using an FTIR spectrometer equipped with a focal plane array detector enabling IR imaging with a lateral pixel resolution of 2.7 μm. IR spectra within the OH-absorption region show two different groups of absorption bands: group 1 (wavenumbers at 3,592 and 3,687 cm^?1) shows strongest absorptions for E||n _β, whereas group 2 (wavenumbers at 3,067 and 3,362 cm^?1) shows strongest absorptions for E||n _γ. The groups are related to different defect types, group 1 to tetrahedral defects (T-site vacancies) and group 2 to octahedral defects (M-site vacancies). The intensity ratio of the bands within one group (i.e. A ₃₀₆₇/A ₃₃₆₂ and A ₃₅₉₂/A ₃₆₈₇) and the intensity ratio of E||n _γ and E||n _α in group 2 bands remain constant within error. In contrast, the intensity ratio of group 2 to group 1 absorption bands [e.g. (A ₃₃₆₂)/(A ₃₆₈₇)] is sensitive to the SiO₂ activity and pressure. On the basis of the results of this and previous studies, a barometer for pure orthoenstatite coexisting with forsterite can be formulated:\( P\,[{\text{GPa}}] = 1.056 + \sqrt {{\frac{{1.025 - A_{{\left( {3362} \right)/\left[ {(3362) + (3687)} \right]}} }}{0.009}}} , \) where A ₍₃₃₆₂₎ and A ₍₃₆₈₇₎ are the integral absorbances of the component E||n _γ of the absorption bands at 3,362 cm^?1 and the component E||n _β of the absorption band at 3,687 cm^?1, respectively.     

Theoretical calculation of equilibrium Mg isotope fractionation between silicate melt and its vapor

Isotope fractionation during the evaporation of silicate melt and condensation of vapor has been widely used to explain various isotope signals observed in lunar soils, cosmic spherules, calcium–aluminum-rich inclusions, and bulk compositions of planetary materials. During evaporation and condensation, the equilibrium isotope fractionation factor (α) between high-temperature silicate melt and vapor is a fundamental parameter that can constrain the melt’s isotopic compositions. However, equilibrium α is difficult to calibrate experimentally. Here we used Mg as an example and calculated equilibrium Mg isotope fractionation in MgSiO₃ and Mg₂SiO₄ melt–vapor systems based on first-principles molecular dynamics and the high-temperature approximation of the Bigeleisen–Mayer equation. We found that, at 2500 K, δ²⁵Mg values in the MgSiO₃ and Mg₂SiO₄ melts were 0.141?±?0.004 and 0.143?±?0.003‰ more positive than in their respective vapors. The corresponding δ²⁶Mg values were 0.270?±?0.008 and 0.274?±?0.006‰ more positive than in vapors, respectively. The general \(\alpha - T\) equations describing the equilibrium Mg α in MgSiO₃ and Mg₂SiO₄ melt–vapor systems were: \(\alpha_{{{\text{Mg}}\left( {\text{l}} \right) - {\text{Mg}}\left( {\text{g}} \right)}} = 1 + \frac{{5.264 \times 10^{5} }}{{T^{2} }}\left( {\frac{1}{m} - \frac{1}{{m^{\prime}}}} \right)\) and \(\alpha_{{{\text{Mg}}\left( {\text{l}} \right) - {\text{Mg}}\left( {\text{g}} \right)}} = 1 + \frac{{5.340 \times 10^{5} }}{{T^{2} }}\left( {\frac{1}{m} - \frac{1}{{m^{\prime}}}} \right)\), respectively, where m is the mass of light isotope ²⁴Mg and m′ is the mass of the heavier isotope, ²⁵Mg or ²⁶Mg. These results offer a necessary parameter for mechanistic understanding of Mg isotope fractionation during evaporation and condensation that commonly occurs during the early stages of planetary formation and evolution.     

Soil erosion assessment on hillslope of GCE using RUSLE model

A new method for obtaining the C factor (i.e., vegetation cover and management factor) of the RUSLE model is proposed. The method focuses on the derivation of the C factor based on the vegetation density to obtain a more reliable erosion prediction. Soil erosion that occurs on the hillslope along the highway is one of the major problems in Malaysia, which is exposed to a relatively high amount of annual rainfall due to the two different monsoon seasons. As vegetation cover is one of the important factors in the RUSLE model, a new method that accounts for a vegetation density is proposed in this study. A hillslope near the Guthrie Corridor Expressway (GCE), Malaysia, is chosen as an experimental site whereby eight square plots with the size of \(8\times 8\) and \(5\times 5\) m are set up. A vegetation density available on these plots is measured by analyzing the taken image followed by linking the C factor with the measured vegetation density using several established formulas. Finally, erosion prediction is computed based on the RUSLE model in the Geographical Information System (GIS) platform. The C factor obtained by the proposed method is compared with that of the soil erosion guideline Malaysia, thereby predicted erosion is determined by both the C values. Result shows that the C value from the proposed method varies from 0.0162 to 0.125, which is lower compared to the C value from the soil erosion guideline, i.e., 0.8. Meanwhile predicted erosion computed from the proposed C value is between 0.410 and \(3.925\, \hbox {t ha}^{-1 }\,\hbox {yr}^{-1}\) compared to 9.367 to \(34.496\, \hbox {t ha}^{-1}\,\hbox {yr}^{-1 }\) range based on the C value of 0.8. It can be concluded that the proposed method of obtaining a reasonable C value is acceptable as the computed predicted erosion is found to be classified as a very low zone, i.e. less than \(10\, \hbox {t ha}^{-1 }\,\hbox {yr}^{-1}\) whereas the predicted erosion based on the guideline has classified the study area as a low zone of erosion, i.e., between 10 and \(50\, \hbox {t ha}^{-1 }\,\hbox {yr}^{-1}\).     

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 β.     

Thermodynamics and phonon dispersion of pyrope and grossular silicate garnets from ab initio simulations

The phonon dispersion and thermodynamic properties of pyrope (\(\hbox {Mg}_3\hbox {Al}_2\hbox {Si}_3\hbox {O}_{12}\)) and grossular (\(\hbox {Ca}_3\hbox {Al}_2\hbox {Si}_3\hbox {O}_{12}\) ) have been computed by using an ab initio quantum mechanical approach, an all-electron variational Gaussian-type basis set and the B3LYP hybrid functional, as implemented in the Crystal program. Dispersion effects in the phonon bands have been simulated by using supercells of increasing size, containing 80, 160, 320, 640, 1280 and 2160 atoms, corresponding to 1, 2, 4, 8, 16 and 27 \(\mathbf {k}\) points in the first Brillouin zone. Phonon band structures, density of states and corresponding inelastic neutron scattering spectra are reported. Full convergence of the various thermodynamic properties, in particular entropy (S) and specific heat at constant volume (\(C_\mathrm{{V}}\)), with the number of \(\mathbf {k}\) points is achieved with 27 \(\mathbf {k}\) points. The very regular behavior of the S(T) and \(C_\mathrm{{V}}(T)\) curves as a function of the number of \(\mathbf {k}\) points, determined by high numerical stability of the code, permits extrapolation to an infinite number of \(\mathbf {k}\) points. The limiting value differs from the 27-\(\mathbf {k}\) case by only 0.40 % at 100 K for S (the difference decreasing to 0.11 % at 1000 K) and by 0.29 % (0.05 % at 1000 K) for \(C_\mathrm{{V}}\). The agreement with the experimental data is rather satisfactory. We also address the problem of the relative entropy of pyrope and grossular, a still debated question. Our lattice dynamical calculations correctly describe the larger entropy of pyrope than grossular by taking into account merely vibrational contributions and without invoking “static disorder” of the Mg ions in dodecahedral sites. However, as the computed entropy difference is found to be larger than the experimental one by a factor of 2–3, present calculations cannot exclude possible thermally induced structural changes, which could lead to further conformational contributions to the entropy.     

The concept of geochemical potential field exemplified by ore formation field

Geochemical potential field is defined as the scope within the earth’s space where a given component in a certain phase of a certain material system is acted upon by a diffusion force, depending on its spatial coordinatesX, Y andZ. The three coordinates follow the relations: $$NF_{ix} = - \frac{{\partial \mu }}{{\partial x}}, NF_{iy} = - \frac{{\partial \mu }}{{\partial y}}, NF_{iz} = - \frac{{\partial \mu }}{{\partial z}}$$ The characteristics of such a field can be summarized as: (1) The summation of geochemical potentials related to the coordinatesX, Y, Z, or pseudo-velocity head, pseudo-pressure head and pseudo-potential head of a certain component in the earth is a constant as given by $$\mu _x + \mu _y + \mu _z = c$$ or $$\mu _{x2} + \mu _{y2} + \mu _{z2} = \mu _{x1} + \mu _{y1} + \mu _{z1} $$ Derived from these relations is the principle of geochemical potential conservation. The following relations have the same physical significance: $$\mu _k + \mu _u + \mu _p = c$$ or $$\mu _{k2} + \mu _{u2} + \mu _{p2} = \mu _{k1} + \mu _{u1} + \mu _{p1} $$ (2) Geochemical potential field is a vector field quantified by geochemical field intensity which is defined as the diffusion force applied to one molecular volume (or one atomic volume) of a certain component moving from its higher concentration phase to lower concentration phase. The geochemical potential field intensity is given by $$\begin{gathered} E = - grad\mu \hfill \\ E = \frac{{RT}}{x}i + \frac{{RT}}{y}j + \frac{{RT}}{z}K \hfill \\ \end{gathered} $$ The present theory has been inferred to interpret the mechanism of formation of some tungsten ore deposits in China.     

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 effect of alkalis and polymerization on the solubility of H<Subscript>2</Subscript>O and CO<Subscript>2</Subscript> in alkali-rich silicate melts

The effect of alkalis on the solubility of H₂O and CO₂ in alkali-rich silicate melts was investigated at 500 MPa and 1,250 °C in the systems with H₂O/(H₂O + CO₂) ratio varying from 0 to 1. Using a synthetic analog of phonotephritic magma from Alban Hills (AH1) as a base composition, the Na/(Na + K) ratio was varied from 0.28 (AH1) to 0.60 (AH2) and 0.85 (AH3) at roughly constant total alkali content. The obtained results were compared with the data for shoshonitic and latitic melts having similar total alkali content but different structural characteristics, e.g., NBO/T parameter (the ratio of non-bridging oxygens over tetrahedrally coordinated cations), as those of the AH compositions. Little variation was observed in H₂O solubility (melt equilibrated with pure H₂O fluid) for the whole compositional range in this study with values ranging between 9.7 and 10.2 wt. As previously shown, the maximum CO₂ content in melts equilibrated with CO₂-rich fluids increases strongly with the NBO/T from 0.29 wt % for latite (NBO/T = 0.17) to 0.45 wt % for shoshonite (NBO/T = 0.38) to 0.90 wt % for AH2 (NBO/T = 0.55). The highest CO₂ contents determined for AH3 and AH1 are 1.18 ± 0.05 wt % and 0.86 ± 0.12 wt %, respectively, indicating that Na is promoting carbonate incorporation stronger than potassium. At near constant NBO/T, CO₂ solubility increases from 0.86 ± 0.12 wt % in AH1 [Na/(Na + K)] = 0.28, to 1.18 ± 0.05 wt % in AH3 [Na/(Na + K)] = 0.85, suggesting that Na favors CO₂ solubility on an equimolar basis. An empirical equation is proposed to predict the maximum CO₂ solubility at 500 MPa and 1,100–1,300 °C in various silicate melts as a function of the NBO/T, (Na + K)/∑cations and Na/(Na + K) parameters: \({\text{wt}}\% \;{\text{CO}}_{2} = - 0.246 + 0.014\exp \left( {6.995 \cdot \frac{\text{NBO}}{T}} \right) + 3.150 \cdot \frac{{{\text{Na}} + {\text{K}}}}{{\varSigma {\text{cations}}}} + 0.222 \cdot \frac{\text{Na}}{{{\text{Na}} + {\text{K}}}}.\) This model is valid for melt compositions with NBO/T between 0.0 and 0.6, (Na + K)/∑cation between 0.08 and 0.36 and Na/(Na + K) ratio from 0.25 to 0.95 at oxygen fugacities around the quartz–fayalite–magnetite buffer and above.