Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geophysical”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Precipitation of low-temperature disordered dolomite induced by extracellular polymeric substances of methanogenic Archaea Methanosarcina barkeri : Implications for sedimentary dolomite formation

Abstract A correlation between methanogenesis and dolomite formation has been reported; however, the mechanism underlying this association is not fully understood. In this study, we conducted forced carbonate precipitation experiments at room temperature in calcite-seeded Ca/Mg carbonate solutions containing either purified non-living biomass or bound extracellular polymeric substances (EPS) of the methanogen Methanosarcina barkeri. Purified non-living biomass and bound EPS was used so as to avoid the possible influence of the complex components of the growing microbial culture on carbonate crystallization. Our results demonstrated that non-living biomass of M. Barkeri can enhance the Mg incorporation into calcitic structure and induce the crystallization of disordered dolomite. In the presence of ~113 mg L–1 of non-living biomass, disordered dolomite with ~41 and 45 mol% of MgCO3 was precipitated in solutions with initial Mg:Ca ratios of 5:1 and 8:1, respectively. A systematic increase in the MgCO3 contents of the precipitated Ca-Mg carbonates was also observed with the increased non-living biomass concentration. Bound EPS was shown to be the component of non-living biomass that catalyzed the precipitation of disordered dolomite. At only ~25 mg L–1 of bound EPS, disordered dolomite with ~47 and 48 mol% of MgCO3 was precipitated in solutions with initial Mg:Ca ratios of 5:1 and 8:1, respectively. We propose that adsorption of bound EPS to growing carbonate surfaces through hydrogen bonding is the key to catalyzing disordered dolomite crystallization, and that this mechanism is also applicable to natural EPS-induced dolomite formation. This study provides significant insight into the formation mechanism of microbial-induced dolomite with high δ13C values.

Geochemistry & Geophysics↗
DOE OSTI2023

Thermodynamic characterization of synthetic lead-arsenate apatites with different halogen substitutions

Abstract Thermodynamic parameters have been measured for synthetic analogs of the mimetite-group minerals Pb5(AsO4)3X (X = OH, Cl, Br, I) belonging to the apatite supergroup. Phases precipitated from aqueous solutions under ambient conditions with well characterized structures and compositions were studied. For each phase, dissolution enthalpy was experimentally determined by oxide melt drop solution calorimetry in a molten solvent of sodium molybdate (3Na2O·4MoO3) at 976 K. The enthalpy of formation from the elements ΔHf,elo was calculated using thermochemical cycles and was −3030.6 ± 11.5, −3026.6 ± 15.8, −2967.6 ± 25.0, and −2993.1 ± 12.2 kJ/mol for Pb5.00(AsO4)3.00OH0.86(CO3)0.07, Pb5.00(AsO4)3.00Cl0.80(CO3)0.10, Pb5.00(AsO4)3.00Br0.80(CO3)0.10, and Pb5.00(AsO4)3.00I0.45OH0.35(CO3)0.10, respectively. These ΔHf,elo values exhibit typical trends for apatites: they increased (were less negative) with the increasing molar mass and ionic radius of X and decreased with the electronegativity and ionization energy of X. The compilation and comparison of data for Ca-, Pb-, P-, and As-apatites revealed correlations indicating that thermodynamic enthalpic stability is largely influenced by chemical factors (e.g., differences in electronegativities of the elements, ionization energy, or ionic characteristics of the bonds) and to a lesser extent by physical and geometric parameters in the crystal structure related to the mass and size of the X anion. Using the correlations, it was possible to estimate the value of hitherto unknown ΔHf,elo for Pb5(AsO4)3F, −3144.3 ± 66.5 kJ/mol. The observed relationships apply to the entire apatite supergroup and can be used to predict the values of ΔHf,elo for phases that have not been studied experimentally. The new data on environmentally significant phases will contribute to the modeling of mineral-water interactions, particularly for potential use in the remediation of soils and wastes contaminated with Pb and As and in the immobilization of radioactive waste containing I-129.

Geochemistry & Geophysics↗
DOE OSTI2021

Seaborgite, LiNa6K2(UO2)(SO4)5(SO3OH)(H2O), the First Uranyl Mineral Containing Lithium

Abstract Seaborgite (IMA2019-087), LiNa6K2(UO2)(SO4)5(SO3OH)(H2O), is a new mineral species from the Blue Lizard mine, Red Canyon, San Juan County, Utah, U.S.A. It is a secondary phase found on gypsum in association with copiapite, ferrinatrite, ivsite, metavoltine, and römerite. Seaborgite occurs in sprays of light-yellow, long flattened prisms or blades, up to about 0.2 mm in length. Crystals are elongated on [100], flattened on {010}, and exhibit the forms {100}, {010}, {001}, and {101}. The mineral is transparent with vitreous luster and very pale-yellow streak. It exhibits bright lime-green fluorescence under a 405 nm laser. The Mohs hardness is ~2½. The mineral has brittle tenacity, curved or conchoidal fracture, and one good cleavage on {100}. The measured density is 2.97(2) g/cm3. The mineral is immediately soluble in H2O at room temperature. The mineral is optically biaxial (–), α = 1.505(2), β = 1.522(2), γ = 1.536(2) (white light); 2Vmeas = 85(1)°; moderate r < ν dispersion; orientation X ^ a ≈ 10°; pleochroic X colorless, Y and Z light green-yellow; X < Y ≈ Z. EPMA and LA-ICP-MS analyses of seaborgite undermeasured its Li, K, and Na. The empirical formula using Li, Na, and K based on the structure refinement is Li1.00Na5.81K2.19(UO2)(SO4)5(SO3OH)(H2O). Seaborgite is triclinic, P1, a = 5.4511(4), b = 14.4870(12), c = 15.8735(15) Å, α = 76.295(5), β = 81.439(6), γ = 85.511(6)°, V = 1203.07(18) Å3, and Z = 2. The structure (R1 = 0.0377 for 1935 I = 2σI) contains [(UO2)2(SO4)8]4– uranyl-sulfate clusters that are linked into a band by bridging LiO4 tetrahedra. The bands are linked through peripheral SO4 tetrahedra forming a thick heteropolyhedral layer. Channels within the layers contain a K site, while an additional K site, six Na sites, and an SO3OH group occupy the space between the heteropolyhedral layers.

Geochemistry & Geophysics↗
DOE OSTI2021

Oxygen isotope ratios in zircon and garnet: A record of assimilation and fractional crystallization in the Dinkey Dome peraluminous granite, Sierra Nevada, California

Abstract The 119 Ma Dinkey Dome pluton in the central Sierra Nevada Batholith is a peraluminous granite and contains magmatic garnet and zircon that are complexly zoned with respect to oxygen isotope ratios. Intracrystalline SIMS analysis tests the relative importance of magmatic differentiation processes vs. partial melting of metasedimentary rocks. Whereas δ18O values of bulk zircon concentrates are uniform across the entire pluton (7.7‰ VSMOW), zircon crystals are zoned in δ18O by up to 1.8‰, and when compared to late garnet, show evidence of changing magma chemistry during multiple interactions of the magma with wall rock during crustal transit. The evolution from an early high-δ18O magma [δ18O(WR) = 9.8‰] toward lower values is shown by high-δ18O zircon cores (7.8‰) and lower δ 18O rims (6.8‰). Garnets from the northwest side of the pluton show a final increase in δ18O with rims reaching 8.1‰. In situ REE measurements show zircon is magmatic and grew before garnets. Additionally, δ18O in garnets from the western side of the pluton are consistently higher (avg = 7.3‰) relative to the west (avg = 5.9‰). These δ18O variations in zircon and garnet record different stages of assimilation and fractional crystallization whereby an initially high-δ18O magma partially melted low-δ18O wallrock and was subsequently contaminated near the current level of emplacement by higher δ18O melts. Collectively, the comparison of δ18O zoning in garnet and zircon shows how a peraluminous pluton can be constructed from multiple batches of variably contaminated melts, especially in early stages of arc magmatism where magmas encounter significant heterogeneity of wall-rock assemblages. Collectively, peraluminous magmas in the Sierran arc are limited to small <100 km2 plutons that are intimately associated with metasedimentary wall rocks and often surrounded by later and larger metaluminous tonalite and granodiorite plutons. The general associations suggest that early-stage arc magmas sample crustal heterogeneities in small melt batches, but that with progressive invigoration of the arc, such compositions are more effectively blended with mantle melts in source regions. Thus, peraluminous magmas provide important details of the nascent Sierran arc and pre-batholithic crustal structure.

Geochemistry & Geophysics↗
DOE OSTI2021

Hydroxylpyromorphite, a mineral important to lead remediation: Modern description and characterization

Abstract Hydroxylpyromorphite, Pb5(PO4)3(OH), has been documented in the literature as a synthetic and naturally occurring phase for some time but has not previously been formally described as a mineral. It is fully described here for the first time using crystals collected underground in the Copps mine, Gogebic County, Michigan. Hydroxylpyromorphite occurs as aggregates of randomly oriented hexagonal prisms, primarily between about 20–35 μm in length and 6–10 μm in diameter. The mineral is colorless and translucent with vitreous luster and white streak. The Mohs hardness is ~3½–4; the tenacity is brittle, the fracture is irregular, and indistinct cleavage was observed on {001}. Electron microprobe analyses provided the empirical formula Pb4.97(PO4)3(OH0.69F0.33Cl0.06)Σ1.08. The calculated density using the measured composition is 7.32 g/cm3. Powder X-ray diffraction data for the type material is compared to data previously reported for hydroxylpyromorphite from the talc mine at Rabenwald, Austria, and from Whytes Cleuch, Wanlockhead, Scotland. Hydroxylpyromorphite is hexagonal, P63/m, at 100 K, a = 9.7872(14), c = 7.3070(10) Å, V = 606.16(19) Å3, and Z = 2. The structure [R1 = 0.0181 for 494 F>4σ(F) reflections] reveals that hydroxylpyromorphite adopts a column anion arrangement distinct from other members of the apatite supergroup due to the presence of fluorine and steric constraints imposed by stereoactive lone-pair electrons of Pb2+ cations. The F– anion sites are displaced slightly from hydroxyl oxygen anions, which allows for stronger hydrogen-bonding interactions that may in turn stabilize the observed column-anion arrangement and overall structure. Our modern characterization of hydroxylpyromorphite provides deeper understanding to a mineral useful for remediation of lead-contaminated water.

Geochemistry & Geophysics↗
DOE OSTI2021

Sound velocities of iron-nickel (Fe90Ni10) alloy up to 8 GPa and 773 K: The effect of nickel on the elastic properties of bcc-iron at high P-T

Sound velocities of iron and iron-based alloys at high pressure and high temperature are crucial for understanding the composition and structure of Earth’s and other telluric planetary cores. In this study, we performed ultrasonic interferometric measurements of both compressional (ν P ) and shear (ν S ) velocities on a polycrystalline body-centered-cubic (bcc)-Fe 90 Ni 10 up to 8 GPa and 773 K. The elastic moduli and their pressure and temperature derivatives are derived from least-square fits to third-order finite strain equations, yielding K S0 = 154.2(8) GPa, G 0 = 73.2(2) GPa, K' S0 = 4.6(2), G' 0 = 1.5(1), ∂K S /∂T = –0.028(1) GPa/K, and ∂G/∂T = –0.023(1) GPa/K. Furthermore, a comparison with literature data on bcc-Fe suggests that nickel not only decreases both P and S wave velocities but also weakens the temperature effects on the elastic moduli of Fe-Ni alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Heat capacity and thermodynamic functions of partially dehydrated sodium and zinc zeolite A (LTA)

Zeolite A (LTA) is an industrially important zeolite that exhibits sorption-induced framework flexibility, the thermodynamics of which are poorly understood. In this work, we report heat capacity measurements on zinc and sodium zeolite A from 1.8 to 300 K and compare the heat capacity of water in sodium zeolite A with that of water in other zeolites. The heat capacity of zeolitic water varies significantly depending on the hydration level and identity of the host zeolite, and more tightly bound water exhibits strong inflections in its heat capacity curve. This suggests a combination of effects, including differences in water-framework binding strength and hydration-dependent flexibility transitions. We also report fits of the heat capacity data using theoretical functions, and we report values for $C_{P,m}^°$,$Δ_0^TS_m^°$,$Δ_0^TH_m^°$, and $Φ_m^°$ from 0 to 300 K. These results contribute to a systematic thermodynamic understanding of the effects of cation exchange, guest molecule confinement, and sorbate-dependent flexibility transitions in zeolites.

Geochemistry & Geophysics↗
DOE OSTI2022

Thermoelastic properties of zircon: Implications for geothermobarometry

Abstract A thermal-pressure equation of state has been determined for zircon (ZrSiO4) that characterizes its thermoelastic behavior at metamorphic conditions. New pressure-volume (P-V) data from a “Mud Tank” zircon have been collected from 1 bar to 8.47(1) GPa using X-ray diffraction, and elastic moduli were measured from room temperature up to 1172 K by resonance ultrasound spectroscopy. These data were fitted simultaneously with temperature-volume (T-V) data from the literature in EosFit7c using a new scaling technique. The parameters of a third-order Birch-Murnaghan EoS with a Mie-Grüneisen-Debye model for thermal pressure have compressional EoS parameters K0T = 224.5(1.2) GPa, K0T' = 4.90(31) with a fixed initial molar volume V0 = 39.26 cm3/mol and thermal parameters γ0 = 0.868(15), q = 2.37(80), and ΘD = 848(38) K. EoS parameters that describe the variation of unit-cell parameters with pressure and temperature were determined using an isothermal-type EoS. This new EoS confirms that zircons are stiffer than garnets and exhibit a much lower thermal expansion. This results in steep isomekes between zircon and garnets, which makes zircon trapped as inclusions in garnets at metamorphic conditions a good piezothermometer.

Geochemistry & Geophysics↗
DOE OSTI2021

Oxygen isotope evidence for input of magmatic fluids and precipitation of Au-Ag-tellurides in an otherwise ordinary adularia-sericite epithermal system in NE China

Abstract Tellurium-rich (Te) adularia-sericite epithermal Au-Ag deposits are an important current and future source of precious and critical metals. However, the source and evolution of ore-forming fluids in these deposits are masked by traditional bulk analysis of quartz oxygen isotope ratios that homogenize fine-scale textures and growth zones. To advance understanding of the source of Te and precious metals, herein, we use petrographic and cathodoluminescence (CL) images of such textures and growth zones to guide high spatial resolution secondary ion mass spectroscopy (SIMS) oxygen isotope analyses (10 µm spot) and spatially correlated fluid inclusion microthermometric measurements on successive quartz bands in contemporary Te-rich and Te-poor adularia-sericite (-quartz) epithermal Au-Ag vein deposits in northeastern China. The results show that large positive oxygen isotope shifts from –7.1 to +7.7‰ in quartz rims are followed by precipitation of Au-Ag telluride minerals in the Te-rich deposit, whereas small oxygen isotope shifts of only 4‰ (–2.2 to +1.6‰) were detected in quartz associated with Au-Ag minerals in the Te-poor deposits. Moreover, fluid-inclusion homogenization temperatures are higher in comb quartz rims (avg. 266.4 to 277.5 °C) followed by Au-Ag telluride minerals than in previous stages (~250 °C) in the Te-rich deposit. The Te-poor deposit has a consistent temperature (~245 °C) in quartz that pre- and postdates Au-Ag minerals. Together, the coupled increase in oxygen isotope ratios and homogenization temperatures followed by precipitation of Au-Ag tellurides strongly supports that inputs of magmatic fluid containing Au, Ag, and Te into barren meteoric water-dominated flow systems are critical to the formation of Te-rich adularia-sericite epithermal Au-Ag deposits. In contrast, Te-poor adularia-sericite epithermal Au-Ag deposits show little or no oxygen isotope or fluid-inclusion evidence for inputs of magmatic fluid.

Geochemistry & Geophysics↗
DOE OSTI2022

Pauling’s rules for oxide-based minerals: A re-examination based on quantum mechanical constraints and modern applications of bond-valence theory to Earth materials

Since their introduction in 1929, Pauling’s five rules have been used by scientists from many disciplines to rationalize and predict stable arrangements of atoms and coordination polyhedra in crystalline solids; amorphous materials such as silicate glasses and melts; nanomaterials, poorly crystalline solids; aqueous cation and anion complexes; and sorption complexes at mineral-aqueous solution interfaces. The predictive power of these simple yet powerful rules was challenged recently by George et al. (2020), who performed a statistical analysis of the performance of Pauling’s five rules for about 5000 oxide crystal structures. They concluded that only 13% of the oxides satisfy the last four rules simultaneously and that the second rule has the most exceptions. They also found that Pauling’s first rule is satisfied for only 66% of the coordination environments tested and concluded that no simple rule linking ionic radius to coordination environment will be predictive due to the variable quality of univalent radii. We address these concerns and discuss quantum mechanical calculations that complement Pauling’s rules, particularly his first (radius sum and radius ratio rule) and second (electrostatic valence rule) rules. We also present a more realistic view of the bonded radii of atoms, derived by determining the local minimum in the electron density distribution measured along trajectories between bonded atoms known as bond paths, i.e., the bond critical point (rc). Electron density at the bond critical point is a quantum mechanical observable that correlates well with Pauling bond strength. Moreover, a metal atom in a polyhedron has as many bonded radii as it has bonded interactions, resulting in metal and O atoms that may not be spherical. O atoms, for example, are not spherical in many oxide-based crystal structures. Instead, the electron density of a bonded oxygen is often highly distorted or polarized, with its bonded radius decreasing systematically from ~1.38 Å when bonded to highly electropositive atoms like sodium to 0.64 Å when bonded to highly electronegative atoms like nitrogen. Bonded radii determined for metal atoms match the Shannon (1976) radii for more electropositive atoms, but the match decreases systematically as the electronegativities of the M atoms increase. As a result, significant departures from the radius ratio rule in the analysis by George et al. (2020) is not surprising. We offer a modified, more fundamental version of Pauling’s first rule and demonstrate that the second rule has a one-to-one connection between the electron density accumulated between the bonded atoms at the bond critical point and the Pauling bond strength of the bonded interaction. Pauling’s second rule implicitly assumes that bond strength is invariant with bond length for a given pair of bonded atoms. Many studies have since shown that this is not the case, and Brown and Shannon (1973) developed an equation and a set of parameters to describe the relation between bond length and bond strength, now redefined as bond valence to avoid confusion with Pauling bond-strength. Brown (1980) used the valence-sum rule, together with the path rule and the valence-matching principle, as the three axioms of bond-valence theory (BVT), a powerful method for understanding many otherwise elusive aspects of crystals and also their participation in dynamic processes. We show how a priori bond-valence calculations can predict unstrained bond-lengths and how bond-valence mapping can locate low-Z atoms in a crystal structure (e.g., Li) or examine possible diffusion pathways for atoms through crystal structures. In addition, we briefly discuss Pauling’s third, fourth, and fifth rules, the first two of which concern the sharing of polyhedron elements (edges and faces) and the common instability associated with structures in which a polyhedron shares an edge or face with another polyhedron and contains high-valence cations. The olivine [α-(MgxFe1–x)2SiO4] crystal structure is used to illustrate the distortions from hexagonal close-packing of O atoms caused by metal-metal repulsion across shared polyhedron edges. We conclude by discussing several applications of BVT to Earth materials, including the use of BVT to: (1) locate H+ ions in crystal structures, including the location of protons in the crystal structures of nominally anhydrous minerals in Earth’s mantle; (2) determine how strongly bonded (usually anionic) structural units interact with weakly bonded (usually cationic) interstitial complexes in complex uranyl-oxide and uranyl-oxysalt minerals using the valence-matching principle; (3) calculate Lewis acid strengths of cations and Lewis base strengths of anions; (4) determine how (H2O) groups can function as bond-valence transformers by dividing one bond into two bonds of half the bond valence; (5) help characterize products of sorption reactions of aqueous cations (e.g., Co2+ and Pb2+) and oxyanions [e.g., selenate (Se6+O4)2- and selenite (Se4+O3)2-] at mineral-aqueous solution interfaces and the important role of protons in these reactions; and (6) help characterize the local coordination environments of highly charged cations (e.g., Zr4+, Ti4+, U4+, U5+, and U6+) in silicate glasses and melts.

Geochemistry & Geophysics↗
DOE OSTI2022

Oxygen diffusion in garnet: Experimental calibration and implications for timescales of metamorphic processes and retention of primary O isotopic signatures

Knowledge of oxygen diffusion in garnet is crucial for a correct interpretation of oxygen isotope signatures in natural samples. A series of experiments was undertaken to determine the diffusivity of oxygen in garnet, which remains poorly constrained. The first suite included high-pressure (HP), nominally dry experiments performed in piston-cylinder apparatus at: (1) T = 1050–1600 °C and P = 1.5 GPa and (2) T = 1500 °C and P = 2.5 GPa using yttrium aluminum garnet (YAG; Y3Al5O12) cubes. Second, HP H2O-saturated experiments were conducted at T = 900 °C and P = 1.0–1.5 GPa, wherein YAG crystals were packed into a YAG + Corundum powder, along with 18O-enriched H2O. Third, 1 atm experiments with YAG cubes were performed in a gas-mixing furnace at T = 1500–1600 °C under Ar flux. Finally, an experiment at T = 900 °C and P = 1.0 GPa was done using a pyrope cube embedded into pyrope powder and 18O-enriched H2O. Experiments using grossular were not successful. Profiles of 18O/(18O+16O) in the experimental charges were analyzed with three different secondary ion mass spectrometers (SIMS): sensitive high-resolution ion microprobe (SHRIMP II and SI), CAMECA IMS-1280, and NanoSIMS. Considering only the measured length of 18O diffusion profiles, similar results were obtained for YAG and pyrope annealed at 900 °C, suggesting limited effects of chemical composition on oxygen diffusivity. However, in both garnet types, several profiles deviate from the error function geometry, suggesting that the behavior of O in garnet cannot be fully described as simple concentration-independent diffusion, certainly in YAG and likely in natural pyrope as well. The experimental results are better described by invoking O diffusion via two distinct pathways with an inter-site reaction allowing O to move between these pathways. Modeling this process yields two diffusion coefficients (D values) for O, one of which is approximately two orders of magnitude higher than the other. Taken together, Arrhenius relationships are:logDm2s-1=-7.2(±1.3)+(-321(±32)kJmol-12.303RT) for the slow pathway, andlogDm2s-1=-5.4(±0.7)+(-321(±20)kJmol-12.303RT) for the fast pathway. We interpret the two pathways as representing diffusion following vacancy and inter-stitial mechanisms, respectively. Regardless, our new data suggest that the slow mechanism is prevalent in garnet with natural compositions, and thus is likely to control the retentivity of oxygen isotopic signatures in natural samples. The diffusivity of oxygen is similar to Fe-Mn diffusivity in garnet at 1000–1100 °C and Ca diffusivity at 850 °C. However, the activation energy for O diffusion is larger, leading to lower diffusivities at P-T conditions characterizing crustal metamorphism. Therefore, original O isotopic signatures can be retained in garnets showing major element zoning partially re-equilibrated by diffusion, with the uncertainty caveat of extrapolating the experimental data to lower temperature conditions.

Geochemistry & Geophysics↗
DOE OSTI2022

Experimental calibration of an Fe3+/Fe2+-in-amphibole oxybarometer and its application to shallow magmatic processes at Shiveluch Volcano, Kamchatka

Abstract Oxygen fugacity is an important but difficult parameter to constrain for primitive arc magmas. In this study, the partitioning behavior of Fe3+/Fe2+ between amphibole and glass synthesized in piston-cylinder and cold-seal apparatus experiments is developed as an oxybarometer, applicable to magmas ranging from basaltic to dacitic composition. The partitioning of Fe2+ is strongly dependent on melt polymerization; the relative compatibility of Fe2+ in amphibole decreases with increasing polymerization. The Fe2+/Mg distribution coefficient between amphibole and melt is a relatively constant value across all compositions and is, on average, 0.27. The amphibole oxybarometer is applied to amphibole in mafic enclaves, cumulates, and basaltic tephra erupted from Shiveluch volcano in Kamchatka with measured Fe3+/FeTotal. An average Fe3+/Fe2+ amphibole-glass distribution coefficient for basalt is used to convert the Fe3+/FeTotal of amphibole in samples from Shiveluch to magmatic oxygen fugacity relative to NNO. The fO2 of primitive melts at the volcano is approximately NNO+2 and is faithfully recorded in amphibole from an amphibole-rich cumulate and the basaltic tephra. Apparently, higher fO2 recorded by amphibole in mafic enclaves likely results from partial dehydrogenation of amphibole during residence in a shallow andesite storage region. We identify three pulses of mafic magma recharge within two weeks of, a month before, and two to three months before the eruption and find that, at each of these times, the host andesite was recharged by at least two magmas at varying stages of differentiation. Application of the amphibole oxybarometer not only gives insight into magmatic fO2 but also potentially details of shallow magmatic processes.

Geochemistry & Geophysics↗
DOE OSTI2022

Sound speed and refractive index of amorphous CaSiO3 upon pressure cycling to 40 GPa

Abstract Brillouin spectroscopy at room temperature and pressures up to 40 GPa documents nearly identical elasticity and refractive index of amorphous CaSiO3 created by two different methods: temperature-quenching the melt at ambient pressure and pressure-amorphizing crystalline wollastonite at room temperature. We find reproducible hysteresis of 0 to 8% on pressure cycling that is small relative to the 30 to 60% changes in shear and longitudinal wave velocities over this pressure range. Together with observed changes in refractive index and previous results from Raman spectroscopy, these measurements reveal a continuous and reversible change in atomic packing induced by pressure. Unlike many other silicate glasses, amorphous CaSiO3 exhibits highly reproducible properties, behaving smoothly and reversibly under pressure cycling and possessing similar structure and elasticity regardless of synthesis paths for the starting material, which suggests that the amorphous solid may mimic the liquid over the pressure range investigated.

Geochemistry & Geophysics↗
DOE OSTI2023

X-ray absorption spectroscopy study of Mn reference compounds for Mn speciation in terrestrial surface environments

Abstract X-ray absorption spectroscopy (XAS) offers great potential to identify and quantify Mn species in surface environments by means of linear combination fit (LCF), fingerprint, and shell-fit analyses of bulk Mn XAS spectra. However, these approaches are complicated by the lack of a comprehensive and accessible spectrum library. Additionally, molecular-level information on Mn coordination in some potentially important Mn species occurring in soils and sediments is missing. Therefore, we investigated a suite of 32 natural and synthetic Mn reference compounds, including Mn oxide, oxyhydroxide, carbonate, phosphate, and silicate minerals, as well as organic and adsorbed Mn species, by Mn K-edge X-ray absorption near edge structure (XANES) and extended X-ray absorption fine structure (EXAFS) spectroscopy. The ability of XAS to infer the average oxidation state (AOS) of Mn was assessed by comparing XANES-derived AOS with the AOS obtained from redox titrations. All reference compounds were studied for their local (<5 Å) Mn coordination environment using EXAFS shell-fit analysis. Statistical analyses were employed to clarify how well and to what extent individual Mn species (groups) can be distinguished by XAS based on spectral uniqueness. Our results show that LCF analysis of normalized XANES spectra can reliably quantify the Mn AOS within ~0.1 v.u. in the range +2 to +4. These spectra are diagnostic for most Mn species investigated, but unsuitable to identify and quantify members of the manganate and Mn(III)-oxyhydroxide groups. First-derivative XANES fingerprinting allows the unique identification of pyrolusite, ramsdellite, and potentially lithiophorite within the manganate group. However, XANES spectra of individual Mn compounds can vary significantly depending on chemical composition and/or crystallinity, which limits the accuracy of XANES-based speciation analyses. In contrast, EXAFS spectra provide a much better discriminatory power to identify and quantify Mn species. Principal component and cluster analyses of k2-weighted EXAFS spectra of Mn reference compounds implied that EXAFS LCF analysis of environmental samples can identify and quantify at least the following primary Mn species groups: (1) Phyllo- and tectomanganates with large tunnel sizes (2 × 2 and larger; hollandite sensu stricto, romanèchite, todorokite); (2) tectomanganates with small tunnel sizes (2 × 2 and smaller; cryptomelane, pyrolusite, ramsdellite); (3) Mn(III)-dominated species (nesosilicates, oxyhydroxides, organic compounds, spinels); (4) Mn(II) species (carbonate, phosphate, and phyllosilicate minerals, adsorbed and organic species); and (5) manganosite. All Mn compounds, except for members of the manganate group (excluding pyrolusite) and adsorbed Mn(II) species, exhibit unique EXAFS spectra that would allow their identification and quantification in mixtures. Therefore, our results highlight the potential of Mn K-edge EXAFS spectroscopy to assess bulk Mn speciation in soils and sediments. A complete XAS-based speciation analysis of bulk Mn in environmental samples should preferably include the determination of Mn valences following the “Combo” method of Manceau et al. (2012), EXAFS LCF analyses based on principal component and target transformation results, as well as EXAFS shell-fit analyses for the validation of LCF results. For this purpose, all 32 XAS reference spectra are provided in the Online Materials1 for further use by the scientific community.

Geochemistry & Geophysics↗
DOE OSTI2023

The influence of Al2O3 on the structural properties of MgSiO3 akimotoite

Abstract Akimotoite, a MgSiO3 polymorph present in the lower transition zone within ultramafic portions of subducting slabs and potentially also in the ambient mantle, will partition some amount of Al, raising the question of how this will affect its crystal structure and properties. In this study, a series of samples along the MgSiO3-Al2O3 (akimotoite-corundum) solid solution have been investigated by means of single-crystal X-ray diffraction to examine their crystal chemistry. Results show a strong nonlinear behavior of the a- and c-axes as a function of Al content, which arises from fundamentally different accommodation mechanisms in the akimotoite and corundum structures. Furthermore, two Al2O3-bearing akimotoite samples were investigated at high pressure to determine the different compression mechanisms associated with Al substitution. Al2O3-bearing akimotoite becomes more compressible at least up to 20 mol% Al2O3, due likely to an increase in compressibility as the Al cation is incorporated into the SiO6 octahedron. This observation is in strong contrast to the stiffer corundum end-member having a KT = 250 GPa, which is larger than that of the akimotoite end-member [KT = 205(1) GPa]. These findings have implications for mineral physics models of elastic properties, which have in the past assumed linear mixing behavior between the MgSiO3 akimotoite and Al2O3 corundum end-members to calculate sound wave velocities for Al-bearing akimotoite at high pressure and temperature.

Geochemistry & Geophysics↗
DOE OSTI2022

Hydrogen solubility in FeSi alloy phases at high pressures and temperatures

Abstract Light elements alloying with metallic Fe can change the properties and therefore play a key role in the structure and dynamics of planetary cores. Hydrogen and silicon are possible light elements in the rocky planets’ cores. However, hydrogen storage in Fe-Si alloy systems remains unclear at high pressures and high temperatures because of experimental difficulties. Taking advantage of pulsed laser heating combined with high-energy synchrotron X-ray diffraction, we studied reactions between FeSi and H in laser-heated diamond-anvil cells (LHDACs) up to 61.9 GPa and 3500 K. We found that under H-saturated conditions the amount of H alloying with FeSi (0.3 and <0.1 wt% for the B20 and B2 structures, respectively) is much smaller than that in pure Fe metal (>1.8 wt%). Our experiments also suggest that H remains in the crystal structure of FeSi alloy when recovered to 1 bar. Further density functional theory (DFT) calculations indicate that the low-H solubility likely results from the highly distorted interstitial sites in the B20 and B2 structures, which are not favorable for H incorporation. The recovery of H in the B20 FeSi crystal structure at ambient conditions could open up possibilities to understand geochemical behaviors of H during core formation in future experiments. The low-H content in FeSi alloys suggests that if a planetary core is Si-rich, Si can limit the ingassing of H into the Fe-rich core.

Geochemistry & Geophysics↗
DOE OSTI2022

Calorimetric study of skutterudite (CoAs2.92) and heazlewoodite (Ni3S2)

Abstract Nickel and cobalt arsenides, sulfarsenides, and sulfides occur in many hydrothermal ore deposits, but their thermodynamic properties are not well known, in some cases not known at all. In this work, we determined a full set of thermodynamic properties for heazlewoodite and skutterudite. Both phases were synthesized in evacuated silica tubes at elevated temperatures, and electron microprobe analyses gave their compositions as Ni3S2 and CoAs2.92, respectively. Enthalpies of formation were measured by high-temperature oxide-melt solution calorimetry. The reference phases were pure elements, thus eliminating any systematic errors related to such phases. The enthalpies of formation at T = 298.15 K and P = 105 Pa are –216.0 ± 8.4(2σ) and –88.2 ± 6.1 kJ·mol−1 for Ni3S2 and CoAs2.92, respectively. Entropies were calculated from low-temperature heat capacity (CP) data from relaxation (PPMS) calorimetry and are 133.8 ± 1.6 and 106.4 ± 1.3 J·mol–1·K–1, respectively. The calculated Gibbs free energies of formation are –210.0 ± 8.4 and –79.9 ± 6.2 kJ·mol−1 for Ni3S2 and CoAs2.92, respectively. The PPMS CP data, together with a set of differential scanning calorimetry measurements, were used to derive CP polynomials up to 700 K with the Kieffer model based on previously published frequencies of acoustic and optic modes. Equilibrium constants for selected reactions with an aqueous phase were calculated up to 700 K. Geochemical modeling in these systems, however, should await until more reliable data for other phases from the system Co-Ni-As-S are available.

Geochemistry & Geophysics↗
DOE OSTI2023

Eu speciation in apatite at 1 bar: An experimental study of valence-state partitioning by XANES, lattice strain, and Eu/Eu* in basaltic systems

Abstract Partition coefficients for rare earth elements (REEs) between apatite and basaltic melt were determined as a function of oxygen fugacity (fO2; iron-wüstite to hematite-magnetite buffers) at 1 bar and between 1110 and 1175 °C. Apatite-melt partitioning data for REE3+ (La, Sm, Gd, Lu) show near constant values at all experimental conditions, while bulk Eu becomes more incompatible (with an increasing negative anomaly) with decreasing fO2. Experiments define three apatite calibrations that can theoretically be used as redox sensors. The first, a XANES calibration that directly measures Eu valence in apatite, requires saturation at similar temperature-composition conditions to experiments and is defined by: ( E u 3 + ∑ E u ) Apatite = 1 1 + 10 - 0.10 ± 0.01 × l o g ⁡ ( f o 2 ) - 1.63 ± 0.16 . The second technique involves analysis of Sm, Eu, and Gd in both apatite and coexisting basaltic melt (glass), and is defined by: ( Eu E u * ) D Sm × Gd = 1 1 + 10 - 0.15 ± 0.03 × l o g ⁡ ( f o 2 ) - 2.46 ± 0.41 . The third technique is based on the lattice strain model and also requires analysis of REE in both apatite and basalt. This calibration is defined by ( Eu E u * ) D lattice strain = 1 1 + 10 - 0.20 ± 0.03 × l o g ⁡ ( f o 2 ) - 3.03 ± 0.42 . The Eu valence-state partitioning techniques based on (Sm×Gd) and lattice strain are virtually indistinguishable, such that either methodology is valid. Application of any of these calibrations is best carried out in systems where both apatite and coexisting glass are present and in direct contact with one another. In holocrystalline rocks, whole rock analyses can be used as a guide to melt composition, but considerations and corrections must be made to either the lattice strain or Sm×Gd techniques to ensure that the effect of plagioclase crystallization either prior to or during apatite growth can be removed. Similarly, if the melt source has an inherited either a positive or negative Eu anomaly, appropriate corrections must also be made to lattice strain or Sm×Gd techniques that are based on whole rock analyses. This being the case, if apatite is primary and saturates from the parent melt early during the crystallization sequence, these corrections may be minimal. The partition coefficients for the REE between apatite and melt range from a maximum DEu3+ = 1.67 ± 0.25 (as determined by lattice strain) to DLu3+ = 0.69 ± 0.10. The REE partition coefficient pattern, as observed in the Onuma diagram, is in a fortuitous situation where the most compatible REE (Eu3+) is also the polyvalent element used to monitor fO2. These experiments provide a quantitative means of assessing Eu anomalies in apatite and how they be used to constrain the oxygen fugacity of silicate melts.

Geochemistry & Geophysics↗