Abstract: Aquaporins (AQPs) are transmembrane proteins that conduct not only water molecules across the cell membrane but also other solutes, such as reactive oxygen and nitrogen species (RONS), produced (among others) by cold atmospheric plasma (CAP). These RONS may induce oxidative stress in the cell interior, which plays a role in cancer treatment. The underlying mechanisms of the transport of RONS across AQPs, however, still remain obscure. We apply molecular dynamics simulations to investigate the permeation of both hydrophilic (H₂O₂and OH) and hydrophobic (NO₂and NO) RONS through AQP1. Our simulations show that these RONS can all penetrate across the pores of AQP1. The permeation free energy barrier of OH and NO is lower than that of H₂O₂and NO₂, indicating that these radicals may have easier access to the pore interior and interact with the amino acid residues of AQP1. We also study the effect of RONS-induced oxidation of both the phospholipids and AQP1 (i.e., sulfenylation of Cys₁₉₁) on the transport of the above-mentioned RONS across AQP1. Both lipid and protein oxidation seem to slightly increase the free energy barrier for H₂O₂and NO₂permeation, while for OH and NO, we do not observe a strong effect of oxidation. The simulation results help to gain insight in the underlying mechanisms of the noticeable rise of CAP-induced RONS in cancer cells, thereby improving our understanding on the role of AQPs in the selective anticancer capacity of CAP.

Abstract: Plasma catalysis is promising for various environmental, energy and chemical synthesis applications, but the underlying mechanisms are far from understood. Modeling can help to obtain a better insight in these mechanisms. Some burning questions relate to the plasma behavior inside packed bed reactors and whether plasma can penetrate into catalyst pores. In this paper, we try to provide answers to these questions, by means of both fluid modeling and particle-in-cell/Monte Carlo collision simulations. We present a short overview of recent findings obtained in our group by means of modeling, i.e., the enhanced electric field near the contact points and the streamer propagation through the packing in packed bed reactors, as well as the plasma behavior in catalyst pores, to determine the minimum pore size in which plasma streamers can penetrate.

Abstract: Magnetron sputtering is a widely used physical vapor deposition technique for deposition and formation of nanocatalyst thin films and clusters. Nevertheless, so far only few studies investigated this formation process at the fundamental level. We here review atomic scale molecular dynamics simulations aimed at elucidating the nanocatalyst growth process through magnetron sputtering. We first introduce the basic magnetron sputtering background and machinery of molecular dynamics simulations, and then describe the studies conducted in this field so far. We also present a perspective view on how the field may be developed further.

Abstract: Utilizing ammonia as a storage medium for hydrogen is currently receiving increased attention. A possible method to retrieve the hydrogen is by plasma-catalytic decomposition. In this work, we combined an experimental study, using a dielectric barrier discharge plasma reactor, with a plasma kinetic model, to get insights into the decomposition mechanism. The experimental results revealed a similar effect on the ammonia conversion when changing the flow rate and power, where increasing the specific energy input (higher power or lower flow rate) gave an increased conversion. A conversion as high as 82 % was achieved at a specific energy input of 18 kJ/Nl. Furthermore, when changing the discharge volume from 31 to 10 cm3, a change in the plasma distribution factor from 0.2 to 0.1 was needed in the model to best describe the conversions of the experiments. This means that a smaller plasma volume caused a higher transfer of energy through micro-discharges (non-uniform plasma), which was found to promote the decomposition of ammonia. These results indicate that it is the collisions between NH3 and the high-energy electrons that initiate the decomposition. Moreover, the rate of ammonia destruction was found by the model to be in the order of 1022 molecules/(cm3 s) during the micro-discharges, which is 5 to 6 orders of magnitude higher than in the afterglows. A considerable re-formation of ammonia was found to take place in the afterglows, limiting the overall conversion. In addition, the model revealed that implementation of packing material in the plasma introduced high concentrations of surface-bound hydrogen atoms, which introduced an additional ammonia re-formation pathway through an Eley-Rideal reaction with gas phase NH2. Furthermore, a more uniform plasma is predicted in the presence of MgAl2O4, which leads to a lower average electron energy during micro-discharges and a lower conversion (37 %) at a comparable residence time for the plasma alone (51 %).

Abstract: Tunnel field-effect transistors based on van der Waals heterostructures are emerging device concepts for low-power applications, auguring sub-60 mV/dec subthreshold swing values. In these devices, the channel is built from a stack of several different two-dimensional materials whose nature allows tailoring the band alignments and enables a good electrostatic control of the device. In this work, we propose a theoretical study of the variability of the performances of a MoS2-ZrS2 tunnel field-effect transistor induced by fluctuations of the relative position or the orientation of the layers. Our results indicate that although a steep subthreshold slope (20 mV/dec) is achievable, fluctuations in the relative orientation of the ZrS2 layer with respect to the MoS2 one lead to a significant variability in the tunneling current by about one decade. This arises from changes in the orbital overlap between the layers and from the modulation of the transport direction.

Abstract: A database-driven approach combined with ab initio density functional theory (DFT) simulations is used to identify and simulate alternative ferroelectric materials beyond Hf(Zr)O-2. The database-driven screening method identifies a class of wurtzite ferroelectric materials. DFT simulations of wurtzite magnesium chalcogenides, including MgS, MgSe, and MgTe, show their potential to achieve improved ferroelectric (FE) stability, simple atomistic unit cell structure, and large FE polarization. Strain engineering can effectively modulate the FE switching barrier height for facilitating FE switching. The effect of the piezoelectric property on the FE switching barrier heights is also examined.

Abstract: We introduce a simple enhanced sampling approach for the calculation of free energy differences and barriers along a one-dimensional reaction coordinate. First, a small number of short nonequilibrium simulations are carried out along the reaction coordinate, and the Jarzynski equality is used to learn an approximate free energy surface from the nonequilibrium work distribution. This free energy estimate is represented in a compact form as an artificial neural network and used as an external bias potential to accelerate rare events in a subsequent molecular dynamics simulation. The final free energy estimate is then obtained by reweighting the equilibrium probability distribution of the reaction coordinate sampled under the influence of the external bias. We apply our reweighted Jarzynski sampling recipe to four processes of varying scales and complexities.spanning chemical reaction in the gas phase, pair association in solution, and droplet nucleation in supersaturated vapor. In all cases, we find reweighted Jarzynski sampling to be a very efficient strategy, resulting in rapid convergence of the free energy to high precision.

Abstract: The hyperdynamics method is a powerful tool to simulate slow processes at the atomic level. However, the construction of an optimal hyperdynamics potential is a task that is far from trivial. Here, we propose a generally applicable implementation of the hyperdynamics algorithm, borrowing two concepts from metadynamics. First, the use of a collective variable (CV) to represent the accelerated dynamics gives the method a very large flexibility and simplicity. Second, a metadynamics procedure can be used to construct a suitable history-dependent bias potential on-the-fly, effectively turning the algorithm into a self-learning accelerated molecular dynamics method. This collective variable-driven hyperdynamics (CVHD) method has a modular design: both the local system properties on which the bias is based, as well as the characteristics of the biasing method itself, can be chosen to match the needs of the considered system. As a result, system-specific details are abstracted from the biasing algorithm itself, making it extremely versatile and transparent. The method is tested on three model systems: diffusion on the Cu(001) surface and nickel-catalyzed methane decomposition, as examples of reactive processes with a bond-length-based CV, and the folding of a long polymer-like chain, using a set of dihedral angles as a CV. Boost factors up to 109, corresponding to a time scale of seconds, could be obtained while still accurately reproducing correct dynamics.

Abstract: The deposition of ultra-thin tungsten films and their related 2D chalcogen compounds on large area dielectric substrates by gas phase reactions is challenging. The lack of nucleation sites complicates the adsorption of W-related precursors and subsequent sulfurization usually requires high temperatures. We propose here a technique in which a thin solid amorphous silicon film is used as reductant for the gas phase precursor WF6 leading to the conversion to metallic W. The selectivity of the W conversion towards the underlying dielectric surfaces is demonstrated. The role of the Si surface preparation, the conversion temperature, and Si thickness on the formation process is investigated. Further, the in situ conversion of the metallic tungsten into thin stoichiometric WS2 is achieved by a cyclic approach based on WF6 and H2S pulses at the moderate temperature of 450 1C, which is much lower than usual oxide sulfurization processes.

Abstract: Catalysis science and technology increased efforts recently to progress beyond conventional “thermal” catalysis and face the challenges of net-zero emissions and electrification of production. Nevertheless, a better gaps and opportunities analysis is necessary. This review analyses four emerging areas of unconventional or less- conventional catalysis which share the common aspect of using directly renewable energy sources: (i) plasma catalysis, (ii) catalysis for flow chemistry and process intensification, (iii) application of electromagnetic (EM) fields to modulate catalytic activity and (iv) nanoscale generation at the catalyst interface of a strong local EM by plasmonic effect. Plasma catalysis has demonstrated synergistic effects, where the outcome is higher than the sum of both processes alone. Still, the underlying mechanisms are complex, and synergy is not always obtained. There is a crucial need for a better understanding to (i) design catalysts tailored to the plasma environment, (ii) design plasma reactors with optimal transport of plasma species to the catalyst surface, and (iii) tune the plasma conditions so they work in optimal synergy with the catalyst. Microfluidic reactors (flow chemistry) is another emerging sector leading to the intensification of catalytic syntheses, particularly in organic chemistry. New unconventional catalysts must be designed to exploit in full the novel possibilities. With a focus on (a) continuous-flow photocatalysis, (b) electrochemical flow catalysis, (c) microwave flow catalysis and (d) ultra­ sound flow activation, a series of examples are discussed, with also indications on scale-up and process indus­ trialisation. The third area discussed regards the effect on catalytic performances of applying oriented EM fields spanning several orders of magnitude. Under well-defined conditions, gas breakdown and, in some cases, plasma formation generates activated gas phase species. The EM field-driven chemical conversion processes depend further on structured electric/magnetic catalysts, which shape the EM field in strength and direction. Different effects influencing chemical conversion have been reported, including reduced activation energy, surface charging, hot spot generation, and selective local heating. The last topic discussed is complementary to the third, focusing on the possibility of tuning the photo- and electro-catalytic properties by creating a strong localised electrical field with a plasmonic effect. The novel possibilities of hot carriers generated by the plasmonic effect are also discussed. This review thus aims to stimulate the reader to make new, creative catalysis to address the challenges of reaching a carbon-neutral world.

Abstract: In this work, the effects of N-doping into the Co-doped single vacancy (Co-SV-G) and di-vacancy graphene flake (Co-dV-G) are investigated and compared toward direct oxidation of methane to methanol (DOMM) employing two different oxidants (N2O and O-2) using density functional theory (DFT) calculation. We found that DOMM on CoN3-G utilizing the N2O molecule as oxygen-donor proceeds via a two-step reaction with low activation energies. In addition, we found that although CoN3-G might be a good catalyst for methane conversion, it can also catalyze the oxidation of methanol to CO2 and H2O due to the required low activation barriers. Moreover, the adsorption behaviors of CHx (x = 0-4) species and dehydrogenation of CHx (x = 1-4) species on CoN3-G are investigated. We concluded that CoN3-G can be used as an efficient catalyst for DOMM and N-2O reduction at ambient conditions which may serve as a guide for fabricating effective C/N catalysts in energy-related devices.

Abstract: Pancreatic ductal adenocarcinoma (PDAC) is a challenging neoplastic disease, mainly due to the development of resistance to radio- and chemotherapy. Cold atmospheric plasma (CAP) is an alternative technology that can eliminate cancer cells through oxidative damage, as shown in vitro, in ovo, and in vivo. However, how CAP affects the pancreatic stellate cells (PSCs), key players in the invasion and metastasis of PDAC, is poorly understood. This study aims to determine the effect of an anti-PDAC CAP treatment on PSCs tissue developed in ovo using mono- and co-cultures of RLT-PSC (PSCs) and Mia PaCa-2 cells (PDAC). We measured tissue reduction upon CAP treatment and mRNA expression of PSC activation markers and extracellular matrix (ECM) remodelling factors via qRT-PCR. Protein expression of selected markers was confirmed via immunohistochemistry. CAP inhibited growth in Mia PaCa-2 and co-cultured tissue, but its effectiveness was reduced in the latter, which correlates with reduced ki67 levels. CAP did not alter the mRNA expression of PSC activation and ECM remodelling markers. No changes in MMP2 and MMP9 expression were observed in RLT-PSCs, but small changes were observed in Mia PaCa-2 cells. Our findings support the ability of CAP to eliminate PDAC cells, without altering the PSCs.

Abstract: Despite modern therapeutic advances, the survival prospects of pancreatic cancer patients remain poor, due to chemoresistance and dysregulated oncogenic kinase signaling networks. We applied a novel kinome activitymapping approach using biological peptide targets as phospho-sensors to identify vulnerable kinase de­pendencies for therapy sensitization by physical plasma. Ser/Thr-kinome specific activity changes were mapped upon induction of ferroptotic cell death in pancreatic tumor cells exposed to reactive oxygen and nitrogen species of plasma-treated water (PTW). This revealed a broad kinome activity response involving the CAMK, the AGC and CMGC family of kinases. This systems-level kinome network response supports stress adaptive switches between chemoresistant anti-oxidant responses of Kelch-like ECH-associated protein 1 (KEAP1)/Heme Oxy­genase 1 (HMOX1) and ferroptotic cell death sensitization upon suppression of Nuclear factor (erythroid derived 2)-like 2 (NRF2) and Glutathione peroxidase 4 (GPX4). This is further supported by ex vivo experiments in the chicken chorioallantoic membrane assay, showing decreased GPX4 and Glutathione (GSH) expression as well as increased lipid peroxidation, along with suppressed BxPC-3 tumor growth in response to PTW. Taken all together, we demonstrate that plasma treated water-derived oxidants sensitize pancreatic cancer cells to fer­roptotic cell death by targeting a NRF2-HMOX1-GPX4 specific kinase signaling network.