Abstract: In this paper, an overview is given of modeling investigations carried out in our research group for a better understanding of plasmas used for medical, environmental and nano applications. The focus is both on modeling the plasma chemistry and the plasma-surface interactions. The plasma chemistry provides the densities and fluxes of the important plasma species. This information can be used as input when modeling the plasma-surface interactions. The combination of plasma simulations and plasma – surface interaction simulations provides a more comprehensive understanding of the underlying processes for these applications.

Abstract: We present here some of our modeling efforts for PECVD growth of nanostructured carbon materials with focus on amorphous hydrogenated carbon. Experimental data from an expanding thermal plasma setup were used as input for the simulations. Attention was focused both on the film growth mechanism, as well as on the hydrocarbon reaction mechanisms during growth of the films. 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. The film growth results are in correspondence with the experiment. Furthermore, it is found that thin a-C:H films can be densified using an additional H-flux towards the substrate.

Abstract: Plasma-based CO2 conversion is gaining increasing interest worldwide. A large research effort is devoted to improving the energy efficiency. For this purpose, it is very important to understand the underlying mechanisms of the CO2 conversion. The latter can be obtained by computer modeling, describing in detail the behavior of the various plasma species and all relevant chemical processes. However, the accuracy of the modeling results critically depends on the accuracy of the assumed input data, like cross sections. This is especially true for the cross section of electron impact dissociation, as the latter process is believed
to proceed through electron impact excitation, but it is not clear from the literature which excitation channels effectively lead to dissociation. Therefore, the present paper discusses the effect of different electron impact dissociation cross sections reported in the literature on the calculated CO2 conversion, for a dielectric barrier discharge (DBD) and a microwave (MW) plasma. Comparison is made to experimental data for the DBD case, to elucidate which cross section might be the most realistic. This comparison reveals that the cross sections proposed
by Itikawa and by Polak and Slovetsky both seem to underestimate the CO2 conversion. The cross sections recommended by Phelps with thresholds of 7 eV and 10.5 eV yield a CO2 conversion only slightly lower than the experimental data, but the sum of both cross sections overestimates the values, indicating that these cross sections represent dissociation, but most probably also include other (pure excitation) channels. Our calculations indicate that the choice of the electron impact dissociation cross section is crucial for the DBD, where this process is the dominant mechanism for CO2 conversion. In the MW plasma, it is only significant at pressures up to 100 mbar, while it is of minor importance for higher pressures, when dissociation proceeds mainly through collisions of CO2 with heavy particles.

Abstract: Modeling results are presented to gain a better insight in the properties of a SiH4/O2/Ar inductively coupled plasma (ICP) and how it interacts with a silicon substrate (wafer), as applied in the microelectronics industry for the fabrication of electronic devices. The SiH4/O2/Ar ICP is used for the filling of microtrenches with isolating material (SiO2), as applied in shallow trench isolation (STI). In this article, a detailed reaction set that describes the plasma chemistry of SiH4/O2/Ar discharges as well as surface processes, such as sputtering, oxidation, and deposition, is presented. Results are presented on the plasma properties during the plasma enhanced chemical vapor deposition process (PECVD) for different gas ratios, as well as on the shape of the filled trenches and the surface compositions of the deposited layers. For the operating conditions under study it is found that the most important species accounting for deposition are SiH2, SiH3O, SiH3 and SiH2O, while SiH+2, SiH+3, O+2 and Ar+ are the dominant species for sputtering of the surface. By diluting the precursor gas (SiH4) in the mixture, the deposition rate versus sputtering rate can be controlled for a desired trench filling process. From the calculation results it is clear that a high deposition rate will result in undesired void formation during the trench filling, while a small deposition rate will result in undesired trench bottom and mask damage by sputtering. By varying the SiH4/O2 ratio, the chemical composition of the deposited layer will be influenced. However, even at the highest SiH4/O2 ratio investigated (i.e., 3.2:1; low oxygen content), the bulk deposited layer consists mainly of SiO2, suggesting that low-volatile silane species deposit first and subsequently become oxidized instead of being oxidized first in the plasma before deposition. Finally, it was found that the top surface of the deposited layer contained less oxygen due to preferential sputtering of O atoms, making the top layer more Si-rich. However, this effect is negligible at a SiH4/O2 ratio of 2:1 or lower.

Abstract: CO2 dissociation and its subsequent conversion into added-value chemicals is a promising strategy for recycling CO2 gas into reusable products. One of the possible methods is direct plasma-induced dissociation. In this work we study the efficiency of CO2 dissociation in pulsed atmospheric-pressure gas discharge between two conducting electrodes by a 0-D numerical plasma model. The purpose of the study is to provide results on the optimal conditions of CO2 conversion with respect to the energy efficiency and dissociation by varying the maximum power density value and the pulse length. The power density is directly related to the discharge current and the reduced electric field in the discharge. We consider pulse lengths in the range from hundreds of nanosecond up to milliseconds. The results obtained show that the dissociation degree and energy efficiency are sensitive to the pulse length (duration) and the power density, so that a considerable improvement of the discharge performance can be achieved by fine-tuning these parameters. The study is intended to provide guidance in designing an experimental set-up and a power supply with the characteristics necessary to achieve optimal conversion.

Abstract: A theoretical study of ultrashort laser-induced electron emission from a negatively biased metallic cathode has been performed. Classical as well as tunneling electron emission mechanisms are considered. It was found that electron emission is governed by an interplay of processes inside as well as above the cathode. A hybrid model is proposed, where the electron distribution within the target is retrieved from Boltzmann scattering integrals, while the charge distribution above it is studied by a Particle-In-Cell simulation. The results indicate that non-equilibrium effects determine the initial emission process, whereas the space charge above the target suppresses the effectively emitted charge.

Abstract: One of the main issues in plasma chemistry modeling is that the cross sections and rate coefficients are subject to uncertainties, which yields uncertainties in the modeling results and hence hinders the predictive capabilities. In this paper, we reveal the impact of these uncertainties on the model predictions of plasma-based dry reforming in a dielectric barrier discharge. For this purpose, we performed a detailed uncertainty analysis and sensitivity study. 2000 different combinations of rate coefficients, based on the uncertainty from a log-normal distribution, are used to predict the uncertainties in the model output. The uncertainties in the electron density and electron temperature are around 11% and 8% at the maximum of the power deposition for a 70% confidence level. Still, this can have a major effect on the electron impact rates and hence on the calculated conversions of CO2 and CH4, as well as on the selectivities of CO and H2. For the CO2 and CH4 conversion, we obtain uncertainties of 24% and 33%, respectively. For the CO and H2 selectivity, the corresponding uncertainties are 28% and 14%, respectively. We also identify which reactions contribute most to the uncertainty in the model predictions. In order to improve the accuracy and reliability of plasma chemistry models, we recommend using only verified rate coefficients, and we point out the need for dedicated verification experiments.

Abstract: Given the ecological problems associated to the CO2 emissions of fossil fuels, plasma technology has gained
interest for conversion of CO2 into value-added products. Microwave plasmas operating at atmospheric pressure
have proven to be especially interesting, due to the high gas temperatures inside the reactor (i.e. up to 6000 K)
allowing for efficient thermal dissociation of CO2 into CO and O2. However, the performance of these high
temperature plasmas is limited by recombination of CO back into CO2 once the gas cools down in the afterglow.
In this work, we computationally investigated several quenching nozzles, developed and experimentally tested
by Hecimovic et al., [1] for their ability to quickly cool the gas after the plasma, thereby quenching the CO
recombination reactions. Using a 3D computational fluid dynamics model and a quasi-1D chemical kinetics
model, we reveal that a reactor without nozzle lacks gas mixing between hot gas in the center and cold gas near
the reactor walls. Especially at low flow rates, where there is an inherent lack of convective cooling due to the
low gas flow velocity, the temperature in the afterglow remains high (between 2000 and 3000 K) for a relatively
long time (in the 0.1 s range). As shown by our quasi-1D chemical kinetics model, this results in a important loss
of CO due to recombination reactions. Attaching a nozzle in the effluent of the reactor induces fast gas quenching
right after the plasma. Indeed, it introduces (i) more convective cooling by forcing cool gas near the walls to mix
with hot gas in the center of the reactor, as well as (ii) more conductive cooling through the water-cooled walls of
the nozzle. Our model shows that gas quenching and the suppression of recombination reactions have more
impact at low flow rates, where recombination is the most limiting factor in the conversion process.

Abstract: A chemical kinetics model was developed to characterise the gas‐phase dynamics of H₂production in nanosecond‐pulsed CH₄plasmas. Pulsed behaviour was observed in the calculated electric field, electron temperature and species densities at all pressures. The model agrees reasonably with experimental results, showing CH₄conversion at 30% and C₂H₂and H₂as major products. The underlying mechanisms in CH₄dissociation and H₂formation were analysed, highlighting the large contribution of vibrationally excited CH₄and H₂to coupling energy from the plasma into gas‐phase heating, and revealing that H₂synthesis is not affected by applied pressure, with selectivity remaining unchanged at ~42% in the 1–5 bar range.