Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “oxygen defects”

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 469 records · Page 26

Ammonolysis with N 2 -diluted NH 3 suppresses Ta( IV ) defects in BaTaO 2 N and enhances photocatalytic water oxidation

BaTaO 2 N stands out among oxynitride photocatalysts because of its ability to capture visible light and to drive the photoelectrochemical water oxidation reaction. However, its solar energy conversion performance is limited by electron–hole recombination at Ta( IV ) defects in the material. These defects are formed by overreduction of the Ta(v) oxide precursor by excess ammonia under the high temperature conditions during ammonolysis. Here we show for the first time that Ta( IV ) defect concentrations can be lowered by conducting the ammonolysis reaction in mixed NH 3 /N 2 gas. The obtained BaTaO 2 N samples crystallize in the cubic CaTiO 3 structure type and form 200–300 nm faceted nanocrystals, based on X-ray diffraction, scanning electron microscopy, and HRTEM. Electron paramagnetic resonance spectra observe the Ta( IV ) defects at g = 1.999 and confirm an 11-fold reduction for the product synthesized in mixed (0.13 : 1.0 vol) NH 3 /N 2 gas, equivalent to 1.14 × 10 16 cm −3 Ta( IV ) ions. This optimized BaTaO 2 N has nearly twice the photocatalytic oxygen evolution activity (AQE of 6.78% at 400 nm) of a reference material made with 1.0 atm ammonia and 78% higher photoelectrochemical water oxidation photocurrent (0.9 mA cm −2 at 1.23 V vs. RHE) under simulated sunlight. According to X-ray photoelectron spectroscopy, remaining Ta( IV ) defects are concentrated in the surface region of the BaTaO 2 N particles, where >50% of all Ta ions are found in the +4 oxidation state. This surface Ta( IV ) population can be directly observed in Vibrating Kelvin Probe Surface Photovoltage Spectra (VK-SPV) via its 1.2–1.4 eV photovoltage onset. Here, it suggests that the surface Ta( IV ) ions contribute empty d-states 0.5–0.7 eV below the BaTaO 2 N conduction band edge. These findings highlight how the energetics and concentrations of Ta( IV ) defects influence the photoelectrochemical water oxidation ability of BaTaO 2 N. Additionally, the work establishes ammonolysis with diluted NH 3 as a new tool to minimize defects in BaTaO 2 N and to raise its solar energy conversion efficiency toward its theoretical limit. Because of its simplicity, the reduced ammonia pressure strategy will likely be applicable to other oxynitrides, which generally suffer from overreduction problems during ammonolysis.

Salmanion, Mahya [University of California, Davis,↗

Ammonolysis Under NH 3 –Limiting Conditions as a Pathway to Improved LaTiO 2 N Water Splitting Photoanodes

LaTiO 2 N is a promising intermediate band gap semiconductor for the water splitting reaction, a pathway to hydrogen fuel from solar energy. However, the photoelectrochemical (PEC) activity of the material is hindered by defects, particularly Ti(III) species, which promote photocarrier recombination. These defects are formed during the high-temperature ammonolysis reaction. Here we show that improved LaTiO 2 N materials can be synthesized under NH 3 -limiting conditions by introducing N 2 to lower the NH 3 partial pressure to0.13atm.This reduces the Ti(III) defect density in the material from 6.06 × 10 16 to ∼4.61 × 10 15 cm −3 , by a factor of 13, based on electron paramagnetic resonance (EPR) spectroscopy. Any remaining Ti(III) defects are localized at the LaTiO 2 N surface, according to X-ray photoelectron spectroscopy (XPS), due to the formation of a depletion layer in the semico. Optical absorption spectra of the improved LaTiO 2 N reveal a blue-shifted band gap absorption edge and a suppressed sub-band gap absorption. Defect removal also reduces a sub-band gap surface photovoltage feature visible in the 1.0 atm reference material. The improved LaTiO 2 N supports a 1.57 mA cm −2 water oxidation photocurrent at 1.23 V RHE under simulated sunlight conditions, and an enhanced quantum efficiency of 4.5% (400 nm) for photocatalytic oxygen evolution from aqueous silver nitrate solution. Stable PEC operation is observed for over 55 min. This confirms that ammonolysis under NH 3 -limiting conditions improves the solar energy conversion properties of LaTiO 2 N. The ability to control metal ion defects in oxynitrides by varying the ammonia partial pressure during ammonolysis might be generally useful for the preparation of metal nitrides and oxynitrides.

defects↗

Computational Fermi Level Engineering and Doping-Type Conversion of Ga2O3 Via Three-Step Processing

Ga2O3 is being actively explored for power electronics, and other applications due to its ultra- wide bandgap and low projected fabrication cost of high-quality crystals. N-type doping of Ga2O3 can be achieved and tuned, but p-type doping faces fundamental obstacles. Successful engineering of Ga2O3 based devices requires critical control of doping density, Fermi level position, providing opportunities for predictive process simulation. We use first-principles defect equilibrium calculations to simulate a 3-step growth- annealing-quench protocol for hydrogen assisted Mg doping in Ga2O3, taking into account the hydrogen- oxygen-water gas phase equilibrium. We predict type conversion to a net p-type regime following O-rich annealing after growth under reducing conditions in the presence of H2. We show that there is an optimal temperature that maximizes the net acceptor density during the equilibrium annealing step for a given Mg doping level. Quenching of non-equilibrium annealed samples then results in a Fermi level EF below mid-gap down to about EV +1.5 eV, creating a significant number of uncompensated neutral MgGa0 acceptors. This type converted Ga2O3 can create significant built-in field in a p-n junction with an adjoining n-type material and could enable vertical MOSFETs.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS,M↗

First-principles thermodynamic assessment of Sr-containing secondary phase formation in strontium-substituted lanthanum manganites for solid oxide cell applications

Sr-secondary phase formation is a potentially significant degradation mode with direct impact upon solid-oxide cell (SOC) commercial viability. A first-principles based thermodynamic study was performed for La 1−x Sr x MnO 3±δ (LSM) perovskites to assess their stability against formation of different Sr-secondary phases, including SrO, SrCrO 4 , SrSO 4 , SrCO 3 , and Sr(OH) 2 , for SOC applications. The Sr-secondary phase formation reaction free energies were determined via a thermodynamic model by combining ab initio lattice dynamics calculations for the solid phases and ab initio thermodynamic data for the gas phases. The current approach expands the previously reported thermodynamic modeling studies by integrating first-principles based point defect equilibria into the thermodynamic analysis. The modeling results obtained using this new approach indicate an increased tendency to form the SrO oxide upon decreasing the oxygen partial pressure. Additionally, the enhancing factors to form the Sr-related secondary phase from the associated SrO activity in LSM are further quantified by considering the equilibrium of SrO reacting with the contaminant gas species as a function of temperature and pressure.

defect thermodynamics modeling↗

Deep and Shallow Gap States in Reduced and n-Type Doped m -ZrO 2

Monoclinic zirconium dioxide (m-ZrO 2 ) is a wide-band-gap functional oxide with many applications. Understanding the role of intrinsic defects and dopants is important for improving the properties of m-ZrO 2 relevant to its applications. In this work, we characterize the electronic states and energy levels of oxygen vacancies (VO’s) and n-type Nb, Ta, Rb, H, and F dopants in the bulk and at the majority (1̅11) surface of m-ZrO 2 , using accurate dielectric-dependent hybrid functional calculations. Our results show that VO’s generally behave as deep donors with the excess electrons localized at the vacant O sites in the form of F centers, in agreement with previous studies. In contrast, surface VO’s at two-fold coordinated oxygen sites show a rather shallow (1+/0) transition level. Shallow levels are also obtained for adsorbed hydrogen (H ads ) and fluorine substituting an oxygen atom (FO) either in the bulk or at the surface, whereas Nb and Ta form deep levels. Altogether, our results suggest that FO is the best candidate for realizing n-type conductivity in bulk m-ZrO 2 , while surface VO’s and Hads provide an efficient way to control the surface chemistry of m-ZrO 2 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Defect-Driven Oxide Transformations and the Electrochemical Interphase

The redox reaction pathway is crucial to the sustainable production of the fuels and chemicals required for a carbon-neutral society. Our society is becoming increasingly dependent on devices using batteries and electrolyzers, all of which rely on a series of redox reactions. The overall properties of oxide materials make them very well suited for such electrochemical and catalytic applications due to their associated cationic redox properties and the static site-adsorbate interactions. As these technologies have matured, it has become apparent that defect-driven redox reactions, defect-coupled diffusion, and structural transformations that are both time- and rate-dependent are also critical materials processes. This change in focus, considering not only redox properties but also more complex, dynamic behaviors, represents a new research frontier in the molecular sciences as they are strongly linked to device operation and degradation and lie at the heart of various phenomena that take place at electrochemical interfaces. Fundamental studies of the structural, electronic, and chemical transformation mechanisms are key to the advancement of materials and technological innovations that could be implemented in various electrochemical systems. In this Account, we focus on recent studies and advances in characterizing and understanding the dynamic redox evolution and structural transformations that take place in model perovskites and layered oxides under reactive conditions and correlate them with degradation mechanisms and operations in electrolyzers and batteries. We show that the dynamic evolution of oxygen vacancies and cationic migration in the surface or bulk occurs at the solid-liquid interface, using a combination of different synchrotron-based Xray spectroscopies and scattering probes. Detailed redox-structure-reactivity correlation studies show how defects and diffusion processes can be tailored to drive various physical and chemical transformations in electrolyzers and batteries. We also highlight a strong correlation between oxygen redox reactivity and structural reorganization in both model thin films and particles, helping to bridge the gap between fundamental studies of the reaction mechanism and device applications. Finally, on the basis of these findings, we discuss strategies to probe and tune the redox reactivity and structural stability of the redox-active oxide interphase toward devising efficient pathways for energy and chemical harvesting.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Irradiation Driven Restructuring of Nanocrystalline ThO 2 and Th 1–x U x O 2 Thin Films

Irradiation induced structural changes of actinide oxide materials is a key consideration in their development and use as nuclear fuels. This study reported on the synthesis of ThO 2 and Th 1–x U x O 2 (x = 0.15, 0.50) thin films, fabricated using electrospray-assisted solution combustion synthesis, and their responses to ion irradiation. Krypton ion irradiations, up to a fluence of 1 × 10 16 ions/cm 2 , were carried out to simulate radiation damage induced by fission products in a reactor environment. Structural and chemical changes induced by irradiation were analyzed using high-resolution scanning transmission electron microscopy (STEM), energy-dispersive X-ray spectroscopy (EDS), and electron energy-loss spectroscopy (EELS). It was determined that the extent and nature of irradiation-induced damage are strongly correlated with the uranium content. ThO 2 films were most susceptible to radiation-induced damage, with significant cavity formation and delamination from the substrate at high fluence. Of the compositions studied, Th 0.85 U 0.15 O 2 films showed the highest stability, characterized by moderate grain growth and the absence of voids or severe defect structures. In contrast, Th 0.5 U 0.5 O 2 films accumulated extensive damage, including the formation of a nanocrystalline central region. EELS analysis indicated that oxygen displacement is the primary driver of structural degradation in Th 0.5 U 0.5 O 2 films. α-particle spectroscopy confirmed minimal actinide loss across all compositions, underscoring the mechanical robustness of the films. These findings provide insight into the irradiation-induced damage mechanisms in Th O2 and Th 1–x U x O 2 systems, supporting their development as potential materials for nuclear fuels and irradiation-tolerant thin film targets in nuclear physics measurements.

Th1−xUxO2↗

Solar Wind Implantation Into the Lunar Regolith: Monte Carlo Simulations of H Retention in a Surface With Defects and the H2 Exosphere

The solar wind implants protons into the top 20–30 nm of lunar regolith grains, and the implanted hydrogen will diffuse out of the regolith but also interact with oxygen in the regolith oxides. We apply a statistical approach to estimate the diffusion of hydrogen in the regolith hindered by forming temporary bonds with regolith oxygen atoms. A Monte Carlo simulation was used to track the temporal evolution of bound OH surface content and the H2 exosphere. The model results are consistent with the interpretation of the Chandrayaan‐1 M3 observations of infrared absorption spectra by surface hydroxyls as discussed in Li and Milliken (2017, https://doi.org/10.1126/sciadv.1701471). The model reproduced the latitudinal concentration of OH by using a Gaussian energy distribution of f(U(sub 0) = 0.5 eV, U(sub W) = 0.078–0.1 eV) to characterize the activation energy barrier to the diffusion of hydrogen in space weathered regolith. In addition, the model results of the exospheric content of H2 are consistent with observations by the Lyman Alpha Mapping Project on the Lunar Reconnaissance Orbiter. Therefore, we provide support for hydroxyl formation by chemically trapped solar wind protons.

Tucker, O. J.↗

Solar Wind Implantation Into the Lunar Regolith: Monte Carlo Simulations of H Retention in a Surface With Defects and the H2 Exosphere

The solar wind implants protons into the top 20–30 nm of lunar regolith grains, and the implanted hydrogen will diffuse out of the regolith but also interact with oxygen in the regolith oxides. We apply a statistical approach to estimate the diffusion of hydrogen in the regolith hindered by forming temporary bonds with regolith oxygen atoms. A Monte Carlo simulation was used to track the temporal evolution of bound OH surface content and the H2 exosphere. The model results are consistent with the interpretation of the Chandrayaan-1 M3 observations of infrared absorption spectra by surface hydroxyls as discussed in Li and Milliken (2017). Themodel reproduced the latitudinal concentration of OH by using a Gaussian energy distribution of f(U(o) = 0.5eV, U(w) = 0.078 – 0.1 eV) to characterize the activation energy barrier to the diffusion of hydrogen in space weathered regolith. In addition, the model results of the exospheric content of H2 are consistent with observations by the Lyman Alpha Mapping Project on the Lunar Reconnaissance Orbiter. Therefore, we provide support for hydroxyl formation by chemically trapped solar wind protons.

O. J. Tucker↗

Solar Wind Implantation Into the Lunar Regolith: Monte Carlo Simulations of H Retention in a Surface With Defects and the H2 Exosphere

The solar wind implants protons into the top 20–30 nm of lunar regolith grains, and the implanted hydrogen will diffuse out of the regolith but also interact with oxygen in the regolith oxides. We apply a statistical approach to estimate the diffusion of hydrogen in the regolith hindered by forming temporary bonds with regolith oxygen atoms. A Monte Carlo simulation was used to track the temporal evolution of bound OH surface content and the H2 exosphere. The model results are consistent with the interpretation of the Chandrayaan‐1 M3 observations of infrared absorption spectra by surface hydroxyls as discussed in Li and Milliken (2017, https://doi.org/10.1126/sciadv.1701471). The model reproduced the latitudinal concentration of OH by using a Gaussian energy distribution of f(U0 = 0.5 eV, UW = 0.078–0.1 eV) to characterize the activation energy barrier to the diffusion of hydrogen in space weathered regolith. In addition, the model results of the exospheric content of H2 are consistent with observations by the Lyman Alpha Mapping Project on the Lunar Reconnaissance Orbiter. Therefore, we provide support for hydroxyl formation by chemically trapped solar wind protons.

O J Tucker↗

Oxygen Vacancies Alter Methanol Oxidation Pathways on NiOOH

A thorough comprehension of the mechanism underlying the methanol oxidation reaction (MOR) on Ni-based catalysts is critical for future electrocatalytic design and development. However, the mechanism of MOR on these materials remains a matter of controversy. Herein, we combine in situ surface-enhanced infrared absorption spectroscopy (SEIRAS) and density functional theory (DFT) calculations to identify the active sites and determine the mechanism of MOR on monometallic Ni-based catalysts in alkaline media. The SEIRAS results show that formate and (bi)carbonate are formed after the commencement of the MOR with potential-dependent relative distributions. These spectroscopic results are in good agreement with the DFT-computed reaction profiles over an oxygen vacancy, suggesting that the MOR mainly proceeds through the formate-involving pathway, in which the early consumption of methanol yields formate as the major product, while increasing potential drives further oxidation of formate to (bi)carbonate. We also find a parallel pathway for the generation of (bi)carbonate at high potentials that bypasses the formation of formate. The two main pathways are thermodynamically more feasible than the one predominantly reported in the literature for MOR on NiOOH that involves CHO and/or CO as key intermediates. These DFT results are supported by spectroscopic evidence showing that no band associated with CHO or CO can be detected by SEIRAS, which is attributed to the nature of the oxygen vacancies as the active sites, suppressing deep dehydrogenation of CH 2 O to CHO. Furthermore, this work thus shows the promising role of defect engineering in promoting the electrocatalytic MOR activity and selectivity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Microscopic relaxation channels in materials for superconducting qubits

Despite mounting evidence that materials imperfections are a major obstacle to practical applications of superconducting qubits, connections between microscopic material properties and qubit coherence are poorly understood. Here, we combine measurements of transmon qubit relaxation times ( T 1 ) with spectroscopy and microscopy of the polycrystalline niobium films used in qubit fabrication. By comparing films deposited using three different techniques, we reveal correlations between T 1 and intrinsic film properties such as grain size, enhanced oxygen diffusion along grain boundaries, and the concentration of suboxides near the surface. Qubit and resonator measurements show signatures of two-level system defects, which we propose to be hosted in the grain boundaries and surface oxides. We also show that the residual resistance ratio of the polycrystalline niobium films can be used as a figure of merit for qubit lifetime. This comprehensive approach to understanding qubit decoherence charts a pathway for materials-driven improvements of superconducting qubit performance.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Fluorine spillover for ceria- vs silica-supported palladium nanoparticles: A MD study using machine learning potentials

Supported metallic nanoparticles play a central role in catalysis. However, predictive modeling is particularly challenging due to the structural and dynamic complexity of the nanoparticle and its interface with the support, given that the sizes of interest are often well beyond those accessible via traditional ab initio methods. With recent advances in machine learning, it is now feasible to perform MD simulations with potentials retaining near-density-functional theory (DFT) accuracy, which can elucidate the growth and relaxation of supported metal nanoparticles, as well as reactions on those catalysts, at temperatures and time scales approaching those relevant to experiments. Furthermore, the surfaces of the support materials can also be modeled realistically through simulated annealing to include effects such as defects and amorphous structures. We study the adsorption of fluorine atoms on ceria and silica supported palladium nanoparticles using machine learning potential trained by DFT data using the DeePMD framework. We show defects on ceria and Pd/ceria interfaces are crucial for the initial adsorption of fluorine, while the interplay between Pd and ceria and the reverse oxygen migration from ceria to Pd control spillover of fluorine from Pd to ceria at later stages. In contrast, silica supports do not induce fluorine spillover from Pd particles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Two-photon excitation cross-section in light and intermediate atoms

The method of explicit summation over the intermediate states is used along with LS coupling to derive an expression for two-photon absorption cross section in light and intermediate atoms in terms of integrals over radial wave functions. Two selection rules, one exact and one approximate, are also derived. In evaluating the radial integrals, for low-lying levels, the Hartree-Fock wave functions, and for high-lying levels, hydrogenic wave functions obtained by the quantum defect method are used. A relationship between the cross section and the oscillator strengths is derived. Cross sections due to selected transitions in nitrogen, oxygen, and chlorine are given. The expression for the cross section is useful in calculating the two-photon absorption in light and intermediate atoms.

Omidvar, K.↗

Rhenium-Oxygen Interactions at High Temperatures

The reaction of pure rhenium metal with dilute oxygen/argon mixtures was studied from 600 to 1400 C. Temperature, oxygen pressure, and flow rates were systematically varied to determine the rate-controlling steps. At lower temperatures the oxygen/rhenium chemical reaction is rate limiting; at higher temperatures gas-phase diffusion of oxygen through the static boundary layer is rate limiting. At all temperatures post-reaction microstructures indicate preferential attack along certain crystallographic planes and defects.

Jacobson, Nathan S.↗

Simulation of etch pit formation through active sites in carbon fiber micro-structures

Erosion of carbon due to oxidation does not occur uniformly but through localized etch pit formation as a result of active surface sites. In this work we demonstrate a newly developed capability to capture pitting of carbon fiber microstructures such as FiberForm, which is commonly used as the base material for NASA’s spacecraft ablative thermal protection systems (TPS). The simulations are performed at the meso-scale in order to capture the pit formation and growth using direct simulation Monte Carlo (DSMC), specifically using the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code. Legacy and latest models both assume uniform reactivity of carbon surface sites with oxygen even at the meso-scale level. However, in reality the carbon surface has locally different reaction rates due to the presence of defects at the atomic scale. The defective nature of these sites enhances their reactivity with atmospheric gases compared to the non-defective sites (2-3 orders of magnitude) and are termed as “active sites”. Thus, these sites tend to be the first to react and eventually get removed through the formation of gasses such as CO, CO2, and CN. Their removal results in all the neighboring atoms becoming defective, thus leading to chain reaction of localized carbon removal and formation of etch pits. Capturing the formation of pits during the ablation simulation of carbon micro-structures is critical to predicting their structural failure. Recently a detailed surface chemistry framework was developed in SPARTA, capable of incorporating various reaction mechanisms such as adsorption, desorption, Eley-Rideal (ER) and Langmuir-Hinshelwood (LH) mechanisms. We have implemented the capability of a single surface having multiple site sets with different reactivities within this framework. We have used this feature to model the presence of active sites on carbon surfaces, whose reactivity is orders of magnitude higher than that of the passive sites due to the presence of defects. The active site fraction is a property of surface elements within SPARTA and is directly proportional to the local reactivity of each surface element. By introducing an initial distribution of the active site fraction across the carbon surface and propagating it in a manner that mimics the evolution of real reacting carbon surfaces, we can capture the formation and growth of etch pits as a result of surface consumption reactions such as oxidation.

DSMC↗

Simulation of Etch Pit Formation in DSMC Through Active Sites in Carbon Fiber Micro-Structures

Erosion of carbon due to oxidation does not occur uniformly but through localized etch pit formation because of active surface sites. In this work we demonstrate a newly developed capability to capture pitting of carbon fiber microstructures such as FiberForm, which is commonly used as the base material for NASA’s spacecraft ablative thermal protection systems (TPS). The simulations are performed at the meso-scale in order to capture the pit formation and growth using direct simulation Monte Carlo (DSMC), specifically using the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code. Legacy and latest models both assume uniform reactivity of carbon surface sites with oxygen even at the meso-scale level. However, in reality the carbon surface has locally different reaction rates due to the presence of defects at the atomic scale. The defective nature of these sites enhances their reactivity with atmospheric gases compared to the non-defective sites (2-3 orders of magnitude) and are termed as “active sites”. Thus, these sites tend to be the first to react and eventually get removed through the formation of gases such as CO, CO2, and CN. Their removal results in all the neighboring atoms becoming defective, thus leading to chain reaction of localized carbon removal and formation of etch pits. Capturing the formation of pits during the ablation simulation of carbon micro-structures is critical to predicting their structural failure. Recently a detailed surface chemistry framework was developed in SPARTA, capable of incorporating various detailed surface reaction mechanisms. We have implemented the capability of a single surface having multiple site sets with different reactivities within this framework. We have used this feature to model the presence of active sites on carbon surfaces, whose reactivity is orders of magnitude higher than that of the passive sites due to the presence of defects. The active site fraction is a property of surface elements within SPARTA and is directly proportional to the local reactivity of each surface element. By introducing an initial distribution of the active site fraction across the carbon surface and propagating it in a manner that mimics the evolution of real reacting carbon surfaces, we can capture the formation and growth of etch pits as a result of surface consumption reactions such as oxidation.

DSMC↗

DSMC Simulation of Etch Pit Formation Through Active Sites in Carbon Fiber Micro-Structures

Erosion of carbon due to oxidation does not occur uniformly but through localized etch pit formation because of active surface sites. In this work we demonstrate a newly developed capability to capture pitting of carbon fiber microstructures such as FiberForm, which is commonly used as the base material for NASA’s spacecraft ablative thermal protection systems (TPS). The simulations are performed at the meso-scale in order to capture the pit formation and growth using direct simulation Monte Carlo (DSMC), specifically using the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code. Legacy and latest models both assume uniform reactivity of carbon surface sites with oxygen even at the meso-scale level. However, in reality the carbon surface has locally different reaction rates due to the presence of defects at the atomic scale. The defective nature of these sites enhances their reactivity with atmospheric gases compared to the non-defective sites (2-3 orders of magnitude) and are termed as “active sites”. Thus, these sites tend to be the first to react and eventually get removed through the formation of gases such as CO, CO2, and CN. Their removal results in all the neighboring atoms becoming defective, thus leading to chain reaction of localized carbon removal and formation of etch pits. Capturing the formation of pits during the ablation simulation of carbon micro-structures is critical to predicting their structural failure. Recently a detailed surface chemistry framework was developed in SPARTA, capable of incorporating various detailed surface reaction mechanisms. We have implemented the capability of a single surface having multiple site sets with different reactivities within this framework. We have used this feature to model the presence of active sites on carbon surfaces, whose reactivity is orders of magnitude higher than that of the passive sites due to the presence of defects. The active site fraction is a property of surface elements within SPARTA and is directly proportional to the local reactivity of each surface element. By introducing an initial distribution of the active site fraction across the carbon surface and propagating it in a manner that mimics the evolution of real reacting carbon surfaces, we can capture the formation and growth of etch pits as a result of surface consumption reactions such as oxidation.

DSMC↗