Abstract: A chemical kinetics model is developed for a CO2/N2 microwave plasma, focusing especially on the vibrational levels of both CO2 and N2. The model is used to calculate the CO2 and N2 conversion as well as the energy efficiency of CO2 conversion for different power densities and for N2 fractions in the CO2/N2 gas mixture ranging from 0 to 90%. The calculation results are compared with measurements, and agreements within 23% and 33% are generally found for the CO2 conversion and N2 conversion, respectively. To explain the observed trends, the destruction and formation processes of both CO2 and N2 are analyzed, as well as the vibrational distribution functions of both CO2 and N2. The results indicate that N2 contributes in populating the lower asymmetric levels of CO2, leading to a higher absolute CO2 conversion upon increasing N2 fraction. However, the effective CO2 conversion drops because there is less CO2 initially present in the gas mixture; thus, the energy efficiency also drops with rising N2 fraction.

Abstract: In this paper, the negative ion behavior in a C4F8 inductively coupled plasma (ICP) is investigated using a hybrid model. The model predicts a non-monotonic variation of the total negative ion density with power at low pressure (1030 mTorr), and this trend agrees well with experiments that were carried out in many fluorocarbon (fc) ICP sources, like C2F6, CHF3, and C4F8. This behavior is explained by the availability of feedstock C4F8 gas as a source of the negative ions, as well as by the presence of low energy electrons due to vibrational excitation at low power. The maximum of the negative ion density shifts to low power values upon decreasing pressure, because of the more pronounced depletion of C4F8 molecules, and at high pressure (∼50 mTorr), the anion density continuously increases with power, which is similar to fc CCP sources. Furthermore, the negative ion composition is identified in this paper. Our work demonstrates that for a clear understanding of the negative ion behavior in radio frequency C4F8 plasma sources, one needs to take into account many factors, like the attachment characteristics, the anion composition, the spatial profiles, and the reactor configuration. Finally, a detailed comparison of our simulation results with experiments is conducted.

Abstract: We perform first-principles calculations to investigate the preferred positions of B and P dopants, both neutral and in their preferred charge state, in Si and Si/SiO2 core-shell nanowires (NWs). In order to understand the observed trends in the formation energy, we isolate the different effects that determine these formation energies. By making the distinction between the unrelaxed and the relaxed formation energy, we separate the impact of the relaxation from that of the chemical environment. The unrelaxed formation energies are determined by three effects: (i) the effect of strain caused by size mismatch between the dopant and the host atoms, (ii) the local position of the band edges, and (iii) a screening effect. In the case of the SiNW (Si/SiO2 NW), these effects result in an increase of the formation energy away from the center (interface). The effect of relaxation depends on the relative size mismatch between the dopant and host atoms. A large size mismatch causes substantial relaxation that reduces the formation energy considerably, with the relaxation being more pronounced towards the edge of the wires. These effects explain the surface segregation of the B dopants in a SiNW, since the atomic relaxation induces a continuous drop of the formation energy towards the edge. However, for the P dopants, the formation energy starts to rise when moving from the center but drops to a minimum just next to the surface, indicating a different type of behavior. It also explains that the preferential location for B dopants in Si/SiO2 core-shell NWs is inside the oxide shell just next to the interface, whereas the P dopants prefer the positions next to the interface inside the Si core, which is in agreement with recent experiments. These preferred locations have an important impact on the electronic properties of these core-shell NWs. Our simulations indicate the possibility of hole gas formation when B segregates into the oxide shell.

Abstract: Because of its excellent properties in nanotechnology applications, atomic layer deposition of ruthenium (Ru) has been the subject of numerous experimental studies. Recently, two different Ru precursors were compared for plasma-enhanced atomic layer deposition (PEALD) of Ru, and their reactivity was found to be different. Inhibition was observed for bis(ethylcyclopentadienyl)ruthenium (Ru(EtCp)(2)), while nearly linear growth behavior was observed for (methylcyclopentadienyl-pyrrolyl)ruthenium (Ru(MeCp)Py). To understand this difference in reactivity, we investigate the adsorption of RuCp, and RuCpPy (i.e., without substituents) on a TiN surface using calculations based on periodic boundary conditions density functional theory (DFT) combined with experiments based on Rutherford backscattering spectroscopy (RBS). The calculations demonstrate that the RuCpPy precursor chemisorbs on the TiN(100) surface while the RuCp2 precursor only physisorbs. We propose a reaction mechanism for the chemisorption of RuCpPy. The area density of the calculated RuCpPy surface species is compared with the experimental values from RBS. The impact of a H-plasma is also investigated. The DFT calculations and experimental results from RBS provide insight into the adsorption processes of the RuCpPy and RuCp2 precursors on the TiN(100) surface.

Abstract: Continued miniaturization of metal-oxide-semiconductor field-effect transistors (MOSFETs) requires an ever-decreasing thickness of the gate oxide. The structure of ultrathin silicon oxide films, however, critically depends on the oxidation mechanism. Using reactive atomistic simulations, we here demonstrate how the oxidation mechanism in hyperthermal oxidation of such structures may be controlled by the oxidation temperature and the oxidant energy. Specifically, we study the interaction of hyperthermal oxygen with energies of 15 eV with thin SiOx (x ≤ 2) films with a native oxide thickness of about 10 Å. We analyze the oxygen penetration depth probability and compare with results of the hyperthermal oxidation of a bare Si(100){2 × 1} (c-Si) surface. The temperature-dependent oxidation mechanisms are discussed in detail. Our results demonstrate that, at low (i.e., room) temperature, the penetrated oxygen mostly resides in the oxide region rather than at the SiOx|c-Si interface. However, at higher temperatures, starting at around 700 K, oxygen atoms are found to penetrate and to diffuse through the oxide layer followed by reaction at the c-Si boundary. We demonstrate that hyperthermal oxidation resembles thermal oxidation, which can be described by the DealGrove model at high temperatures. Furthermore, defect creation mechanisms that occur during the oxidation process are also analyzed. This study is useful for the fabrication of ultrathin silicon oxide gate oxides for metal-oxide-semiconductor devices as it links parameters that can be straightforwardly controlled in experiment (oxygen temperature, velocity) with the silicon oxide structure.

Abstract: Nonthermal atmospheric pressure plasmas are gaining increasing attention for biomedical applications. However, very little fundamental information on the interaction mechanisms between the plasma species and biological cells is currently available. We investigate the interaction of important plasma species, such as OH, H2O2, O, O3, as well as O2 and H2O, with bacterial peptidoglycan by means of reactive molecular dynamics simulations, aiming for a better understanding of plasma disinfection. Our results show that OH, O, O3, and H2O2 can break structurally important bonds of peptidoglycan (i.e., CO, CN, or CC bonds), which consequently leads to the destruction of the bacterial cell wall. The mechanisms behind these breakups are, however, dependent on the impinging plasma species, and this also determines the effectiveness of the cell wall destruction.

Abstract: A database-driven approach combined with ab initio density functional theory (DFT) simulations is used to identify and simulate alternative ferroelectric materials beyond Hf(Zr)O-2. The database-driven screening method identifies a class of wurtzite ferroelectric materials. DFT simulations of wurtzite magnesium chalcogenides, including MgS, MgSe, and MgTe, show their potential to achieve improved ferroelectric (FE) stability, simple atomistic unit cell structure, and large FE polarization. Strain engineering can effectively modulate the FE switching barrier height for facilitating FE switching. The effect of the piezoelectric property on the FE switching barrier heights is also examined.

Abstract: In this paper, the splitting of CO2 in a pulsed plasma system, such as a dielectric barrier discharge (DBD), is evaluated from a chemical point of view by means of numerical modeling. For this purpose, a chemical reaction set of CO2 in an atmospheric pressure plasma is developed, including the vibrational states of CO2, O2, and CO. The simulated pulses are matched to the conditions of a filament (or microdischarge) and repeated with intervals of 1 μs. The influence of vibrationally excited CO2 as well as other neutral species, ions, and electrons on the CO2 splitting is discussed. Our calculations predict that the electrons have the largest contribution to the CO2 splitting at the conditions under study, by electron impact dissociation. The contribution of vibrationally excited CO2 levels in the splitting of CO2 is found be 6.4%, when only considering one microdischarge pulse and its afterglow, but it can be much higher for consecutive discharge pulses, as is typical for a filamentary DBD, when the interpulse time is short enough and accumulation effects in the vibrationally excited CO2 densities can occur.

Abstract: The adsorption probability and reaction behavior of CHx plasma species on various nickel catalyst surfaces is investigated by means of reactive molecular dynamics (MD) simulations using the ReaxFF potential. Such catalysts are used in the reforming of hydrocarbons and in the growth of carbon nanotubes, and further insight in the underlying mechanisms of these processes is needed to increase their applicability. Single and consecutive impacts of CHx radicals (x={1,2,3}) were performed on four different Ni surfaces, at a temperature of 400 K. The adsorption probability is shown to be related to the number of free electrons, i.e. a higher number leads to more adsorptions, and the steric hindrance caused by the hydrogen atoms bonded to the impacting CHx species. Furthermore, some of the CH bonds break after adsorption, which generally leads to diffusion of the hydrogen atom over the surface. Additionally, these adsorbed H-atoms can be used in reactions to form new molecules, such as CH4 and C2Hx, although this is dependent on the precise morphology of the surface. New molecules are also formed by subtraction of H-atoms from adsorbed radicals, leading to occasional formation of H2 and C2Hx molecules.

Abstract: We apply an implicit particle-in-cell Monte-Carlo (PIC-MC) method to study a radio-frequency argon microdischarge at steady state in the glow discharge limit, in which the microdischarge is sustained by secondary electron emission from the electrodes. The plasma density, electron energy distribution function (EEDF), and electron temperature are calculated in a wide range of operating conditions, including driving voltage, microdischarge gap, and pressure. Also, the effect of gap size scaling (in the range of 50-1000 μm) on the plasma sustaining voltage and peak electron density at atmospheric pressure is examined, which has not been explored before. In our simulations, three different EEDFs, i.e., a so-called three temperature hybrid mode, a two temperature α mode, and a two temperature γ mode distribution, are identified at different gaps and voltages. The maximum sustaining voltage to avoid a transition from the glow mode to an arc is predicted, as well as the minimum sustaining voltage for a steady glow discharge. Our calculations elucidate that secondary electrons play an essential role in sustaining the discharge, and as a result the relationship between breakdown voltage and gap spacing is far away from the Paschen law at atmospheric pressure.

Abstract: Using reactive molecular dynamics simulations, we have investigated the effect of single-impact, low-energy (thermal-100 eV) bombardment of a Si(100){2 × 1} surface by atomic and molecular oxygen. Penetration probability distributions, as well as defect formation distributions, are presented as a function of the impact energy for both species. It is found that at low impact energy, defects are created chemically due to the chemisorption process in the top layers of the surface, while at high impact energy, additional defects are created by a knock-on displacement of Si. These results are of particular importance for understanding device performances of silica-based microelectronic and photovoltaic devices.

Abstract: Particle-in-cell/Monte Carlo collision simulations, coupled with an external circuit, are used to investigate the behavior of pulsed dual-frequency (DF) capacitively coupled plasmas (CCPs). It is found that the phase shift θ between the high (or low) frequency source and the pulse modulation has a great influence on the ion density and the ionization rate. By pulsing the high frequency source, the time-averaged ion density shows a maximum when θ = 90∘. The time-averaged ion energy distribution functions (IEDFs) at the driven electrode, however, keep almost unchanged, illustrating the potential of pulsed DF-CCP for independent control of ion density (and flux) and ion energy. A detailed investigation of the temporal evolution of the plasma characteristics indicates that several high frequency harmonics can be excited at the initial stage of a pulse period by tuning the phase shift θ, and this gives rise to strong sheath oscillations, and therefore high ionization rates. For comparison, the pulsing of the low frequency source is also studied. In this case, the ion density changes slightly as a function of time, and the time-averaged ion density shows the same trend as in the HF modulation for different phase shifts θ. Moreover, the time-averaged IEDFs at the driven electrode can be modulated, showing the potential to reduce the maximum ion bombardment energy.

Abstract: This contribution presents the effects of a mild O2 plasma treatment on the structural, optical, and electrical properties of single-layer (SLG) and bilayer graphene (BLG). Unexpectedly, we observe only photoluminescence in the SLG parts of a graphene flake composed of regions of various thickness upon O2 plasma treatment, whereas the BLG and few-layer graphene (FLG) parts remain optically unchanged. Confirmed with X-ray photoelectron spectroscopy (XPS) that O2 plasma induces epoxide and hydroxyl-like groups in graphene, density functional theory (DFT) calculations are carried out on representative epoxidized and hydroxylated SLG and BLG models to predict density of states (DOS) and band structures. Sufficiently oxidized SLG shows a bandgap and thus loss of semimetallic behavior, while oxidized BLG maintains its semimetallic behavior even at high oxygen density in agreement with the results of the photoluminescence spectroscopy (PL) experiments. DFT calculations confirm that the Fermi velocity in epoxidized BLG is remarkably comparable with that of pristine SLG, pointing to a similarity of electronic band structure. The similarity is also experimentally demonstrated by the electrical characterization of a plasma-treated BLG-FET. As expected from the electronegative oxygen adatoms in the graphene, epoxidized BLG presents conductive features typical of hole doping. Moreover, the electrical characteristics suggest band structures closely related to that of epoxidized graphene while deviating from that of hydroxylated graphene. Finally, upon O2 plasma treatment of BLG, the four-component 2D peak around 2700 cm1 in the Raman spectrum evolves into a single Lorentzian line, very like the 2D peak of pristine SLG. Summarizing, the data in this contribution recommend that a controlled O2 plasma treatment, which is compatible with CMOS process flow in contrast to wet chemical oxidation methods, provides an efficient and valuable technique to exploit the transport properties of the bottom layer of BLG.

Abstract: The phase transition processes induced by ultrashort, 100 fs pulsed laser irradiation of Au, Cu, and Ni are studied by means of a combined atomistic-continuum approach. A moderately low absorbed laser fluence range, from 200 to 600 J/m2 is considered to study phase transitions by means of a local and a nonlocal order parameter. At low laser fluences, the occurrence of layer-by-layer evaporation has been observed, which suggests a direct solid to vapor transition. The calculated amount of molten material remains very limited under the conditions studied, especially for Ni. Therefore, our results show that a kinetic equation that describes a direct solid to vapor transition might be the best approach to model laser-induced phase transitions by continuum models. Furthermore, the results provide more insight into the applicability of analytical superheating theories that were implemented in continuum models and help the understanding of nonequilibrium phase transitions.

Abstract: Cold atmospheric pressure plasma jets (plasma) operated in ambient air provide a rich source of reactive oxygen and nitrogen species (RONS), which are known to influence biological processes important in disease. In the plasma treatment of diseased tissue such as subcutaneous cancer tumors, plasma RONS need to first traverse an interface between the plasma-skin surface and second be transported to millimeter depths in order to reach deep-seated diseased cells. However, the mechanisms in the plasma generation of RONS within soft tissues are not understood. In this study, we track the plasma jet delivery of RONS into a tissue model target and we delineate two processes: through target delivery of RONS generated (primarily) in the plasma jet and in situ RONS generation by UV photolysis within the target. We demonstrate that UV photolysis promotes the rapid generation of RONS in the tissue model target’s surface after which the RONS are transported to millimeter depths via a slower molecular process. Our results imply that the flux of UV photons from plasma jets is important for delivering RONS through seemingly impenetrable barriers such as skin. The findings have implications not only in treatments of living tissues but also in the functionalization of soft hydrated biomaterials such as hydrogels and extracellular matrix derived tissue scaffolds.