Abstract: Secondary electron emission is an important process that plays a significant role in several plasma‐related applications. As measuring the secondary electron yield experimentally is very challenging, quantitative modelling of this process to obtain reliable yield data is critical as input for higher‐scale simulations. Here, we build upon our previous work combining density functional theory calculations with a model originally developed by Hagstrum to extend its application to metallic surfaces. As plasmonic effects play a much more important role in the secondary electron emission mechanism for metals, we introduce an approach based on Poisson point processes to include both surface and bulk plasmon excitations to the process. The resulting model is able to reproduce the yield spectra of several available experimental results quite well but requires the introduction of global fitting parameters, which describe the strength of the plasmon interactions. Finally, we use an in‐house developed workflow to calculate the electron yield for a list of elemental surfaces spanning the periodic table to produce an extensive data set for the community and compare our results with more simplified approaches from the literature.

Abstract: Plasma-based dry reforming of methane (DRM) into
high-value-added oxygenates is an appealing approach to enable
otherwise thermodynamically unfavorable chemical reactions at
ambient pressure and near room temperature. However, it suffers
from coke deposition due to the deep decomposition of CH4. In this
work, we assess the DRM performance upon O2 addition, as well as
varying temperature, CO2/CH4 ratio, discharge power, and gas
residence time, for optimizing oxygenate production. By adding O2,
the main products can be shifted from syngas (CO + H2) toward
oxygenates. Chemical kinetics modeling shows that the improved
oxygenate production is due to the increased concentration of
oxygen-containing radicals, e.g., O, OH, and HO2, formed by electron
impact dissociation [e + O2 → e + O + O/O(1D)] and subsequent
reactions with H atoms. Our study reveals the crucial role of oxygen-coupling in DRM aimed at oxygenates, providing practical
solutions to suppress carbon deposition and at the same time enhance the oxygenates production in plasma-assisted DRM.

Abstract: The DFT-1/2 method is a band-gap correction with GW precision at a density functional theory (DFT) computational cost. The method was also extended to correct the gap between defect levels, allowing for the calculation of optical transitions. However, this method fails when the atomic character of the occupied and unoccupied defect levels is similar as we illustrate by two examples, the tetrahedral hydrogen interstitial and the negatively charged vacancy in diamond. We solve this problem by decoupling the effect of the occupied and unoccupied defect levels and call this the decoupled DFT-1/2 method for defects.

Abstract: technologies have been expanding, and one of the most exciting and rapidly growing
applications is in biology and medicine. Most biomedical studies with DBD plasma systems are performed in vitro, which include cells grown on the surface of plastic well plates, or in vivo, which include animal research models (e.g. mice, pigs). Since many DBD systems use the biological target as the secondary electrode for direct plasma generation and treatment, they are sensitive to the surface properties of the target, and thus can be altered based on the in vitro or in vivo system used. This could consequently affect biological response from plasma treatment. Therefore, in this study, we investigated the DBD plasma behavior both in vitro (i.e. 96-well flat bottom plates, 96-well U-bottom plates, and 24-well flat bottom plates), and in vivo (i.e. mouse skin). Intensified charge coupled device (ICCD) imaging was performed and the plasma discharges were visually distinguishable between the different systems. The geometry of the wells did not affect DBD plasma generation for low application distances (≤ 2 mm), but differentially affected plasma uniformity on the bottom of the well at greater distances. Since DBD plasma treatment in vitro is rarely performed in dry wells for plasma medicine experiments, the effect of well wetness was also investigated. In all in vitro cases, the uniformity of the DBD plasma was affected when comparing wet versus dry wells, with the plasma in the wide-bottom wells appearing the most similar to plasma generated on mouse skin. Interestingly, based on quantification of ICCD images, the DBD plasma intensity per surface area demonstrated an exponential one-phase decay with increasing application distance, regardless of the in vitro or in vivo system. This trend is similar to that of the energy per pulse of plasma, which is used to determine the total plasma treatment energy for biological systems. Optical emission spectroscopy performed on the plasma revealed similar trends in radical species generation between the plastic well plates and mouse skin. Therefore, taken together, DBD plasma intensity per surface area may be a valuable parameter to be used as a simple method for in situ monitoring during biological treatment and active plasma treatment control, which can be applied for in vitro and in vivo systems.

Abstract: In this study, the multi-layer self-assembled molecular network formation of an alkylated tetrathiafulvalene compound is studied at the liquid-solid interface between 1-phenyloctane and graphite. A combined theoretical/experimental approach associating force-field and quantum-chemical calculations with scanning tunnelling microscopy is used to determine the two-dimensional self-assembly beyond the monolayer, but also to further the understanding of the molecular adsorption conformation and its impact on the molecular packing within the assemblies at the monolayer and bilayer level.

Abstract: Chirality in gold nanostructures offers an exciting opportunity to tune their differential optical response to left- and right-handed circularly polarized light, as well as their interactions with biomolecules and living matter. However, tuning and understanding such interactions demands quantification of the structural features that are responsible for the chiral behavior. Electron tomography (ET) enables structural characterization at the single-particle level and has been used to quantify the helicity of complex chiral nanorods. However, the technique is time-consuming and consequently lacks statistical value. To address this issue, we introduce herein a high-throughput methodology that combines images acquired by secondary electron-based electron beam-induced current (SEEBIC) with quantitative image analysis. As a result, the geometric chirality of hundreds of nanoparticles can be quantified in less than 1 h. When combining the drastic gain in data collection efficiency of SEEBIC with a limited number of ET data sets, a better understanding of how the chiral structure of individual chiral nanoparticles translates into the ensemble chiroptical response can be reached.

Abstract: Recent experiments have revealed that the macroscopic Kelvin equation remains surprisingly accurate even for nanoscale capillaries. This phenomenon was so far explained by the oscillatory behavior of the solid−liquid interfacial free energy. We here demonstrate thermodynamic and capillarity inconsistencies with this explanation. After revising the Kelvin equation, we ascribe its validity at nanoscale confinement to the effect of disjoining pressure.
To substantiate our hypothesis, we employed molecular dynamics simulations to evaluate interfacial heat transfer and wetting properties. Our assessments unveil a breakdown in a previously established proportionality between the work of adhesion and the Kapitza conductance at capillary heights below 1.3 nm, where the dominance of the work of adhesion shifts primarily from energy to entropy. Alternatively, the peak density of the initial water layer can effectively probe the work of adhesion. Unlike under bulk conditions, high confinement renders the work of adhesion entropically unfavorable.

Abstract: We have developed a comprehensive kinetic model to study the O atom kinetics in an O₂plasma and its afterglow. By adopting a pseudo-1D plug-flow formalism within the kinetic model, our aim is to assess how far the O atoms travel in the plasma afterglow, evaluating its potential as a source of O atoms for post-plasma gas conversion applications. Since we could not find experimental data for pure O₂plasma at atmospheric pressure, we first validated our model at low pressure (1–10 Torr) where very good experimental data are available. Good agreement between our model and experiments was achieved for the reduced electric field, gas temperature and the densities of the dominant neutral species, i.e. O₂(a), O₂(b) and O. Subsequently, we confirmed that the chemistry set is consistent with thermodynamic equilibrium calculations at atmospheric pressure. Finally, we investigated the O atom densities in the O₂plasma and its afterglow, for which we considered a microwave O₂plasma torch, operating at a pressure between 0.1 and 1 atm, for a flow rate of 20 slm and an specific energy input of 1656 kJ mol⁻¹. Our results show that for both pressure conditions, a high dissociation degree of ca. 92% is reached within the discharge. However, the O atoms travel much further in the plasma afterglow for<italic>p</italic>= 0.1 atm (9.7 cm) than for<italic>p</italic>= 1 atm (1.4 cm), attributed to the longer lifetime (3.8 ms at 0.1 atm vs 1.8 ms at 1 atm) resulting from slower three-body recombination kinetics, as well as a higher volumetric flow rate.

Abstract: Plasma‐based CO₂conversion has attracted increasing interest. However, to understand the impact of plasma operation on post‐plasma processes, we studied the effect of adding N₂, N₂/CH₄and N₂/CH₄/H₂O to a CO₂gliding arc plasmatron (GAP) to obtain valuable insights into their impact on exhaust stream composition and temperature, which will serve as feed gas and heat for post‐plasma catalysis (PPC). Adding N₂improves the CO₂conversion from 4 % to 13 %, and CH₄addition further promotes it to 44 %, and even to 61 % at lower gas flow rate (6 L/min), allowing a higher yield of CO and hydrogen for PPC. The addition of H₂O, however, reduces the CO₂conversion from 55 % to 22 %, but it also lowers the energy cost, from 5.8 to 3 kJ/L. Regarding the temperature at 4.9 cm post‐plasma, N₂addition increases the temperature, while the CO₂/CH₄ratio has no significant effect on temperature. We also calculated the temperature distribution with computational fluid dynamics simulations. The obtained temperature profiles (both experimental and calculated) show a decreasing trend with distance to the exhaust and provide insights in where to position a PPC bed.

Abstract: A combination of a gliding arc plasmatron (GAP) reactor and a newly designed tubular catalyst bed (N-bed) was applied to investigate the post-plasma catalytic (PPC) effect for dry reforming of methane (DRM). As comparison, a traditional plasma catalyst bed (T-bed) was also utilized. The post-plasma catalytic effect of a Ni-based mixed oxide (Ni/MO) catalyst with a thermal catalytic performance of 77% CO2 and 86% CH4 conversion at 700 ℃ was studied. Although applying the T-bed had little effect on plasma based CO2 and CH4 conversion, an increase in selectivity to H2 was obtained with a maximum value of 89% at a distance of 2 cm. However, even when only α-Al2O3 packing material was used in the N-bed configuration, compared to the plasma alone and the T-bed, an increase of the CO2 and CH4 conversion from 53% and 53% to 69% and 69% to 83% was achieved. Addition of the Ni/MO catalyst further enhanced the DRM reaction, resulting in conversions of 79% for CO2 and 91% for
CH4. Hence, although no insulation nor external heating was applied to the N-bed post plasma, it provides a slightly better conversion than the thermal catalytic performance with the same catalyst, while being fully electrically driven. In addition, an enhanced CO selectivity to 96% was obtained and the energy cost was reduced from ~ 6 kJ/L (plasma alone) to 4.3 kJ/L. To our knowledge, it is the first time that a post-plasma catalytic system achieves this excellent catalytic performance for DRM without extra external heating or insulation.

Abstract: The basal plane of graphene can function as a selective barrier that is permeable to protons but impermeable to all ions and gases, stimulating its use in applications such as membranes, catalysis and isotope separation. Protons can chemically adsorb on graphene and hydrogenate it, inducing a conductor–insulator transition that has been explored intensively in graphene electronic devices. However, both processes face energy barriersand various strategies have been proposed to accelerate proton transport, for example by introducing vacancies, incorporating catalytic metalsor chemically functionalizing the lattice. But these techniques can compromise other properties, such as ion selectivity or mechanical stability. Here we show that independent control of the electric field,<italic>E</italic>, at around 1 V nm⁻¹, and charge-carrier density,<italic>n</italic>, at around 1 × 10¹⁴ cm⁻², in double-gated graphene allows the decoupling of proton transport from lattice hydrogenation and can thereby accelerate proton transport such that it approaches the limiting electrolyte current for our devices. Proton transport and hydrogenation can be driven selectively with precision and robustness, enabling proton-based logic and memory graphene devices that have on–off ratios spanning orders of magnitude. Our results show that field effects can accelerate and decouple electrochemical processes in double-gated 2D crystals and demonstrate the possibility of mapping such processes as a function of<italic>E</italic>and<italic>n</italic>, which is a new technique for the study of 2D electrode–electrolyte interfaces.

Abstract: In order to model correctly the corrosion of spent nuclear fuel under disposal conditions, it is important to understand its behavior in the presence of oxidants. To advance in this direction, we consider the oxidation of UO2. We investigate computationally the adsorption of various species on its three most stable surfaces: (111), (110), and (100), with emphasis on incorporating a full non-collinear PBE+U approach. Various species, namely O, O2, H2O and H2O2 are considered due to their relevance for the oxidation of UO2. The dissociation energy and an estimate for the dissociation barrier for O2 were obtained, using the preferred adsorption configurations of O and O2. The adsorption configurations for H2O in our study compare well with previous studies that used collinear approximations, both in terms of relative stability of configurations and bond lengths. Differences in adsorption energies were found, which may be important for reaction kinetics. Dissociative reactions in which the water molecule splits in hydrogen and hydroxyl occur only on one of the three surfaces. The hydrogen further reacts with a surface oxygen to also form a hydroxyl group. Not surprisingly, we find that H2O2 binds more strongly to the three surfaces than water (lower formation energy), and similar to H2O adsorption, dissociative reactions may occur. The dissociated hydrogen reacts with a surface oxygen to form a hydroxyl group and the hydroperoxyl molecule binds with a surface uranium. Our study, which includes a detailed study of electron transfer, magnetic structure and the preferred adsorption configurations, gives insight into the uranium oxidation states and the influence of surface geometry on adsorption. The findings contribute to a more comprehensive understanding of the early stages of UO2 oxidation.

Abstract: Plasma catalysis is an emerging technology, but a lot of questions about the underlying surface mechanisms remain unanswered. One of these questions is how important Eley−Rideal (ER) reactions are, next to Langmuir−Hinshelwood reactions. Most plasma catalysis kinetic models predict ER reactions to be important and sometimes even vital for the surface chemistry. In this work, we take a critical look at how ER reactions involving H radicals are incorporated in kinetic models describing CO2 hydrogenation and NH3 synthesis. To this end, we construct potential energy surface (PES) intersections, similar to elbow plots constructed for dissociative chemisorption. The results of the PES intersections are in agreement with ab initio molecular dynamics (AIMD) findings in literature while being computationally much cheaper. We find that, for the reactions studied here, adsorption is more probable than a reaction via the hot atom (HA) mechanism, which in turn is more probable than a reaction via the ER mechanism. We also conclude that kinetic models of plasma-catalytic systems tend to overestimate the importance if ER reactions. Furthermore, as opposed to what is often assumed in kinetic models, the choice of catalyst will influence the ER reaction probability. Overall, the description of ER reactions is too much “ideal” in models. Based on our indings, we make a number of recommendations on how to incorporate ER reactions in kinetic models to avoid overestimation of their importance.

Abstract: A gliding arc (GA) plasma, operating at atmospheric pressure in a gas mixture of 50% N₂and 50% O₂, is studied using laser-induced fluorescence spectroscopy. The main goal is to determine the two-dimensional distribution of both the gas temperature and the NO ground state density in the afterglow. As GA plasma discharges at atmospheric pressure normally produce rather high NO_{<italic>x</italic>}densities, the high concentration of relevant absorbers, such as NO, may impose essential restrictions for the use of ‘classical’ laser-induced fluorescence methods (dealing with excitation in the bandhead vicinity), as the laser beam would be strongly absorbed along its propagation in the afterglow. Since this was indeed the case for the studied discharge, an approach dealing with laser-based excitation of separate rotational lines is proposed. In this case, due to a non-saturated absorption regime, simultaneous and reliable measurements of both the NO density and the gas temperature (using a reference fitting spectrum) are possible. The proposed method is applied to provide a two-dimensional map for both the NO density and the gas temperature at different plasma conditions. The results show that the input gas flow rate strongly alters the plasma shape, which appears as an elongated column at low input gas flow rate and spreads laterally as the flow rate increases. Finally, based on temperature map analysis, a clear correlation between the gas temperature and NO concentration is found. The proposed method may be interesting for the plasma-chemical analysis of discharges with high molecular production yields, where knowledge of both molecular concentration and gas temperature is required.