Abstract: We investigate the performance of an atmospheric pressure glow discharge (APGD) reactor for CO2 conversion in three different configurations, through experiments and simulations. The first (basic) configuration utilizes the well-known pin-to-plate design, which offers a limited conversion. The second configuration improves the reactor performance by employing a vortex-flow generator. The third, “confined” configuration is a complete redesign of the reactor, which encloses the discharge in a limited volume, significantly surpassing the conversion rate of the other two designs. The plasma properties are investigated using an advanced plasma model.

Abstract: We performed computer simulations to investigate the effect of oxidation on the extracellular cystine (CYC) uptake by the xC− antiporter. The latter is important for killing of cancer cells. Specifically, applying molecular dynamics (MD) simulations we studied the transport of CYC across xCT, i.e., the light subunit of the xC− antiporter, in charge of bidirectional transport of CYC and glutamate. We considered the outward facing (OF) configuration of xCT, and to study the effect of oxidation, we modified the Cys327 residue, located in the vicinity of the extracellular milieu, to cysteic acid (CYO327). Our computational results showed that oxidation of Cys327 results in a free energy barrier for CYC translocation, thereby blocking the access of CYC to the substrate binding site of the OF system. The formation of the energy barrier was found to be due to the conformational changes in the channel. Analysis of the MD trajectories revealed that the reorganization of the side chains of the Tyr244 and CYO327 residues play a critical role in the OF channel blocking. Indeed, the calculated distance between Tyr244 and either Cys327 or CYO327 showed a narrowing of the channel after oxidation. The obtained free energy barrier for CYC translocation was found to be 33.9kJmol−1, indicating that oxidation of Cys327, by e.g., cold atmospheric plasma, is more effective in inhibiting the xC− antiporter than in the mutation of this amino acid to Ala (yielding a barrier of 32.4kJmol−1). The inhibition of the xC− antiporter may lead to Cys starvation in some cancer cells, eventually resulting in cancer cell death.

Abstract: The mechanisms of plasma in medicine are broadly attributed to plasma-derived reactive oxygen and nitrogen species (RONS). In order to exert any intracellular effects, these plasma-derived RONS must first traverse a major barrier in the cell membrane. The cell membrane lipid composition, and thereby the magnitude of this barrier, is highly variable between cells depending on type and state (e.g. it is widely accepted that healthy and cancerous cells have different membrane lipid compositions). In this study, we investigate how plasma-derived RONS interactions with lipid membrane components can potentially be exploited in the future for treatment of diseases. We couple phospholipid vesicle experiments, used as simple cell models, with molecular dynamics (MD) simulations of the lipid membrane to provide new insights into how the interplay between phospholipids and cholesterol may influence the response of healthy and diseased cell membranes to plasma-derived RONS. We focus on the (i) lipid tail saturation degree, (ii) lipid head group type, and (iii) membrane cholesterol fraction. Using encapsulated molecular probes, we study the influence of the above membrane components on the ingress of RONS into the vesicles, and subsequent DNA damage. Our results indicate that all of the above membrane components can enhance or suppress RONS uptake, depending on their relative concentration within the membrane. Further, we show that higher RONS uptake into the vesicles does not always correlate with increased DNA damage, which is attributed to ROS reactivity and lifetime. The MD simulations indicate the multifactorial chemical and physical processes at play, including (i) lipid oxidation, (ii) lipid packing, and (iii) lipid rafts formation. The methods and findings presented here provide a platform of knowledge that could be leveraged in the development of therapies relying on the action of plasma, in which the cell membrane and oxidative stress response in cells is targeted.

Abstract: Plasma-treated liquids have great potential for biomedical applications. However, insight into the underlying mechanisms and the exact chemistry is still scarce. In this study, we present the combination of a 0D chemical kinetics and a 2D fluid dynamics model to investigate the plasma treatment of a buffered water solution with the kINPen (R) plasma jet. Using this model, we calculated the gas and liquid flow profiles and the transport and chemistry of all species in the gas and the liquid phase. Moreover, we evaluated the stability of the reactive oxygen and nitrogen species after plasma treatment. We found that of all species, only H2O2, HNO2/NO2-, and HNO3/NO3- are stable in the buffered solution after plasma treatment. This is because both their production and loss processes in the liquid phase are dependent on short-lived radicals (e.g. OH, NO, and NO2). Apart from some discrepancy in the absolute values of the concentrations, which can be explained by the model, all general trends and observations in our model are in qualitative agreement with experimental data and literature.

Abstract: The combination of catalysis with non-thermal plasma is a promising alternative to thermal catalysis. A dielectric-barrier discharge reactor was used to study plasma-catalytic dry reforming of methane at ambient pressure and temperature and a fixed plasma power of 45 W. The effect of different catalytic packing materials was evaluated in terms of conversion, product selectivity, and energy efficiency. The conversion of CO2 (~22%) and CH4 (~33%) were found to be similar in plasma-only and when introducing packing materials in plasma. The main reason is the shorter residence time of the gas due to packing geometry, when compared at identical flow rates. H2, CO, C2-C4 hydrocarbons, and oxygenates were identified in the product gas. High selectivity towards H2 and CO were found for all catalysts and plasma-only, with a H2/CO molar ratio of ~0.9. The lowest syngas selectivity was obtained with Cu/Al2O3 (~66%), which instead, had the highest alcohol selectivity (~3.6%).

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: 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: 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: 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: CO2 conversion by plasma technology is gaining increasing interest. We present a carbon (charcoal) bed placed after a Gliding Arc Plasmatron (GAP) reactor, to enhance the CO2 conversion, promote O/O2 removal and in­ crease the CO fraction in the exhaust mixture. By means of an innovative (silo) system, the carbon is constantly supplied, to avoid carbon depletion upon reaction with O/O2. Using this carbon bed, the CO2 conversion is enhanced by almost a factor of two (from 7.6 to 12.6%), while the CO concentration even increases by a factor of three (from 7.2 to 21.9%), and O2 is completely removed from the exhaust mixture. Moreover, the energy ef­ ficiency of the conversion process drastically increases from 27.9 to 45.4%, and the energy cost significantly drops from 41.9 to 25.4 kJ.L− 1. We also present the temperature as a function of distance from the reactor outlet, as well as the CO2, CO and O2 concentrations and the temperature in the carbon bed as a function of time, which is important for understanding the underlying mechanisms. Indeed, these time-resolved measurements reveal that the initial enhancements in CO2 conversion and in CO concentration are not maintained in our current setup. Therefore, we present a model to study the gasification of carbon with different feed gases (i.e., O2, CO and CO2 separately), from which we can conclude that the oxygen coverage at the surface plays a key role in determining the product composition and the rate of carbon consumption. Indeed, our model insights indicate that the drop in CO2 conversion and in CO concentration after a few minutes is attributed to deactivation of the carbon bed, due to rapid formation of oxygen complexes at the surface.

Abstract: This study aims at sensing in situ reactive oxygen and nitrogen species (RONS) and specifically superoxide anion (O-2(center dot-)) in aqueous buffer solutions exposed to cold atmospheric plasmas (CAPs). CAPs were generated by ionizing He gas shielded with variable N-2/O-2 mixtures. Thanks to ultramicroelectrodes protected against the high electric fields transported by the ionization waves of CAPs, the production of superoxide and several RONS was electrochemically directly detected in liquids during their plasma exposure. Complementarily, optical emissive spectroscopy (OES) was used to study the plasma phase composition and its correlation with the chemistry in the exposed liquid. The specific production of O-2(center dot-), a biologically reactive redox species, was analyzed by cyclic voltammetry (CV), in both alkaline (pH 11), where the species is fairly stable, and physiological (pH 7.4) conditions, where it is unstable. To understand its generation with respect to the plasma chemistry, we varied the shielding gas composition of CAPs to directly impact on the RONS composition at the plasma-liquid interface. We observed that the production and accumulation of RONS in liquids, including O(2)(center dot-)depends on the plasma composition, with N-2-based shieldings providing the highest superoxide concentrations (few 10s of micromolar at most) and of its derivatives (hundreds of micromolar). In situ spectroscopic and electrochemical analyses provide a high resolution kinetic and quantitative understanding of the interactions between CAPs and physiological solutions for biomedical applications.

Abstract: Utilizing ammonia as a storage medium for hydrogen is currently receiving increased attention. A possible method to retrieve the hydrogen is by plasma-catalytic decomposition. In this work, we combined an experimental study, using a dielectric barrier discharge plasma reactor, with a plasma kinetic model, to get insights into the decomposition mechanism. The experimental results revealed a similar effect on the ammonia conversion when changing the flow rate and power, where increasing the specific energy input (higher power or lower flow rate) gave an increased conversion. A conversion as high as 82 % was achieved at a specific energy input of 18 kJ/Nl. Furthermore, when changing the discharge volume from 31 to 10 cm3, a change in the plasma distribution factor from 0.2 to 0.1 was needed in the model to best describe the conversions of the experiments. This means that a smaller plasma volume caused a higher transfer of energy through micro-discharges (non-uniform plasma), which was found to promote the decomposition of ammonia. These results indicate that it is the collisions between NH3 and the high-energy electrons that initiate the decomposition. Moreover, the rate of ammonia destruction was found by the model to be in the order of 1022 molecules/(cm3 s) during the micro-discharges, which is 5 to 6 orders of magnitude higher than in the afterglows. A considerable re-formation of ammonia was found to take place in the afterglows, limiting the overall conversion. In addition, the model revealed that implementation of packing material in the plasma introduced high concentrations of surface-bound hydrogen atoms, which introduced an additional ammonia re-formation pathway through an Eley-Rideal reaction with gas phase NH2. Furthermore, a more uniform plasma is predicted in the presence of MgAl2O4, which leads to a lower average electron energy during micro-discharges and a lower conversion (37 %) at a comparable residence time for the plasma alone (51 %).

Abstract: We present a zero-dimensional kinetic model to characterise specifically the gas-phase dynamics of methane
conversion in a nanosecond pulsed discharge (NPD) plasma reactor. The model includes a systematic approach to
capture the nanoscale power discharges and the rapid ensuing changes in electric field, gas and electron temperature,
as well as species densities. The effects of gas temperature and reactor pressure on gas conversion and
product selectivity are extensively investigated and validated against experimental work. We discuss the
important reaction pathways and provide an analysis of the dynamics of the heating and cooling mechanisms. H
radicals are found to be the most populous plasma species and they participate in hydrogenation and dehydrogenation
reactions, which are the dominant recombination reactions leading to C2H4 and C2H2 as main
products (depending on the pressure).

Abstract: By combining dielectric barrier discharge plasma and external heating, we exploit a two-stage hybrid plasmathermal
system (HPTS), i.e., a plasma stage followed by a thermal stage, for direct non-oxidative coupling of
CH4 to C2H4 and H2, yielding a CH4 conversion of ca. 17 %. In the two-stage HPTS, the plasma first converts CH4
into C2H6 and C3H8, which in the thermal stage leads to a high C2H4 selectivity of ca. 63 % by pyrolysis, with H2
selectivity of ca. 64 %.

Abstract: We calculate bubble nucleation rates in a Lennard-Jones fluid through explicit molecular dynamics simulations. Our approach-based on a recent free energy method (dubbed reweighted Jarzynski sampling), transition state theory, and a simple recrossing correction-allows us to probe a fairly wide range of rates in several superheated and cavitation regimes in a consistent manner. Rate predictions from this approach bridge disparate independent literature studies on the same model system. As such, we find that rate predictions based on classical nucleation theory, direct brute force molecular dynamics simulations, and seeding are consistent with our approach and one another. Published rates derived from forward flux sampling simulations are, however, found to be outliers. This study serves two purposes: First, we validate the reliability of common modeling techniques and extrapolation approaches on a paradigmatic problem in materials science and chemical physics. Second, we further test our highly generic recipe for rate calculations, and establish its applicability to nucleation processes.