Trace element partitioning and melt structure: An experimental study at 1 atm pressure

Diopside-melt and forsterite-melt rare earth (REE) and Ni partition coefficients have been determined as a function of bulk compositions of the melt. Available Raman spectroscopic data have been used to determine the structures of the melts coexisting with diopside and forsterite. The compositional dependence of the partition coefficients is then related to the structural changes of the melt.The melts in all experiments have a ratio of nonbridging oxygens to tetrahedral cations ($\text{NBO}\text{T}$) between 1 and 0. The quenched melts consist of structural units that have, on the average, 2 (chain), 1 (sheet) and 0 (three-dimensional network) nonbridging oxygens per tetrahedral cation. The proportions of these structural units in the melts, as well as the overall $\text{NBO}\text{T}$, change as a function of the bulk composition of the melt.It has been found that Ce, Sm, Tm and Ni crystal-liquid partition coefficients ($\text{K}{}^{\mathrm{crystal?liq}}{}_{\mathrm{i}}\text{=}\text{C}{}^{\mathrm{crystal}}{}_{\mathrm{i}}\text{C}{}^{\mathrm{liq}}{}_{\mathrm{i}}$) decrease linearly with increasing $\text{NBO}\text{T}$. The values of the individual REE crystal-liquid trace element partition coefficients have different functional relations to $\text{NBO}\text{T}$, so that the degree of light REE enrichment of the melts would depend on their $\text{NBO}\text{T}$.The solution mechanisms of minor oxides such as CO₂, H₂O, TiO₂, P₂O₅ and Fe₂O₃ in silicate melts are known. These data have been recast as changes of $\text{NBO}\text{T}$ of the melts with regard to the type of oxide and its concentration in the melt. From such data the dependence of crystal-liquid partition coefficients on concentration and type of minor oxide in melt solution has been calculated.     

Density calculations for silicate liquids. I. Revised method for aluminosilicate compositions

Equations are developed for calculating the density of aluminosilicate liquids as a function of composition and temperature. The mean molar volume at reference temperature T_r, is given by $\text{V}{}_{\mathrm{r}}\text{= ∑X}{}_{\mathrm{i}}\text{V}\text{?}{}^{\mathrm{o}}{}_{\mathrm{i}}\text{+ X}{}_{\mathrm{A}}\text{V}\text{?}{}^{\mathrm{o}}{}_{\mathrm{A}}$, where the summation is taken over all oxide components except A1₂O₃, X stands for mole fraction, terms are constants derived independently from an analysis of volume-composition relations in alumina-free silicate liquids, and $\text{V}\text{?}{}^{\mathrm{o}}{}_{\mathrm{A}}$ is the composition-dependent apparent partial molar volume of Al₂O₃. The thermal expansion coefficient of aluminosilicate liquids is given by $\text{α = ∑X}{}_{\mathrm{i}}\text{}\text{?}\text{ga}{}_{\mathrm{i}}{}^{\mathrm{o}}\text{+ X}{}_{\mathrm{A}}\text{}\text{?}\text{ga}{}_{\mathrm{A}}{}^{\mathrm{o}}$, where $\text{}\text{?}\text{ga}{}_{\mathrm{i}}{}^{\mathrm{o}}$ terms are constants independent of temperature and composition, and $\text{}\text{?}\text{ga}{}^{\mathrm{o}}{}_{\mathrm{A}}$ is a composition-dependent term representing the effect of Al₂O₃ on the thermal expansion. Parameters necessary to calculate the volume of silicate liquids at any temperature T according to V(T) = V_rexp[α(T-T_r)], where T_r = 1400°C have been evaluated by least-square analysis of selected density measurements in aluminosilicate melts. Mean molar volumes of aluminosilicate liquids calculated according to the model equation conform to experimentally measured volumes with a root mean square difference of 0.28 $\text{cc}\text{mole}$ and an average absolute difference of 0.90% for 248 experimental observations. The compositional dependence of $\text{V}\text{?}{}^{\mathrm{o}}{}_{\mathrm{A}}$ is discussed in terms of several possible interpretations of the structural role of Al³⁺ in aluminosilicate melts.     

Diffusivity of oxygen in jadeite and diopside melts at high pressures

The diffusivity of oxygen was determined in melts of Jadeite (NaAlSi₂O₆) and diopside (CaMgSi₂O₆) compositions using diffusion couples with ¹⁸O as a tracer. In the Jadeite melt, the diffusivity of oxygen increases from $\text{6.87}{}_{\mathrm{?0.25}}{}^{\mathrm{+0.28}}\text{× 10}{}^{\mathrm{?10}}\text{cm}{}^2\text{/}\text{sec}$ at 5 Kb to $\text{1.32 ± 0.08 × 10}{}^{\mathrm{?9}}\text{cm}{}^2\text{/}\text{sec}$ at 20 Kb at constant temperature (1400°C), whereas in the diopside melt at 1650°C, the diffusivity decreases from $\text{7.30}{}_{\mathrm{?0.18}}{}^{0.29}\text{× 10}{}^{\mathrm{?7}}\text{cm}{}^2\text{/}\text{sec}$ at 10 Kb to $\text{5.28}{}_{\mathrm{?0.55}}{}^{\mathrm{+0.60}}\text{× 10}{}^{\mathrm{?7}}\text{cm}{}^2\text{/}\text{sec}$ at 17 Kb. These results demonstrate that the diffusivity is inversely correlated with the viscosity of the melt. For the jadeite melt, in particular, the inverse correlation is very well approximated by the Eyring equation using the diameter of oxygen ions as a unit distance of translation, suggesting that the viscous flow is rate-limited by the diffusion of individual oxygen ions. In the diopside melt, the activation volume is slightly greater than the molar volume of oxygen ion, indicating that the individual oxygen ion is the diffusion unit. The negative activation volume obtained for the jadeite melt is interpreted as the volume decrease associated with a diffusive jump of an oxygen ion due to local collapse of the network structure.     

PbO-SiO2 melts: structure and thermodynamics of mixing

Thermodynamic properties of PbO-SiO₂ melts, obtained from published data and calculated from freezing point depressions, reflect the gradual polymerization of silicate anions in the melt as the $\text{SiO}{}_2\text{PbO}$ ratio is increased. The free energy of mixing curve at 1000°C has a minimum at 40 mole % SiO₂ and is convex-upward between 72 and 98 mole % SiO₂. The latter is an indication of metastable liquid immiscibility. The free energy minimum is correlated with the maximum in the distribution of nonbridging oxygens in the melt. In SiO₂-poor melts, the activities of PbO and SiO₂ (pure liquid standard states) show sharp negative deviations from ideality. The PbO activity reflects the paucity of free oxygen species in the melt whereas the SiO₂ activity reflects the depolymerized state of the silicate anions. In more SiO₂-rich melts, the activity of SiO₂ shows a positive deviation from ideality which is qualitatively correlated to a polymerization parameter. The heat of mixing term has a minimum of ?2000 cal at 35 mole % SiO₂ and a maximum of +200 cal at 90 mole % SiO₂. The minimum is associated with the exothermic heat effect obtained during the reaction (O⁰) + (O^2?) = 2(O^?), whereas the maximum corresponds to the endothermic heat effect obtained when coordination polyhedra of oxygens form around the Pb cation. The entropy of mixing curve has the same form but is systematically smaller than a theoretical curve calculated on the assumption of random mixing of oxygen species. The discrepancy is due to the entropy loss obtained by the clustering of oxygen species to form complex silicate species.     

Plagioclase—melt equilibria

The crystallization of plagioclase feldspar from magmatic liquid has been investigated experimentally under equilibrium conditions at 1 atm total pressure in the temperature range 1400-1095°C. Natural and synthetic melts of composition basalt to rhyolite were used, crystallizing plagioclase of composition An₈₉-An₃₂.The experimental results are analyzed initially in terms of elemental plagioclase/melt partition coefficients (D). D_Si is always less than unity and is invariant with temperature. D_A1 is always greater than unity and is relatively insensitive to temperature. D_Na is less than unity above 1200°C and is strongly dependent upon temperature. D_Ca is greater than unity below 1430°C and is strongly dependent upon temperature.Analysis of the temperature-dependence of equilibrium constants for plagioclase-melt formation and exchange reactions in which several mixing models for the melt are considered, leads to the conclusion that, with appropriate choice of melt-components, the melt-components mix quasi-ideally. At fixed temperature in the absence of H₂O, the equilibrium constant for the equilibrium of albite with the melt is insensitive to changes in melt-composition, and is insensitive to changes in pressure up to at least 10 kbars. As a consequence the composition of plagioclase crystallizing at known temperature and at low total pressure from a dry melt of known composition may be predicted []. However, the equilibrium constant is sensitive to changes in water pressure.The analysis further suggests that Na is intimately associated with tetrahedrally-coordinated Al in the melt, while Ca appears to be partitioned between at least two distinct melt-sites.     

Distribution coefficients of Eu and Sr for plagioclase-liquid and clinopyroxene-liquid equilibria in oceanic ridge basalt: an experimental study

The distribution coefficients of Eu and Sr for plagioclase-liquid and clinopyroxene-liquid pairs as a function of temperature and oxygen fugacity were experimentally investigated using an oceanic ridge basalt enriched with Eu and Sr as the starting material. Experiments were conducted between 1190° and 1140°C over a range of oxygen fugacities between 10^?8 and 10^?14 atm.The molar distribution coefficients are given by the equations: log^PL_Sr = 7320/T ? 4.62logK^CPX_Sr = 18020/T ? 13.10. Similarly, the weight fraction distribution coefficients are given by the equations: logD^PL_Sr = 6570/T ? 4.30logD^CPX_Sr = 18434/T ? 13.62.Although the mole fraction distribution coefficients have a smaller dependence on bulk composition than do the weight fraction distribution coefficients, they are not independent of bulk composition, thereby restricting the application of these experimental results to rocks similar to oceanic ridge basalts in bulk composition.Because the Sr distribution coefficients are independent of oxygen fugacity, they may be used as geothermometers. If the temperature can be determined independently — for example, with the Sr distribution coefficients, the Eu distribution coefficients may be used as oxygen geobarometers. Throughout the range of oxygen fugacities ascribed to terrestrial and lunar basalts, plagioclase concentrates Eu but clinopyroxene rejects Eu.     

The system pargasite-H2O-CO2: a model for melting of a hydrous mineral with a mixed-volatile fluid—I. Experimental results to 8 kbar

The stability of the amphibole pargasite [NaCa₂Mg₄Al(Al₂Si₆))O₂₂(OH)₂] in the melting range has been determined at total pressures (P) of 1.2 to 8 kbar. The activity of H₂O was controlled independently of P by using mixtures of H₂O + CO₂ in the fluid phase. The mole fraction of H₂O in the fluid ($\text{X}{}_{\mathrm{H}}{}_2\text{O}{}^{\mathrm{1fl}}$) ranged from 1.0 to 0.2.At P < 4 kbar the stability temperature (T) of pargasite decreases with decreasing $\text{X}{}_{\mathrm{H}}{}_2\text{O}{}^{\mathrm{1fl}}$ at constant P. Above P ? 4 kbar stability T increases as $\text{X}{}_{\mathrm{H}}{}_2\text{O}{}^{\mathrm{1fl}}$ is decreased below one, passes through a T maximum and then decreases with a further decrease in $\text{X}{}_{\mathrm{H}}{}_2\text{O}{}^{\mathrm{1fl}}$. This behavior is due to a decrease in the H₂O content of the silicate liquid as $\text{X}{}_{\mathrm{H}}{}_2\text{O}{}^{\mathrm{1fl}}$ decreases. The magnitude of the T maximum increases from about 10°C (relative to the stability T for $\text{X}{}_{\mathrm{H}}{}_2\text{O}{}^{\mathrm{1fl}}\text{= 1}$) at P = 5 kbar to about 30°C at P = 8 kbar, and the position of the maximum shifts from $\text{X}{}_{\mathrm{H}}{}_2\text{O}{}^{\mathrm{1fl}}\text{? 0.6}$ at P = 5 kbar to $\text{X}{}_{\mathrm{H}}{}_2\text{O}{}^{\mathrm{1fl}}\text{? 0.4}$ at P = 8 kbar.The H₂O content of liquid coexisting with pargasite has been estimated as a function of at 5 and 8 kbar P, and can be used to estimate the H₂O content of magmas. Because pargasite is stable at low values of $\text{X}{}_{\mathrm{H}}{}_2\text{O}{}^{\mathrm{1fl}}$ at high P and T, hornblende can be an important phase in igneous processes even at relatively low H₂O fugacities.     

Revised values for the Gibbs free energy of formation of [Al(OH)4 aq−], diaspore,boehmite and bayerite at 298.15 K and 1 bar,the thermodynamic properties of kaolinite to 800 K and 1 bar,and the heats of solution of several gibbsite samples

Solution calorimetric measurements compared with solubility determinations from the literature for the same samples of gibbsite have provided a direct thermochemical cycle through which the Gibbs free energy of formation of [Al(OH)_{4 aq}^?] can be determined. The Gibbs free energy of formation of [Al(OH)_{4 aq}^?] at 298.15 K is ?1305 ± 1 kJ/mol. These heat-of-solution results show no significant difference in the thermodynamic properties of gibbsite particles in the range from 50 to 0.05 μm.The Gibbs free energies of formation at 298.15 K and 1 bar pressure of diaspore, boehmite and bayerite are ?9210 ± 5.0, ?918.4 ± 2.1 and ?1153 ± 2 kJ/mol based upon the Gibbs free energy of [A1(OH)_{4 aq}^?] calculated in this paper and the acceptance of ?1582.2 ± 1.3 and ?1154.9 ± 1.2 kJ/mol for the Gibbs free energy of formation of corundum and gibbsite, respectively.Values for the Gibbs free energy formation of [Al(OH)_{2 aq}⁺] and [AlO_{2 aq}^?] were also calculated as ?914.2 ± 2.1 and ?830.9 ± 2.1 kJ/mol, respectively. The use of [AlC_{2 aq}^?] as a chemical species is discouraged.A revised Gibbs free energy of formation for [H₄SiO_4aq⁰] was recalculated from calorimetric data yielding a value of ?1307.5 ± 1.7 kJ/mol which is in good agreement with the results obtained from several solubility studies.Smoothed values for the thermodynamic functions C_P⁰, ($\text{H}{}_{\mathrm{T}}{}^0\text{- H}{}_{298}{}^0\text{)}\text{T}$, $\text{(}\text{G}{}_{\mathrm{T}}{}^0\text{- H}{}_{298}{}^0\text{)}\text{T}$, S_T⁰ - S₀⁰, $\text{ΔH}{}_{\mathrm{?,298}}{}^0$ kaolinite are listed at integral temperatures between 298.15 and 800 K. The heat capacity of kaolinite at temperatures between 250 and 800 K may be calculated from the following equation: C_P⁰ = 1430.26 ? 0.78850 T + 3.0340 × 10^?4T² ?1.85158 × 10^?4T²$\text{1}\text{2}\text{+ 8.3341 × 10}{}^6\text{T}{}^{\mathrm{?2}}$.The thermodynamic properties of most of the geologically important Al-bearing phases have been referenced to the same reference state for Al, namely gibbsite.     

Liquidus relationships in the system CaCO3Ca(OH)2CaS and the solubility of sulfur in carbonatite magmas

CaCO₃Ca(OH)₂CaS serves as a model system for sulfide solubility in carbonatite magmas. Experiments at 1 kbar delineate fields for primary crystallization of CaCO₃, Ca(OH)₂ and CaS. The three fields meet at a ternary eutectic at 652°C with liquid composition (wt%): CaCO₃ = 46.1%, Ca(OH)₂ = 51.9%, CaS = 2.0%. Two crystallization sequences are possible for liquids that precipitate calcite, depending upon whether the liquid is on the low-CaS side, or the high-CaS side of the line connecting CaCO₃ to the eutectic liquid. Low-CaS liquids precipitate no sulfide until the eutectic temperature is reached leading to sulfide enrichment. The higher-CaS liquids precipitate some sulfide above the eutectic temperature, but the sulfide content of the melt is not greatly depleted as the eutectic temperature is approached. Theoretical considerations indicate that sulfide solubility in carbonate melts will be directly proportional to and inversely proportional to ; it also is likely to be directly proportional to melt basicity, defined here by . A strong similarity exists in the processes which control sulfide solubility in carbonate and in silicate melts. By analogy with silicates, ferrous iron, which was absent in our experiments, may also exert an important influence on sulfide solubility in natural carbonatite magmas.     

Olivine/silicate melt partitioning of germanium: an example of a nearly constant partition coefficient

The partitioning of germanium between forsterite (Fo) and liquids in the diopside-anorthiteforsterite join was investigated by electron microprobe analysis of Ge-doped samples equilibrated at 1300°–1450°C. Germanium is somewhat incompatible in Fo relative to the haplobasaltic melts, with a grand mean for all simple partition coefficients (D_Fo-l^Ge) of 0.68 ± 0.06. For the melt composition range studied, D_Fo-l^Ge is virtually constant in isothermal series of experiments, and shows only minor overall temperature dependence. The exchange reaction partition coefficient $\text{K}{}_{\mathrm{D}}\text{= (}\text{Mg}{}_2\text{GeO}{}_4\text{)}{}_{\mathrm{Fo}}\text{(SiO}{}_2\text{)}{}_{\mathrm{l}}\text{(Mg}{}_2\text{SiO}{}_4\text{)}{}_{\mathrm{Fo}}\text{(GeO}{}_2\text{)}{}_{\mathrm{l}}$] is near unity in all cases, with a grand mean of 0.93 ± 0.11. One exploratory run at 20 kbar yielded a distinctly lower partition coefficient (D_Fo-l^Ge = 0.54 ± 0.04), which confirms the negative pressure dependence predicted by the thermodynamics of Ge ai Si exchange.These new data indicate that absolute Ge enrichment must occur in terrestrial magmas undergoing olivine fractionation, while $\text{Ge}\text{Si}$ remains nearly constant.     

A valid Margules formulation for an asymmetric ternary solution: revision of the olivine-ilmenite thermometer,with applications

The olivine-ilmenite thermometer of Andersen and Lindsley (1979) was based on an incorrect formulation for the excess free energy of an asymmetric ternary solution. A valid formulation is derived and used to revise the parameters of the olivine-ilmenite thermometer. For olivine and ilmenite that have equilibrated above 700°C, temperature can be calculated from: $\text{T(°C) = ?273 +¦-12549 + P[0.03X}{}_{\mathrm{fa}}\text{+ 0.01099(X}{}_{\mathrm{gk}}\text{?X}{}_{\mathrm{il}}\text{)?0.062] + 10496 X}{}_{\mathrm{fa}}\text{+ 5767(X}{}_{\mathrm{gk}}\text{?X}{}_{\mathrm{il}}\text{) + X}{}_{\mathrm{hem}}\text{(38602?141550X}{}_{\mathrm{il}}\text{?47183X}{}_{\mathrm{gk}}\text{)|/[5.67?R ln K}{}_{\mathrm{D}}\text{+ 6.52X}{}_{\mathrm{fa}}\text{+ 3.09(X}{}_{\mathrm{gk}}\text{?X}{}_{\mathrm{il}}\text{) + X}{}_{\mathrm{hem}}\text{(16.49?109.46 X}{}_{\mathrm{il}}\text{?36.49X}{}_{\mathrm{gk}}\text{)]}$ with $\text{K}{}_{\mathrm{d}}\text{=}\text{(X}{}_{\mathrm{il}}\text{X}{}_{\mathrm{fo}}\text{)}\text{(X}{}_{\mathrm{gk}}\text{X}{}_{\mathrm{fa}}\text{)}$. The revised model gives W^il·gk_G = 5767?3.09T + 0.011P and ΔG_exch = 7301 ? 8.9T ? 0.047P (T in K, P in bars). Applications include Apollo 17 breccias and kimberlites.     

Diffusion of ions in sea water and in deep-sea sediments

The tracer-diffusion coefficient of ions in water, D_j⁰, and in sea water, $\text{D}{}_{\mathrm{j}}{}^1$, differ by no more than zero to 8 per cent. When sea water diffuses into a dilute solution of water, in order to maintain the electro-neutrality, the average diffusion coefficients of major cations become greater but of major anions smaller than their respective $\text{D}{}_{\mathrm{j}}{}^1$ or D_j⁰ values. The tracer diffusion coefficients of ions in deep-sea sediments, D_j,sed., can be related to $\text{D}{}_{\mathrm{j}}{}^1$ by $\text{D}{}_{\mathrm{j,sed.}}\text{= D}{}_{\mathrm{j}}{}^1\text{·}\text{α}\text{θ}{}^2$, where θ is the tortuosity of the bulk sediment and a a constant close to one.     

The effect of pressure on the solubility of minerals in water and seawater

The effect of presure on the solubility of minerals in water and seawater can be estimated from In

$\text{(}\text{K}{}^{\mathrm{P}}{}_{\mathrm{sp}}\text{K}{}^0{}_{\mathrm{sp}}\text{) +}\text{(?ΔVP + 0.5ΔKP}{}^2\text{)}\text{RT}$

where the volume (ΔV) and compressibility (ΔK) changes at atmospheric pressure (P = 0) are given by

$\text{ΔV =}\text{V}\text{?}\text{(M}{}^{\mathrm{+}}\text{, X}{}^{\mathrm{?}}\text{) ?}\text{V}\text{?}\text{[MX(s)]ΔK =}\text{K}\text{?}\text{(M}{}^{\mathrm{+}}\text{, X}{}^{\mathrm{?}}\text{) ?}\text{K}\text{?}\text{[MX(s)]}$

Values of the partial molal volume ($\text{V}\text{?}$) and compressibilty ($\text{K}\text{?}$) in water and seawater have been tabulated for some ions from 0 to 50°C. The compressibility change is quite large (~10 × 10^?3 cm³ bar^?1 mol^?1) for the solubility of most minerals. This large compressibility change accounts for the large differences observed between values of ΔV obtained from linear plots of In K_sp versus P and molal volume data (Macdonald and North, 1974; North, 1974). Calculated values of $\text{K}{}^{\mathrm{P}}{}_{\mathrm{sp}}\text{K}{}^{\mathrm{o}}{}_{\mathrm{sp}}$ for the solubility of CaCO₃, SrSO₄ and CaF₂ in water were found to be in good agreement with direct measurements (Macdonald and North, 1974). Similar calculations for the solubility of minerals in seawater are also in good agreement with direct measurements (Ingle, 1975) providing that the surface of the solid phase is not appreciably altered.     

Effect of pressure on the diffusivity of network-forming cations in melts of jadeitic compositions

The diffusivities of network-forming cations (Si⁴⁺, Al³⁺, Ge⁴⁺ and Ga³⁺) in melts of the jadeitic composition NaAl(Si, Ge)₂O₆ and Na(Al, Ga)Si₂O₆ have been measured at pressures between 6 and 20 kbar at 1400°C. The rates of interdiffusion of Si⁴⁺-Ge⁴⁺ and Al³⁺-Ge³⁺ increase with increasing pressure at constant temperature. The results are consistent with the ion-dynamics computer simulations of Jadeite melt by Angellet al. (1982, 1983). The coefficient measured for the Si⁴⁺-Ge⁴⁺ interdiffusion is between 8 × 10^?10 and 2.5 × 10^?8$\text{cm}{}^2\text{sec}$ at 6 kbar, depending on the composition of the melt, whereas at 20 kbar it is between 7 × 10^?9 and 2 × 10^?7$\text{cm}{}^2\text{sec}$. The effect of pressure is greater for more Si-rich compositions (i.e., closer to NaAlSi₂O₆ composition). The coefficient measured for the Al³⁺-Ga³⁺ inter- diffusion is between 9 × 10^?10 and 3 × 10^?9 cm²/sec at 6 kbar and between 3 × 10^?9 and 1 × 10^?8$\text{cm}{}^2\text{sec}$ at 20 kbar. The rate of increase in diffusivity with pressure of Al³⁺-Ga³⁺ (a factor of 3–4) is smaller than that of Si⁴⁺-Ge⁴⁺ (a factor of 7–17).The Si⁴⁺-Ge⁴⁺ interdiffusion in melts of Na₂O · 4(Si, Ge)O₂ composition has also been measured at 8 and 15 kbar for comparison. The effect of pressure on the diffusivity in this melt is significantly smaller than that for the jadeitic melts. The increase in diffusivity of the network-forming cations in jadeitic melts with increasing pressure may be related to the decrease in viscosity of the same melt. The present results, as well as the ion-dynamics simulations, suggest that the homogenization of partial melts and mixing of magmas would be more efficient at greater depths.     

Water and magmas; a mixing model

A model for the mixing of H₂O and silicate melts has been derived from the experimentally determined effects of H₂O on the viscosity (fluidity), volumes, electrical conductivities, and especially the thermodynamic properties of hydrous aluminosilicate melts. It involves primarily the reaction of H₂O with those O^?2 ions of the melt that are shared (bridging) between adjacent (Al, Si)O₄ tetrahedra to produce OH^? ions. However, in those melts that contain trivalent ions in tetrahedral coordination, such as the Al³⁺ ion in feldspathic melts, the model further involves exchange of a proton from H₂O with a non-tetrahedrally coordinated cation that must be present to balance the net charge on the AlO₄ group. This cation exchange reaction, which goes essentially to completion, results in dissociation of the H₂O and is limited only by the availability of H₂O and the number of exchangeable cations per mole of aluminosilicate.In the system NaAlSi₃O₈-H₂O, upon which this thermodynamic model is based, there is 1 mole of exchangeable cations (Na⁺) per mole (GFW) of NaAlSi₃O₈, consequently ion exchange occurs for H₂O contents up to a 1:1 mole ratio (X^m_w = mole fraction H₂O = 0.5). For mole fractions of H₂O greater than 0.5, no further exchange can occur and the reaction with additional bridging oxygens of the melt produces 2 moles of associated OH^? ions per mole of H₂O dissolved. These reactions lead to a linear dependence of the thermodynamic activity of H₂O (a^m_w) on the square of its mole fraction (X^m_w) for values of X^m_w, up to 0.5 and an exponential dependence on X^m_w at higher H₂O contents. Thus, for values of $\text{X}{}^{\mathrm{m}}{}_{\mathrm{w}}\text{? 0.5, a}{}^{\mathrm{m}}{}_{\mathrm{w}}\text{= k(X}{}^{\mathrm{m}}{}_{\mathrm{w}}\text{)}{}^2$, where k is a Henry's law constant for the dissociated solute.Extension of the thermodynamic model for NaAlSi₃O₈-H₂O to predict H₂O solubilities and other behavior of compositionally more complex aluminosilicate melts (magmas) requires placing these melts on an equimolal basis with NaAlSi₃O₈. This is readily accomplished using chemical analyses of quenched glasses by normalizing to the stoichiometric requirements of NaAlSi₃O₈, first in terms of equal numbers of exchangeable cations for mole fractions of H₂O up to 0.5 and secondly in terms of 8 moles of oxygen for higher H₂O contents. Chemical analyses of three igneous-rock glasses, ranging in composition from tholeiitic basalt to lithium-rich pegmatite, were thus recast and the experimental H₂O solubilities were computed on this equimolal basis. The resulting equimolal solubilities are all the same, within experimental error, as the solubility of H₂O in NaAlSi₃O₈ melt calculated from the thermodynamic relations.The equivalence of equimolal solubilities implies that the Henry's law constant (k), which is a function of temperature and pressure, is independent of aluminosilicate composition over a wide range. Moreover, as a consequence of the Gibbs-Duhem relation and the properties of exact differentials, it is clear that the silicate components of the melt, properly defined, mix ideally. Thus, a relatively simple mixing model for H₂O in silicate melts has led to a quantitative thermodynamic model for magmas that has far-reaching consequences in igneous petrogenesis.     

Effet de la température sur les distributions de Na,Rb et Cs entre la sanidine,la muscovite,la phlogopite et une solution hydrothermale sous une pression de 1 kbar

The distribution of trace amounts of Na, Rb and Cs, between muscovite, phlogopite, sanidine and hydrothermal solution have been studied by ion exchange in a temperature range from 400 to 800°C.These distributions have been expressed with a partition ratio P^aq?m_x = $\text{(}\text{X}\text{K}\text{)}{}^{\mathrm{aq}}\text{(}\text{X}\text{K}\text{)}{}^{\mathrm{m}}$ (where X is Na, Rb or Cs).In the case of Na and Cs in muscovite, even for the dilute solutions, the ratio P^aq?m_x is not the equilibrium constant k_x of exchange reactions. In other cases, P^aq?m_x does not depend on the trace alkali ion concentration in silicates (X) and is equal to k_x. Variations of P_x or k_x with T are greater for Na and Cs than for Rb. Generally, k_x decreases with increase in T. The function log $\text{P}{}_{\mathrm{x}}\text{= f(}\text{1}\text{T}\text{)}$ is not linear for Na or Cs, but in the case of Rb, $\text{f(}\text{1}\text{T}\text{)}$ is linear and the standard enthalpy and entropy of exchange reactions have been estimated by applying the Arrhenius relation.The distribution relations obtained between silicate and vapour phase permit the determination of distributions of Na, Rb and Cs between two minerals mI and mII, relative to K. These have been expressed with the partition ratio Q_x =$\text{(}\text{X}\text{K}\text{)}{}^{\mathrm{mI}}\text{(}\text{X}\text{K}\text{)}{}^{\mathrm{mII}}$. Variations of Q_x with T are not remarkable, and even for Rb between phlogopite and feldspar are negligible. Nevertheless, one may use the distributions of Rb and Cs between muscovite and feldspar for geothermometry. Experimental results have been applied to some rocks by effecting corrections from the major element composition of the natural minerals. Estimated temperatures are near to 400°C in the granites and pegmatite studied here.     

Dissociation constants for alkali earth and sodium borate ion pairs from 10 to 50°C

Potentiometric measurements in dilute sodium borate solutions with added alkali earth chlordie salts yield the following expressions for the dissociation constants of alkali earth borate ion pairs from 10 to 50°C:

$\text{p}\text{K(}\text{MgH}{}_2\text{BO}{}_3{}^{\mathrm{+}}\text{=1.266+0.001204 T}$

$\text{p}\text{K(}\text{CaH}{}_2\text{BO}{}_3{}^{\mathrm{+}}\text{=1.154+0.002170 T}$

$\text{p}\text{K(}\text{SrH}{}_2\text{BO}{}_3{}^{\mathrm{+}}\text{=1.033+0.001738 T}$

$\text{p}\text{K(}\text{BaH}{}_2\text{BO}{}_3{}^{\mathrm{+}}\text{=1.942+0.001850 T}$

where T is in °K. Enthalpies for the dissociation reactions at 25°C are less than 1 kcal./mole for all the alkali earth borate ion pairs.Values for pK(NaH₂BO₃°) from 5 to 55°C computed from the experimental data of Owen and King are in good agreement with those determined potentiometrically. The average value from both methods is 0.22 ± 0.1 at 25°C.Application to seawater of computed pK's for MgH₂BO₃⁺, CaH₂BO₃⁺ and NaH₂BO₃⁰ yields an apparent dissociation constant for boric acid of 8.73 vs. 8.70 measured by Lyman, 8.68 by Buch and 8.73 by Byrne and Kester.     

On calculating activity coefficients and other excess functions from the intracrystalline exchange properties of a double-site phase

For a phase at equilibrium in which two cation species are partitioned ideally between two sub-lattice sites, the excess functions of mixing (free energy, enthalpy and entropy) are directly related to the bulk composition of the phase and ΔG_E°(T, P), the standard-state intra- crystalline exchange free energy. If the phase is not at equilibrium internally, an additional ordering parameter is necessary to fix the excess free energy of mixing, G_mix^EX, unambiguously. Conversely, for any fixed G_mix^EX there exists an infinity of possible intracrystalline cation dis- tributions, only one of which is the equilibrium distribution for the specified temperature and pressure. As ideal intraphase cation ordering becomes more pronounced, G_mix^EX decreases. In response, the total free energy of mixing for the phase decreases progressively for non-end member compositions, approaching, at the limits of ordering, values appropriate for stabilizing compounds of intermediate composition.The model-dependent activity coefficient for component A in the phase, γ_A^T, can be calculated for any bulk composition, X_A^T, either from G_mix^EX directly or from more basic equations involving the interrelation of chemical potentials at equilibrium. A general form for γ_A^T is ln $\text{γ}{}_{\mathrm{A}}{}^{\mathrm{T}}\text{= 1n[}\text{2(X}{}_{\mathrm{A}}{}^{\mathrm{\alpha }}\text{X}{}_{\mathrm{A}}{}^{\mathrm{\beta }}\text{)}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}\text{/}\text{(X}{}_{\mathrm{A}}{}^{\mathrm{\alpha }}\text{+X}{}_{\mathrm{A}}{}^{\mathrm{\beta }}\text{)}\text{]+Y}$, where Xj^κ denotes the mole fraction of species j in site κ. The first term on the right-hand side of this equation is the contribution to γ_A^T from ideal intracrystalline partitioning, and is common to the several theories lately presented to model intraphase cation partitioning. It can be shown rigorously that this term contributes to a negative deviation from ideality for the bulk phase. The second term is the contribution to the macroscopic activity coefficient from non-ideal intraphase partitioning, and is related to an enthalpy of mixing, H_mix^N in excess of that resulting from ideal inter-site cation ordering. While the expression represented by Y can take several functional forms, the additional enthalpy can be evaluated explicitly for specific non-ideal partitioning models from the relation $\text{H}{}_{\mathrm{mix}}{}^{\mathrm{N}}\text{= 2RT(1? X}{}_{\mathrm{A}}{}^{\mathrm{T}}\text{) ∝}\text{Y}\text{(1 ? X}{}_{\mathrm{A}}{}^{\mathrm{T}}\text{)}{}^2\text{dX}{}_{\mathrm{A}}{}^{\mathrm{T}}$.In those cases, G_mix^EX can also be determined exactly.