Abstract: Molecular dynamics simulations using the Brenner potential have been performed to investigate reaction mechanisms of various hydrocarbon radicals with low kinetic energies on amorphous hydrogenated carbon (a-C:H) surfaces and to simulate thin a-C:H film growth. Experimental data from an expanding thermal plasma setup were used as input for the simulations. The hydrocarbon reaction mechanisms were studied both during growth of the films and on a set of surface sites specific for a-C:H surfaces. Thin film growth was studied using experimentally detected growth species. It is found that the reaction mechanisms and sticking coefficients are dependent on the specific surface sites, and the structural properties of the growth radicals. Furthermore, it is found that thin a-C:H films can be densified using an additional H-flux towards the substrate.

Abstract: In this theoretical study, the electron emission from a copper surface under ultrashort pulsed laser irradiation is investigated using a one-dimensional particle in cell model. Thermionic emission as well as multiphoton photoelectron emission were taken into account. The emitted electrons create a negative space charge above the target; consequently the generated electric field reduces the electron emission by several orders of magnitude. The simulations indicate that the space charge effect should be considered when investigating electron emission related phenomena in materials under ultrashort pulsed laser irradiation of metals.

Abstract: The plasma behavior in a parallel-plate dielectric barrier discharge (DBD) is simulated by a two-dimensional particle-in-cell/Monte Carlo collision model, comparing for the first time an unpacked (empty) DBD with a packed bed DBD, i.e., a DBD filled with dielectric spheres in the gas gap. The calculations are performed in air, at atmospheric pressure. The discharge is powered by a pulse with a voltage amplitude of −20 kV. When comparing the packed and unpacked DBD reactors with the same dielectric barriers, it is clear that the presence of the dielectric packing leads to a transition in discharge behavior from a combination of negative streamers and unlimited surface streamers on the bottom dielectric surface to a combination of predominant positive streamers and limited surface discharges on the dielectric surfaces of the beads and plates. Furthermore, in the packed bed DBD, the electric field is locally enhanced inside the dielectric material, near the contact points between the beads and the plates, and therefore also in the plasma between the packing beads and between a bead and the dielectric wall, leading to values of $4\times {10}

Abstract: Packed bed plasma reactors (PBPRs) are gaining increasing interest for use in environmental applications, such as greenhouse gas conversion into value-added chemicals or renewable fuels and volatile pollutant removal (e.g. NOx, VOC, K), as they enhance the conversion and energy efficiency of the process compared to a non-packed reactor. However, the plasma behaviour in a PBPR is not well understood. In this paper we demonstrate, by means of a fluid model, that the discharge behaviour changes considerably when changing the size of the packing beads and their dielectric constant, while keeping the interelectrode spacing constant. At low dielectric constant, the plasma is spread out over the full discharge gap, showing significant density in the voids as well as in the connecting void channels. The electric current profile shows a strong peak during each half cycle. When the dielectric constant increases, the plasma becomes localised in the voids, with a current profile consisting of many smaller peaks during each half cycle. For large bead sizes, the shift from full gap discharge to localised discharges takes place at a higher dielectric constant than for smaller beads. Furthermore, smaller beads or beads with a lower dielectric constant require a higher breakdown voltage to cause plasma formation.

Abstract: Computer simulations have been performed to describe the effect of small admixtures of hydrogen to an argon glow discharge in the Grimm-type configuration. The two-dimensional density profiles of the various plasma species (i.e., electrons, Ar+, ArH+, H+, H2+ and H3+ ions, H atoms and H2 molecules, Ar metastable atoms and sputtered Cu atoms) are presented for 1% H2 added to the argon glow discharge, and the effect of different H2 additions (varying between 0.1 and 10%) on the species densities, the hydrogen dissociation degree, and the sputtering process, are investigated. Finally, the relative contributions of various production and loss processes for the different plasma species are calculated.

Abstract: In this paper, a simulation method is described to predict the etching behaviour of Cl2/O2/Ar inductively coupled plasmas on a Si substrate, as used in shallow trench isolation for the production of electronic devices. The hybrid plasma equipment model (HPEM) developed by Kushner et al is applied to calculate the plasma characteristics in the reactor chamber and two additional Monte Carlo simulations are performed to predict the fluxes, angles and energy of the plasma species bombarding the Si substrate, as well as the resulting surface processes such as etching and deposition. The simulations are performed for a wide variety of operating conditions such as gas composition, chamber pressure, power deposition and substrate bias. It is predicted by the simulations that when the fraction of oxygen in the gas mixture is too high, the oxidation of the Si substrate is superior to the etching of Si by chlorine species, resulting in an etch rate close to zero as is also observed in the experiments.

Abstract: The physical processes in an Ar/N2 magnetron discharge used for the reactive sputter deposition of TiNx thin films were simulated with a 2d3v particle-in-cell/Monte Carlo collisions (PIC/MCC) model. Cathode currents and voltages were calculated self-consistently and compared with experiments. Also, ion fractions were calculated and validated with mass spectrometric measurements. With this PIC/MCC model, the influence of N2/Ar gas ratio on the particle densities and fluxes was investigated, taking into account the effect of the poisoned state of the target.

Abstract: We study the evolution of the structural and electronic properties of crystalline indium gallium zinc oxide (IGZO) upon amorphization by first-principles calculation. The bottom of the conduction band (BCB) is found to be constituted of a pseudo-band of molecular orbitals that resonate at the same energy on different atomic sites. They display a bonding character between the s orbitals of the metal sites and an anti-bonding character arising from the interaction between the oxygen and metal s orbitals. The energy level of the BCB shifts upon breaking of the crystal symmetry during the amorphization process, which may be attributed to the reduction of the coordination of the cationic centers. The top of the valence band (TVB) is constructed from anti-bonding oxygen p orbitals. In the amorphous state, they have random orientation, in contrast to the crystalline state. This results in the appearance of localized tail states in the forbidden gap above the TVB. Zinc is found to play a predominant role in the generation of these tail states, while gallium hinders their formation. Last, we study the dependence of the fundamental gap and effective mass of IGZO on mechanical strain. The variation of the gap under strain arises from the enhancement of the anti-bonding interaction in the BCB due to the modification of the length of the oxygen-metal bonds and/or to a variation of the cation coordination. This effect is less pronounced for the amorphous material compared to the crystalline material, making amorphous IGZO a semiconductor of choice for flexible electronics. Finally, the effective mass is found to increase upon strain, in contrast to regular materials. This counterintuitive variation is due to the reduction of the electrostatic shielding of the cationic centers by oxygen, leading to an increase of the overlaps between the metal orbitals at the origin of the delocalization of the BCB. For the range of strain typically met in flexible electronics, the induced variation in the effective mass is found to be negligible (less than 1%).

Abstract: A packed bed dielectric barrier discharge (DBD) was studied by means of fluid modelling, to investigate the influence of the dielectric constant of the packing on the plasma characteristics, for two different gap sizes. The electric field strength and electron temperature are much more enhanced in a microgap reactor than
in a mm-gap reactor, leading to more current peaks per half-cycle, but also to non-quasineutral plasma. Increasing the dielectric constant enhances the electric field further, but only up to a certain value of dielectric constant, being 9 for a microgap and 100 for a mm-gap reactor. The enhanced electric field results in a higher electron temperature, but also lower electron density. This last one strongly affects the reaction rate.

Abstract: The vast biomedical potential of cold atmospheric pressure plasmas (CAPs) is governed by the formation of reactive species. These biologically active species are formed upon the interaction of CAPs with the surroundings. In biological milieu, water plays an essential role. The development of biomedical CAPs thus requires understanding of the sources of the reactive species in aqueous media exposed to the plasma. This is especially important in case of the COST RF plasma jet, which is developed as a reference microplasma system. In this work, we investigated the formation of the OH radicals, H atoms and H2O2 in aqueous solutions exposed to the COST plasma jet. This was done by combining experimental and modelling approaches. The liquid phase species were analysed using UV-Vis spectroscopy and spin trapping with hydrogen isotopes and electron paramagnetic resonance (EPR) spectroscopy. The discrimination between the species formed from the liquid phase and the gas phase molecules was performed by EPR and 1H-NMR analyses of the liquid samples. The concentrations of the reactive species in the gas phase plasma were obtained using a zero-dimensional (0D) chemical kinetics computational model. A three-dimensional (3D) fluid dynamics model was developed to provide information on the induced humidity in the plasma effluent. The comparison of the experimentally obtained trends for the formation of the species as a function of the feed gas and effluent humidity with the modelling results suggest that all reactive species detected in our system are mostly formed in the gas phase plasma inside the COST jet, with minor amounts arising from the plasma effluent humidity.

Abstract: A numerical model for NxOy remediation in humid air plasma produced with a dielectric barrier discharge at atmospheric pressure is presented. Special emphasis is given to NO2 and N2O reduction with the decrease of O2 content in the feedstock gas. A detailed reaction mechanism including electronic and ionic processes, as well as the contribution of radicals and excited atomic/molecular species is proposed. The temporal evolution of the densities of NO, NO2 and N2O species, and some other by-products, is analyzed, and the major pathways for the NxOy remediation are discussed for one pulse. Subsequently, simulations are presented for a multi-pulses case, where three O2 contents are tested for optimization of the remediation process. It is found that when the gas mixture O2/N2/H2O/NOx has no initial O2 content, the best NOx and N2O remediation is achieved.

Abstract: The physical processes in an Ar/O2 magnetron discharge used for the reactive sputter deposition of TiOx thin films were simulated with a 2d3v particle-in-cell/Monte Carlo collisions (PIC/MCC) model. The plasma species taken into account are electrons, Ar+ ions, fast Arf atoms, metastable Arm* atoms, Ti+ ions, Ti atoms, O+ ions, O2+ ions, O− ions and O atoms. This model accounts for plasmatarget interactions, such as secondary electron emission and target sputtering, and the effects of target poisoning. Furthermore, the deposition process is described by an analytical surface model. The influence of the O2/Ar gas ratio on the plasma potential and on the species densities and fluxes is investigated. Among others, it is shown that a higher O2 pressure causes the region of positive plasma potential and the O− density to be more spread, and the latter to decrease. On the other hand, the deposition rates of Ti and O are not much affected by the O2/Ar proportion. Indeed, the predicted stoichiometry of the deposited TiOx film approaches x=2 for nearly all the investigated O2/Ar proportions.