Abstract: There is currently a growing interest in the realisation and optimization of hybrid plasma/catalyst systems for a multitude of applications, ranging from nanotechnology to environmental chemistry. In spite of this interest, there is, however, a lack in fundamental understanding of the underlying processes in such systems. While a lot of experimental research is already being carried out to gain this understanding, only recently the first simulations have appeared in the literature. In this contribution, an overview is presented on atomic scale simulations of plasma catalytic processes as carried out in our group. In particular, this contribution focusses on plasma-assisted catalyzed carbon nanostructure growth, and plasma catalysis for greenhouse gas conversion. Attention is paid to what can routinely be done, and where challenges persist.

Abstract: Amorphous semiconductors are usually characterized by a low charge carrier mobility, essentially related to their lack of long-range order. The development of such material with higher charge carrier mobility is hence challenging. Part of the issue comes from the difficulty encountered by first-principles simulations to evaluate concepts such as the electron effective mass for disordered systems since the absence of periodicity induced by the disorder precludes the use of common concepts derived from condensed matter physics. In this paper, we propose a methodology based on first-principles simulations that partially solves this problem, by quantifying the degree of delocalization of a wave function and of the connectivity between the atomic sites within this electronic state. We validate the robustness of the proposed formalism on crystalline and molecular systems and extend the insights gained to disordered/amorphous InGaZnO4 and Si. We also explore the properties of p-type oxide semiconductor candidates recently reported to have a low effective mass in their crystalline phases [G. Hautier et al., Nat. Commun. 4, 2292 (2013)]. Although in their amorphous phase none of the candidates present a valence band with delocalization properties matching those found in the conduction band of amorphous InGaZnO4, three of the seven analyzed materials show some potential. The most promising candidate, K2Sn2O3, is expected to possess in its amorphous phase a slightly higher hole mobility than the electron mobility in amorphous silicon.

Abstract: Experiments have demonstrated the potential selective anticancer capacity of cold atmospheric plasmas (CAPs), but the underlying mechanisms remain unclear. Using computer simulations, we try to shed light on the mechanism of selectivity, based on aquaporins (AQPs), i.e. transmembrane protein channels transferring external H 2 O 2 and other reactive oxygen species, created e.g. by CAPs, to the cell interior. Specifically, we perform molecular dynamics simulations for the permeation of H 2 O 2 through AQP1 (one of the members of the AQP family) and the palmitoyl-oleoyl-phosphatidylcholine (POPC) phospholipid bilayer (PLB). The free energy barrier of H 2 O 2 across AQP1 is lower than for the POPC PLB, while the permeability coefficient, calculated using the free energy and diffusion rate profiles, is two orders of magnitude higher. This indicates that the delivery of H 2 O 2 into the cell interior should be through AQP. Our study gives a better insight into the role of AQPs in the selectivity of CAPs for treating cancer cells.

Abstract: The interaction between cold atmospheric pressure plasma and liquids is receiving increasing attention for various applications. In particular, the use of plasma-treated liquids (PTL) for biomedical applications is of growing importance, in particular for sterilization and cancer treatment. However, insight into the
underlying mechanisms of plasma–liquid interactions is still scarce. Here, we present a 2D fluid dynamics model for the interaction between a plasma jet and liquid water. Our results indicate that the formed reactive species originate from either the gas phase (with further solvation) or are formed at the liquid interface. A clear increase in the aqueous density of H2O2, HNO2/NO2- and NO3-
is observed as a function of time, while the densities of O3, HO2/O2- and ONOOH/ONOO- are found to quickly reach a maximum due to chemical reactions in solution. The trends observed in our model correlate well with experimental observations from the literature.

Abstract: We perform molecular dynamics simulations to study the flip-flop motion of phosphatidylserine (PS) across the plasma membrane upon increasing oxidation degree of the membrane. Our computational results show that an increase of the oxidation degree in the lipids leads to a decrease of the free energy barrier for translocation of PS through the membrane. In other words, oxidation of the lipids facilitates PS flip-flop motion across the membrane, because in native phospholipid bilayers this is only a “rare event” due to the high energy barriers for the translocation of PS. The present study provides an atomic-scale insight into the mechanisms of the PS flip-flop upon oxidation of lipids, as produced for example by cold atmospheric plasma, in living cells.

Abstract: Epitaxially grown ultra-thin Si layers are often used to passivate Ge surfaces in the high-k gate module of (strained) Ge FinFET devices. We use Si4H10 as Si precursor as it enables epitaxial Si growth at temperatures down to 330 degrees C. C-V characteristics of blanket capacitors made on Ge virtual substrates point to the presence of an optimal Si thickness. In case of compressively strained Ge fin structures, the Si growth results in non-uniform and high strain levels in the strained Ge fin. These strain levels have been calculated for different shapes of the Ge fin and in function of the grown Si thickness. The high strain is the driving force for potential (unwanted) Ge surface reflow during the Si deposition. The Ge surface reflow is strongly affected by the strength of the H-passivation during Si-capping and can be avoided by carefully selected process conditions.

Abstract: A number of recent computational material design studies based on density functional theory (DFT) calculations have put forward a new class of materials with electrically switchable chemical characteristics that can be exploited in the development of tunable gas storage and electrocatalytic applications. We find systematic flaws in almost every computational study of gas adsorption on polarized or charged surfaces, stemming from an improper and unreproducible treatment of periodicity, leading to very large errors of up to 3 eV in some cases. Two simple corrective procedures that lead to consistent results are proposed, constituting a crucial course correction to the research in the field.

Abstract: In this study, we performed density functional theory (DFT) calculations to investigate different reaction mechanisms of CO oxidation catalyzed by the Si atom embedded defective BC2N nanostructures as well as the analysis of the structural and electronic properties. The structures of all the complexes are optimized and characterized by frequency calculations at the M062X/6-31G* computational level. Also, The electronic structures and thermodynamic parameters of adsorbed CO and O-2 molecules over Si-doped BC2N nanostructures are examined in detail. Moreover, to investigate the curvature effect on the CO oxidation reaction, all the adsorption and CO oxidation reactions on a finite-sized armchair (6,6) Si-BC2NNT are also studied. Our results indicate that there can be two possible pathways for the CO oxidation with O-2 molecule: O-2(g) + CO(g) -> O-2(ads) + CO(ads) -> CO2(g) + O-(ads) and O-(ads) + CO(g) -> CO2(g). The first reaction proceeds via the Langmuir-Hinshelwood (LH) mechanism while the second goes through the Eley-Rideal (ER) mechanism. On the other hand, by increasing the tube diameter, the energy barrier increases due to the strong adsorption energy of the O-2 molecule which is related to its dissociation over the tube surface. Our calculations indicate that the two step energy barrier of the oxidation reaction over Si-BC2NNS is less than that over the Si-BC2NNT. Hence, Si-BC2NNS may serve as an efficient and highly activated substrate to CO oxidation rather than (4,4) Si-BC2NNT. (C) 2018 Elsevier B.V. All rights reserved.

Abstract: Microwave (MW) plasmas represent a promising solution for efficient CO2 dissociation. MW discharges are also very versatile and can be sustained at various pressure and gas flow regimes. To identify the most favorable conditions for the further scale-up of the CO2 decomposition reaction, a MW plasma reactor operating in pure CO2 in a wide pressure range (200 mbar–1 bar) is studied. Three different gas flow configurations are explored: a direct, reverse and a vortex regime. The CO2 conversion and energy efficiency drop almost linearly with increasing pressure, regardless of the gas flow regime. The results obtained in the direct flow configuration underline the importance of post-discharge cooling, as the exhaust of the MW plasma reactor in this regime expanded into the vacuum chamber without additional quenching. As a result, this system yields exhaust temperatures of up to 1000 K, which explains the lowest conversion (∼3.5% at 200 mbar and 2% at 1 bar). A post-discharge cooling step is introduced for the reverse gas inlet regime and allows the highest conversion to be achieved (∼38% at 200 mbar and 6.2% at 1 bar, with energy efficiencies of 23% and 3.7%). Finally, a tangential gas inlet is utilized in the vortex configuration to generate a swirl flow pattern. This results in the generation of a stable discharge in a broader range of CO2 flows (15–30 SLM) and the highest energy efficiencies obtained in this study (∼25% at 300 mbar and ∼13% at 1 bar, at conversions of 21% and 12%). The experimental results are complemented with computational fluid dynamics simulations and with the analysis of the latest literature to identify the further research directions.

Abstract: Catalyst preparation with plasmas is increasingly attracting interest. A plasma is a partially ionized gas, consisting of electrons, ions, molecules, radicals, photons, and excited species, which are all active species for catalyst preparation and treatment. Under the influence of plasma, nucleation and crystal growth in catalyst preparation can be very different from those in the conventional thermal approach. Some thermodynamically unfavorable reactions can easily take place with plasmas. Compounds such as sulfides, nitrides, and phosphides that are produced under harsh conditions can be synthesized by plasma under mild conditions. Plasmas can produce catalysts with smaller particle sizes and controllable structure. Plasma is also a facile tool for reduction, oxidation, doping, etching, coating, alloy formation, surface treatment, and surface cleaning in a simple and direct way. A rapid and convenient plasma template removal has thus been established for zeolite synthesis. It can operate at room temperature and allows the catalyst preparation on temperature-sensitive supporting materials. Plasma is typically effective for the production of various catalysts on metallic substrates. In addition, plasma-prepared transition-metal catalysts show enhanced low-temperature activity with improved stability. This provides a useful model catalyst for further improvement of industrial catalysts. In this review, we aim to summarize the recent advances in catalyst preparation with plasmas. The present understanding of plasma-based catalyst preparation is discussed. The challenges and future development are addressed.

Abstract: Because of the unique properties of plasma technology, its use in gas conversion applications is gaining significant interest around the globe. Plasma-based CO2 and CH4 conversion has become a major research area. Many investigations have already been performed regarding the single-component gases, that is, CO2 splitting and CH4 reforming, as well as for two-component mixtures, that is, dry reforming of methane
(CO2/CH4), partial oxidation of methane (CH4/O2), artificial photosynthesis (CO2/H2O), CO2 hydrogenation (CO2/H2), and even first steps toward the influence of N2 impurities have been taken, that is, CO2/N2 and CH4/N2. In this Feature Article we briefly discuss the advances made in literature for these different steps from a plasma chemistry modeling point of view. Subsequently, we present a comprehensive plasma chemistry set, combining the knowledge gathered in this field so far and supported with extensive experimental data. This set can be used for chemical kinetics plasma modeling for all possible combinations of CO2, CH4, N2, O2, and H2O to investigate the bigger picture of the underlying plasmachemical pathways for these mixtures in a dielectric barrier discharge plasma. This is extremely valuable
for the optimization of existing plasma-based CO2 conversion and CH4 reforming processes as well as for investigating the influence of N2, O2, and H2O on these processes and even to support plasma-based multireforming processes.

Abstract: Reactive oxygen and nitrogen species released by cold physical plasma are being proposed as effectors in various clinical conditions connected to inflammatory processes. As these plasmas can be tailored in a wide range, models to compare and control their biochemical footprint are desired to infer on the molecular mechanisms underlying the observed effects and to enable the discrimination between different plasma sources. Here, an improved model to trace short-lived reactive species is presented. Using FTIR, high-resolution mass spectrometry, and molecular dynamics computational simulation, covalent modifications of cysteine treated with different plasmas were deciphered and the respective product pattern used to generate a fingerprint of each plasma source. Such, our experimental model allows a fast and reliable grading of the chemical potential of plasmas used for medical purposes. Major reaction products were identified to be cysteine sulfonic acid, cystine, and cysteine fragments. Less abundant products, such as oxidized cystine derivatives or S-nitrosylated cysteines, were unique to different plasma sources or operating conditions. The data collected point at hydroxyl radicals, atomic O, and singlet oxygen as major contributing species that enable an impact on cellular thiol groups when applying cold plasma in vitro or in vivo.

Abstract: In this work, a chemical kinetics study on methane activation for hydrogen production in a warm plasma, i.e., N-2 rotating gliding arc (RGA), was performed for the first time to get new insights into the underlying reaction mechanisms and pathways. A zero-dimensional chemical kinetics model was developed, which showed a good agreement with the experimental results in terms of the conversion of CH4 and product selectivities, allowing us to get a better understanding of the relative significance of various important species and their related reactions to the formation and loss of CH4, H-2, and C2H2 etc. An overall reaction scheme was obtained to provide a realistic picture of the plasma chemistry. The results reveal that the electrons and excited nitrogen species (mainly N-2(A)) play a dominant role in the initial dissociation of CH4. However, the H atom induced reaction CH4+ H -> CH3+ H-2, which has an enhanced reaction rate due to the high gas temperature (over 1200 K), is the major contributor to both the conversion of CH4 and H-2 production, with its relative contributions of > 90% and > 85%, respectively, when only considering the forward reactions. The coexistence and interaction of thermochemical and plasma chemical processes in the rotating gliding arc warm plasma significantly enhance the process performance. The formation of C-2 hydrocarbons follows a nearly one-way path of C2H6 -> C2H4 -> C2H2, explaining why the selectivities of C-2 products decreased in the order of C2H2 > C2H4 > C2H6.

Abstract: Based on a rational classification of defects in amorphous materials, we propose a simplified model to describe intrinsic defects and hydrogen impurities in amorphous indium gallium zinc oxide (a-IGZO). The proposed approach consists of organizing defects into two categories: point defects, generating structural anomalies such as metal-metal or oxygen-oxygen bonds, and defects emerging from changes in the material stoichiometry, such as vacancies and interstitial atoms. Based on first-principles simulations, it is argued that the defects originating from the second group always act as perfect donors or perfect acceptors. This classification simplifies and rationalizes the nature of defects in amorphous phases. In a-IGZO, the most important point defects are metal-metal bonds (or small metal clusters) and peroxides (O-O single bonds). Electrons are captured by metal-metal bonds and released by the formation of peroxides. The presence of hydrogen can lead to two additional types of defects: metal-hydrogen defects, acting as acceptors, and oxygen-hydrogen defects, acting as donors. The impact of these defects is linked to different instabilities observed in a-IGZO. Specifically, the diffusion of hydrogen and oxygen is connected to positive-and negative-bias stresses, while negative-bias illumination stress originates from the formation of peroxides.

Abstract: Plasma generation inside catalyst pores is of utmost importance for plasma catalysis, as the existence of plasma species inside the pores affects the active surface area of the catalyst available to the plasma species for catalytic reactions. In this paper, the electric field enhancement, and thus the plasma production inside catalyst pores with different pore shapes is studied with a two-dimensional fluid model. The results indicate that the electric field will be significantly enhanced near tip-like structures. In a conical pore with small opening, the strongest electric field appears at the opening and bottom corners of the pore, giving rise to a prominent ionization rate throughout the pore. For a cylindrical pore, the electric field is only enhanced at the bottom corners of the pore, with lower absolute value, and thus the ionization rate inside the pore is only slightly enhanced. Finally, in a conical pore with large opening, the electric field is characterized by a maximum at the bottom of the pore, yielding a similar behavior for the ionization rate. These results demonstrate that the shape of the pore has a significantly influence on the electric field enhancement, and thus modifies the plasma properties.