Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “oxygen diffusion”

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,225 records · Page 68

Exploring an ‘ideal hill': how lithology and transport mechanisms affect the possibility of a steady state during weathering and erosion

We introduce a model of chemical reaction within hills to explore how evolving porosity (and by inference, permeability) affects flow pathways and weathering. The model consists of hydrologic and reactive-transport equations that describe alteration of ferrous minerals and feldspar. These reactions were chosen because previous work emphasized that oxygen- and acid-driven weathering affects porosity differently in felsic and mafic rocks. A parameter controlling the order of the fronts is presented. In the absence of erosion, the two reaction fronts move at different velocities and the relative depths depend on geochemical conditions and starting composition. In turn, the fronts and associated changes in porosity drastically affect both the vertical and lateral velocities of water flow. For these instances, estimates of weathering advance rates based on simple models that posit unidirectional constant-velocity advection do not apply. In the model hills, weathering advance rates diminish with time as the Darcy velocities decrease with depth. The system can thus attain a dynamical steady state at any erosion rate where the regolith thickness is constant in time and velocities of both fronts become equal to one another and to the erosion rate. The slower the advection velocities in a system, the faster it attains a steady state. For example, a massive rock with relatively fast-dissolving minerals such as diabase -- where solute transport across the reaction front mainly occurs by diffusion -- can reach a steady state more quickly than granitoid rocks in which advection contributes to solute transport. The attainment of a steady state is controlled by coupling between weathering and hydrologic processes that force water to flow horizontally above reaction fronts where permeability changes rapidly with depth and acts as a partial barrier to fluid flow.

54 ENVIRONMENTAL SCIENCES↗

Dividing the Concentrator Target From the Genesis Mission

The Genesis spacecraft, launched in 2001, traveled to a Lagrangian point between the Earth and Sun to collect particles from the solar wind and return them to Earth. However, during the return of the spacecraft in 2004, the parachute failed to open during descent, and the Genesis spacecraft crashed into the Utah desert. Many of the solar wind collectors were broken into smaller pieces, and the field team rapidly collected the capsule and collector pieces for later assessment. On each of the next few days, the team discovered that various collectors had survived intact, including three of four concentrator targets. Within a month, the team had imaged more than 10,000 fragments and packed them for transport to the Astromaterials Acquisition and Curation Office within the ARES Directorate at JSC. Currently, the Genesis samples are curated along with the other extraterrestrial sample collections within ARES. Although they were broken and dirty, the Genesis solar wind collectors still offered the science community the opportunity to better understand our Sun and the solar system as a whole. One of the more highly prized concentrator collectors survived the crash almost completely intact. The Genesis Concentrator was designed to concentrate the solar wind by a factor of at least 20 so that solar oxygen and nitrogen isotopes could be measured. One of these materials was the Diamond-on-Silicon (DoS) concentrator target. Unfortunately, the DoS concentrator broke on impact. Nevertheless, the scientific value of the DoS concentrator target was high. The Genesis Allocation Committee received a request for approximately 1 cm(sup 2) of the DoS specimen taken near the focal point of the concentrator for the analysis of solar wind nitrogen isotopes. The largest fragment, Genesis sample 60000, was designated for this allocation and needed to be precisely cut. The requirement was to subdivide the designated sample in a manner that prevented contamination of the sample and minimized the risk of losing or breaking the precious requested sample fragment. The Genesis curator determined that the use of laser scribing techniques to "cut" a precise line and subsequently cleave the sample (in a controlled break of the sample along that line) was the best method for accomplishing the sample subdivision. However, there were risks, including excess heating of the sample, that could cause some of the implanted solar wind to be lost via thermal diffusion. Accidentally breaking the sample during the handling and cleaving process was an additional risk. Early in fiscal year 2013, to address this delicate, complicated task, the ARES Directorate assembled its top scientists to develop a cutting plan that would ensure success when applied to the actual concentrator target wafer; i.e., to produce an approximately 1 cm(sup 2) piece from the requested area of the wafer. The team, subsequently referred to as the JSC Genesis Tiger Team, spent months researching and testing parameters and techniques related to scribing, cleaving, transporting, handling, and holding (i.e., mounting) the specimen. The investigation required considerable "thinking outside the box," and many, many trials using nonflight wafer analogs. After all preliminary testing, the following method was adopted as the final cutting plan. It was used in two final end-to-end practice runs before being used on the actual flight target wafer. The wafer was oriented on the laser cutting stage with the 100 and 010 directions of the wafer parallel to the corresponding X and Y directions of the cutting stage. The laser was programed to scribe 31 lines of the appropriate length along the Y stage direction. The programed scribe lines were separated by 5 micron in the X direction. The laser parameters were set as follows: (1) The laser power was 0.5 watts; (2) each line consisted of 50 passes, with the Z position being advanced 5 micron per pass; and (3) 30 s would elapse before the next line was scribed to allow for wafer cool down from any possible heating via the laser. The ablated material that "stuck" in the "scribe-cut" was removed from the "cut" using an ultrasonic micro-tool. After all the ablated silicon was removed from the wafer, the wafer was repositioned in exactly the same orientation on the laser stage. The laser was focused using the bottom of the wafer channel, and the 31-line scribing pattern described above was reprogrammed using the Z position of the groove bottom as the starting Z value instead of the top wafer surface, which was used previously. Upon completion of the second set of scribes, the ultrasonic micro-tool was again used to clean out the cut. The wafer was remounted on the stage in exactly the same orientation as before. The laser was again focused on the bottom of the groove. This time, however, the laser was.programed to scribe only one line down the exact center of the channel. The final scribe line consisted of 100 passes with a Z advance of 5 micron per pass and with the laser power set at 0.5 watts. As mentioned above, the final cutting plan was practiced in two end-to-end trials using non-flight, triangular-shaped silicon wafers similar in size and orientation to the actual DOS 60000 target sample. The actual scribing of the triangular-shaped wafers required scribing two lines and cleaving (i.e. scribe-cleave, then scribe-cleave) to obtain the piece requested for allocation. Early in December 2012, after many months of experiments and practicing and perfecting the techniques and procedures, the team successfully subdivided the Genesis DoS 60000 target sample, one of the most scientifically important samples from the Genesis mission (figure 2). On December 17, 2012, the allocated piece of concentrator target sample was delivered to the requesting principal investigator.The cutting plan developed for the subdivision of this sample will be used as the model for subdividing future requested Genesis flight wafers (appropriately modified for different wafer types).

Lauer, H. V., Jr.↗

Structure of confined laminar spray diffusion flames: Numerical investigation

The structure of confined laminar spray diffusion flames is investigated numerically by solving the gas-phase conservation equations for mass species, continuity, momentum, and energy and the liquid-phase equations for droplet position, velocity, size, and temperature. A one-step global reaction scheme along with six equilibrium reactions are employed to model the flame chemistry. Monodisperse as well as polydisperse sprays are considered. The numerical results demonstrate that liquid spray flames substantially differ from gaseous flames in their structure, i.e., temperature, concentration, and velocity fields, shape, and dimensions under the same conditions. Spray flames are predicted to be taller and narrower than their counterpart gaseous ones and their shapes are almost cylindrical. This is in agreement with experimental observations. The numerical computations also show that the use of the equilibrium reactions with the one-step reaction scheme decreases the flame temperature compared to the one-step reaction scheme without the equilibrium reactions and more importantly increases the surface area of the flame zone due to a phenomenon termed 'equilibrium broadening.' The spray flames also possess a finite thickness with minimal overlap of the fuel and oxygen species. A case for which a fuel-mixture consisting of 20 to 80 percent gas-liquid by mass is introduced into the combustor is also investigated and compared with predictions using only gaseous or liquid fuel.

Mawid, M. A.↗

Windows(Registered Trademark)-Based Software Models Cyclic Oxidation Behavior

Oxidation of high-temperature aerospace materials is a universal issue for combustion-path components in turbine or rocket engines. In addition to the question of the consumption of material due to growth of protective scale at use temperatures, there is also the question of cyclic effects and spallation of scale on cooldown. The spallation results in the removal of part of the protective oxide in a discontinuous step and thereby opens the way for more rapid oxidation upon reheating. In experiments, cyclic oxidation behavior is most commonly characterized by measuring changes in weight during extended time intervals that include hundreds or thousands of heating and cooling cycles. Weight gains occurring during isothermal scale-growth processes have been well characterized as being parabolic or nearly parabolic functions of time because diffusion controls reaction rates. In contrast, the net weight change in cyclic oxidation is the sum of the effects of the growth and spallation of scale. Typically, the net weight gain in cyclic oxidation is determined only empirically (that is, by measurement), with no unique or straightforward mathematical connection to either the rate of growth or the amount of metal consumed. Thus, there is a need for mathematical modeling to infer spallation mechanisms. COSP is a computer program that models the growth and spallation processes of cyclic oxidation on the basis of a few elementary assumptions that were discussed in COSP: A Computer Model of Cyclic Oxidation, Oxidation of Metals, vol. 36, numbers 1 and 2, 1991, pages 81-112. Inputs to the model include the selection of an oxidation-growth law and a spalling geometry, plus oxide-phase, growth-rate, cycle-duration, and spall-constant parameters. (The spalling fraction is often shown to be a constant factor times the existing amount of scale.) The output of COSP includes the net change in weight, the amounts of retained and spalled oxide, the total amounts of oxygen and metal consumed, and the terminal rates of weight loss and metal consumption.

Smialek, J. L.↗

Microgravity effects on water supply and substrate properties in porous matrix root support systems

The control of water content and water movement in granular substrate-based plant root systems in microgravity is a complex problem. Improper water and oxygen delivery to plant roots has delayed studies of the effects of microgravity on plant development and the use of plants in physical and mental life support systems. Our international effort (USA, Russia and Bulgaria) has upgraded the plant growth facilities on the Mir Orbital Station (OS) and used them to study the full life cycle of plants. The Bulgarian-Russian-developed Svet Space Greenhouse (SG) system was upgraded on the Mir OS in 1996. The US developed Gas Exchange Measurement System (GEMS) greatly extends the range of environmental parameters monitored. The Svet-GEMS complex was used to grow a fully developed wheat crop during 1996. The growth rate and development of these plants compared well with earth grown plants indicating that the root zone water and oxygen stresses that have limited plant development in previous long-duration experiments have been overcome. However, management of the root environment during this experiment involved several significant changes in control settings as the relationship between the water delivery system, water status sensors, and the substrate changed during the growth cycles. c 2001 Published by Elsevier Science Ltd. All rights reserved.

NASA Discipline Life Support Systems↗

Applications of Nickelate perovskites for neuromorphic computing from electronic structure and Machine Learning

While the limit of Moore's law is presently being reached with current microelectronic technologies, we need to develop new paradigms that overcome this limitation. In that respect, neuromorphic computing is a concept that emulates the neural behavior and response of the human brain, and it has been recognized as a promising alternative approach. In this research project, we will perform multi-fidelity scale bridging to explore the potential use of materials with metal to insulator transition for neuromorphic applications. In particular, rare earth nickelates are promising for such purposes, as the transition in these materials is quite sensitive to a broad set of different external stimuli. Our multi-fidelity approach will bridge the high-fidelity electronic structure calculations with classical potentials. We will bridge dynamical mean field theory with a classical atomistic representation via a deep learning force field. The neural network is trained with energies, charges, and forces obtained by accurate electronic structure theories based on Dynamical Mean Field Theory. The configurational space is generated from known crystal phases, ab initio molecular dynamics with exchange-correlation functionals corrected with the Hubbard model, disordered phases with different concentrations of oxygen vacancies, and nonsymmetrical positions and induced strain by grain interfaces or contact with a substrate. Strategies to train the model with a reduced number of training examples are obtained from active learning methods, and new structures for improving the learning process are generated by using machine learning autoencoders. This classical potential will be validated through a diversity of electronic structure methods and represents an important step to combine the flexibility and accuracy of first-principles with the speed of classical potentials. The generated multi-fidelity surrogate model will be used to understand the role of strain, oxygen vacancies, proton doping, the variation of the crystal phase, substrate effects, vibrational effects as the octahedral rotation, grain boundaries and defect effects on the response of a Metal to Insulator Transition (MIT) in correlated materials. Long time and large-scale simulations will help understand the role of different stimuli to control the hysteresis of the MIT, as it has been experimentally suggested. Selected configurations will be analyzed with higher-level theories to provide an accurate electronic description and to study how the orbitals and charges are rearranged under different conditions.

36 MATERIALS SCIENCE↗

CO2 Acquisition Membrane (CAM)

The objective of CAM is to develop, test, and analyze thin film membrane materials for separation and purification of carbon dioxide (CO2) from mixtures of gases, such as those found in the Martian atmosphere. The membranes are targeted toward In Situ Resource Utilization (ISRU) applications that will operate in extraterrestrial environments and support future unmanned and human space missions. A primary application is the Sabatier Electrolysis process that uses Mars atmosphere CO2 as raw material for producing water, oxygen, and methane for rocket fuel and habitat support. Other applications include use as an inlet filter to collect and concentrate Mars atmospheric argon and nitrogen gases for habitat pressurization, and to remove CO2 from breathing gases in Closed Environment Life Support Systems (CELSS). CAM membrane materials include crystalline faujasite (FAU) zeolite and rubbery polymers such as silicone rubber (PDMS) that have been shown in the literature and via molecular simulation to favor adsorption and permeation of CO2 over nitrogen and argon. Pure gas permeation tests using commercial PDMS membranes have shown that both CO2 permeance and the separation factor relative to other gases increase as the temperature decreases, and low (Delta)P(Sub CO2) favors higher separation factors. The ideal CO2/N2 separation factor increases from 7.5 to 17.5 as temperature decreases from 22 C to -30 C. For gas mixtures containing CO2, N2, and Ar, plasticization decreased the separation factors from 4.5 to 6 over the same temperature range. We currently synthesize and test our own Na(+) FAU zeolite membranes using standard formulations and secondary growth methods on porous alumina. Preliminary tests with a Na(+) FAU membrane at 22 C show a He/SF6 ideal separation factor of 62, exceeding the Knudsen diffusion selectivity by an order of magnitude. This shows that the membrane is relatively free from large defects and associated non-selective (viscous flow) transport mechanisms. The Membrane Test Facility (MTF) has been developed to measure membrane permeance over a wide range of temperature and pressure. The facility uses two volume compartments separated by the membrane that are instrumented to measure temperature, delta pressure across the membrane, and gas composition. A thermal shroud supports and encloses the membrane, and provides temperature control. Methods were developed to determine membrane permeance using the first order decay of the pressure difference between the sealed compartments, using the total pressure for pure gases, and partial pressure of each species in gas mixtures. The technique provides an end-to-end measurement of gas permeance that includes concentration polarization effects. Experiments have shown that in addition to membrane permeance properties, the geometry and design of associated structures play an important role in how membrane systems will function on Mars.

Mason, Larry W.↗

Effect of metal encapsulation on bulk superconducting properties of niobium thin films used in qubits

Niobium metal occupies nearly 100% of the volume of a typical 2D transmon device. While the aluminum Josephson junction is of utmost importance, maintaining quantum coherence across the entire device means that pair-breaking in Nb leads, capacitive pads, and readout resonators can be a major source of decoherence. The established contributors are surface oxides and hydroxides, as well as absorbed hydrogen and oxygen. Metal encapsulation of freshly grown surfaces with non-oxidizing metals, preferably without breaking the vacuum, is a successful strategy to mitigate these issues. While the positive effects of encapsulation are undeniable, it is important to understand its impact on the macroscopic behavior of niobium films. We present a comprehensive study of the bulk superconducting properties of Nb thin films encapsulated with gold and palladium/gold, and compare them to those of bare Nb films. Magneto-optical imaging, magnetization, resistivity, and London and Campbell penetration depth measurements reveal significant differences in encapsulated samples. Both sputtered, and epitaxial Au-capped films exhibit the highest residual resistivity ratio and superconducting transition temperature, as well as the lowest upper critical field, London penetration depth, and critical current. These results are in good agreement with the microscopic theory of anisotropic normal and superconducting states of Nb. We conclude that pair-breaking in the bulk of niobium films, driven by disorder throughout the film rather than just at the surface, is a significant source of quantum decoherence in transmons. We also conclude that gold capping not only passivates the surface but also affects the properties of the entire film, significantly reducing the scattering rate due to defects likely induced by surface diffusion if the film is not protected immediately after fabrication.

Datta, Amlan [Ames Lab; Iowa State U.]↗

Mechanical Properties of Carbon Fiber Reinforced Composites Exposed to Cryogenic Conditions and Space Radiation via Simulation and Testing

As NASA missions extend beyond low Earth orbit, increasing reliance is placed on carbon fiber reinforced polymer (CFRP) composites for spacecraft structures where mass efficiency, durability, and long-term reliability are critical. In service, these materials are subjected to a combination of ultraviolet radiation, vacuum, ionizing radiation, atomic oxygen, and extreme thermal excursions under sustained mechanical loading. Flight systems such as the Boeing Starliner and SpaceX Dragon employ external composite structures that will experience these environments for extended durations. Although prior spaceflight and ground studies have reported limited changes in bulk mechanical properties, the synergistic effects of these environments on composite microstructure, particularly at the fiber matrix interphase, remain insufficiently characterized and represent a potential qualification and reliability risk. This study investigates the effects of short-term cryogenic exposure on a radiation shielding carbon epoxy composite, SC2020, as a ground-based analog for space relevant thermal extremes. The SC2020 material system has previously flown on the International Space Station under the Materials International Space Station Experiment (MISSE) program. Composite specimens were exposed to liquid nitrogen for 6 and 24 hours and evaluated using a multiscale characterization framework that combined ASTM D3039 tensile testing, Atomic Force Microscopy (AFM) based interphase analysis, and helium gas permeability measurements. Tensile testing showed no statistically significant or permanent degradation in global strength or modulus following cryogenic exposure. In contrast, AFM measurements revealed reductions in interphase modulus, weakened adhesion, and increased nanoscale heterogeneity, indicating localized degradation mechanisms not captured by conventional bulk testing. Gas permeability measurements showed a progressive increase in helium diffusion with exposure duration, consistent with micro-void formation or partial interfacial debonding. The results indicate that cryogenic exposure initiates degradation at the fiber matrix interphase while leaving global mechanical properties largely unchanged over short durations. These findings underscore the importance of multiscale diagnostics for identifying early-stage damage mechanisms that may influence long term performance and qualification margins for spaceflight composite structures. The data presented establish a cryogenic baseline for comparison with forthcoming MISSE flight exposure results and support ongoing NASA Established Program to Stimulate Competitive Research (EPSCoR) efforts aimed at improving composite qualification methodologies, risk assessment, and reliability prediction for space environments.

composite reliability↗

A Burning Rate Emulator (BRE) for Study in Microgravity

A gas-fueled burner, the Burning Rate Emulator (BRE), is used to emulate condensed-phase fuel flames. The design has been validated to easily measure the burning behavior of condensed-phase fuels by igniting a controlled stream of gas fuel and diluent. Four properties, including the heat of combustion, the heat of gasification, the surface temperature, and the laminar smoke point, are assumed to be sufficient to define the steady burning rate of a condensed-phase fuel. The heat of gasification of the fuel is determined by measuring the heat flux and the fuel flow rate. Microgravity BRE tests in the NASA 5.2 s drop facility have examined the burning of pure methane and ethylene (pure and 50 in N2 balance). Fuel flow rates, chamber oxygen concentration and initial pressure have been varied. Two burner sizes, 25 and 50 mm respectively, are chosen to examine the nature of initial microgravity burning. The tests reveal bubble-like flames that increase within the 5.2s drop but the heat flux received from the flame appears to asymptotically approach steady state. Portions of the methane flames appear to locally detach and extinguish at center, while its shape remains fixed, but growing. The effective heat of gasification is computed from the final measured net heat flux and the fuel flow rate under the assumption of an achieved steady burning. Heat flux (or mass flux) and flame position are compared with stagnant layer burning theory. The analysis offers the prospect of more complete findings from future longer duration ISS experiments.

diffusion flames↗

The use of superoxide mixtures as air-revitalization chemicals in hyperbaric, self-contained, closed-circuit breathing apparatus

In portable breathing apparatus applications at 1 atm, potassium superoxide (KO2) has exhibited low-utilization efficiency of the available oxygen (O2) and diminished carbon dioxide-(CO2) scrubbing capacity caused by the formation of a fused, hydrated-hydroxide/carbonate product coating on the superoxide granules. In earlier work, it was discovered that granules fabricated from an intimate mixture of KO2 and calcium superoxide, Ca(O2)2, did not exhibit formation of a fused product coating and the utilization efficiency with respect to both O2 release and CO2 absorption was superior to KO2 granules when both types of granules were reacted with humidified CO2 under identified conditions. In the work described here, single pellets of KO2, KO2/Ca(O2), mixtures and commercially available KO2 tables and granules were reacted with a flow of humidified CO2 in helium at 1- and 10-atm total pressure and at an initial temperature of 40 C. In the 1-atm flow tests, the reaction rates and utilization efficiency of the KO2/Ca(O2)2 pellets were markedly superior to the KO2 pellets, tablets, and granules when the samples were reacted under identical conditions. However, at 10 atm, the rates of O2 release and CO2 absorption, as well as the utilization efficiencies of all the superoxide samples, were one-third to one-eighth of the values observed at 1 atm. The decrease in reaction performance at 10 atm compared to that at 1 atm has been attributed principally to the lower bulk diffusivity of the CO2 and H2O reactants in helium at the higher pressure and secondarily to the moderation of the reaction temperature caused by the higher heat capacity of the 10-atm helium.

Wood, P. C.↗

Modeling Electrolytic O2 Recovery from Metabolic CO2 for Advanced Closed Loop Life Support Systems in Extraterrestrial Human Missions

The International Space Station (ISS) is currently equipped with a complex, heavy, and power consuming system that recovers approximately 50% of O2 from metabolic CO2. Future long duration human missions to the Moon and Mars will necessitate a sustainable and highly efficient metabolic oxygen recovery system capable of yielding a minimum of 75% O2 recovery. A Macrofluidic Electrochemical Reactor (MFECR) technology development effort is currently underway at NASA Marshall Space Flight Center (MSFC) to significantly increase current metabolic O2 recovery efficiency, expand mission sustainability, and reduce complexity of the system. The novel design combines CO2 conversion to O2 along with C2H4 as byproduct and water electrolysis (currently conducted in two separate units) into a single compact unit that runs at standard conditions and is theoretically capable of generating O2 with a theoretical maximum metabolic CO2 conversion of 73% while consuming less than metabolic water. This paper presents a comprehensive multi-physic 3D model developed at MSFC on CO2 conversion to O2 and C2H4 at standard conditions via MFECR. The 3D spatial domain of the model is a replica of the actual MFECR’s 3D drawing generated for the MFECR fabrication and operated to recover O2 from CO2 yielding C2H4 as byproduct. Electrochemical (EC) physics that includes EC multicomponent reaction mechanisms, mass transport, and electrical current density distributions is coupled in the model with all the other physics phenomena involved in the MFECR’s process, such as two-phase flow, free and porous fluid regimes, multicomponent mass transfer, heat transfer, and DC electrical current generation along with Joule heating effect. The EC reaction sections of the MFECR consists of two porous gas diffusion electrodes (GDE) and an electrolyte serpentine channel sandwiched in the middle. The CO2 feeds the cathode serpentine channel and part of the O2 product is fed back to the anode serpentine chamber. An alkaline solution feeds the electrolyte serpentine chamber wetting the GDEs of both, the anode and cathode allowing the OH- ionic transport between them. The EC reactions in the cathode’s GDE yield C2H4 from CO2 and H2 from water while the EC reaction in the anode’s GDE yields O2. The authors will present in this paper the validation of the model using experimental data and the utilization of the validated model in building a reliable simulator that will not only assist the authors on the MFECR design but also the optimization of its operation in the ISS and future spatial human missions.

Jesus A Dominguez↗

Limits to Creation of Oxygen-Rich Atmospheres on Planets in the Outer Reaches of the Conventional Habitable Zone

Abundant free oxygen appears to be a requirement for macroflora and macrofauna. To the best of our knowledge, a general discussion of which habitable planets are conducive to oxygen has not taken place. Theories for the rise of oxygen fall into 4 categories: (i) It is governed by an intrinsic rate of biological innovation, independent of environmental factors. (ii) It is caused by mantle evolution, probably consequent to secular cooling. (iii) It is caused by hydrogen escape, which irreversibly oxidizes the Earth. (iv) It is Gaia's response to the brightening Sun, its rise prevented until reduced greenhouse gases were no longer needed to maintain a clement climate. All but the first of these make implicit astronomical predictions that can be quantified and made explicit. Here we address the third hypothesis. In this hypothesis hydrogen escape acts like an hourglass that continues until all relevant reduced mineral buffers have been oxidized (titrated, as it were) and the surface made safe for O2. The hypothesis predicts that abundant free O2 will be absent from habitable planets that have not experienced significant hydrogen escape. Where hydrogen escape is modest or insignificant, the atmosphere can be approximated as hydrostatic, which makes assessing radiative cooling by embedded molecules, atoms, and ions such as CO2 and H3+ straightforward. In particular, H2 is efficient at exciting non-LTE CO2 15 micron emission, which makes radiative cooling very effective when H2 is abundant. We can therefore map out the region of phase space in which habitable planets do not lose hydrogen, and therefore do not develop O2 atmospheres. A related matter is the power of radiative cooling by embedded molecules to enforce the diffusion limit to hydrogen escape. This matter in particular is relevant to addressing the empirical observation that rocky planets with thin or negligible atmospheres are rarely or never bigger than approx.1.6 Earth radii.

Zahnle, Kevin↗

A Straightforward Model for Quantifying Local pH Gradients Governing the Oxygen Evolution Reaction

The production and consumption of protons by an electrocatalyst will, under certain conditions, generate localized microenvironments with properties distinct from those of the bulk solution. These local properties are particularly impactful for reactions involving proton-coupled electron transfer, where the generation of locally basic or acidic environments may significantly influence the energy efficiency and reaction selectivity of the electrocatalyst. Whereas local pH environments have been observed and characterized in reductive half-reactions, including the CO 2 reduction and hydrogen evolution reactions, the incompatibility of conventional techniques and materials has limited studies in oxidative half-reactions, including the oxygen evolution reaction (OER), which provides the reducing equivalents for solar-to-fuels electrolysis. With the straightforward parameters bulk pH, buffer composition and pK a , and mass transport, we develop a model for describing local pH as a function of current density regardless of the microscopic details of the mechanism. Using an acid-stable PbO x OER catalyst, we observe the formation and dissipation of pH gradients during the OER and validate the model with voltammetric and potentiometric studies. Here, the model predicts how local acidic environments can develop over a narrow OER current density window, thus providing further motivation for the development of OER catalysts that are stable to acid, even when operating in basic aqueous conditions. More generally, the model is not restricted to the OER and is useful for determining the onset of local pH gradients for other electrocatalytic reactions that involve the consumption or generation of protons in energy conversion reactions.

Anions↗

An automated cluster surface scanning method for exploring reaction paths on metal-cluster surfaces

Metal-cluster surfaces present a wide variety of unique coordination environments. This complexity makes it difficult to manually probe the surface reactivity of such clusters. Here, we present a simple and automated method to systematically discover reaction pathways on cluster surfaces, based on the automated cluster surface scanning (ACSS) technique for mapping out potential energy surfaces. We showcase our method on 55-atom icosahedral Cu and Ag clusters, where we determine the activation energies of four elementary steps common in heterogeneous catalysis – hydrogen recombination (H* + H* → H 2 * + *), oxygen recombination (O* + O* → O 2 * + *), water formation (OH* + H* → H 2 O* + *), and CO oxidation (CO* + O* → CO 2 * + *) – with density functional theory calculations (DFT-PBE + D3). We show that the ACSS method requires significantly less human effort than the established manually performed climbing-image nudged elastic band (MP + CI-NEB) technique and locates transition states with comparable accuracy (root-mean-squared error of 0.10 eV) and similar computational cost. Rigorous sampling of the potential energy surface with the ACSS method allows one to locate all lowest-energy reaction pathways obtained via the MP+CI-NEB approach, as well as alternative pathways that one may have missed with the MP+CI-NEB approach due to the many possible pathways available on these clusters. The accuracy and efficiency afforded by the ACSS method could enable high-throughput exploration of the diverse reactivity of metal clusters.

36 MATERIALS SCIENCE↗

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↗

Smoothed particle hydrodynamics modeling and analysis of oxide reduction process for uranium oxides

A common kinetic feature for oxide reduction chemical/electrochemical processes is oxygen transport via a porous metallic layer, which has been considered as a rate-determining step for reducing uranium oxides to metallic uranium. Accounting this kinetic behavior must involve the resolution of the moving reactive interface between shrinking oxide and expanding metal phases. This study presents a numerical model using smoothed particles hydrodynamics (SPH) to effectively deal with the evolution of the shrinking core reaction interface and oxygen transport via mass transfer of lithium oxide (Li 2 O) species in multiple mass transfer domains. We successfully validated the proposed model against a theoretical derivation for the oxide reduction process handling a shrinking oxide core with molten salt and metal ash medium on a simple planar geometry. Armed with successful validation results, the model examined a realistic reactant geometry to extend the arguments beyond the one-dimensional analyses, allowing the proposed model to apply to general application scenarios with multi-dimensional geometry. Here, this study demonstrated that the proposed model could simulate and evaluate an arbitrary reaction basket design without iterative experimental trials, which is prohibitive for a scaled high-temperature molten salt study in an inert environment. The construct potentially provides not only deep insights on multiphysics behaviors governing the process dynamics but also a robust framework for evaluating and screening candidate basket designs in the most cost-effective manner.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Variations in the Peak Position of the 6.2 micron Interstellar Emission Feature: A Tracer of N in the Interstellar Polycyclic Aromatic Hydrocarbon Population

This paper presents the results of an investigation of the molecular characteristics that underlie the observed peak position and profile of the nominal 6.2 micron interstellar emission band generally attributed to the CC stretching vibrations of polycyclic aromatic hydrocarbons (PAHs). It begins with a summary of recent experimental and theoretical studies ofthe spectroscopic properties of large (>30 carbon atoms) PAH cations as they relate to this aspect of the astrophysical problem. It then continues with an examination of the spectroscopic properties of a number of PAH variants within the context of the interstellar 6.2 micron emission, beginning with a class of compounds known as polycyclic aromatic nitrogen heterocycles (PANHs; PAHs with one or more nitrogen atoms substituted into their carbon skeleton). In this regard, we summarize the results of recent relevant experimental studies involving a limited set of small PANHs and their cations and then report the results of a comprehensive computational study that extends that work to larger PANH cations including many nitrogen-substituted variants of coronene(+) (C24H12(+)), ovalene(+) (C32H14(+)), circumcoronene(+) (C54H18(+)), and circum-circumcoronene(+) (C96H24(+)). Finally, we report the results of more focused computational studies of selected representatives from a number of other classes of PAH variants that share one or more of the key attributes of the PANH species studied. These alternative classes of PAH variants include (1) oxygen- and silicon-substituted PAH cations; (2) PAH-metal ion complexes (metallocenes) involving the cosmically abundant elements magnesium and iron; and (3) large, asymmetric PAH cations. Overall, the studies reported here demonstrate that increasing PAH size alone is insuEcient to account for the position of the shortest wavelength interstellar 6.2 micron emission bands, as had been suggested by earlier studies. On the other hand, this work reveals that substitution of one or more nitrogen atoms within the interior of the carbon skeleton of a PAH cation induces a significant blueshift in the position of the dominant CC stretching feature of these compounds that is suf6cient to account for the position of the interstellar bands. Subsequent studies of the effects of substitution by other heteroatoms (O and Si), metal ion complexation (Fe(+), Mg(+), and Mg(2+)), and molecular symmetry variation-all of which fail to reproduce the blueshift observed in the PANH cations-indicate that N appears to be unique in its ability to accommodate the position of the interstellar 6.2 micron bands while simultaneously satisfying the other constraints of the astrophysical problem. This result implies that the peak position of the interstellar feature near 6.2 micron traces the degree of nitrogen substitution in the population, that most of the PAHs responsible for the interstellar IR emission features incorporate nitrogen within their aromatic networks, and that a lower limit of 1%-2% of the cosmic nitrogen is sequestered within the interstellar PAH population. Finally, in view of the ubiquity and abundance of interstellar PAHs and the permanent dipoles and distinctive electronic structures of these nitrogen-substituted variants, this work impacts a wide range of observational phenomena outside of the infrared region of the spectrum including the forest of unidentified molecular rotational features and the anomalous Galactic foreground emission in the microwave, and the diffuse interstellar bands (DIBs) and other structure in the interstellar extinction curve in the ulhviolet/visible. These astrophysical ramifications are discussed, and the dipole moments and rotational constants are tabulated to facilitate further investigations of the astrophysical role of nitrogen-substituted aromatic compounds.

Hudgins, Douglas M.↗