Abstract: In this study we report on a 2D fluid model of a gliding arc discharge in argon. Despite the 3D nature of the discharge, 2D models are found to be capable of providing very useful information about the operation of the discharge. We employ two modelsan axisymmetric and a Cartesian one. We show that for the considered experiment and the conditions of a low current arc (around 30 mA) in argon, there is no significant heating of the cathode surface and the discharge is sustained by field electron emission from the cathode accompanied by the formation of a cathode spot. The obtained discharge power and voltage are relatively sensitive to the surface properties and particularly to the surface roughness, causing effectively an amplification of the normal electric field. The arc body and anode region are not influenced by this and depend mainly on the current value. The gliding of the arc is modelled by means of a 2D Cartesian model. The arcelectrode contact points are analysed and the gliding mechanism along the electrode surface is discussed. Following experimental observations, the cathode spot is simulated as jumping from one point to another. A complete arc cycle is modelled from initial ignition to arc decay. The results show that there is no interaction between the successive gliding arcs.

Abstract: A combined Monte Carlo (MC)/analytical surface model is developed to study the plasma processes occurring during the reactive sputter deposition of TiOx thin films. This model describes the important plasma species with a MC approach (i.e. electrons, Ar+ ions, {\rm O}_2

Abstract: In this paper, an overview is given of different modeling approaches used for describing gas discharge plasmas, as well as plasma-surface interactions. A fluid model is illustrated for describing the detailed plasma chemistry in capacitively coupled rf discharges. The strengths and limitations of Monte Carlo simulations and of a particle-in-cell-Monte Carlo collisions model are explained for a magnetron discharge, whereas the capabilities of a hybrid Monte Carlo-fluid approach are illustrated for a direct current glow discharge used for spectrochemical analysis of materials. Finally, some examples of molecular dynamics simulations, for the purpose of plasma-deposition, are given.

Abstract: In this paper, an O2 inductively coupled plasma used for plasma enhanced atomic layer deposition of Al2O3 thin films is investigated by means of modeling. This work intends to provide more information about basic plasma properties such as species densities and species fluxes to the substrate as a function of power and pressure, which might be hard to measure experimentally. For this purpose, a hybrid model developed by Kushner et al is applied to calculate the plasma characteristics in the reactor volume for different chamber pressures ranging from 1 to 10 mTorr and different coil powers ranging from 50 to 500 W. Density profiles of the various oxygen containing plasma species are reported as well as fluxes to the substrate under various operating conditions. Furthermore, different orientations of the substrate, which can be placed vertically or horizontally in the reactor, are taken into account. In addition, special attention is paid to the recombination process of atomic oxygen on the different reactor walls under the stated operating conditions. From this work it can be concluded that the plasma properties change significantly in different locations of the reactor. The plasma density near the cylindrical coil is high, while it is almost negligible in the neighborhood of the substrate. Ion and excited species fluxes to the substrate are found to be very low and negligible. Finally, the orientation of the substrate has a minor effect on the flux of O2, while it has a significant effect on the flux of O. In the horizontal configuration, the flux of atomic oxygen can be up to one order of magnitude lower than in the vertical configuration.

Abstract: The initial stage of nanoparticle formation and growth in radiofrequency acetylene (C2H2) plasmas is investigated by means of a self-consistent one-dimensional fluid model. A detailed chemical kinetic scheme, containing electron impact, ion-neutral, and neutral-neutral reactions, has been developed in order to predict the underlying dust growth mechanisms and the most important dust precursors. The model considers 41 different species (neutrals, radicals, ions, and electrons) describing hydrocarbons (CnHm) containing up to 12 carbon atoms. Possible routes for particle growth are discussed. Both positive and negative ion reaction pathways are considered, as consecutive anion- and cation-molecule reactions seem to lead to a fast build up of the carbon skeleton.

Abstract: The conversion of methane to value-added chemicals and fuels is considered to be one of the challenges of the 21st century. In this paper we study, by means of fluid modeling, the conversion of methane to higher hydrocarbons or oxygenates by partial oxidation with CO2 or O2 in a dielectric barrier discharge. Sixty-nine different plasma species (electrons, ions, molecules, radicals) are included in the model, as well as a comprehensive set of chemical reactions. The calculation results presented in this paper include the conversion of the reactants and the yields of the reaction products as a function of residence time in the reactor, for different gas mixing ratios. Syngas (i.e. H2 + CO) and higher hydrocarbons (C2Hx) are typically found to be important reaction products.

Abstract: Using Brownian dynamics simulations, we investigate the dynamics of colloids confined in two-dimensional narrow channels driven by a nonuniform force Fdr(y). We considered linear-gradient, parabolic, and deltalike driving-force profiles. This driving force induces melting of the colloidal solid (i.e., shear-induced melting), and the colloidal motion experiences a transition from elastic to plastic regime with increasing Fdr. For intermediate Fdr (i.e., in the transition region) the response of the system, i.e., the distribution of the velocities of the colloidal chains õi(y), in general does not coincide with the profile of the driving force Fdr(y), and depends on the magnitude of Fdr, the width of the channel, and the density of colloids. For example, we show that the onset of plasticity is first observed near the boundaries while the motion in the central region is elastic. This is explained by: (i) (in)commensurability between the chains due to the larger density of colloids near the boundaries, and (ii) the gradient in Fdr. Our study provides a deeper understanding of the dynamics of colloids in channels and could be accessed in experiments on colloids (or in dusty plasma) with, e.g., asymmetric channels or in the presence of a gradient potential field.

Abstract: This paper demonstrates that the CO2 conversion in a dielectric barrier discharge rises drastically upon addition of Ar or He, and the effect is more pronounced for Ar than for He. The effective CO2 conversion, on the other hand, drops upon addition of Ar or He, which is logical due to the lower CO2 content in the gas mixture, and the same is true for the energy efficiency, because a considerable fraction of the energy is then consumed into ionization/excitation of Ar or He atoms. The higher absolute CO2 conversion upon addition of Ar or He can be explained by studying in detail the Lissajous plots and the current profiles. The breakdown voltage is lower in the CO2/Ar and CO2/He mixtures, and the discharge gap is more filled with plasma, which enhances the possibility for CO2 conversion. The rates of electron impact excitationdissociation of CO2, estimated from the electron densities and mean electron energies, are indeed higher in the CO2/Ar and (to a lower extent) in the CO2/He mixtures, compared to the pure CO2 plasma. Moreover, charge transfer between Ar+ or Ar2+ ions and CO2, followed by electron-ion dissociative recombination of the CO2+ ions, might also contribute to, or even be dominant for the CO2 dissociation. All these effects can explain the higher CO2 conversion, especially upon addition of Ar, but also upon addition of He.

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: Electron acceleration from the interaction of an intense short-pulse laser with low density plasma is considered. The relation between direct electron acceleration within the laser pulse and that in the wake is investigated analytically. The magnitude and location of the ponderomotive-force-caused charge separation field with respect to that of the pulse determine the relative effectiveness of the two acceleration mechanisms. It is shown that there is an optimum condition for acceleration in the wake. Electron acceleration within the pulse dominates as the pulse becomes sufficiently short, and the latter directly drives and even traps the electrons. The latter can reach ultrahigh energies and can be extracted by impinging the pulse on a solid target. (C) 2003 American Institute of Physics.

Abstract: In this paper, we investigate the etch process of so-called low-k organic material by means of a N2/H2 capacitively coupled plasma, as applied in the micro-electronics industry for the manufacturing of computer chips. In recent years, such an organic material has emerged as a possible alternative for replacing bulk SiO2 as a dielectric material in the back-end-of-line, because of the smaller parasitic capacity between adjacent conducting lines, and thus a faster propagation of the electrical signals throughout the chip. Numerical simulations with a hybrid plasma model, using an extensive plasma and surface chemistry set, as well as experiments are performed, focusing on the plasma properties as well as the actual etching process, to obtain a better insight into the underlying mechanisms. Furthermore, the effects of gas pressure, applied power and gas composition are investigated to try to optimize the etch process. In general, the plasma density reaches a maximum near the wafer edge due to the so-called 'edge effect'. As a result, the etch rate is not uniform but will also reach its maximum near the wafer edge. The pressure seems not to have a big effect. A higher power increases the etch rate, but the uniformity becomes (slightly) worse. The gas mixing ratio has no significant effect on the etch process, except when a pure H2 or N2 plasma is used, illustrating the synergistic effects of a N2/H2 plasma. In fact, our calculations reveal that the N2/H2 plasma entails an ion-enhanced etch process. The simulation results are in reasonable agreement with the experimental values. The microscopic etch profile shows the desired anisotropic shape under all conditions under study.

Abstract: A fluid model is self-consistently established to investigate the harmonic effects in an inductively coupled plasma, where the electromagnetic field is solved by the finite difference time domain technique. The spatiotemporal distribution of harmonic current density, harmonic potential, and other plasma quantities, such as radio frequency power deposition, plasma density, and electron temperature, have been investigated. Distinct differences in current density have been observed when calculated with and without Lorentz force, which indicates that the nonlinear Lorentz force plays an important role in the harmonic effects, especially at low frequencies. Moreover, the even harmonics are larger than the odd harmonics both in the current density and the potential. Finally, the dependence of various plasma quantities with and without the Lorentz force on various driving frequencies is also examined. It is shown that the deposited power density decreases and the depth of penetration increases slightly because of the Lorentz force. The electron density increases distinctly while the electron temperature remains almost the same when the Lorentz force is taken into account.