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: Despite the recent promising potential of plasma-based nitrogen fixation, the technology faces significant challenges in efficient upscaling. To tackle this challenge, we investigate two reactors, i.e., a small one, operating in a flow rate range of 5-20 ln min-1 and current range of 200-500 mA, and a larger one, operating at higher flow rate (100-300 ln min-1) and current (400-1000 mA). Both reactors operate in a pin-to-pin configuration and are powered by direct current (DC) from the same power supply unit, to allow easy comparison and evaluate the effect of upscaling. In the small reactor, we achieve the lowest energy cost (EC) of 2.8 MJ mol-1, for a NOx concentration of 1.72%, at a flow rate of 20 ln min-1, yielding a production rate (PR) of 33 g h-1. These values are obtained in air; in oxygen-enriched air, the results are typically better, at the cost of producing oxygen-enriched air. In the large reactor, the higher flow rates reduce the NOx concentration due to lower SEI, while maintaining a similar EC. This stresses the important effect of the geometrical configuration of the arc, which is typically concentrated in the center of the reactor, resulting in limited coverage of the reacting gas flow, and this is identified as the limiting factor for upscaling. However, our experiments reveal that by changing the reactor configuration, and thus the plasma geometry and power deposition mechanisms, the amount of gas treated by the plasma can be enhanced, leading to successful upscaling. To obtain more insights in our experiments, we performed thermodynamic equilibrium calculations. First of all, they show that our measured lowest EC closely aligns with the calculated minimum thermodynamic equilibrium at atmospheric pressure. In addition, they reveal that the limited NOx production in the large reactor results from the contracted nature of the plasma. To solve this limitation, we let the large reactor operate in so-called torch configuration. Indeed, the latter enhances the NOx concentrations compared to the pin-to-pin configuration, yielding a PR of 80 g h-1 at an EC of 2.9 MJ mol-1 and NOx concentration of 0.31%. This illustrates the importance of reactor design in upscaling. With the focus on feasibility evaluation of scaling-up plasma-based nitrogen fixation by combined experiments and thermodynamic modelling, we aim to tackle the challenge of design and development of an energy-efficient and scaled-up plasma reactor.

Abstract: Plasma medicine has attracted tremendous interest in a variety of medical conditions, ranging from wound healing to antimicrobial applications, even in cancer treatment, through the interactions of cold atmospheric plasma (CAP) and various biological tissues directly or indirectly. The underlying mechanisms of CAP treatment are still poorly understood although the oxidative effects of CAP with amino acids, peptides, and proteins have been explored experimentally. In this study, machine learning (ML) technology is introduced to efficiently unveil the interaction mechanisms of amino acids and reactive oxygen species (ROS) in seconds based on the data obtained from the reactive molecular dynamics (MD) simulations, which are performed to probe the interaction of five types of amino acids with various ROS on the timescale of hundreds of picoseconds but with the huge computational load of several days. The oxidative reactions typically start with H-abstraction, and the details of the breaking and formation of chemical bonds are revealed; the modification types, such as nitrosylation, hydroxylation, and carbonylation, can be observed. The dose effects of ROS are also investigated by varying the number of ROS in the simulation box, indicating agreement with the experimental observation. To overcome the limits of timescales and the size of molecular systems in reactive MD simulations, a deep neural network (DNN) with five hidden layers is constructed according to the reaction data and employed to predict the type of oxidative modification and the probability of occurrence only in seconds as the dose of ROS varies. The well-trained DNN can effectively and accurately predict the oxidative processes and productions, which greatly improves the computational efficiency by almost ten orders of magnitude compared with the reactive MD simulation. This study shows the great potential of ML technology to efficiently unveil the underpinning mechanisms in plasma medicine based on the data from reactive MD simulations or experimental measurements. In this study, since reactive molecular dynamics simulation can currently only describe interactions between a few hundred atoms in a few hundred picoseconds, deep neural networks (DNN) are introduced to enhance the simulation results by predicting more data efficiently. image

Abstract: We propose a plasma chemical looping CO2 splitting (PCLCS) approach that enables highly efficient CO2 conversion into O-2-free CO at mild temperatures. PCLCS achieves an impressive 84% CO2 conversion and a 1.3 mmol g(-1) CO yield, with no O-2 detected. Crucially, this strategy significantly lowers the temperature required for conventional chemical looping processes from 650 to 1000 degrees C to only 320 degrees C, demonstrating a robust synergy between plasma and the Ce0.7Zr0.3O2 oxygen carrier (OC). Systematic experiments and density functional theory (DFT) calculations unveil the pivotal role of plasma in activating and partially decomposing CO2, yielding a mixture of CO, O-2/O, and electronically/vibrationally excited CO2*. Notably, these excited CO2* species then efficiently decompose over the oxygen vacancies of the OCs, with a substantially reduced activation barrier (0.86 eV) compared to ground-state CO2 (1.63 eV), contributing to the synergy. This work offers a promising and energy-efficient pathway for producing O-2-free CO from inert CO2 through the tailored interplay of plasma and OCs.

Abstract: A multi-species Monte Carlo (MC) model, combined with an analytical surface model, has been developed in order to investigate the general plasma processes occurring during the sputter deposition of complex oxide films in a dual-magnetron sputter deposition system. The important plasma species, such as electrons, Ar+ ions, fast Ar atoms and sputtered metal atoms (i.e. Mg and Al atoms) are described with the so-called multi-species MC model, whereas the deposition of MgxAlyOz films is treated by an analytical surface model. Targetsubstrate distances for both magnetrons in the dual-magnetron setup are varied for the purpose of growing stoichiometric complex oxide thin films. The metal atoms are sputtered from pure metallic targets, whereas the oxygen flux is only directed toward the substrate and is high enough to obtain fully oxidized thin films but low enough to avoid target poisoning. The calculations correspond to typical experimental conditions applied to grow these complex oxide films. In this paper, some calculation results are shown, such as the densities of various plasma species, their fluxes toward the targets and substrate, the deposition rates, as well as the film stoichiometry. Moreover, some results of the combined model are compared with experimental observations. Note that this is the first complete model, which can be applied for large and complicated magnetron reactor geometries, such as dual-magnetron configurations. With this model, we are able to describe all important plasma species as well as the deposition process. It can also be used to predict film stoichiometries of complex oxide films on the substrate.

Abstract: The electron bounce resonance heating (BRH) in dual-frequency capacitively coupled plasmas operated in oxygen and argon has been studied by different experimental methods. In comparison with the electropositive argon discharge, the BRH in an electronegative discharge occurs at larger electrode gaps. Kinetic particle simulations reveal that in the oxygen discharge, the bulk electric field becomes quite strong and is out of phase with the sheath field. Therefore, it retards the resonant electrons when traversing the bulk, resulting in a suppressed BRH. This effect becomes more pronounced at lower high-frequency power, when the discharge mode changes from electropositive to electronegative.

Abstract: A Langmuir probe and an ICCD are employed to study the discharge mode transition in Ar inductively coupled plasma. Electron density and plasma emission intensity are measured during the E (capacitive discharge) to H (inductive discharge) mode transitions at different pressures. It is found that plasma exists with a low electron density and a weak emission intensity in the E mode, while it has a high electron density and a strong emission intensity in the H mode. Meanwhile, the plasma emission intensity spatial (2D image) profile is symmetrical in the H mode, but the 2D image is an asymmetric profile in the E mode. Moreover, the electron density and emission intensity jump up discontinuously at high pressure, but increase almost continuously at the E to H mode transition under low pressure.

Abstract: In recent years there has been growing interest in the use of low-temperature atmospheric pressure plasmas for biomedical applications. Currently, however, there is very little fundamental knowledge regarding the relevant interaction mechanisms of plasma species with living cells. In this paper, we investigate the interaction of important plasma species, such as O3, O2 and O atoms, with bacterial peptidoglycan (or murein) by means of reactive molecular dynamics simulations. Specifically, we use the peptidoglycan structure to model the gram-positive bacterium Staphylococcus aureus murein. Peptidoglycan is the outer protective barrier in bacteria and can therefore interact directly with plasma species. Our results demonstrate that among the species mentioned above, O3 molecules and especially O atoms can break important bonds of the peptidoglycan structure (i.e. CO, CN and CC bonds), which subsequently leads to the destruction of the bacterial cell wall. This study is important for gaining a fundamental insight into the chemical damaging mechanisms of the bacterial peptidoglycan structure on the atomic scale.

Abstract: An inductively coupled plasma, connected to a mass spectrometer interface, is computationally investigated. The effect of pressure behind the sampler, injector gas flow rate, auxiliary gas flow rate, and applied power is studied. There seems to be an optimum range of injector gas flow rate for each setup which guaranties the presence and also a proper length of the central channel in the torch. Moreover, our modeling results show that for any specific purpose, it is possible to control that either only the central gas flow passes through the sampler orifice or that it is accompanied by the auxiliary gas flow. It was also found that depending on geometry, the variation of outgoing gas flow rate is much less than the variation of the injector gas flow rate and this causes a slightly higher pressure inside the torch. The general effect of increasing the applied power is a rise in the plasma temperature, which results in a higher ionization in the coil region. However, the negative effect is reducing the length of the cool central channel which is important to transfer the sample substances to the sampler. Using a proper applied power can enhance the efficiency of the system. Indeed, by changing the gas path lines, the power can control which flow (i.e., only from injector gas or also from the auxiliary gas) goes to the sampler orifice. Finally, as also reported from experiments in literature, the pressure behind the sampler has no dramatic effect on the plasma characteristics.

Abstract: The adsorption probability and reaction behavior of CHx plasma species on various nickel catalyst surfaces is investigated by means of reactive molecular dynamics (MD) simulations using the ReaxFF potential. Such catalysts are used in the reforming of hydrocarbons and in the growth of carbon nanotubes, and further insight in the underlying mechanisms of these processes is needed to increase their applicability. Single and consecutive impacts of CHx radicals (x={1,2,3}) were performed on four different Ni surfaces, at a temperature of 400 K. The adsorption probability is shown to be related to the number of free electrons, i.e. a higher number leads to more adsorptions, and the steric hindrance caused by the hydrogen atoms bonded to the impacting CHx species. Furthermore, some of the CH bonds break after adsorption, which generally leads to diffusion of the hydrogen atom over the surface. Additionally, these adsorbed H-atoms can be used in reactions to form new molecules, such as CH4 and C2Hx, although this is dependent on the precise morphology of the surface. New molecules are also formed by subtraction of H-atoms from adsorbed radicals, leading to occasional formation of H2 and C2Hx molecules.

Abstract: In this paper, the splitting of CO2 in a pulsed plasma system, such as a dielectric barrier discharge (DBD), is evaluated from a chemical point of view by means of numerical modeling. For this purpose, a chemical reaction set of CO2 in an atmospheric pressure plasma is developed, including the vibrational states of CO2, O2, and CO. The simulated pulses are matched to the conditions of a filament (or microdischarge) and repeated with intervals of 1 μs. The influence of vibrationally excited CO2 as well as other neutral species, ions, and electrons on the CO2 splitting is discussed. Our calculations predict that the electrons have the largest contribution to the CO2 splitting at the conditions under study, by electron impact dissociation. The contribution of vibrationally excited CO2 levels in the splitting of CO2 is found be 6.4%, when only considering one microdischarge pulse and its afterglow, but it can be much higher for consecutive discharge pulses, as is typical for a filamentary DBD, when the interpulse time is short enough and accumulation effects in the vibrationally excited CO2 densities can occur.

Abstract: The exact structure and properties of the Si vertical bar SiO2 interface are very important in microelectronics and photovoltaic devices such as metal-oxide-semiconductor field-effect transistors (MOSFETs) and solar cells. Whereas Si vertical bar SiO2 structures are traditionally produced by thermal oxidation, hyperthermal oxidation shows a number of promising advantages. However, the Si vertical bar SiO2 interface induced in hyperthermal Si oxidation has not been properly investigated yet. Therefore, in this work, the interface morphology and interfacial stresses during hyperthermal oxidation at room temperature are studied using reactive molecular dynamics simulations based on the ReaxFF potential. Interface thickness and roughness, as well as the bond length and bond angle distributions in the interface are discussed and compared with other models developed for the interfaces induced by traditional thermal oxidation. The formation of a compressive stress is observed. This compressive stress, which at the interface amounts about 2 GPa, significantly slows down the inward silica growth. This value is close to the experimental value in the Si vertical bar SiO2 interface obtained in traditional thermal oxidation.