Abstract: Microdischarges are an important type of plasma discharge that possess several unique characteristics, such as the presence of a stable glow discharge, high plasma density and intense excimer radiation, leading to several potential applications. The intense and controllable gas heating within the extremely small dimensions of microdischarges has been exploited in microthruster technologies by incorporating a micro-nozzle to generate the thrust. This kind of microthruster has a significantly improved specific impulse performance compared to conventional cold gas thrusters, and can meet the requirements arising from the emerging development and application of micro-spacecraft. In this paper, we performed a self-consistent 2D particle-in-cell simulation, with a Monte Carlo collision model, of a microdischarge operating in a prototype micro-plasma thruster with a hollow cylinder geometry and a divergent micro-nozzle. The model takes into account the thermionic electron emission including the Schottky effect, the secondary electron emission due to cathode bombardment by the plasma ions, several different collision processes, and a non-uniform argon background gas density in the cathode-anode gap. Results in the high-pressure (several hundreds of Torr), high-current (mA) operating regime showing the behavior of the plasma density, potential distribution, and energy flux towards the hollow cathode and anode are presented and discussed. In addition, the results of simulations showing the effect of different argon gas pressures, cathode material work function and discharge voltage on the operation of the microdischarge thruster are presented. Our calculated properties are compared with experimental data under similar conditions and qualitative and quantitative agreements are reached.

Abstract: A number of recent computational material design studies based on density functional theory (DFT) calculations have put forward a new class of materials with electrically switchable chemical characteristics that can be exploited in the development of tunable gas storage and electrocatalytic applications. We find systematic flaws in almost every computational study of gas adsorption on polarized or charged surfaces, stemming from an improper and unreproducible treatment of periodicity, leading to very large errors of up to 3 eV in some cases. Two simple corrective procedures that lead to consistent results are proposed, constituting a crucial course correction to the research in the field.

Abstract: In this study, we performed density functional theory (DFT) calculations to investigate different reaction mechanisms of CO oxidation catalyzed by the Si atom embedded defective BC2N nanostructures as well as the analysis of the structural and electronic properties. The structures of all the complexes are optimized and characterized by frequency calculations at the M062X/6-31G* computational level. Also, The electronic structures and thermodynamic parameters of adsorbed CO and O-2 molecules over Si-doped BC2N nanostructures are examined in detail. Moreover, to investigate the curvature effect on the CO oxidation reaction, all the adsorption and CO oxidation reactions on a finite-sized armchair (6,6) Si-BC2NNT are also studied. Our results indicate that there can be two possible pathways for the CO oxidation with O-2 molecule: O-2(g) + CO(g) -> O-2(ads) + CO(ads) -> CO2(g) + O-(ads) and O-(ads) + CO(g) -> CO2(g). The first reaction proceeds via the Langmuir-Hinshelwood (LH) mechanism while the second goes through the Eley-Rideal (ER) mechanism. On the other hand, by increasing the tube diameter, the energy barrier increases due to the strong adsorption energy of the O-2 molecule which is related to its dissociation over the tube surface. Our calculations indicate that the two step energy barrier of the oxidation reaction over Si-BC2NNS is less than that over the Si-BC2NNT. Hence, Si-BC2NNS may serve as an efficient and highly activated substrate to CO oxidation rather than (4,4) Si-BC2NNT. (C) 2018 Elsevier B.V. All rights reserved.

Abstract: A hybrid Monte Carlofluid model is applied to simulate the wafer-temperature-dependent etching of silicon with SF6 inductively coupled plasmas (ICP). The bulk plasma within the ICP reactor volume as well as the surface reactions occurring at the wafer are self-consistently described. The calculated etch rates are validated by experiments. The calculations and experiments are performed at two different wafer temperatures, i.e. 300 and 173 K, resembling conventional etching and cryoetching, respectively. In the case of cryoetching, a physisorbed SFx layer (x = 06) is formed on the wafer, which is negligible at room temperature, because of fast thermal desorption, However, even in the case of cryoetching, this layer can easily be disintegrated by low-energy ions, so it does not affect the etch rates. In the investigated pressure range of 19 Pa, the etch rate is always slightly higher at cryogenic conditions, both in the experiments and in the model, and this could be explained in the model due to a local cooling of the gas above the wafer, making the gas denser and increasing the flux of reactive neutrals, like F and F2, towards the wafer.

Abstract: In this work particle-in-cell/Monte Carlo collision simulations are performed to study the heating mechanism and plasma characteristics in direct current (dc) superposed radio-frequency (RF) capacitively coupled plasmas, operated both in single-frequency (SF) and dual-frequency (DF) regimes. An RF (60/2 MHz) source is applied on the bottom electrode to sustain the discharge, and a dc source is fixed on the top electrode. The heating mechanism appears to be very different in dc superposed SF and DF discharges. When only a single source of 60 MHz is applied, the plasma bulk region is reduced by the dc source, thus the ionization rate and hence the electron density decrease with rising dc voltage. However, when a DF source of 60 and 2 MHz is applied, the electron density can increase upon addition of a dc voltage, depending on the gap length and applied dc voltage. This is explained from the spatiotemporal ionization rates in the DF discharge. In fact, a completely different behavior is observed for the ionization rate in the two half-periods of the LF source. In the first LF half-period, the situation resembles the dc superposed SF discharge, and the reduced plasma bulk region due to the negative dc bias results in a very small effective discharge area and a low ionization rate. On the other hand, in the second half-period, the negative dc bias is to some extent counteracted by the LF voltage, and the sheath close to the dc electrode becomes particularly thin. Consequently, the amplitude of the high-frequency sheath oscillations at the top electrode is largely enhanced, while the LF sheath at the bottom electrode is in its expanding phase and can thus well confine the high-energy electrons. Therefore, the ionization rate increases considerably in this second LF half-period. Furthermore, in addition to the comparison between SF and DF discharges and the effect of gap length and dc voltage, the effect of secondary electrons is examined.

Abstract: In this work, we study the abatement of dilute trichloroethylene (TCE) in air with a negative direct current corona discharge. A numerical model is used to theoretically investigate the underlying plasma chemistry for the removal of TCE, and a reaction pathway for the abatement of TCE is proposed. The Cl atom, mainly produced by dissociation of COCl, is one of the controlling species in the TCE destruction chemistry and contributes to the production of chlorine containing by-products. The effect of humidity on the removal efficiency is studied and a good agreement is found between experiments and the model for both dry (5 % relative humidity (RH)) and humid air (50 % RH). An increase of the relative humidity from 5 % to 50 % has a negative effect on the removal efficiency, decreasing by ±15 % in humid air. The main loss reactions for TCE are with ClO·, O· and CHCl2. Finally, the by-products and energy cost of TCE abatement are discussed.

Abstract: A Monte Carlo model is utilized for studying the behavior of electrons in the afterglow of an analytical microsecond dc pulsed glow discharge. This model uses several quantities as input data, such as electric field and potential, ion flux at the cathode, the fast argon ion and atom impact ionization rates, slow electron density, the electrical characterization of the pulse (voltage and current profiles) and temperature profile. These quantities were obtained by earlier Monte Carlo fluid calculations for a pulsed discharge. Our goal is to study the behavior of the so-called Monte Carlo electrons (i.e., those electrons created at the cathode or by ionization collisions in the plasma which are followed by using the Monte Carlo model) from their origin to the moment when they are absorbed at the cell walls or when they have lost their energy by collisions (being transferred to the group of slow electrons) in the afterglow of the pulsed discharge. The thermalization of the electrons is a phenomenon where the electron-electron Coulomb collisions acquire a special importance. Indeed, in the afterglow the cross sections of the other electron reactions taken into account in the model are very low, because of the very low electron energy. We study the electron energy distributions at several times during and after the pulse and at several positions in the plasma cell, focusing on the thermalization and on the behavior of the electrons in the afterglow. Also, the time evolution of the rates of the various collision processes, the average electron energy, the densities of Monte Carlo and slow electrons and the ionization degree are investigated.

Abstract: Strong terahertz (1THz=1012Hz) radiation can be generated by the electron oscillation in fs-laser-induced wake fields. The interaction of a fs-laser pulse with a low-density plasma layer is studied in detail using numerical simulations. The spatial distribution and temporal evolution of terahertz electron current developed in a low-density plasma layer are presented, which enables us to calculate the intensity distribution of THz radiation. It is shown that laser and plasma parameters, such as laser intensity, pulse width, and background plasma density, are of key importance to the process. The optimum condition for wake-field excitation and terahertz emission is discussed upon the simulation results. Radiation peaked at 6.4 THz, with 900 fs duration and 9% bandwidth, can be generated in a plasma of density 5×1017cm−3. It turns out that the maximum radiation intensity scales as n03a04 when wake field is resonantly excited, where n0 and a0 are, respectively, the plasma density and the normalized field amplitude of the laser pulse.

Abstract: In situ doped epitaxial Si: P films with P concentrations > 1 x 10(21) at./cm(3) are suitable for source-drain stressors of n-FinFETs. These films combine the advantages of high conductivity derived from the high P doping with the creation of tensile strain in the Si channel. It has been suggested that the tensile strain developed in the Si: P films is due to the presence of local Si3P4 clusters, which however do not contribute to the electrical conductivity. During laser annealing, the Si3P4 clusters are expected to disperse resulting in an increased conductivity while the strain reduces slightly. However, the existence of Si3P4 is not proven. Based on first-principles simulations, we demonstrate that the formation of vacancy centered Si3P4 clusters, in the form of four P atoms bonded to a Si vacancy, is thermodynamically favorable at such high P concentrations. We suggest that during post epi-growth annealing, a fraction of the P atoms from these clusters are activated, while the remaining part goes into interstitial sites, thereby reducing strain. We corroborate our conjecture experimentally using positron annihilation spectroscopy, electron spin resonance, and Rutherford backscattering ion channeling studies. (C) 2016 AIP Publishing LLC.

Abstract: We present a 3D and 2D Cartesian quasi-neutral plasma model for a low current argon gliding arc discharge, including strong interactions between the gas flow and arc plasma column.
The 3D model is applied only for a short time of 0.2 ms due to its huge computational cost. It mainly serves to verify the reliability of the 2D model. As the results in 2D compare well with those in 3D, they can be used for a better understanding of the gliding arc basic characteristics. More specifically, we investigate the back-breakdown phenomenon induced by an artificially controlled plasma channel, and we discuss its effect on the gliding arc characteristics. The
back-breakdown phenomenon, or backward-jump motion of the arc, as observed in the experiments, results in a drop of the gas temperature, as well as in a delay of the arc velocity with respect to the gas flow velocity, allowing more gas to pass through the arc, and thus increasing the efficiency of the gliding arc for gas treatment applications.

Abstract: Doped Cr2O3 has been shown to be a p-type transparent conducting oxide (TCO). Its conductivity, however, is low. As for most p-type TCOs, the main problem is the high effective hole mass due to flat valence bands. We use first-principles methods to investigate whether one can increase the valence band dispersion (i.e. reduce the hole mass) by anion alloying with sulfur, while keeping the band gap large enough for transparency. The alloying concentrations considered are given by Cr(4)SxO(6-x), with x = 1-5. To be able to describe the electronic properties of these materials accurately, we first study Cr2O3, examining critically the accuracy of different density functionals and methods, including PBE, PBE+U, HSE06, as well as perturbative approaches within the GW approximation. Our results demonstrate that Cr4S2O4 has an optical band gap of 3.08 eV and an effective hole mass of 1.8 m(e). This suggests Cr4S2O4 as a new p-type TCO host candidate.

Abstract: We perform molecular dynamics simulations to study the flip-flop motion of phosphatidylserine (PS) across the plasma membrane upon increasing oxidation degree of the membrane. Our computational results show that an increase of the oxidation degree in the lipids leads to a decrease of the free energy barrier for translocation of PS through the membrane. In other words, oxidation of the lipids facilitates PS flip-flop motion across the membrane, because in native phospholipid bilayers this is only a “rare event” due to the high energy barriers for the translocation of PS. The present study provides an atomic-scale insight into the mechanisms of the PS flip-flop upon oxidation of lipids, as produced for example by cold atmospheric plasma, in living cells.