Abstract: Oligomerized [Cu-O-Cu] species are reported to be efficient in promoting plasma catalytic one-step steam reforming of methane to methanol and hydrogen, achieving 6.8% CH4 conversion and 73.1% CH3OH selectivity without CO2.

Abstract: With the growth of the world population, the agricultural sector is required to meet an increasing demand for nutrients and currently relies on industrially produced fertilizers. Among them, nitrogen-based fertilizers are the most common choice and require N2 to be converted into more reactive molecules in a process called “nitrogen fixation”. This is mainly performed through the Haber-Bosch process, which, is not ideal since it requires large-scale facilities to be economical and is associated with a high energy cost and high CO2 emissions, resulting in an environmental impact that is pushing for the study of greener alternatives. Among these, plasma-based nitrogen fixation into NOx is promising, and gliding arc plasma, specifically, proved to be suitable for nitrogen fixation. This thesis aims to study plasma-based nitrogen fixation focusing on an atmospheric pressure gliding arc plasma on three different levels. On a fundamental level, an approach dealing with laser-based excitation of separate rotational lines was successfully developed. This method can be implemented on atmospheric discharges that produce rather high NOx densities and, thus, can impose essential restrictions for the use of “classical” laser-induced fluorescence methods. The approach is then implemented, providing a discussion on the two-dimensional distributions of both the gas temperature and the NO ground state density. A clear correlation between these quantities is found and the effects of both the gas temperature and the plasma power on NO and NO2 concentrations are discussed, revealing how the NO oxidation is already significant in the plasma afterglow region and how the gas flow rate is a crucial parameter affecting the temperature gradients. >From a technological level, the conventional approach of introducing external resistors to stabilize the arc is challenged by studying both its performance and its stability replacing the external resistor with an inductor. We conclude that similar stabilization results can be obtained while significantly lowering the overall energy cost, which decreased from up to a maximum of 7.9 MJ/mol N to 3 MJ/mol N. Finally, we study whether a small-scale fertilizer production facility based on a gliding arc plasma can be a local competitive alternative. This is done by proposing a comparative model to understand how capital, operative expenditures and transport costs affect the production costs. The model highlights how, with the current best available technology, plasma-based nitrogen fixation, while being an interesting alternative for NOx synthesis, still requires a more efficient use of H2 for direct NH3 production.

Abstract: Atmospheric pressure plasmas have shifted in recent years from being a burgeoning research field in the academic setting to an actively investigated technology in the chemical, oil, and environmental industries. This is largely driven by the climate change mitigation efforts, as well as the evident pathways of value creation by converting greenhouse gases (such as CO2) into useful chemical feedstock. Currently, most high technology readiness level (TRL) plasma-based technologies are based on volumetric and power-based scaling of thermal plasma systems, which results in large capital investment and regular maintenance costs. This work investigates bringing a quasi-thermal (so-called “warm”) plasma setup, namely, a gliding arc plasmatron, from a lab-scale to a pilot-scale capacity with an increase in throughput capacity by a factor of 10. The method of scaling is the parallelization of plasmatron reactors within a single housing, with the aim of maintaining a warm plasma regime while simultaneously improving build cost and efficiency (compared to separate reactors operating in parallel). Special attention is also given to the safety and control features implemented in the setup, a key component required for integration into industrial systems. The performance of the multi-reactor gliding arc plasmatron (MRGAP) reactor is investigated, focusing on the influence of flow rate and the number of active reactors. The location of active reactors was deemed to have a negligible effect on the monitored metrics of conversion, energy efficiency, and energy cost. The optimum operating conditions were found to be with the most active reactors (five) at the highest investigated flow rate (80 L/min). Analysis of results suggests that an optimum conversion (9%) and plug power-based energy efficiency (19%) can be maintained at a specific energy input (SEI) around 5.3 kJ/L (or 1 eV/molecule). The concept of parallelization of plasmatron reactors within a singular housing was demonstrated to be a viable method for scaling up from a lab-scale to a prototype-scale device, with performance analysis suggesting that increasing the power (through adding more reactor channels) and total flow rate, while maintaining an SEI around 5.3 or 4.2 kJ/L, i.e., 1.3 or 1 eV/molecule (based on plug power and plasma-deposited power, respectively), can result in increased conversion rate without sacrificing absolute conversion or energy efficiency.

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: We study the properties of a capacitive 13.56 MHz discharge with a mixture of Ar/C₂H₂ taking into account the plasmochemistry and growth of heavy hydrocarbons. A hybrid model was developed to combine the kinetic description for electron motion and the fluid approach for negative and positive ion transports and plasmochemical processes. A significant change in plasma parameters related to injection of 5.8% portion of acetylene in argon was observed and analyzed. We found that the electronegativity of the mixture is about 30%. The densities of negatively and positively charged heavy hydrocarbons are sufficiently large to be precursors for the formation of nanoparticles in the discharge volume.

Abstract: Cold atmospheric pressure plasmas have proven to provide an alternative treatment of cancer by targeting tumorous cells while leaving their healthy counterparts unharmed. However, the underlying mechanisms of the plasma–cell interactions are not yet fully understood. Reactive oxygen species, and in particular hydroxyl radicals (OH), are known to play a crucial role in plasma driven apoptosis of
malignant cells. In this paper we investigate the interaction of OH radicals, as well as H2O2 molecules and HO2 radicals, with DNA by means of reactive molecular dynamics simulations using the ReaxFF force field. Our results provide atomic-scale insight into the dynamics of oxidative stress on DNA caused by the OH radicals, while H2O2 molecules appear not reactive within the considered timescale. Among the observed processes are the formation of 8-OH-adduct radicals, forming the first stages towards the formation of 8-oxoGua and 8-oxoAde, H-abstraction reactions of the amines, and the partial opening of loose DNA ends in aqueous solution.

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: A packed bed dielectric barrier discharge plasma reactor is computationally studied with a fluid model. Two different complementary axisymmetric 2D geometries are used to mimic the intrinsic 3D problem. It is found that a packing enhances the electric field strength and electron temperature at the contact points of the dielectric material due to polarization of the beads by the applied potential. As a result, these contact points prove to be of direct importance to initiate the plasma. At low applied potential, the discharge stays at the contact points, and shows the properties of a Townsend discharge. When a high enough potential is applied, the plasma will be able to travel through the gaps in between the beads from wall to wall, forming a kind of glow discharge. Therefore, the inclusion of a so-called ‘channel of voids’ is indispensable in any type of packed bed modelling.

Abstract: Modelling at the molecular or atomic scale can be very useful for obtaining a better insight in plasma medicine. This paper gives an overview of different atomic/molecular scale modelling approaches that can be used to study the direct interaction of plasma species with biomolecules or the consequences of these interactions for the biomolecules on a somewhat longer time-scale. These approaches include density functional theory (DFT), density functional based tight binding (DFTB), classical reactive and non-reactive molecular dynamics (MD) and united-atom or coarse-grained MD, as well as hybrid quantum mechanics/molecular mechanics (QM/MM) methods. Specific examples will be given for three important types of biomolecules, present in human cells, i.e. proteins, DNA and phospholipids found in the cell membrane. The results show that each of these modelling approaches has its specific strengths and limitations, and is particularly useful for certain applications. A multi-level approach is therefore most suitable for obtaining a global picture of the plasma–biomolecule interactions.

Abstract: The hyperdynamics method is a powerful tool to simulate slow processes at the atomic level. However, the construction of an optimal hyperdynamics potential is a task that is far from trivial. Here, we propose a generally applicable implementation of the hyperdynamics algorithm, borrowing two concepts from metadynamics. First, the use of a collective variable (CV) to represent the accelerated dynamics gives the method a very large flexibility and simplicity. Second, a metadynamics procedure can be used to construct a suitable history-dependent bias potential on-the-fly, effectively turning the algorithm into a self-learning accelerated molecular dynamics method. This collective variable-driven hyperdynamics (CVHD) method has a modular design: both the local system properties on which the bias is based, as well as the characteristics of the biasing method itself, can be chosen to match the needs of the considered system. As a result, system-specific details are abstracted from the biasing algorithm itself, making it extremely versatile and transparent. The method is tested on three model systems: diffusion on the Cu(001) surface and nickel-catalyzed methane decomposition, as examples of reactive processes with a bond-length-based CV, and the folding of a long polymer-like chain, using a set of dihedral angles as a CV. Boost factors up to 109, corresponding to a time scale of seconds, could be obtained while still accurately reproducing correct dynamics.

Abstract: Device performance of Si|SiOx core-shell based nanowires critically depends on the exact control over the oxide thickness. Low-temperature plasma oxidation is a highly promising alternative to thermal oxidation allowing for improved control over the oxidation process, in particular for ultrasmall Si nanowires. We here elucidate the room temperature plasma oxidation mechanisms of ultrasmall Si nanowires using hybrid molecular dynamics / force-bias Monte Carlo simulations. We demonstrate how the oxidation and concurrent water formation mechanisms are a function of the oxidizing plasma species and we demonstrate how the resulting core-shell oxide thickness can be controlled through these species. A new mechanism of water formation is discussed in detail. The results provide a detailed atomic level explanation of the oxidation process of highly curved Si surfaces. These results point out a route toward plasma-based formation of ultrathin core-shell Si|SiOx nanowires at room temperature.

Abstract: In this paper the various elementary plasma—surface interaction processes occurring in plasma catalysis are critically evaluated. Specifically, plasma catalysis at atmospheric pressure is considered. The importance of the various processes is analyzed for the most common plasma catalysis sources, viz. the dielectric barrier discharge and the gliding arc. The role and importance of surface chemical reactions (including adsorption, surface-mediated association and dissociation reactions, and desorption), plasma-induced surface modification, photocatalyst activation, heating, charging, surface discharge formation and electric field enhancement are discussed in the context of plasma catalysis. Numerous examples are provided to demonstrate the importance of the various processes.

Abstract: Carbon dioxide conversion and utilization has gained significant interest over the years. A novel gas conversion technique with great potential in this area is plasma technology. A lot of research has already been performed, but mostly on pure gases. In reality, N2 will always be an important impurity in effluent
gases. Therefore, we performed an extensive combined experimental and computational study on the effect of N2 in the range of 1–98% on CO2 splitting in dielectric barrier discharge (DBD) plasma. The presence of up to 50% N2 in the mixture barely influences the effective (or overall) CO2 conversion and energy efficiency, because the N2 metastable molecules enhance the absolute CO2 conversion, and this compensates for the lower CO2 fraction in the mixture. Higher N2 fractions, however, cause a drop in the CO2 conversion and energy efficiency. Moreover, in the entire CO2/N2 mixing ratio, several harmful compounds, i.e., N2O and NOx compounds, are produced in the range of several 100 ppm. The reaction pathways for the formation of these compounds are explained based on a kinetic analysis, which allows proposing solutions on how to prevent the formation of these harmful compounds.

Abstract: We studied the transport of copper droplets through an inductively coupled plasma, connected to the sampling cone of a mass spectrometer, by means of a computational model. The sample droplets are followed until they become evaporated. They are inserted as liquid particles from the central inlet and the effects of injection position (i.e. “on” and “off” axis), droplet diameter, as well as mass loading flow rate are investigated. It is shown that more “on-axis” injection of the droplets leads to a more straight path line, so that the droplets move less in the radial direction and are evaporated more on the central axis, enabling a better sample transfer efficiency to the sampler cone. Furthermore, there are optimum ranges of diameters and flow rates, which guarantee the proper position of evaporation along the torch, i.e. not too early, so that the sample can get lost in the torch, and not too late, which reduces the chance of becoming ionized before reaching the sampler.

Abstract: A modelling study is performed of high-density low-pressure inductively coupled CF4/CHF3/H2/Cl2/O2/HBr plasmas under different gas mixing ratios. A reaction set describing the complete plasma chemistry is presented and discussed. The gas fraction of each component in this mixture is varied to investigate the sensitivity of the plasma properties, like electron density, plasma potential and species densities, towards the gas mixing ratios. This research is of great interest for microelectronics applications because these gases are often combined in two (or more)-component mixtures, and mixing gases or changing the fraction of a gas can sometimes yield unwanted reaction products or unexpected changes in the overall plasma properties due to the increased chemical complexity of the system. Increasing the CF4 fraction produces more F atoms for chemical etching as expected, but also more prominently lowers the density of Cl atoms, resulting in an actual drop in the etch rate under certain conditions. Furthermore, CF4 decreases the free electron density when mixed with Cl2. However, depending on the other gas components, CF4 gas can also sometimes enhance free electron density. This is the case when HBr is added to the mixture. The addition of H2 to the gas mixture will lower the sputtering process, not only due to the lower overall positive ion density at higher H2 fractions, but also because more H+, H2 + and H3 + are present and they have very low sputter yields. In contrast, a larger Cl2 fraction results in more chemical etching but also in less physical sputtering due to a smaller abundance of positive ions. Increasing the O2 fraction in the plasma will always lower the etch rate due to more oxidation of the wafer surface and due to a lower plasma density. However, it is also observed that the density of F atoms can actually increase with rising O2 gas fraction. This is relevant to note because the exact balance between fluorination and oxidation is important for fine-tuning the overall etch rate and for control of the sidewall profile. Finally, HBr is often used as a chemical etcher, but when mixed with F- or Cl-containing gases, HBr creates the same diluting effects as Ar or He, because a
higher fraction results in less chemical etching but more (physical) sputtering.