Abstract: Atomistic simulation methods for the quantification of free energies are in wide use. These methods operate by sampling the probability density of a system along a small set of suitable collective variables (CVs), which is, in turn, expressed in the form of a free energy surface (FES). This definition of the FES can capture the relative stability of metastable states but not that of the transition state because the barrier height is not invariant to the choice of CVs. Free energy barriers therefore cannot be consistently computed from the FES. Here, we present a simple approach to calculate the gauge correction necessary to eliminate this inconsistency. Using our procedure, the standard FES as well as its gauge-corrected counterpart can be obtained by reweighing the same simulated trajectory at little additional cost. We apply the method to a number of systems—a particle solvated in a Lennard-Jones fluid, a Diels–Alder reaction, and crystallization of liquid sodium—to demonstrate its ability to produce consistent free energy barriers that correctly capture the kinetics of chemical or physical transformations, and discuss the additional demands it puts on the chosen CVs. Because the FES can be converged at relatively short (sub-ns) time scales, a free energy-based description of reaction kinetics is a particularly attractive option to study chemical processes at more expensive quantum mechanical levels of theory.

Abstract: Plasma has been studied for several years to convert CO2 into value-added products. If CO2 could be converted in the presence of H2O as a cheap H-source for making syngas and oxygenates, it would mimic natural photosynthesis. However, CO2/H2O plasmas have not yet been extensively studied, not by experiments, and certainly not computationally. Therefore, we present here a kinetic modelling study to obtain a greater understanding of the vibrational kinetics of a CO2/H2O microwave plasma. For this purpose, we first created an electron impact cross section set for H2O, using a swarm-derived method. We added the new cross section set and CO2/H2O-related chemistry to a pure CO2 model. While it was expected that H2O addition mainly causes quenching of the CO2 asymmetric mode vibrational levels due to the additional CO2/H2O vibrational-translational relaxation, our model shows that the modifications in the vibrational kinetics are mainly induced by the strong electron dissociative attachment to H2O molecules, causing a reduction in electron density, and the corresponding changes in the input of energy into the CO2 vibrational levels by electron impact processes.

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: The dielectric properties of confined water is of fundamental interest and is still controversial. For water confined in channels with height smaller than h = 8 angstrom, we found a commensurability effect and an extraordinary decrease in the out-of-plane dielectric constant down to the limit of the dielectric constant of optical water. Spatial resolved polarization density data obtained from molecular dynamics simulations are found to be antisymmetric across the channel and are used as input in a mean-field model for the dielectric constant as a function of the height of the channel for h > 15 angstrom. Our results are in excellent agreement with a recent experiment [L. Fumagalli et al., Science 360, 1339 (2018)].

Abstract: We study the photoionization of the s states in the systems bound by sufficiently weak central fields V(r) for the large photon energies corresponding to the relativistic photoelectrons. We demonstrate that the energy dependence of the photoionization cross section can be obtained without solving the wave equation. We show that the shape of the energy dependence of the cross section is determined by analytical properties of the binding potential V(r). We find the cross sections for the potentials V(r) which have singularities in the origin, on the real axis and in the complex plane.

Abstract: There is increasing interest in plasma technology for CO2 conversion because it can operate at mild conditions and it can store fluctuating renewable electricity into
value-added compounds and renewable fuels. This perspective paper aims to provide a view on the future for non-specialists who want to understand the role of plasma
technology in the new scenario for sustainable and low-carbon energy and chemistry. Thus, it is prepared to give a personal view on future opportunities and challenges. First, we introduce the current state-of-the-art and the potential of plasma-based CO2 conversion. Subsequently, we discuss the challenges to overcome the current limitations and to apply plasma technology on a large scale. The final section discusses the general context and the potential benefits of plasma-based CO2 conversion for our life and the impact on climate change. It also includes a brief analysis on the future scenario for energy and chemical production, and how plasma technology may realize new paths for CO2 utilization.

Abstract: In this work, the effects of N-doping into the Co-doped single vacancy (Co-SV-G) and di-vacancy graphene flake (Co-dV-G) are investigated and compared toward direct oxidation of methane to methanol (DOMM) employing two different oxidants (N2O and O-2) using density functional theory (DFT) calculation. We found that DOMM on CoN3-G utilizing the N2O molecule as oxygen-donor proceeds via a two-step reaction with low activation energies. In addition, we found that although CoN3-G might be a good catalyst for methane conversion, it can also catalyze the oxidation of methanol to CO2 and H2O due to the required low activation barriers. Moreover, the adsorption behaviors of CHx (x = 0-4) species and dehydrogenation of CHx (x = 1-4) species on CoN3-G are investigated. We concluded that CoN3-G can be used as an efficient catalyst for DOMM and N-2O reduction at ambient conditions which may serve as a guide for fabricating effective C/N catalysts in energy-related devices.

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: 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: Positron annihilation spectroscopy in the Doppler and coincidence Doppler mode was applied on Ge1 xSnx epitaxial layers, grown by chemical vapor deposition with different total As concentrations (1019-1021 cm3), high active As concentrations (1019 cm3), and similar Sn concentrations (5.9%-6.4%). Positron traps are identified as mono-vacancy complexes. Vacancy-As complexes, V-Asi, formed during the growth were studied to deepen the understanding of the electrical passivation of the Ge1 xSnx:As epilayers. Larger monovacancy complexes, V-Asi (i 2), are formed as the As doping increases. The total As concentration shows a significant impact on the saturation of the number of As atoms (i 1/4 4) around the vacancies in the sample epilayers. The presence of V-Asi complexes decreases the dopant activation in the Ge1 xSnx:As epilayers. Furthermore, the presence of Sn failed to hinder the formation of larger V-Asi complexes and thus failed to reduce the donor-deactivation.

Abstract: Natural deep eutectic solvents (NADES) are mixtures of naturally derived compounds with a significantly decreased melting point due to the specific interactions among the constituents. NADES have benign properties (low volatility, flammability, toxicity, cost) and tailorable physicochemical properties (by altering the type and molar ratio of constituents), hence they are often considered as a green alternative to common organic solvents. Modeling the relation between their composition and properties is crucial though, both for understanding and predicting their behavior. Several efforts were done to this end, yet this review aims at structuring the present knowledge as an outline for future research. First, we reviewed the key properties of NADES and relate them to their structure based on the available experimental data. Second, we reviewed available modeling methods applicable to NADES. At the molecular level, density functional theory and molecular dynamics allow interpreting density differences and vibrational spectra, and computation of interaction energies. Additionally, properties at the level of the bulk media can be explained and predicted by semi-empirical methods based on ab initio methods (COSMO-RS) and equation of state models (PC-SAFT). Finally, methods based on large datasets are discussed; models based on group contribution methods and machine learning. A combination of bulk media and dataset modeling allows qualitative prediction and interpretation of phase equilibria properties on the one hand, and quantitative prediction of melting point, density, viscosity, surface tension and refractive indices on the other hand. In our view, multiscale modeling, combining the molecular and macroscale methods, will strongly enhance the predictability of NADES properties and their interaction with solutes, yielding truly tailorable solvents to accommodate (bio)chemical reactions.

Abstract: This paper benchmarks various epitaxial growth schemes based on n-type group-IV materials as viable source/drain candidates for Ge nMOS devices. Si:P grown at low temperature on Ge, gives an active carrier concentration as high as 3.5 x 10(20) cm(-3) and a contact resistivity down to 7.5 x 10(-9) Omega.cm(2). However, Si:P growth is highly defective due to large lattice mismatch between Si and Ge. Within the material stacks assessed, one option for Ge nMOS source/drain stressors would be to stack Si:P, deposited at contact level, on top of a selectively grown n-SiyGe1-x-ySnx at source/drain level, in line with the concept of Si passivation of n-Ge surfaces to achieve low contact resistivities as reported in literature (Martens et al. 2011 Appl. Phys. Lett., 98, 013 504). The saturation in active carrier concentration with increasing P (or As)-doping is the major bottleneck in achieving low contact resistivities for as-grown Ge or SiyGe1-x-ySnx. We focus on understanding various dopant deactivation mechanisms in P-doped Ge and Ge1-xSnx alloys. First principles simulation results suggest that P deactivation in Ge and Ge1-xSnx can be explained both by P-clustering and donor-vacancy complexes. Positron annihilation spectroscopy analysis, suggests that dopant deactivation in P-doped Ge and Ge1-xSnx is primarily due to the formation of P-n-V and SnmPn-V clusters. (C) 2020 The Author(s). Published on behalf of The Electrochemical Society by IOP Publishing Limited.

Abstract: This study focuses on a magnetically stabilized gliding arc (MGA) plasma. Two fully coupled flow-plasma models (in 3D and 2D) are presented. The 3D model is applied to compare the arc dynamics of the MGA with a traditional gas-driven gliding arc. The 2D model is used for a detailed parametric study on the effect of the external magnetic field. The results show that the relative velocity between the plasma and feed gas is generated due to the Lorentz force, which can increase the plasma-treated gas fraction. The magnetic field also helps to decrease the gas temperature by enhancing heat transfer and to increase the electron number density. This work shows the potential of an external magnetic field to control the gliding arc behavior, for enhanced gas conversion at low gas flow rates.

Abstract: Plasma is gaining interest for CH4 conversion into higher hydrocarbons and H2. However, the performance in terms of conversion and selectivity toward different hydrocarbons is different for different plasma types, and the underlying mechanisms are not yet fully understood. Therefore, we study here these mechanisms in different plasma sources, by means of a chemical kinetics model. The model is first validated by comparing the calculated conversions and hydrocarbon/H2 selectivities with experimental results in these different plasma types and over a wide range of specific energy input (SEI) values. Our model predicts that vibrational−translational nonequilibrium is negligible in all CH4 plasmas investigated, and instead, thermal conversion is important. Higher gas temperatures also lead to a more selective production of unsaturated hydrocarbons (mainly C2H2) due to neutral dissociation of CH4 and subsequent dehydrogenation processes, while three-body recombination reactions into saturated hydrocarbons (mainly C2H6, but also higher hydrocarbons) are dominant in low temperature plasmas.

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.