Abstract: This paper focuses on the investigation of CF4 decomposition in a low-pressure inductively coupled plasma by means of a global model. The influence of O2 on the CF4 decomposition process is studied for conditions used in semiconductor manufacturing processes. The model is applied for different powers and O2 contents ranging between 2% and 98% in the CF4/O2 gas mixture. The model includes the reaction mechanisms in the gas phase coupled with the surface reactions and sticking probabilities of the species at the walls. The calculation results are first compared with experimental results from the literature (for the electron density, temperature and F atom density) at a specific power, in the entire range of CF4/O2 gas mixture ratios, and the obtained agreements indicate the validity of the model. The main products of the gas mixture, obtained from this model, include CO, CO2 and COF2 together with a low fraction of F2. The most effective reactions for the formation and loss of the various species in this process are also determined in detail. Decomposition of CF4 produces mostly CF3 and F radicals. These radicals also contribute to the backward reactions, forming again CF4. This study reveals that the maximum decomposition efficiency of CF4 is achieved at a CF4/O2 ratio equal to 1, at the applied power of 300 W.

Abstract: A chemical kinetics model is developed for a CO2/N2 microwave plasma, focusing especially on the vibrational levels of both CO2 and N2. The model is used to calculate the CO2 and N2 conversion as well as the energy efficiency of CO2 conversion for different power densities and for N2 fractions in the CO2/N2 gas mixture ranging from 0 to 90%. The calculation results are compared with measurements, and agreements within 23% and 33% are generally found for the CO2 conversion and N2 conversion, respectively. To explain the observed trends, the destruction and formation processes of both CO2 and N2 are analyzed, as well as the vibrational distribution functions of both CO2 and N2. The results indicate that N2 contributes in populating the lower asymmetric levels of CO2, leading to a higher absolute CO2 conversion upon increasing N2 fraction. However, the effective CO2 conversion drops because there is less CO2 initially present in the gas mixture; thus, the energy efficiency also drops with rising N2 fraction.

Abstract: In this paper, the effects of prolonged bias-enhanced nucleation (prolonged BEN) on the growth mechanisms of diamond are investigated by molecular dynamics (MD) and combined MD-Metropolis Monte Carlo (MD-MMC) simulations. First, cumulative impacts of CxHy+ and Hx+ on an a-C:H/nanodiamond composite were simulated; second, nonconsecutive impacts of the dominant ions were simulated in order to understand the observed phenomena in more detail. As stated in the existing literature, the growth of diamond structures during prolonged BEN is a process that takes place below the surface of the growing film. The investigation of the penetration behavior of CxHy+ and Hx+ species shows that the carbon-containing ions remain trapped within this amorphous phase where they dominate mechanisms like precipitation of sp3 carbon clusters. The H+ ions, however, penetrate into the crystalline phase at high bias voltages (>100 V), destroying the perfect diamond structure. The experimentally measured reduction of grain sizes at high bias voltage, reported in the literature, might thus be related to penetrating H+ ions. Furthermore, the CxHy+ ions are found to be the most efficient sputtering agents, preventing the build up of defective material.

Abstract: Using reactive molecular dynamics simulations, the melting behavior of nickel-carbon 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 nickel-carbon 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: It is known that film structure may change dramatically with the extent of surface diffusion during the film growth process. In the present work, surface diffusion, induced thermally or activated by energetic impacts, is investigated theoretically under conditions appropriate for magnetron sputter-deposition of MgAlO thin films with varying stoichiometry. The distribution of surface diffusion energy barriers available to the system was determined for each stoichiometry, which allowed assessing in a qualitative way how much surface diffusion will take place on the time scale available between deposition events. The activation energy barriers increase with the Al concentration in the film, and therefore, the surface diffusion rates in the time frame of typical deposition rates drop, which can explain the decrease in crystallinity in the film structure and the transition to amorphous structure. The deposition process and the immediate surface diffusion enhanced by the energetic adatoms are simulated by means of a molecular dynamics model. The longer-time thermal surface diffusion and the energy landscape are studied by the temperature accelerated dynamics method, applied in an approximate way. The surface diffusion enhanced by the energetic impacts appears to be very important for the film structure in the low-temperature deposition regime.

Abstract: The vibrational level populations of the electronic ground state of the nitrogen molecule have been calculated for typical glow discharge conditions in argonnitrogen mixtures with nitrogen concentrations between 0.1 and 1%. Stationary solutions of the master equations of the vibrational levels have been obtained using numerical methods. The main mechanisms responsible for the population and depopulation of the vibrational levels, and for the overall shape of the vibrational distribution function are pointed out. It has been found that vibrationvibration collisions play only a minor role and therefore the population of the vibrational levels is basically determined by the electron temperature.

Abstract: The study presents a 2D fluid plasma model of a gliding arc discharge for dissociation of CO2 which allows its subsequent conversion into value-added chemicals. The model is based on the balance equations of charged and neutral particles, the electron energy balance equation, the gas thermal balance equation and the current continuity equation. By choosing the modeling domain to be the plane perpendicular to the arc current, the numerical calculations are significantly simplified. Thus, the model allows us to explore the influence of the gas instabilities (turbulences) on the energy efficiency of CO2 conversion. This paper presents results for plasma parameters at different values of the effective turbulent thermal conductivity leading to enhanced energy transport.

Abstract: The mechanisms of plasma in medicine are broadly attributed to plasma-derived reactive oxygen and nitrogen species (RONS). In order to exert any intracellular effects, these plasma-derived RONS must first traverse a major barrier in the cell membrane. The cell membrane lipid composition, and thereby the magnitude of this barrier, is highly variable between cells depending on type and state (e.g. it is widely accepted that healthy and cancerous cells have different membrane lipid compositions). In this study, we investigate how plasma-derived RONS interactions with lipid membrane components can potentially be exploited in the future for treatment of diseases. We couple phospholipid vesicle experiments, used as simple cell models, with molecular dynamics (MD) simulations of the lipid membrane to provide new insights into how the interplay between phospholipids and cholesterol may influence the response of healthy and diseased cell membranes to plasma-derived RONS. We focus on the (i) lipid tail saturation degree, (ii) lipid head group type, and (iii) membrane cholesterol fraction. Using encapsulated molecular probes, we study the influence of the above membrane components on the ingress of RONS into the vesicles, and subsequent DNA damage. Our results indicate that all of the above membrane components can enhance or suppress RONS uptake, depending on their relative concentration within the membrane. Further, we show that higher RONS uptake into the vesicles does not always correlate with increased DNA damage, which is attributed to ROS reactivity and lifetime. The MD simulations indicate the multifactorial chemical and physical processes at play, including (i) lipid oxidation, (ii) lipid packing, and (iii) lipid rafts formation. The methods and findings presented here provide a platform of knowledge that could be leveraged in the development of therapies relying on the action of plasma, in which the cell membrane and oxidative stress response in cells is targeted.

Abstract: Plasma-treated liquids have great potential for biomedical applications. However, insight into the underlying mechanisms and the exact chemistry is still scarce. In this study, we present the combination of a 0D chemical kinetics and a 2D fluid dynamics model to investigate the plasma treatment of a buffered water solution with the kINPen (R) plasma jet. Using this model, we calculated the gas and liquid flow profiles and the transport and chemistry of all species in the gas and the liquid phase. Moreover, we evaluated the stability of the reactive oxygen and nitrogen species after plasma treatment. We found that of all species, only H2O2, HNO2/NO2-, and HNO3/NO3- are stable in the buffered solution after plasma treatment. This is because both their production and loss processes in the liquid phase are dependent on short-lived radicals (e.g. OH, NO, and NO2). Apart from some discrepancy in the absolute values of the concentrations, which can be explained by the model, all general trends and observations in our model are in qualitative agreement with experimental data and literature.