Abstract: In this study, we investigated the diffusion of H-atoms to the subsurface and their further diffusion into the bulk of a Ni(111) crystal by means of density functional theory calculations in the context of thermal and plasma-assisted catalysis. The H-atoms at the surface can originate from the dissociative adsorption of H2 or CH4 molecules, determining the surface H-coverage. When a threshold H-coverage is passed, corresponding to 1.00 ML for the crystalline Ni(111) surface, the surface-bound H-atoms start to diffuse to the subsurface. A similar threshold coverage is observed for the interstitial H-coverage. Once the interstitial sites are filled up with a coverage above 1.00 ML of H, dissolution of interstitial H-atoms to the layer below the interstitial sites will be initiated. Hence, by applying a high pressure or inducing a reactive plasma and high temperature, increasing the H-flux to the surface, a large amount of hydrogen can diffuse in a crystalline metal like Ni and can be absorbed. The formation of metal hydride may modify the entire reaction kinetics of the system. Equivalently, the H-atoms in the bulk can easily go back to the surface and release a large amount of heat. In a plasma process, H-atoms are formed in the plasma, and therefore the energy barrier for dissociative adsorption is dismissed, thus allowing achievement of the threshold coverage without applying a high pressure as in a thermal process. As a result, depending on the crystal plane and type of metal, a large number of H-atoms can be dissolved (absorbed) in the metal catalyst, explaining the high efficiency of plasma-assisted catalytic reactions. Here, the mechanism of H-dissolution is established as a chemical identifier for the investigation of the reaction kinetics of a chemical process.

Abstract: The interaction between cold atmospheric pressure plasma and liquids is receiving increasing attention for various applications. In particular, the use of plasma-treated liquids (PTL) for biomedical applications is of growing importance, in particular for sterilization and cancer treatment. However, insight into the
underlying mechanisms of plasma–liquid interactions is still scarce. Here, we present a 2D fluid dynamics model for the interaction between a plasma jet and liquid water. Our results indicate that the formed reactive species originate from either the gas phase (with further solvation) or are formed at the liquid interface. A clear increase in the aqueous density of H2O2, HNO2/NO2- and NO3-
is observed as a function of time, while the densities of O3, HO2/O2- and ONOOH/ONOO- are found to quickly reach a maximum due to chemical reactions in solution. The trends observed in our model correlate well with experimental observations from the literature.

Abstract: Numerous applications have required the study of CO2 plasmas since the 1960s, from CO2 lasers to spacecraft heat shields. However, in recent years, intense research activities on the subject have restarted because of environmental problems associated with CO2 emissions. The present review provides a synthesis of the current state of knowledge on the physical chemistry of cold CO2 plasmas. In particular, the different modeling approaches implemented to address specific aspects of CO2 plasmas are presented. Throughout the paper, the importance of conducting joint experimental, theoretical and modeling studies to elucidate the complex couplings at play in CO2 plasmas is emphasized. Therefore, the experimental data that are likely to bring relevant constraints to the different modeling approaches are first reviewed. Second, the calculation of some key elementary processes obtained with semi-empirical, classical and quantum methods is presented. In order to describe the electron kinetics, the latest coherent sets of cross section satisfying the constraints of “electron swarm” analyses are introduced, and the need for self-consistent calculations for determining accurate electron energy distribution function (EEDF) is evidenced. The main findings of the latest zero-dimensional (0D) global models about the complex chemistry of CO2 and its dissociation products in different plasma discharges are then given, and full state-to-state (STS) models of only the vibrational-dissociation kinetics developed for studies of spacecraft shields are described. Finally, two important points for all applications using CO2 containing plasma are discussed: the role of surfaces in contact with the plasma, and the need for 2D/3D models to capture the main features of complex reactor geometries including effects induced by fluid dynamics on the plasma properties. In addition to bringing together the latest advances in the description of CO2 non-equilibrium plasmas, the results presented here also highlight the fundamental data that are still missing and the possible routes that still need to be investigated.

Abstract: Ammonia is an industrial large-volume chemical, with its main application in fertilizer production. It also attracts increasing attention as a green-energy vector. Over the past century, ammonia production has been dominated by the Haber-Bosch process, in which a mixture of nitrogen and hydrogen gas is converted to ammonia at high temperatures and pressures. Haber-Bosch processes with natural gas as the source of hydrogen are responsible for a significant share of the global CO(2)emissions. Processes involving plasma are currently being investigated as an alternative for decentralized ammonia production powered by renewable energy sources. In this work, we present the PNOCRA process (plasma nitrogen oxidation and catalytic reduction to ammonia), combining plasma-assisted nitrogen oxidation and lean NO(x)trap technology, adopted from diesel-engine exhaust gas aftertreatment technology. PNOCRA achieves an energy requirement of 4.6 MJ mol(-1)NH(3), which is more than four times less than the state-of-the-art plasma-enabled ammonia synthesis from N(2)and H(2)with reasonable yield (>1 %).

Abstract: In this contribution, we present an overview of the various techniques for combining atomistic molecular dynamics with Monte Carlo simulations, mainly in the context of condensed matter systems, as well as a brief summary of the main accelerated dynamics techniques. Special attention is given to the force bias Monte Carlo technique and its combination with molecular dynamics, in view of promising recent developments, including a definable timescale. Various examples of the application of combined molecular dynamics / Monte Carlo simulations are given, in order to demonstrate the enhanced simulation efficiency with respect to either pure molecular dynamics or Monte Carlo.

Abstract: Pancreatic ductal adenocarcinoma (PDAC) is a challenging neoplastic disease, mainly due to the development of resistance to radio- and chemotherapy. Cold atmospheric plasma (CAP) is an alternative technology that can eliminate cancer cells through oxidative damage, as shown in vitro, in ovo, and in vivo. However, how CAP affects the pancreatic stellate cells (PSCs), key players in the invasion and metastasis of PDAC, is poorly understood. This study aims to determine the effect of an anti-PDAC CAP treatment on PSCs tissue developed in ovo using mono- and co-cultures of RLT-PSC (PSCs) and Mia PaCa-2 cells (PDAC). We measured tissue reduction upon CAP treatment and mRNA expression of PSC activation markers and extracellular matrix (ECM) remodelling factors via qRT-PCR. Protein expression of selected markers was confirmed via immunohistochemistry. CAP inhibited growth in Mia PaCa-2 and co-cultured tissue, but its effectiveness was reduced in the latter, which correlates with reduced ki67 levels. CAP did not alter the mRNA expression of PSC activation and ECM remodelling markers. No changes in MMP2 and MMP9 expression were observed in RLT-PSCs, but small changes were observed in Mia PaCa-2 cells. Our findings support the ability of CAP to eliminate PDAC cells, without altering the PSCs.

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: The effect of pore size on plasma catalysis is crucial but still unclear. Studies have shown plasma cannot enter micropores and mesopores, so catalysts for traditional thermocatalysis may not fit plasma catalysis. Here, 3D porous Cu and CuO with different pore sizes were prepared using uniform silica particles (10–2000 nm) as templates, and compared in plasma-catalytic dry reforming. In most cases, the smaller the pore size, the higher the conversion of CH4 and CO2. Large pores reachable by more electrons did not improve the reaction efficiency. We attribute this to the small surface area and large crystallite size, as indicated by N2-sorption, mercury intrusion and XRD. While the smaller pores might not be reachable by electrons, due to the sheath formed in front of them, as predicted by modeling, they can still be reached by radicals formed in the plasma, and ions can even be attracted into these pores. An exception are the samples synthesized from 1 μm silica, which show better performance. We believe this is due to the electric field enhancement for pore sizes close to the Debye length. The performances of CuO and Cu with different pore sizes can provide references for future research on oxide supports and metal components of plasma catalysts.

Abstract: A gliding arc plasma is a potential way to convert CO2 into CO and O2, due to its non-equilibrium character, but little is known about the underlying mechanisms. In this paper, a self-consistent two-dimensional (2D) gliding arc model is developed, with a detailed non-equilibrium CO2 plasma chemistry, and validated with experiments. Our calculated values of the electron number density in the plasma, the CO2 conversion and energy efficiency show reasonable agreement with the experiments, indicating that the model can provide a realistic picture of the plasma chemistry. Comparison of the results with classical thermal conversion, as well as other plasma-based technologies for CO2 conversion reported in literature, demonstrates the non-equilibrium character of the gliding arc, and indicates that the gliding arc is a promising plasma reactor for CO2 conversion. However, some process modifications should be exploited to further improve its performance. As the model provides a realistic picture of the plasma behaviour, we use it first to investigate the plasma characteristics in a whole gliding arc cycle, which is necessary to understand the underlying mechanisms. Subsequently, we perform a chemical kinetics analysis, to investigate the different pathways for CO2 loss and formation. Based on the revealed discharge properties and the underlying CO2 plasma chemistry, the model allows us to propose solutions on how to further improve the
CO2 conversion and energy efficiency by a gliding arc plasma.

Abstract: A packed bed dielectric barrier discharge (DBD) is widely used for plasma catalysis applications, but the exact plasma characteristics in between the packing beads are far from understood. Therefore, we study here these plasma characteristics by means of fluid modelling and experimental observations using ICCD imaging, for packing materials with different dielectric constants. Our study reveals that a packed bed DBD reactor in dry air at atmospheric pressure may show three types of discharges, i.e. positive restrikes, filamentary microdischarges, which can also be localized between two packing beads, and surface discharges (so-called surface ionization
waves). Restrikes between the dielectric surfaces result in the formation of filamentary microdischarges, while surface charging creates electric field components parallel to the dielectric surfaces, leading to the formation of surface discharges. A transition in discharge mode occurs from surface discharges to local filamentary discharges between the packing beads when the dielectric constant of the packing rises from 5 to 1000. This may have implications for the efficiency of plasma catalytic gas treatment, because the catalyst activation may be limited by constraining the discharge to the contact points of the beads. The production of reactive species occurs most in the positive restrikes, the surface discharges and the local microdischarges in between the beads, and is less significant in the longer filamentary microdischarges. The faster streamer propagation and discharge development with higher dielectric constant of the packing beads leads to a faster production of reactive species. This study is of great interest for plasma catalysis, where packing beads with different dielectric constants are often used as supports for the catalytic materials. It allows us to better understand how different packing materials can influence the performance of packed bed plasma reactors for environmental applications.

Abstract: We propose a novel Power to Synthesis Gas (P2SG) approach, composed of two high-efficiency and renewable electricity-driven units, i.e., plasma catalytic reforming (PCR) and water electrolysis (WE), to produce high quality syngas from CH4, CO2 and H2O. As WE technology is already commercial, we mainly focus on the PCR unit, consisting of gliding arc plasma and Ni-based catalyst, for oxidative dry reforming of methane. An energy efficiency of 78.9% and energy cost of 1.0 kWh/Nm3 at a CH4 conversion of 99% and a CO2 conversion of 79% are obtained. Considering an energy efficiency of 80% for WE, the P2SG system yields an overall energy efficiency of 79.3% and energy cost of 1.8 kWh/Nm3. High-quality syngas is produced without the need for posttreatment units, featuring the ideal stoichiometric number of 2, with concentration of 94.6 vol%, and a desired CO2 fraction of 1.9 vol% for methanol synthesis. The PCR unit has the advantage of fast response to adapting to fluctuation of renewable electricity, avoiding local hot spots in the catalyst bed and coking, in contrast to conventional catalytic processes. Moreover, pure O2 from the WE unit is directly utilized by the PCR unit for oxidative dry reforming of methane, and thus, no air separation unit, like in conventional processes, is required. This work demonstrates the viability of the P2SG approach for large-scale energy storage of renewable electricity via electricity-to-fuel conversion.

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: 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: 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: Direct conversion of CH4 and CO2 to liquid fuels and chemicals under mild conditions is appealing for biogas conversion and utilization but challenging due to the inert nature of both gases. Herein, we report a promising plasma process for the catalyst-free single-step conversion of CH4 and CO2 into higher value oxygenates (i.e., methanol, acetic acid, ethanol, and acetone) at ambient pressure and room temperature using a water-cooled dielectric barrier discharge (DBD) reactor, with methanol being the main liquid product. The distribution of liquid products could be tailored by tuning the discharge power, reaction temperature and residence time. Lower discharge powers (10–15 W) and reaction temperatures (5–20 ◦ C) were favourable for the production of liquid products, achieving the highest methanol selectivity of 43% at 5 ◦ C and 15 W. A higher discharge power and reaction temperature, on the other hand, produced more gaseous products, particularly H2 (up to 26% selec­tivity) and CO (up to 33% selectivity). In addition, varying these process parameters (discharge power, reaction temperature and residence time) resulted in a simultaneous change in key discharge properties, such as mean electron energy (Ee), electron density (ne) and specific energy input (SEI), all of which are essential determiners of plasma chemical reactions. According to the results of artificial neural network (ANN) models, the relative importance of these process parameters and key discharge indicators on reaction performance follows the order: discharge power > reaction temperature > residence time, and SEI > ne > Ee, respectively. This work provides new insights into the contributions and tuning mechanism of multiple parameters for optimizing the reaction performance (e.g., liquid production) in the plasma gas conversion process.