Dynamics of Na in sodium aluminosilicate glasses and liquids

²³Na NMR measurements on Na₂Si₃O₇, Na₃AlSi₆O₁₅, and NaAlSi₃O₈ glasses from room temperature to 1200°C show that the dynamics and local structure of sodium in silicate/aluminosilicate glasses and melts vary with composition and temperature.The peak positions decrease in frequency between room temperature and 200°C indicating that the Na sees a larger average site as temperature is increased. Between 200°–300° and 700°C, line widths, nutation frequencies and peak positions are consistent with motional averaging of quadrupolar satellites. Above 700°C there is little or no change in the peak positions with temperature. Chemical shifts of the materials at 1000°C (Na₂Si₃O₇: 3.6; Na₃AlSi₆O₁₅:-1.3; NaAlSi₃O₈:-6.4 ppm) indicate a slight change in the average Na coordination number from 6–7 for the silicate to 7–8 for the aluminosilicates.     

Surface tension of melts in the system CaMgSi2O6-CaAl2Si2O8 and its structural significance

The surface tension between silicate melts and air has been measured for melt compositions lying on the diopside-anorthite (Di-An) join from 1300° C to 1580° C. It ranges from 300 dyne/cm to 400 dyne/cm, and decreases with increasing temperature, except for a pure diopside composition. At relatively high temperatures, the surface tension decreases as the anorthite content increases, whereas at lower temperature it is almost constant. These results suggest that diopside melt has a more discrete structure at higher temperatures, whereas, anorthite-bearing melts do not dissociate in the temperature range studied. They also suggest that the structure of both the surface and interior parts of the melt are almost identical at lower temperatures, but at higher temperatures, the surface part has a more polymerized structure with Al₂O₃ enrichment. The surface energy, obtained from the relationship between surface tension and temperature, increases from 294 erg/cm² (Di composition) to 1013 erg/cm² (Di40An60) with increasing anorthite content.     

Viscosity of silicate melts in CaMgSi2O6–NaAlSi2O6 system at high pressure

In situ X-ray viscometry of the silicate melts was carried out at high pressure and at high temperature. The viscosity of the silicate melts in the diopside(Di)–jadeite(Jd) system was determined in the pressure range from 1.88 GPa to 7.9 GPa and in the temperature range from 2,003 K to 2,173 K. The viscosity of the Di 25%–Jd 75% melt decreases continuously to 5.0 GPa, whereas the viscosity of the Di 50%–Jd 50% melt increases over 3.5 GPa. The viscosity of the Di50%–Jd 50% melt reaches a minimum around 3.5 GPa. Since the amounts of silicon in the two melts are the same, the difference in the pressure dependence of the viscosity may be controlled by another network-forming element, i.e., aluminum. The difference in the pressure dependence of the viscosities in the melts with two intermediate compositions in the Di–Jd system is estimated to be due to the difference in the melt structures at high pressures and high temperatures.     

Liquidus karrooite stability and composition at reducing conditions

We have determined the stability of rutile and karrooite on the liquidus of pseudobinary silicate melts of anorthite–diopside eutectic composition with, in addition, of up to 32 wt% TiO₂ at one atm total pressure and at wide range of oxygen fugacities. At 1,300°C and at an fO₂ below 10^–11.2 atm rutile (TiO₂) is replaced as liquidus phase by the pseudobrookite-type MgTi₂O₅–Ti₃O₅ solid solution with some Al in the crystal structure. The composition of karrooite was found to be strongly dependent on oxygen fugacity. Crystalline phases were identified by transmission electron microscopy (TEM). The results obtained in this study are relevant for understanding the chemistry of lunar armalcolites. Using excess rutile solubility at reducing conditions, we estimated the Ti⁴⁺ /Ti³⁺ ratio in silicate melts at 1,300°C as function of fO₂.     

Rare earth partitioning between immiscible carbonate and silicate liquids and CO2 vapor: Results and implications for the formation of light rare earth-enriched rocks

Melting relations at 5 and 20 kbar on the composition join sanidine-potassium carbonate are dominated by a two-liquid region that covers over 60% of the join at 1,300 ° C. At this temperature, the silicate melt contains approximately 19 wt% carbonate component at 5 kbar and 32 wt% carbonate component at 20 kbar. The conjugate carbonate melt contains less than 5 wt% silicate component, and it varies less as a function of temperature than does the silicate melt.Partition coefficients for Ce, Sm, and Tm between the immiscible carbonate and silicate melts at 1,200 ° and 1,300 ° C at 5 and 20 kbar are in favor of the carbonate melt by a factor of 2–3 for light REE and 5–8 for heavy REE. The effect of pressure on partitioning cannot be evaluated independently because of complementary changes in melt compositions.Minimum REE partition coefficients for CO₂ vapor/carbonate melt and CO₂ vapor/silicate melt can be calculated from the carbonate melt/silicate melt partition coefficients, the known proportions of melt, and maximum estimates of the proportion of CO₂ vapor. The vapor phase is enriched in light REE relative to both melts at 20 kbar and enriched in all REE, especially the light elements, at 5 kbar. The enrichment of REE in CO₂ vapor relative to both melts is 3–4 orders of magnitude in excess of that in water vapor (Mysen, 1979) at 5 kbar and is approximately the same as that in water vapor at 20 kbar.Mantle metasomatism by a CO₂-rich vapor enriched in light REE, occurring as a precursor to magma genesis, may explain the enhanced REE contents and light REE enrichment of carbonatites, alkali-rich silicate melts, and kimberlites. Light REE enrichment in fenites and the granular suite of nodules from kimberlites attests to the mobility of REE in CO₂-rich fluids under both mantle and crustal conditions.     

The onset of non-Newtonian rheology of silicate melts

The viscoelastic behavior of silicate melts has been measured for a range of compositions (NaAlSi₃O₈, NaCaAlSi₂O₇, CaMgSi₂O₆, Li₂Si₄O₉, Na₂Si₄O₉, K₂Si₄O₉, Na₂Si₃O₇, K₂Si₃O₇ and Na₂Si₂O₅) using the fiber elongation method. A1l compositions exhibit Newtonian behavior at low strain-rates, but non-Newtonian behavior at higher strain-rates, with strain-rate increasing faster than the applied stress. The decrease in shear viscosity observed at the high strain-rates ranges from 0.3 to 1.6 log₁₀ units (Pa s). The relaxation strain-rates, _relax, of these melts have been estimated from the low strain-rate, Newtonian, shear viscosity, using the Maxwell relationship; _relax= ^–1=(_s/G)^–1. For all compositions investigated, the onset of non-Newtonian rheology is observed at strain-rates 2.5+0.5 orders of magnitude less than the calculated relaxation strain-rate. This difference between the non-Newtonian onset and the relaxation strain-rate is larger than that predicted by the single relaxation time Maxwell model. Normalization of the experimental strain-rates to the relaxation strain-rate predicted from the Maxwell relation, eliminates the composition. and temperature-dependence of the onset of non-Newtonian behavior. The distribution of relaxation in the viscoelastic region appears to be unrelated to melt chemistry. This conclusion is consistent with the torsional, frequency domain study of Mills (1974) which illustrated a composition-invariance of the distribution of the imaginary component of the shear modulus in melts on the Na₂O-SiO₂ join. The present, time domain study of viscoelasticity contrasts with frequency domain studies in terms of the absolute strains employed. The present study employs relatively large total strains (up to 2). This compares with typical strains of 10^–8 in ultrasonic (frequency domain) studies. The stresses used to achieve the strain-rates required to observe viscoelastic behavior in this study approach the tensile strength of the fibers with the result that some of our experiments were terminated by fiber breakage. Although the breakage is unrelated to the observation of non-Newtonian viscosity, their close proximity in this and earlier studies suggests that brittle failure of igneous melts, may, in general, be preceded by a period of non-Newtonian rheology.     

Oxygen diffusion in basalt and andesite melts: experimental results and discussion of chemical versus tracer diffusion

Chemical diffusion coefficients for oxygen in melts of Columbia River basalt (Ice Harbor Dam flow) and Mt. Hood andesite have been determined at 1 atm. The diffusion model is that of sorption or desorption of oxygen into a sphere of uniform initial concentration from a constant and semi-infinite atmosphere. The experimental design utilizes a thermogravimetric balance to monitor the rate of weight change arising from the response of the sample redox state to an imposed fO₂. Oxygen diffusion coefficients are approximately an order-ofmagnitude greater for basaltic melt than for andesitic melt. At 1260° C, the oxygen diffusion coefficients are: D=1.65×10^–6cm²/s and D=1.43×10^–7cm²/s for the basalt and andesite melts, respectively. The high oxygen diffusivity in basaltic melt correlates with a high ratio of nonbridging oxygen/tetrahedrally coordinated cations, low melt viscosity, and high contents of network-modifying cations. The dependence of the oxygen diffusion coefficient on temperature is: D=36.4exp(–51,600±3200/RT)cm²/s for the basalt and D=52.5exp(–60,060±4900/RT)cm²/s for the andesite (R in cal/deg-mol; T in Kelvin). Diffusion coefficients are independent of the direction of oxygen diffusion (equilibrium can be approached from extremely oxidizing or reducing conditions) and thus, melt redox state. Characteristic diffusion distances for oxygen at 1260° C vary from 10^-2 to 10² m over the time interval of 1 to 10⁶ years. A compensation diagram shows two distinct trends for oxygen chemical diffusion and oxygen tracer diffusion. These different linear relationships are interpreted as supporting distinct oxygen transport mechanisms. Because oxygen chemical diffusivities are generally greater than tracer diffusivities and their Arrhenius activation energies are less, transport mechanisms involving either molecular oxygen or vacancy diffusion are favored.     

Carbon isotopic fractionation between CO2 vapour,silicate and carbonate melts: an experimental study to 30 kbar

The carbon isotopic fractionation between CO₂ vapour and sodamelilite (NaCaAlSi₂O₇) melt over a range of pressures and temperatures has been investigated using solid-media piston-cylinder high pressure apparatus. Ag₂C₂O₄ was the source of CO₂ and experimental oxygen fugacity was buffered at hematite-magnetite by the double capsule technique. The abundance and isotopic composition of carbon dissolved in sodamelilite (SM) glass were determined by stepped heating and the ¹³C of coexisting vapour was determined directly by capsule piercing. CO₂ solubility in SM displays a complex behavior with temperature. At pressures up to 10 kbars CO₂ dissolves in SM to form carbonate ion complexes and the solubility data suggest slight negative temperature dependence. Above 20 kbars CO₂ reacts with SM to form immiscible Na-rich silicate and Ca-rich carbonate melts and CO₂ solubility in Na-enriched silicate melt rises with increasing temperature above the liquidus. Measured values for carbon isotopic fractionation between CO₂ vapour and carbonate ions dissoived in sodamelilite melt at 1200°–1400° C and 5–30 kbars average 2.4±0.2, favouring¹³C enrichment in CO₂ vapour. The results are maxima and are independent of pressure and temperature. Similar values of 2 are obtained for the carbon isotopic fractionation between CO₂ vapour and carbonate melts at 1300°–1400° C and 20–30 kbars.     

Molecular dynamics simulations of pressure and temperature effects on MgSiO3 and Mg2SiO4 melts and glasses

Ab-initio interionic potentials for Mg²⁺, Si⁴⁺, and O^2– have been used in molecular dynamics (MD) simulations to investigate diffusivity changes, pressure-induced structural transitions, and temperature effects on polymerization in MgSiO₃ and Mg₂SiO₄ melts and glasses. The potential gives reasonable agreement with the 0.1 MPa radial distribution function of MgSiO₃ glass. Maxima in the diffusion coefficients of Si⁴⁺ and O^2– occur as pressure is increased on the MgSiO₃ melt. The controlling structural mechanism for this behavior is the Q¹ species of SiO₄ tetrahedra. Mg²⁺ diffusion coefficients decrease monotonically with pressure in both melt compositions. Increasing Mg²⁺ coordination number and population of 3- and 4-membered SiO₄ rings with pressure combine to hinder translation of the Mg²⁺ ions. The dominant changes in structure with pressure are a decrease in the intertetrahedral (Si-O--Si) angle up to approximately 4 g/cm³ and coordination changes of the ions above this density. Temperature effects on viscosity in these simulated melts are indirectly studied by analyzing polymerization changes with temperature. Polymerization and coordination numbers increase with decreasing temperature and a small quench rate effect is observed. Fair agreement is found between the MD simulations and experimental equation of state for Mg₂SiO₄, but the equation of state predictions for MgSiO₃ melts are much less accurate. The zero pressure volume, V ₀, is significantly higher and K ₀ is lower in the simulations than empirical values. The inadequacies reflect error in using the ionic approximation for polymerized systems and a need to collect more data for a variety of molecular configurations in the development of ab-initio potentials.     

Solution of H2O in diopside melts: a thermodynamic model

Structural similarities between dry diopside melt and superhydrous albite melt (X _w >0.5) — both lack three-dimensional silicate units — suggest that thermodynamic relations may be similar. A model based on that assumption successfully predicts diopside melting relations and H₂O solubilities. For the model, the three partial differential equations describing solution of H₂O in albite melt for X _w >0.5 have been integrated for diopside melt from X _w =0 to X _w at least as large as 0.76, with two exceptions: an alternative partial differential equation for Henrian solution of H₂O in dilute melts was applied for X _w <0.20, and an alternative differential equation for the pressure dependence of a _w at pressures below 2 kbar was developed. The latter alternative equation yields relatively small ¯V_w's at low pressures rather than the large ¯V_w's calculated from the equation from the albite system. Available experimental solubility data are not precise enough to offer a choice between the small-¯V_w and large-¯V_w equations. Integration of all the partial differential equations was constrained solely by the P and T of a single experimentally-determined point on the H₂O-saturated solidus.Solubilities calculated by a Henrian-analogue solution model (a _di=X _di ² ) from the experimental H₂O saturated solidus lie outside experimental solubility constraints for dilute melts. On the other hand, a Henrian model (a _di=X_di) successfully predicts solubilities in dilute melts. The formulation of the Henrian model and magnitudes of model molar entropies of solution are consistent with the hypothesis that H₂O dissolves in diopside melt as an essentially undissociated species with little ordering on melt structural sites. That species could in turn be consistently, if not uniquely, interpreted to be molecular H₂O or a hydroxylation (OH^–) complex formed from nonbridging oxygens.     

Concentrations and behavior of oxygen and oxide ion in melts of composition CaO · MgO · xSiO2

The concentrations and behavior of oxygen and oxide ion were studied in silicate melts of composition CaO · MgO · xSiO₂ (1.25 ≤ x ≤ 3) in the temperature range 1425 to 1575°C by cyclic voltammetry and chronopotentiometry. Electroreduction of oxygen is a reversible, 2 electron process involving dissociated oxygen atoms. The Henry's Law constant for O₂ in molten diopside (CaO · MgO · 2SiO₂) is 0.023 ± 0.004 mole/l atm at 1450°C. The diffusion coefficient for molecular oxygen in diopside melt is 4.5 ± .5 × 10^?6 cm²/sec at 1450°C and the activation energy of diffusion is 80 ± 2 kcal/mole. Oxide ions produced by electroreduction of oxygen, rapidly dissociate silicate polymers, causing the concentration of free oxide ions in diopside melt to be buffered at a low level (4.7 ± .8 × 10^?5 mole/l). The concentration of free oxide ion increases at higher proportions of metal oxides but remains at this value in more silicic melts. The rate of formation of oxide ions by polymerization in diopside melt is 0.021 ± .007 mole/l sec. Thermodynamic parameters (the standard free energy, enthalpy and entropy) for the oxidation of Ni, Co, and Zn in diopside melt in equilibrium with gaseous oxygen agree with those for solid oxide systems. The platinum reference electrode in molten diopside is a reversible, oxygen electrode.     

Experimental determination of oxygen isotope fractionations between CO2 vapor and soda-melilite melt

We report results of experiments constraining oxygen isotope fractionations between CO₂ vapor and Na-rich melilitic melt at 1 bar and 1250 and 1400°C. The fractionation factor constrained by bracketed experiments, 1000^.lnα_{CO2-Na melilitic melt}, is 2.65±0.25 ‰ (±2σ; n=92) at 1250°C and 2.16±0.16 ‰ (2σ; n=16) at 1400°C. These values are independent of Na content over the range investigated (7.5 to 13.0 wt. % Na₂O). We combine these data with the known reduced partition function ratio of CO₂ to obtain an equation describing the reduced partition function ratio of Na-rich melilite melt as a function of temperature. We also fit previously measured CO₂-melt or -glass fractionations to obtain temperature-dependent reduced partition function ratios for all experimentally studied melts and glasses (including silica, rhyolite, albite, anorthite, Na-rich melilite, and basalt). The systematics of these data suggest that reduced partition function ratios of silicate melts can be approximated either by using the Garlick index (a measure of the polymerization of the melt) or by describing melts as mixtures of normative minerals or equivalent melt compositions. These systematics suggest oxygen isotope fractionation between basalt and olivine at 1300°C of approximately 0.4 to 0.5‰, consistent with most (but not all) basalt glass-olivine fractionations measured in terrestrial and lunar basalts.     

The role of phosphorus in rhyolitic liquids as determined from the homogeneous iron redox equilibrium

The redox ratio of iron is used as an indicator of solution properties of silicate liquids in the system (SiO–Al₂O₃–K₂O–FeO–Fe₂O₃–P₂O₅). Glasses containing 80–85 mol% SiO₂ with 1 mol% Fe₂O₃ and compositions covering a range of K₂O/Al₂O₃ were synthesized at 1400°C in air (fixed fO₂). Variations in the ratio FeO/FeO_1.5 resulting from the addition of P₂O₅ are used to determine the solution behavior of phosphorus and its interactions with other cations in the silicate melt. In 80 mol% SiO₂ peralkaline melts the redox ratio, expressed as FeO/FeO_1.5, is unchanged relative to the reference curve with the addition of 3 mol% P₂O₅. Yet, the iron redox ratio in the 85 mol% SiO₂ potassium aluminosilicate melts is decreased relative to phosphorus-free liquids even for small amounts of P₂O₅ (0.5 mol%). The redox ratio in peraluminous melts is decreased relative to phosphorus- free liquids at P₂O₅ concentrations of 3 mol%. In peraluminous liquids, complexing of both Fe⁺³–O–P⁺⁵ and Al⁺³–O–P⁺⁵ occur. The activity coefficient of Fe⁺³ is decreased because more ferric iron can be accommodated than in phosphorus-free liquids. In peralkaline melts, there is no evidence that P⁺⁵ is removing K⁺ from either Al⁺³ or Fe⁺³ species. In chargebalanced melts with 3 mol% Fe₂O₃ and very high P₂O₅ concentrations, phosphorus removes K⁺ from K–O–Fe⁺³ complexes resulting in a redox increase. P₂O₅ should be accommodated easily in peraluminous rhyolitic liquids and phosphate saturation may be suppressed relative to metaluminous rhyolites. In peralkaline melts, phosphate solubility may increase as a result of phosphorus complexing with alkalis. The complexing stoichiometry may be variable, however, and the relative influence of peralkalinity versus temperature on phosphate solubility in rhyolitic melts deserves greater attention.     

Experimental study on the phase equilibria in the join CaMgSi2O6-CaCrCrSiO6 with special reference to the blue diopside

The behaviour of tetrahedrally coordinated and octahedrally coordinated Cr³⁺ ions in diopside is discussed from the study on the join CaMg-Si₂O₆-CaCrCrSiO₆. The molecule CaCrCrSiO₆ decomposes into uvarovite+eskolaite and its maximum solubility in diopside is 6.7 wt percent at 940 ° C. Crystalline phases are diopside ss (ss is abbreviation of solid solution), uvarovite ss, wollastonite ss, spinel and eskolaite. The diopside ss is blue in colour. Its optical spectra were measured in the wavelenght range of 325–2600 nm, and assigned after tetrahedral configuration Td and octahedral configuration Oh. It is estimated that octahedral Cr³⁺ ions are in high spin state, while tetrahedral Cr³⁺ ions may be probably in low spin state. The _t and B are 10,300–10,370 cm^–1 and 429–432 cm^–1. The CFSE for tetrahedral low spin Cr³⁺ ions is nearly the same as that for octahedral high spin Cr³⁺ ions. The ionic radii of tetrahedral low spin Cr³⁺ ions calculated are 0.47–0.53 Å, shrinked from the radius of octahedral high spin Cr³⁺ ion (0.615 Å) as much as 14–24 percent. Petrologic implications of the result are also discussed.The first half of the D. Sc. dissertation of K. Ikeda presented to Hokkaido University in June, 1976     

Equilibrium density and expansivity of silicate melts in the glass transition range

Equilibrium volumes and expansivities of three liquids in the system anorthite (CaAl₂Si₂O₈)–diopside (CaMgSi₂O₆) have been derived from dilatometric measurements of the equilibrium length of samples in the glass transition range. The typical temperature range of 40 K for the measurements is limited at low temperature by the very long times necessary to reach structural equilibrium and at high temperature by the penetration of the rod used to measure sample dilatation. Despite such narrow intervals, the expansivities are determined to better than 3% thanks to the high precision with which length changes are measured. The coefficient of volume thermal expansion (1/V dV/dT) of the fully relaxed liquid just above the glass transition is found to decrease linearly from diopside composition (139 ± 4 × 10⁻⁶ K⁻¹) to anorthite composition (59 ± 2 × 10⁻⁶ K⁻¹). These values are greater than those determined for the same liquids at superliquidus temperatures, demonstrating that expansivities of silicate melts may decrease markedly with increasing temperature. A predictive model based upon partial molar volumes which vary as a linear function of the logarithm of temperature is proposed. Received: 25 February 2000 / Accepted: 29 May 2000     

Solution mechanisms of H2O in depolymerized peralkaline melts

Solubility mechanisms of water in depolymerized silicate melts quenched from high temperature (1000°-1300°C) at high pressure (0.8-2.0 GPa) have been examined in peralkaline melts in the system Na₂O-SiO₂-H₂O with Raman and NMR spectroscopy. The Na/Si ratio of the melts ranged from 0.25 to 1. Water contents were varied from ∼3 mol% and ∼40 mol% (based on O = 1). Solution of water results in melt depolymerization where the rate of depolymerization with water content, ∂(NBO/Si)/∂X_H2O, decreases with increasing total water content. At low water contents, the influence of H₂O on the melt structure resembles that of adding alkali oxide. In water-rich melts, alkali oxides are more efficient melt depolymerizers than water. In highly polymerized melts, Si-OH bonds are formed by water reacting with bridging oxygen in Q⁴-species to form Q³ and Q² species. In less polymerized melts, Si-OH bonds are formed when bridging oxygen in Q³-species react with water to form Q²-species. In addition, the presence of Na-OH complexes is inferred. Their importance appears to increase with Na/Si. This apparent increase in importance of Na-OH complexes with increasing Na/Si (which causes increasing degree of depolymerization of the anhydrous silicate melt) suggests that water is a less efficient depolymerizer of silicate melts, the more depolymerized the melt. This conclusion is consistent with recently published ¹H and ²⁹Si MAS NMR and ¹H-²⁹Si cross polarization NMR data.     

Simulation of the molecular mass distributions of Si anions on the liquidus of the Na<Subscript>2</Subscript>O-SiO<Subscript>2</Subscript> system with regard for the disproportionation of <Emphasis Type="Italic">Q</Emphasis> structons

A new version of the STRUCTON (2009) computer model is proposed for the simulation of the molecular mass distributions (MMD) characterizing the diversity of anions in silicate melts depending on their polymerization and temperature. In contrast to earlier versions, the new version of the model accounts for disproportionation reactions of Q ⁿ species and makes use of their proportions in the statistical simulations of the origin of real Si-O complexes. The new potentialities of the STRUCTON program package are illustrated by its application to studying the structural-chemical characteristics of melts in the Na₂O-SiO₂ system along its liquidus line, including the points of eutectics and phase transitions at 0.333 ≤ $ N_{SiO_2 } $ N_{SiO_2 } < 0.500. This problem is solved with the use of a temperature-composition dependence of polymerization constants K _p ^Na in the Toop-Samis approximation. The variations in K _p ^Na were proved to be as large as three orders of magnitude due to both the temperature effect at a constant composition and the composition effect at a constant temperature. The results of the MMD simulations on the liquidus show that the concentration of the SiO₄⁴⁻ ion strongly decreases, and the proportion of chain species increases compared to those at a stochastic distribution. The concentration of the Si₂O₇⁶⁻ anion reaches its maximum (∼42%) at 40 mol % in the liquid, i.e., the composition of Na₆Si₂O₇. At $ N_{SiO_2 } $ N_{SiO_2 } > 0.40, this ion dominates over the SiO₄⁴⁻ monomer. More silicic melts with $ N_{SiO_2 } $ N_{SiO_2 } ≥ 0.45, are dominated by (Si _n O_3n )³ⁿ⁻ ring species, and the concentrations of these species are related as (Si₃O₉)⁶⁻ > (Si₄O₁₂)⁸⁻ > (Si₅O₁₅)¹⁰⁻. The maximum concentration of these flat rings also occurs near the composition of stoichiometric metasilicate with Si/O = 0.333. The comparison of the dependence of the average size of anions i _av and the average number of their species on depolymerization indicates that a change in the proportion of Q ⁿ species in melt at decreasing temperature results in structural restyling and an increase in the average size of Si-O complexes. The average number of anion species thereby decreases compared to that in a stochastic MMD. The results presented in this publication direct the progress in the thermodynamic theory of silicate melts to a new avenue that makes use of the capabilities and advantages of the ion-polymer model, the theory of associated solutions, spectroscopic data, and the experimental study of variations in oxide activities depending on composition and temperature.     

Water Content of Basalt Erupted on the ocean floor

Deep sea pillow basalts dredged from the ocean floor show that vesicularity changes with composition as well as with depth. Alkalic basalts are more vesicular than tholeiitic basalts erupted at the same depth. The vesicularity data, when related to experimentally determined solubility of water in basalt, indicate that K-poor oceanic tholeiites originally contained about 0.25 percent water, Hawaiian tholeiites of intermediate K-content, about 0.5 percent water, and alkali-rich basalts, about 0.9 percent water. Analyses of fresh basalt pillows show a systematic increase of H₂O⁺ as the rocks become more alkalic. K-poor oceanic tholeiites contain 0.06–0.42 percent H₂O⁺, Hawaiian tholeiites, 0.31–0.60 percent H₂O⁺, and alkali rich basalts 0.49–0.98 percent H₂O⁺. The contents of K₂O, P₂O₅, F, and Cl increase directly with an increase in H₂O⁺ content such that at 1.0 weight percent H₂O⁺, K₂O is 1.58 percent, P₂O₅ is 0.55 percent, F is 0.07 percent, and Cl is 0.1 percent. The measured weight percent of deuterium on the rim of one Hawaiian pillow is –6.0 (relative to SMOW); this value, which is similar to other indications of magmatic water, suggests that no appreciable sea water was absorbed by the pillow during or subsequent to eruption on the ocean floor.Concentrations of volatile constituents in the alkali basalt melts relative to tholeiitic melts can be explained by varying degrees of partial melting of mantle material or by fractional crystallization of a magma batch.Publication authorized by the Director, U.S. Geological Survey.     

Structure and Composition Effects on the Oxygen Isotope Fractionation in Silicate Melts

The influence of melt composition and structure on the oxygen isotope fractionation was studied for the multicomponent (SiO₂ ± TiO₂ + Al₂O₃ ± Fe₂O₃ + MgO ± CaO) system at 1500°C and 1 atm. The experiments show that significant oxygen isotope effects can be observed in silicate melts even at such high temperature. It is shown that the ability of silicate melt to concentrate ¹⁸O isotope is mainly determined by its structure. In particular, an increase of the NBO/T ratio in the experimental glasses from 0.11 to 1.34 is accompanied by a systematic change of oxygen isotope difference between melt and internal standard by values from–0.85 to +1.29‰. The obtained data are described by the model based on mass-balance equations and the inferred existence of O⁰, O^–, and O^2– (bridging, non-bridging, and free oxygen) ions in the melts. An application of the model requires the intra-structure isotope fractionation between bridging and non-bridging oxygens. Calculations show that the intra-structure isotope fractionation in our experiments is equal to 4.2 ± 1.0‰. To describe the obtained oxygen isotope effects at the melts relatively to temperature and fraction of non-bridging oxygen a general equation was proposed.     

Structure and properties of fluorine-bearing aluminosilicate melts: the system Na2O-Al2O3-SiO2-F at 1 atm

The chemical interaction between fluorine and highly polymerized sodium aluminosilicate melts [Al/(Al+Si)= 0.125–0.250 on the join NaAlO₂-SiO₂] has been studied with Raman spectroscopy. Fluorine is dissolved to form F^– ions that are electrically neutralized with Na⁺ or Al³⁺. There is no evidence for association of fluorine with either Si⁴⁺ or Al³⁺ in four-fold coordination and no evidence of fluorine in six-fold coordination with Si⁴⁺ in these melt compositions. Upon solution of fluorine nonbridging oxygens are formed and are a part of structural units with nonbridging oxygen per tetrahedral cations (NBO/T) about 2 and 1. The proportions of these two depolymerized units in the melts increase systematically with increasing F/(F+O) at constant Al/(Al+Si) and with decreasing Al/(Al+Si) at constant F/(F+O). Depolymerization (increasing NBO/T) of silicate melts results from a fraction of aluminum and alkalies (in the present study; Na⁺) reacting to form fluoride complexes. In this process an equivalent amount of Na⁺ (orginally required for Al-³⁺charge-balance) or Al³⁺ (originally required Na⁺ to exist in tetrahedral coordination) become network-modifiers.The structural data have been used to develop a method for calculating the viscosity of fluorine-bearing sodium aluminosilicate melts at 1 atm. Where experimental viscosity data are available, the calculated and measured values are within 5% of each other.A method is also suggested by which the liquidus phase equilibria of fluorine-bearing aluminosilicate melts may be predicted. In accord with published experimental data it is suggested, for example, that — on the basis of the determined solubility mechanism of fluorine in aluminosilicate melts — with increasing fluorine content of feldspar-quartz systems, the liquidus boundaries between aluminosilicate minerals (e.g., feldspars) and quartz shift away from silica.