Abstract: Ion irradiation of carbon nanotubes very often leads to defect formation. However, we have recently shown that Ar ion irradiation in a limited energy window of 1025 eV may enhance the initial cap nucleation process, when the carbon network is in contact with the metal nanocatalyst. Here, we employ reactive molecular dynamics simulations to demonstrate that ion irradiation in a higher energy window of 1035 eV may also heal network defects after the nucleation stage through a non-metal-mediated mechanism, when the carbon network is no longer in contact with the metal nanocatalyst. The results demonstrate the possibility of beneficially utilizing ions in e.g. plasma-enhanced chemical vapour deposition of carbon nanotubes.

Abstract: Controlling and steering the growth of single walled carbon nanotubes is often believed to require controlling of the nucleation stage. Yet, little is known about the microscopic mechanisms governing the nucleation from hydrocarbon molecules. Specifically, we address here the dehydrogenation of hydrocarbon molecules and the formation of all-carbon graphitic islands on metallic nanoclusters from hydrocarbon molecules under conditions typical for carbon nanotube growth. Employing reactive molecular dynamics simulations, we demonstrate for the first time that the formation of a graphitic network occurs through the intermediate formation of vertically oriented, not fully dehydrogenated graphitic islands. Upon dehydrogenation of these vertical graphenes, the islands curve over the surface, thereby forming a carbon network covering the nanoparticle. The results indicate that controlling the extent of dehydrogenation offers an additional parameter to control the nucleation of carbon nanotubes.

Abstract: In metallized film capacitors, there exists an air gap of about 0.2 μm between the films, with a pressure ranging generally from 130 atm. Because of the created potential difference between the two films, a microdischarge is formed in this gap. In this paper, we use an implicit particle-in-cell Monte Carlo collision simulation method to study the discharge properties in this direct-current microdischarge with 0.2 μm gap in a range of different voltages and pressures. The discharge process is significantly different from a conventional high pressure discharge. Indeed, the high electric field due to the small gap sustains the discharge by field emission. At low applied voltage (~15 V), only the electrons are generated by field emission, while both electrons and ions are generated as a stable glow discharge at medium applied voltage (~50 V). At still higher applied voltage (~100 V), the number of electrons and ions rapidly multiplies, the electric field reverses, and the discharge changes from a glow to an arc regime.

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: In this review paper, we present several examples of reactive molecular dynamics simulations, which contribute to a better understanding of the underlying mechanisms in plasma medicine on the atomic scale. This includes the interaction of important reactive oxygen plasma species with the outer cell wall of both gram-positive and gram-negative bacteria, and with lipids present in human skin. Moreover, as most biomolecules are surrounded by a liquid biofilm, the behavior of these plasma species in a liquid (water) layer is presented as well. Finally, a perspective for future atomic scale modeling studies is given, in the field of plasma medicine in general, and for cancer treatment in particular.

Abstract: We present a zero-dimensional kinetic model of CO2 splitting in non-equilibrium plasmas. The model includes a description of the CO2 vibrational kinetics (25 vibrational levels up to the dissociation limit of the molecule), taking into account state-specific VT and VV relaxation reactions and the effect of vibrational excitation on other chemical reactions. The model is applied to study the reaction kinetics of CO2 splitting in an atmospheric-pressure dielectric barrier discharge (DBD) and in a moderate-pressure microwave discharge. The model results are in qualitative agreement with published experimental works. We show that the CO2 conversion and its energy efficiency are very different in these two types of discharges, which reflects the important dissociation mechanisms involved. In the microwave discharge, excitation of the vibrational levels promotes efficient dissociation when the specific energy input is higher than a critical value (2.0 eV/molecule under the conditions examined). The calculated energy efficiency of the process has a maximum of 23%. In the DBD, vibrationally excited levels do not contribute significantly to the dissociation of CO2 and the calculated energy efficiency of the process is much lower (5%).

Abstract: Using reactive molecular dynamics simulations, the melting behavior of nickelcarbon nanoclusters is examined. The phase diagrams of icosahedral and Wulff polyhedron clusters are determined using both the Lindemann index and the potential energy. Formulae are derived for calculating the equilibrium constants and the solid and liquid fractions during a phase transition, allowing more rational determination of the melting temperature with respect to the arbitrary Lindemann value. These results give more insight into the properties of nickelcarbon nanoclusters in general and can specifically be very useful for a better understanding of the synthesis of carbon nanotubes using the catalytic chemical vapor deposition method.

Abstract: Plasma catalysis holds great promise for environmental applications, provided that the process viability can be maximized in terms of energy efficiency and product selectivity. This requires a fundamental understanding of the various processes taking place and especially the mutual interactions between plasma and catalyst. In this review, we therefore first examine the various effects of the plasma on the catalyst and of the catalyst on the plasma that have been described in the literature. Most of these studies are purely experimental. The urgently needed fundamental understanding of the mechanisms underpinning plasma catalysis, however, may also be obtained through modelling and simulation. Therefore, we also provide here an overview of the modelling efforts that have been developed already, on both the atomistic and the macroscale, and we identify the data that can be obtained with these models to illustrate how modelling and simulation may contribute to this field. Last but not least, we also identify future modelling opportunities to obtain a more complete understanding of the various underlying plasma catalytic effects, which is needed to provide a comprehensive picture of plasma catalysis.

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: Plasma technology is gaining increasing interest for the splitting of CO2 into CO and O2. We have performed experiments to study this process in a dielectric barrier discharge (DBD) plasma with a wide range of parameters. The frequency and dielectric material did not affect the CO2 conversion and energy efficiency, but the discharge gap can have a considerable effect. The specific energy input has the most important effect on the CO2 conversion and energy efficiency. We have also presented a plasma chemistry model for CO2 splitting, which shows reasonable agreement with the experimental conversion and energy efficiency. This model is used to elucidate the critical reactions that are mostly responsible for the CO2 conversion. Finally, we have compared our results with other CO2 splitting techniques and we identified the limitations as well as the benefits and future possibilities in terms of modifications of DBD plasmas for greenhouse gas conversion in general.

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: 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: In traditional capacitively coupled plasmas, the discharge can be described by an electrostatic model, in which the Poisson equation is employed to determine the electrostatic electric field. However, current plasma reactors are much larger and driven at a much higher frequency. If the excitation wavelength k in the plasma becomes comparable to the electrode radius, and the plasma skin depth d becomes comparable to the electrode spacing, the electromagnetic (EM) effects will become significant and compromise the plasma uniformity. In this regime, capacitive discharges have to be described by an EM model, i.e., the full set of Maxwells equations should be solved to address the EM effects. This paper gives an overview of the theory, simulation and experiments that have recently been carried out to understand these effects, which cause major uniformity problems in plasma processing for microelectronics and flat panel display industries. Furthermore, some methods for improving the plasma uniformity are also described and compared.

Abstract: Computer simulations are performed for an argon inductively coupled plasma (ICP) with a capacitive radio-frequency bias power, to investigate the bias effect on the discharge mode transition and on the plasma characteristics at various ICP currents, bias voltages, and bias frequencies. When the bias frequency is fixed at 13.56 MHz and the ICP current is low, e.g., 6A, the spatiotemporal averaged plasma density increases monotonically with bias voltage, and the bias effect is already prominent at a bias voltage of 90 V. The maximum of the ionization rate moves toward the bottom electrode, which indicates clearly the discharge mode transition in inductive/capacitive discharges. At higher ICP currents, i.e., 11 and 13 A, the plasma density decreases first and then increases with bias voltage, due to the competing mechanisms between the ion acceleration power dissipation and the capacitive power deposition. At 11 A, the bias effect is still important, but it is noticeable only at higher bias voltages. At 13 A, the ionization rate is characterized by a maximum at the reactor center near the dielectric window at all selected bias voltages, which indicates that the ICP power, instead of the bias power, plays a dominant role under this condition, and no mode transition is observed. Indeed, the ratio of the bias power to the total power is lower than 0.4 over a wide range of bias voltages, i.e., 0300V. Besides the effect of ICP current, also the effect of various bias frequencies is investigated. It is found that the modulation of the bias power to the spatiotemporal distributions of the ionization rate at 2MHz is strikingly different from the behavior observed at higher bias frequencies. Furthermore, the minimum of the plasma density appears at different bias voltages, i.e., 120V at 2MHz and 90V at 27.12 MHz.