Abstract: The resistance of a TiO2 anatase (001) surface to coke formation was studied in the context of dry reforming of methane using density functional theory (DFT) calculations. As carbon atoms act as precursors for coke formation, the resistance to coke formation can be measured by the carbon coverage of the surface. This is related to the stability of different CHx (x = 0−3) species and their rate of hydrogenation and dehydrogenation on the TiO2 surface. Therefore, we studied the reaction mechanisms and their corresponding rates as a function of the temperature for the dehydrogenation of the species on the surface. We found that the stabilities of C and CH are significantly lower than those of CH3 and CH2. The hydrogenation rates of the different species are significantly higher than the dehydrogenation rates in a temperature range of 300−1000 K. Furthermore, we found that dehydrogenation of CH3, CH2, and CH will only occur at appreciable rates starting from 600, 900, and 900 K, respectively. On the basis of these results, it is clear that the anatase (001) surface has a high coke resistance, and it is thus not likely that the surface will become poisoned by coke during dry reforming of methane. As the rate limiting step in dry reforming is the dissociative adsorption of CH4, we studied an alternative approach to thermal catalysis. We found that the temperature threshold for dry reforming is at least 700 K. This threshold temperature may be lowered by the use of plasma-catalysis, where the appreciable rates of adsorption of plasma-generated CHx radicals result in bypassing the rate limiting step of the reaction.

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: Nonequilibrium plasmas offer a pathway for energy-efficient CO2 conversion through vibrationally induced dissociation. However, the efficiency of this pathway is limited by a rise in gas temperature, which increases vibrational−translational (VT) relaxation and quenches the vibrational levels. Therefore, we investigate here the effect of plasma pulsing on the VT nonequilibrium and on the CO2 conversion by means of a zerodimensional chemical kinetics model, with self-consistent gas temperature calculation. Specifically, we show that higher energy efficiencies can be reached by correctly tuning the plasma pulse and interpulse times. The ideal plasma pulse time corresponds to the time needed to reach the highest vibrational temperature. In addition, the highest energy efficiencies are obtained with long interpulse times, that is, ≥0.1 s, in which the gas temperature can entirely drop to room temperature. Furthermore, additional cooling of the reactor walls can give higher energy efficiencies at shorter interpulse times of 1 ms. Finally, our model shows that plasma pulsing can significantly improve the energy efficiency at low reduced electric fields (50 and 100 Td, typical for microwave and gliding arc plasmas) and intermediate ionization degrees (5 × 10−7 and 10−6).

Abstract: Molecular dynamics simulations were used to investigate the size-dependent melting mechanism of nickel nanoclusters of various sizes. The melting process was monitored by the caloric curve, the overall cluster Lindemann index, and the atomic Lindemann index. Size-dependent melting temperatures were determined, and the correct linear dependence on inverse diameter was recovered. We found that the melting mechanism gradually changes from dynamic coexistence melting to surface melting with increasing cluster size. These findings are of importance in better understanding carbon nanotube growth by catalytic chemical vapor deposition as the phase state of the catalyst nanoparticle codetermines the growth mechanism.

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: Continued miniaturization of metal-oxide-semiconductor field-effect transistors (MOSFETs) requires an ever-decreasing thickness of the gate oxide. The structure of ultrathin silicon oxide films, however, critically depends on the oxidation mechanism. Using reactive atomistic simulations, we here demonstrate how the oxidation mechanism in hyperthermal oxidation of such structures may be controlled by the oxidation temperature and the oxidant energy. Specifically, we study the interaction of hyperthermal oxygen with energies of 15 eV with thin SiOx (x ≤ 2) films with a native oxide thickness of about 10 Å. We analyze the oxygen penetration depth probability and compare with results of the hyperthermal oxidation of a bare Si(100){2 × 1} (c-Si) surface. The temperature-dependent oxidation mechanisms are discussed in detail. Our results demonstrate that, at low (i.e., room) temperature, the penetrated oxygen mostly resides in the oxide region rather than at the SiOx|c-Si interface. However, at higher temperatures, starting at around 700 K, oxygen atoms are found to penetrate and to diffuse through the oxide layer followed by reaction at the c-Si boundary. We demonstrate that hyperthermal oxidation resembles thermal oxidation, which can be described by the DealGrove model at high temperatures. Furthermore, defect creation mechanisms that occur during the oxidation process are also analyzed. This study is useful for the fabrication of ultrathin silicon oxide gate oxides for metal-oxide-semiconductor devices as it links parameters that can be straightforwardly controlled in experiment (oxygen temperature, velocity) with the silicon oxide structure.

Abstract: Nonthermal atmospheric pressure plasmas are gaining increasing attention for biomedical applications. However, very little fundamental information on the interaction mechanisms between the plasma species and biological cells is currently available. We investigate the interaction of important plasma species, such as OH, H2O2, O, O3, as well as O2 and H2O, with bacterial peptidoglycan by means of reactive molecular dynamics simulations, aiming for a better understanding of plasma disinfection. Our results show that OH, O, O3, and H2O2 can break structurally important bonds of peptidoglycan (i.e., CO, CN, or CC bonds), which consequently leads to the destruction of the bacterial cell wall. The mechanisms behind these breakups are, however, dependent on the impinging plasma species, and this also determines the effectiveness of the cell wall destruction.

Abstract: A one-dimensional fluid model for a dielectric barrier discharge in CH4/O2 and CH4/CO2 gas mixtures is developed. The model describes the gas-phase chemistry for partial oxidation and for dry reforming of methane. The spatially averaged densities of the various plasma species are presented as a function of time and initial gas mixing ratio. Besides, the conversion of the inlet gases and the selectivities of the reaction products are calculated. Syngas, higher hydrocarbons, and higher oxygenates are typically found to be important reaction products. Furthermore, the main underlying reaction pathways for the formation of syngas, methanol, formaldehyde, and other higher oxygenates are determined.

Abstract: Aquaporins (AQPs) are transmembrane proteins that conduct not only water molecules across the cell membrane but also other solutes, such as reactive oxygen and nitrogen species (RONS), produced (among others) by cold atmospheric plasma (CAP). These RONS may induce oxidative stress in the cell interior, which plays a role in cancer treatment. The underlying mechanisms of the transport of RONS across AQPs, however, still remain obscure. We apply molecular dynamics simulations to investigate the permeation of both hydrophilic (H₂O₂and OH) and hydrophobic (NO₂and NO) RONS through AQP1. Our simulations show that these RONS can all penetrate across the pores of AQP1. The permeation free energy barrier of OH and NO is lower than that of H₂O₂and NO₂, indicating that these radicals may have easier access to the pore interior and interact with the amino acid residues of AQP1. We also study the effect of RONS-induced oxidation of both the phospholipids and AQP1 (i.e., sulfenylation of Cys₁₉₁) on the transport of the above-mentioned RONS across AQP1. Both lipid and protein oxidation seem to slightly increase the free energy barrier for H₂O₂and NO₂permeation, while for OH and NO, we do not observe a strong effect of oxidation. The simulation results help to gain insight in the underlying mechanisms of the noticeable rise of CAP-induced RONS in cancer cells, thereby improving our understanding on the role of AQPs in the selective anticancer capacity of CAP.

Abstract: Plasma catalysis is promising for various environmental, energy and chemical synthesis applications, but the underlying mechanisms are far from understood. Modeling can help to obtain a better insight in these mechanisms. Some burning questions relate to the plasma behavior inside packed bed reactors and whether plasma can penetrate into catalyst pores. In this paper, we try to provide answers to these questions, by means of both fluid modeling and particle-in-cell/Monte Carlo collision simulations. We present a short overview of recent findings obtained in our group by means of modeling, i.e., the enhanced electric field near the contact points and the streamer propagation through the packing in packed bed reactors, as well as the plasma behavior in catalyst pores, to determine the minimum pore size in which plasma streamers can penetrate.

Abstract: Zeogrids and Zeotiles are hierarchical materials built from assembled MFI zeolite precursor units. Permanent secondary porosity in these materials is obtained through self assembly of nanoparticles encountered in MFI zeolite synthesis in the presence of supramolecular templates. Hereon, the aggregated species are termed nanoslabs. Zeogrids are layered materials with lateral spacings between nanoslabs creating galleries qualifying as supermicropores. Zeotiles present a diversity of tridimensional nanoslab assemblies with mesopores. Zeotile-1, -4 and -6 are hexagonal mesostructures. Zeotile-1 has triangular and hexagonal channels; Zeotile-4 has hexagonal channels interconnected via slits. Zeotile-2 has a cubic structure with gyroid type mesoporosity. The behavior of Zeogrids and Zeotiles in adsorption, membrane and chromatographic separation and catalysis has been characterized and compared with zeolites and mesoporous materials derived from unstructured silica sources. Shape selectivity was detected via adsorption of n- and iso-alkanes. The mesoporosity of Zeotiles can be exploited in chromatographic separation of biomolecules. Zeotiles present attractive separation properties relevant to CO2 sequestration. Because of its facile synthesis procedure without hydrothermal steps Zeogrid is convenient for membrane synthesis. The performance of Zeogrid membrane in gas separation, nanofiltration and pervaporation is reported. In the Beckmann rearrangement of cyclohexanone oxime Zeogrids and Zeotiles display a catalytic activity characteristic of silicalite-1 zeolites. Introduction of acidity and redox catalytic activity can be achieved via incorporation of Al and Ti atoms in the nanoslabs during synthesis. Zeogrids are active in hydrocracking, catalytic cracking, alkylation and epoxidation reactions. Zeogrids and Zeotiles often behave differently from ordered mesoporous materials as well as from zeolites and present a valuable extension of the family of hierarchical silicate based materials.

Abstract: Background: Strong electric fields are knownto affect cell membrane permeability,which can be applied for therapeutic purposes, e.g., in cancer therapy. A synergistic enhancement of this effect may be accomplished by the presence of reactive oxygen species (ROS), as generated in cold atmospheric plasmas. Little is known about the synergy between lipid oxidation by ROS and the electric field, nor on howthis affects the cell membrane permeability.
Method: We here conduct molecular dynamics simulations to elucidate the dynamics of the permeation process under the influence of combined lipid oxidation and electroporation. A phospholipid bilayer (PLB), consisting of di-oleoyl-phosphatidylcholine molecules covered with water layers, is used as a model system for the plasma membrane.
Results and conclusions:Weshow howoxidation of the lipids in the PLB leads to an increase of the permeability of the bilayer to ROS, although the permeation free energy barriers still remain relatively high. More importantly, oxidation of the lipids results in a drop of the electric field threshold needed for pore formation (i.e., electroporation) in the PLB. The created pores in the membrane facilitate the penetration of reactive plasma species deep into the cell interior, eventually causing oxidative damage.
General significance: This study is of particular interest for plasma medicine, as plasma generates both ROS and electric fields, but it is also of more general interest for applications where strong electric fields and ROS both come into play.

Abstract: Plasma catalysis has emerged as a promising approach for driving thermodynamically unfavorable chemical
reactions. Nevertheless, comprehending the mechanisms involved remains a challenge, leading to uncertainty
about whether the optimal catalyst in plasma catalysis aligns with that in thermal catalysis. In this research, we
explore this question by studying monometallic catalysts (Fe, Co, Ni and Mo) and bimetallic catalysts (Fe-Co, Mo-
Co, Fe-Ni and Mo-Ni) in both thermal catalytic and plasma catalytic NH3 decomposition. Our findings reveal that
the Fe-Co bimetallic catalyst exhibits the highest activity in thermal catalysis, the Fe-Ni bimetallic catalyst
outperforms others in plasma catalysis, indicating a discrepancy between the optimal catalysts for the two
catalytic modes in NH3 decomposition. Comprehensive catalyst characterization, kinetic analysis, temperature
program surface reaction experiments and plasma diagnosis are employed to discuss the key factors influencing
NH3 decomposition performance.

Abstract: Thermal-catalytic gas processing is integral to many current industrial processes. Ever-increasing demands on conversion and energy efficiencies are a strong driving force for the development of alternative approaches. Similarly, synthesis of several functional materials (such as nanowires and nanotubes) demands special processing conditions. Plasma catalysis provides such an alternative, where the catalytic process is complemented by the use of plasmas that activate the source gas. This combination is often observed to result in a synergy between plasma and catalyst. This Review introduces the current state-of-the-art in plasma catalysis, including numerous examples where plasma catalysis has demonstrated its benefits or shows future potential, including CO2 conversion, hydrocarbon reforming, synthesis of nanomaterials, ammonia production, and abatement of toxic waste gases. The underlying mechanisms governing these applications, as resulting from the interaction between the plasma and the catalyst, render the process highly complex, and little is known about the factors leading to the often-observed synergy. This Review critically examines the catalytic mechanisms relevant to each specific application.