Search NASA⌕ Search

SEARCH · Search NASA

Results for “Experimental evolution”

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,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

PROPERTIES OF CARBON–OXYGEN WHITE DWARFS FROM MONTE CARLO STELLAR MODELS

We investigate properties of carbon–oxygen white dwarfs with respect to the composite uncertainties in the reaction rates using the stellar evolution toolkit, Modules for Experiments in Stellar Astrophysics (MESA) and the probability density functions in the reaction rate library STARLIB. These are the first Monte Carlo stellar evolution studies that use complete stellar models. Focusing on 3 M{sub ⊙} models evolved from the pre main-sequence to the first thermal pulse, we survey the remnant core mass, composition, and structure properties as a function of 26 STARLIB reaction rates covering hydrogen and helium burning using a Principal Component Analysis and Spearman Rank-Order Correlation. Relative to the arithmetic mean value, we find the width of the 95% confidence interval to be ΔM{sub 1TP} ≈ 0.019 M{sub ⊙} for the core mass at the first thermal pulse, Δt{sub 1TP} ≈ 12.50 Myr for the age, Δlog(T{sub c}/K) ≈ 0.013 for the central temperature, Δlog(ρ{sub c}/g cm{sup −3}) ≈ 0.060 for the central density, ΔY{sub e,c} ≈ 2.6 × 10{sup −5} for the central electron fraction, ΔX{sub c}({sup 22}Ne) ≈ 5.8 × 10{sup −4}, ΔX{sub c}({sup 12}C) ≈ 0.392, and ΔX{sub c}({sup 16}O) ≈ 0.392. Uncertainties in the experimental {sup 12}C(α,γ){sup 16}O, triple-α, and {sup 14}N(p,γ){sup 15}O reaction rates dominate these variations. We also consider a grid of 1–6 M{sub ⊙} models evolved from the pre main-sequence to the final white dwarf to probe the sensitivity of the initial–final mass relation to experimental uncertainties in the hydrogen and helium reaction rates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Inference finds consistency between a neutrino flavor evolution model and Earth-based solar neutrino measurements

We continue examining statistical data assimilation (SDA), an inference methodology, to infer solutions to neutrino flavor evolution, for the first time using real - rather than simulated - data. The model represents neutrinos streaming from the Sun's center and undergoing a Mikheyev-Smirnov-Wolfenstein (MSW) resonance in flavor space, due to the radially-varying electron number density. The model neutrino energies are chosen to correspond to experimental bins in the Sudbury Neutrino Observatory (SNO) and Borexino experiments, which measure electron-flavor survival probability at Earth. In conclusion, the procedure successfully finds consistency between the observed fluxes and the model, if the MSW resonance - that is, flavor evolution due to solar electrons - is included in the dynamical equations representing the model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Coupling scales for modelling heavy metal vaporization from municipal solid waste incineration in a fluid bed by CFD

Highlights: • A CFD two-scale model is formulated to simulate heavy metal vaporization from waste incineration in fluidized beds. • MSW particle is modelled with the macroscopic particle model. • Influence of bed dynamics on HM vaporization is included. • CFD predicted results agree well with experimental data reported in literature. • This approach may be helpful for fluidized bed reactor modelling purposes. - Abstract: Municipal Solid Waste Incineration (MSWI) in fluidized bed is a very interesting technology mainly due to high combustion efficiency, great flexibility for treating several types of waste fuels and reduction in pollutants emitted with the flue gas. However, there is a great concern with respect to the fate of heavy metals (HM) contained in MSW and their environmental impact. In this study, a coupled two-scale CFD model was developed for MSWI in a bubbling fluidized bed. It presents an original scheme that combines a single particle model and a global fluidized bed model in order to represent the HM vaporization during MSW combustion. Two of the most representative HM (Cd and Pb) with bed temperatures ranging between 923 and 1073 K have been considered. This new approach uses ANSYS FLUENT 14.0 as the modelling platform for the simulations along with a complete set of self-developed user-defined functions (UDFs). The simulation results are compared to the experimental data obtained previously by the research group in a lab-scale fluid bed incinerator. The comparison indicates that the proposed CFD model predicts well the evolution of the HM release for the bed temperatures analyzed. It shows that both bed temperature and bed dynamics have influence on the HM vaporization rate. It can be concluded that CFD is a rigorous tool that provides valuable information about HM vaporization and that the original two-scale simulation scheme adopted allows to better represent the actual particle behavior in a fluid bed incinerator.

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

Alteration and Erosion of Rock Matrix Bordering a Carbonate-Rich Shale Fracture

A novel reactive transport model has been developed to examine the processes that affect fracture evolution in a carbonate-rich shale. An in situ synchrotron X-ray microtomography experiment, flowing CO 2 saturated water through a single fracture mini-core of Niobrara Shale provided the experimental observations for the development and testing of the model. The phenomena observed included the development of a porous altered layer, flow channeling, and increasingly limited calcite dissolution. The experimental observations cannot be explained by models that consider only mineral dissolution and development of an altered layer. The difference between the fracture volume change recorded by the microtomography images and what would be expected from mineral dissolution alone suggest that there is erosion of the altered layer as it develops. The numerical model includes this additional mechanism, with the erosion rate based on the thickness of the altered layer, and successfully captures the evolution of the geochemical reactions and morphology of the fracture. The findings imply that the abundance (with a threshold of approximately 35%) and reactivity of the rapidly reacting mineral control the development and erodibility of the altered layer on the fracture surfaces, and therefore fracture opening.

58 GEOSCIENCES↗
NASA NTRS1999

Synergistic Effects of Physical Aging and Damage on Long-Term Behavior of Polymer Matrix Composites

The research consisted of two major parts, first modeling and simulation of the combined effects of aging and damage on polymer composites and secondly an experimental phase examining composite response at elevated temperatures, again activating both aging and damage. For the simulation, a damage model for polymeric composite laminates operating at elevated temperatures was developed. Viscoelastic behavior of the material is accounted for via the correspondence principle and a variational approach is adopted to compute the temporal stresses within the laminate. Also, the effect of physical aging on ply level stress and on overall laminate behavior is included. An important feature of the model is that damage evolution predictions for viscoelastic laminates can be made. This allows us to track the mechanical response of the laminate up to large load levels though within the confines of linear viscoelastic constitutive behavior. An experimental investigation of microcracking and physical aging effects in polymer matrix composites was also pursued. The goal of the study was to assess the impact of aging on damage accumulation, in ten-ns of microcracking, and the impact of damage on aging and viscoelastic behavior. The testing was performed both at room and elevated temperatures on [+/- 45/903](sub s) and [02/903](sub s) laminates, both containing a set of 90 deg plies centrally located to facilitate investigation of microcracking. Edge replication and X-ray-radiography were utilized to quantify damage. Sequenced creep tests were performed to characterize viscoelastic and aging parameters. Results indicate that while the aging times studied have limited ]Influence on damage evolution, elevated temperature and viscoelastic effects have a profound effect on the damage mode seen. Some results are counterintuitive, including the lower strain to failure for elevated temperature tests and the catastrophic failure mode observed for the [+/- 45/9O3](sub s), specimens. The fracture toughness for transverse cracks increases with increasing temperature for both systems: transverse cracking was completely absent prior to failure in [+/- 45/903](sub s), and was suppressed for [02/903](sub s). No significant effect of damage on aging or viscoelastic parameters was observed.

Brinson, L. Cate↗
DOE OSTI2025

Ensemble cure kinetics network (ECK-Net): A method to derive cure kinetics of thermosetting resin

This paper introduces an Ensemble Cure Kinetics Network (ECK-Net), a neural network (NN)–based framework for modeling the cure kinetics of thermosetting resins within a phenomenological context. ECK-Net replaces traditional analytic models, which require extensive chemical insight and multiple isothermal/non-isothermal experiments, with a data-driven surrogate that maps nonlinear relationships between temperature, degree of cure, and reaction rate from differential scanning calorimetry data. The proposed approach predicts input-dependent kinetic coefficients of a generalized nth-order reaction equation rather than reaction rates directly, enabling a single unified model to represent various epoxy systems without relying on iso-conversional analysis or predefined functional forms. To ensure robustness, multiple independently trained networks under different random initializations are blended through an ensemble strategy, effectively mitigating the stochastic variability inherent to neural networks. The framework is validated using experimental datasets from multiple resin systems, including aerospace-grade materials (Toray 3900-2, Cycom 5320-1, and Hexcel 8552) and a windmill-grade resin (RIMR 035c). The model accurately reproduces the temporal evolution of the degree of cure under manufacturers’ recommended cure cycles across all tested resins systems, yielding Pearson’s correlation coefficients of 0.992, 0.994, 0.993, 0.997, respectively. To demonstrate process-level applicability, the trained network was implemented within the Abaqus environment to simulate out-of-autoclave (OOA) curing process of the CFRP panel composed of Toray T830H-6K/3900-2D prepreg. The simulation results showed excellent agreement with experimental temperature response (maximum peak temperature, simulation: 189.6 °C, experiment: 188.5 °C) and the final degree of cure (simulation: 0.948, experiment: 0.960 ± 0.013), confirming ECK-Net’s capability as a reliable alternative to conventional cure kinetics modeling methods.

Composite curing↗
DOE OSTI2015

Deformation mechanisms in a precipitation-strengthened ferritic super alloy revealed by in situ neutron dffraction studies at elevated temperatures

The ferritic superalloy Fe–10Ni–6.5Al–10Cr–3.4Mo strengthened by ordered (Ni,Fe)AlB2-type precipitates is a candidate material for ultra-supercritical steam turbine applications above 923 K. Despite earlier success in improving its room-temperature ductility, the creep resistance of this material at high temperatures needs to be further improved, which requires a fundamental understanding of the high-temperature deformation mechanisms at the scales of individual phases and grains. In situ neutron diffraction has been utilized to investigate the lattice strain evolution and the microscopic load-sharing mechanisms during tensile deformation of this ferritic superalloy at elevated temperatures. Finite-element simulations based on the crystal plasticity theory are employed and compared with the experimental results, both qualitatively and quantitatively. Based on these interphase and intergranular load-partitioning studies, it is found that the deformation mechanisms change from dislocation slip to those related to dislocation climb, diffusional flow and possibly grain boundary sliding, below and above 873 K, respectively. Insights into microstructural design for enhancing creep resistance are also discussed.

nuclear (including radiation effects), defects, me↗
DOE OSTI2020

Modeling the role of local crystallographic correlations in microstructures of Ti-6Al-4V using a correlated structure visco-plastic self-consistent polycrystal plasticity formulation

This paper presents a multi-level crystal plasticity-based simulation framework for modeling mechanical response and microstructure evolution of Ti-6Al-4V with α-lath/lamellar microstructures. Here, the model is a correlated structure visco-plastic self-consistent (CS-VPSC) formulation linking three scales: a single crystals micro-scale, a lath/lamellar colony meso-scale, and a lath/lamellar aggregate macro-scale. A selected hardening law for the evolution of critical resolved shear stress per slip system used in CS-VPSC is phenomenological. However, it adjusts the resistances of basal and prismatic slip systems based on the geometry of slip transfer between adjacent lamellae. Consistent with experimental evidences, the resolved shear stress on the pyramidal slip planes is dependent not only on the stress in the direction of slip but also on the two orthogonal shear stress components and the three normal stress components (non-Schmid effects). Electron backscatter diffraction (EBSD) data in conjunction with a procedure relying on α→β phase transformation is used to construct paired variants of α-laths/lamellae satisfying their local crystallographic correlations. The procedure fits volume fractions of individual laths/lamellae with the experimental EBSD data and selects a distribution of habit planes between adjacent variants with respect to the loading direction. The simulation framework is applied to interpret the deformation behavior in tension and compression along two sample directions of Ti-6Al-4V fabricated via laser powder bed fusion. Moreover, the model is used to simulate texture evolution during rolling of the material to large strains. It is demonstrated that the model is capable of predicting plastic anisotropy/asymmetry and the concomitant texture evolution. While the model reveals a significant effect of habit plane inclination with respect to the loading direction on yield stress, the comparison of the data and model predictions shows that a random distribution of habit planes fits the flow response. It is further inferred that the tension-compression asymmetry arises from the non-Schmid effects.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Formation of buffer-gas-trap based positron beams

Presented here are experimental measurements, analytic expressions, and simulation results for pulsed, magnetically guided positron beams formed using a Penning-Malmberg style buffer gas trap. In the relevant limit, particle motion can be separated into motion along the magnetic field and gyro-motion in the plane perpendicular to the field. Analytic expressions are developed which describe the evolution of the beam energy distributions, both parallel and perpendicular to the magnetic field, as the beam propagates through regions of varying magnetic field. Simulations of the beam formation process are presented, with the parameters chosen to accurately replicate experimental conditions. The initial conditions and ejection parameters are varied systematically in both experiment and simulation, allowing the relevant processes involved in beam formation to be explored. These studies provide new insights into the underlying physics, including significant adiabatic cooling, due to the time-dependent beam-formation potential. Methods to improve the beam energy and temporal resolution are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Creep ductility limiting mechanisms in an additively manufactured Al-Ce-Ni-Mn-Zr alloy

Tensile creep response and cavitation damage evolution in an additively manufactured Al-7.5Ce-4.5Ni-0.4Mn-0.7Zr (wt%) alloy with peak-aging and overaging treatments were investigated in the 300–400 ºC range. Microstructural heterogeneity and its response to heat treatment and subsequent creep deformation were studied to understand the interplay between cavity formation, creep lifetime and ductility. Increasing the applied stress activated the nucleation of more cavities, an experimental observation that is well described using the vacancy accumulation model. Cavities nucleated prematurely due to localized plasticity in the denuded zones that formed at/near melt-pool or grain boundaries. Microstructure/deformation heterogeneity with consequent evolution of stress triaxiality, especially at lower stresses, causes accelerated cavitation, thus producing low creep ductility (∼ 0.2–2.4 %), compared to (∼12–21 %) ductility of the alloy measured by regular tensile tests at equivalent temperatures. A constrained diffusional cavity growth mechanism with continuous cavity nucleation during creep is established as the dominant mechanism, implying that cavitation involves vacancy diffusion, yet its growth rate is dictated by the minimum creep rate. In conclusion, the ductility-limiting creep and cavitation mechanisms discussed here provide new insight into the creep behavior of 3D-printed metallic alloys.

Additive manufacturing↗
DOE OSTI2019

Reversible and Selective CO 2 to HCO 2 - Electrocatalysis near the Thermodynamic Potential

Reversible catalysis is an emblem of energy-efficient chemical transformations, but can only be achieved if the changes in free energy of intermediate steps are minimized and the catalytic cycle is devoid of high transition-state barriers. Using these criteria, we show reversible CO 2 /HCO 2 - conversion catalyzed by [Pt(depe) 2 ] 2+ (depe=1,2- bis (diethylphosphino)ethane). Direct measurement of the free energies associated with each catalytic step correctly predicts a slight bias towards CO 2 reduction. We demonstrate how the experimentally measured free energy of each step directly contributes to the <50mV overpotential. We also find that for CO 2 reduction, H 2 evolution is negligible and the Faradaic efficiency for HCO 2 - production is nearly quantitative. A free-energy analysis reveals H 2 evolution is endergonic, providing a thermodynamic basis for highly selective CO 2 reduction.

10 SYNTHETIC FUELS↗
DOE OSTI2026

Predicting fracture behavior in single crystal nickel using a coupled crystal plasticity phase field damage approach

Understanding the fracture behavior of single crystal metals is critical for predicting material performance under mechanical loading. Here, in this study, we investigate the fracture characteristics of single crystal nickel tensile bars using a crystal plasticity coupled phase field damage (CP-PFD) model. Experimental tensile tests were conducted on 15 specimens spanning five crystallographic orientations and three thickness variants per orientation. The results revealed two distinct fracture modes: brittle fractures with 45-degree angled surfaces and ductile fractures characterized by significant necking. The CP-PFD model successfully replicated these fracture behaviors, demonstrating strong agreement with experimental observations. The model effectively predicted the strain at which necking and fracture occurred, as well as the orientation-dependent fracture mechanisms. By comparing experimental and simulated fracture surfaces, we establish the CP-PFD model as a robust tool for predicting single crystal behavior and damage evolution. This work provides insight into the microstructural dependence of fracture behavior and establishes a predictive framework for modeling orientation-dependent damage evolution in single-crystal nickel.

Crystal plasticity↗
DOE OSTI2022

Explosion behavior investigation and safety assessment of large-format lithium-ion pouch cells

We report large-format lithium-ion (Li-ion) batteries with high energy density for electric vehicles are prone to thermal runaway (or even explosion) under abusive conditions. In this study, overcharge induced explosion behaviors of large-format Li-ion pouch cells with Li[Ni 0.8 Co 0.1 Mn 0.1 ]O 2 cathode at different current rates (C-rates) (0.5C, 1C, 2C) were investigated. The explosion characteristics of the cells were elucidated by discussing the evolution of the cell voltage, the surface temperature and the shock wave pressure. Generally, the whole overcharge process could be divided into four stages according to the evolution of several key parameters and the overcharge behaviors; the overcharge C-rate has a great influence on cells’ thermal behaviors. The experimental results showed that the thermal runaway process of Li-ion cells caused by overcharging consisted of two kinds of explosions, physical explosion and chemical explosion. The existence of observable negative pressure zone in the pressure curves indicated that the Li-ion cells are not a self-supplying oxygen system during the explosion. Further, the explosion dynamics parameters were matched. An explosion TNT-equivalent conversion strategy that depended on the pressure of the shock wave was utilized to evaluate the released energy and its hazards. In addition, with respect to the overcharge of Li-ion pouch cells, a safety assessment method and a safety management method were proposed based on the explosion behaviors. From the perspective of battery safety, this study is of great significance for the safety design of Li-ion cells and can provide guidance for engineers to optimize the safety function of battery packs.

25 ENERGY STORAGE↗
DOE OSTI2018

On visualizing material flow and precipitate evolution during probeless friction stir spot welding of an Al-Li alloy

Highlights: • Severe thermal cycle and deformation refined grains and intermetallic particles in SZ. • T{sub 1} phases fully dissolved in SZ followed by re-precipitation of GP zone and δ′ phases. • The predominant flow direction remained in spiral-like pattern during P-FSSW. • Local abrasion existed in the lap-weld interface due to flow difference, forming joining. - Abstract: The precipitate evolution and material flow visualization during probeless friction stir spot welding (P-FSSW) of an Al-Li alloy has been investigated experimentally and numerically in this work. Due to severe heat input during welding, the main strengthening phases of T{sub 1} fully dissolved in the stir zone followed by the re-precipitation of GP zone and δ′ phases, assisted with the refinement of the intermetallic particles and grains. The history of material flow was three-dimensionally visualized with a tracer, studied with X-ray tomography. The experiments and the simulation results provided a consistent interpretation of a circumferential material flow which moved inwards radially near the top surface. Upward flow in the bottom sheet improved with dwell time, causing the interface to bend. Local abrasion between the two sheets due to flow difference led to the disruption of oxide layer, and formed a connection at the interface.

36 MATERIALS SCIENCE↗
DOE OSTI2018

A brittle fracture mechanism in thermally aged duplex stainless steels revealed by in situ high-energy X-ray diffraction

A direct relationship between the microstructural evolution and macroscopic fracture behaviors of a thermally aged (at 475 °C for 400 h) duplex stainless steel (DSS) has been established with experimental results from atom probe tomography (APT), nanoindentation and in situ synchrotron-based high-energy X-ray diffraction (HE-XRD). Furthermore, the APT experiments demonstrate that ferrite in DSS spinodally decomposes into Cr-enriched and Cr-depleted domains during thermal aging, which leads to a severe hardening effect in ferrite. The lattice strain development during deformation acquired with the in situ HE-XRD measurements confirms the cleavage fracture of α{110} and α{200} ferrite grains aligned perpendicular to loading direction. The thermally aged DSS can be readily fractured by connecting the cleavage cracks in ferrite. Conclusively, the spinodal-decomposition-induced hardening in ferrite and the premature failure of ferrite control the final brittle fracture in the aged DSS.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Investigation of Structural Evolution of Li 1.1 V 3 O 8 by In Situ X-ray Diffraction and Density Functional Theory Calculations

Li 1.1 V 3 O 8 is a promising cathode material for Li-ion batteries, because of its high theoretical capacity (362 mAh g-1) and good rate performance. In this study, the structural evolution of Li 1.1 V 3 O 8 material during electrochemical dis(charge) processes was investigated using a combination of theoretical calculations and experimental data. Density functional theory (DFT) was used to predict the intermediate structures at various lithiation states, as well as the stability of major phases. In order to validate these predictions, in situ X-ray diffraction (XRD) data was collected operando, allowing for the phase transformations to be monitored under current load and eliminating the possibilities of structural relaxation processes and environmental oxidation. Rietveld refinement was performed to fit the diffraction data with the DFT-derived structures and to analyze the fractions of major phases as a function of dis(charge). The DFT calculations identified three stable states that were validated by the in situ XRD result: a Li-poor α-phase (Li1), a Li-rich α-phase (Li2.5), and a β-phase (Li4). The DFT-predicted particle shape based on the surface energy of the (100), (001), and (010) planes rationalized the preferential orientation of Li 1.1 V 3 O 8 particles along the [010] direction in the electrode. Furthermore, the onset and offset of the α → β transition, as well as the phase fractions of α and β determined via in situ XRD, related well with the DFT-derived relative stability of each phase. Thus, by integrating DFT calculations with experimental work, this work provides a thorough understanding of the structural transformations in Li 1.1 V 3 O 8 during electrochemical dis(charge).

25 ENERGY STORAGE↗
DOE OSTI2018

Physical Origins of the Transient Absorption Spectra and Dynamics in Thin-Film Semiconductors: The Case of BiVO 4

Transient absorption (TA) spectroscopy is uniquely suited for understanding kinetic processes initiated by light over vast ranges of time. In combination with white light probes, the recorded differential absorption spectra can contain spectroscopic signatures characteristic of specific charge carrier population densities. However, disentangling the often-complex and convoluted spectra is made challenging without robust analysis methods relating the underlying physical mechanisms to the spectral components. In this work, we address the origin of the transient spectra of a model system of emerging solar energy harvesting materials using a monoclinic BiVO 4 thin film. Using ground-state optical properties of the semiconductor, we find the main derivative-like spectral response to be related to shifting and broadening of oscillators, rather than specific carrier-related transitions. However, by using the Drude optical model of free carriers, we also identify the transient response related to free-hole density. Importantly, sample heating from the optical pumping, which begins at ~10 ps and plateaus by ~200 ps, dominates the overall spectral response at longer times. On the basis of a physical model of the spectral response, a kinetic model is developed that describes the pump power dependence of the free-hole density, as well as the temporal evolution of the spectral changes associated with shifting and broadening of oscillators. First-principles density functional theory calculations are used to rationalize experimental measurements. Here, this comprehensive approach to analyzing and modeling the TA spectra offers a generalizable basis for understanding the complex pump-probe data, reveals thermal heating artifacts that are frequently erroneously assigned to long-lived photocarriers, and offers a path to eliminating ambiguity in analysis of photocarrier dynamics in solid-state systems.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Quantum Mechanical Screening of Single-Atom Bimetallic Alloys for the Selective Reduction of CO 2 to C 1 Hydrocarbons

Electrocatalytic reduction of CO 2 to energy-rich hydrocarbons such as alkanes, alkenes, and alcohols is a very challenging task. So far, only copper has proven to be capable of such a conversion. Here, we report density functional theory (DFT) calculations combined with the Poisson-Boltzmann implicit solvation model to show that single-atom alloys (SAAs) are promising electrocatalysts for CO 2 reduction to C 1 hydrocarbons in aqueous solution. The majority component of the SAAs studied is either gold or silver, in combination with isolated single atoms, M (M = Cu, Ni, Pd, Pt, Co, Rh, and Ir), replacing surface atoms. We envision that the SAA behaves as a one-pot tandem catalyst: first gold (or silver) reduces CO 2 to CO, and the newly formed CO is then captured by M and is further reduced to C 1 hydrocarbons such as methane or methanol. We studied 28 SAAs, and found about half of them selectively favor the CO 2 reduction reaction over the competing hydrogen evolution reaction. Most of those promising SAAs contain M = Co, Rh, or Ir. The reaction mechanism of two SAAs, Rh@Au(100) and Rh@Ag(100), is explored in detail. Both of them reduce CO 2 to methane but via different pathways. For Rh@Au(100), reduction occurs through the pathway *CO → *CHO → *CHOH → *CH + H 2 O (l) → *CH 2 + H 2 O (l) → *CH 3 + H 2 O (l) → * + H 2 O (l) + CH 4(g) ; whereas, for Rh@Ag(100), the pathway is *CO → *CHO → *CH 2 O→ *OCH 3 → *O + CH 4(g) → *OH + CH 4(g) →* + H 2 O (l) + CH 4(g) . The minimum applied voltages to drive the two electrocatalytic systems are -1.01 and -1.12 V RHE for Rh@Au(100) and Rh@Ag(100), respectively, at which the Faradaic efficiencies for CO 2 reduction to CO are 60% for gold and 90% for silver. This suggests that SAA can efficiently reduce CO 2 to methane with as small as 40% loss to the hydrogen evolution reaction for Rh@Au(100) and as small as 10% for Rh@Ag(100). Lastly, we hope these computational results can stimulate experimental efforts to explore the use of SAA to catalyze CO 2 electrochemical reduction to hydrocarbons.

36 MATERIALS SCIENCE↗