Abstract: Experiments have demonstrated the potential selective anticancer capacity of cold atmospheric plasmas (CAPs), but the underlying mechanisms remain unclear. Using computer simulations, we try to shed light on the mechanism of selectivity, based on aquaporins (AQPs), i.e. transmembrane protein channels transferring external H 2 O 2 and other reactive oxygen species, created e.g. by CAPs, to the cell interior. Specifically, we perform molecular dynamics simulations for the permeation of H 2 O 2 through AQP1 (one of the members of the AQP family) and the palmitoyl-oleoyl-phosphatidylcholine (POPC) phospholipid bilayer (PLB). The free energy barrier of H 2 O 2 across AQP1 is lower than for the POPC PLB, while the permeability coefficient, calculated using the free energy and diffusion rate profiles, is two orders of magnitude higher. This indicates that the delivery of H 2 O 2 into the cell interior should be through AQP. Our study gives a better insight into the role of AQPs in the selectivity of CAPs for treating cancer cells.

Abstract: Cold atmospheric plasmas (CAPs) have attracted significant interest for their potential benefits in medical applications, including cancer therapy. The therapeutic effects of CAPs are related to reactive oxygen and nitrogen species (ROS and RNS) present in the plasma. The impact of ROS has been extensively studied, but the role of RNS in CAP-treatment remains poorly understood at the molecular level. Here, we investigate the permeation of RNS and ROS across native and oxidized phospholipid bilayers (PLBs) by means of computer simulations. The results reveal significantly lower free energy barriers for RNS (i.e. NO, NO2, N2O4) and O3 compared to hydrophilic ROS, such as OH, HO2 and H2O2. This suggests that the investigated RNS and O3 can permeate more easily through both native and oxidized PLBs in comparison to hydrophilic ROS, indicating their potentially important role in plasma medicine.

Abstract: A thorough understanding of the skyrmion motion through nanotracks is a prerequisite to realize the full potential of spintronic applications like the skyrmion racetrack memory. One of the challenges is to place the data, i.e. skyrmions, on discrete fixed positions, e.g. below a read or write head. In the domain-wall racetrack memory, one proposed solution to this problem was patterning the nanotrack with notches. Following this approach, this paper reports on the skyrmion mobility through a nanotrack with periodic notches (constrictions) made using variations in the chiral Dzyaloshinskii-Moriya interaction. We observe that such notches induce a coupling between the mobility and the skyrmion breathing mode, which manifests itself as velocity-dependent oscillations of the skyrmion diameter and plateaus in which the velocity is independent of the driving force. Despite the fact that domain walls are far more rigid objects than skyrmions, we were able to perform an analogous study and, surprisingly, found even larger plateaus of constant velocity. For both systems it is straightforward to tune the velocity at these plateaus by changing the design of the notched nanotrack geometry, e.g. by varying the distance between the notches. Therefore, the notch-induced coupling between the excited modes and the mobility could offer a strategy to stabilize the velocity against unwanted perturbations in racetrack-like applications. In the last part of the paper we focus on the low-current mobility regimes, whose very rich dynamics at nonzero temperatures are very similar to the operating principle of recently developed probabilistic logic devices. This proves that the mobility of nanomagnetic structures through a periodically modulated track is not only interesting from a fundamental point of view, but has a future in many spintronic applications.

Abstract: Plasma catalysis is gaining increasing interest for various gas conversion applications, such as CO2 conversion into value-added chemicals and fuels, CH4 activation into hydrogen, higher hydrocarbons or oxygenates, and NH3 synthesis. Other applications are already more established, such as for air pollution control, e.g. volatile organic compound remediation, particulate matter and NOx removal. In addition, plasma is also very promising for catalyst synthesis and treatment. Plasma catalysis clearly has benefits over ‘conventional’ catalysis, as outlined in the Introduction. However, a better insight into the underlying physical and chemical processes is crucial. This can be obtained by experiments applying diagnostics, studying both the chemical processes at the catalyst surface and the physicochemical mechanisms of plasma-catalyst interactions, as well as by computer modeling. The key challenge is to design cost-effective, highly active and stable catalysts tailored to the plasma environment. Therefore, insight from thermal catalysis as well as electro- and photocatalysis is crucial. All these aspects are covered in this Roadmap paper, written by specialists in their field, presenting the state-of-the-art, the current and future challenges, as well as the advances in science and technology needed to meet these challenges.

Abstract: In this work, we extensively investigate the structural and electronic properties of van der Waals heterostructures (HTs) constructed by MoS${2}$/$BC₃$, MoS${2}$/$C₃N$, MoS${2}$/$C₃N₄$, MoS${2}$/$C₄N₃$ and those using Janus MoSSe instead of MoS$2$ by performing density functional theory calculations. The electronic band structure calculations and the corresponding partial density of states reveal that the significant changes are driven by quite strong layer-layer interaction between the constitutive layers. Our results show that although all monolayers are semiconductors as free-standing layers, the MoS${2}$/$C₃N$ and MoS${2}$/$C₄N₃$ bilayer HTs display metallic behavior as a consequence of transfer of charge carriers between two constituent layers. In addition, it is found that in MoSSe/$C₃N$ bilayer HT, the degree of metallicity is affected by the interface chalcogen atom type when Se atoms are facing to $C₃N$ layer, the overlap of the bands around the Fermi level is smaller. Moreover, the half-metallic magnetic $C₄N₃$ is shown to form magnetic half-metallic trilayer HT with MoS$2$ independent of the stacking sequence, i.e. whether it is sandwiched or two $C₄N₃$ layer encapsulate MoS$2$ layer. We further analyze the trilayer HTs in which MoS$2$ is encapsulated by two different monolayers and it is revealed that at least with one magnetic monolayer, it is possible to construct a magnetic trilayer. While the trilayer of $C₄N₃$/MoS${2}$/$BC₃$ and $C₄N₃$/MoS${2}$/$C₃N₄$ exhibit half-metallic characteristics, $C₄N₃$/MoS${_2}$/$C₃$N possesses a magnetic metallic ground state. Overall, our results reveal that holly structures of BCN crystals are suitable for heterostructure formation even over van der Waals type interaction which significantly changes electronic nature of the constituent layers.

Abstract: The rate, selectivity and efficiency of plasma-based conversion processes is strongly affected by nonequilibrium phenomena. High concentrations of vibrationally excited molecules are such a plasma-induced effect. It is frequently assumed that vibrationally excited molecules are important in plasma catalysis because their presence lowers the apparent activation energy of dissociative chemisorption reactions and thus increases the conversion rate. A detailed atomic-level understanding of vibrationally stimulated catalytic reactions in the context of plasma catalysis is however lacking. Here, we couple a recently developed statistical model of a plasma-induced vibrational nonequilibrium to molecular dynamics simulations, enhanced sampling methods, and machine learning techniques. We quantify the impact of a vibrational nonequilibrium on the dissociative chemisorption barrier of H2 and CH4 on nickel catalysts over a wide range of vibrational temperatures. We investigate the effect of surface structure and compare the role of different vibrational modes of methane in the dissociation process. For low vibrational temperatures, very high vibrational efficacies are found, and energy in bend vibrations appears to dominate the dissociation of methane. The relative impact of vibrational nonequilibrium is much higher on terrace sites than on surface steps. We then show how our simulations can help to interpret recent experimental results, and suggest new paths to a better understanding of plasma catalysis.

Abstract: The 2022 Roadmap is the next update in the series of Plasma Roadmaps published by<italic>Journal of Physics</italic>D with the intent to identify important outstanding challenges in the field of low-temperature plasma (LTP) physics and technology. The format of the Roadmap is the same as the previous Roadmaps representing the visions of 41 leading experts representing 21 countries and five continents in the various sub-fields of LTP science and technology. In recognition of the evolution in the field, several new topics have been introduced or given more prominence. These new topics and emphasis highlight increased interests in plasma-enabled additive manufacturing, soft materials, electrification of chemical conversions, plasma propulsion, extreme plasma regimes, plasmas in hypersonics, data-driven plasma science and technology and the contribution of LTP to combat COVID-19. In the last few decades, LTP science and technology has made a tremendously positive impact on our society. It is our hope that this roadmap will help continue this excellent track record over the next 5–10 years.

Abstract: The crystal structure of the newly discovered double-Hg-layer copper oxide superconductor (Hg, Pr)(2)Ba-2(Y, Ca)Cu2O8-delta was studied as a function of chemical doping using neutron and electron diffraction and high-resolution transmission electron microscopy (HREM). Rietveld refinements of the structural parameters from neutron powder diffraction data indicate that the oxygen site O3 on the mercury plane is highly defective, being both partially occupied and displaced from the high-symmetry position. The variable concentration of oxygen vacancies partially compensates for the cation doping and, together with the O3 displacement field, makes some of the Hg atoms acquire an unusual pyramidal coordination. HREM images confirm that the structure is of the '2212' type, with very few defects. In some grains, faint superstructure reflections were evidenced by electron diffraction, suggesting that both the oxygen vacancies and the O3 displacement field may order at least on a local scale.

Abstract: The calculation of defect levels in semiconductors within a density functional theory approach suffers greatly from the band gap problem. We propose a band gap correction scheme that is based on the separation of energy differences in electron addition and relaxation energies. We show that it can predict defect levels with a reasonable accuracy, particularly in the case of defects with conduction band character, and yet is simple and computationally economical. We apply this method to ZnO doped with group III elements (Al, Ga, In). As expected from experiment, the results indicate that Zn substitutional doping is preferred over interstitial doping in Al, Ga, and In-doped ZnO, under both zinc-rich and oxygen-rich conditions. Further, all three dopants act as shallow donors, with the +1 charge state having the most advantageous formation energy. Also, doping effects on the electronic structure of ZnO are sufficiently mild so as to affect little the fundamental band gap and lowest conduction bands dispersion, which secures their n-type transparent conducting behavior. A comparison with the extrapolation method based on LDA+U calculations and with the HeydScuseriaErnzerhof hybrid functional (HSE) shows the reliability of the proposed scheme in predicting the thermodynamic transition levels in shallow donor systems.

Abstract: Seed-mediated growth is the most efficient methodology to control the size and shape of colloidal metal nanoparticles. In this process, the final nanocrystal shape is defined by the crystalline structure of the initial seed as well as by the presence of ligands and other additives that help to stabilize certain crystallographic facets. We analyze here the growth mechanism in aqueous solution of silver shells on presynthesized gold nanoparticles displaying various well-defined crystalline structures and morphologies. A thorough three-dimensional electron microscopy characterization of the morphology and internal structure of the resulting core-shell nanocrystals indicates that {100} facets are preferred for the outer silver shell, regardless of the morphology and crystallinity of the gold cores. These results are in agreement with theoretical analysis based on the relative surface energies of the exposed facets in the presence of halide ions.

Abstract: Recent findings revealed that surface oxygen can participate in the oxygen evolution reaction (OER) for the most active catalysts, which eventually triggers a new mechanism for which the deprotonation of surface intermediates limits the OER activity. We propose in this work a “dual strategy” in which tuning the electronic properties of the oxide, such as La1-xSrxCoO3-delta, can be dissociated from the use of surface functionalization with phosphate ion groups (P-i) that enhances the interfacial proton transfer. Results show that the P-i functionalized La0.5Sr0.5CoO3-delta gives rise to a significant enhancement of the OER activity when compared to La0.5Sr0.5Co3-delta and LaCoO3. We further demonstrate that the P-i surface functionalization selectivity enhances the activity when the OER kinetics is limited by the proton transfer. Finally, this work suggests that tuning the catalytic activity by such a “dual approach” may be a new and largely unexplored avenue for the design of novel high-performance catalysts.

Abstract: We present an approach to incorporate the effect of vibrational nonequilibrium in molecular dynamics (MD) simulations. A perturbed canonical ensemble, in which selected modes are excited to higher temperature while all others remain equilibrated at low temperature, is simulated by applying a specifically tailored bias potential. Our method can be readily applied to any (classical or quantum mechanical) MD setup at virtually no additional computational cost and allows the study of reactions of vibrationally excited molecules in nonequilibrium environments such as plasmas. In combination with enhanced sampling methods, the vibrational efficacy and mode selectivity of vibrationally stimulated reactions can then be quantified in terms of chemically relevant observables, such as reaction rates and apparent free energy barriers. We first validate our method for the prototypical hydrogen exchange reaction and then show how it can capture the effect of vibrational excitation on a symmetric SN2 reaction and radical addition on CO2.

Abstract: Water permeation between stacked layers of hBN sheets forming 2D nanochannels is investigated using large-scale ab initio-quality molecular dynamics simulations. A high-dimensional neural network potential trained on density-functional theory calculations is employed. We simulate water in van der Waals nanocapillaries and study the impact of nanometric confinement on the structure and dynamics of water using both equilibrium and nonequilibrium methods. At an interlayer distance of 10.2 A confinement induces a first-order phase transition resulting in a well-defined AA-stacked bilayer of hexagonal ice. In contrast, for h < 9 A, the 2D water monolayer consists of a mixture of different locally ordered patterns of squares, pentagons, and hexagons. We found a significant change in the transport properties of confined water, particularly for monolayer water where the water-solid friction coefficient decreases to half and the diffusion coefficient increases by a factor of 4 as compared to bulk water. Accordingly, the slip-velocity is found to increase under confinement and we found that the overall permeation is dominated by monolayer water adjacent to the hBN membranes at extreme confinements. We conclude that monolayer water in addition to bilayer ice has a major contribution to water transport through 2D nanochannels.

Abstract: A multiscale modeling and simulation approach, including first-principles calculations, ab initio molecular dynamics simulations, and a tight binding approach, is employed to study band flattening of the electronic band structure of oxidized monolayer graphene. The width offlat bands can be tuned by strain, the external electric field, and the density of functional groups and their distribution. A transition to a conducting state is found for monolayer graphene with impurities when it is subjected to an electric field of similar to 1.0 V/angstrom. Several parallel impurity-induced flat bands appear in the low-energy spectrum of monolayer graphene when the number of epoxy groups is changed. The width of the flat band decreases with an increase in tensile strain but is independent of the electric field strength. Here an alternative and easy route for obtaining band flattening in thermodynamically stable functionalized monolayer graphene is introduced. Our work discloses a new avenue for research on band flattening in monolayer graphene.