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: Understanding the folding and stability of membrane proteins is of great importance in protein science. Recently, osmolytes and ionic liquids (ILs) are increasingly being used as drug delivery systems in the biopharmaceutical industry. However, the stability of membrane proteins in the presence of osmolytes and ILs is not yet fully understood. Besides, the effect of oxidative stress on membrane proteins with osmolytes or ILs has not been investigated. Therefore, we studied the influence of osmolytes and ILs as co-solvents on the stability of a model membrane protein (i.e., Bacteriorhodopsin in purple membrane of Halobacterium salinarum), using UV–Vis spectroscopy and molecular dynamics (MD) simulations. The MD simulations allowed us to determine the flexibility and solvent accessible surface area (SASA) of Bacteriorhodopsin protein in the presence and/or absence of cosolvents, as well as to carry out principal component analysis (PCA) to identify the most important movements in this protein. In addition, by means of UV–Vis spectroscopy we studied the effect of oxidative stress generated by cold atmospheric plasma on the stability of Bacteriorhodopsin in the presence and/or absence of co-solvents. This study is important for a better understanding of the stability of proteins in the presence of oxidative stress.

Abstract: Ammonia is a crucial nutrient used for plant growth and as a building block in pharmaceutical and chemical industry, produced via nitrogen fixation of the ubiquitous atmospheric N2. Current industrial ammonia production relies heavily on fossil resources, but a lot of work is put into developing non-fossil based pathways. Among these is the use of nonequilibrium plasma. In this work, we investigated water vapor as H source for nitrogen fixation into NH3 by non-equilibrium plasma. The highest selectivity towards NH3 was observed with low amounts of added H2O vapor, but the highest production rate was reached at high H2O vapor.

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: Atmospheric pressure gliding arc (GA) discharges are gaining increasing interest for CO2 conversion and other gas conversion applications, due to their simplicity and high energy efficiency. However, they are characterized by some drawbacks, such as non-uniform gas treatment, limiting the conversion, as well as the development of a hot cathode spot, resulting in severe electrode degradation. In this work, we built a dual-vortex plasmatron, which is a GA plasma reactor with innovative electrode configuration, to solve the above problems. The design aims to improve the CO2 conversion capability of the GA reactor by elongating the arc in two directions, to increase the residence time of the gas inside the arc, and to actively cool the cathode spot by rotation of the arc and gas convection. The measured CO2 conversion and corresponding energy efficiency indeed look very promising. In addition, we developed a fluid dynamics non-thermal plasma model with argon chemistry, to study the arc behavior in the reactor and to explain the experimental results.

Abstract: We present a zero-dimensional plasma kinetics model, including both surface and gas phase kinetics, to determine the role of vibrationally excited states in plasma-catalytic ammonia synthesis. We defined a new method to systematically capture the conditions of dielectric barrier discharges (DBDs), including those found in packed bed DBDs. We included the spatial and temporal nature of such discharges by special consideration of the number of micro-discharges in the model. We introduce a parameter that assigns only a part of the plasma power to the microdischarges, to scale the model conditions from filamentary to uniform plasma. Because of the spatial and temporal behaviour of the micro-discharges, not all micro-discharges occurring in the plasma reactor during a certain gas residence time are affecting the molecules. The fraction of power considered in the model ranges from 0.005 %, for filamentary plasma, to 100 %, for uniform plasma. If vibrational excitation is included in the plasma chemistry, these different conditions, however, yield an ammonia density that is only varying within one order of magnitude. At only 0.05 % of the power put into the uniform plasma component, a model neglecting vibrational excitation clearly does not result in adequate amounts of ammonia. Thus, our new model, which accounts for the concept in which not all the power is deposited by the micro-discharges, but some part may also be distributed in between them, suggests that vibrational kinetic processes are really important in (packed bed) DBDs. Indeed, vibrational excitation takes place in both the uniform plasma between the micro-discharges and in the strong micro-discharges, and is responsible for an increased N2 dissociation rate. This is shown here for plasma-catalytic ammonia synthesis, but might also be valid for other gas conversion applications.

Abstract: In this work, we use density functional theory calculations to study the combined effect of external electric fields, surface morphology, and surface charge on CO2 activation over Cu(111), Cu(211), Cu(110), and Cu(001) surfaces. We observe that the binding energy of the CO2 molecule on Cu surfaces increases significantly upon increasing the applied electric field strength. In addition, rougher surfaces respond more effectively to the presence of the external electric field toward facilitating the formation of a carbonate-like CO2 structure and the transformation of the most stable adsorption mode from physisorption to chemisorption. The presence of surface charges further strengthens the electric field effect and consequently causes an improved bending of the CO2 molecule and C−O bond length elongation. On the other hand, a net charge in the absence of an externally applied electric field shows only a marginal effect on CO2 binding. The chemisorbed CO2 is more stable and further activated when the effects of an external electric field, rough surface, and surface charge are combined. These results can help to elucidate the underlying factors that control CO2 activation in heterogeneous and plasma catalysis, as well as in electrochemical processes.

Abstract: Mcl1 is an important anti-apoptotic member of the Bcl2 family proteins that are upregulated in several cancer malignancies. The canonical binding groove (CBG) located at the surface of Mcl1 exhibits a critical role in binding partners selectively via the BH3-domain of pro-apoptotic Bcl2 family members that trigger the downregulation of Mcl1 function. There are several crystal structures of point-mutated pro-apoptotic Bim peptides in complex with Mcl1. However, the mechanistic effects of such point-mutations towards peptide binding and complex stability still remain unexplored. Here, the effects of the reported point mutations in Bim peptides and their binding mechanisms to Mcl1 were computationally evaluated using atomistic-level steered molecular dynamics (SMD) simulations. A range of external-forces and constant-velocities were applied to the Bim peptides to uncover the mechanistic basis of peptide dissociation from the CBG of Mcl1. Although the peptides showed similarities in their dissociation pathways, the peak rupture forces varied significantly. According to simulations results, the disruption of the conserved polar contacts at the complex interface causes a sequential release of the peptides from the CBG of Mcl1. Overall, the results obtained from the current study may provide valuable insights for the development of novel anti-cancer peptide-inhibitors that can downregulate Mcl1’s function.

Abstract: The penetration of reactive oxygen and nitrogen species (RONS) across the stratum corneum (SC) is a necessary and crucial process in many skin‐related plasma medical applications. To gain more insights into this penetration behavior, we combined experimental measurements of the permeability of dry and moist SC layers with computer simulations of model lipid membranes. We measured the permeation of relatively stable molecules, which are typically generated by plasma, namely H2O2, NO3−, and NO2−. Furthermore, we calculated the permeation free energy profiles of the major plasma‐generated RONS and their derivatives (i.e., H2O2, OH, HO2, O2, O3, NO, NO2, N2O4, HNO2, HNO3, NO2−, and NO3−) across native and oxidized SC lipid bilayers, to understand the mechanisms of RONS permeation across the SC. Our results indicate that hydrophobic RONS (i.e., NO, NO2, O2, O3, and N2O4) can translocate more easily across the SC lipid bilayer than hydrophilic RONS (i.e., H2O2, OH, HO2, HNO2, and HNO3) and ions (i.e., NO2− and NO3−) that experience much higher permeation barriers. The permeability of RONS through the SC skin lipids is enhanced when the skin is moist and the lipids are oxidized. These findings may help to understand the underlying mechanisms of plasma interaction with a biomaterial and to optimize the environmental parameters in practice in plasma medical applications.

Abstract: Plasma technology provides a sustainable, fossil-free method for N2 fixation, i.e., the conversion of inert atmospheric N2 into valuable substances, such as NOx or ammonia. In this work, we present a novel gliding arc plasmatron at atmospheric pressure for NOx production at different N2/O2 gas feed ratios, offering a promising NOx yield of 1.5% with an energy cost of 3.6 MJ/mol NOx produced. To explain the underlying mechanisms, we present a chemical kinetics model, validated by experiments, which provides insight into the NOx formation pathways and into the ambivalent role of the vibrational kinetics. This allows us to pinpoint the factors limiting the yield and energy cost, which can help to further improve the process.

Abstract: We explore the consequences of nonthermal plasma-activation on product yields in catalytic ammonia synthesis, a reaction that is equilibrium-limited at elevated temperatures. We employ a minimal microkinetic model that incorporates the influence of plasma-activation on N2 dissociation rates to predict NH3 yields into and across the equilibrium-limited regime. NH3 yields are predicted to exceed bulk thermodynamic equilibrium limits on materials that are thermal-rate-limited by N2 dissociation. In all cases, yields revert to bulk equilibrium at temperatures at which thermal reaction rates exceed plasma-activated ones. Beyond-equilibrium NH3 yields are observed in a packed bed dielectric barrier discharge reactor and exhibit sensitivity to catalytic material choice in a way consistent with model predictions. The approach and results highlight the opportunity to exploit synergies between nonthermal plasmas and catalysts to affect transformations at conditions inaccessible through thermal routes.

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: We studied hydrogen sulfide (H2S) decomposition into hydrogen (H2) and sulfur (S2) in a gliding arc plasmatron (GAP) and microwave (MW) plasma by a combination of 0D and 2D models. The conversion, energy efficiency, and plasma distribution are examined for different discharge conditions, and validated with available experiments from literature. Furthermore, a comparison is made between GAP and MW plasma. The GAP operates at atmospheric pressure, while the MW plasma experiments to which comparison is made were performed at reduced pressure. Indeed, the MW discharge region becomes very much contracted near atmospheric pressure, at the conditions under study, as revealed by our 2D model. The models predict that thermal reactions play the most important role in H2S decomposition in both plasma types. The GAP has a higher energy efficiency but lower conversion than the MW plasma at their typical conditions. When compared at the same conversion, the GAP exhibits a higher energy efficiency and lower energy cost than the MW plasma.

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: Cancer is the second leading cause of death worldwide, and while science has advanced significantly to improve the treatment outcome and quality of life in cancer patients, there are still many issues with the current therapies, such as toxicity and the development of resistance to treatment [...]