Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Aqueous solution”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Green Electrode Processing Enabled by Fluoro-Free Multifunctional Binders for Lithium-Ion Batteries

The eco-friendly processing of conjugated polymer binder for lithium-ion batteries demands improved polymer solubility by introducing functional moieties, while this strategy will concurrently sacrifice polymer conductivity. Employing the polyfluorene-based binder poly(2,7-9,9 (di(oxy-2,5,8-trioxadecane))fluorene) (PFO), soluble in water-ethanol mixtures, a novel approach is presented to solve this trade-off, which features integration of aqueous solution processing with subsequent controlled thermal-induced cleavage of solubilizing side chains, to produce hierarchically ordered structures (HOS). The thermal processing can enhance the intermolecular π–π stacking of polyfluorene backbone for better electrochemical performance. Notably, HOS-PFO demonstrated a substantial 6–7 orders of magnitude enhancement in electronic conductivity, showcasing its potential as a functional binder for lithium-ion batteries. As an illustration, HOS-PFO protected SiOx anodes, utilizing in situ side chain decomposition of PFO surrounding SiOx particles after aqueous processing are fabricated. HOS-PFO contributed to the stable cycling and high-capacity retention of practical SiOx anodes (3.0 mAh cm -2 ), without the use of any conducting carbon additives or fluorinated electrolyte additives. It is proposed that this technique represents a universal approach for fabricating electrodes with conjugated polymer binders from aqueous solutions without compromising conductivity.

25 ENERGY STORAGE↗
DOE OSTI2025

Minimal implicit-solvent coarse-grained simulation of Pluronic block copolymers with ionic liquids

Pluronic block copolymers, composed of poly(ethylene oxide) (PEO) and poly(propylene oxide) (PPO) in a triblock structure (PEO–PPO–PEO), are well known for their amphiphilic character and ability to self‐assemble into micelles in aqueous solution. The addition of ionic liquids (ILs) can further modulate the core–shell structures of these copolymers, influencing their stability, critical micellization temperature, and size. However, fully atomistic simulations often become prohibitively expensive due to the size and complexity of these systems. In this work, coarse‐grained simulations using a minimal implicit‐solvent model were performed to examine how two classes of ILs, namely, 1‐alkyl‐3‐methylimidazolium ([C n C 1 im]) and 1‐alkyl‐3‐methylpyrrolidinium ([C n C 1 pyrr]), change the micellization of Pluronic block copolymers in aqueous solution. The effects of IL concentration and alkyl group length were investigated, and the model greatly improved the efficiency of simulating large‐scale micelle systems. Furthermore, the numerical simulations are qualitatively compared with experimental investigations. Our results show that adding ILs expands the micelle core by embedding IL tails among the PPO blocks, thereby increasing overall micelle size. Less polar ILs generally induce more pronounced micellar growth. However, the effect of IL tail length on conformation and micellar packing is non‐monotonic. Up to moderate chain lengths (around C8–C10), the IL tails can extend sufficiently to increase local separation within the micelle; at longer tail lengths, enhanced hydrophobic clustering and steric hindrance cause the tails to bend or fold, capping further expansion. In addition, although block copolymer chains tend to pack more closely in the presence of longer‐tailed ILs, the random coil size of an individual polymer chain does not necessarily shrink. Meanwhile, these insights provide a deeper understanding of how Pluronic/IL systems interact, informing applications in drug delivery, cosmetics, food, and environmental engineering. Finally, our minimal implicit‐solvent model can be applied to larger systems and longer timescales, substantially reducing computational cost while reproducing key structural trends observed experimentally.

Atomistic simulations↗
DOE OSTI2020

Mineral-specific heterogeneous neptunium sorption onto geological repository rocks in oxic and anoxic conditions and various temperatures

Neptunium poses challenges for radioactive wastes management due to its long half-life and ingrowth from americium decay, as well as its high mobility in water in its oxidized forms. Np(V)-bearing aqueous solutions were contacted with polished and crushed gneiss and dolerite from the Yeniseisky site, Russia, which is under consideration as a location for an underground repository, both in oxic and anoxic conditions at 20 and 90 °C. Neptunium sorption on crushed gneiss and dolerite in oxic conditions at 20 °C was quite low, but increased substantially in anoxic conditions. Autoradiography on samples treated at 20 °C revealed that sorbed Np was associated mostly with dark colored iron-bearing minerals, and solvent extraction indicated partial reduction of Np(V) to Np(IV). Sorption of Np from aqueous solutions contacted with dolerite and gneiss at 90 °C was in all cases higher than at the lower temperature. Autoradiography of polished rock samples treated in anoxic conditions at 90 °C revealed that Fe-bearing minerals had 3–4 orders of magnitude more Np sorbed on them than the quartz and feldspar grains. X-ray photoelectron spectroscopy indicated this sorbed Np was present as NpO 2 . This study supports the concept that the geologic conditions present at the planned repository site will retard Np migration.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Vaterite Optical Petrography in Lake Sturgeon Otoliths

Vaterite occurring in fish otoliths exhibits plumose and spherulitic habits, the latter being like those grown from aqueous solutions. Vaterite spherulites are known to have sheaf-like texture and can be up to 100 μm in size. In thin section, the spherulites typically show uniaxial positive interference figures between cross-polarizers using a polarizing light microscope but without conoscopic or Bertrand lenses because the spherulites mimic the effect of these lenses. Plumose vaterite can show both uniaxial (+) and biaxial (+) interference figures with a small 2V, which suggests that crystal structural models for vaterite should have the plane of the carbonate group parallel to the optic axis or acute bisectrix. Finally, vaterite spherulites grown from aqueous solution exhibit similar structure and growth habits.

58 GEOSCIENCES↗
DOE OSTI2024

Modeling and Optimization of Zeolites for Contaminant Removal from Coal Combustion Impoundment Leachates

Coal Combustion Residual (CCR) landfills in the U.S. can lead to elevated contaminant concentrations in groundwater and runoff, particularly of arsenic, boron, and selenium. Zeolites can be readily synthesized from materials contained within the coal combustion ash impoundments and can be cation-exchanged to enhance adsorption capacity, selectivity, and reactivity. However, the optimization of zeolites becomes a daunting task when accounting for the variety of Si:Al ratios, the species of extra-framework cations present, and the zeolite pore structure. Molecular simulations provide methods to study and guide the design of zeolites for the sorption of contaminants from aqueous solution. In this work, models that predict the sorption of arsenic, boron, and selenium from water using cation-exchanged zeolites. Because of the lack of experimental adsorption data for these species of contaminants, models were parametrized to reproduce the results of first-principles calculations and then used to predict sorption for zeolites for a dataset containing 6000 combinations of zeolites and sorption conditions. Machine learning was used to train a model to predict sorption for materials in this database based on the results of the molecular simulations. Next, a genetic algorithm was used to optimize zeolites for the removal of each contaminant from aqueous solution for individual impoundment sites based on contaminant concentrations reported by the Electric Power Research Institute.

Findley, John↗
DOE OSTI2020

Leveraging Actinide Hydrolysis Chemistry for Targeted Th and U Separations using Amidoxime-Functionalized Poly(HIPE)s

Polymerized high internal phase emulsions (poly(HIPE)s) are porous polymer monoliths whose synthesis can easily be tailored to allow incorporation of functional units. In this work, nitrile containing poly(HIPE)s have been prepared with either acrylonitrile (AN) or 4-cyanostyrene (4CS) comonomers. Post-synthetic modification of these nitrile-containing poly(HIPE)s yields their corresponding amidoximated analogues, which were studied for actinide uptake. These amidoxime-functionalized, porous polymers were shown to adsorb 95 % Th4+ species from aqueous solution within 30 minutes. In contrast to other amidoxime containing polymers the uptake of UO22+ in these poly(HIPE)s is lower under similar conditions. A critical analysis of actinide separations and high-energy X-ray scattering data provides insight into the polymers’ selectivity, enabled by the uptake of multinuclear Th clusters.

actinides↗
DOE OSTI2021

Potential Dependent Mn Oxidation and Its Role in Passivation of Ni 38 Fe 20 Cr 22 Mn 10 Co 10 Multi-Principal Element Alloy Using Multi-Element Resolved Atomic Emission Spectroelectrochemistry

The identity of passivating oxides on multi-principal element alloys is of great interest as their optimization offers the potential for exceptional corrosion resistance in aqueous solutions over a broad range of potential and pH. This study focuses on a non-equiatomic Ni 38 Fe 20 Cr 22 Mn 10 Co 10 solid solution alloy and tracks the fate of each alloying element during linear sweep voltammetry, low and intermediate potential holds in the passive potential domain as well as during open circuit relaxation after anodic polarization in slightly acidified Cl − solution. Ni dissolves at all potentials investigated in this work, Fe and Co are incorporated into oxides or hydroxides in low concentrations whilst Cr and Mn are enriched at passive potentials. At low passivating potentials, Mn(II) dissolves and is incorporated in minor amounts in oxides containing large concentrations of Cr(III). Considerable enrichment in Mn(II)-species occurs relative to Cr(III) in the oxide at 0.1 V vs SCE. Electrochemical impedance spectroscopy suggests the presence of layered oxides with marginal passivation at high Mn(II) levels. The formation of these oxides depends on a combination of thermodynamic and kinetic factors as well as the sequence of passivation.

Han, Junsoo (ORCID:0000000192952954)↗
DOE OSTI2023

On electron beam-induced degradation of vinyl ester thermosets

We have demonstrated that electron beam radiolysis induces scissions of the C-O-C bonds along the backbone of the chains of unsaturated polyester thermosets of different compositions based on dicyclopentadiene, isophthalic acid, epoxy vinyl ester, and terephthalic acid. The radiolysis is imminent irrespective of the degree of crosslinking in the thermosets both in neat resins and in the presence of solvents. Electron Paramagnetic Resonance (EPR) results show the formation of the alkoxyl radicals and C-centered radicals as the primary intermediate products of the C-O-C scissions. While the alkoxyl radicals of these resins exhibit very good stability even six months after the irradiation, the C-centered radicals decay very rapidly via their reactions with oxygen that is available either as adsorbed in the resins, dissolved in solvents or from the environment. The radiolytically produced •OH radicals in the unsaturated-ester aqueous solutions play a major role in inducing scissions on the backbone of the polymer chains. The solvated electrons (e s ) from organic solvents, such as dimethyl sulfoxide and isopropyl alcohol, also induce direct scission of the C-O-C bonds, giving rise to the formation of alkoxyl radicals and C-centered radicals. However, a considerable fraction of es is scavenged by the dissolved O 2 to produce O 2 - . Despite the radiation-induced scissions, irradiation of the resins at a dose level of 1000 kGy results in an increase of the glass transition temperature, Tg. This is due to the simultaneous radiation-induced polymerization of the vinyl monomers, toluene, and styrene that are present in the resins. The increase in Tg is observed in all the resins except for the terephthalic polyester resin in aqueous solution, in which the absence of these monomers results in Tg decreasing sharply with irradiation. This work demonstrates that ionizing radiation triggers continuous free radical-chain reactions that lead to the formation of recyclable oligomers.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

Probing the influence of ion-pairing on ligand-field excited-state dynamics

Exploration of the photophysical and photochemical properties of transition metal complexes has driven ground-breaking advancements in solar energy conversion technologies, including photoredox catalysis. While significant research has been devoted to understanding excited state properties of second- and third-row transition metal complexes, earth-abundant first-row metal complexes have received comparatively little attention in this context until very recently. In particular, the role of ion-pairing – which has been identified as a potentially significant factor for Ir(III)-based photosensitizers – has not been examined with regard to its influence on the ligand-field excited states that dominate much of first-row photophysics. A key challenge in studying ion-pair interactions lies in quantifying the extent and nature of ion-pairing, particularly in non-aqueous media where the vast majority of photophysical studies are performed. Cobalt(III) polypyridyl complexes provide an attractive platform to address such questions due to their demonstrated potential for applications in photoredox catalysis involving ligand-field excited states. In the present study, we prepared a cobalt(III) polypyridyl complex, [Co(4,4′-OMebpy) 3 ](BAr F 4 ) 3 (where 4,4′-OMebpy is 4,4′-dimethoxy-2,2′-bipyridine and BAr F 4 is tetrakis(3,5-bis(trifluoromethyl)-phenyl)borate) to probe ion-pairing in non-aqueous solutions. Specifically, analysis of data acquired from both variable-temperature diffusion ordered spectroscopy (DOSY) NMR and 1-D rotating-frame nuclear Overhauser effect (ROE) experiments allowed us to identify and differentiate between solvent-separated ion pairs in high-dielectric media and contact ion pairs in a low-dielectric solvent. Time-resolved absorption spectroscopy was then used to measure ground-state recovery dynamics under these varying conditions of ion-pairing, the results of which revealed an increase in excited-state lifetime for contact ion-pairs that we suggest arises from a reduction in outer-sphere reorganization energy relative to conditions which favored solvent-separated ion pairs. We believe this study demonstrates that one can leverage broadly available NMR-based methods to understand ion-pairing in non-aqueous solutions, which in turn can provide a microscopic picture of intermolecular interactions that can impact the photophysical properties of transition metal-based chromophores.

Ghosh, Atanu [Michigan State University, East Lans↗
DOE OSTI2024

Insights into the structure and dynamics of K + ions at the muscovite–water interface from machine learning potential simulations

The surfaces of many minerals are covered by naturally occurring cations that become partially hydrated and can be replaced by hydronium or other cations when the surface is exposed to water or an aqueous solution. These ion exchange processes are relevant to various chemical and transport phenomena, yet elucidating their microscopic details is challenging for both experiments and simulations. Here, in this work, we make a first step in this direction by investigating the behavior of the native K + ions at the interface between neat water and the muscovite mica (001) surface with ab-initio -based machine learning molecular dynamics and enhanced sampling simulations. Our results show that the desorption of the surface K + ions in pure ion-free water has a significant free energy barrier irrespective of their local surface arrangement. In contrast, facile K + diffusion between mica’s ditrigonal cavities characterized by different Al/Si orderings is observed. This behavior suggests that the K + ions may favor a dynamic disordered surface arrangement rather than complete desorption when exposed to deionized water.

Ab-initio molecular dynamics↗
DOE OSTI2025

Ambient mechanosynthesis of flexible two-dimensional covalent organic frameworks

We present the first ambient mechanosynthesis of 16 flexible covalent organic frameworks (COFs) within an hour. Notably, one representative COF exhibited a high iodine uptake capacity of ∼4.3 g g −1 from aqueous solutions and 5.97 g g −1 from vapor. Flexible two-dimensional covalent organic frameworks (2D COFs) constructed from nonplanar building blocks represent an emerging paradigm in COF design. Nevertheless, the prevailing solvothermal synthesis suffers from low time efficiency, environmental unfriendliness, and cumbersome protocols. Here, we address these challenges by developing the first ambient mechanosynthesis of a diverse library of flexible 2D COFs. Sixteen distinct triazine-cored Schiff-base COFs, including five as-yet-unreported ones, were rapidly synthesized via ball milling using 2,4,6-tris(4-aminophenoxy)-1,3,5-triazine (TPT-NH 2 ) and 2,4,6-tris(4-formylphenoxy)-1,3,5-triazine (TPT-CHO) as building blocks. Notably, the representative COF, MC-flexible-COF-1, was synthesized in as little as one hour under mechanosynthesis conditions, whereas it remained unattainable via the traditional solvothermal method, despite prolonged heating and extensive solvent screening. The highly dynamic nature of the imine linkage was unequivocally demonstrated through mechanochemical “scrambling” experiments using molecular model compounds. Furthermore, MC-flexible-COF-1 exhibited a high iodine uptake capacity of ∼4.30 g g −1 from aqueous solutions and 5.97 g g −1 from the vapor phase. This work underscores the immense potential of mechanochemistry as a powerful and sustainable tool for the rapid synthesis of advanced 2D COFs, including those inaccessible via conventional solution-based methods.

Nailwal, Yogendra↗
DOE OSTI2025

Illuminating the mechanistic impacts of an Fe-quaterpyridine functionalized crystalline poly(triazine imide) semiconductor for photocatalytic CO 2 reduction

The strategy of incorporating earth-abundant catalytic centers into light-absorbing architectures is desirable from the viewpoint of low cost, low toxicity, and versatility for activating small molecules to produce solar-based fuels. Herein, we show that an Fe-quaterpyridine molecular catalyst can be anchored to a light-absorbing, crystalline, carbon nitride (PTI), to yield a molecular-catalyst/material hybrid, Fe-qpy-PTI, capable of facilitating CO 2 reduction to CO selectively (up to ∼97–98%) in aqueous solution under low-intensity light irradiation. This hybrid material leverages the ability of the Fe-qpy catalyst to bind CO 2 upon a one-electron reduction, as achieved by transfer of excited electrons from the carbon-nitride semiconductor. At a low incident power density of only 50 mW cm −2 , the catalytic activity of the hybrid material was measured across a range of catalyst loadings from 0.1–3.8 wt%, yielding CO rates of up to 596 μmol g −1 h −1 for a 3.8 wt% loading during a 3 h experiment. Over the course of 8 h, the hybrid material attained a CO evolution rate of 608 μmol g −1 h −1 and 305 turnovers for a TOF of ∼38 h −1 and an apparent quantum yield of 2.6%. Higher light intensities provided an initial increase in activity but negatively impacted photocatalytic rates with time, with an AQY of 0.6% at 150 mW cm −2 and 0.4% at 250 mW cm −2 . Transient absorption spectroscopy results showed electron survival probabilities consistent with the trends in observed product rates. Computational modeling was also used to evaluate and understand the mechanistic pathway of the high product selectivity for CO versus H 2 . These results thus help unveil key factors for leveraging the mechanistic understanding of molecular catalysts for CO 2 reduction for pairing with light absorbing semiconductors and establishing optimal conditions to attain maximal rates in aqueous solution.

McGuigan, Scott [Baylor University, Waco, TX (Unit↗
DOE OSTI2021

Oxidation and Associated Pore Structure Modification During Experimental Alteration of Granite

Weathering plays a crucial role in a number of environmental processes, and the microstructure and evolution of multi-scale pore space is a critically important factor in the weathering of geological formations. In igneous rocks the infiltration of meteoric water into initially relatively dry material can cause the rock to begin to disaggregate, increasing porosity and surface area, and allowing further disaggregation and weathering to occur. These processes, in turn, allow biota to colonize the rock, further enhancing the weathering rate. In some igneous rocks this may be driven by oxidation of primary minerals. One such mineral, biotite, has been repeatedly mentioned as a cause of cracking during oxidation. However, the scale-dependence of the processes by which this occurs are poorly understood. We cannot, therefore, accurately extrapolate laboratory reaction rates to the field in predictive numerical models. In order to better understand the effects of oxidation and test the hypothesis that fracture and disaggregation are initiated by swelling of oxidizing biotites, we reacted granite cores in a selenic acid-rich aqueous solution at 200°C for up to 438 days. Elevated temperatures and selenic acid were used to provide relatively fast reaction rates and highly oxidizing conditions in sealed reaction vessels. These experiments were analyzed using a combination of imaging, X-ray diffraction, Mössbauer spectroscopy, and small- and ultra-small angle neutron scattering to interrogate porosity and microfracture formation. The experimental results show little observable biotite swelling, but significant transport and growth of iron oxides and/or clays along grain boundaries throughout the sample. Significant increases in porosity were also observed at the sample rim, likely associated with feldspar alteration. Fractures and transport were observed throughout the core, suggesting that stresses due to crystallization pressures caused by the growing iron phases may be the initiating factors in granite weathering, possibly followed by biotite swelling after sufficient permeability is achieved.

Anovitz, Lawrence M.↗
DOE OSTI2024

Liquid Interfacial Electron Microscopy Identifies Nanogalvanic Corrosion in Pearlitic Steel

The nanoscale mechanisms of localized corrosion in low carbon steels have remained elusive due to the complexity of studying the degradative material behavior at nanoscale solid-liquid interfaces. We identified various steps in the nanogalvanic corrosion processes using in-situ liquid-cell scanning transmission electron microscopy (STEM) using a microfluidic holder by Hummingbird Scientific. Initial work, performed at low magnification, identified the initiation point on a 1018 low-carbon steel surface. This initiation point was determined to be a triple junction of two ferrite grains bridging a cementite grain in contact with a baseline electrolyte of 6 uM CO2 dissolved in a buffered (2.78 uM Na2SO4) aqueous solution, pH 6.1. The pre-etched low-carbon steel surface was prepared using focused ion beam lift-out procedures to extract a cross-section of the low-carbon steel surface, which then was thinned to about 150 nm and transferred to a SiN membrane microfluidic window. The transfer was made using a lift-out needle to attach the low-carbon steel lamella to the corner of the SiN window, and then Pt/C deposition held the lamella in contact with the window while it was released from the lift out needle. To identify the triple point on the low carbon steel lamella, prior to attachment on the SiN window, the sample was characterized for compositional variations with energy dispersive x-ray spectroscopy mapping, grain orientation and phase mapping with precession electron diffraction, and thickness mapping with energy filtered transmission electron microscopy. This pre-characterization prior to the in-situ experiment provided a map of the multiphase and multigrain structure, where the in-situ liquid cell imaging provided a clear understanding of the initiation point on the sample. These data were cross-correlated to paint a holistic picture of the triple junction site, enabling low electron-fluence in-situ snapshot imaging to avoid dominating the native corrosion reactions with effects from the incident electron beam. This initial result identified that localized, nanogalvanic corrosion at the phase interface was the dominant corrosion process in the low-carbon steel, so we next targeted the observation of an array of these nanogalvanic features phase boundaries in a pearlite grain. Near-surface ferrite/cementite phase interfaces that typify pearlitic low-carbon steel were extracted, pre-characterized, and imaged for the in-situ corrosion processes. The sample was a cross-section from a pearlite grain, with alternating ferrite and cementite grains that extended microns down from the pre-etched low-carbon steel pipe surface. After contact with a buffered aqueous solution, the phase boundaries between the ferrite and cementite began to dissolve, with observable material loss and thickness changes in the dark-field and bright-field STEM images. Within minutes, the corrosion front proceeded deeper into the material, claiming a thin layer of ferrite around all exposed phase boundaries before progressing laterally into the ferrite matrix, converting the ferrite to corrosion product normal to each buried cementite grain. Formation of the corrosion product causes a volumetric expansion, creating a lateral wedging force that mechanically ejects the cementite grains from their grooves and leaves behind percolation channels into the steel substructure. Rapid and deleterious, this nanogalvanic corrosion pathway represents an important target for understanding and preventing run-away degradation in this common building material. Observation of this corrosion mechanism was enabled by the combination of pre-characterization using standard structural, grain, and compositional analysis in the TEM, which provides maps for understanding the reaction propagation captured in low-dose, in-situ, liquid-cell STEM.

corrosion↗
DOE OSTI2020

Influence of Proton Activity Gaps between Electrodes on Open-Circuit Potential of H 2 /H 2 and H 2 /Air Cells

Polymer -electrolyte-fuel-cell open-circuit voltages (OCVs) are exactly defined by equation (1), where cathode and anode proton activities [(aH+)cathode and (aH+)anode, respectively] usually are identical, so the third term in the right-hand side of equation (2) is ignored. OCV=E0+RT/2F*ln(a1/2 O2*(a2 H+)cathode/aH2O)-RT/2F*ln((a2 H+)anode/aH2) (1) =E0+RT/2F*ln(a1/2 O2*aH2/aH2O)+RT/2F*ln((a2 H+)cathode/(a2 H+)anode) (2) Water vapor pressure is a colligative property that fundamentally correlates to electrolyte concentrations in aqueous solutions. Proton activity is a function of acid concentration, such as pH, when electrolytes are acids. In polymer-electrolyte membranes, water vapor pressure and acid concentration are understood as relative humidity (RH) and water uptake (λ), respectively, where λ represents number of water molecules per sulfonic acid molecule. Several investigations have reported the relationship between RH and λ, meaning that proton activities and associated water uptakes are intimately related to RH. In actual fuel cell operation, cathode RH is determined by ambient-atmosphere and/or humidifier RH(s), and anode RH depends on hydrogen-circulator RH. Therefore, RH is not always identical at both electrodes, and the difference between electrode RHs is considerable during dry operation of polymer electrolyte fuel cells. Therefore, the third term in the right-hand side of equation (2) may be significant for dry operation. We measured OCVs when hydrogen was supplied to both electrodes at 80°C. One electrode (A) was fixed at 30% RH, while RH at the other electrode (B) was varied (0, 5, 10, 20, and 30%). Measured OCVs varied from 0 to 75 mV. For fuel cell tests, electrode A was supplied with hydrogen at 30% RH; electrode B, oxygen at 0, 5, 10, 20, and 30% RH. OCVs deviated from that measured when RH at electrode B was 30%, increasing from 0 to 60 mV with decreasing RH at electrode B. Results are also shown in Figure 1. Proton activities of both electrodes were thermodynamically calculated. The Gibbs–Duhem relation was applied to obtain molar Gibbs free energies of water and sulfonic acid, and proton activity coefficient was calculated using the Gibbs free energy of sulfonic acid and the relationship between RH and λ1–4, assuming that protons and sulfonic anions show identical ionic-activity coefficients. OCVs were estimated using the third term in the right-hand side of equation (2). Results are shown in Figure 1. Fuel-cell current–voltage performance was poor when RHs at the anode and cathode were 30 and 20%, respectively. To determine kinetic current, we measured the oxygen-reduction reaction (ORR) using a rotating-disk electrode (RDE) in concentrated-acid aqueous solutions, which modeled catalyst-layer ionomers. Kinetic currents decreased with acid concentrations. References T. A. Zawodzinski, Jr., C. Derouin, S. Radzinski, R. J. Sherman, V. T. Smith, T. E. Springer and S. Gottesfeld , J. Electrochem. Soc., 140,1041 (1993) P. K. Das and A. Z. Weber, Proceedings of the ASME 2013 11th Fuel Cell Science, Engineering and Technology Conference, Fuel Cell 18010 (2013) V. A. Sethuraman, J. W. Weidner, A. T. Haug, S. Motupally,b and L. V. Protsailo, J. Electrochem. Soc., 155, B50 (2008) A. Kusoglu and A. Z. Weber, Chem. Rev., 117, 987 (2017) Figure 1

Yoshida, Toshihiko↗
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↗