Search NASA⌕ Search

SEARCH · Search NASA

Results for “MD simulation”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Machine learning enabled quantification of the hydrogen bonds inside the polyelectrolyte brush layer probed using all-atom molecular dynamics simulations

The configuration of densely grafted charged polyelectrolyte (PE) brushes is strongly dictated by the properties and behavior of the counterions that screen the PE brush charges and the solvent molecules (typically water) that solvate the brush molecules and these screening counterions. Only recently, efforts have been made to study the PE brushes atomistically, thereby shedding light on the properties of brush-supported ions and water molecules. However, even for such efforts, there are limitations associated with using a generic definition to estimate certain properties of water and ions inside the brush layer. For example, water–water hydrogen bonds (HBs) will behave differently for locations outside and inside the brush layer, given the fact that the densely closely grafted PE brush molecules create a soft nanoconfinement where the water connectivity becomes highly disrupted: therefore, using the same definition to quantify the HBs inside and outside the brush layer will be unwise. In this paper, we address this limitation by employing an unsupervised machine learning (ML) approach to predict the water–water hydrogen bonding inside a cationic PE brush layer modeled using all-atom molecular dynamics (MD) simulations. Here, the ML method, which relies on a clustering approach and uses the equilibrium coordinates of the water molecules (obtained from the all-atom MD simulations) as the input, is capable of identifying the structural modification of water–water HBs (revealed through appropriate clustering of the data) inside the PE brush layer induced soft nanoconfinement. Such capabilities would not have been possible by using a generic definition of the HBs. Our calculations lead to four key findings: (1) the clusters formed inside and outside the brush layer are structurally similar; (2) the margin of the cluster is shorter inside the PE brush layer confirming the possible disruption of the HBs inside the PE brush layer; (3) the average “hydrogen–acceptor-oxygen–donor-oxygen” angle that defines the HB is reduced for the HBs formed inside the brush layer; (4) the use of the generic definition (definition usable for characterizing the HBs in brush-free bulk) leads to an overprediction of the number of HBs formed inside the PE brush layer.

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

Tight-binding study of beryllium

In this work, we have applied the Naval Research Laboratory Tight-binding (NRL-TB) Method to the hcp alkaline-earth metal Beryllium. This approach consists of fitting to energy band and total energy data generated by the Linearized Augmented Plane Wave (LAPW) for the fcc, bcc, sc, and hcp structures as a function of volume. First, we found that using the Generalized Gradient Approximation (GGA) for the input data secures a better agreement to experimental volumes than the LDA. Second, we found that including small volumes to the fit is needed to get TB parameters that are transferable enough to successfully predict the energies of structures that were not fitted, and to perform molecular dynamics (MD) simulations. The non-orthogonal Hamiltonian obtained, is successful in correctly predicting the hcp lattice as the ground state of Be and finding eight other structures, that were not fitted to LAPW, positioned at higher energies. In addition, the NRL-TB produces accurate energy bands, densities of states, elastic constants, the Bain path and several quantities derived from MD simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Relating Geometric Nanoconfinement and Local Molecular Environment to Diffusion in Ionic Polymer Membranes

Understanding structure-diffusion relationships in ionic polymer membranes not only illuminates fundamental mechanisms for membrane transport, but also informs further membrane development. In this study, we employ NMR diffusometry and molecular dynamics (MD) simulations to investigate diffusion coefficients and activation energies of diffusion (Ea) for both a lithium-form perfluorosulfonic acid membrane and aqueous solutions of lithium triflate. NMR shows that, at high water/ion molar ratio (?), Ea for water diffusion is slightly lower in the membrane than in solution. Conversely, membrane transport exhibits much higher Ea at low hydration as compared to solutions. MD simulations of a model system consisting of carbon nanotubes with varying diameter reveal that Ea of diffusion clearly relates to both the geometric nanoconfinement of the hydrophilic pathways in the membrane as well as the local molecular environment. These results demonstrate that Ea of diffusion can serve as a revelatory tool for the study of molecular transport processes and the coupling of morphology with transport at the nanoscale. This study thus provides new insights and new experimental and computational models for understanding transport in a wide range of polymer membranes, such as those used in molecular separations applications.

Zhang, Rui↗
DOE OSTI2025

Prediction of carbon nanostructure mechanical properties and the role of defects using machine learning

Graphene-based nanostructures hold immense potential as strong and lightweight materials, however, their mechanical properties such as modulus and strength are difficult to fully exploit due to challenges in atomic-scale engineering. This study presents a database of over 2,000 pristine and defective nanoscale CNT bundles and other graphitic assemblies, inspired by microscopy, with associated stress–strain curves from reactive molecular dynamics (MD) simulations using the reactive INTERFACE force field (IFF-R). These 3D structures, containing up to 80,000 atoms, enable detailed analyses of structure-stiffness-failure relationships. By leveraging the database and physics- and chemistry-informed machine learning (ML), accurate predictions of elastic moduli and tensile strength are demonstrated at speeds 1,000 to 10,000 times faster than efficient MD simulations. Hierarchical Graph Neural Networks with Spatial Information (HS-GNNs) are introduced, which integrate chemistry knowledge. HS-GNNs as well as extreme gradient boosted trees (XGBoost) achieve forecasts of mechanical properties of arbitrary carbon nanostructures with only 3 to 6% mean relative error. The reliability equals experimental accuracy and is up to 20 times higher than other ML methods. Predictions maintain 8 to 18% accuracy for large CNT bundles, CNT junctions, and carbon fiber cross-sections outside the training distribution. The physics- and chemistry-informed HS-GNN works remarkably well for data outside the training range while XGBoost works well with limited training data inside the training range. The carbon nanostructure database is designed for integration with multimodal experimental and simulation data, scalable beyond 100 nm size, and extendable to chemically similar compounds and broader property ranges. The ML approaches have potential for applications in structural materials, nanoelectronics, and carbon-based catalysts.

Winetrout, Jordan J.↗
DOE OSTI2023

Multiscale molecular simulations for the solvation of lignin in ionic liquids

Lignin, the second most abundant biopolymer found in nature, has emerged as a potential source of sustainable fuels, chemicals, and materials. Finding suitable solvents, as well as technologies for efficient and affordable lignin dissolution and depolymerization, are major obstacles in the conversion of lignin to value-added products. Certain ionic liquids (ILs) are capable of dissolving and depolymerizing lignin but designing and developing an effective IL for lignin dissolution remains quite challenging. To address this issue, the COnductor-like Screening MOdel for Real Solvents (COSMO-RS) model was used to screen 5670 ILs by computing logarithmic activity coefficients (ln(γ)) and excess enthalpies (H E ) of lignin, respectively. Based on the COSMO-RS computed thermodynamic properties (ln(γ) and H E ) of lignin, anions such as acetate, methyl carbonate, octanoate, glycinate, alaninate, and lysinate in combination with cations like tetraalkylammonium, tetraalkylphosphonium, and pyridinium are predicted to be suitable solvents for lignin dissolution. The dissolution properties such as interaction energy between anion and cation, viscosity, Hansen solubility parameters, dissociation constants, and Kamlet–Taft parameters of selected ILs were evaluated to assess their propensity for lignin dissolution. Furthermore, molecular dynamics (MD) simulations were performed to understand the structural and dynamic properties of tetrabutylammonium [TBA] + -based ILs and lignin mixtures and to shed light on the mechanisms involved in lignin dissolution. MD simulation results suggested [TBA] + -based ILs have the potential to dissolve lignin because of their higher contact probability and interaction energies with lignin when compared to cholinium lysinate.

09 BIOMASS FUELS↗
DOE OSTI2022

Dual Atomic Coherence in the Self-Assembly of Patchy Heterostructural Nanocrystals

Advances in the synthesis and self-assembly of nanocrystals have enabled researchers to create a plethora of different nanoparticle superlattices. But while many super lattices with complex types of translational order have been realized, rotational order of nanoparticle building blocks within the lattice is more difficult to achieve. Self-assembled superstructures with atomically coherent nanocrystal lattices, which are desirable due to their exceptional electronic and optical properties, have been fabricated only for a few selected systems. Here, we combine experiments with molecular dynamics (MD) simulations to study the self-assembly of heterostructural nanocrystals (HNCs), consisting of a near-spherical quantum dot (QD) host decorated with a small number of epitaxially grown gold nanocrystal (Au NC) "patches ". Self-assembly of these HNCs results in face-centered-cubic (fcc) superlattices with well-defined orientational relationships between the atomic lattices of both QD hosts and Au patches. MD simulations indicate that the observed dual atomic coherence is linked to the number, size, and relative positions of gold patches. Here this study provides a strategy for the design and fabrication of NC superlattices with large structural complexity and delicate orientational order.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

A spall and diffraction study of nanosecond pressure release across the iron ε-α phase boundary

The extreme response of polycrystalline iron at high pressures and high strain rates is revealed by means of high-power laser pulses. The compression portion of the pulse coupled with x-ray diffraction identifies the expected body-centered cubic (α) to hexagonal close packed (ε) displacive transformation. Upon release, observation shows that the complete reverse transformation takes approximately 8 ns and that the structure returns to its initial microstructural configuration, in a reversible transformation path. This is in good agreement with molecular dynamics (MD) simulations which predict an inverse dependence between transformation time and strain rate. The grain size is reduced from μm to nm range during compression and begins increasing back to the original grain size on decompression. The kinetics of the transition is dictated by heterogenous nucleation as it follows the Johnson-Mehl-Avrami-Kolmogorov equation with the appropriate time exponent of ~1. This is confirmed by MD simulations which also identify profuse twinning and dislocation generation. Here, the tensile pulse generated upon reflection at the free surface is captured by time-resolved free surface velocity measurements from which a peak tensile stress of 7 GPa is obtained, in stark contrast with its quasi-static value of ~200 MPa. At these strain rates, the strength of grain interiors, which is determined by twinning and slip exceeds the strength of the boundaries, and failure initiates preferentially in the latter.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Lytic polysaccharide monooxygenase increases cellobiohydrolases activity by promoting decrystallization of cellulose surface

Efficient depolymerization of crystalline cellulose requires cooperation between multiple cellulolytic enzymes. Through biochemical approaches, molecular dynamics (MD) simulation, and single-molecule observations using high-speed atomic force microscopy (HS-AFM), we quantify and track synergistic activity for cellobiohydrolases (CBHs) with a lytic polysaccharide monooxygenase (LPMO) from Phanerochaete chrysosporium. Increasing concentrations of LPMO (AA9D) increased the activity of a glycoside hydrolase family 6 CBH, Cel6A, whereas the activity of a family 7 CBH (Cel7D) was enhanced only at lower concentrations of AA9D. MD simulation suggests that the result of AA9D action to produce chain breaks in crystalline cellulose can oxidatively disturb the crystalline surface by disrupting hydrogen bonds. HS-AFM observations showed that AA9D increased the number of Cel7D molecules moving on the substrate surface and increased the processivity of Cel7D, thereby increasing the depolymerization performance, suggesting that AA9D not only generates chain ends but also amorphizes the crystalline surface, thereby increasing the activity of CBHs.

09 BIOMASS FUELS↗
DOE OSTI2021

Structure of ice confined in silica nanopores

Observed anomalous thermodynamic properties of confined water such as deviations in the melting point and freezing point motivate the determination of the structure of confined water as a function of pore size and temperature. Here, we investigate the dynamic evolution of the structure of confined ice in SBA-15 porous materials with pore diameters of 4 nm, 6 nm, and 8 nm at temperatures ranging from 183 K to 300 K using in operando Wide-Angle X-Ray Scattering (WAXS) measurements, X-Ray Partial Distribution Function (PDF) measurements, and classical Molecular Dynamics (MD) simulations. Formation of hexagonal ice structures is noted in all the three pore sizes. In silica nanopores with diameters of 4 nm, cubic ice formation is noted in addition to hexagonal ice. Longer lasting hydrogen bonds and longer residence times of the water molecules in the first coordination shell contribute to observed crystalline organization of ice in confinement. Self-diffusion coefficients of confined liquid water, predicted from classical MD simulations, are four orders of magnitude higher compared to ice formed in confinement. These experimental and simulation results provide comprehensive insights underlying the organization of confined water and ice in silica nanopores and the underlying physico-chemical interactions that contribute to the observed structures.

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

A Revisited Mechanism of the Graphite-to-Diamond Transition at High Temperature

The graphite-diamond transition, under high-pressure and high-temperature conditions, has been a central subject in physical science. However, its atomistic mechanism remains under debate. Employing large-scale molecular dynamics (MD) simulations, we report a mechanism whereby the diamond nuclei in the graphite matrix propagate in two preferred directions, among which the graphite [120] is about 2.5 times faster than [001]. Consequently, cubic diamond (CD) is the kinetically favorable product, while only a few hexagonal diamonds (HDs) can exist as the twins of CDs. The coherent interface of t-(100)gr//(11-1)cd + [010]gr//[1-10]cd observed in MD simulation was confirmed by our high-resolution transmission electron microscopy experiment. This proposed mechanism not only clarifies the role of HD in graphite-diamond transition but also yields atomistic insight into strengthening synthetic diamond via microstructure engineering.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Predicting thermophysical properties of molten salts in the MgCl 2 -NaCl-KCl-LiCl system with a shell-model potential

Ternary eutectic salts composed of MgCl 2 , NaCl, and KCl, referred to as MNK salts, have recently emerged as promising candidates as high-temperature heat transfer fluids and thermal energy storage media. Here we performed classical molecular dynamics (MD) simulations to predict the densities, specific heat capacities, viscosities, and ionic self-diffusivities for MNK salts over a wide temperature range. The impact of LiCl additive on their thermophysical properties was also investigated. To capture the electronic polarization of Cl anions by neighboring cations, we developed a novel shell-model potential using the force-matching method and a dataset of ab initio calculated interatomic forces. Our extensive MD simulations predict structure and properties for pure salts and binary/ternary salt mixtures in the MgCl 2 -NaCl-KCl-LiCl system in overall good agreement with available experimental and theoretical data, which corroborates the accuracy and reliability of our developed potential.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Multiscale investigation of the mechanism of biomass deconstruction in the dimethyl isosorbide/water Co-solvent pretreatment system

In the context of promoting a circular bioeconomy, the development of green and efficient lignocellulosic biomass pretreatment technologies so as to realize high value-added biomass utilization is of intense interest. We demonstrated the potential of the bio-based green solvent dimethyl isosorbide (DMI) for the fractionation of Eucalyptus biomass with excellent performance. Here, to investigate the mechanisms involved in biomass fractionation, microimaging and microspectroscopic techniques were employed together with molecular dynamics (MD) simulation and COSMO-RS quantum chemical calculations to derive multiscale information. Both the microstructure and regional chemistry of the cell wall vary significantly with the volume ratio of DMI/H 2 O. The strongest effects were found at DMI/H 2 O = 9:1 and showed visible cell wall tearing cracks and cell wall deformation and collapse as well as the lowest values of cell wall thickness and circularity. From the MD simulations, lignin exhibits collapsed-like structure in pure H 2 O with low solvent accessibility surface area (SASA) and radius of gyration (R g ). In contrast, lignin in DMI/H 2 O shows extended structure with high SASA and solvent interactions dominated by van der Waals forces, with maximal contact in the 9:1 (v/v) system. Further, the COSMO-RS calculated sigma (σ-) potential suggests the intermolecular interactions in DMI and DMI/H 2 O co-solvent are weak, leading to stronger interaction with lignin and correspondingly higher lignin dissolution. The radial distribution functions and σ-potential all show that again DMI/H 2 O at 9:1 is an optimal volume ratio for high lignin dissolution. Finally, this study provides a solvent-ratio dependent mechanism for the action of polar aprotic solvents in the deconstruction of biomass.

09 BIOMASS FUELS↗
DOE OSTI2024

Accelerated Sequence Design of Star Block Copolymers: An Unbiased Exploration Strategy via Fusion of Molecular Dynamics Simulations and Machine Learning

Star block copolymers (s-BCPs) have potential applications as novel surfactants or amphiphiles for emulsification, compatibilization, chemical transformations, and separations. s-BCPs have chain architectures where three or more linear diblock copolymer arms comprised of two chemically distinct linear polymers, e.g., solvophobic and solvophilic chains, are covalently joined at one point. The chemical composition of each of the subunit polymer chains comprising the arms, their molecular weights, and the number of arms can be varied to tailor the surface and interfacial activity of these architecturally unique molecules. Further, this makes identification of the optimal s-BCP design nontrivial as the total number of plausible s-BCP architectures is experimentally or computationally intractable. In this work, we use molecular dynamics (MD) simulations coupled with a reinforcement learning-based Monte Carlo tree search (MCTS) to identify s-BCP designs that minimize the interfacial tension between polar and nonpolar solvents. We first validate the MCTS approach for the design of small- and medium-sized s-BCPs and then use it to efficiently identify sequences of copolymer blocks for large-sized s-BCPs. The structural origins of interfacial tension in these systems are also identified by using the configurations obtained from MD simulations. Chemical insights into the arrangement of copolymer blocks that promote lower interfacial tension were mined using machine learning (ML) techniques. Overall, this work provides an efficient approach to solve design problems via fusion of simulations and ML and provides important groundwork for future experimental investigation of s-BCPs for various applications.

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

Elucidating the Atomic Structures of the Gel Layer Formed during Aluminoborosilicate Glass Dissolution: An Integrated Experimental and Simulation Study

The altered glasses produced during aqueous dissolution of silicate and borosilicate glasses are among the most complex structures to understand at the atomic level due to their amorphous nature, random porosity and various levels of hydration. In this study, we gained insights of the complex atomic structures of altered aluminoborosilicate glasses by combining a range of experimental and computational approaches. The altered glasses were prepared by dissolution of three glasses with varying level of alumina in acid for 7 days. A comprehensive set of experimental [elemental analysis, high-energy X-ray diffraction, 29 Si and 27 Al solid-state nuclear magnetic resonance (NMR), and O 1s X-ray photoelectron spectroscopy (XPS)] and modeling (molecular dynamics (MD) simulations using non-reactive and reactive force fields) were used to study the atomic structures of these altered glasses. Elemental analysis showed that most of the B in the pristine glasses was leached into the solution and was not contained in the altered glass. The 29 Si and 27 Al solid-state NMR spectra revealed that the altered glasses have more polymerized silicate network as compared to those in the pristine glasses due to the reformation of linkages among Si and Al oxygen polyhedral in the altered glasses. The bridging and non-bridging (or hydroxyl O) atoms in the altered glasses were also quantified from their O 1s XPS spectra. Atomic structure models of the altered glasses were constructed by using MD simulations using the reactive force field based on the compositional information obtained from experiments. Various pore structures were generated using the charge-scaling (CS) method by using different initial densities and CS temperatures; the best CS parameters of each alternated glass were then determined by comparing with the experimental structure factors obtained from high energy X-ray diffraction. Pore and the atomic structures and vibrational properties around these pore surfaces were analyzed. Furthermore, these results from this comprehensive study thus provides a realistic insight into the pore morphology, atomic structure, and vibrational properties of altered glasses.

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

Enhanced Ion Transport and Molecular Packing Stability in Asymmetric 2D Nanostructured π‐Conjugated Thieno[3,2‐b]Thiophene‐Based Liquid Crystal

Organic semiconductors based on liquid crystal (LC) molecules have attracted increasing interest. Here, in this work, two linear LCs based on 2,5‐bis(thien‐2‐yl)thieno[3,2‐b]thiophene (BTTT) mesogen are designed and synthesized, including BTTT/dEO3 with two symmetrically attached tri(ethylene oxide) groups and BTTT/mEO6 with one asymmetrically attached hexa(ethylene oxide) group. These two molecules have comparable functional‐group compositions but different molecular geometries, leading to their moderately different material performances. Both LCs show smectic mesophases with relatively low transition temperatures as confirmed by differential scanning calorimetry and polarized optical microscopy. A combination of experimental grazing incidence wide‐angle X‐ray scattering and molecular dynamics (MD) simulations reveals a herringbone packing motif of BTTT segments in both LCs while a smaller molecular tilt angle in BTTT/mEO6. Ionic conductivities are measured by doping LCs with different amounts of ionic dopants, lithium bis(trifluoromethanesulfonyl)imide (LiTFSI). BTTT/mEO6 shows better smectic phase stability to higher LiTFSI doping ratios. Both LCs exhibit similar ionic conductivities in the smectic phases, but BTTT/mEO6 outperforms BTTT/dEO3 by a factor of three in the amorphous phase at higher temperatures. MD simulations, performed to examine the ion solvation environment, reveal that BTTT/mEO6 is more efficient in coordinating Li‐ions and screening their interactions with TFSI‐ions which further promote ionic transport.

grazing incidence wide-angle X-ray scattering↗
DOE OSTI2020

Effect of Hydrophilicity and Interfacial Water Structure on Particle Attachment

Oriented attachment (OA) of nanoparticles is an important and common crystallization pathway. Efforts to isolate and quantify the interparticle forces that underlie OA have only recently been initiated and the role played by the interfacial solvent structure in the confined region between nanocrystal surfaces in enabling and controlling OA remains an open question. In this work, classical molecular dynamics (MD) simulations were performed to compute the free energy landscape of nanoparticles approaching mineral substrates in both vacuum and pure water. Periclase (MgO), muscovite (KAl 3 Si 3 O 10 (OH) 2 ), and pyrophyllite (Al 2 Si 4 O 10 (OH) 2 ) nanoparticles exposing the (100), (001), and (001) surfaces, respectively, were brought to contact with substrates exposing the same surfaces. These mineral surfaces were selected to investigate a range of hydrophilicities (periclase > muscovite > pyrophyllite). The MD simulations revealed that the interfacial water structure created free energy minima for nanoparticle approach that corresponded to integer numbers of intervening water layers. The depth of the free energy minima and the free energy barriers between minima correlated positively with the hydrophilicity of the mineral surface. While the range of solvent-mediated interactions was predicted to be short (a few intervening water layers at most for MgO), the free energy barriers were calculated to be large (fraction of an eV per nm 2 ) and the effect of misalignment diminished rapidly with the number of intervening water layers. The results help explain the formation of mesocrystals and the observed ability of nanoparticles to rotate before attachment while in a long-lived solvent-separated state.

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

Rational Design of High-Performing Electrodes in Energy Storage Devices

Our approach in this project is to control conductivities of different ions by modulating the pore sizes in MOFs as sulfur hosts in cathodes. To validate our hypothesis that controlling the pore size of MOFs can effectively prevent LiPSs shuttling while allowing sufficient Li ion transport, which will ensure electrochemical cycling performance, we performed molecular dynamics (MD) simulations using a series of UiO MOFs with different linker sizes including LiPSs and electrolyte components (1M LiTFSI in DME/DOL 1:1 volume) in pores. Specifically, UiO-66, UiO-67, and UiO-68 are Zr-based MOFs with a linker length of 5.67 Å, 10.07 Å, and 14.35 Å, respectively, leading to different pore diameters of 8.4 Å, 12.8 Å, and 17.8 Å. We calculated the mean-squared displacements of Li and LiPS from MD simulations at 300 K for 50 ns or over to assess their mobility in MOFs.

25 ENERGY STORAGE↗
DOE OSTI2024

A transient site balance model for atomic layer etching

We present a transient site balance model of plasma-assisted atomic layer etching of silicon (Si) with alternating exposure to chlorine gas (Cl 2 ) and argon ions (Ar + ). Molecular dynamics (MD) simulation results are used to provide parameters for the model. The model couples the dynamics of a top monolayer surface region ('top layer') and a perfectly mixed subsurface region ('mixed layer'). The differential equations describing the rates of change of the Cl coverage in the two layers are transient mass balances. Model predictions include Cl coverages and rates of etching of various species from the surface as a function of Cl 2 or Ar + fluence. The simplified phenomenological model reproduces the MD simulation results well over a range of conditions. Comparing model predictions directly to experimental optical emission spectroscopy data, as reported in a previous paper, provides further evidence of the accuracy of the model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗