Abstract: Plasma catalysis is gaining increasing interest for various environmental applications, but the crucial question is whether plasma can be created inside catalyst pores and under which conditions. In practice, various catalytic support materials are used, with various dielectric constants. We investigate here the influence of the dielectric constant on the plasma properties inside catalyst pores and in the sheath in front of the pores, for various pore sizes. The calculations are performed by a two-dimensional fluid model for an atmospheric pressure dielectric barrier discharge in helium. The electron impact ionization rate, electron temperature, electron and ion density, as well as the potential distribution and surface charge density, are analyzed for a better understanding of the discharge behavior inside catalyst pores. The results indicate that, in a 100 μm pore, the electron impact ionization in the pore, which is characteristic for the plasma generation inside the pore, is greatly enhanced for dielectric constants below 300. Smaller pore sizes only yield enhanced ionization for smaller dielectric constants, i.e., up to εr = 200, 150, and 50 for pore sizes of 50, 30, and 10 μm. Thus, the most common catalyst supports, i.e., Al2O3 and SiO2, which have dielectric constants around εr = 8−11 and 4.2, respectively, should allow more easily that microdischarges can be formed inside catalyst pores, even for smaller pore sizes. On the other hand, ferroelectric materials with dielectric constants above 300 never seem to yield plasma enhancement inside catalyst pores, not even for 100 μm pore sizes. Furthermore, it is clear that the dielectric constant of the material has a large effect on the extent of plasma enhancement inside the catalyst pores, especially in the range between εr = 4 and εr = 200. The obtained results are explained in detail based on the surface charge density at the pore walls,
and the potential distribution and electron temperature inside and above the pores. The results obtained with this model are
important for plasma catalysis, as the production plasma species in catalyst pores might affect the catalyst properties, and thus
improve the applications of plasma catalysis.

Abstract: The application of nonthermal atmospheric pressure plasma is emerging as an alternative and efficient technique for the inactivation of bacterial biofilms. In this study, reactive molecular dynamics simulations were used to examine the reaction mechanisms of hydroxyl radicals, as key reactive oxygen plasma species in biological systems, with several organic molecules (i.e., alkane, alcohol, carboxylic acid, and amine), as prototypical components of biomolecules in the biofilm. Our results demonstrate that organic molecules containing hydroxyl and carboxyl groups may act as trapping agents for the OH radicals. Moreover, the impact of OH radicals on N-acetyl-glucosamine, as constituent component of staphylococcus epidermidis biofilms, was investigated. The results show how impacts of OH radicals lead to hydrogen abstraction and subsequent molecular damage. This study thus provides new data on the reaction mechanisms of plasma species, and particularly the OH radicals, with fundamental components of bacterial biofilms.

Abstract: Nanostructured carbon materials show a tremendous variety in atomic structure, morphology, properties, and applications. As all properties are ultimately determined by the structure of the material, a thorough understanding of the growth mechanisms that give rise to the particular structure is critical. On many occasions, it has been shown that plasma enhanced growth can be strongly beneficial. This review will describe the authors current understanding of plasma enhanced growth of carbon nanotubes, the prototypical example of nanostructured carbon materials, as obtained from experiments, simulations, and modeling. Specific emphasis is put on where experiments and computational approaches correspond, and where they differ. Also, the current status on simulating PECVD growth of some other carbon nanomaterials is reviewed, including amorphous carbon, graphene, and metallofullerenes. Finally, computational challenges with respect to the simulation of PECVD growth are identified.

Abstract: Molecular dynamics simulations are used to study the effects of different cations on the permeation of charged polymers through flat capillaries with heights below 2 nm. Interestingly, we found that, despite being monovalent, Li+ , Na+ , and K+ cations have different effects on polymer permeation, which consequently affects their transmission speed throughout those capillaries. We attribute this phenomenon to the interplay of the cations' hydration free energies and the hydrodynamic drag in front of the polymer when it enters the capillary. Different alkali cations exhibit different surface versus bulk preferences in small clusters of water under the influence of an external electric field. This paper presents a tool to control the speed of charged polymers in confined spaces using cations.

Abstract: Heterostructures of low-dimensional semiconducting materials, such as transition metal dichalcogenides (MX2), are promising building blocks for future electronic and optoelectronic devices. The patterning of one MX2 material on top of another one is challenging due to their structural similarity. This prevents an intrinsic etch stop when conventional anisotropic dry etching processes are used. An alternative approach consist in a two-step process, where a sacrificial silicon layer is pre-patterned with a low damage plasma process, stopping on the underlying MoS2 film. The pre-patterned layer is used as sacrificial template for the formation of the top WS2 film. This study describes the optimization of a cyclic Ar/Cl-2 atomic layer etch process applied to etch silicon on top of MoS2, with minimal damage, followed by a selective conversion of the patterned Si into WS2. The impact of the Si atomic layer etch towards the MoS2 is evaluated: in the ion energy range used for this study, MoS2 removal occurs in the over-etch step over 1-2 layers, leading to the appearance of MoOx but without significant lattice distortions to the remaining layers. The combination of Si atomic layer etch, on top of MoS2, and subsequent Si-to-WS2 selective conversion, allows to create a WS2/MoS2 heterostructure, with clear Raman signals and horizontal lattice alignment. These results demonstrate a scalable, transfer free method to achieve horizontally individually patterned heterostacks and open the route towards wafer-level processing of 2D materials.

Abstract: In recent years, the ability of cold atmospheric pressure plasmas (CAPS) to selectively induce cell death in cancer cells has been widely established. This selectivity has been assigned to the reactive oxygen and nitrogen species (RONS) created in CAPs. To provide new insights in the search for an explanation
for the observed selectivity, we calculate the transfer free energy of multiple ROS across membranes containing a varying amount of cholesterol. The cholesterol fraction is investigated as a selectivity parameter because membranes of cancer cells are known to contain lower fractions of cholesterol compared to healthy cells. We find that cholesterol has a significant effect on the permeation of
reactive species across a membrane. Indeed, depending on the specific reactive species, an increasing cholesterol fraction can lead to (i) an increase of the transfer free energy barrier height and width, (ii) the formation of a local free energy minimum in the center of the membrane and (iii) the creation of extra free energy barriers due to the bulky sterol rings. In the context of plasma oncology, these observations suggest that the increased ingress of RONS in cancer cells can be explained by the decreased cholesterol fraction of their cell membrane.

Abstract: Modelling at the molecular or atomic scale can be very useful for obtaining a better insight in plasma medicine. This paper gives an overview of different atomic/molecular scale modelling approaches that can be used to study the direct interaction of plasma species with biomolecules or the consequences of these interactions for the biomolecules on a somewhat longer time-scale. These approaches include density functional theory (DFT), density functional based tight binding (DFTB), classical reactive and non-reactive molecular dynamics (MD) and united-atom or coarse-grained MD, as well as hybrid quantum mechanics/molecular mechanics (QM/MM) methods. Specific examples will be given for three important types of biomolecules, present in human cells, i.e. proteins, DNA and phospholipids found in the cell membrane. The results show that each of these modelling approaches has its specific strengths and limitations, and is particularly useful for certain applications. A multi-level approach is therefore most suitable for obtaining a global picture of the plasma–biomolecule interactions.

Abstract: Plasma generation inside catalyst pores is of utmost importance for plasma catalysis, as the existence of plasma species inside the pores affects the active surface area of the catalyst available to the plasma species for catalytic reactions. In this paper, the electric field enhancement, and thus the plasma production inside catalyst pores with different pore shapes is studied with a two-dimensional fluid model. The results indicate that the electric field will be significantly enhanced near tip-like structures. In a conical pore with small opening, the strongest electric field appears at the opening and bottom corners of the pore, giving rise to a prominent ionization rate throughout the pore. For a cylindrical pore, the electric field is only enhanced at the bottom corners of the pore, with lower absolute value, and thus the ionization rate inside the pore is only slightly enhanced. Finally, in a conical pore with large opening, the electric field is characterized by a maximum at the bottom of the pore, yielding a similar behavior for the ionization rate. These results demonstrate that the shape of the pore has a significantly influence on the electric field enhancement, and thus modifies the plasma properties.

Abstract: In recent years there has been growing interest in the use of low-temperature atmospheric pressure plasmas for biomedical applications. Currently, however, there is very little fundamental knowledge regarding the relevant interaction mechanisms of plasma species with living cells. In this paper, we investigate the interaction of important plasma species, such as O3, O2 and O atoms, with bacterial peptidoglycan (or murein) by means of reactive molecular dynamics simulations. Specifically, we use the peptidoglycan structure to model the gram-positive bacterium Staphylococcus aureus murein. Peptidoglycan is the outer protective barrier in bacteria and can therefore interact directly with plasma species. Our results demonstrate that among the species mentioned above, O3 molecules and especially O atoms can break important bonds of the peptidoglycan structure (i.e. CO, CN and CC bonds), which subsequently leads to the destruction of the bacterial cell wall. This study is important for gaining a fundamental insight into the chemical damaging mechanisms of the bacterial peptidoglycan structure on the atomic scale.

Abstract: Piezoelectricity is considered to be one of the key functionalities in biomaterials to boost bone tissue regeneration, however, integrating biocompatibility, biodegradability and 3D structure with pronounced piezoresponse remains a material challenge. Herein, novel hybrid biocompatible 3D scaffolds based on biodegradable poly(3-hydroxybutyrate) (PHB) and reduced graphene oxide (rGO) flakes have been developed. Nanoscale insights revealed a more homogenous distribution and superior surface potential values of PHB fibers (33 +/- 29 mV) with increasing rGO content up to 1.0 wt% (314 +/- 31 mV). The maximum effective piezoresponse was detected at 0.7 wt% rGO content, demonstrating 2.5 and 1.7 times higher out-of-plane and in-plane values, respectively, than that for pure PHB fibers. The rGO addition led to enhanced zigzag chain formation between paired lamellae in PHB fibers. In contrast, a further increase in rGO content reduced the alpha-crystal size and prevented zigzag chain conformation. A corresponding model explaining structural and molecular changes caused by rGO addition in electrospun PHB fibers is proposed. In addition, finite element analysis revealed a negligible vertical piezoresponse compared to lateral piezoresponse in uniaxially oriented PHB fibers based on alpha-phase (P2(1)2(1)2(1) space group). Thus, the present study demonstrates promising results for the development of biodegradable hybrid 3D scaffolds with an enhanced piezoresponse for various tissue engineering applications.

Abstract: Molecular dynamics simulations using improved TersoffBrenner potential parameters were performed to investigate Cl+ etching of a {2×1} reconstructed Si(100) surface. Steady-state Si etching accompanying the Cl coverage of the surface is observed. Furthermore, a steady-state chlorinated reaction layer is formed. The thickness of this reaction layer is found to increase with increasing energy. The stoichiometry of SiClx species in the reaction layer is found to be SiCl:SiCl2:SiCl3 = 1.0:0.14:0.008 at 50 eV. These results are in excellent agreement with available experimental data. While elemental Si products are created by physical sputtering, most SiClx (0<x<4) etch products are produced by chemical-enhanced physical sputtering.

Abstract: Natural deep eutectic solvents (NADES) represent a green alternative to common organic solvents in the biochemical industry due to their benign behavior and tailorable properties, in particular as media for enzymatic reactions. However, to fully exploit their potential in enzymatic reactions, there is a need for a more fundamental understanding of how these neoteric solvents influence the course of these reac-tions. Thus, the aim of this study is to investigate the influence of NADES with various molar composi-tions on the stability and structure of enzymes, applying molecular dynamics simulations. This can help to better understand the effect of individual compounds of NADES, in addition to eutectic mixtures. More specifically, we simulate the behavior of Candida antarctica lipase B (CALB) enzyme in NADES com-posed of choline chloride with either urea, ethylene glycol or glycerol. Hereto, we monitor the NADES microstructure, the general stability of the enzyme and changes in the structure of its active sites and sur-face residues. Our simulations show that none of the studied NADES systems significantly disrupt the microstructure of the solvent or the stability of the CALB enzyme within the time scales of the simula-tions. The enzyme preserves its initial structure, size and intra-chain hydrogen bonds in all investigated compositions and, for the first time reported, also in NADES with increased hydrogen bond donating com-pound ratios. As the main novelty, our results indicate that, in addition to the composition, the molar ratio can be an additional variable to fine-tune the physicochemical properties of NADES without altering the enzyme characteristics. These findings could facilitate the development and application of task -tailored NADES media for biocatalytic processes. (c) 2022 Elsevier B.V. All rights reserved.

Abstract: Monte Carlo (MC) methods have a long-standing history as partners of molecular dynamics (MD) to simulate the evolution of materials at the atomic scale. Among these techniques, the uniform-acceptance force-bias Monte Carlo (UFMC) method [ G. Dereli Mol. Simul. 8 351 (1992)] has recently attracted attention [ M. Timonova et al. Phys. Rev. B 81 144107 (2010)] thanks to its apparent capacity of being able to simulate physical processes in a reduced number of iterations compared to classical MD methods. The origin of this efficiency remains, however, unclear. In this work we derive a UFMC method starting from basic thermodynamic principles, which leads to an intuitive and unambiguous formalism. The approach includes a statistically relevant time step per Monte Carlo iteration, showing a significant speed-up compared to MD simulations. This time-stamped force-bias Monte Carlo (tfMC) formalism is tested on both simple one-dimensional and three-dimensional systems. Both test-cases give excellent results in agreement with analytical solutions and literature reports. The inclusion of a time scale, the simplicity of the method, and the enhancement of the time step compared to classical MD methods make this method very appealing for studying the dynamics of many-particle systems.

Abstract: The ever increasing global production and dispersion of methane requires novel chemistry to transform it into easily condensable energy carriers that can be integrated into the chemical infrastructure. In this context, single atom catalysts have attracted considerable interest due to their outstanding catalytic activity. We here use density functional theory (DFT) computations to compare the reaction and activation energies of M and MN4 embedded graphene (M = Ni and Si) on the methane-to-methanol conversion near room temperature. Thermodynamically, conversion of methane to methanol is energetically favorable at ambient conditions. Both singlet and triplet spin state of the studied systems are considered in all of the calculations. The DFT results show that the barriers are significantly lower when the complexes are in the triplet state than in the singlet state. In particular, Si-G with the preferred spin multiplicity of triplet seems to be viable catalysts for methane oxidation thanks to the corresponding lower energy barriers and higher stability of the obtained configurations. Our results provide insights into the nature of methane conversion and may serve as guidance for fabricating cost-effective graphene-based single atom catalysts.