Abstract: e investigate atmospheric pressure microwave (MW) plasma (2.45 GHz) conversion in CO2 and CH4 mixtures (i.e., dry reforming of methane, DRM) focusing on reaction performance and carbon formation. Promising energy costs of ~2.8–3.0 eV/molecule or ~11.1–11.9 kJ/L are amongst the best performance to date considering the current state-of-the-art for plasma-based DRM for all types of plasma. The conversion is in the range of ~46–49% and ~55–67% for CO2 and CH4, respectively, producing primarily syngas (i.e., H2 and CO) with H2/CO ratios of ~0.6–1 at CH4 fractions ranging from 30% to 45%. Water is the largest byproduct with levels ranging ~7–14% in the exhaust. Carbon particles visibly impact the plasma at higher CH4 fractions (> 30%), where they become heated and incandescent. Particle luminosity increases with increasing CH4 fractions, with the plasma becoming unstable near a 1:1 mixture (i.e., > 45% CH4). Electron microscopy of the carbon material reveals an agglomerated morphology of pure carbon nanoparticles. The mean particle size is determined as ~20 nm, free of any metal contamination, consistent with the electrode-less MW design.

Abstract: Uniform acceptance force biased Monte Carlo (UFMC) simulations have previously been shown to be a powerful tool to simulate atomic scale processes, enabling one to follow the dynamical path during the simulation. In this contribution, we present a simple proof to demonstrate that this uniform acceptance still complies with the condition of detailed balance, on the condition that the characteristic parameter lambda = 1/2 and that the maximum allowed step size is chosen to be sufficiently small. Furthermore, the relation to Metropolis Monte Carlo (MMC) is also established, and it is shown that UFMC reduces to MMC by choosing the characteristic parameter lambda = 0 [Rao, M. et al. Mol. Phys. 1979, 37, 1773]. Finally, a simple example compares the UFMC and MMC methods.

Abstract: Accurate barriers for rate controlling elementary reactions on metal surfaces are key to understanding, controlling, and predicting the rate of heterogeneously catalyzed processes. While barrier heights for gas phase reactions have been extensively benchmarked, dissociative chemisorption barriers for the reactions of molecules on metal surfaces have received much less attention. The first database called SBH10 and containing 10 entries was recently constructed based on the specific reaction parameter approach to density functional theory (SRP-DFT) and experimental results. We have now constructed a new and improved database (SBH17) containing 17 entries based on SRP-DFT and experiments. For this new SBH17 benchmark study, we have tested three algorithms (high, medium, and light) for calculating barrier heights for dissociative chemisorption on metals, which we have named for the amount of computational effort involved in their use. We test the performance of 14 density functionals at the GGA, GGA+vdW-DF, and meta-GGA rungs. Our results show that, in contrast with the previous SBH10 study where the BEEF-vdW-DF2 functional seemed to be most accurate, the workhorse functional PBE and the MS2 density functional are the most accurate of the GGA and meta-GGA functionals tested. Of the GGA+vdW functionals tested, the SRP32-vdW-DF1 functional is the most accurate. Additionally, we found that the medium algorithm is accurate enough for assessing the performance of the density functionals tested, while it avoids geometry optimizations of minimum barrier geometries for each density functional tested. The medium algorithm does require metal lattice constants and interlayer distances that are optimized separately for each functional. While these are avoided in the light algorithm, this algorithm is found not to give a reliable description of functional performance. The combination of relative ease of use and demonstrated reliability of the medium algorithm will likely pave the way for incorporation of the SBH17 database in larger databases used for testing new density functionals and electronic structure methods.

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: The evaluation of nucleation rates from molecular dynamics trajectories is hampered by the slow nucleation time scale and impact of finite size effects. Here, we show that accurate nucleation rates can be obtained in a very general fashion relying only on the free energy barrier, transition state theory, and a simple dynamical correction for diffusive recrossing. In this setup, the time scale problem is overcome by using enhanced sampling methods, in casu metadynamics, whereas the impact of finite size effects can be naturally circumvented by reconstructing the free energy surface from an appropriate ensemble. Approximations from classical nucleation theory are avoided. We demonstrate the accuracy of the approach by calculating macroscopic rates of droplet nucleation from argon vapor, spanning 16 orders of magnitude and in excellent agreement with literature results, all from simulations of very small (512 atom) systems.

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: The effect of gas heating in laboratory sputter magnetrons is investigated by means of numerical modeling. The model is two-dimensional in the coordinate space and three-dimensional in the velocity space based on the particle-in-cellMonte Carlo collisions technique. It is expanded in a way that allows the inclusion of the neutral plasma particles (fast gas atoms and sputtered atoms), which makes it possible to calculate the gas temperature and its influence on the discharge behavior in a completely self-consistent way. The results of the model are compared to experimental measurements and to other existing simulation results. The results show that gas heating is pressure dependent (rising with the increase in the gas pressure) and should be taken into consideration at pressures above 10 mTorr.

Abstract: The influence of the composition on the crystallographic properties of deposited Mg(M)O (with M=Al, Cr, Ti, Y, and Zr) films is studied. For a flexible control of the composition, dual reactive magnetron sputtering was used as deposition technique. Two different approaches to predict the composition are discussed. The first is an experimental way based on the simple relationship between the deposition rate and the target-substrate distance. The second is a route using a Monte Carlo based particle trajectory code. Both methods require a minimal experimental input and enable the user to quickly predict the composition of complex thin films. Good control and flexibility allow us to study the compositional effects on the growth of Mg(M)O films. Pure MgO thin films were grown with a (111) preferential out-of-plane orientation. When adding M to MgO, two trends were noticed. The first trend is a change in the MgO lattice parameters compared to pure MgO. The second tendency is a decrease in the crystallinity of the MgO phase. The experimentally determined crystallographic properties are shown to be in correspondence with the predicted properties from molecular dynamics simulations.

Abstract: In this paper, the negative ion behavior in a C4F8 inductively coupled plasma (ICP) is investigated using a hybrid model. The model predicts a non-monotonic variation of the total negative ion density with power at low pressure (1030 mTorr), and this trend agrees well with experiments that were carried out in many fluorocarbon (fc) ICP sources, like C2F6, CHF3, and C4F8. This behavior is explained by the availability of feedstock C4F8 gas as a source of the negative ions, as well as by the presence of low energy electrons due to vibrational excitation at low power. The maximum of the negative ion density shifts to low power values upon decreasing pressure, because of the more pronounced depletion of C4F8 molecules, and at high pressure (∼50 mTorr), the anion density continuously increases with power, which is similar to fc CCP sources. Furthermore, the negative ion composition is identified in this paper. Our work demonstrates that for a clear understanding of the negative ion behavior in radio frequency C4F8 plasma sources, one needs to take into account many factors, like the attachment characteristics, the anion composition, the spatial profiles, and the reactor configuration. Finally, a detailed comparison of our simulation results with experiments is conducted.