Abstract: Molecular dynamics (MD) simulations were performed to provide atomic scale insight in the initial interaction between hydroxyl radicals (OH) and peptide systems in solution. These OH radicals are representative reactive oxygen species produced by cold atmospheric plasmas. The use of plasma for biomedical applications is gaining increasing interest, but the fundamental mechanisms behind the plasma modifications still remain largely elusive. This study helps to gain more insight in the underlying mechanisms of plasma medicine but is also more generally applicable to peptide oxidation, of interest for other applications. Combining both reactive and nonreactive MD simulations, we are able to elucidate the reactivity of the amino acids inside the peptide systems and their effect on their structure up to 1 μs. Additionally, experiments were performed, treating the simulated peptides with a plasma jet. The computational results presented here correlate well with the obtained experimental data and highlight the importance of the chemical environment for the reactivity of the individual amino acids, so that specific amino acids are attacked in higher numbers than expected. Furthermore, the long time scale simulations suggest that a single oxidation has an effect on the 3D conformation due to an increase in hydrophilicity and intra- and intermolecular interactions.

Abstract: We study the mechanisms of CO2 conversion in a nanosecond repetitively pulsed (NRP) discharge, by means of a chemical kinetics model. The calculated conversions and energy efficiencies are in reasonable agreement with experimental results over a wide range of specific energy input values, and the same applies to the evolution of gas temperature and CO2 conversion as a function of time in the afterglow, indicating that our model provides a realistic picture of the underlying mechanisms in the NRP discharge and can be used to identify its limitations and thus to suggest further improvements. Our model predicts that vibrational excitation is very important in the NRP discharge, explaining why this type of plasma yields energy-efficient CO2 conversion. A significant part of the CO2 dissociation occurs by electronic excitation from the lower vibrational levels toward repulsive electronic states, thus resulting in dissociation. However, vibration−translation (VT) relaxation (depopulating the higher vibrational levels) and CO + O recombination (CO + O + M → CO2 + M), as well as mixing of the converted gas with fresh gas entering the plasma in between the pulses, are limiting factors for the conversion and energy efficiency. Our model predicts that extra cooling, slowing down the rate of VT relaxation and of the above recombination reaction, thus enhancing the contribution of the highest vibrational levels to the overall CO2 dissociation, can further improve the performance of the NRP discharge for energy-efficient CO2 conversion.

Abstract: Plasma technology is gaining increasing interest for CO2 conversion. The gas temperature in (and after) the plasma reactor largely affects the performance. Therefore, we examine the effect of cooling/quenching, during and after the plasma, on the CO2 conversion and energy efficiency, for typical “warm” plasmas, by means of chemical kinetics modeling. For plasmas at low specific energy input (SEI ∼ 0.5 eV/molecule), it is best to quench at the plasma end, while for high-SEI plasmas (SEI ∼ 4 eV/molecule), quenching at maximum conversion is better. For low-SEI plasmas, quenching can even increase the conversion beyond the dissociation in the plasma, known as superideal quenching. To better understand the effects of quenching at different plasma conditions, we study the dissociation and recombination rates, as well as the vibrational distribution functions (VDFs) of CO2, CO, and O2. When a high vibrational−translational (VT) nonequilibrium exists at the moment of quenching, the dissociation and recombination reaction rates both increase. Depending on the conversion degree at the moment of quenching, this can lead to a net increase or decrease of CO2 conversion. In general, however, and certainly for equilibrium plasmas at high temperature, quenching after the plasma helps prevent recombination reactions and clearly enhances the final CO2 conversion. We also investigate the effect of different quenching cooling rates on the CO2 conversion and energy efficiency. Finally, we compare plasma-based conversion to purely thermal conversion. For warm plasmas with typical temperatures of 3000−4000 K, the conversion is roughly thermal.

Abstract: Supersonic flows provide a high thermodynamic
nonequilibrium, which is crucial for energy-efficient conversion of
CO 2 in microwave plasmas and are therefore of great interest.
However, the effect of the flow on the chemical reactions is poorly
understood. In this work, we present a combined flow and plasma
chemical kinetics model of a microwave CO 2 plasma in a Laval
nozzle setup. The effects of the flow field on the different dissociation
and recombination mechanisms, the vibrational distribution, and the
vibrational transfer mechanism are discussed. In addition, the effect
of experimental parameters, like position of power deposition, outlet
pressure, and specific energy input, on the CO 2 conversion and
energy efficiency is examined. The short residence time of the gas in
the plasma region, the shockwave, and the maximum critical heat,
and thus power, that can be added to the flow to avoid thermal
choking are the main obstacles to reaching high energy efficiencies.

Abstract: We describe the plasma chemistry in a helium flowing atmospheric pressure afterglow (FAPA) used for analytical spectrometry, by means of a quasione-dimensional (1D) plasma chemical kinetics model. We study the effect of typical impurities present in the feed gas, as well as the afterglow in ambient humid air. The model provides the species density profiles in the discharge and afterglow regions and the chemical pathways. We demonstrate that H, N, and O atoms are formed in the discharge region, while the dominant reactive neutral species in the afterglow are O3 and NO. He* and He2* are responsible for Penning ionization of O2, N2, H2O, H2, and N, and especially O and H atoms. Besides, He2+ also contributes to ionization of N2, O2, H2O, and O through charge transfer reactions. From the pool of ions created in the discharge, NO+ and (H2O)3H+ are the dominant ions in the afterglow. Moreover, negatively charged clusters, such as NO3H2O− and NO2H2O−, are formed and their pathway is discussed as well. Our model predictions are in line with earlier observations in the literature about the important reagent ions and provide a comprehensive overview of the underlying pathways. The model explains in detail why helium provides a high analytical sensitivity because of high reagent ion formation by both Penning ionization and charge transfer. Such insights are very valuable for improving the analytical performance of this (and other) ambient desorption/ionization source(s).

Abstract: Plasma catalysis is emerging for plasma-assisted gas conversion processes. However, the underlying mechanisms of plasma catalysis are poorly understood. In this work, we present a 1D heterogeneous catalysis model with axial dispersion (i.e., accounting for back-mixing and molecular diffusion of fluid elements in the process stream in the axial direction), for plasma-catalytic NO production from N2/O2 mixtures. We investigate the concentration and reaction rates of each species formed as a function of time and position across the catalyst, in order to determine the underlying mechanisms. To obtain insights into how the performance of the process can be further improved, we also study how changes in the postplasma gas flow composition entering the catalyst bed and in the operation conditions of the catalytic stage affect the performance of NO production.

Abstract: Plasma catalysis is emerging for plasma-assisted gas conversion
processes. However, the underlying mechanisms of plasma catalysis are poorly
understood. In this work, we present a 1D heterogeneous catalysis model with axial
dispersion (i.e., accounting for back-mixing and molecular diffusion of fluid elements in
the process stream in the axial direction), for plasma-catalytic NO production from
N2/O2 mixtures. We investigate the concentration and reaction rates of each species
formed as a function of time and position across the catalyst, in order to determine the
underlying mechanisms. To obtain insights into how the performance of the process
can be further improved, we also study how changes in the postplasma gas flow
composition entering the catalyst bed and in the operation conditions of the catalytic
stage affect the performance of NO production.

Abstract: The continuous miniaturization of nanodevices, such as transistors, solar cells, and optical fibers, requires the controlled synthesis of (ultra)thin gate oxides (<10 nm), including Si gate-oxide (SiO2) with high quality at the atomic scale. Traditional thermal growth of SiO2 on planar Si surfaces, however, does not allow one to obtain such ultrathin oxide due to either the high oxygen diffusivity at high temperature or the very low sticking ability of incident oxygen at low temperature. Two recent techniques, both operative at low (room) temperature, have been put forward to overcome these obstacles: (i) hyperthermal oxidation of planar Si surfaces and (ii) thermal or plasma-assisted oxidation of nonplanar Si surfaces, including Si nanowires (SiNWs). These nanooxidation processes are, however, often difficult to study experimentally, due to the key intermediate processes taking place on the nanosecond time scale.
In this Account, these Si nano-oxidation techniques are discussed from a computational point of view and compared to both hyperthermal and thermal oxidation experiments, as well as to well-known models of thermal oxidation, including the Deal−Grove, Cabrera−Mott, and Kao models and several alternative mechanisms. In our studies, we use reactive molecular dynamics (MD) and hybrid MD/Monte Carlo simulation techniques, applying the Reax force field. The incident energy of oxygen species is chosen in the range of 1−5 eV in hyperthermal oxidation of planar Si surfaces in order to prevent energy-induced damage. It turns out that hyperthermal growth allows for two growth modes, where the ultrathin oxide thickness depends on either (1) only the kinetic energy of the incident oxygen species at a growth temperature below Ttrans = 600 K, or (2) both the incident energy and the growth temperature at a growth temperature above Ttrans. These modes are specific to such ultrathin oxides, and are not observed in traditional thermal oxidation, nor theoretically considered by already existing models. In the case of thermal or plasma-assisted oxidation of small Si nanowires, on the other hand, the thickness of the ultrathin oxide is a function of the growth temperature and the nanowire diameter. Below Ttrans, which varies with the nanowire diameter, partially oxidized SiNW are formed, whereas complete oxidation to a SiO2 nanowire occurs only above Ttrans. In both nano-oxidation processes at lower temperature (T < Ttrans), final sandwich c-Si|SiOx|a-SiO2 structures are obtained due to a competition between overcoming the energy barrier to penetrate into Si subsurface layers and the compressive stress (∼2−3 GPa) at the Si crystal/oxide interface. The overall atomic-simulation results strongly indicate that the thickness of the intermediate SiOx (x < 2) region is very limited (∼0.5 nm) and constant irrespective of oxidation parameters. Thus, control over the ultrathin SiO2 thickness with good quality is indeed possible by accurately tuning the oxidant energy, oxidation temperature and surface curvature.
In general, we discuss and put in perspective these two oxidation mechanisms for obtaining controllable ultrathin gate-oxide films, offering a new route toward the fabrication of nanodevices via selective nano-oxidation.

Abstract: Lipid oxidation is associated with several inflammatory and neurodegenerative diseases, but many questions to unravel its effects on biomembranes are still open due to the complexity of the topic. For instance, recent studies indicated that phase-separated domains can have a significant effect on membrane function. It is reported that domain interfaces are “hot spots” for pore formation, but the underlying mechanisms and the effect of oxidation-induced phase separation on membranes remain elusive. Thus, to evaluate the permeability of the membrane coexisting of liquid-ordered (Lo) and liquid-disordered (Ld) domains, we performed atomistic molecular dynamics simulations. Specifically, we studied the membrane permeability of nonoxidized or oxidized homogeneous membranes (single-phase) and at the Lo/Ld domain interfaces of heterogeneous membranes, where the Ld domain is composed of either oxidized or nonoxidized lipids. Our simulation results reveal that the addition of only 1.5% of lipid aldehyde molecules at the Lo/Ld domain interfaces of heterogeneous membranes increases the membrane permeability, whereas their addition at homogeneous membranes does not have any effect. This study is of interest for a better understanding of cancer treatment methods based on oxidative stress (causing among others lipid oxidation), such as plasma medicine and photodynamic therapy.

Abstract: We perform computer simulations supported by experiments to investigate the oxidation of an important signaling protein, that is, human epidermal growth factor (hEGF), caused by cold atmospheric plasma (CAP) treatment. Specifically, we study the conformational changes of hEGF with different degrees of oxidation, to mimic short and long CAP treatment times. Our results indicate that the oxidized structures become more flexible, due to their conformational changes and breakage of the disulfide bonds, especially at higher oxidation degrees. MM/GBSA calculations reveal that an increasing oxidation level leads to a lower binding free energy of hEGF with its receptor. These results help to understand the fundamentals of the use of CAP for wound healing versus cancer treatment at short and longer treatment times.

Abstract: Plasma has been studied for several years to convert CO2 into value-added products. If CO2 could be converted in the presence of H2O as a cheap H-source for making syngas and oxygenates, it would mimic natural photosynthesis. However, CO2/H2O plasmas have not yet been extensively studied, not by experiments, and certainly not computationally. Therefore, we present here a kinetic modelling study to obtain a greater understanding of the vibrational kinetics of a CO2/H2O microwave plasma. For this purpose, we first created an electron impact cross section set for H2O, using a swarm-derived method. We added the new cross section set and CO2/H2O-related chemistry to a pure CO2 model. While it was expected that H2O addition mainly causes quenching of the CO2 asymmetric mode vibrational levels due to the additional CO2/H2O vibrational-translational relaxation, our model shows that the modifications in the vibrational kinetics are mainly induced by the strong electron dissociative attachment to H2O molecules, causing a reduction in electron density, and the corresponding changes in the input of energy into the CO2 vibrational levels by electron impact processes.

Abstract: A pulsed microwave surfaguide discharge operating at 2.45 GHz was used for the conversion of molecular nitrogen into valuable compounds in several gas mixtures: N2 :O2 , N2 :O2 :CO2 and N2 :CO2 . The ro-vibrational absorption bands of the molecular species were monitored by a Fourier transform infrared apparatus in the post-discharge region in order to evaluate the relative number density of species, specifically NO production. The effects of specific energy input, pulse frequency, gas flow fraction, gas admixture and gas flow rate were studied for better understanding and optimization of the NO production yield and the corresponding energy cost (EC). By both the experiment and modelling, a highest NO yield is obtained at N2 :O2 (1:1) gas ratio in N2 :O2 mixture. The NO yield reveals a small growth followed by saturation when pulse repetition frequency increases. The energy efficiency start decreasing after the energy input reaches about 5 eV/molec, whereas the NO yield rises steadily at the same time. The lowest EC of about 8 MJ mol−1 corresponding to the yield and the energy efficiency of about 7% and 1% are found, respectively, in an optimum discharge condition in our case.

Abstract: In this study, we investigated the diffusion of H-atoms to the subsurface and their further diffusion into the bulk of a Ni(111) crystal by means of density functional theory calculations in the context of thermal and plasma-assisted catalysis. The H-atoms at the surface can originate from the dissociative adsorption of H2 or CH4 molecules, determining the surface H-coverage. When a threshold H-coverage is passed, corresponding to 1.00 ML for the crystalline Ni(111) surface, the surface-bound H-atoms start to diffuse to the subsurface. A similar threshold coverage is observed for the interstitial H-coverage. Once the interstitial sites are filled up with a coverage above 1.00 ML of H, dissolution of interstitial H-atoms to the layer below the interstitial sites will be initiated. Hence, by applying a high pressure or inducing a reactive plasma and high temperature, increasing the H-flux to the surface, a large amount of hydrogen can diffuse in a crystalline metal like Ni and can be absorbed. The formation of metal hydride may modify the entire reaction kinetics of the system. Equivalently, the H-atoms in the bulk can easily go back to the surface and release a large amount of heat. In a plasma process, H-atoms are formed in the plasma, and therefore the energy barrier for dissociative adsorption is dismissed, thus allowing achievement of the threshold coverage without applying a high pressure as in a thermal process. As a result, depending on the crystal plane and type of metal, a large number of H-atoms can be dissolved (absorbed) in the metal catalyst, explaining the high efficiency of plasma-assisted catalytic reactions. Here, the mechanism of H-dissolution is established as a chemical identifier for the investigation of the reaction kinetics of a chemical process.

Abstract: The hydrogenation of carbon dioxide in a dielectric barrier discharge plasma is studied with a one-dimensional fluid model. The spatially averaged densities of the most important end products formed in the CO2/H2 mixture are determined as a function of the initial gas mixing ratio. CO and H2O are found to be present at the highest densities and to a lower content also CH4, C2H6, CH2O, CH3OH, O2, and some other higher hydrocarbons and oxygenates. The main underlying reaction
pathways for the conversion of the inlet gases and the formation of CO, CH4, CH2O, and CH3OH are pointed out for various gas mixing ratios. The CO2 conversion and the production of value added products is found to be quite low, also in comparison to a CO2/CH4 mixture, and this can be explained by the model.

Abstract: An Ar plasma sustained by a surfaguide wave launcher is investigated at intermediate pressure (200–2667 Pa). Two 2D self-consistent models (quasi-neutral and plasma bulk-sheath) are developed and benchmarked. The complete set of electromagnetic and fluid equations and the boundary conditions are presented. The transformation of fluid equations from a local reference frame, that is, moving with plasma or when the gas flow is zero, to a laboratory reference frame, that is,
accounting for the gas flow, is discussed. The pressure range is extended down to 80 Pa by experimental measurements. The electron temperature decreases with pressure. The electron density depends linearly on power, and changes its behavior with pressure depending on the product of pressure and radial plasma size.