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: In this paper we study dry reforming of methane (DRM) in a gliding arc plasmatron (GAP) in the presence of N₂and O₂. N₂is added to create a stable plasma at equal fractions of CO₂and CH₄, and because emissions from industrial plants typically contain N₂, while O₂is added to enhance the process. We test different gas mixing ratios to evaluate the conversion and energy cost. We obtain conversions between 31 and 52% for CO₂and between 55 and 99% for CH₄, with total energy costs between 3.4 and 5.0 eV per molecule, depending on the gas mixture. This is very competitive when benchmarked with the literature. In addition, we present a chemical kinetics model to obtain deeper insight in the underlying plasma chemistry. This allows determination of the major reaction pathways to convert CO₂and CH₄, in the presence of O₂and N₂, into CO and H₂. We show that N₂assists in the CO₂conversion, but part of the applied energy is also wasted in N₂excitation. Adding O₂enhances the CH₄conversion, and lowers the energy cost, while the CO₂conversion remains constant, and only slightly drops at the highest O₂fractions studied, when CH₄is fully oxidized into CO₂.

Abstract: This paper gives an overview of the modeling work developed in our group in the last 25 years for various plasmas used in analytical spectrochemistry, i.e., glow discharges (GDs), inductively coupled plasmas (ICPs), and laser ablation (LA) for sample introduction in the ICP and for laser-induced breakdown spectroscopy (LIBS). The modeling approaches are briefly presented, which are different for each case, and some characteristic results are illustrated. These plasmas are used not only in analytical chemistry but also in other applications, and the insights obtained in these other fields were quite helpful for us to develop models for the analytical plasmas. Likewise, there is now a huge interest in plasma–liquid interaction, atmospheric pressure glow discharges (APGDs), and dielectric barrier discharges (DBDs) for environmental, medical, and materials applications of plasmas. The insights obtained in these fields are also very relevant for ambient desorption/ionization sources and for liquid sampling, which are nowadays very popular in analytical chemistry, and they could be very helpful in developing models for these sources as well.

Abstract: The combination of catalysis and nonthermal plasma holds promise for enabling difficult chemical conversions. The possible synergy between both depends strongly on the nature of the reactive plasma species and the catalyst material. In this paper, we show how vibrationally excited species and plasma-generated radicals interact with transition metal catalysts and how changing the catalyst material can improve the conversion rates and product selectivity. We developed a microkinetic model to investigate the impact of vibrational excitations and plasma-generated radicals on the nonoxidative coupling of methane over transition metal surfaces. We predict a significant increase in ethylene formation for vibrationally excited methane. Plasma-generated radicals have a stronger impact on the turnover frequencies with high selectivity toward ethylene on noble catalysts and mixed selectivity on non-noble catalysts. In general, we show how the optimal catalyst material depends on the desired products as well as the plasma conditions.

Abstract: DBD plasma reactors are commonly used in a static ‘one inlet – one outlet’ design that goes against reactor design principles for multi-component reactions, such as dry reforming of methane (DRM). Therefore, in this paper we have developed a novel reactor design, and investigated how the shape and size of the reaction zone, as well as gradual gas addition, and the method of mixing CO2 and CH4 can influence the conversion and product com­ position of DRM. Even in the standard ‘one inlet – one outlet’ design, the direction of the gas flow (i.e. short or long path through the reactor, which defines the gas velocity at fixed residence time), as well as the dimensions of the reaction zone and the power delivery to the reactor, largely affect the performance. Using gradual gas addition and separate plasma activation zones for the individual gases give increased conversions within the same operational parameters, by optimising mixing ratios and kinetics. The choice of the main (pre-activated) gas and the direction of gas flow largely affect the conversion and energy cost, while the gas inlet position during separate addition only influences the product distribution.

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: Plasma reactors are interesting for gas-based chemical conversion but the fundamental relation between the plasma chemistry and selected conditions remains poorly understood. Apparent kinetic parameters for the loss and formation processes of individual components of gas conversion processes, can however be extracted by performing experiments in an extended residence time range (2–75 s) and fitting the gas composition to a firstorder kinetic model of the evolution towards partial chemical equilibrium (PCE). We specifically investigated the differences in kinetic characteristics and PCE state of the CO2 dissociation and CH4 reforming reactions in a dielectric barrier discharge reactor (DBD), how these are mutually affected when combining both gases in the dry reforming of methane (DRM) reaction, and how they change when a packing material (non-porous SiO2) is added to the reactor. We find that CO2 dissociation is characterized by a comparatively high reaction rate of 0.120 s−1 compared to CH4 reforming at 0.041 s−1; whereas CH4 reforming reaches higher equilibrium conversions, 82% compared to 53.6% for CO2 dissociation. Combining both feed gases makes the DRM reaction to proceed at a relatively high rate (0.088 s−1), and high conversion (75.4%) compared to CO2 dissociation, through accessing new chemical pathways between the products of CO2 and CH4. The addition of the packing material can also distinctly influence the conversion rate and position of the equilibrium, but its precise effect depends strongly on the gas composition. Comparing different CO2:CH4 ratios reveals the delicate balance of the combined chemistry. CO2 drives the loss reactions in DRM, whereas CH4 in the mixture suppresses back reactions. As a result, our methodology provides some of the insight necessary to systematically tune the conversion process.

Abstract: Ammonia is one of the most produced chemicals, mainly synthesized from fossil fuels for fertilizer applications. Furthermore, ammonia may be one of the energy carriers of the future, when it is produced from renewable electricity. This has spurred research on alternative technologies for green ammonia production. Research on plasma-driven ammonia synthesis has recently gained traction in academic literature. In the current review, we summarize the literature on plasma-driven ammonia synthesis. We distinguish between mechanisms for ammonia synthesis in the presence of a plasma, with and without a catalyst, for different plasma conditions. Strategies for catalyst design are discussed, as well as the current understanding regarding the potential plasma-catalyst synergies as function of the plasma conditions and their implications on energy efficiency. Finally, we discuss the limitations in currently reported models and experiments, as an outlook for research opportunities for further unravelling the complexities of plasma-catalytic ammonia synthesis, in order to bridge the gap between the currently reported models and experimental results.

Abstract: Biological membranes are under constant attack of free radicals, which may lead to lipid nitro-oxidation, pro­ ducing a complex mixture of nitro-oxidized lipids that are responsible for structural and dynamic changes on the membrane. Despite the latter, nitro-oxidized lipids are also associated with several inflammatory and neuro­ degenerative diseases, the underlying mechanisms of which remain elusive. We perform atomistic molecular dynamics simulations using several isomers of nitro-oxidized lipids to study their effect on the structure and permeability of the membrane, as well as the interaction between the mixture of these products in the phos­pholipid membrane environment. Our results show that the stereo- and positional isomers have a stronger effect on the properties of the membrane composed of oxidized lipids compared to that containing nitrated lipids. Nevertheless, nitrated lipids lead to three-fold increase in water permeability compared to oxidized lipids. In addition, we show that in a membrane consisting of combined nitro-oxidized lipid products, the presence of oxidized lipids protects the membrane from transient pores. Is well stablished that plasma application and photodynamic therapy produces a number of oxidative species used to kill cancer cells, through membrane damage induced by nitro-oxidative stress. This study is important to elucidate the mechanisms and the molecular level properties involving the reactive species produced during that cancer therapies.

Abstract: We use microkinetic modeling to examine the potential of plasma-catalytic partial oxidation (POX) of CH4 as a promising new approach to produce oxygenates. We study how different plasma species affect POX of CH4 on the Pt(111) surface, and we discuss the associated kinetic and mechanistic changes. We discuss the effect of vibrationally excited CH4 and O2, as well as plasma-generated radicals and stable intermediates. Our results show that vibrational excitation enhances the turnover frequency (TOF) of catalytic CH4 dissociation and has good potential for improving the selectivities toward CH3OH, HCOOH, and C2 hydrocarbons. Nevertheless, when also considering plasma-generated radicals, we find that these species mainly govern the surface chemistry. Additionally, we find that plasma-generated radicals and stable intermediates enhance the TOFs of COx and oxygenates, increase the selectivity toward oxygenates, and make the formation of HCOOH more significant on Pt(111). We also briefly illustrate the potential impact of Eley−Rideal reactions that involve plasma-generated radicals. Finally, we reveal how various radicals affect the catalyst surface chemistry and we link this to the formation of different products. This allows us to make suggestions on how the plasma composition should be altered to improve the formation of desired products.

Abstract: While spectroscopic data on small hydrocarbons in interstellar media in combination with crossed molecular beam (CMB) experiments have provided a wealth of information on astrochemically relevant species, much of the underlying mechanistic pathways of their formation remain elusive. Therefore, in this work, the chemical reaction mechanisms of C(³P_J) + C₆H₆and C⁺(²P) + C₆H₆systems using the quantum mechanical molecular dynamics (QMMD) technique at the PBE0-D3(BJ) level of theory is investigated, mimicking a CMB experiment. Both the dynamics of the reactions as well as the electronic structure for the purpose of the reaction network are evaluated. The method is validated for the first reaction by comparison to the available experimental data. The reaction scheme for the C(³P_J) + C₆H₆system covers the literature data,<italic>e.g.</italic>the major products are the 1,2-didehydrocycloheptatrienyl radical (C₇H₅) and benzocyclopropenyl radical (C₆H₅–CH), and it reveals the existence of less common pathways for the first time. The chemistry of the C⁺(²P_J) + C₆H₆system is found to be much richer, and we have found that this is because of more exothermic reactions in this system in comparison to those in the C(³P_J) + C₆H₆system. Moreover, using the QMMD simulation, a number of reaction paths have been revealed that produce three distinct classes of reaction products with different ring sizes. All in all, at all the collision energies and orientations, the major product is the heptagon molecular ion for the ionic system. It is also revealed that the collision orientation has a dominant effect on the reaction products in both systems, while the collision energy mostly affects the charged system. These simulations both prove the applicability of this approach to simulate crossed molecular beams, and provide fundamental information on reactions relevant for the interstellar medium.

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: N-Based fertilisers are paramount to support our still-growing world population. Current industrial N₂fixation is heavily fossil fuel-dependent, therefore, a lot of work is put into the development of fossil-free pathways. Plasma technology offers a fossil-free and flexible method for N₂fixation that is compatible with renewable energy sources. We present here a pulsed plasma jet for direct NO_{<italic>x</italic>}production from air. The pulsed power allows for a record-low energy consumption (EC) of 0.42 MJ (mol N)⁻¹. This is the lowest reported EC in plasma-based N₂fixation at atmospheric pressure thus far. We compare our experimental data with plasma chemistry modelling, and obtain very good agreement. Hence, we can use our model to explain the underlying mechanisms responsible for this low EC. The pulsed power and the corresponding pulsed gas temperature are the reason for the very low EC: they provide a strong vibrational–translational non-equilibrium and promote the non-thermal Zeldovich mechanism. This insight is important for the development of the next generation of plasma sources for energy-efficient NO_{<italic>x</italic>}production.

Abstract: It is well established that lipid aldehydes (LAs) are able to increase the permeability of cell membranes and induce their rupture. However, it is not yet clear how LAs are distributed in phase-separated membranes (PSMs), which are responsible for the transport of selected molecules and intracellular signaling. Thus, we investigate here the distribution of LAs in a PSM by coarse-grained molecular dynamics simulations. Our results reveal that LAs derived from mono-unsaturated lipids tend to accumulate at the interface between the liquid-ordered/liquiddisordered domains, whereas those derived from poly-unsaturated lipids remain in the liquid-disordered domain. These results are important for understanding the effects caused by oxidized lipids in membrane structure, properties and organization.

Abstract: Plasma-based NOx production is of interest for sustainable N2 fixation, but more research is needed to improve its performance. One of the current limitations is recombination of NO back into N2 and O2 molecules immediately after the plasma reactor. Therefore, we developed a novel so-called “effusion nozzle”, to improve the perfor­mance of a rotating gliding arc plasma reactor for NOx production, but the same principle can also be applied to other plasma types. Experiments in a wide range of applied power, gas flow rates and N2/O2 ratios demonstrate an enhancement in NOx concentration by about 8%, and a reduction in energy cost by 22.5%. In absolute terms, we obtain NOx concentrations up to 5.9%, at an energy cost down to 2.1 MJ/mol, which are the best values reported to date in literature. In addition, we developed four complementary models to describe the gas flow, plasma temperature and plasma chemistry, aiming to reveal why the effusion nozzle yields better performance. Our simulations reveal that the effusion nozzle acts as very efficient heat sink, causing a fast drop in gas tem­perature when the gas molecules leave the plasma, hence limiting the recombination of NO back into N2 and O2. This yields an overall higher NOx concentration than without the effusion nozzle. This immediate quenching right at the end of the plasma makes our effusion nozzle superior to more conventional cooling options, like water cooling In addition, this higher NOx concentration can be obtained at a slightly lower power, because the effusion nozzle allows for the ignition and sustainment of the plasma at somewhat lower power. Hence, this also explains the lower energy cost. Overall, our experimental results and detailed modeling analysis will be useful to improve plasma-based NOx production in other plasma reactors as well.