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: 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: 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: 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: 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: 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: 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 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: 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: 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: 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: 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: In this work, we use density functional theory calculations to study the combined effect of external electric fields, surface morphology, and surface charge on CO2 activation over Cu(111), Cu(211), Cu(110), and Cu(001) surfaces. We observe that the binding energy of the CO2 molecule on Cu surfaces increases significantly upon increasing the applied electric field strength. In addition, rougher surfaces respond more effectively to the presence of the external electric field toward facilitating the formation of a carbonate-like CO2 structure and the transformation of the most stable adsorption mode from physisorption to chemisorption. The presence of surface charges further strengthens the electric field effect and consequently causes an improved bending of the CO2 molecule and C−O bond length elongation. On the other hand, a net charge in the absence of an externally applied electric field shows only a marginal effect on CO2 binding. The chemisorbed CO2 is more stable and further activated when the effects of an external electric field, rough surface, and surface charge are combined. These results can help to elucidate the underlying factors that control CO2 activation in heterogeneous and plasma catalysis, as well as in electrochemical processes.

Abstract: Plasma-catalytic ammonia synthesis is receiving ever increasing attention, especially in packed bed dielectric barrier discharge (DBD) reactors. The latter typically operate in the filamentary regime when used for gas conversion applications. While DBDs are in principle well understood and already applied in the industry, the incorporation of packing materials and catalytic surfaces considerably adds to the complexity of the plasma physics and chemistry governing the ammonia formation. We employ a plasma kinetics model to gain insights into the ammonia formation mechanisms, paying special attention to the role of filamentary microdischarges and their afterglows. During the microdischarges, the synthesized ammonia is actually decomposed, but the radicals created upon electron impact dissociation of N2 and H2 and the subsequent catalytic reactions cause a net ammonia gain in the afterglows of the microdischarges. Under our plasma conditions, electron impact dissociation of N2 in the gas phase followed by the adsorption of N atoms is identified as a rate-limiting step, instead of dissociative adsorption of N2 on the catalyst surface. Both elementary Eley−Rideal and Langmuir−Hinshelwood reaction steps can be found important in plasma-catalytic NH3 synthesis.

Abstract: To gain full control over the growth of carbon nanotubes (CNTs) using plasma-enhanced chemical vapor deposition (PECVD), a thorough understanding of the underlying plasma-catalyst mechanisms is required. Oxygen-containing species are often used as or added to the growth precursor gas, but these species also yield various radicals and ions, which may simultaneously etch the CNT during the growth. At present, the effect of these reactive species on the growth onset has not yet been thoroughly investigated. We here report on the etching mechanism of incipient CNT structures from OH and O radicals as derived from combined (reactive) molecular dynamics (MD) and force-bias Monte Carlo (tfMC) simulations. Our results indicate that the oxygen-containing radicals initiate a dissociation process. In particular, we show how the oxygen species weaken the interaction between the CNT and the nanocluster. As a result of this weakened interaction, the CNT closes off and dissociates from the cluster in the form of a fullerene. Beyond the specific systems studied in this work, these results are generically important in the context of PECVD-based growth of CNTs using oxygen-containing precursors.