Abstract: Nitrogen oxyanions and oxyacids are important agents in atmospheric chemistry and medical biology. Although their chemical behavior in solution is relatively well understood, they may behave very differently at the water/air interface of atmospheric aerosols or at the membrane/water interface of cells. Here, we developed a fully classical model for molecular dynamics simulations of NO3−, NO2−, HNO3, and HNO2 in the framework of the GROMOS 53A6 and 54A7 force field versions. The model successfully accounted for the poorly structured solvation shell and ion pairing tendency of NO3−. Accurate pure-liquid properties and hydration free energies were obtained for the oxyacids. Simulations at the water/air interface showed a local enrichment of HNO3 and depletion of NO3−. The effect was discussed in light of earlier spectroscopic data and ab initio calculations, suggesting that HNO3 behaves as a weaker acid at the surface of water. Our model will hopefully allow for efficient and accurate simulations of nitrogen oxyanions and oxyacids in solution and at microheterogeneous interface environments.

Abstract: The hyperdynamics method is a powerful tool to simulate slow processes at the atomic level. However, the construction of an optimal hyperdynamics potential is a task that is far from trivial. Here, we propose a generally applicable implementation of the hyperdynamics algorithm, borrowing two concepts from metadynamics. First, the use of a collective variable (CV) to represent the accelerated dynamics gives the method a very large flexibility and simplicity. Second, a metadynamics procedure can be used to construct a suitable history-dependent bias potential on-the-fly, effectively turning the algorithm into a self-learning accelerated molecular dynamics method. This collective variable-driven hyperdynamics (CVHD) method has a modular design: both the local system properties on which the bias is based, as well as the characteristics of the biasing method itself, can be chosen to match the needs of the considered system. As a result, system-specific details are abstracted from the biasing algorithm itself, making it extremely versatile and transparent. The method is tested on three model systems: diffusion on the Cu(001) surface and nickel-catalyzed methane decomposition, as examples of reactive processes with a bond-length-based CV, and the folding of a long polymer-like chain, using a set of dihedral angles as a CV. Boost factors up to 109, corresponding to a time scale of seconds, could be obtained while still accurately reproducing correct dynamics.

Abstract: We introduce a simple enhanced sampling approach for the calculation of free energy differences and barriers along a one-dimensional reaction coordinate. First, a small number of short nonequilibrium simulations are carried out along the reaction coordinate, and the Jarzynski equality is used to learn an approximate free energy surface from the nonequilibrium work distribution. This free energy estimate is represented in a compact form as an artificial neural network and used as an external bias potential to accelerate rare events in a subsequent molecular dynamics simulation. The final free energy estimate is then obtained by reweighting the equilibrium probability distribution of the reaction coordinate sampled under the influence of the external bias. We apply our reweighted Jarzynski sampling recipe to four processes of varying scales and complexities.spanning chemical reaction in the gas phase, pair association in solution, and droplet nucleation in supersaturated vapor. In all cases, we find reweighted Jarzynski sampling to be a very efficient strategy, resulting in rapid convergence of the free energy to high precision.

Abstract: Thermal-catalytic gas processing is integral to many current industrial processes. Ever-increasing demands on conversion and energy efficiencies are a strong driving force for the development of alternative approaches. Similarly, synthesis of several functional materials (such as nanowires and nanotubes) demands special processing conditions. Plasma catalysis provides such an alternative, where the catalytic process is complemented by the use of plasmas that activate the source gas. This combination is often observed to result in a synergy between plasma and catalyst. This Review introduces the current state-of-the-art in plasma catalysis, including numerous examples where plasma catalysis has demonstrated its benefits or shows future potential, including CO2 conversion, hydrocarbon reforming, synthesis of nanomaterials, ammonia production, and abatement of toxic waste gases. The underlying mechanisms governing these applications, as resulting from the interaction between the plasma and the catalyst, render the process highly complex, and little is known about the factors leading to the often-observed synergy. This Review critically examines the catalytic mechanisms relevant to each specific application.

Abstract: The unique capabilities of microsecond dwell time (DT) single-particle inductively coupled plasma mass spectrometry (spICPMS) were utilized to characterize the cloud of ions generated from the introduction of suspensions of gold nanoparticles (AuNPs) into the plasma. A set of narrowly distributed particles with diameters ranging from 15.4 to 100.1 nm was synthesized and characterized according to established protocols. Statistically significant numbers of the short transient spICPMS events were evaluated by using 50 μs DT for their summed intensity, maximum intensity, and duration, of which all three were found to depend on the particle diameter. The summed intensity increases from 10 to 1661 counts and the maximum intensity from 6 to 309 counts for AuNPs with diameters from 15.4 to 83.2 nm. The event duration rises from 322 to 1007 μs upon increasing AuNP diameter. These numbers represent a comprehensive set of key data points of the ion clouds generated in ICPMS from AuNPs. The extension of event duration is of high interest to appoint the maximum possible particle number concentration at which separation of consecutive events in spICPMS can still be achieved. Moreover, the combined evaluation of all above-mentioned ion cloud characteristics can explain the regularly observed prolonged single-particle events. The transport and ionization behavior of AuNPs in the ICP was also computationally modeled to gain insight into the size-dependent signal generation. The simulated data reveals that the plasma temperature, and therefore the point of ionization of the particles, is the same for all diameters. However, the maximum number density of Au+, as well as the extent of the ion cloud, depends on the particle diameter, in agreement with the experimental data, and it provides an adequate explanation for the observed ion cloud characteristics.

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: The continuous miniaturization of nanodevices, such as transistors, solar cells, and optical fibers, requires the controlled synthesis of (ultra)thin gate oxides (<10 nm), including Si gate-oxide (SiO2) with high quality at the atomic scale. Traditional thermal growth of SiO2 on planar Si surfaces, however, does not allow one to obtain such ultrathin oxide due to either the high oxygen diffusivity at high temperature or the very low sticking ability of incident oxygen at low temperature. Two recent techniques, both operative at low (room) temperature, have been put forward to overcome these obstacles: (i) hyperthermal oxidation of planar Si surfaces and (ii) thermal or plasma-assisted oxidation of nonplanar Si surfaces, including Si nanowires (SiNWs). These nanooxidation processes are, however, often difficult to study experimentally, due to the key intermediate processes taking place on the nanosecond time scale.
In this Account, these Si nano-oxidation techniques are discussed from a computational point of view and compared to both hyperthermal and thermal oxidation experiments, as well as to well-known models of thermal oxidation, including the Deal−Grove, Cabrera−Mott, and Kao models and several alternative mechanisms. In our studies, we use reactive molecular dynamics (MD) and hybrid MD/Monte Carlo simulation techniques, applying the Reax force field. The incident energy of oxygen species is chosen in the range of 1−5 eV in hyperthermal oxidation of planar Si surfaces in order to prevent energy-induced damage. It turns out that hyperthermal growth allows for two growth modes, where the ultrathin oxide thickness depends on either (1) only the kinetic energy of the incident oxygen species at a growth temperature below Ttrans = 600 K, or (2) both the incident energy and the growth temperature at a growth temperature above Ttrans. These modes are specific to such ultrathin oxides, and are not observed in traditional thermal oxidation, nor theoretically considered by already existing models. In the case of thermal or plasma-assisted oxidation of small Si nanowires, on the other hand, the thickness of the ultrathin oxide is a function of the growth temperature and the nanowire diameter. Below Ttrans, which varies with the nanowire diameter, partially oxidized SiNW are formed, whereas complete oxidation to a SiO2 nanowire occurs only above Ttrans. In both nano-oxidation processes at lower temperature (T < Ttrans), final sandwich c-Si|SiOx|a-SiO2 structures are obtained due to a competition between overcoming the energy barrier to penetrate into Si subsurface layers and the compressive stress (∼2−3 GPa) at the Si crystal/oxide interface. The overall atomic-simulation results strongly indicate that the thickness of the intermediate SiOx (x < 2) region is very limited (∼0.5 nm) and constant irrespective of oxidation parameters. Thus, control over the ultrathin SiO2 thickness with good quality is indeed possible by accurately tuning the oxidant energy, oxidation temperature and surface curvature.
In general, we discuss and put in perspective these two oxidation mechanisms for obtaining controllable ultrathin gate-oxide films, offering a new route toward the fabrication of nanodevices via selective nano-oxidation.

Abstract: Nonthermal plasma is a promising alternative for ammonia synthesis at gentle conditions. Metal meshes of Fe, Cu, Pd, Ag, and Au were employed as catalysts in radio frequency plasma for ammonia synthesis. The energy yield for all these transition metal catalysts ranged between 0.12 and 0.19 g-NH3/kWh at 300 W and, thus, needs further improvement. In addition, a semimetal, pure gallium, was used for the first time as catalyst for ammonia synthesis, with energy yield of 0.22 g-NH3/kWh and with a maximum yield of ∼10% at 150 W. The emission spectra, as well as computer simulations, revealed hydrogen recombination as a primary governing parameter, which depends on the concentration or flux of H atoms in the plasma and on the catalyst surface. The simulations helped to elucidate the underlying mechanism, implicating the dominance of surface reactions and surface adsorbed species. The rate limiting step appears to be NH2 formation on the surface of the reactor wall and on the catalyst surface, which is different from classical catalysis.

Abstract: Binding of the SARS-CoV-2 S-glycoprotein to cell receptors is vital for the entry of the virus into cells and subsequent infection. ACE2 is the main cell receptor for SARS-CoV-2, which can attach to the C-terminal receptor-binding domain (RBD) of the SARS-CoV-2 S-glycoprotein. The GRP78 receptor plays an anchoring role, which attaches to the RBD and increases the chance of other RBDs binding to ACE2. Although high levels of reactive oxygen and nitrogen species (RONS) are produced during viral infections, it is not clear how they affect the RBD structure and its binding to ACE2 and GRP78. In this research, we apply molecular dynamics simulations to study the effect of oxidation of the highly reactive cysteine (Cys) amino acids of the RBD on its binding to ACE2 and GRP78. The interaction energy of both ACE2 and GRP78 with the whole RBD, as well as with the RBD main regions, is compared in both the native and oxidized RBDs. Our results show that the interaction energy between the oxidized RBD and ACE2 is strengthened by 155 kJ/mol, increasing the binding of the RBD to ACE2 after oxidation. In addition, the interaction energy between the RBD and GRP78 is slightly increased by 8 kJ/mol after oxidation, but this difference is not significant. Overall, these findings highlight the role of RONS in the binding of the SARS-CoV-2 S-glycoprotein to host cell receptors and suggest an alternative mechanism by which RONS could modulate the entrance of viral particles into the cells.

Abstract: Plasma degradation is one of the most effective methods for the abatement of greenhouse gas sulfur hexafluoride
(SF6). To evaluate the potential of MgO as a catalyst in plasma degradation, we investigate the catalytic properties
of MgO on SF6 adsorption and activation by density functional theory (DFT) where the O-defective and
hydroxylated surfaces are considered as two typical plasma-generated surfaces. Our results show that perfect
MgO (001) and (111) surfaces cannot interact with SF6 and only physical adsorption happens. In case of Odefective
MgO surfaces, the O vacancy is the most stable adsorption site. SF6 undergoes a decomposition to SF5
and F over the O-defective MgO (001) surface and undergoes an elongation of the bottom S-F bond over the Odefective
(111) surface. Besides, SF6 shows a physically adsorption at the stepsite of the MgO (001) surface,
accompanied by small changes in its bond angle and length. Furthermore, SF6 is found to be physically and
chemically adsorbed over 0.5 and 1.0 ML (monolayer) H-covered O-terminated MgO (111) surfaces, respectively.
The SF6 molecule undergoes a self-decomposition on the 1.0 ML hydroxylated surface via a surface bonding
process. This study shows that defective and hydroxylated MgO surfaces have the surface capacities for SF6
activation, which shows that MgO has potential as packing material in SF6 waste treatment in packed-bed
plasmas.

Abstract: Nowadays, lasers are commonly applied in spectrochemical analysis methods, for sampling, plasma formation or a combination of both. Despite the numerous investigations that have been performed on these applications, the underlying processes are still insufficiently understood. In order to fasten progress in the field and in honor of the lifework of professor Rick Russo, we here provide a brief overview of the fundamental mechanisms in lasermatter interaction as proposed in literature, and throw the spotlight on some aspects that have not received much attention yet. For an organized discussion, we choose laser ablation, laser desorption and the associated gaseous plasma formation as the central processes in this perspective article, based on a classification of the laserbased spectrochemical analysis techniques and the corresponding laser-matter interaction regimes. First, we put the looking glass over the excitation and thermalization mechanisms in the laser-irradiated condensed phase, for which we propose the so-called multi-plasma model. This novel model can be understood as an extension of the well-known two-temperature model, featuring multiple thermodynamic dimensions, each of which corresponds to a quasi-particle type. Next, the focus is placed on the mass transfer and ionization mechanisms, after which we shortly highlight the possible role of anisotropic and magnetic effects in the laser-excited material. We hope this perspective article motivates more fundamental research on laser-matter interaction, as a continuation of the lifework of Rick Russo.

Abstract: Melanoma arises from pigment-producing cells called melanocytes located in the basal layers of the epidermis of the skin. Cytoglobin (CYGB) is a ubiquitously expressed hexacoordinated globin that is highly enriched in me­lanocytes and frequently downregulated during melanomagenesis. Previously, we showed that non-thermal plasma (NTP)-produced reactive oxygen and nitrogen species (RONS) lead to the formation of an intra­ molecular disulfide bridge that would allow CYGB to function as a redox-sensitive protein. Here, we investigate the cytotoxic effect of indirect NTP treatment in two melanoma cell lines with divergent endogenous CYGB expression levels, and we explore the role of CYGB in determining treatment outcome. Our findings are consistent with previous studies supporting that NTP cytotoxicity is mediated through the production of RONS and leads to apoptotic cell death in melanoma cells. Furthermore, we show that NTP-treated solutions elicit an antioxidant response through the activation of nuclear factor erythroid 2–related factor 2 (NRF2). The knock­ down and overexpression of CYGB respectively sensitizes and protects melanoma cells from RONS-induced apoptotic cell death. The presence of CYGB enhances heme-oxygenase 1 (HO-1) and NRF2 protein expression levels, whereas the absence impairs their expression. Moreover, analysis of the CYGB-dependent transcriptome demonstrates the tumor suppressor long non-coding RNA maternally expressed 3 (MEG3) as a hitherto unde­ scribed link between CYGB and NRF2. Thus, the presence of CYGB, at least in melanoma cells, seems to play a central role in determining the therapeutic outcome of RONS-inducing anticancer therapies, like NTP-treated solutions, possessing both tumor-suppressive and oncogenic features. Hence, CYGB expression could be of in­ terest either as a biomarker or as a candidate for future targeted therapies in melanoma.

Abstract: Many current anti-cancer therapies rely on increasing the intracellular reactive oxygen and nitrogen species (RONS) contents with the aim to induce irreparable damage, which subsequently results in tumor cell death. A novel tool in cancer therapy is the use of cold atmospheric plasma (CAP), which has been found to be very effective in the treatment of many different cancer cell types in vitro as well as in vivo, mainly through the vast generation of RONS. One of the key determinants of the cell's fate will be the interaction of RONS, generated by CAP, with important proteins, i.e. redox-regulatory proteins. One such protein is cytoglobin (CYGB), a recently discovered globin proposed to be involved in the protection of the cell against oxidative stress. In this study, the effect of plasma-produced RONS on CYGB was investigated through the treatment of CYGB with CAP for different treatment times. Spectroscopic analysis of CYGB showed that although chemical modifications occur, its secondary structure remains intact. Mass spectrometry experiments identified these modifications as oxidations of mainly sulfur-containing and aromatic amino acids. With longer treatment time, the treatment was also found to induce nitration of the heme. Furthermore, the two surface-exposed cysteine residues of CYGB were oxidized upon treatment, leading to the formation of intermolecular disulfide bridges, and potentially also intramolecular disulfide bridges. In addition, molecular dynamics and docking simulations confirmed, and further show, that the formation of an intramolecular disulfide bond, due to oxidative conditions, affects the CYGB 3D structure, thereby opening the access to the heme group, through gate functioning of His117. Altogether, the results obtained in this study (1) show that plasma-produced RONS can extensively oxidize proteins and (2) that the oxidation status of two redox-active cysteines lead to different conformations of CYGB.

Abstract: Plasma-based gas conversion has great potential for enabling carbon-free fertilizer production powered by renewable electricity. Sustaining an energy-efficient plasma process without eroding the containment vessel is currently a significant challenge, limiting scaling to higher powers and throughputs. Isolation of the plasma from contact with any solid surfaces is an advantage, which both limits energy loss to the walls and prevents material erosion that could lead to disastrous soil contamination. This paper presents highly energy-efficient nitrogen fixation from air into NOx by microwave plasma, with the plasma filament isolated at the center of a quartz tube using a vortex gas flow. NOx production is found to scale very efficiently when increasing both gas flow rate and absorbed power. The lowest energy cost recorded of ~2 MJ/mol, for a total NOx production of ~3.8%, is the lowest reported up to now for atmospheric pressure plasmas.

Abstract: A falling film dielectric barrier discharge (DBD) plasma reactor combined with adsorption on activated carbon textile material was optimized to minimize the formation of hazardous oxidation by-products from the treatment of persistent pesticides (alachlor, diuron and isoproturon) in water. The formation of by-products and the reaction mechanism was investigated by HPLC-TOF-MS. The maximum concentration of each by-product was at least two orders of magnitude below the initial pesticide concentration, during the first 10 min of treatment. After 30 min of treatment, the individual by-product concentrations had decreased to values of at least three orders of magnitude below the initial pesticide concentration. The proposed oxidation pathways revealed five main oxidation steps: dechlorination, dealkylation, hydroxylation, addition of a double-bonded oxygen and nitrification. The latter is one of the main oxidation mechanisms of diuron and isoproturon for air plasma treatment. To our knowledge, this is the first time that the formation of nitrificated intermediates is reported for the plasma treatment of non-phenolic compounds.