Abstract: Plasma catalysis is an emerging new technology for the electrification and downscaling of NH3 synthesis. Increasing attention is being paid to the optimization of plasma catalysis with respect to the plasma conditions, the catalyst material, and their mutual interaction. In this work we use microkinetic models to study how the total conversion process is impacted by the combination of different plasma conditions and transition metal catalysts. We study how plasma-generated radicals and vibrationally excited N2 (present in a dielectric barrier discharge plasma) interact with the catalyst and impact the NH3 turnover frequencies (TOFs). Both filamentary and uniform plasmas are studied, based on plasma chemistry models that provided plasma phase speciation and vibrational distribution functions. The Langmuir−Hinshelwood reaction rate coefficients (i.e., adsorption reactions and subsequent reactions among adsorbates) are determined using conventional scaling relations. An additional set of Eley−Rideal reactions (i.e., direct reactions of plasma radicals with adsorbates) was added and a sensitivity analysis on the assumed reaction rate coefficients was performed. We first show the impact of different vibrational distribution functions on the catalytic dissociation of N2 and subsequent production of NH3, and we gradually include more radical reactions, to illustrate the contribution of these species and their corresponding reaction pathways. Analysis over a large range of catalysts indicates that different transition metals (metals such as Rh, Ni, Pt, and Pd) optimize the NH3TOFs depending on the population of the vibrational levels of N2. At higher concentrations of plasma-generated radicals, the NH3 TOFs become less dependent on the catalyst material, due to radical adsorptions on the more noble catalysts and Eley−Rideal reactions on the less noble catalysts.

Abstract: We present an approach to incorporate the effect of vibrational nonequilibrium in molecular dynamics (MD) simulations. A perturbed canonical ensemble, in which selected modes are excited to higher temperature while all others remain equilibrated at low temperature, is simulated by applying a specifically tailored bias potential. Our method can be readily applied to any (classical or quantum mechanical) MD setup at virtually no additional computational cost and allows the study of reactions of vibrationally excited molecules in nonequilibrium environments such as plasmas. In combination with enhanced sampling methods, the vibrational efficacy and mode selectivity of vibrationally stimulated reactions can then be quantified in terms of chemically relevant observables, such as reaction rates and apparent free energy barriers. We first validate our method for the prototypical hydrogen exchange reaction and then show how it can capture the effect of vibrational excitation on a symmetric SN2 reaction and radical addition on CO2.

Abstract: A two-dimensional model is presented that describes the behavior of argon metastable atoms, copper atoms, and copper ions in an argon direct. current glow discharge, in the standard cell of the VG9000 glow discharge mass spectrometer for analyzing flat samples. The model is combined with a previously developed model for the electrons, argon ions, and atoms in the same cell to obtain an overall picture of the glow discharge, The results of the present model comprise the number densities of the described plasma species, the relative contributions of different production and loss processes for the argon metastable atoms, the thermalization profile of the sputtered copper atoms, the relative importance of the different ionization mechanisms for the copper atoms, the ionization degree of copper, the copper ion-to-argon ion density ratio, and the relative roles of copper ions, argon ions, and atoms in the sputtering process. All these quantities are calculated for a range of voltages and pressures, Moreover, since the sticking coefficient of copper atoms on solid surfaces is not well-known in the literature, the influence of this parameter on the results is briefly discussed.

Abstract: Despite modern therapeutic advances, the survival prospects of pancreatic cancer patients remain poor, due to chemoresistance and dysregulated oncogenic kinase signaling networks. We applied a novel kinome activitymapping approach using biological peptide targets as phospho-sensors to identify vulnerable kinase de­pendencies for therapy sensitization by physical plasma. Ser/Thr-kinome specific activity changes were mapped upon induction of ferroptotic cell death in pancreatic tumor cells exposed to reactive oxygen and nitrogen species of plasma-treated water (PTW). This revealed a broad kinome activity response involving the CAMK, the AGC and CMGC family of kinases. This systems-level kinome network response supports stress adaptive switches between chemoresistant anti-oxidant responses of Kelch-like ECH-associated protein 1 (KEAP1)/Heme Oxy­genase 1 (HMOX1) and ferroptotic cell death sensitization upon suppression of Nuclear factor (erythroid derived 2)-like 2 (NRF2) and Glutathione peroxidase 4 (GPX4). This is further supported by ex vivo experiments in the chicken chorioallantoic membrane assay, showing decreased GPX4 and Glutathione (GSH) expression as well as increased lipid peroxidation, along with suppressed BxPC-3 tumor growth in response to PTW. Taken all together, we demonstrate that plasma treated water-derived oxidants sensitize pancreatic cancer cells to fer­roptotic cell death by targeting a NRF2-HMOX1-GPX4 specific kinase signaling network.

Abstract: Binding of the SARS-CoV-2 S-glycoprotein to cell receptors is vital for the entry of the virus into cells and subsequent infection. ACE2 is the main cell receptor for SARS-CoV-2, which can attach to the C-terminal receptor-binding domain (RBD) of the SARS-CoV-2 S-glycoprotein. The GRP78 receptor plays an anchoring role, which attaches to the RBD and increases the chance of other RBDs binding to ACE2. Although high levels of reactive oxygen and nitrogen species (RONS) are produced during viral infections, it is not clear how they affect the RBD structure and its binding to ACE2 and GRP78. In this research, we apply molecular dynamics simulations to study the effect of oxidation of the highly reactive cysteine (Cys) amino acids of the RBD on its binding to ACE2 and GRP78. The interaction energy of both ACE2 and GRP78 with the whole RBD, as well as with the RBD main regions, is compared in both the native and oxidized RBDs. Our results show that the interaction energy between the oxidized RBD and ACE2 is strengthened by 155 kJ/mol, increasing the binding of the RBD to ACE2 after oxidation. In addition, the interaction energy between the RBD and GRP78 is slightly increased by 8 kJ/mol after oxidation, but this difference is not significant. Overall, these findings highlight the role of RONS in the binding of the SARS-CoV-2 S-glycoprotein to host cell receptors and suggest an alternative mechanism by which RONS could modulate the entrance of viral particles into the cells.

Abstract: Catalysis science and technology increased efforts recently to progress beyond conventional “thermal” catalysis and face the challenges of net-zero emissions and electrification of production. Nevertheless, a better gaps and opportunities analysis is necessary. This review analyses four emerging areas of unconventional or less- conventional catalysis which share the common aspect of using directly renewable energy sources: (i) plasma catalysis, (ii) catalysis for flow chemistry and process intensification, (iii) application of electromagnetic (EM) fields to modulate catalytic activity and (iv) nanoscale generation at the catalyst interface of a strong local EM by plasmonic effect. Plasma catalysis has demonstrated synergistic effects, where the outcome is higher than the sum of both processes alone. Still, the underlying mechanisms are complex, and synergy is not always obtained. There is a crucial need for a better understanding to (i) design catalysts tailored to the plasma environment, (ii) design plasma reactors with optimal transport of plasma species to the catalyst surface, and (iii) tune the plasma conditions so they work in optimal synergy with the catalyst. Microfluidic reactors (flow chemistry) is another emerging sector leading to the intensification of catalytic syntheses, particularly in organic chemistry. New unconventional catalysts must be designed to exploit in full the novel possibilities. With a focus on (a) continuous-flow photocatalysis, (b) electrochemical flow catalysis, (c) microwave flow catalysis and (d) ultra­ sound flow activation, a series of examples are discussed, with also indications on scale-up and process indus­ trialisation. The third area discussed regards the effect on catalytic performances of applying oriented EM fields spanning several orders of magnitude. Under well-defined conditions, gas breakdown and, in some cases, plasma formation generates activated gas phase species. The EM field-driven chemical conversion processes depend further on structured electric/magnetic catalysts, which shape the EM field in strength and direction. Different effects influencing chemical conversion have been reported, including reduced activation energy, surface charging, hot spot generation, and selective local heating. The last topic discussed is complementary to the third, focusing on the possibility of tuning the photo- and electro-catalytic properties by creating a strong localised electrical field with a plasmonic effect. The novel possibilities of hot carriers generated by the plasmonic effect are also discussed. This review thus aims to stimulate the reader to make new, creative catalysis to address the challenges of reaching a carbon-neutral world.

Abstract: We developed a novel in situ (i.e. inside plasma and during operation) IR dielectric barrier discharge cell allowing investigation of plasma catalysis in transmission mode, atmospheric pressure, flow conditions (WHSV similar to 0-50 000 mL g(-1) h(-1)), at relevant discharge voltages (similar to 0-50 kV) and frequencies (similar to 0-5 kHz). We applied it to study the IR-active surface species formed on a SiO2 support and on a 3 wt% Ru/SiO2 catalyst, which can help to reveal the important surface reaction mechanisms during the plasma-catalytic dry reforming of methane (DRM). Moreover, we present a technique for the challenging task of estimating the temperature of a catalyst sample in a plasma-catalytic system in situ and during plasma operation. We found that during the reaction, water is immediately formed at the SiO2 surface, and physisorbed formic acid is formed with a delay. As Ru/SiO2 is subject to greater plasma-induced heating than SiO2 (with a surface temperature increase in the range of 70-120 degrees C, with peaks up to 150 degrees C), we observe lower amounts of physisorbed water on Ru/SiO2, and less physisorbed formic acid formation. Importantly, the formation of surface species on the catalyst sample in our plasma-catalytic setup, as well as the observed conversions and selectivities in plasma conditions, can not be explained by plasma-induced heating of the catalyst surface, but must be attributed to other plasma effects, such as the adsorption of plasma-generated radicals and molecules, or the occurrence of Eley-Rideal reactions.

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: Mesoporous carbon, which can be replicated from mesoporous silica and whose surface is hydrophobic, can be an ideal catalyst for the esterification of fatty acids. Here we report an easy and low cost way to prepare sulfonic acid group-functionalized mesoporous carbon. A sample of calcined mesoporous silica SBA-15 was added to an aqueous sucrose solution followed by drying and calcination at different temperatures. In contrast to existing procedures, the obtained hybrid Si/C material was then first sulfonated in H2SO4, before the final removal of the silica template in order to stabilize the porous structure towards the liquid phase sulfonation treatment. Thus the silicacarbon composites, instead of the mesoporous carbon, were successfully sulfonated to introduce SO3H groups, while keeping the ordered mesoporous structure intact. The influence of carbonization temperature was investigated, suggesting an optimum temperature of 873 K. The SO3H group-functionalized mesoporous carbon, denoted as CMK-3-873-SO3H, was characterized by means of XRD, N2 physisorption, SEM, FT-IR, elemental analysis and TEM. It followed that a uniform mesoporous carbon was obtained with an average pore size of 3.89 nm, a specific surface of 807 m2/g and a SO3H group loading of 0.39 meq/g of dry material. Compared with other solid acid catalysts, the resulting material shows enhanced activity in the acid-catalyzed esterification of oleic acid with methanol, and can be used repeatedly. The increased catalytic performance is attributed to the hydrophobic surface and larger pore size of the new catalyst. It can effectively accommodate long chain fatty acids and reject formed water, making the active sites easily accessible.

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: Physical plasmas generate unique mixes of reactive oxygen and nitrogen species (RONS or ROS). Only a bit more than a decade ago, these plasmas, operating at body temperature, started to be considered for medical therapy with considerably little mechanistic redox chemistry or biomedical research existing on that topic at that time. Today, a vast body of evidence is available on physical plasma-derived ROS, from their spatiotemporal resolution in the plasma gas phase to sophisticated chemical and biochemical analysis of these species once dissolved in liquids. Data from<italic>in silico</italic>analysis dissected potential reaction pathways of plasma-derived reactive species with biological membranes, and<italic>in vitro</italic>and<italic>in vivo</italic>experiments in cell and animal disease models identified molecular mechanisms and potential therapeutic benefits of physical plasmas. In 2013, the first medical plasma systems entered the European market as class IIa devices and have proven to be a valuable resource in dermatology, especially for supporting the healing of chronic wounds. The first results in cancer patients treated with plasma are promising, too. Due to the many potentials of this blooming new field ahead, there is a need to highlight the main concepts distilled from plasma research in chemistry and biology that serve as a mechanistic link between plasma physics (how and which plasma-derived ROS are produced) and therapy (what is the medical benefit). This inevitably puts cellular membranes in focus, as these are the natural interphase between ROS produced by plasmas and translation of their chemical reactivity into distinct biological responses.

Abstract: The phase transition processes induced by ultrashort, 100 fs pulsed laser irradiation of Au, Cu, and Ni are studied by means of a combined atomistic-continuum approach. A moderately low absorbed laser fluence range, from 200 to 600 J/m2 is considered to study phase transitions by means of a local and a nonlocal order parameter. At low laser fluences, the occurrence of layer-by-layer evaporation has been observed, which suggests a direct solid to vapor transition. The calculated amount of molten material remains very limited under the conditions studied, especially for Ni. Therefore, our results show that a kinetic equation that describes a direct solid to vapor transition might be the best approach to model laser-induced phase transitions by continuum models. Furthermore, the results provide more insight into the applicability of analytical superheating theories that were implemented in continuum models and help the understanding of nonequilibrium phase transitions.

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: 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.