Abstract: Numerous applications have required the study of CO2 plasmas since the 1960s, from CO2 lasers to spacecraft heat shields. However, in recent years, intense research activities on the subject have restarted because of environmental problems associated with CO2 emissions. The present review provides a synthesis of the current state of knowledge on the physical chemistry of cold CO2 plasmas. In particular, the different modeling approaches implemented to address specific aspects of CO2 plasmas are presented. Throughout the paper, the importance of conducting joint experimental, theoretical and modeling studies to elucidate the complex couplings at play in CO2 plasmas is emphasized. Therefore, the experimental data that are likely to bring relevant constraints to the different modeling approaches are first reviewed. Second, the calculation of some key elementary processes obtained with semi-empirical, classical and quantum methods is presented. In order to describe the electron kinetics, the latest coherent sets of cross section satisfying the constraints of “electron swarm” analyses are introduced, and the need for self-consistent calculations for determining accurate electron energy distribution function (EEDF) is evidenced. The main findings of the latest zero-dimensional (0D) global models about the complex chemistry of CO2 and its dissociation products in different plasma discharges are then given, and full state-to-state (STS) models of only the vibrational-dissociation kinetics developed for studies of spacecraft shields are described. Finally, two important points for all applications using CO2 containing plasma are discussed: the role of surfaces in contact with the plasma, and the need for 2D/3D models to capture the main features of complex reactor geometries including effects induced by fluid dynamics on the plasma properties. In addition to bringing together the latest advances in the description of CO2 non-equilibrium plasmas, the results presented here also highlight the fundamental data that are still missing and the possible routes that still need to be investigated.

Abstract: The reaction behavior of species that may affect the growth of ultrananocrystal line and nanocrystalline diamond ((U)NCD) films is investigated by means of molecular dynamics simulations. Impacts of CHx (x = 0 – 4), C2Hx (x=0-6), C3Hx (x=0-2), C4Hx (x = 0 – 2), H, and H-2 on clean and hydrogenated diamond (100)2 x 1 and (111) 1 x 1 surfaces at two different substrate temperatures are simulated. We find that the different bonding structures of the two surfaces cause different temperature effects on the sticking efficiency. These results predict a temperature-dependent ratio of diamond (100) and (111) growth. Furthermore, predictions of which are the most important hydrocarbon species for (U)NCD growth are made.

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: Plasma is gaining increasing interest for cancer
treatment, but the underlying mechanisms are not yet fully
understood. Using computer simulations at the molecular
level, we try to gain better insight in how plasma-generated
reactive oxygen and nitrogen species (RONS) can
penetrate through the cell membrane. Specifically, we
compare the permeability of various (hydrophilic and
hydrophobic) RONS across both oxidized and nonoxidized cell membranes. We also study pore formation,
and how it is hampered by higher concentrations of
cholesterol in the cell membrane, and we illustrate the
much higher permeability of H2O2 through aquaporin
channels. Both mechanisms may explain the selective
cytotoxic effect of plasma towards cancer cells. Finally, we
also discuss the synergistic effect of plasma-induced
oxidation and electric fields towards pore formation.
Keywords plasma medicine, cancer treatment, computer
modelling, cell membrane, reactive oxygen and nitrogen
species

Abstract: In this work, a detailed structural and spectroscopic study of nanocrystalline diamond (NCD) thin films grown by a continuous bias assisted CVD growth technique is reported. This technique allows the tuning of grain size and phase purity in the deposited material. The crystalline properties of the films are characterized by SEM, TEM, EELS, and Raman spectroscopy. A clear improvement of the crystalline structure of the nanograined diamond film is observed for low negative bias voltages, while high bias voltages lead to thin films consisting of diamond grains of only ∼10 nm nanometer in size, showing remarkable similarities with so-called ultrananocrystalline diamond. These layers arecharacterized by an increasing amount of sp2-bonded carbon content of the matrix in which the diamond grains are embedded. Classical molecular dynamics simulations support the observed experimental data, giving insight in the underlying mechanism for the observed increase in deposition rate with bias voltage. Furthermore, a high atomic concentration of hydrogen has been determined in these films. Finally, Raman scattering analyses confirm that the Raman line observed at ∼1150 cm−1 cannot be attributed to trans-poly-acetylene, which continues to be reported in literature, reassigning it to a deformation mode of CHx bonds in NCD.

Abstract: A molecular dynamics (MD) study of MgxAlyOz thin films grown by magnetron sputtering is presented using an ionic model and comparing two potential sets with formal and partial charges. The applicability of the model and the reliability of the potential sets for the simulation of thin film growth are discussed. The formal charge potential set was found to reproduce the thin film structure in close agreement with the structure of the experimentally grown thin films. Graphical abstract A molecular dynamics study of growth of MgxAlyOz thin films is presented using an ionic model and comparing two potential sets with formal and partial charges. The simulation results with the formal charge potential set showed a transition in the film from a crystalline to an amorphous structure, when the Mg metal content decreases below 50% in very close agreement with the structure of the experimentally deposited films.

Abstract: The compositional influence of Cr and Y on the microstructure of Mg―Cr―O, and Mg―Y―O films synthesized by reactive magnetron sputtering has been investigated by transmission electron microscopy, X-ray diffraction and molecular dynamics simulations. A decrease in crystallinity is observed in these films as the M (Cr or Y) content is increased. It is found that M forms a solid solution with MgO for metal ratios up to ~ 70% and ~ 50% for Cr and Y respectively. Above ~ 70% Cr metal ratio the Mg―Cr―O films are found to be completely amorphous. The Mg―Y―O films are composed of Mg(Y)O and Y2O3 nano crystallites, up to ~ 50% Y metal ratio. Above this ratio, only Y2O3 nano crystallites are found. The preferential < 111> MgO grain alignment is strongly affected by the increase in M content. For M metal ratios up to ~ 50%, there is a selective promotion of the < 100> MgO grain alignments and a decline in the < 111> grain alignments.

Abstract: In this contribution, we present an overview of the various techniques for combining atomistic molecular dynamics with Monte Carlo simulations, mainly in the context of condensed matter systems, as well as a brief summary of the main accelerated dynamics techniques. Special attention is given to the force bias Monte Carlo technique and its combination with molecular dynamics, in view of promising recent developments, including a definable timescale. Various examples of the application of combined molecular dynamics / Monte Carlo simulations are given, in order to demonstrate the enhanced simulation efficiency with respect to either pure molecular dynamics or Monte Carlo.