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: In this review paper, an overview is given of different modelling efforts for plasmas used for the formation and growth of nanostructured materials. This includes both the plasma chemistry, providing information on the precursors for nanostructure formation, as well as the growth processes itself. We limit ourselves to carbon (and silicon) nanostructures. Examples of the plasma modelling comprise nanoparticle formation in silane and hydrocarbon plasmas, as well as the plasma chemistry giving rise to carbon nanostructure formation, such as (ultra)nanocrystalline diamond ((U)NCD) and carbon nanotubes (CNTs). The second part of the paper deals with the simulation of the (plasma-based) growth mechanisms of the same carbon nanostructures, i.e. (U)NCD and CNTs, both by mechanistic modelling and detailed atomistic simulations.

Abstract: We investigate the theoretical possibility of single layer graphene formation on a nickel surface at different substrate temperatures under far-from-equilibrium high precursor flux conditions, employing state-of-the-art hybrid reactive molecular dynamics/uniform acceptance force bias Monte Carlo simulations. It is predicted that under these conditions, the formation of a single layer graphene-like film may proceed through a combined depositionsegregation mechanism on a nickel substrate, rather than by pure surface segregation as is typically observed for metals with high carbon solubility. At 900 K and above, nearly continuous graphene layers are obtained. These simulations suggest that single layer graphene deposition is theoretically possible on Ni under high flux conditions.

Abstract: This work examines the properties of a dielectric barrier discharge (DBD) reactor, built for CO2 decomposition, by means of electrical characterization, optical emission spectroscopy and gas chromatography. The discharge, formed in an electronegative gas (such as CO2, but also O2), exhibits clearly different electrical characteristics, depending on the surface conductivity of the reactor walls. An asymmetric current waveform is observed in the metaldielectric (MD) configuration, with sparse high-current pulses in the positive half-cycle (HC) and a more uniform regime in the negative HC. This indicates that the discharge is operating in two alternating regimes with rather different properties. At high CO2 conversion regimes, a conductive coating is deposited on the dielectric. This so-called coated MD configuration yields a symmetric current waveform, with current peaks in both the positive and negative HCs. In a double-dielectric (DD) configuration, the current waveform is also symmetric, but without current peaks in both the positive and negative HC. Finally, the DD configuration with conductive coating on the inner surface of the outer dielectric, i.e. so-called coated DD, yields again an asymmetric current waveform, with current peaks in the negative HC. These different electrical characteristics are related to the presence of the conductive coating on the dielectric wall of the reactor and can be explained by an increase of the local barrier capacitance available for charge transfer. The different discharge regimes affect the CO2 conversion, more specifically, the CO2 conversion is lowest in the clean DD configuration. It is somewhat higher in the coated DD configuration, and still higher in the MD configuration. The clean and coated MD configuration, however, gave similar CO2 conversion. These results indicate that the conductivity of the dielectric reactor walls can highly promote the development of the high-amplitude discharge current pulses and subsequently the CO2 conversion.

Abstract: In this work, a chemical kinetics study on methane activation for hydrogen production in a warm plasma, i.e., N-2 rotating gliding arc (RGA), was performed for the first time to get new insights into the underlying reaction mechanisms and pathways. A zero-dimensional chemical kinetics model was developed, which showed a good agreement with the experimental results in terms of the conversion of CH4 and product selectivities, allowing us to get a better understanding of the relative significance of various important species and their related reactions to the formation and loss of CH4, H-2, and C2H2 etc. An overall reaction scheme was obtained to provide a realistic picture of the plasma chemistry. The results reveal that the electrons and excited nitrogen species (mainly N-2(A)) play a dominant role in the initial dissociation of CH4. However, the H atom induced reaction CH4+ H -> CH3+ H-2, which has an enhanced reaction rate due to the high gas temperature (over 1200 K), is the major contributor to both the conversion of CH4 and H-2 production, with its relative contributions of > 90% and > 85%, respectively, when only considering the forward reactions. The coexistence and interaction of thermochemical and plasma chemical processes in the rotating gliding arc warm plasma significantly enhance the process performance. The formation of C-2 hydrocarbons follows a nearly one-way path of C2H6 -> C2H4 -> C2H2, explaining why the selectivities of C-2 products decreased in the order of C2H2 > C2H4 > C2H6.

Abstract: The electron bounce resonance heating (BRH) in dual-frequency capacitively coupled plasmas operated in oxygen and argon has been studied by different experimental methods. In comparison with the electropositive argon discharge, the BRH in an electronegative discharge occurs at larger electrode gaps. Kinetic particle simulations reveal that in the oxygen discharge, the bulk electric field becomes quite strong and is out of phase with the sheath field. Therefore, it retards the resonant electrons when traversing the bulk, resulting in a suppressed BRH. This effect becomes more pronounced at lower high-frequency power, when the discharge mode changes from electropositive to electronegative.

Abstract: Uniform-acceptance force-bias Monte Carlo (fbMC) methods have been shown to be a powerful technique to access longer timescales in atomistic simulations allowing, for example, phase transitions and growth. Recently, a new fbMC method, the time-stamped force-bias Monte Carlo (tfMC) method, was derived with inclusion of an estimated effective timescale; this timescale, however, does not seem able to explain some of the successes the method. In this contribution, we therefore explicitly quantify the effective timescale tfMC is able to access for a variety of systems, namely a simple single-particle, one-dimensional model system, the Lennard-Jones liquid, an adatom on the Cu(100) surface, a silicon crystal with point defects and a highly defected graphene sheet, in order to gain new insights into the mechanisms by which tfMC operates. It is found that considerable boosts, up to three orders of magnitude compared to molecular dynamics, can be achieved for solid state systems by lowering of the apparent activation barrier of occurring processes, while not requiring any system-specific input or modifications of the method. We furthermore address the pitfalls of using the method as a replacement or complement of molecular dynamics simulations, its ability to explicitly describe correct dynamics and reaction mechanisms, and the association of timescales to MC simulations in general.

Abstract: In this paper, two types of melanoma and glioblastoma cancer cell lines are treated with cold atmospheric plasma to assess the effect of several parameters on the cell viability. The cell viability decreases with treatment duration and time until analysis in all cell lines with varying sensitivity. The majority of dead cells stains both AnnexinV (AnnV) and propidium iodide, indicating that the plasma-treated non-viable cells are mostly late apoptotic or necrotic. Genetic mutations might be involved in the response to plasma. Comparing the effects of two gas mixtures, as well as indirect plasma-activated medium versus direct treatment, gives different results per cell line. In conclusion, this study confirms the potential of plasma for cancer therapy and emphasizes the influence of experimental parameters on therapeutic outcome.

Abstract: The laser ablation-induced plasma was used as a composition-con trolled source for ion implantation in Si crystals. Then, laser ablation in combination with inductively coupled plasma mass spectrometry was used for the elemental depth profiling of the implanted samples. Monte Carlo simulations permitted us to conclude that a depth resolution of tens of nm would be necessary to define the shape of the implantation profiles, as is obtained using XPS and RBS, whereas a hundred nm depth resolution is sufficient to determine the total implanted dose. The detection power of LA-ICP-MS would routinely allow rapid analytical control on the trace level implanted dose. Nevertheless, this technique is limited in terms of depth profiling resolution due to pulse mixing and signal tailing induced during the aerosol transport. Raw signal processing procedures were developed for the minimization of shapeline dispersion, deconvolution of pulse mixing and more appropriate assessment of the implanted profiles. Shapeline dispersion could be corrected for by determining the signal waning constant and implementing this information for a non-affine alibi transformation of the LA-ICP-MS signal traces. Pulse mixing deconvolution was attained with an algorithm that considered accumulated signal intensity due to pulse-on-pulse stacking, i.e., the latest pulse on top of all antecedent individual pulses' exponential tails proportionally.

Abstract: A combined experimental and first-principles study is performed to study the origin of conductivity in
ZnO:Al nanoparticles synthesized under controlled conditions via a reflux route using benzylamine as a
solvent. The experimental characterization of the samples by Raman, nuclear magnetic resonance (NMR)
and conductivity measurements indicates that upon annealing in nitrogen, the Al atoms at interstitial
positions migrate to the substitutional positions, creating at the same time Zn interstitials. We provide
evidence for the fact that the formed complex of AlZn and Zni corresponds to the origin of the Knight
shifted peak (KS) we observe in 27Al NMR. As far as we know, the role of this complex has not been
discussed in the literature to date. However, our first-principles calculations show that such a complex is
indeed energetically favoured over the isolated Al interstitial positions. In our calculations we also
address the charge state of the Al interstitials. Further, Zn interstitials can migrate from Al_Zn and possibly
also form Zn clusters, leading to the observed increased conductivity.

Abstract: The exact structure and properties of the Si vertical bar SiO2 interface are very important in microelectronics and photovoltaic devices such as metal-oxide-semiconductor field-effect transistors (MOSFETs) and solar cells. Whereas Si vertical bar SiO2 structures are traditionally produced by thermal oxidation, hyperthermal oxidation shows a number of promising advantages. However, the Si vertical bar SiO2 interface induced in hyperthermal Si oxidation has not been properly investigated yet. Therefore, in this work, the interface morphology and interfacial stresses during hyperthermal oxidation at room temperature are studied using reactive molecular dynamics simulations based on the ReaxFF potential. Interface thickness and roughness, as well as the bond length and bond angle distributions in the interface are discussed and compared with other models developed for the interfaces induced by traditional thermal oxidation. The formation of a compressive stress is observed. This compressive stress, which at the interface amounts about 2 GPa, significantly slows down the inward silica growth. This value is close to the experimental value in the Si vertical bar SiO2 interface obtained in traditional thermal oxidation.