Abstract: The atmospheric pressure dielectric barrier discharge in helium is a pulsed discharge in nature. If during the electrical current pulse a glow discharge is reached, then this pulse will last only a few microseconds in operating periods of sinusoidal voltage with lengths of about 10 to 100 µs. In this paper we demonstrate that right before a glow discharge is reached, the discharge very closely resembles the commonly assumed Townsend discharge structure, but actually contains some significant differing features and hence should not be considered as a Townsend discharge. In order to clarify this, we present calculation results of high time and space resolution of the pulse formation. The results indicate that indeed a maximum of ionization is formed at the anode, but that the level of ionization remains high and that the electric field at that time is significantly disturbed. Our results also show where this intermediate structure comes from.

Abstract: Plasma modelling is an exciting subject in which virtually all physical disciplines are represented. Plasma models combine the electromagnetic, statistical and fluid dynamical theories that have their roots in the 19th century with the modern insights concerning the structure of matter that were developed throughout the 20th century. The present cluster issue consists of 20 invited contributions, which are representative of the state of the art in plasma modelling and numerical simulation. These contributions provide an in-depth discussion of the major theories and modelling and simulation strategies, and their applications to contemporary plasma-based technologies. In this editorial review, we introduce and complement those papers by providing a bird's eye perspective on plasma modelling and discussing the historical context in which it has surfaced.

Abstract: Current-voltage analysis using different optical band pass filters has been performed on Cu2ZnSnSe4 and Cu2ZnSn(S, Se)(4) thin-film solar cells. When using red or orange light (i.e. wavelengths above 600 nm), a distortion appears in the I-V curve of the Cu2ZnSnSe4 solar cell, indicating an additional potential barrier to the current flow in the device for these conditions of illumination. This barrier is reduced when using a Cu2ZnSn(S, Se)(4) absorber. Numerical simulations demonstrate that the barrier visible under red light could be explained by a positive conduction band offset at the front interface coupled with compensating defects in the buffer layer.

Abstract: While plasma treatment of skin diseases and wound healing has been proven highly effective, the underlying mechanisms, and more generally the effect of plasma radicals on skin tissue, are not yet completely understood. In this paper, we perform ReaxFF-based reactive molecular dynamics simulations to investigate the interaction of plasma generated OH radicals with a model system composed of free fatty acids, ceramides, and cholesterol molecules. This model system is an approximation of the upper layer of the skin (stratum corneum). All interaction mechanisms observed in our simulations are initiated by H-abstraction from one of the ceramides. This reaction, in turn, often starts a cascade of other reactions, which eventually lead to the formation of aldehydes, the dissociation of ceramides or the elimination of formaldehyde, and thus eventually to the degradation of the skin barrier function.

Abstract: The size and lattice constant evolution of Pb nanoparticles (NPs) synthesized by high fluence implantation in crystalline Si have been studied with a variety of experimental techniques. Results obtained from small-angle x-ray scattering showed that the Pb NPs grow with increasing implantation fluence and annealing duration. The theory of NP growth kinetics can be applied to qualitatively explain the size evolution of the Pb NPs during the implantation and annealing processes. Moreover, the lattice constant of the Pb NPs was evaluated by conventional x-ray diffraction. The lattice dilatation was observed to decrease with increasing size of the Pb NPs. Such lattice constant tuning can be attributed to the pseudomorphism caused by the lattice mismatch between the Pb NPs and the Si matrix.

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.