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.     

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.     

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.     

Thermochemistry of the high structural state plagioclases

The enthalpies of solution of a suite of 19 high-structural state synthetic plagioclases were measured in a Pb₂B₂O₅ melt at 970 K. The samples were crystallized from analyzed glasses at 1200°C and 20 kbar pressure in a piston-cylinder apparatus. A number of runs were also made on Amelia albite and Amelia albite synthetically disordered at 1050–1080°C and one bar for one month and at 1200°C and 20 kbar for 10 hr. The component oxides of anorthite, CaO, Al₂O₃ and SiO₂, were remeasured.The ΔH of disorder of albite inferred in the present study from albite crystallized from glass is 3.23 kcal, which agrees with the 3.4 found by Holm and Kleppa (1968). It is not certain whether this value includes the ΔH of a reversible displacive transition to monoclinic symmetry, as suggested by Helgesonet al. (1978) for the Holm-Kleppa results. The enthalpy of solution value for albite accepted for the solid solution series is based on the heat-treated Amelia albite and is 2.86 kcal less than for untreated Amelia albite.The enthalpy of formation from the oxides at 970 K of synthetic anorthite is ?24.06 ± 0.31 kcal, significantly higher than the ?23.16 kcal found by Charluet al. (1978), and in good agreement with the value of ?23.89 ± 0.82 given by Robieet al. (1979), based on acid calorimetry.The excess enthalpy of mixing in high plagioclase can be represented by the expression, valid at 970 K: ΔH^ex(±0.16 kcal) = 6.7461 X_abX²_An + 2.0247 X_AnX²_Ab where X_Ab and X_An are, respectively, the mole fractions of NaAlSi₃O₈ and CaAl₂Si₂O₈. This ΔH^ex, together with the mixing entropy of Kerrick and Darken's (1975) Al-avoidance model, reproduces almost perfectly the free energy of mixing found by Orville (1972) in aqueous cation-exchange experiments at 700°C. It is likely that Al-avoidance is the significant stabilizing factor in the high plagioclase series, at least for X_An≥ 0.3. At high temperatures the plagioclases have nearly the free energies of ideal one-site solid solutions. The Al-avoidance model leads to the following Gibbs energy of mixing for the high plagioclase series: $\text{ΔG}{}^{\mathrm{mix}}\text{= ΔH}{}^{\mathrm{ex}}\text{+ RT X}{}_{\mathrm{Ab}}\text{ln[X}{}^2{}_{\mathrm{Ab}}\text{(2 ? X}{}_{\mathrm{Ab}}\text{)]+ X}{}_{\mathrm{An}}\text{ln}\text{[X}{}_{\mathrm{An}}\text{(1+X}{}_{\mathrm{An}}\text{)}{}^2\text{]}\text{4}$. The entropy and enthalpy of mixing should be very nearly independent of temperature because of the unlikelihood of excess heat capacity in the albite-anorthite join.     

A detailed Pb isotopic study of crustal contamination/assimilation: The Edgecumbe volcanic field,SE Alaska

To better understand the process of crustal contamination/assimilation, 23 Pb isotopic compositions and 12 concentrations have been measured on lavas and basement rocks from the Edgecumbe volcanic field, SE Alaska. Measured isotopic ratios have the following ranges: ${}^{206}\text{Pb}{}^{204}\text{Pb}\text{= 18.477–19.161}$; ${}^{207}\text{Pb}{}^{204}\text{Pb}\text{= 15.562–15.679}$; ${}^{208}\text{Pb}{}^{204}\text{Pb}\text{= 38.17–38.85}$. While the data form well-constrained linear arrays on Pb-Pb diagrams, no simple correlation exists with major element composition. Basaltic lavas (≤ 51 wt% SiO₂) are characterized by two isotopic groups. The olivine basalt (≤ 48% SiO₂) is more radiogenic than the plagioclase basalt (48–51%) which also shows more heterogeneity. In the silica range 52–55%, Pb isotopic ratios increase significantly but remain fairly constant in the range 55–70% SiO₂. Lead concentrations vary from 1 ppm in the basalts to 7 ppm in the rhyodacites. Analyzed basement rocks are more radiogenic than any of the lavas $\text{(}{}^{206}\text{Pb}{}^{204}\text{Pb}\text{= 19.20}$; ${}^{207}\text{Pb}{}^{204}\text{Pb}\text{= 15.65}$; ${}^{208}\text{Pb}{}^{204}\text{Pb}\text{= 38.86}$. The Pb isotopic data are qualitatively consistent with the contamination process described by Myerset al. (1984). However, because of fundamental differences in the mixing relations between the Sr system studied earlier and the Pb system, the new Pb data have revealed details of the process not apparent from the Sr data alone. In particular, it has been shown that the parent magma was more primitive than originally assumed, and that two contamination events are recorded in the lavas. The first event, involving a mafic parent and different crustal contaminants, produced the intermediate and siliceous hybrids in cupolas located above the main basaltic chamber. The types of country rock intruded as well as the degree of partial fusion achieved in individual cupolas controlled the range of hybrid compositions produced while the eruption sequence was determined by the order in which the cupolas were tapped. The second contamination event produced the plagioclase basalt, the most voluminous basaltic unit, by mixing the mafic parent with the olivine basalt, an independent, primary magma. Our results suggest crustal contamination models that assume bulk assimilation of crustal end members may be too simplistic.     

The kinetics of silica-water reactions

A differential rate equation for silica-water reactions from 0–300°C has been derived based on stoichiometry and activities of the reactants in the reaction SiO₂(s) + 2H₂O(l) = H₄SiO₄(aq)

where ($\text{A}\text{M}$) = (the relative interfacial area between the solid and aqueous phases/the relative mass of water in the system), and k₊ and k_? are the rate constants for, respectively, dissolution and precipitation. The rate constant for precipitation of all silica phases is $\text{log k}{}_{\mathrm{?}}\text{= ? 0.707 ?}\text{2598}\text{T}\text{(T, K)}$ and E_act for this reaction is 49.8 kJ mol^?1. Corresponding equilibrium constants for this reaction with quartz, cristobalite, or amorphous silica were expressed as $\text{log K = a + bT +}\text{c}\text{T}$. Using $\text{K =}\text{k}{}_{\mathrm{+}}\text{k}{}_{\mathrm{?}}$, k was expressed as $\text{log k}{}_{\mathrm{+}}\text{= a + bT +}\text{c}\text{T}$ and a corresponding activation energy calculated:

     

8.

Liquidus relationships in the system CaCO₃Ca(OH)₂CaS and the solubility of sulfur in carbonatite magmas     

George R. Helz  Peter J. Wyllie 《Geochimica et cosmochimica acta》1979,43(2):259-265

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.     

9.

Calcium diffusion in a simple silicate melt to 30 kbar     

E.Bruce Watson 《Geochimica et cosmochimica acta》1979,43(3):313-322

Calcium-45 was used as a radiotracer to measure self-diffusion coefficients for Ca in a sodium-calcium-aluminosilicate melt (29% Na₂O, 5% CaO, 10% Al₂O₃, 56% SiO₂) at temperatures in the range 1100–1400°C and pressures to 30 kbar. Calcium diffusivity (D_Ca) was found to depend upon both temperature and pressure in a complex but systematic manner: $\text{(}\text{?D}{}_{\mathrm{Ca}}\text{?P}\text{)}{}_{\mathrm{T}}$ is always negative and has a larger absolute value at lower temperatures; $\text{(}\text{?D}{}_{\mathrm{Ca}}\text{?T}\text{)}{}_{\mathrm{P}}$ is positive and increases with increasing pressure. The overall dependence of D_Ca upon T and P is given approximately by $\text{D}{}_{\mathrm{caT.P}}\text{= [0.0025 exp(}\text{-23,107}\text{RT}\text{)] exp [P}\text{(0.7297T ? 1261.32)}\text{RT}\text{]}$. When expressed in terms of volume (V_a) and energy (E) of activation, the results are as follows: V_a ranges from 2.2 cm³/mole at 1400°C to 11.9 cm³/mole at 1100°C. and E ranges from 25.4 kcal/mole (1 kban to 49.8 kcal/mole (20 kbar).From the systematic dependence of D_Ca upon T and P, it is concluded that diffusion of Ca²⁺ in silicate melts does not take place by means of a vacant site mechanism, but is controlled instead by the amount and distribution of free volume in the melt structure.If it is assumed that the viscosity of the melt used in this study decreases with increasing pressure (Kushiro, 1976, J. Geophys. Res.81, 6351–6356) as D_Ca does, then the Stokes-Einstein inverse relation between viscosity and diffusivity is clearly violated, and its validity for silicate melts must be questioned. Thus, it appears that in silicate melts, unlike many liquids, viscous flow and diffusion are fundamentally different transport processes, involving different structural units.The effect of pressure on calcium diffusion is too small to invalidate kinetic models of upper mantle processes that have been based upon diffusivity values measured at 1 atm. Pressure may, however, induce significant reductions in the diffusion rates of large ions such as Rb⁺ or SiO^4?₄ in silicate melts.     

10.

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

John Geover 《Geochimica et cosmochimica acta》1974,38(10):1527-1548

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.     

11.

Serpentinization of cumulate ultramafic rocks from the North Arm Mountain massif of the Bay of Islands ophiolite   总被引：1，自引：0，他引：1  

Stephen C. Komor  Don Elthon  John F. Casey 《Geochimica et cosmochimica acta》1985,49(11):2331-2338

Partially serpentinized dunites and wehrlites comprise the bulk of the cumulate ultramafic unit at the North Arm Mountain massif of the Bay of Islands ophiolite complex, Newfoundland. In a suite of 59 dunites and werhlites from the base of the unit, the serpentinized portions consist of lizardite + chrysotile + brucite + (accessory) magnetite. The ratio of $\text{(}\text{lizardite}\text{+}\text{chrysotile}\text{)}\text{to brucite}\text{= ~8:2 (}\text{weight percent}\text{)}$. Petrographic observations show that most serpentinization occurred at the expense of olivine; only limited amounts of clinopyroxene were serpentized. An estimated volume increase of 32% accompanied serpentinization of the peridotites. Reconstructions of the primary modal proportions of wehrlites (made taking this volume increase into account) contain an average of 6% more clinopyroxene and 6% less olivine than do modal reconstructions that ignore the volume increase. Mass balance calculations provide no clear evidence for appreciable metasomatism of Al₂O₃, CaO, FeO, MgO, or SiO₂ during Serpentinization. The presence of brucite, the evidence that most serpentinization occurred at the expense of olivine, and the lack of appreciable metasomatism, suggest that the primary reaction that controlled serpentinization of the peridotites is: $\text{2Mg}{}_2\text{SiO}{}_2\text{+ 3H}{}_2\text{O ? Mg}{}_3\text{Si}{}_2\text{O}{}_5\text{(OH)}{}_4\text{+ Mg(OH)}{}_2$. olivine added serpentine brucite     

12.

The partial molal volume of silicic acid in 0.725 M NaCl at 1°C determined by the neutralization of Na₂SiO₃     

J.Peter Hershey  Iver W Duedall 《Geochimica et cosmochimica acta》1983,47(11):1931-1936

The partial molal volume of silicic acid ($\text{V}\text{?}\text{(Si(OH)}{}_4\text{)}$) in 0.725 M NaCl at 1°C was calculated from the measured volume change (Δ$\text{V}\text{?}{}_{\mathrm{n}}$) due to the neutralization of anhydrous sodium metasilicate with HCl and the $\text{V}\text{?}\text{(HCl)}$ and $\text{V}\text{?}\text{(NaCl)}$ obtained from the literature. $\text{V}\text{?}\text{(Si(OH)}{}_4\text{) = 59.0 cm}{}^3\text{mol}\text{? 1}$, determined under experimental conditions of pH = 2.2, compares favorably with $\text{V}\text{?}\text{(Si(OH)}{}_4\text{) = 58.9 cm}{}^3\text{mol}{}^{\mathrm{?1}}$ calculated from the measured volume change due to the hydrolysis of the meta-silicate salt at pH = 11 and from the partial molal volume due to electrostriction ($\text{V}\text{?}{}_{\mathrm{elect}}$) of water by charged Si species present in the solution at the high pH. This agreement lends support to a semiempirical model for calculating $\text{V}\text{?}{}_{\mathrm{elect}}$ in developed by Millero (1969). $\text{V}\text{?}\text{(NaOH) = ? 5.45 cm}{}^3\text{mol}{}^{\mathrm{?1}}$ in 0.725 M NaCl needed for this calculation was also determined in this work. The rate of polymerization of Si(OH)₄ at 1°C was monitored to insure that the monomer Si(OH)₄ was the main Si species present during the determination of $\text{V}\text{?}\text{(Si(OH)}{}_4\text{)}$ by neutralization of the alkali silicate. $\text{V}\text{?}\text{(Si(OH)}{}_4\text{)}$ determined in this study compares favorably with the value calculated from high pressure solubility measurements.     

13.

Lead orthophosphates—II. Stability of cholopyromophite at 25°c     

Jerome O. Nriagu 《Geochimica et cosmochimica acta》1973,37(3):367-377

The conversion of secondary lead orthophosphate [PbHPO₄] into chloropyromorphite [Pb₅(PO₄)₃Cl] in ca. 10^?1 M NaCl solutions has been investigated at 25°C. From the composition of the supernatant solutions, the solubility product constant for Pb₅(PO₄)₃Cl has been calculated to be 10^?84.4±0.1, corresponding to $\text{ΔG}{}_{\mathrm{?}}\text{°}$ of ?906.2 kcal mol^?1. The solution equilibria and phase relationships in the system PbCl₂-PbO-P₂O₈-H₂O are discussed along with the geological implications.     

14.

Use of Pitzer's equations to determine the media effect on the formation of lead chloro complexes     

Frank.J. Millero  Robert H. Byrne 《Geochimica et cosmochimica acta》1984,48(5):1145-1150

The spectrophotometric measurements of chloro complexes of lead in aqueous HCl, NaCl, MgCl₂ and CaCl₂ solutions at 25°C have been analyzed using Pitzer's specific interaction equations. Parameters for activity coefficients of the complexes PbCl⁺, PbCl₂⁰ and PbCl₃^? have been determined for the various media. Values of $\text{K}{}_1\text{= 30.0 ± 0.6}$, $\text{K}{}_2\text{= 106.7 ± 2.1}$ and $\text{K}{}_3\text{= 73.0 ± 1.5}$ were obtained for the cumulative formation constants. [$\text{Pb}{}^{\mathrm{2+}}\text{+ nCl}{}^{\mathrm{?}}\text{→ PbCl}{}_{\mathrm{n}}{}^{\mathrm{2?n}}\text{)}$]. These values are in reasonable agreement with literature data. The Pitzer parameters for the PbCl ion pairs in various media were used to calculate the speciation of Pb²⁺ in an artificial seawater solution.     

15.

Mineral-solution equilibria—III. The system Na₂OAl₂O₃SiO₂H₂OHCl     

Robert K. Popp  John D. Frantz 《Geochimica et cosmochimica acta》1980,44(7):1029-1037

Chemical equilibrium between sodium-aluminum silicate minerals and chloride bearing fluid has been experimentally determined in the range 500–700°C at 1 kbar, using rapid-quench hydrothermal methods and two modifications of the Ag + AgCl acid buffer technique. The temperature dependence of the thermodynamic equilibrium constant (K) for the reaction $\text{NaAlSi}{}_3\text{O}{}_8\text{+ HCl}{}^{\mathrm{o}}\text{= NaCl}{}^{\mathrm{o}}\text{+}\text{1}\text{2Al}{}_2\text{SiO}{}_5\text{, +}\text{5}\text{2SiO}{}_2\text{+}\text{1}\text{2H}{}_2\text{O}$ Albite Andalusite Qtz. can be described by the following equation: log k = ?4.437 + 5205.6/T(K) The data from this study are consistent with experimental results reported by Montoya and Hemley (1975) for lower temperature equilibria defined by the assemblages albite + paragonite + quartz + fluid and paragonite + andalusite + quartz + fluid. Values of the equilibrium constants for the above reactions were used to estimate the difference in Gibbs free energy of formation between NaCl^o and HCl^o in the range 400–700°C and 1–2 kbar. Similar calculations using data from phase equilibrium studies reported in the literature were made to determine the difference in Gibbs free energy of formation between KCl^o and HCl^o. These data permit modelling of the chemical interaction between muscovite + kspar + paragonite + albite + quartz assemblages and chloride-bearing hydrothermal fluids.     

16.

Mineral chemistry and crystallization sequences in kimberlite and lamproite dikes from the Sisimiut area,central West Greenland     

《Lithos》1987,20(5):391-417

     

17.

Rutile solubility and titanium coordination in silicate melts     

James E Dickinson  Paul C Hess 《Geochimica et cosmochimica acta》1985,49(11):2289-2296

The solubility of rutile has been determined in a series of compositions in the K₂O-Al₂O₃-SiO₂ system ($\text{K}{}^1\text{=}\text{K}{}_2\text{O}\text{(K}{}_2\text{O}$ + Al₂O₃) = 0.38–0.90), and the CaO-Al₂O₃-SiO₂ system ($\text{C}{}^1\text{=}\text{CaO}\text{(CaO + Al}{}_2\text{O}{}_3\text{)}\text{= 0.47–0.59}$). Isothermal results in the KAS system at 1325°C, 1400°C, and 1475°C show rutile solubility to be a strong function of the $\text{K}{}^1$ ratio. For example, at 1475°C the amount of TiO₂ required for rutile saturation varies from 9.5 wt% ($\text{K}{}^1\text{= 0.38}$) to 11.5 wt% ($\text{K}{}^1\text{= 0.48}$) to 41.2 wt% ($\text{K}{}^1\text{= 0.90}$). In the CAS system at 1475°C, rutile solubility is not a strong function of $\text{C}{}^1$. The amount of TiO₂ required for saturation varies from 14 wt% ($\text{C}{}^1\text{= 0.48}$) to 16.2 wt% ($\text{C}{}^1\text{= 0.59}$).The solubility changes in KAS melts are interpreted to be due to the formation of strong complexes between Ti and K⁺ in excess of that needed to charge balance Al³⁺. The suggested stoichiometry of this complex is K₂Ti₂O₅ or K₂Ti₃O₇. In CAS melts, the data suggest that Ca²⁺ in excess of A1³⁺ is not as effective at complexing with Ti as is K⁺. The greater solubility of rutile in CAS melts when $\text{C}{}^1$ is less than 0.54 compared to KAS melts of equal $\text{K}{}^1$ ratio results primarily from competition between Ti and Al for complexing cations (Ca vs. K).TiK_β x-ray emission spectra of KAS glasses ($\text{K}{}^1\text{= 0.43–0.60}$) with 7 mole% added TiO₂, rutile, and Ba₂TiO₄, demonstrate that the average Ti-O bond length in these glasses is equal to that of rutile rather than Ba₂TiO₄, implying that Ti in these compositions is 6-fold rather than 4-fold coordinated. Re-examination of published spectroscopic data in light of these results and the solubility data, suggests that the 6-fold coordination polyhedron of Ti is highly distorted, with at least one Ti-O bond grossly undersatisfied in terms of Pauling's rules.     

18.

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

Y Bottinga  D Weill  P Richet 《Geochimica et cosmochimica acta》1982,46(6):909-919

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.     

19.

Viscosity and configurational entropy of silicate melts     

Pascal Richet 《Geochimica et cosmochimica acta》1984,48(3):471-483

With the configurational entropy theory of relaxation processes of Adam and Gibbs (1965), one predicts that the viscosity depends on temperature according to log $\text{η = A}{}_{\mathrm{e}}\text{+}\text{B}{}_{\mathrm{e}}\text{TS}{}_{\mathrm{<mtext>conf</mtext>}}$, where S_conf is the configurational entropy of the liquid. Thermochemical calculations of S_conf performed for some mineral compositions show the importance of non-configurational contributions to the entropy differences between amorphous and crystalline phases. Except for the case of SiO₂, the available thermodynamic data indicate that the above equation for viscosity accounts quantitatively for the experimentally determined temperature dependence of the viscosity of silicate melts. The Adam and Gibbs theory also provides a simple rationale for the non linear variation of the logarithmic viscosity with composition in mixed alkali silicate liquids at low temperatures, the minimum of viscosity resulting from the contribution of the entropy of mixing to S_conf.     

20.

Etude des inclusions fluides du système N₂-CO₂ de dolomites et de quartz de Tunisie septentrionale. Données de la microcryoscopie et de l'analyse à la microsonde à effet Raman     

N Guilhaumou  P Dhamelincourt  J-C Touray  J Touret 《Geochimica et cosmochimica acta》1981,45(5):657-673

Optical and analytical studies were performed on 400 N₂ + CO₂ gas bearing inclusions in dolomites and quartz from Triassic outcrops in northern Tunisia. Other fluids present include brines (NaCl and KCl bearing inclusions) and rare liquid hydrocarbons. At the time of trapping, such fluids were heterogeneous gas + brine mixtures. In hydrocarbon free inclusions the $\text{N}{}_2\text{(N}{}_2\text{+ CO}{}_2\text{)}$ mole ratio was determined using two different non-destructive and punctual techniques: Raman microprobe analysis, and optical estimation of the volume ratios of the different phases selected at low temperatures. In the observed range of compositions, the two methods agree reasonably well.The N₂ + CO₂ inclusions are divided into three classes of composition: (a) $\text{N}{}_2\text{(N}{}_2\text{+ CO}{}_2\text{)}\text{> 0,57}$: Liquid nitrogen is always visible at very low temperature and homogenisation occurs in the range ?151°C to ? 147°C (nitrogen critical temperature) dry ice (solid CO₂) sublimates between ?75°C and ?60°C; (b) $\text{0,20 <}\text{N}{}_2\text{(N}{}_2\text{+ CO}{}_2\text{)}\text{? 0,57}$: liquid nitrogen is visible at very low temperature but dry ice melts on heating; liquid and gas CO₂ homogenise to liquid phase between ?51°C to ?22°C; (c) $\text{N}{}_2\text{(N}{}_2\text{+ CO}{}_2\text{)}\text{? 0,20}$: liquid nitrogen is not visible even at very low temperature (?195°C) and liquid and gas CO₂ homogenise to liquid phase between ?22°C and ?15°C. The observed phases changes are used to propose a preliminary phase diagram for the system CO₂-N₂ at low temperatures.Assuming additivity of partial pressures, isochores for the CO₂-N₂ inclusions have been computed. The intersection of these isochores with those for brine inclusions in the same samples may give the P and T of trapping of the fluids.     

	a	b	c	Eact(kJ mol -1)
Quarts	1.174	-2.028 x 103	-4158	67.4–76.6
α-Cristobalite	-0.739	0	-3586	68.7
β-Cristobalite	-0.936	0	-3392	65.0
Amorphous silica	-0.369	-7.890 x 10-4	3438	60.9–64.9

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.     

Calcium diffusion in a simple silicate melt to 30 kbar

Calcium-45 was used as a radiotracer to measure self-diffusion coefficients for Ca in a sodium-calcium-aluminosilicate melt (29% Na₂O, 5% CaO, 10% Al₂O₃, 56% SiO₂) at temperatures in the range 1100–1400°C and pressures to 30 kbar. Calcium diffusivity (D_Ca) was found to depend upon both temperature and pressure in a complex but systematic manner: $\text{(}\text{?D}{}_{\mathrm{Ca}}\text{?P}\text{)}{}_{\mathrm{T}}$ is always negative and has a larger absolute value at lower temperatures; $\text{(}\text{?D}{}_{\mathrm{Ca}}\text{?T}\text{)}{}_{\mathrm{P}}$ is positive and increases with increasing pressure. The overall dependence of D_Ca upon T and P is given approximately by $\text{D}{}_{\mathrm{caT.P}}\text{= [0.0025 exp(}\text{-23,107}\text{RT}\text{)] exp [P}\text{(0.7297T ? 1261.32)}\text{RT}\text{]}$. When expressed in terms of volume (V_a) and energy (E) of activation, the results are as follows: V_a ranges from 2.2 cm³/mole at 1400°C to 11.9 cm³/mole at 1100°C. and E ranges from 25.4 kcal/mole (1 kban to 49.8 kcal/mole (20 kbar).From the systematic dependence of D_Ca upon T and P, it is concluded that diffusion of Ca²⁺ in silicate melts does not take place by means of a vacant site mechanism, but is controlled instead by the amount and distribution of free volume in the melt structure.If it is assumed that the viscosity of the melt used in this study decreases with increasing pressure (Kushiro, 1976, J. Geophys. Res.81, 6351–6356) as D_Ca does, then the Stokes-Einstein inverse relation between viscosity and diffusivity is clearly violated, and its validity for silicate melts must be questioned. Thus, it appears that in silicate melts, unlike many liquids, viscous flow and diffusion are fundamentally different transport processes, involving different structural units.The effect of pressure on calcium diffusion is too small to invalidate kinetic models of upper mantle processes that have been based upon diffusivity values measured at 1 atm. Pressure may, however, induce significant reductions in the diffusion rates of large ions such as Rb⁺ or SiO^4?₄ in silicate melts.     

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.     

Serpentinization of cumulate ultramafic rocks from the North Arm Mountain massif of the Bay of Islands ophiolite

Partially serpentinized dunites and wehrlites comprise the bulk of the cumulate ultramafic unit at the North Arm Mountain massif of the Bay of Islands ophiolite complex, Newfoundland. In a suite of 59 dunites and werhlites from the base of the unit, the serpentinized portions consist of lizardite + chrysotile + brucite + (accessory) magnetite. The ratio of $\text{(}\text{lizardite}\text{+}\text{chrysotile}\text{)}\text{to brucite}\text{= ~8:2 (}\text{weight percent}\text{)}$. Petrographic observations show that most serpentinization occurred at the expense of olivine; only limited amounts of clinopyroxene were serpentized. An estimated volume increase of 32% accompanied serpentinization of the peridotites. Reconstructions of the primary modal proportions of wehrlites (made taking this volume increase into account) contain an average of 6% more clinopyroxene and 6% less olivine than do modal reconstructions that ignore the volume increase. Mass balance calculations provide no clear evidence for appreciable metasomatism of Al₂O₃, CaO, FeO, MgO, or SiO₂ during Serpentinization. The presence of brucite, the evidence that most serpentinization occurred at the expense of olivine, and the lack of appreciable metasomatism, suggest that the primary reaction that controlled serpentinization of the peridotites is: $\text{2Mg}{}_2\text{SiO}{}_2\text{+ 3H}{}_2\text{O ? Mg}{}_3\text{Si}{}_2\text{O}{}_5\text{(OH)}{}_4\text{+ Mg(OH)}{}_2$. olivine added serpentine brucite     

The partial molal volume of silicic acid in 0.725 M NaCl at 1°C determined by the neutralization of Na2SiO3

The partial molal volume of silicic acid ($\text{V}\text{?}\text{(Si(OH)}{}_4\text{)}$) in 0.725 M NaCl at 1°C was calculated from the measured volume change (Δ$\text{V}\text{?}{}_{\mathrm{n}}$) due to the neutralization of anhydrous sodium metasilicate with HCl and the $\text{V}\text{?}\text{(HCl)}$ and $\text{V}\text{?}\text{(NaCl)}$ obtained from the literature. $\text{V}\text{?}\text{(Si(OH)}{}_4\text{) = 59.0 cm}{}^3\text{mol}\text{? 1}$, determined under experimental conditions of pH = 2.2, compares favorably with $\text{V}\text{?}\text{(Si(OH)}{}_4\text{) = 58.9 cm}{}^3\text{mol}{}^{\mathrm{?1}}$ calculated from the measured volume change due to the hydrolysis of the meta-silicate salt at pH = 11 and from the partial molal volume due to electrostriction ($\text{V}\text{?}{}_{\mathrm{elect}}$) of water by charged Si species present in the solution at the high pH. This agreement lends support to a semiempirical model for calculating $\text{V}\text{?}{}_{\mathrm{elect}}$ in developed by Millero (1969). $\text{V}\text{?}\text{(NaOH) = ? 5.45 cm}{}^3\text{mol}{}^{\mathrm{?1}}$ in 0.725 M NaCl needed for this calculation was also determined in this work. The rate of polymerization of Si(OH)₄ at 1°C was monitored to insure that the monomer Si(OH)₄ was the main Si species present during the determination of $\text{V}\text{?}\text{(Si(OH)}{}_4\text{)}$ by neutralization of the alkali silicate. $\text{V}\text{?}\text{(Si(OH)}{}_4\text{)}$ determined in this study compares favorably with the value calculated from high pressure solubility measurements.     

Lead orthophosphates—II. Stability of cholopyromophite at 25°c

The conversion of secondary lead orthophosphate [PbHPO₄] into chloropyromorphite [Pb₅(PO₄)₃Cl] in ca. 10^?1 M NaCl solutions has been investigated at 25°C. From the composition of the supernatant solutions, the solubility product constant for Pb₅(PO₄)₃Cl has been calculated to be 10^?84.4±0.1, corresponding to $\text{ΔG}{}_{\mathrm{?}}\text{°}$ of ?906.2 kcal mol^?1. The solution equilibria and phase relationships in the system PbCl₂-PbO-P₂O₈-H₂O are discussed along with the geological implications.     

Use of Pitzer's equations to determine the media effect on the formation of lead chloro complexes

The spectrophotometric measurements of chloro complexes of lead in aqueous HCl, NaCl, MgCl₂ and CaCl₂ solutions at 25°C have been analyzed using Pitzer's specific interaction equations. Parameters for activity coefficients of the complexes PbCl⁺, PbCl₂⁰ and PbCl₃^? have been determined for the various media. Values of $\text{K}{}_1\text{= 30.0 ± 0.6}$, $\text{K}{}_2\text{= 106.7 ± 2.1}$ and $\text{K}{}_3\text{= 73.0 ± 1.5}$ were obtained for the cumulative formation constants. [$\text{Pb}{}^{\mathrm{2+}}\text{+ nCl}{}^{\mathrm{?}}\text{→ PbCl}{}_{\mathrm{n}}{}^{\mathrm{2?n}}\text{)}$]. These values are in reasonable agreement with literature data. The Pitzer parameters for the PbCl ion pairs in various media were used to calculate the speciation of Pb²⁺ in an artificial seawater solution.     

Mineral-solution equilibria—III. The system Na2OAl2O3SiO2H2OHCl

Chemical equilibrium between sodium-aluminum silicate minerals and chloride bearing fluid has been experimentally determined in the range 500–700°C at 1 kbar, using rapid-quench hydrothermal methods and two modifications of the Ag + AgCl acid buffer technique. The temperature dependence of the thermodynamic equilibrium constant (K) for the reaction $\text{NaAlSi}{}_3\text{O}{}_8\text{+ HCl}{}^{\mathrm{o}}\text{= NaCl}{}^{\mathrm{o}}\text{+}\text{1}\text{2Al}{}_2\text{SiO}{}_5\text{, +}\text{5}\text{2SiO}{}_2\text{+}\text{1}\text{2H}{}_2\text{O}$ Albite Andalusite Qtz. can be described by the following equation: log k = ?4.437 + 5205.6/T(K) The data from this study are consistent with experimental results reported by Montoya and Hemley (1975) for lower temperature equilibria defined by the assemblages albite + paragonite + quartz + fluid and paragonite + andalusite + quartz + fluid. Values of the equilibrium constants for the above reactions were used to estimate the difference in Gibbs free energy of formation between NaCl^o and HCl^o in the range 400–700°C and 1–2 kbar. Similar calculations using data from phase equilibrium studies reported in the literature were made to determine the difference in Gibbs free energy of formation between KCl^o and HCl^o. These data permit modelling of the chemical interaction between muscovite + kspar + paragonite + albite + quartz assemblages and chloride-bearing hydrothermal fluids.     

Rutile solubility and titanium coordination in silicate melts

The solubility of rutile has been determined in a series of compositions in the K₂O-Al₂O₃-SiO₂ system ($\text{K}{}^1\text{=}\text{K}{}_2\text{O}\text{(K}{}_2\text{O}$ + Al₂O₃) = 0.38–0.90), and the CaO-Al₂O₃-SiO₂ system ($\text{C}{}^1\text{=}\text{CaO}\text{(CaO + Al}{}_2\text{O}{}_3\text{)}\text{= 0.47–0.59}$). Isothermal results in the KAS system at 1325°C, 1400°C, and 1475°C show rutile solubility to be a strong function of the $\text{K}{}^1$ ratio. For example, at 1475°C the amount of TiO₂ required for rutile saturation varies from 9.5 wt% ($\text{K}{}^1\text{= 0.38}$) to 11.5 wt% ($\text{K}{}^1\text{= 0.48}$) to 41.2 wt% ($\text{K}{}^1\text{= 0.90}$). In the CAS system at 1475°C, rutile solubility is not a strong function of $\text{C}{}^1$. The amount of TiO₂ required for saturation varies from 14 wt% ($\text{C}{}^1\text{= 0.48}$) to 16.2 wt% ($\text{C}{}^1\text{= 0.59}$).The solubility changes in KAS melts are interpreted to be due to the formation of strong complexes between Ti and K⁺ in excess of that needed to charge balance Al³⁺. The suggested stoichiometry of this complex is K₂Ti₂O₅ or K₂Ti₃O₇. In CAS melts, the data suggest that Ca²⁺ in excess of A1³⁺ is not as effective at complexing with Ti as is K⁺. The greater solubility of rutile in CAS melts when $\text{C}{}^1$ is less than 0.54 compared to KAS melts of equal $\text{K}{}^1$ ratio results primarily from competition between Ti and Al for complexing cations (Ca vs. K).TiK_β x-ray emission spectra of KAS glasses ($\text{K}{}^1\text{= 0.43–0.60}$) with 7 mole% added TiO₂, rutile, and Ba₂TiO₄, demonstrate that the average Ti-O bond length in these glasses is equal to that of rutile rather than Ba₂TiO₄, implying that Ti in these compositions is 6-fold rather than 4-fold coordinated. Re-examination of published spectroscopic data in light of these results and the solubility data, suggests that the 6-fold coordination polyhedron of Ti is highly distorted, with at least one Ti-O bond grossly undersatisfied in terms of Pauling's rules.     

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.     

Viscosity and configurational entropy of silicate melts

With the configurational entropy theory of relaxation processes of Adam and Gibbs (1965), one predicts that the viscosity depends on temperature according to log $\text{η = A}{}_{\mathrm{e}}\text{+}\text{B}{}_{\mathrm{e}}\text{TS}{}_{\mathrm{<mtext>conf</mtext>}}$, where S_conf is the configurational entropy of the liquid. Thermochemical calculations of S_conf performed for some mineral compositions show the importance of non-configurational contributions to the entropy differences between amorphous and crystalline phases. Except for the case of SiO₂, the available thermodynamic data indicate that the above equation for viscosity accounts quantitatively for the experimentally determined temperature dependence of the viscosity of silicate melts. The Adam and Gibbs theory also provides a simple rationale for the non linear variation of the logarithmic viscosity with composition in mixed alkali silicate liquids at low temperatures, the minimum of viscosity resulting from the contribution of the entropy of mixing to S_conf.     

Etude des inclusions fluides du système N2-CO2 de dolomites et de quartz de Tunisie septentrionale. Données de la microcryoscopie et de l'analyse à la microsonde à effet Raman

Optical and analytical studies were performed on 400 N₂ + CO₂ gas bearing inclusions in dolomites and quartz from Triassic outcrops in northern Tunisia. Other fluids present include brines (NaCl and KCl bearing inclusions) and rare liquid hydrocarbons. At the time of trapping, such fluids were heterogeneous gas + brine mixtures. In hydrocarbon free inclusions the $\text{N}{}_2\text{(N}{}_2\text{+ CO}{}_2\text{)}$ mole ratio was determined using two different non-destructive and punctual techniques: Raman microprobe analysis, and optical estimation of the volume ratios of the different phases selected at low temperatures. In the observed range of compositions, the two methods agree reasonably well.The N₂ + CO₂ inclusions are divided into three classes of composition: (a) $\text{N}{}_2\text{(N}{}_2\text{+ CO}{}_2\text{)}\text{> 0,57}$: Liquid nitrogen is always visible at very low temperature and homogenisation occurs in the range ?151°C to ? 147°C (nitrogen critical temperature) dry ice (solid CO₂) sublimates between ?75°C and ?60°C; (b) $\text{0,20 <}\text{N}{}_2\text{(N}{}_2\text{+ CO}{}_2\text{)}\text{? 0,57}$: liquid nitrogen is visible at very low temperature but dry ice melts on heating; liquid and gas CO₂ homogenise to liquid phase between ?51°C to ?22°C; (c) $\text{N}{}_2\text{(N}{}_2\text{+ CO}{}_2\text{)}\text{? 0,20}$: liquid nitrogen is not visible even at very low temperature (?195°C) and liquid and gas CO₂ homogenise to liquid phase between ?22°C and ?15°C. The observed phases changes are used to propose a preliminary phase diagram for the system CO₂-N₂ at low temperatures.Assuming additivity of partial pressures, isochores for the CO₂-N₂ inclusions have been computed. The intersection of these isochores with those for brine inclusions in the same samples may give the P and T of trapping of the fluids.