Abstract: 200 years ago, humanity started the industrial revolution by discovering fossil fuels, which lead to unprecedented technological advancements. However it has become alarmingly clear that the major environmental concerns associated with fossil fuels require a short-term transition from a carbon-based energy economy to a sustainable one based on green electricity. A key step concerning this transition exists in developing electricity-driven alternatives for chemical processes that rely on fossil fuels as a raw material. A technology that is gaining increasing interest to achieve this, is plasma technology. Using plasmas to induce chemical reactions by selectively heating electrons in a gas has already delivered promising results for gas conversion applications like CO2 conversion and N2 fixation, but plasma reactors still require optimization to be considered industrially competitive to existing fossil fuel-based processes and emerging other electricity-based technologies. In this thesis I develop computational models to describe plasma reactors and identify key mechanisms in three different plasma reactors for three different gas conversion applications, i.e. N2 fixation, combined CO2-CH4 conversion and CO2 splitting. I first developed models to describe a new rotating gliding arc (GA) reactor operating in two arc modes, which, as revealed by my model, are characterized by distinct plasma chemistry pathways. Subsequently, my colleague and I study the quenching effect of an effusion nozzle to this rotating GA reactor, reaching the best results to date for N2 fixation into NOx at atmospheric pressure, i.e., NOx concentrations up to 5.9%, at an energy cost down to 2.1 MJ/mol. Afterwards, I investigate the possible improvement of N2 admixtures in plasma-based CO2 and CH4 conversion, as significant amounts of N2 are often found in industrial CO2 waste streams, and gas separations are financially costly. Through combining my models with the experiment from a fellow PhD student, we reveal that moderate amounts of N2 (i.e. around 20%) increase both the electron density and the gas temperature to yield an overall energy cost reduction of 21%. Finally, I model quenching nozzles for plasma-based CO2 conversion in a microwave reactor, to explain the enhancements in CO2 conversion that were demonstrated in experiments. Through computational modelling I reveal that the nozzle introduces fast gas quenching resulting in the suppression of recombination reactions, which have more impact at low flow rates, where recombination is the most limiting factor in the conversion process.

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: Research in plasma reactor designs is developing rapidly as plasma technology is gaining increasing interest for sustainable gas conversion applications, like the conversion of greenhouse gases into value-added chemicals and renewable fuels, and fixation of N₂from air into precursors of mineral fertilizer. As plasma is generated by electric power and can easily be switched on/off, these applications allows for efficient conversion and energy storage of intermittent renewable electricity. In this paper, we present a new comprehensive modelling approach for the design and development of gliding arc plasma reactors, which reveals the fluid dynamics, the arc behaviour and the plasma chemistry by solving a unique combination of five complementary models. This results in a complete description of the plasma process, which allows one to efficiently evaluate the performance of a reactor and indicate possible design improvements before actually building it. We demonstrate the capabilities of this method for an experimentally validated study of plasma-based NO_xformation in a rotating gliding arc reactor, which is gaining increasing interest as a flexible, electricity-driven alternative for the Haber–Bosch process. The model demonstrates the importance of the vortex flow and the presence of a recirculation zone in the reactor, as well as the formation of hot spots in the plasma near the cathode pin and the anode wall that are responsible for most of the NO_xformation. The model also reveals the underlying plasma chemistry and the vibrational non-equilibrium that exists due to the fast cooling during each arc rotation. Good agreement with experimental measurements on the studied reactor design proves the predictive capabilities of our modelling approach.

Abstract: Plasma is gaining increasing interest for N2 fixation, being a flexible, electricity-driven alternative for the current conventional fossil fuel-based N2 fixation processes. As the vibrational-induced dissociation of N2 is found to be an energy-efficient pathway to acquire atomic N for the fixation processes, plasmas that are in vibrational nonequilibrium seem promising for this application. However, an important challenge in using nonequilibrium plasmas lies in preventing vibrational−translational (VT) relaxation processes, in which vibrational energy crucial for N2 dissociation is lost to gas heating. We present here both experimental and modeling results for the vibrational and gas temperature in a microsecond-pulsed microwave (MW) N2 plasma, showing how power pulsing can suppress this unfavorable VT relaxation and achieve a maximal vibrational nonequilibrium. By means of our kinetic model, we demonstrate that pulsed plasmas take advantage of the long time scale on which VT processes occur, yielding a very pronounced nonequilibrium over the whole N2 vibrational ladder. Additionally, the effect of pulse parameters like the pulse frequency and pulse width are investigated, demonstrating that the advantage of pulsing to inhibit VT relaxation diminishes for high pulse frequencies (around 7000 kHz) and long power pulses (above 400 μs). Nevertheless, all regimes studied here demonstrate a clear vibrational nonequilibrium while only requiring a limited power-on time, and thus, we may conclude that a pulsed plasma seems very interesting for energyefficient vibrational excitation.

Abstract: A quantitative HPLC method was applied to determine the amounts of C-60 and C-70 present in extracts of soot produced in the electric arc reactor and in flames. The combustion method was found to yield a higher C-70/C-60 ratio (0.67) compared with the evaporation experiment where the C-70/C-60 ratio amounts to 0.27.

Abstract: To gain full control over the growth of carbon nanotubes (CNTs) using plasma-enhanced chemical vapor deposition (PECVD), a thorough understanding of the underlying plasma-catalyst mechanisms is required. Oxygen-containing species are often used as or added to the growth precursor gas, but these species also yield various radicals and ions, which may simultaneously etch the CNT during the growth. At present, the effect of these reactive species on the growth onset has not yet been thoroughly investigated. We here report on the etching mechanism of incipient CNT structures from OH and O radicals as derived from combined (reactive) molecular dynamics (MD) and force-bias Monte Carlo (tfMC) simulations. Our results indicate that the oxygen-containing radicals initiate a dissociation process. In particular, we show how the oxygen species weaken the interaction between the CNT and the nanocluster. As a result of this weakened interaction, the CNT closes off and dissociates from the cluster in the form of a fullerene. Beyond the specific systems studied in this work, these results are generically important in the context of PECVD-based growth of CNTs using oxygen-containing precursors.

Abstract: The idea of stacking multiple monolayers of different two-dimensional materials has become a global pursuit. In this work, a silicene armchair nanoribbon of width W and van der Waals-bonded to different transition-metal dichalcogenide (TMD) bilayer substrates MoX2 and WX2, where X = S, Se, Te is considered. The orbital resolved electronic structure and ballistic transport properties of these systems are simulated by employing van der Waals-corrected density functional theory and nonequilibrium Green's functions. We find that the lattice mismatch with the underlying substrate determines the electronic structure, correlated with the silicene buckling distortion and ultimately with the contact resistance of the two-terminal system. The smallest lattice mismatch, obtained with the MoTe2 substrate, results in the silicene ribbon properties coming close to those of a freestanding one. With the TMD bilayer acting as a dielectric layer, the electronic structure is tunable from a direct to an indirect semiconducting layer, and subsequently to a metallic electronic dispersion layer, with a moderate applied perpendicular electric field.

Abstract: This study presents first-principles results on the electronic functionalization of silicene and germanene monolayers by means of chemisorption of adatom species H, Li, F, Sc, Ti, V. Three general adatom-monolayer configurations are considered, each having its distinct effect on the electronic structure, yielding metallic or semiconducting dispersions depending on the adatom species and configuration. The induced bandgap is a (in)direct F gap ranging from 0.2 to 2.3 eV for both silicene and germanene. In general the alternating configuration was found to be the most energetically stable. The boatlike and chairlike conformers are degenerate with the former having anisotropic effective carrier masses. The top configuration leads to the planar monolayer and predominately to a gapped dispersion. The hollow configuration with V adatoms retains the Dirac cone, but with strong orbital planar hybridization at the Fermi level. We also observe a planar surface state the Fermi level for the latter systems. (C) 2013 Elsevier B.V. All rights reserved.

Abstract: Plasma medicine has been claimed to provide a novel route to heal wounds and regenerate skin, although very little is currently known about the elementary processes taking place. We carried out a series of ReaxFF-based reactive molecular dynamics simulations to investigate the interaction of O and OH radicals with lipids, more specifically with α-linolenic acid as a model for the free fatty acids present in the upper skin layer. Our calculations predict that the O and OH radicals most typically abstract a H atom from the fatty acids, which can lead to the formation of a conjugated double bond, but also to the incorporation of alcohol or aldehyde groups, thereby increasing the hydrophilic character of the fatty acids and changing the general lipid composition of the skin. Within the limitations of the investigated model, no formation of possibly toxic products was observed.

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: In recent years, the ability of cold atmospheric pressure plasmas (CAPS) to selectively induce cell death in cancer cells has been widely established. This selectivity has been assigned to the reactive oxygen and nitrogen species (RONS) created in CAPs. To provide new insights in the search for an explanation
for the observed selectivity, we calculate the transfer free energy of multiple ROS across membranes containing a varying amount of cholesterol. The cholesterol fraction is investigated as a selectivity parameter because membranes of cancer cells are known to contain lower fractions of cholesterol compared to healthy cells. We find that cholesterol has a significant effect on the permeation of
reactive species across a membrane. Indeed, depending on the specific reactive species, an increasing cholesterol fraction can lead to (i) an increase of the transfer free energy barrier height and width, (ii) the formation of a local free energy minimum in the center of the membrane and (iii) the creation of extra free energy barriers due to the bulky sterol rings. In the context of plasma oncology, these observations suggest that the increased ingress of RONS in cancer cells can be explained by the decreased cholesterol fraction of their cell membrane.