Search NASA⌕ Search

SEARCH · Search NASA

Results for “DFT simulations”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Atomic cluster expansion potential for large scale simulations of hydrocarbons under shock compression

We present an Atomic Cluster Expansion (ACE) machine learned potential developed for high-fidelity atomistic simulations of hydrocarbons, targeting pressures and temperatures near and above supercritical fluid regimes for molecular fluids. A diverse set of stoichiometries were covered in training, including 1:0 (pure carbon), 1:4 (methane), and 1:1 (benzene), and rich bonding environments sampled at supercritical temperatures, hydrogen rich, reactive mixtures where metastable stoichiometries arise, including 1:2 (ethylene) and 1:3 (ethane). A high-fidelity training database was constructed by performing large-scale quantum molecular dynamic simulations [density functional theory (DFT) MD] of diamond, graphite, methane, and benzene. A novel approach to selecting structures from DFT MD is also presented, which allows for the rapid selection of unique DFT MD frames from complex trajectories. Comparisons to DFT and experimental data demonstrate that the presented ACE potential accurately reproduces isotherms, carbon melting curves, radial distribution functions, and shock Hugoniots for carbon and hydrocarbon systems for pressures up to 100 GPa and temperatures up to 6000 K for hydrocarbon systems and up to 9000 K for pure carbon systems. This work delivers a potential that can be used for accurate, large-scale simulations of shocked hydrocarbons and demonstrates a methodology for fitting and validating machine learning interatomic potentials to complex molecular environments, which can be applied to energetic materials in future works.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Cesium-Induced Active Sites for C–C Coupling and Ethanol Synthesis from CO 2 Hydrogenation on Cu/ZnO(000$\bar{1}$) Surfaces

The efficient conversion of carbon dioxide, a major air pollutant, into ethanol or higher alcohols is a big challenge in heterogeneous catalysis, generating great interest in both basic scientific research and commercial applications. Here, we report the facilitated methanol synthesis and the enabled ethanol synthesis from carbon dioxide hydrogenation on a catalyst generated by codepositing Cs and Cu on a ZnO(000$\bar{1}$) substrate. A combination of catalytic testing, X-ray photoelectron spectroscopy (XPS) measurements, and calculations based on density functional theory (DFT) and kinetic Monte Carlo (KMC) simulation was used. The results of XPS showed a clear change in the reaction mechanism when going from Cs/Cu(111) to a Cs/Cu/ZnO(000$\bar{1}$) catalyst. The Cs-promoting effect on C–C coupling is a result of a synergy among Cs, Cu, and ZnO components that leads to the presence of CH x and CH y O species on the surface. Furthermore, according to the DFT-based KMC simulations, the deposition of Cs introduces multifunctional sites with a unique structure at the Cu–Cs–ZnO interface, particularly being able to promote the interaction with CO 2 and thus the methanol synthesis predominantly via the formate pathway. More importantly, it tunes the CHO binding strongly enough to facilitate the HCOOH decomposition to CHO via the formate pathway, but weakly enough to allow further hydrogenation to methanol. The fine-tuning of CHO binding also enables a close alignment of a CHO pair to facilitate the C–C coupling and eventually ethanol synthesis. Our study opens new possibilities to allow the highly active and selective conversion of carbon dioxide to higher alcohols on widely used and low-cost Cu-based catalysts.

36 MATERIALS SCIENCE↗
DOE OSTI2020

STRUCTURE AND PROPERTIES STUDY ON ENERGY MATERIALS: THERMOELECTRIC MATERIAL TETRAHEDRITE AND LITHIUM ION CONDUCTOR

Development of efficient energy materials is critical in order to ease the energy demand and reduce our dependence of fossil fuel. Thermoelectric materials are promising due to their capability of generating electrical power by recovering waste heat. The performance of thermoelectric materials is quantified by a dimensionless figure of merit zT, which depends on their properties such as electrical conductivity, Seebeck coefficient and thermal conductivity. Tetrahedrites, a copper antimony sulfosalt mineral, typified by Cu 12-x MxSb 4 S 13 , where M is a transition metal element such as Ni, Zn, Fe or Mn, have great potential for thermoelectric application due to their relatively high zT (close to 1 at 700 K), earth-abundance, environmental friendliness, favorable electrical properties, and most importantly intrinsic low lattice thermal conductivity (less than 1 W m -1 K -1 ) in wide temperature. In addition to energy recovery, reliable energy storage devices are also emerging to relieve the energy demand and improve the efficiency of consuming energy resources. Lithium-ion batteries are known to be reliable and successful electrochemical energy storage devices and appliable in various aspects, including laptops, smartphones and electrical vehicles. Lithium phosphorous oxynitride (LiPON) are widely used as thin-film solid-state electrolytes in Li-ion battery, which is the only demonstrated solid-state electrolyte that is quite stable in direct contact with Li metal at potentials from 0-5 V. However, the structure of LiPON, the effects of N doping, and the origin of its good electrochemical stability remains inconclusive. In this thesis, reliable modeling techniques accompanied with experimental tools, are applied to study the thermoelectric material tetrahedrite and the ionic conductor LiPON, in order to study their structural and dynamical properties. Accurate and efficient density-functional theory (DFT) and density-functional tight-binding (DFTB) methods, combined with molecular dynamics (MD) simulations are utilized in order to investigate the structures and properties of these energy materials. The incoherent and coherent atomic dynamics study of tetrahedrite Cu 10.5 NiZn 0.5 Sb 4 S 13 provides the origin of softening upon cooling by investigate the motion of Cu12e at different temperatures. The dynamic structure factors in the longitudinal and transverse direction will also be discussed. The Cu movement of Cu-rich tetrahedrite Cu 14 Sb 4 S 13 is revealed by Cu self-diffusivity, nuclear density map and “nudged elastic band” (NEB). Moreover, we investigate the effect of simulation cell size and basis sets on the DFT-based MD simulation results using tetrahedrite Cu1 0 Zn 2 Sb 4 S 13 thermoelectric as a model material, showing the advantage of larger cell by accessing smaller Q range. In addition, the low-temperature structural properties of Cu 12 Sb 4 S 13 is measured by neutron diffraction, which indicates that no cubic to tetragonal transition occurs at metal-semiconductor transition (MST) temperature. Thermoelectric properties such as Seebeck coefficient, electrical resistivity and electrical thermal conductivity will also be investigated. DFTB method is implemented to study the structure and transport properties of Li 3 PO 4 and LiPON, while the exploration of N doping effect is included. Lastly, the LiPON/Li interphase will be revealed in order to study the origin of electrochemical stability.

Li, Junchao↗
DOE OSTI2020

STRUCTURE AND PROPERTIES STUDY ON ENERGY MATERIALS: THERMOELECTRIC MATERIAL TETRAHEDRITE AND LITHIUM ION CONDUCTOR LiPON

Development of efficient energy materials is critical in order to ease the energy demand and reduce our dependence of fossil fuel. Thermoelectric materials are promising due to their capability of generating electrical power by recovering waste heat. The performance of thermoelectric materials is quantified by a dimensionless figure of merit zT, which depends on their properties such as electrical conductivity, Seebeck coefficient and thermal conductivity. Tetrahedrites, a copper antimony sulfosalt mineral, typified by Cu12-xMxSb4S13, where M is a transition metal element such as Ni, Zn, Fe or Mn, have great potential for thermoelectric application due to their relatively high zT (close to 1 at 700 K), earth-abundance, environmental friendliness, favorable electrical properties, and most importantly intrinsic low lattice thermal conductivity (less than 1 W m-1 K-1) in wide temperature. In addition to energy recovery, reliable energy storage devices are also emerging to relieve the energy demand and improve the efficiency of consuming energy resources. Lithium-ion batteries are known to be reliable and successful electrochemical energy storage devices and appliable in various aspects, including laptops, smartphones and electrical vehicles. Lithium phosphorous oxynitride (LiPON) are widely used as thin-film solid-state electrolytes in Li-ion battery, which is the only demonstrated solid-state electrolyte that is quite stable in direct contact with Li metal at potentials from 0-5 V. However, the structure of LiPON, the effects of N doping, and the origin of its good electrochemical stability remains inconclusive.In this thesis, reliable modeling techniques accompanied with experimental tools, are applied to study the thermoelectric material tetrahedrite and the ionic conductor LiPON, in order to study their structural and dynamical properties. Accurate and efficient density-functional theory (DFT) and density-functional tight-binding (DFTB) methods, combined with molecular dynamics (MD) simulations are utilized in order to investigate the structures and properties of these energy materials. The incoherent and coherent atomic dynamics study of tetrahedrite Cu10.5NiZn0.5Sb4S13 provides the origin of softening upon cooling by investigate the motion of Cu12e at different temperatures. The dynamic structure factors in the longitudinal and transverse direction will also be discussed. The Cu movement of Cu-rich tetrahedrite Cu14Sb4S13 is revealed by Cu self-diffusivity, nuclear density map and "nudged elastic band" (NEB). Moreover, we investigate the effect of simulation cell size and basis sets on the DFT-based MD simulation results using tetrahedrite Cu10Zn2Sb4S13 thermoelectric as a model material, showing the advantage of larger cell by accessing smaller Q range. In addition, the low-temperature structural properties of Cu12Sb4S13 is measured by neutron diffraction, which indicates that no cubic to tetragonal transition occurs at metal-semiconductor transition (MST) temperature. Thermoelectric properties such as Seebeck coefficient, electrical resistivity and electrical thermal conductivity will also be investigated. DFTB method is implemented to study the structure and transport properties of Li3PO4 and LiPON, while the exploration of N doping effect is included. Lastly, the LiPON/Li interphase will be revealed in order to study the origin of electrochemical stability

25 ENERGY STORAGE↗
DOE OSTI2024

Temperature and density dependent pair potential for deuterium under shock

Large-scale classical molecular dynamics (CMD) simulations naturally include the microscopic physics necessary for atomistic modeling of shock release at the ablator-fuel interface in an inertial confinement fusion (ICF) capsule. Here, the multi-megabar shocks utilized in ICF experiments can drive the deuterium fuel from ambient to electron volt temperatures (T) and multi-fold compression. Modeling interatomic interactions over such an extreme range of conditions is challenging for empirical bond order potentials. We generate a pair potential for deuterium with explicit temperature and mass density dependence from ab initio density functional theory molecular dynamics using the iterative Boltzmann inversion method. This potential accurately reproduces the radial distribution functions and pressures from DFT in CMD equilibrium simulations across a wide range of thermodynamic conditions, yet fails to return the expected Hugoniot relations when used in direct CMD shock simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Structural and Spectroscopic Properties of Butanediol-Modified Boehmite Materials

Glycoboehmite (GB) materials are synthesized by a solvothermal reaction to form layered aluminum oxyhydroxide (boehmite) modified by intercalated butanediol molecules. These hybrid materials offer a platform to design materials with potentially novel sorption, wetting, and catalytic properties. Several synthetic methods have been used, resulting in different structural and spectroscopic properties, but atomistic detail is needed to determine the interlayer structure to explore the synthetic control of GB materials. Here, in this study, we use classical molecular dynamics (MD) simulations to compare the structural properties of GB interlayers containing chemisorbed butanediol molecules as a function of diol loading. Accompanying quantum (density functional theory, DFT) static calculations and MD simulations are used to validate the classical model and compute the infrared spectra of various models. Classical MD results reveal the existence of two unique interlayer environments at higher butanediol loading, corresponding to smaller (cross-linked) and expanded interlayers. DFT-computed infrared spectra reveal the sensitivity of the aluminol O–H stretch frequencies to the interlayer environment, consistent with the spectrum of the synthesized material. Insight from these simulations will aid in the characterization of the newly synthesized GB materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Experimental and Theoretical Studies of the Surface Oxidation Process of Rare–Earth Tritellurides

Recent studies have established Van der Waals (vdW) layered and 2D rare-earth tritellurides (RTe 3 ) as superconductors and near room-temperature charge density wave (CDW) materials. Their environmental stability raises natural concern owing to aging/stability effects observed in other tellurium-based layered crystals. Here, the results establish the stability and environmental aging characteristics of these RTe 3 systems involving a variety of metals such as La, Nd, Sm, Gd, Dy, and Ho. The atomic force microscopy (AFM) and scanning electron microscopy (SEM) results show that all the RTe 3 sheets oxidize to form thin TeO x layers that are primarily confined to the surface, edges, and grain boundaries. Time-resolved in situ Raman spectroscopy measurements are used to understand the kinetics of the oxidization process for different lanthanide metal cations and establish their relative stability/resilience to oxidization. Overall results indicate that the vdW layers show higher air stability as the 4f electron number decreases going from Ho to La, resulting in the most stable LaTe 3 compared to the least stable HoTe 3 . Comprehensive quantum mechanical simulations reveal that environmental degradation originates from a strong oxidizing reaction with O 2 molecules, while humidity (H 2 O) plays a negligible role unless Te vacancies are present. Moreover, the simulations explain the effects of 4f electrons on the work function and Te vacancies formation, which directly impact the aging characteristics of RTe 3 layers. Interestingly, optical and electrical measurements show that the CDW response is still observed in aged RTe 3 layers owing to the presence of underlying pristine/nonoxidized RTe 3 layers, except CDW transition temperatures increase due to the thickness effect. Overall results offer the first in-depth environmental aging studies on these materials, which can be applied to engineer and design their chemical stability, surface properties, and overall CDW characteristics.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Understanding Solid-Electrolyte Interphase Formation at the Lithium Metal Anode of Lithium-Sulfur Batteries

Lithium-sulfur (Li-S) batteries are one of the most promising chemistries for the next generation of energy storage devices. In a standard Li-S cell, the metallic lithium anode is an essential component of the battery due to its low density, extremely high theoretical specific capacity, and very low negative electrochemical potential. However, several challenges related to the use of Li metal have prevented this battery technology from becoming commercially available. For instance, Li-anodes are highly reactive, which results in the continuous decomposition of the electrolyte and the formation of the solid-electrolyte interphase (SEI) layer. This, combined with the non-uniform deposition of Liions during plating at the anode surface, can seriously affect the performance, cycling, and safety of the battery. In addition, it has been suggested that a controlled SEI formation at the metallic Li anode can yield enhanced battery performance. Hence, a comprehensive understanding of how the SEI layer is formed can help elucidate improvements in this battery technology. In this work, the formation of the SEI layer in Li-S batteries is investigated using density functional theory (DFT) and ab initio molecular dynamics (AIMD) simulations. The stability of liquid organic-based electrolyte components, salt concentration, electron-rich environments, and the use of inorganic solid-state electrolyte (SSE) materials are explored in order to provide molecular-level fundamental insights into how the nature and composition of the electrolyte can alter the initial stages of the SEI formation. The products of electrolyte decomposition and electro-/chemical reaction pathways of liquid electrolyte component such as solvent and salts are predicted from AIMD simulations of electrolyte-anode interfaces. Energetics of reactions from DFT ground-state optimizations are also presented to confirm the decomposition mechanisms. Some significant differences are then drawn regarding the use of low and high salt concentrations. Simulations under electron-rich environments show additional multielectron electrochemical reactions of solvent and salt decomposition taking place due to the excess of electrons and the presence of radical anions in the solution. Finally, the initial stages of formation, stability, and main constituents of interfaces between Li-metal and Sbased electrodes with sulfide-based SSEs are characterized in detail.

25 ENERGY STORAGE↗
DOE OSTI2021

Effect of excess Mg to control corrosion in molten MgCl 2 and KCl eutectic salt mixture

Structural alloys may experience corrosion when exposed to molten chloride salts due to selective dissolution of active alloying elements. One way to prevent this is to make the molten salt reducing. For the KCl + MgCl 2 eutectic salt mixture, pure Mg can be added to achieve this. However, Mg can form intermetallic compounds with nickel at high temperatures, which may cause alloy embrittlement. This work shows that an optimum level of excess Mg could be added to the molten salt which will prevent corrosion of alloys like 316 H, while not forming any detectable Ni-Mg intermetallic phases on Ni-rich alloy surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Probing the Mechanism of Cadmium Sulfide Cluster Nucleation and Growth via Sequential Infiltration Synthesis

Few-atom metal chalcogenide clusters may be realized through the sequential infiltration of metal–organic precursors in polymer films. However, the underlying nucleation and growth mechanisms that allow for cluster synthesis with near atomic-scale precision are not fully resolved. The kinetics of the sequential infiltration synthesis (SIS) method that control the nucleation and growth mechanisms of primarily Cd 4 S 4 -core clusters within a poly­(4-vinylpyridine) (P4VP) matrix are probed with in situ UV–visible absorbance spectroscopy. Density functional theory (DFT) calculations that allow simulation of the optical properties of cluster fragments further reveal the thermodynamics that guide cluster growth within the P4VP matrix. Here, we conclude that a reactive capture mechanism for cluster nucleation and growth is dominant, although transient dimethyl cadmium adduction to the polymer backbone may contribute to cluster nucleation under shorter metal–organic purge process conditions. Grazing incidence X-ray diffraction (GI-XRD) and X-ray absorption spectroscopy (XAS) analyses further corroborate the cluster size and atom connectivity throughout the stepwise synthesis.

Jayaweera, Nuwanthaka P. [Argonne National Laborat↗
DOE OSTI2023

A Membrane Contactor Enabling Energy-Efficient CO 2 Capture from Point Sources with Deep Eutectic Solvents

We demonstrate a scalable and energy-efficient hollow fiber membrane contactor (HFMC)-based process using a green solvent for CO 2 capture. This process uses a deep eutectic solvent (DES) in an HFMC to provide close interfacial interactions and contact between the DES and CO 2 . This approach overcomes disadvantages associated with direct absorption in DES and could potentially be applied to a variety of solvent-based CO 2 capture methods. Commercial low-cost polymer hollow fiber membranes (e.g., microporous polypropylene) were evaluated for CO 2 capture with reline, a prototypical DES. Single-gas measurements showed that the DES-based polypropylene HFMC can capture and separate CO 2 while rejecting N2. From a mixed gas containing 50 mol % N 2 and 50 mol % CO ν , the DES-based HFMC separated CO 2 with a purity of 96.9 mol %. The effect of several process parameters including solvent flow rate, pressure, and temperature on the CO 2 separation performance was studied. The flux of the recovered CO 2 was 67.43 mmole/m 2 /h at a feed pressure of 4 bar. In situ Fourier transform infrared (FTIR) measurements combined with density functional theory (DFT)-based molecular dynamics simulations revealed that reline absorbs CO 2 by physical absorption without forming a new chemical compound, and CO 2 separation by reline occurs via the pressure swing mechanism. This research provides fundamental insights about physical solvent-based separation processes and a pathway toward practical deployment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Unraveling Adsorbate-Induced Structural Evolution of Iron Carbide Nanoparticles

Iron carbide (Fe x C y ) nanoparticles (NPs) are promising candidates for replacing platinum group metals in industrial applications, such as high-temperature Fischer–Tropsch synthesis. However, due to their amorphous nature, characterization of the active sites has been challenging experimentally and computationally. Here, using a combined density functional theory (DFT), neural network interatomic potential-assisted global optimization, and ensemble learning study, we evaluate dynamic surface changes associated with syngas (H and CO) interactions. For this purpose, we have developed a general procedure that we use to model an experimentally relevant 270-atom Fe 182 C 88 NP using the neural network-assisted stochastic surface walk global optimization algorithm (SSW-NN). Once generated, the Fe 182 C 88 NP active sites and particle morphology are thoroughly characterized before the effects of syngas adsorbate interactions are explored by using DFT and molecular dynamics simulations. Lastly, we explore correlations between geometric and electronic features of the active sites and the adsorption of H (H ads ), using a regularized random forest machine learning algorithm. In doing so, we identified the Fe–C coordination number and p orbital occupancy as the most important descriptors affecting H ads . Furthermore, using a combined ML and quantum chemistry approach, our work demonstrates a general and efficient procedure for generating and probing complex surface phenomena on binary nanoparticles.

Adsorption↗
DOE OSTI2025

Operando Characterization of Fe in Doped Ni x (Fe 1– x )O y H z Catalysts for Electrochemical Oxygen Evolution

Iron-doped nickel oxyhydroxides, Ni x (Fe 1– x )O y H z , are among the most promising oxygen evolution reaction (OER) electrocatalysts in alkaline environments. Although iron (Fe) significantly enhances the catalytic activity, there is still no clear consensus on whether Fe directly participates in the reaction or merely acts as a promoter. To elucidate the Fe’s role, we performed operando X-ray spectroscopy studies supported by DFT on Ni x (Fe 1– x )O y H z electrocatalysts. We probed the reversible changes in the structure and electronic character of Ni x (Fe 1– x )O y H z as the electrode potential is cycled between the resting (here at 1.10 V RHE ) and operational states (1.66 V RHE ). DFT calculations and XAS simulations on a library of Fe structures in various NiO y H z environments are in favor of a distorted local octahedral Fe(III)O 3 (OH) 3 configuration at the resting state with the NiO y H z scaffold going from α-Ni(OH) 2 to γ-NiOOH as the potential is increased. Under catalytic conditions, EXAFS and HERFD spectra reveal changes in p-d mixing (covalency) relative to the resting state between O/OH ligands and Fe leading to a shift from octahedral to square pyramidal coordination at the Fe site. XES measurements and theoretical simulations further support that the Fe equilibrium structure remains in a formal Fe(III) state under both resting and operational conditions. These spectral changes are attributed to potential dependent structural rearrangements around Fe. The results suggest that ligand dissociation leads to the C 4v symmetry as the most stable intermediate of the Fe during OER. This implies that Fe has a weakly coordinated or easily dissociable ligand that could serve to coordinate the O–O bond formation and, tentatively, play an active role in the Ni x (Fe 1– x )O y H z electrocatalyst.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Breaking up CO 2 through pressure-induced redox reaction with rhenium

We present results from density-functional-theory (DFT) based molecular dynamics simulations of fluid mixtures of carbon dioxide (CO 2 ) and rhenium (Re) metal at pressures up to 45 GPa and a temperature of 4000 K. Covalently-bonded carbon chains form as a result of a redox reaction between CO 2 and Re, whereas a pure CO 2 fluid is stable against decomposition at the same conditions. Furthermore, these findings have implications in reducing CO 2 by pressure-induced redox reactions with transition metals, which may play a role in CO 2 destruction and carbon storage as related to possible solutions to the current climate crisis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Multiphase equation of state for Ta 2 O 5

A new equation of state for Ta 2 O 5 is presented. The EOS is constructed using the OpenSesame software and is referred to as SESAME 3530. The EOS uses a combination of density functional theory (DFT) calculations and experimental data. DFT calculations include cold curves and phonons of the solid phases, as well as DFT-based molecular dynamics simulations of the liquid phase. Experimental data includes isobaric, diamond anvil cell, and porous shock Hugoniot data. To fit the data, we create a multiphase EOS consisting of two solid phases and the liquid. Overall agreement with experimental data is shown, and we provide some suggestions for future experiments that could improve our knowledge of the phase diagram.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRSDate not supplied

Building and Breaking Carbon Composites with REACTER

Carbon composites have become indispensable for aerospace and other high-performance applications, and a detailed picture of their morphology and failure mechanisms remains difficult to obtain through experiment. REACTER is a versatile computational modeling tool for atomistic molecular dynamics designed to model chemical reactions at the speed and length scales of classical force fields.1 In this work, several recent features of REACTER were applied to the creation and subsequent mechanical testing of two classes of carbon composites, carbon-fiber reinforced polymers (CFRP) and carbon nanotube (CNT) composites. Carbon fiber core morphologies were created by the method of Desai et al.,2 but using the advanced reaction constraints framework of REACTER, their proposed multistep procedure was reduced to a single uninterrupted molecular dynamics simulation. The carbon fiber filler was embedded into a polymer matrix by simulated in situ polymerization of several thermosetting resins, including bismaleimide and polyarylacetylene, to obtain the final CFRP model. To generate the second class of carbon composite, CNT networks were grown dynamically using the new ‘create atoms’ feature of REACTER, and similarly infiltrated with resin to obtain CNT composites. The resulting models were compared directly to experiment using simulated high resolution transmission electron microscopy and x-ray diffraction. Failure mechanisms were elucidated by simulating mechanically induced bond breaking, as characterized by third-order DFT-based tight-binding (DFTB3) simulations, via a reaction constraint on the total potential energy of the involved atoms.

Molecular Dynamics↗
DOE OSTI2023

Understanding Microscopic Mechanisms of LeTID and LID and their Unifying Features by Electron Paramagnetic Resonance

LID degradation involves not only creation of ~ 10^12 cm^-3 recombination centers, but also ~ 10^16 cm^-3 shallow negative-U traps. In Ga-doped Si, LID EPR defects don't appear, but some traps are still created. LeTID: Si DB and H-hyperfine EPR signatures. We postulate that the defect responsible for LeTID is a partially hydrogenated (multivacancy) with a Si dangling bond and H in the vicinity. O involvement is possible yet unclear. We prove that H is related to the structure of the LeTID defect with isotope experiments and its EPR signal is comparable and linear with the Si DB signal upon LeTID degradation. Working on simulating these results with DFT to obtain more detailed defect structure.

Cz Si↗