Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “geoscience”

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,765 records · Page 98
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

The Stability of Nd Hydroxyl Complexes at Near Neutral to Alkaline pH from 25 to 75 °C: Implications for Monazite Solubility in Hydrothermal Aqueous Fluids

Speciation of rare earth elements (REE) in aqueous fluids controls their mobilization during fluid-rock interaction. Thermodynamic modeling provides important insights into the factors controlling REE mobility and fractionation in aqueous fluids, in which aqueous REE hydroxyl complexes are particularly important across a broad pH range. However, the stability of these REE hydroxyl complexes remains poorly determined due to limited experimental data above ambient temperature. In this work, UV–vis spectrophotometry experiments were conducted to determine the cumulative hydrolysis constants for neodymium (Nd) from 25 to 75 °C in near-neutral to alkaline solutions. The color indicator m -cresol purple was used to determine in situ pH. Alkaline NaOH-bearing aqueous solutions were doped with varying initial Nd concentrations (0 to ∼0.155 mmol/kg), which resulted in the release of protons (H + ) and pH decrease during hydrolysis of Nd 3+ according to Nd 3+ + n H 2 O = Nd(OH) n 3– n + H + . The average OH – ligand number coordinated to Nd 3+ was found to increase from 1.0 to 1.6 at 25 °C and from 1.2 to 2.8 at 75 °C over a pH range from 6.3 to 9.0. The measured speciation shows an increased predominance of Nd(OH) 3 0 over the Nd(OH) 2 + and Nd(OH) 2+ species with increased temperature and pH. The increasing cumulative formation constants (log*β° n , n = 1 to 3) retrieved from 25 to 75 °C differ by 0.1 to 1.9 logarithmic units in comparison to existing thermodynamic databases. These updated thermodynamic data have important implications for geochemical modeling of the speciation of Nd hydroxyl complexes and the solubility of monazite as a function of pH.

58 GEOSCIENCES↗
DOE OSTI2026

Spectrophotometric study of erbium (III) speciation in chloride solutions from 25 to 75 °C

In this study, UV-Vis spectrophotometric experiments were conducted from 25 to 75 °C using Er-HClO4-bearing solutions with varying NaCl concentrations. The molar absorbance coefficient of Er3+ decreases with temperature, and the absorbance intensity gradually decreases with increasing mCl/mEr ratios. The number of absorbing species include Er3+ and ErCl2+ even at the highest mCl/mEr ratio. The formation constant (β 1 ) for the ErCl 2+ species was derived and fitted between 50 and 250 °C with existing literature values, yielding the following equation: logβ 1 0 = −22.4315 + 0.0345T + 3.9379·10 3 /T, where T is temperature in Kevin.

58 GEOSCIENCES↗
DOE OSTI2019

Relationship between Atmospheric Aerosol Mineral Surface Area and Iron Solubility

Size-fractionated dust aerosols (>7.2, 7.2-3, 3-1.5, 1.5-0.95, 0.95-0.49, and <0.49 mu m) originating from source regions, including the Sahara, North America, and the local marine environment, were collected shipboard in the Atlantic Ocean. Iron mineralogy, oxidation state, and solubility were examined using wet chemistry and synchrotron-based iron K-edge X-ray absorption near edge structure (XANES) spectroscopy. The soluble iron concentration showed a strong positive correlation with calculated aerosol an surface area (R 2 = 0.81). The iron XANES spectra, for all size ranges, were nearly identical for Saharan samples. Spectra from North American and local atmosphere (marine) provenances exhibited more variation between size fractions, with no obvious trends with the particle size. Saharan spectra were distinct from North American spectra, which confirms the different compositions of the two sources. Marine spectra displayed a range of both North American and Saharan spectral characteristics, likely reflecting that the marine samples are ultimately derived from these and potentially other sources. However, despite the notable variation in compositions, the surface area is the predominant factor affecting total iron solubility. Aerosol iron solubilization through acidic surface reactions would be most effective in samples with the highest surface area, which likely results in the strong positive correlation with soluble iron concentration in distilled water.

58 GEOSCIENCES↗
DOE OSTI2026

Mitigating Hydrogen-Induced Degradation of Iridium Anodes in Proton Exchange Membrane Water Electrolyzers

One promising method for reducing precious metal usage in proton exchange membrane water electrolysis is lowering iridium (Ir) loading at the anode. However, low-loading Ir catalysts often suffer from poor stability under high current densities. In this study, hydrogen (H2) crossover from the cathode to the anode is identified as a key degradation pathway. Temperature-programmed reduction confirms the reduction of IrO2 at 80 °C in a hydrogen environment, highlighting the vulnerability of IrO2-based catalysts to H2 exposure. To mitigate this effect, palladium (Pd) is introduced as an anode additive, acting as an H2 oxidation catalyst and mitigating IrO2 reduction. This protective role is verified by inductively coupled plasma optical emission spectroscopy and in situ X-ray absorption spectroscopy, showing significantly suppressed Ir dissolution at 80 °C under H2 flow when an O-covered Pd surface is present at oxygen evolution reaction potentials. Results from the current study identify a new strategy in improving activity and durability of catalysts in electrolyzers.

58 GEOSCIENCES↗
DOE OSTI2023

Functional Predictor Variables for the Leaching Potential of Arsenic and Selenium from Coal Fly Ash

The release of leachates from intact coal ash impoundments is a concern due to the enrichment and mobilization of toxic elements such as arsenic (As) and selenium (Se). This study aims to explore the intrinsic properties of coal fly ash that correlate with the relative leachability of As and Se. We performed leaching experiments with 52 fly ash samples collected from 15 different U.S. power plants and representing coal feedstocks from the three major domestic regions. We assessed the mobilization potential of As and Se in fly ash based on standardized leaching protocols and performed multivariate and lasso regression analyses to explore correlations of leachable As and Se contents with characteristics such as major element contents, loss on ignition, and pH. The results of regression models indicated that major elements (Fe, Ca, and Al) for a wide range of fly ashes can serve as predictor variables for the leaching potential of As but not for Se. LOI and pH were not important predictive variables in the models. Both regression approaches resulted in relatively strong fits for leachable As (correlation coefficient R 2 = 0.78 for both models) compared to models for leachable Se (R 2 = 0.49). Overall, these results suggest that correlation models combined with on-site elemental analysis with portable analyzers may enable a screening method for leachable As in coal ash.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2025

Molecular Insights into Geochemical Reactions of Iron-Bearing Minerals: Implications for Hydrogen Geo-Storage

Here, this study investigates the reaction of hydrogen (H 2 ) with pyrite (FeS 2 ), focusing on how temperature and the presence of water influence the reaction pathways and kinetics. Utilizing computational molecular simulations and kinetic analyses, we explore the impact of these factors on the formation of hydrogen sulfide (H 2 S) and related species. First, grand canonical Monte Carlo/molecular dynamics (GCMC/MD) simulations reveal that physical H 2 adsorption occurs in distinct layers on the pyrite surface. In addition, increased temperatures reduce the absolute adsorption capacities. Reactive MD simulations demonstrate that H 2 interacts differently with pyrite under varying conditions. At 298 K, H 2 reacts with pyrite to form HS – , leading to the formation of HS – through covalent bonding with sulfur of pyrite. However, no H 2 S is produced at this temperature, suggesting that a kinetic barrier (i.e., activation energy) may prevent this reaction. At higher temperatures, H 2 S production significantly increases. The presence of water introduces additional complexity to the reaction mechanism. Unlike dry conditions, water enhances H 2 S generation, even at low temperatures. Water also facilitates the formation of additional products, such as SOH, indicating a more intricate chemical environment on the pyrite surface. Our findings identify the association of HS – ions to form H 2 S as the rate-limiting step, with temperature influencing this process. This finding suggests that while the presence of water can create a more dynamic reaction environment, the overall mechanisms leading to H 2 S formation remain consistent. These outcomes suggest the need for developing targeted strategies to manage and control H 2 S emissions within the context of underground hydrogen storage.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Mine Tailings Valorization by Electrochemically Stimulated Mineralization from Mildly Acidic Conditions

This study presents a high-efficiency electrochemical process for the mineralization of calcium carbonate (CaCO 3 ) from mildly acidic mine tailing supernatant water. Electrochemical pH control was used to promote carbonate speciation and precipitation from CO 2 -saturated solutions, achieving >1 mol CaCO 3 precipitated per mol e - at applied potentials between -1.4 and -1.6 V vs Ag/AgCl. Product morphology and polymorph selectivity were tunable via applied potential, yielding calcite and vaterite phases. Experiments using real mine tailings water confirmed selective calcite precipitation, despite the presence of sulfate and other trace elements. These results highlight a viable route for coupling CO 2 utilization with mine tailings valorization.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Achieving Higher Activity of Acidic Oxygen Evolution Reaction Using an Atomically Thin Layer of IrO x over Co 3 O 4

The development of electrocatalysts with reduced iridium (Ir) loading for the oxygen evolution reaction (OER) is essential to produce low-cost green hydrogen from water electrolysis under acidic conditions. Herein, an atomically thin layer of iridium oxide (IrO x ) has been uniformly dispersed onto cobalt oxide (Co 3 O 4 ) nanocrystals to improve the efficient use of Ir for acidic OER. In situ characterization and theoretical calculations reveal that compared to the conventional IrOx cluster, the atomically thin layer of IrO x shows stronger interaction with the Co 3 O 4 and consequently higher OER activity due to the Ir-O-Co bond formation at the interface. Equally important, the facile synthetic method and the promising activity in the proton exchange membrane water electrolyzer, reaching 1 A cm -2 at 1.7 V with remarkable durability, enable potential scale-up applications. In conclusion, these findings provide a mechanistic understanding for designing active, stable and lower-cost electrocatalysts with well-defined structures for acidic OER.

58 GEOSCIENCES↗
DOE OSTI2025

Dynamic Atomistic Polar Structure Underpins Ultrahigh Linear Electro-Optic Coefficient in Transparent Ferroelectric Ceramics

Transparent ferroelectrics with high linear electro-optic (EO) coefficients are critical for advanced electro-optical devices. However, achieving optical transparency in ferroelectric ceramics remains challenging due to visible light scattering caused by defects such as domain walls, grain boundaries, and pores. Here, we report the successful fabrication of transparent ferroelectric ceramics through innovative chemical composition design and an advanced two-step sintering process in the La-doped Pb(Mg 1/3 Nb 2/3 )O 3 –PbTiO 3 system. The optical transparency, which is near the theoretical upper limit, can be attributed to the wide band gap and the minimization of light scattering of defects. By minimizing porosity and engineering grain/domain sizes to differ significantly from the wavelengths of visible light, we suppress scattering, achieving optical transparency near the theoretical upper limit. Strikingly, these ceramics exhibit an ultrahigh linear EO coefficient of ∼1417 pm/V, over 65 times greater than that of LiNbO 3 single crystals, the current industry standard. We attribute this exceptional performance to dynamic atomistic polar structures within switchable, thermally stable domains, which enhance electronic polarization sensitivity. This mechanism is corroborated by dielectric spectroscopy, high-resolution transmission electron microscopy and simulation. Our findings offer insights into the design of cost-effective transparent materials with exceptional EO properties, paving the way for next-generation electro-optical devices.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Probing Terra Incognita of Ni–P Catalysts: Operando Explorations during Hydrogen Evolution Reaction

We have developed two Ni phosphide preparation methods allowing operando XAS surface-sensitive studies of well-defined bulk systems. For Ni K-edge XAS, a Ni 2 P phase-pure powder was sintered into a high-density pellet and polished for grazing incidence XAS. Ni sites were mildly affected by the acidic electrolyte prior to the HER, while the applied cathodic potential caused the reduction of Ni surface states beyond the states of as-prepared Ni 2 P. The computed fully H-covered Ni 2 P [0001] model describes the difference in the operando Ni K-edge GIXAS spectrum well. Upon turning the applied bias off, the Ni sites became immediately oxidized, forming NiO on the surface. Thus, the active phase during the HER is covalent Ni 0 close to that in the intermetallic phosphides, and Ni 2+ oxides formed after, and not during, the HER. For P K-edge XAS, Ni foam was phosphorized to form a thin Ni 3 P layer while preserving its high surface area. Upon immersion in the acidic electrolyte, the P sites underwent removal of P 5+ phosphates and formed new P coordination, possibly due to the adsorption of protons from the electrolyte. These new P surface states were not affected by turning the cathodic current on and off as soon as the sample was immersed in the acidic electrolyte. However, the removal of the sample from the electrochemical cell and drying in air resulted in substantial depletion and oxidation of surface P. Echoing the observed Ni site chemistry during HER, XAS and XPS suggest that the in situ active P sites are different from the oxidized P states observed under ex situ conditions.

Kong, Seongyoung [Iowa State Univ., Ames, IA (Unit↗
DOE OSTI2026

Activationof Oxygen Evolution Electrocatalysis viaReduced Ruthenium–Oxygen–Ruthenium Coordination

Noble metal oxides such as RuO2 are the state-of-the-art electrocatalysts for anodic reactions in acidic electrolytes, but their scarcity and moderate activity greatly limit emerging renewable energy technologies. Here, we show that oxidized overlayers of ruthenium on earth-abundant manganese oxide (MnO2/o-RuOx) nanocrystal supports exhibit Ru chemical states associated with reduced Ru–O–Ru coordination that enable dynamic switching of hydrogen bonding, with *OH intermediates hydrogen bonding to surface O and *OOH intermediates bonding to protruding RuOx clusters. The resulting electrocatalysts exhibit an overpotential of 218.9 ± 0.3 mV at 10 mA cm–2 for the oxygen evolution reaction in acid, corresponding to a 2425% increase in Ru mass activity compared to RuO2, enabling the construction of electrolyzers that achieved 3 A cm–2 at 1.646 V, 5.54 A cm–2 at 1.8 V, and exhibited over 3000-h stability at 100 mA cm–2. These findings motivate further efforts to develop nanomaterials that harness reduced Ru–O–Ru coordination to enable emerging renewable energy technologies.

58 GEOSCIENCES↗
DOE OSTI2019

Effects of Fault Roughness on Coseismic Slip and Earthquake Locations

We report that fault zone structure is well known to exert strong controls on earthquake properties including coseismic slip distribution, rupture propagation direction, and hypocenter location. It has also been well established that the principal slip surface, which accommodates the majority of earthquake displacement, exhibits roughness at all scales following self-affine fractal distributions. Here we explore the relationship between fault roughness and specific earthquake properties including coseismic slip distribution and hypocenter location based on long-term simulations of earthquake catalogs on fractally rough faults. We begin by using the von Kármán autocorrelation function to procedurally generate single faults with different fractal roughness properties, which we place in a homogeneous elastic solid and apply pure right-lateral shear at a constant back slip rate with the earthquake simulator RSQSim. Running the simulations for 10,000 years each, we generate millions of earthquakes including thousands of events with Mw > 6.0, which rupture the surface. We show that the patterns of surface rupture in these large events follow self-affine fractal distributions with consistent fractal dimension related distinctly to the fractal dimension of the fault. In addition, the hypocenters of these large events occur in very specific predictable locations where the longest wavelength structure produces a stress asperity (i.e., restraining bend). The resulting patterns can explain many features observed on real fault systems, including clustered hypocenter locations, spatially variable coseismic slip distributions, and characteristic slip recurrence behavior. Lastly, these results demonstrate a quantitative link between a directly measurable fault property—roughness—and the properties of future earthquakes.

58 GEOSCIENCES↗
DOE OSTI2020

Timescales and Processes of Methane Hydrate Formation and Breakdown, With Application to Geologic Systems

Abstract Gas hydrate is an ice‐like form of water and low molecular weight gas stable at temperatures of roughly −10°C to 25°C and pressures of ~3 to 30 MPa in geologic systems. Natural gas hydrates sequester an estimated one sixth of Earth's methane and are found primarily in deepwater marine sediments on continental margins, but also in permafrost areas and under continental ice sheets. When gas hydrate is removed from its stability field, its breakdown has implications for the global carbon cycle, ocean chemistry, marine geohazards, and interactions between the geosphere and the ocean‐atmosphere system. Gas hydrate breakdown can also be artificially driven as a component of studies assessing the resource potential of these deposits. Furthermore, geologic processes and perturbations to the ocean‐atmosphere system (e.g., warming temperatures) can cause not only dissociation, but also more widespread dissolution of hydrate or even formation of new hydrate in reservoirs. Linkages between gas hydrate and disparate aspects of Earth's near‐surface physical, chemical, and biological systems render an assessment of the rates and processes affecting the persistence of gas hydrate an appropriate Centennial Grand Challenge. This paper reviews the thermodynamic controls on methane hydrate stability and then describes the relative importance of kinetic, mass transfer, and heat transfer processes in the formation and breakdown (dissociation and dissolution) of gas hydrate. Results from numerical modeling, laboratory, and some field studies are used to summarize the rates of hydrate formation and breakdown, followed by an extensive treatment of hydrate dynamics in marine and cryospheric gas hydrate systems.

58 GEOSCIENCES↗
DOE OSTI2019

The Effects of Earthquakes and Fluids on the Metamorphism of the Lower Continental Crust

Rock rheology and density have first-order effects on the lithosphere's response to plate tectonic forces at plate boundaries. Changes in these rock properties are controlled by metamorphic transformation processes that are critically dependent on the presence of fluids. At the onset of a continental collision, the lower crust is in most cases dry and strong. However, if exposed to internally produced or externally supplied fluids, the thickened crust will react and be converted into a mechanically weaker lithology by fluid-driven metamorphic reactions. Fluid introduction is often associated with deep crustal earthquakes. Microstructural evidence, suggest that in strong highly stressed rocks, seismic slip may be initiated by brittle deformation and that wall-rock damage caused by dynamic ruptures plays a very important role in allowing fluids to enter into contact with dry and highly reactive lower crustal rocks. The resulting metamorphism produces weaker rocks which subsequently deform by viscous creep. Volumes of weak rocks contained in a highly stressed environment of strong rocks may experience significant excursions toward higher pressure without any associated burial. Slow and highly localized creep processes in a velocity strengthening regime may produce mylonitic shear zones along faults initially characterized by earthquake-generated frictional melting and wall rock damage. However, stress pulses from earthquakes in the shallower brittle regime may kick start new episodes of seismic slip at velocity weakening conditions. These processes indicate that the evolution of the lower crust during continental collisions is controlled by the transient interplay between brittle deformation, fluid-rock interactions, and creep flow.

58 GEOSCIENCES↗
DOE OSTI2019

Three–Dimensional Acoustic Multipole Waveform Inversion at Yasur Volcano, Vanuatu

Acoustic waveform inversions can provide estimates of volume flow rate and erupted mass, enhancing the ability to estimate volcanic emissions. Previous studies have generally assumed a simple acoustic source (monopole); however, more complex and accurate source reconstructions are possible with a combination of equivalent sources (multipole). We deployed a high–density acoustic network around Yasur volcano, Vanuatu, including acoustic sensors on a tethered aerostat that was moved every ~15–60 min. Using this unique data set we invert for the acoustic multipole source mechanism using a grid search approach for 80 events to examine volume flow rates and dipole strengths. Our method utilizes finite–difference time–domain modeling to obtain the full 3–D Green's functions that account for topography. Inversion results are compared using a monopole–only, multipole (monopole and dipole), simulations that do not include topography, and those that use a subset of sensors. We find that the monopole source is a good approximation when topography is considered. However, initial compression amplitude is not fully captured by a monopole source so source directionality cannot be ruled out. The monopole solution is stable regardless of whether a monopole–only or multipole inversion is performed. Inversions for the dipole components produce estimates consistent with observed source directionality, though these inversions are somewhat unstable given station configurations of typical deployments. Furthermore, our results suggest that infrasound waveform inversion shows promise for realistic quantitative source estimates, but additional work is necessary to fully explore inversion stability, uncertainty, and robustness.

58 GEOSCIENCES↗