Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic methods”

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 505 records · Page 28

A neural master equation framework for multiscale modeling of molecular processes: application to atomic-scale plasma processes

Plasma-surface interactions (PSI) play a crucial role in microelectronics fabrication; however, their multiscale nature and array of complex, often unknown interactions make computational modeling of PSIs extremely difficult. To this end, we propose a general neural master equation (NME) framework that uses master equations to describe the dynamics of a molecular process, wherein neural networks learned from atomistic simulations represent unknown transitions between different system states. By leveraging the physics-based structure of master equations and data-driven state transitions, the NME framework promotes generalizability and physics interpretability, and can bridge disparate length and time scales. The framework is demonstrated for multiscale modeling of Si atomic layer etching and reactive ion etching, where the learned NME-based surface kinetic models exhibit good predictive and extrapolative capabilities for predicting experimentally relevant observables as a function of process parameters. The NME-based surface kinetic models obey physical constraints, which are violated in models based on neural ordinary differential equations. The proposed NME framework for multiscale modeling of molecular processes can pave the way for the discovery of new chemistries and materials in atomic-scale plasma processes.

Chemical engineering↗

A deep learning approach to fast analysis of collective Thomson scattering spectra

Fast analysis of collective Thomson scattering ion acoustic wave features using a deep convolutional neural network model is presented. The network was trained from spectra to predict the plasma parameters, including ion velocities, population fractions, and ion and electron temperatures. A fully kinetic particle-in-cell simulation was used to model a laboratory astrophysics experiment and simulate a diagnostic image of the ion acoustic wave feature. Network predictions were compared with Bayesian inference of the plasma model parameters for both the simulated and experimentally measured images. Both approaches were fairly accurate predicting the simulated image and the network predictions matched a good portion of the Bayesian results for the experimentally measured image. The Bayesian approach is more robust to noise and motivates future work to train deep learning models with realistic noise. The advantage of the deep learning model is making thousands of predictions in a few hundred milliseconds, compared to a few seconds to minutes per prediction for the optimization and Bayesian approaches presented here. The results demonstrate promising capabilities of deep learning models to analyze Thomson data orders of magnitude faster than conventional methods when using the neural network for standalone analysis. If more rigorous analysis is needed, neural network predictions can be used to quickly initialize other optimization methods and increase chances of success. This is especially useful when the dataset becomes very large or highly dimensional and manually refining initial conditions for the entire dataset are no longer tractable.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Control of core–shell nanoparticles properties through plasma synthesis: a computational study

The improved properties of core–shell nanoparticles (CSNPs) over homogeneous nanoparticles (NPs) have expanded and diversified the applications of these nanomaterials. However, controlling the properties of CSNPs can be a challenging task. Low temperature plasmas have proven to be an effective method of producing NPs with uniform size and morphology, and high yield. That said, NP transport and growth dynamics are sensitive to LTP properties. We report on a computational investigation of the evolution of Ge–Si CSNP properties as a function of operating conditions through the modeling of a flowing, two-zone inductively coupled plasma (ICP) reactor. Ar/GeH 4 and Ar/SiH 4 gas mixtures were supplied to separate plasma zones at a pressure of 1 Torr to promote growth of Ge cores and Si shells. The negatively charged CSNPs are trapped electrostatically in the vicinity of the antennas where the plasma is generated and where the majority of particle growth occurs. Particles that grow to a critical size are then de-trapped by fluid drag due to neutral gas flow. A two-dimensional hybrid plasma model coupled with a three-dimensional kinetic NP transport model were utilized to resolve plasma chemistry and NP growth processes that take place on distinct timescales. The trends in CSNP properties and trapping mechanisms associated with flow rate, applied ICP power and inlet precursor fraction are discussed. While the spatial distribution of plasma produced radical species can have significant impact on the NP growth process, the NP transport dynamics are what ultimately dictates the growth environment that is unique to each particle and so determines their final dimension and composition. The key to optimizing reactor conditions involves controlling the spatial density of growth species and plasma profile as a means to tailor particle trapping dynamics suitable to produce CSNPs for a specific application.

36 MATERIALS SCIENCE↗

Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM SciDAC) (Technical Final Report)

Runaway electrons can severely damage the plasma facing components on ITER during a major disruption and pose a major risk for tokamak fusion. It has been recognized that an adequate disruption mitigation system (DMS) is essential for the safe operation of ITER. The United States is responsible for the design and implementation of the disruption mitigation system on ITER, and in July 2016 the Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM) was launched by DOE, in a joint Fusion Energy Sciences (FES) and Advanced Scientific Computing Research (ASCR) collaboration. SCREAM was a comprehensive theory and simulation SciDAC center that provided physics guidance in the avoidance and mitigation of runaway electrons, and in tandem with domestic and international experiments, helped establish the qualitative and quantitative bases for safe operational scenarios and viable mitigation techniques. The SCREAM center assembled a national team of experts in runaway electron physics, tokamak disruptions, magnetohydrodynamic (MHD) simulation, and advanced algorithms and computing. The team combined advanced simulation and analysis capability facilitated by direct participation of ASCR SciDAC institutes with theoretical models and code development by FES scientists to focus on the runaway risk for ITER and tokamaks in general. The research scope was focussed on integrated simulations of kinetic runaway electrons, including MHD and fluid models of impurity transport, within a research plan guided by theory. The specific research tasks were (1) establish the fundamental physics of runaway generation, saturation, and dynamical evolution in a tokamak; (2) examine the critical path toward runaway avoidance; and (3) investigate the viability and effectiveness of the leading candidate schemes for runaway mitigation. In all three areas, members of the team carried out scoping studies that established the readiness for rapid and critical advances, especially in the deployment and further development of large-to extreme-scale simulation tools. Our multi-pronged computational approach included (1) relativistic Fokker-Planck solvers with discretization in phase space, (2) self-consistent particle-in-cell techniques, (3) particle-based Monte-Carlo, and (4) MHD-particle hybrid simulations. Cross-check between these different methods provided an additional means for verification and further bolstered the fidelity of our physics prediction. Validation against experimental results brings confidence to the predictive capability for ITER and frequently leads to new ideas for understanding and mitigating the thermal quench driven runaway electron phenomenon.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Solutions to Droop and the Green Gap by Novel Carrier Injection

Long-wavelength (green, yellow, red) visible light-emitting diodes (LEDs) have historically been inefficient compared to blue light-emitters in the III-nitrides. This is due to many factors, including poor material quality due to the low growth temperatures required for high Incomposition InGaN quantum wells (QWs), and the polar nature of the III-nitride crystal structure. Spontaneous and piezoelectric polarization in the III-nitrides causes enhanced quantum-confined Stark effect (QCSE), which spatially separates electrons and holes and reduces radiative recombination (ie: light output) from the devices. Polarization in the crystal also induces additional potential barriers that electrons and holes must cross over before entering the QWs where they can recombine to emit light. Both of these polarization effects worsen as emission wavelength increases, making blue III-nitride LEDs are the most efficient and red III-nitride LEDs the least efficient. The focus of this project has been solving the latter polarization-related issue through the implementation of V-defects in long wavelength LEDs. V-defects are morphological defects which are commonly observed in c-plane III-nitrides. They are observed as hexagonal pyramid-shaped depressions on the c-plane surface, with six semipolar sidewalls. They typically form at the apex of threading dislocations (TDs) under conditions of kinetically-limited growth and low growth-temperature. They were initially thought to be detrimental to LED performance and much early work focused on eliminating them entirely from III-nitride devices. However, over the past decade work has emerged that indicate that they can improve LED performance by allowing electrons and holes to bypass the polarization-induced barriers present in the c-plane and directly enter the QWs of an LED. This is due to the semipolar nature of the V-defect sidewall: these sidewalls are thin and lack the polarization-induced barriers which prevent carriers easily moving between layers. V-defects have since been determined to be an efficient avenue by which to inject electrons and holes into the c-plane QWs where they can recombine to emit light. Throughout this project we have explored lateral injection through a variety of methods: simulation (Task 2), epitaxial growth of V-defect and non-V-defect LEDs (Tasks 1, 3, 4, 5, 6), and advanced characterization methods (Task 7). All tasks have been completed. A description of each task completed follows from this section.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Evaluating the Origins of Aerobic Oxidation Catalysis with TAM-3, a MOF with Accessible Co(II) Sites and Large Pores

Metal-organic frameworks (MOFs) are attractive platforms that merge concepts of homogeneous and heterogeneous catalysis. Catalyst design and optimization are enabled by an array of synthetic methods that offer independent control over the local chemical structure of lattice-embedded metal ions (i.e., ligand identity and geometry) and the long-range materials properties (i.e., porosity). Establishing the origin of catalytic activity in MOF-promoted reactions remains a significant challenge: The relative rates of catalyst turnover and substrate diffusion dictate the extent to which interstitial sites are accessible and operational in catalysis. To minimize the contributions of surface sites in catalysis, materials with large pore dimensions are often sought, however, the impact of pore expansion on the origins of catalytic activity is similarly challenging to establish. Here, we describe TAM-3, a Co(II) based MOF with accessible metal sites supported by a facially coordinating tris-tetrazole ligand set. TAM-3 features large channel-like pores (17 × 23 Å) and promotes aerobic C−H oxidation and olefin epoxidation. Using a set of simple kinetics experiments, based on the analysis of kinetic isotope effects and olefin oxidation diastereoselectivities, we demonstrate that despite the large pores, interstitial metal ions do not significantly contribute to the observed substrate oxidation. This study highlights the importance of conducting kinetic experiments to assess the origin of apparent catalytic activity with MOFs and the challenge of harnessing reactive oxidants with microporous catalyst materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Experimental electronic structures of the Fe IV =O bond in S=1 heme vs. nonheme sites: Effect of the porphyrin ligand

High-valent Fe IV =O species are common intermediates in biological and artificial catalysts. Heme and nonheme S=1 Fe IV =O sites have been synthesized and studied for decades but little quantitative experimental comparison of their electronic structures has been available, due to the lack of direct methods focused on the iron. This study allows a rigorous determination of the electronic structure of a nonheme Fe IV =O center and its comparison to an Fe IV =O heme site using 1s2p resonant inelastic X-ray scattering (RIXS) and Fe L-edge X-ray absorption spectroscopy (XAS). Further, variable temperature magnetic circular dichroism (VT-MCD) of the ligand field transitions, combined with nuclear resonance vibrational spectroscopy of the two S=1 Fe IV =O systems show that the equatorial ligand field decreases from a nonheme to a heme Fe IV =O site. Alternatively, RIXS and Fe L-edge XAS combined with MCD show that the Fe dπ orbitals are unperturbed in the Fe IV =O heme relative to the nonheme site because the strong axial Fe-O bond uncouples the Fe dπ orbitals from the porphyrin π-system. As a consequence, the thermodynamics and kinetics of the H-atom abstraction reactions are actually very similar for heme compound II and nonheme Fe IV =O active sites.

bioinorganic chemistry↗

Large-scale Linear Magnetic Holes with Magnetic Mirror Properties in Hybrid Simulations of Solar Wind Turbulence

Magnetic holes (MHs) are coherent magnetic field dips whose size ranges from fluid to kinetic scale, ubiquitously observed in the heliosphere and in planetary environments. Despite the long-standing effort in interpreting the abundance of observations, the origin and properties of MHs are still debated. In this Letter, we investigate the interplay between plasma turbulence and MHs, using a 2D hybrid simulation initialized with solar wind parameters. We show that fully developed turbulence exhibits localized elongated magnetic depressions, whose properties are consistent with linear MHs frequently encountered in space. The observed MHs develop self-consistently from the initial magnetic field perturbations by trapping hot ions with large pitch angles. Ion trapping produces an enhanced perpendicular temperature anisotropy that makes MHs stable for hundreds of ion gyroperiods, despite the surrounding turbulence. We introduce a new quantity, based on local magnetic field and ion temperature values, to measure the efficiency of ion trapping, with potential applications to the detection of MHs in satellite measurements. We complement this method by analyzing the ion velocity distribution functions inside MHs. Our diagnostics reveal the presence of trapped gyrotropic ion populations, whose velocity distribution is consistent with a loss cone, as expected for the motion of particles inside a magnetic mirror. Our results have potential implications for the theoretical and numerical modeling of MHs.

79 ASTRONOMY AND ASTROPHYSICS↗

Laser-to-proton conversion efficiency studies for proton fast ignition

We assess the conversion efficiency from intense picosecond laser pulses to multi-MeV ion beams for a wide range of laser and target parameters, using 2D kinetic particle-in-cell simulations. Scalings are addressed in a quasi-one-dimensional geometry, leaving out beam divergence. Then, we study the conversion efficiency into a narrow spatial band along the laser axis for flat and hemispherical targets in large-scale 2D simulations. Combining these findings allows us to calculate the energy requirements for ignition of a compressed ICF target with an intense proton beam in a fast-ignition scenario.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quantification of [MTBD][beti] loading and its effects on Pt/HSC electrode performance in hydrogen fuel cells

High surface area carbon-supported platinum catalysts (Pt/HSC) are widely used in polymer electrolyte membrane fuel cells but often suffer from limited proton and oxygen transport within porous domains. To address these challenges, we integrate the ionic liquid [MTBD][beti] into Pt/HSC catalyst layers using two deposition methods-one-pot and sequential deposition-to tune IL distribution within micropores and mesopores. A combination of ex situ and in operando techniques were employed to elucidate electrode structure-performance relationships across a range of IL loadings. Sequential deposition achieved more efficient pore filling and higher IL retention at lower IL:C ratios, enabling enhanced proton conductivity and protection of active sites. Compared to the IL-free Pt/HSC, IL-modified electrodes demonstrated up to 28% improvement in mass activity and enhanced high-current-density performance under low relative humidity, while maintaining comparable performance under humidified conditions. Electrochemical impedance spectroscopy and CO displacement experiments reveal that improvements are linked to reduced ionic resistance and lower sulfonate adsorption on Pt sites, rather than changes in electrochemical surface area. However, excessive IL loading leads to mass transport losses. These results highlight the importance of selective pore filling and efficient IL distribution in achieving kinetic gains without compromising ionic and gas transport.

08 HYDROGEN↗

Interactions of clathrate hydrate promoters sodium dodecyl sulfate and tetrahydrofuran investigated using 1H diffusion nuclear magnetic resonance at hydrate-forming conditions

Thermodynamic hydrate promoters and kinetic hydrate promoters can be used to reduce the P–T conditions for clathrate hydrate synthesis to decrease the nucleation induction time while increasing growth rates. Two commonly used promoters for hydrate research are tetrahydrofuran (THF) and sodium dodecyl sulfate (SDS), which can increase the overall hydrate promotion when used in tandem as compared to individually. There are several molecular theories regarding how SDS promotes hydrate growth. This study explores the micellular theory, for which hydrate formation depends on surfactant aggregates (micelles) at a critical micelle concentration (CMC) to increase the interfacial surface area. The micellular theory is the most investigated and criticized surfactant hydrate promotion theory. To address questions related to micellar behavior, this study investigates the intermolecular behavior between SDS and THF for the identification of micelles at hydrate-forming conditions. The systems explored contained THF at 3 and 5 wt. % with varying concentrations of SDS below and above the CMC. Several methods including a qualitative visual method, conductivity, interfacial tensiometry, 13C Liquid-state Nuclear Magnetic Resonance (NMR) spectroscopy, and 1H diffusion NMR spectroscopy were evaluated at temperatures below the Krafft point of SDS and above 0 °C. The presence of THF at low concentrations decreased the critical temperature for the formation of SDS micelles, where SDS is solubilized in THF/water solution at hydrate-forming temperatures without precipitation. The CMC of SDS was decreased significantly even at hydrate-forming conditions. Mixed surfactant–cosolvent micellular behavior of SDS in the presence of low concentrations of THF was confirmed at hydrate-forming conditions above 0 °C.

Chemistry↗

Continuous Wet Air Oxidation of the Hydrothermal Liquefaction Aqueous Product from Various Wet Wastes

Wet air oxidation (WAO) offers an effective method for treating waste streams, converting pollutants into benign substances, and holds significant potential for processing the aqueous product from the hydrothermal liquefaction (HTL-AP) of wet wastes, a promising renewable fuel technology. Here, we conducted a comprehensive study of the WAO of HTL-AP from four different wet wastes. Through continuous testing under various conditions, we produced samples with different chemical oxygen demand (COD) levels, enhancing understanding of reaction parameters necessary for substantial COD reduction (>95%). Chemical analysis revealed that alcohols and ketones in the HTL-AP rapidly oxidized to acetic acid through aldehyde intermediates, while acetic acid, other carboxylic acids, and phenols oxidized relatively slowly. The light N-containing compounds were found to exhibit a change in concentration only after the whole sample reaches an 80% COD reduction, indicating their refractory nature under applied conditions. Energy released in the WAO reaction was calculated, and anaerobic toxicity assay demonstrated that WAO treatment enhanced methane production kinetics due to reduced inhibitory effects, suggesting partial oxidative transformation of inhibitory compounds into less toxic derivatives. These findings provide insights into designing effective WAO processes for valorizing HTL aqueous products, addressing key barriers to HTL process commercialization.

anaerobic digestion↗

A spatially-resolved model of neutron-irradiated tungsten coupling stochastic cluster dynamics and finite deformation plasticity

Structural materials used in nuclear reactors face severe degradation in mechanical properties, such as hardening and embrittlement. At the microscopic scale, this occurs due to creation and accumulation of irradiation-induced defects and their interaction with system dislocations. Although techniques exist which can model evolution of irradiation defects, for instance kinetic transport theory-based models, their interaction with mechanical deformation of the bulk material has not been investigated extensively. In this work, we demonstrate a novel spatially-resolved multiscale coupling between microscopic irradiation defect evolution, modeled using Stochastic Cluster Dynamics (SCD) and macroscopic mechanical deformation modeled using a finite-deformation plasticity model. SCD is used to determine the statistically averaged defect cluster spacing, dependent on operating conditions such as irradiation dose and temperature. This acts as an initial condition that governs the critical resolved shear stress of dislocation glide in the macroscopic plasticity model. This framework is used to predict mechanical behavior in post-mortem test of irradiated Tungsten samples, which has found its importance as structural material used in nuclear reactors. The results obtained using the coupled approach are in good agreement with experimental data of uniaxial tension tests. The model is able to capture the effect of temperature and irradiation dose on the material hardening. Two methods are proposed to estimate hardness – using Tabor's Law relating uniaxial yield stress to hardness and from flat-punch simulations. The results are in reasonable agreement with hardness data from micro-indentation experiments of irradiated Tungsten samples. Finally, the model is also able to reveal microstructural details such as spatial variation in defect density and local stress.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Scalable mechanochemical synthesis of high-quality Prussian blue analogues for high-energy and durable potassium-ion batteries

Prussian blue analogues (PBAs) are recognized as promising cathode materials for potassium-ion batteries (PIBs), particularly the low-cost and high-energy K 2 Mn[Fe(CN) 6 ](KMnF). However, conventional solution-based synthesis inevitably introduces [Fe(CN) 6 ] 4− defects and lattice water while suffering low synthesis efficiency, unfavorable to the improvement of electrochemical performance and scalability. Here, in this work, we report a simple solvent-free mechanochemical strategy for the synthesis of a wide variety of K 2 M[Fe(CN) 6 ] (M = Mn, Mg, Ca, etc.) with negligible defects and water, and it is unprecedented to achieve kilogram-level products of high-quality KMnF within just 10 minutes. The as-prepared KMnF delivers a high energy density of 590 Wh kg −1 at 0.2 C and exhibits an astonishing stability over 10 000 cycles and rate ability up to 50 C in a potassium metal half-cell. Encouragingly, a high-areal-capacity pouch cell with 2.2 mAh cm −2 (16.5 mg cm −2 ) exhibits a capacity retention of 80.7% after 500 cycles. Furthermore, systematic in situ characterization reveals underlying mechanism insights into structure–performance relationships. Specifically, the fully coordinated Mn–N 6 octahedral configuration effectively suppresses Mn 3+ Jahn–Teller distortion, enabling reversible phase transitions under both high-voltage and long-term cycling conditions. In addition, minimal defects provide sufficient redox centers, while the continuous three-dimensional framework facilitates rapid K + diffusion kinetics. This work provides a new opportunity for the ultrafast, universal and scalable synthesis of high-quality PBAs, facilitating the practical application of PIBs while enabling precise structural and compositional design of novel PBAs.

Mechanochemical method↗

Contact with Electrically Conductive Inert Solids Alters Intrinsic Heterogeneous Brønsted Acid Catalysis

Interfacial electric fields at heterogeneous catalyst surfaces have been demonstrated to alter kinetics of liquid-phase reactions. In these systems, electric fields are generated from applying a potential to the catalyst through connection to a potentiostat or through electron transfer from redox-active species in solution. Here, we demonstrate that catalyst polarization can also occur by simply contacting electrically conductive inert solids, leading to the counterintuitive conclusion that a catalyst particle touching an inert solid can alter intrinsic reaction rates. Using dehydration of 1-methylcyclopentanol to 1-methylcyclopentene catalyzed by Brønsted-acidic carboxylic acid groups on carbon nanotubes as a proof-of-concept probe reaction, we show that catalyst contact with inert, thermally reduced carbon nanotubes leads to order-of-magnitude changes in reaction rate. Furthermore, we demonstrate that these contact-induced effects can also be observed under standard laboratory reaction conditions, where particle-to-particle contact in stirred catalyst powder suspensions is sufficient to demote rates by ∼8-fold. This work provides the foundation for a new method of reaction rate control, which could have implications whenever heterogeneous catalyst particles are in contact with inert materials for liquid-phase reactions in the presence of electrolyte.

Alcohols↗

Machine learning approaches for intentional materials engineering

In this article, the development of nanoporous metals and metallic composites through dealloying processes presents significant opportunities in materials engineering. However, designing multicomponent precursor alloys and establishing corresponding processing methods that yield predictable compositions and nanostructures remain a complex challenge. This article explores how machine learning (ML)-augmented computational and experimental methodologies can tackle these challenges by predicting precursor alloy compositions, final nanoporous structures, and mechanical properties, while integrating ML-enabled autonomous experimentation for material design and quantification. We highlight recent advancements in applying ML to nanostructured materials design via dealloying and discuss how techniques from other nanomaterial designs can be adapted for improved control over morphological and compositional outcomes in nanoporous and nanocomposite materials. Furthermore, we explore the role of ML in autonomous synchrotron x-ray experimentation, enabling real-time feedback between modeling and experimental setups. ML-driven approaches to microstructure characterization and mechanical property prediction are also examined, with a focus on modeling and advanced imaging techniques such as three-dimensional nanotomography. Finally, this article outlines future directions for ML-enhanced materials science, emphasizing the exploration of high-dimensional parameter spaces and the incorporation of materials kinetics into processing and property evaluation, ultimately advancing the design of nanoporous structures and materials science.

36 MATERIALS SCIENCE↗

Machine learning force field model for kinetic Monte Carlo simulations of itinerant Ising magnets

Here, we present a scalable machine learning (ML) framework for large-scale kinetic Monte Carlo (kMC) simulations of itinerant electron Ising systems. As the effective interactions between Ising spins in such itinerant magnets are mediated by conducting electrons, the calculation of energy change due to a local spin update requires solving an electronic structure problem. Such repeated electronic structure calculations could be overwhelmingly prohibitive for large systems. Assuming the locality principle, a convolutional neural network (CNN) model is developed to directly predict the effective local field and the corresponding energy change associated with a given spin update based on Ising configuration in a finite neighborhood. As the kernel size of the CNN is fixed at a constant, the model can be directly scalable to kMC simulations of large lattices. Our approach is reminiscent of the ML force field models widely used in first-principles molecular dynamics simulations. Applying our ML framework to a square-lattice double-exchange Ising model, we uncover unusual coarsening of ferromagnetic domains at low temperatures. Our work highlights the potential of ML methods for large-scale modeling of similar itinerant systems with discrete dynamical variables.

machine learning↗

Development of a deep potential model for F and CF 2 etching of Si and SiO 2

An understanding of plasma-surface interactions at increasingly smaller scales is invaluable for the development of novel technologies and processing techniques. Molecular dynamics (MD) simulations can provide insights into atomic-scale interactions, though they are restricted by the availability of interatomic potentials. Machine learning methods, such as Deep Potential Molecular Dynamics (DeepMD), provide a systematic framework for the development of accurate and flexible ab initio-based models. In this work, we develop DeepMD models for the ion-enhanced etching of Si and SiO 2 by F and CF 2 radicals. We employ an active learning process to expand the data set on which the model is trained and demonstrate its effect on the model accuracy. The DeepMD results are compared to data from classical MD simulations and experiments. Physical sputtering yields of SiO 2 by Ar + ions show good agreement with previous simulation results using conventional interatomic potentials, though the predicted depth profiles are different. Etching yields are calculated as a function of ion energy and neutral to ion flux ratio for the Ar + ion-enhanced etching of SiO 2 and Si by F atoms, as well as for etching of SiO 2 by CF 2 radicals, showing reasonable agreement with experimental data. Finally, an ion-enhanced surface kinetic model is fitted to the DeepMD etch yields, and the fitted parameters are compared to quantities computed directly from DeepMD simulations. This study illustrates how molecular dynamics simulations using machine learning potentials can provide an accurate model of etching processes relevant to device manufacturing.

Kounis-Melas, Andreas [Princeton Univ., NJ (United↗