Search NASA⌕ Search

SEARCH · Search NASA

Results for “atomic chains”

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 145 records · Page 8

Observation of Individual Fluorine Atom from Highly Oriented Poly (tetrafluoroethylene) Films by Atomic Force Microscopy

Direct observation of the film thickness, molecular structure and individual fluorine atoms from highly oriented poly(tetrafluoroethylene) (PTFE) films were achieved using atomic force microscopy (AFM). A thin PTFE film is mechanically deposited onto a smooth glass substrate at specific temperatures by a friction transfer technique. Atomic resolution images of these films show that the chain-like helical structures of the PTFE macromolecules are aligned parallel to each other with an intermolecular spacing of 5.72 A, and individual fluorine atoms are clearly observed along these twisted molecular chains with an interatomic spacing of 2.75 A. Furthermore, the first direct AFM measurements for the radius of the fluorine-helix, and of the carbon-helix in sub-angstrom scale are reported as 1.70 A and 0.54 A respectively.

Lee, Jonathan A.,↗

Observation of Individual Fluorine Atoms from Highly Oriented Poly(Tetrafluoroethylene) Films by Atomic Force Microscopy

Direct observation of the film thickness, molecular structure, and individual fluorine atoms from highly oriented poly(tetrafluoroethylene) (PTFE) films were achieved using atomic force microscopy (AFM). A thin PTFE film is mechanically deposited onto a smooth glass substrate at specific temperatures by a friction-transfer technique. Atomic resolution images of these films show that the chain-like helical structures of the PTFE macromolecules are aligned parallel to each other with an intermolecular spacing of 5.72 A, and individual fluorine atoms are clearly observed along these twisted molecular chains with an interatomic spacing of 2.75 A. Furthermore, the first direct AFM measurements for the radius of the fluorine-helix, and of the carbon-helix in sub-angstrom scale are reported as 1.7 and 0.54 A respectively.

Lee, J. A.↗

La 4 Co 4 X ( X = Pb , Bi , Sb ) : A demonstration of antagonistic pairs as a route to quasi-low-dimensional ternary compounds

We outline how pairs of strongly immiscible elements, referred to here as antagonistic pairs, can be used to synthesize ternary compounds with low or quasi-reduced-dimensional motifs intrinsically built into their crystal structures. By identifying third elements that are mutually compatible with a given antagonistic pair, ternary compounds can be formed in which the third element segregates the immiscible atoms into spatially separated substructures. Quasi-low-dimensional structural units, such as sheets, chains, or clusters are a natural consequence of the immiscible atoms seeking to avoid close contact in the solid state. Further, as proof of principle, we present the discovery, crystal growth, and basic physical properties of La 4 ⁢Co 4 ⁢$\mathrm{X}$ (X = Pb, Bi, Sb), a family of intermetallic compounds based on the antagonistic pairs Co-Pb and Co-Bi. La 4 ⁢Co 4 ⁢$\mathrm{X}$ adopts an orthorhombic crystal structure (space group Pbam) containing quasi-two-dimensional Co slabs and La-X polyhedra that stack in an alternating manner along the α axis. Consistent with our proposal, the La atoms separate the Co and X substructures, ensuring there are no direct contacts between the members of the immiscible (antagonistic) pair. Within the Co slabs, the atoms occupy the vertices of corner sharing tetrahedra and triangles, and this bonding motif produces narrow electronic bands near the Fermi level that favor magnetism. The Co is moment bearing in each La 4 ⁢Co 4 $\mathrm{X}$ compound studied, and we show that whereas La 4 ⁢Co 4 ⁢Pb behaves as a three-dimensional antiferromagnet with T N =220K, La 4 ⁢Co 4 ⁢Bi and La 4⁢ Co 4 ⁢Sb have behavior consistent with low-dimensional magnetic coupling and ordering, with T N =153K and 143 K, respectively. In addition to the Pb-, Bi-, and Sb-based La 4 ⁢Co 4 ⁢$\mathrm{X}$ compounds, we also were likely able to produce an analogous La 4 ⁢Co 4 ⁢Sn in polycrystalline form, although we were unable to isolate single crystals. We anticipate that identifying and using mutually compatible third elements together with an antagonistic pair represents a generalizable design principle for discovering new materials and new structure types containing low-dimensional substructures.

36 MATERIALS SCIENCE↗

Unveiling the Finite Temperature Physics of Hydrogen Chains via Auxiliary Field Quantum Monte Carlo

The ability to accurately predict the finite temperature properties and phase diagrams of realistic quantum solids is central to uncovering new phases and engineering materials with novel properties ripe for device applications. Nonetheless, there remain comparatively few many-body techniques capable of elucidating the finite temperature physics of solids from first principles. Here, we take a significant step toward developing such a technique by generalizing our previous, exact fully ab initio finite temperature Auxiliary Field Quantum Monte Carlo (FT-AFQMC) method to model periodic solids and employing it to uncover the finite temperature physics of periodic hydrogen chains. Our chains’ unit cells consist of 10 hydrogen atoms modeled in a minimal basis, and we sample 5 k-points from the first Brillouin zone to arrive at a supercell consisting of 50 orbitals and 50 electrons. Based upon our calculations of these chains’ many-body energies, free energies, entropies, heat capacities, double and natural occupancies, and charge and spin correlation functions, we outline their metal–insulator and magnetic ordering as a function of both H–H bond distance and temperature. At low temperatures approaching the ground state, we observe both metal–insulator and ferromagnetic–antiferromagnetic crossovers at bond lengths between 0.5 and 0.75 Å. We then demonstrate how this low-temperature ordering evolves into a metallic phase with decreasing magnetic order at higher temperatures. In order to contextualize our results, we compare the features we observe to those previously seen in one-dimensional, half-filled Hubbard models at finite temperature and in ground state hydrogen chains. Interestingly, we identify signatures of the Pomeranchuk effect in hydrogen chains for the first time and show that spin and charge excitations that typically arise at distinct temperatures in the Hubbard model are indistinguishably coupled in these systems. Beyond qualitatively revealing the many-body phase behavior of hydrogen chains in a numerically exact manner without invoking the phaseless approximation, our efforts shed light on the further theoretical developments that will be required to construct the phase diagrams of the more complex transition metal, lanthanide, and actinide solids of longstanding interest to physicists.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

In Situ Atomic Tracking on the Interfacial Etching and Reconfiguration of Cu-ReSe 2 Contact during Thermal Annealing

The Schottky barrier height can be greatly affected by the metal diffusion, reaction, and covalent bonding formation at the contact. Exploring novel methods and revealing the fundamental mechanisms for contact engineering are of vital importance for microelectronic devices. Here, in this study, the annealing induced interfacial reactions at Cu-ReSe 2 contact are dynamically revealed from the atomic scale. Accompanied by the diffusion of Se to Cu, ReSe 2 is gradually decomposed to a thin Re interlayer through a “chain-by-chain” manner. Theoretical calculations show that the Cu atoms can facilitate the chemical bond breaking of ReSe 2 , significantly lowering the Se diffusion energy barrier toward Cu. The formed Re/ReSe 2 heterostructure presents a metal-like band structure, which underscores the critical role of Cu in altering the interfacial chemistry and promoting carrier transport across the interface. Our results can provide vital insights into the contact properties of ReSe 2 and provide a possible method for fabricating high-performance ReSe 2 -based devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Fermi-Löwdin-orbital self-interaction correction using the optimized-effective-potential method within the Krieger-Li-Iafrate approximation

Perdew-Zunger self-interaction correction (PZ-SIC) offers a route to remove self interaction errors on an orbital-by-orbital basis. A recent formulation of PZ-SIC by Pederson, Ruzsinszky and Perdew proposes restricting the unitary transformation to localized orbitals called Fermi-Löwdin orbitals. This formulation, called the FLOSIC method, simplifies PZ-SIC calculations and was implemented self-consistently using a Jacobi-like (FLOSIC-Jacobi) iteration scheme. In this work we implement the FLOSIC approach using the Krieger-Li-Iafrate (KLI) approximation to the optimized effective potential (OEP). We compare the results of present FLOSIC-KLI approach with FLOSIC-Jacobi scheme for atomic energies, atomization energies, ionization energies, barrier heights, polarizability of chains of hydrogen molecules etc. to validate the FLOSIC-KLI approach. The FLOSICKLI approach, which is within the realm of Kohn-Sham theory, predicts smaller energy gaps between frontier orbitals due to the lowering of eigenvalues of the lowest unoccupied orbitals. Results show that atomic energies, atomization energies, ionization energy as an absolute of highest occupied orbital eigenvalue, and polarizability of chains of hydrogen molecules between the two methods agree within 2%. Lastly the FLOSIC-KLI approach is used to determine the vertical ionization energies of water clusters.

74 ATOMIC AND MOLECULAR PHYSICS↗

Spectroscopic camera analysis of the roles of molecularly assisted reaction chains during detachment in JET L-mode plasmas

The roles of the molecularly assisted ionization (MAI), recombination (MAR) and dissociation (MAD) reaction chains with respect to the purely atomic ionization and recombination processes were studied experimentally during detachment in low-confinement mode (L-mode) plasmas in JET with the help of experimentally inferred divertor plasma and neutral conditions, extracted previously from filtered camera observations of deuterium Balmer emission, and the reaction coefficients provided by the ADAS, AMJUEL and H2VIBR atomic and molecular databases. The direct contribution of MAI and MAR in the outer divertor particle balance was found to be inferior to the electron–atom ionization (EAI) and electron–ion recombination (EIR). Near the outer strike point, a strong atom source due to the D$^+_2$ -driven MAD was, however, observed to correlate with the onset of detachment at outer strike point temperatures of T e,osp = 0.9 –2.0 eV via increased plasma-neutral interactions before the increasing dominance of EIR at T e,osp < 0.9 eV, followed by increasing degree of detachment. The analysis was supported by predictions from EDGE2D-EIRENE simulations which were in qualitative agreement with the experimental observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimental and Kinetic Modeling Study of 3-Methyl-2-butenol (Prenol) Oxidation

Longer chain alcohols with 4–5 carbon atoms are attractive alternative fuels as they can be derived from biological sources and since their combustion leads to lower exhaust gas levels of NO x and soot compared to commercial fossil fuels. The auto-ignition behavior of fuels that contain both a hydroxyl group and a C=C double bond in their molecular structure is not well established in the literature. Understanding the influence of these functional groups on the ignition behavior of fuels is critical in the development of tailor-made fuels for advanced combustion engines. In this study, ignition delay times of an unsaturated alcohol, 3-methyl-2-butenol (prenol), are measured using a high-pressure shock tube and a rapid compression machine at pressures of 15 and 30 bar at equivalence ratios of 0.5, 1.0, and 2.0 in “air” in the temperature range 600–1400 K. Furthermore, a detailed kinetic model is developed and validated using the new experimental data in this study. In addition, speciation data in a jet-stirred reactor, ignition delay times, and laminar burning velocities available in the literature were also used to validate the new kinetic model. Fuel flux and sensitivity analyses are performed using this new model to determine the important fuel consumption pathways and critical reactions that affect the reactivity of prenol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Volatile Products of the Autoxidation of Poly(ethylenimine) in CO 2 Sorbents

Poly(ethyleneimine) (PEI)-based CO 2 sorbents are a promising material class for use in CO 2 capture applications, particularly direct air capture (DAC), due to their high amine content and CO 2 capacities. The sorbent lifetime is a key uncertainty in the deployment of these materials in such applications, as they oxidize under conditions relevant to candidate DAC process cycles. Here, in this study, we utilize thermogravimetric analysis/differential scanning calorimetry/FTIR spectroscopy to characterize the nature, rate, and quantity of volatile species formed from the oxidation of PEI when supported on mesoporous Al 2 O 3 . We show that NH 3 , CO 2 , and H 2 O are primary volatile species formed from PEI oxidation and that their production rate matches the overall oxidation rate as measured via heat flow from a differential scanning calorimeter. We show that the total quantity of NH 3 evolved is consistent with terminal primary amines being cleaved from the polymer chain and that approximately one hydrogen atom per PEI repeat unit reacts to form H 2 O during the reaction. Finally, metadynamics and quantum chemical simulations performed on a simplified system are used to highlight the likely importance of radicals in the elementary reactions that comprise the oxidation reaction set. The findings here represent the first characterization of the volatile products from the oxidation of PEI in CO 2 sorbents and lead to a greater understanding of the mechanisms behind the oxidative degradation of these materials.

36 MATERIALS SCIENCE↗

Side chain engineering control of mixed conduction in oligoethylene glycol-substituted polythiophenes

A major limitation for polymeric mixed ionic/electronic conductors (MIECs) is the trade-off between ionic and electronic conductivity; changes made that improve one typically hinder the other. In order to address this fundamental problem, this work provides insight into ways that we could improve one type of conduction without hindering the other. Here, we investigated a common oligoethylene glycol side chain polymer by adjusting the oxygen atom content and position, providing structural insights for materials that better balanced the two conduction pathways. The investigated polymer series showed the prototypical conflict between ionic and electronic conduction for oxygen atom content, with increasing oxygen atom content increasing ionic conductivity, but decreasing electronic conductivity; however, by increasing the oxygen atom distance from the polymer backbone, both ionic and electronic conductivity could be improved. Following these rules, we show that poly(3-(methoxyethoxybutyl)thiophene), when blended with lithium bistrifluoromethanesulfonimide (LiTFSI), matches the ionic conductivity of a comparable MIEC [poly(3-(methoxyethoxyethoxymethyl)thiophene)], while simultaneously showing higher electronic conductivity, highlighting the potential of this design strategy. We also provide strategies for tuning the MIEC performance to fit a desired application, depending on if electronic, ionic, or balanced conduction is most important. These results have implications beyond just polythiophene-based MIECs, as these strategies for balancing backbone crystallization and coordinating group interconnectivity apply for all semicrystalline conjugated polymers.

36 MATERIALS SCIENCE↗

Skeletides: A Modular, Simplified Physical Model of Protein Secondary Structure

Three-dimensional (3D) models are essential for visualization and conceptual understanding of complex architectures such as protein structure. Although there is a plethora of software platforms that allow digital depictions of protein structure at an atomic level in silico, physical models are needed to convey an intuitive understanding of biomolecular architecture. However, it is a challenge to represent all the relevant features of proteins in a single physical model due to their sheer structural complexity. Here, we describe a modular protein model that focuses only on representation of the secondary structure—the underlying structural skeleton. The simplified model consists of amino acid units, which can be linked together to reproduce the two most fundamental structural features of protein secondary structure: the relative positions of the amino acid alpha carbon atoms, and the intra-main chain hydrogen bonding pattern. We use 3D printing and magnets to create a set of three modular amino acid building blocks, which when linked together into a chain, can faithfully represent alpha helices, beta sheets, and turns. These simple to make models can be used to quickly assemble the alpha carbon trace of an entire protein domain, which conveys a tactile experience of the complexity of protein skeletal architecture. These models highlight the modularity of protein structure: where a single structural unit, the amino acid, can be linked together to form a larger, regular secondary structure. These models also have a propensity to spontaneously organize into alpha helices and beta sheets, as demonstrated by their ability to autonomously assemble when placed in a circulating water tank. These models provide a missing educational tool to expand knowledge of protein structure, foster deeper insight into protein folding, and inspire greater interest in biomacromolecular architecture.

36 MATERIALS SCIENCE↗

PDBspheres: a method for finding 3D similarities in local regions in proteins

Abstract We present a structure-based method for finding and evaluating structural similarities in protein regions relevant to ligand binding. PDBspheres comprises an exhaustive library of protein structure regions (‘spheres’) adjacent to complexed ligands derived from the Protein Data Bank (PDB), along with methods to find and evaluate structural matches between a protein of interest and spheres in the library. PDBspheres uses the LGA (Local–Global Alignment) structure alignment algorithm as the main engine for detecting structural similarities between the protein of interest and template spheres from the library, which currently contains >2 million spheres. To assess confidence in structural matches, an all-atom-based similarity metric takes side chain placement into account. Here, we describe the PDBspheres method, demonstrate its ability to detect and characterize binding sites in protein structures, show how PDBspheres—a strictly structure-based method—performs on a curated dataset of 2528 ligand-bound and ligand-free crystal structures, and use PDBspheres to cluster pockets and assess structural similarities among protein binding sites of 4876 structures in the ‘refined set’ of the PDBbind 2019 dataset.

59 BASIC BIOLOGICAL SCIENCES↗

Metallization with generic metallo-organic inks

The use and fabrication of metallo-organic films are discussed. Metallo-organic compounds are ones in which a metal is linked to a long chain carbon ligand through a hetero atom such as O, S, N, P or As. Films formed by the thermal decomposition of these metallo-organics are called MOD films. In order that the products of decomposition contain only CO2, H2O, and in rare cases nitrogen compounds, and to avoid S containing products, the use of a set of metallo-organic compounds for ink fabrication where the linking hetero atom was oxygen was pioneered. These links were made from commercially available carboxylates, or synthesized from commonly available reagents. The processing is described and the molecular design criteria are given. The particular carboxylates or amine carboxylates selected were the octoates or neodecanoates, and they are described.

Vest, G. M.↗

Tholins as Coloring Agents on Solar System Bodies

Pre-biotic organic materials appear to be common on many small bodies in the outer Solar System, as evidenced by the color properties of these objects. We report on our continuing study of color properties in connection with the presence of complex organic solids (tholins) among the planets and their satellites, the asteroids, and the trans- Neptunian objects (Kuiper Belt objects). Most small, icy bodies in the Solar System, whether they have high or low surface reflectance (albedo), show a pronounced downward slope in reflectance at wavelengths shorter than approx. 1 micron. This increasing absorption of sunlight toward shorter wavelengths is characteristic of pi-bonds in hydrocarbons having chains or rings of conjugated C atoms. Tholins, which contain polycyclic aromatic and aliphatic hydrocarbons, exhibit these color properties. Using the complex refractive indices of tholins in models of the reflectance spectra of icy bodies in the Solar System, we find that these complex organic materials satisfactorily account for the coloration so widely observed. The new results presented here show that the wide variety of colors of Kuiper Belt objects can be fit very well with tholins, as can the colors of Pluto and Triton. The implications of these fits of Kuiper Belt objects is that complex organic material is created on their surfaces by energetic particle bombardment of native ices, and also may be accreted from external sources. In the cases of Pluto and Triton, photochemistry of their weak N2 + CH4 + CO atmospheres produces complex organic molecules that precipitate to the surface, providing local color.

Cruikshank, D. P.↗

Working Fluids for Increasing Capacities of Heat Pipes

A theoretical and experimental investigation has shown that the capacities of heat pipes can be increased through suitable reformulation of their working fluids. The surface tensions of all of the working fluids heretofore used in heat pipes decrease with temperature. As explained in more detail below, the limits on the performance of a heat pipe are associated with the decrease in the surface tension of the working fluid with temperature, and so one can enhance performance by reformulating the working fluid so that its surface tension increases with temperature. This improvement is applicable to almost any kind of heat pipe in almost any environment. The heat-transfer capacity of a heat pipe in its normal operating-temperature range is subject to a capillary limit and a boiling limit. Both of these limits are associated with the temperature dependence of surface tension of the working fluid. In the case of a traditional working fluid, the decrease in surface tension with temperature causes a body of the liquid phase of the working fluid to move toward a region of lower temperature, thus preventing the desired spreading of the liquid in the heated portion of the heat pipe. As a result, the available capillary-pressure pumping head decreases as the temperature of the evaporator end of the heat pipe increases, and operation becomes unstable. Water has widely been used as a working fluid in heat pipes. Because the surface tension of water decreases with increasing temperature, the heat loads and other aspects of performance of heat pipes that contain water are limited. Dilute aqueous solutions of long-chain alcohols have shown promise as substitutes for water that can offer improved performance, because these solutions exhibit unusual surface-tension characteristics: Experiments have shown that in the cases of an aqueous solution of an alcohol, the molecules of which contain chains of more than four carbon atoms, the surface tension increases with temperature when the temperature exceeds a certain value. There are also other liquids that have surface tensions that increase with temperature and could be used as working fluids in heat pipes. For example, as a substitute for ammonia, which is the working fluid in some heat pipes, one could use a solution of ammonia and an ionic surfactant.

Chao, David F.↗

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↗