Abstract: The global relaxation process after pulsed laser induced photodetachment in a rf electronegative SIH4 discharge is studied by a self-consistent kinetic one-dimensional particle-in-cell-Monte Carlo model. Our results reveal a comprehensive physical picture of the relaxation process, including the main plasma variables, after a perturbation up to the full recovery of the steady state. A strong influence of the photodetachment on the discharge is found, which results from an increase of the electron density, leading to a weaker bulk field, and hence to a drop in the high energy tail of the electron energy distribution function (EEDF), a reduction of the reaction rates of electron impact attachment and ionization, and a subsequent decrease of the positive and negative ion densities. All the plasma quantities related to electrons recover synchronously. The recovery time of the ion densities is about 1-2 orders of magnitude longer than that of the electrons due to different recovery mechanisms. The modeled behavior of all the charged particles agrees very well with experimental results from the literature. In addition, our work clarifies some unclear processes assumed in the literature, such as the relaxation of the EEDF, the evolution of the electric field, and the recovery of negative ions.

Abstract: The He-Ar-Cu+ IR laser operates in a hollow-cathode discharge, typically in a mixture of helium with a few-% Ar. The population inversion of the Cu+ ion levels, responsible for laser action, is attributed to asymmetric charge transfer between He+ ions and sputtered Cu atoms. The Ar gas is added to promote sputtering of the Cu cathode. In this paper, a hybrid modeling network consisting of several different models for the various plasma species present in a He-Ar-Cu hollow-cathode discharge is applied to investigate the effect of Ar concentration in the gas mixture on the discharge behavior, and to find the optimum He/Ar gas ratio for laser operation. It is found that the densities of electrons, Ar+ ions, Ar-m* metastable atoms, sputtered Cu atoms and Cu+ ions increase upon the addition of more Ar gas, whereas the densities of He+ ions, He-2(+) ions and He-m* metastable atoms drop considerably. The product of the calculated Cu atom and He+ ion densities, which determines the production rate of the upper laser levels, and hence probably also the laser output power, is found to reach a maximum around 1-5% Ar addition. This calculation result is compared to experimental measurements, and reasonable agreement has been reached.

Abstract: The influence of the composition on the crystallographic properties of deposited Mg(M)O (with M=Al, Cr, Ti, Y, and Zr) films is studied. For a flexible control of the composition, dual reactive magnetron sputtering was used as deposition technique. Two different approaches to predict the composition are discussed. The first is an experimental way based on the simple relationship between the deposition rate and the target-substrate distance. The second is a route using a Monte Carlo based particle trajectory code. Both methods require a minimal experimental input and enable the user to quickly predict the composition of complex thin films. Good control and flexibility allow us to study the compositional effects on the growth of Mg(M)O films. Pure MgO thin films were grown with a (111) preferential out-of-plane orientation. When adding M to MgO, two trends were noticed. The first trend is a change in the MgO lattice parameters compared to pure MgO. The second tendency is a decrease in the crystallinity of the MgO phase. The experimentally determined crystallographic properties are shown to be in correspondence with the predicted properties from molecular dynamics simulations.

Abstract: The interaction between cold atmospheric pressure plasma and liquids is receiving increasing attention for various applications. In particular, the use of plasma-treated liquids (PTL) for biomedical applications is of growing importance, in particular for sterilization and cancer treatment. However, insight into the
underlying mechanisms of plasma–liquid interactions is still scarce. Here, we present a 2D fluid dynamics model for the interaction between a plasma jet and liquid water. Our results indicate that the formed reactive species originate from either the gas phase (with further solvation) or are formed at the liquid interface. A clear increase in the aqueous density of H2O2, HNO2/NO2- and NO3-
is observed as a function of time, while the densities of O3, HO2/O2- and ONOOH/ONOO- are found to quickly reach a maximum due to chemical reactions in solution. The trends observed in our model correlate well with experimental observations from the literature.

Abstract: The interaction of thermal and hyperthermal Ni ions with gas-phase C60 fullerene was investigated at two temperatures with classical molecular dynamics simulations using a recently developed interatomic many-body potential. The interaction between Ni and C60 is characterized in terms of the NiC60 binding sites, complex formation, and the collision and temperature induced deformation of the C60 cage structure. The simulations show how ion implantation theoretically allows the synthesis of both endohedral Ni@C60 and exohedral NiC60 metallofullerene complexes.

Abstract: The mechanisms of plasma in medicine are broadly attributed to plasma-derived reactive oxygen and nitrogen species (RONS). In order to exert any intracellular effects, these plasma-derived RONS must first traverse a major barrier in the cell membrane. The cell membrane lipid composition, and thereby the magnitude of this barrier, is highly variable between cells depending on type and state (e.g. it is widely accepted that healthy and cancerous cells have different membrane lipid compositions). In this study, we investigate how plasma-derived RONS interactions with lipid membrane components can potentially be exploited in the future for treatment of diseases. We couple phospholipid vesicle experiments, used as simple cell models, with molecular dynamics (MD) simulations of the lipid membrane to provide new insights into how the interplay between phospholipids and cholesterol may influence the response of healthy and diseased cell membranes to plasma-derived RONS. We focus on the (i) lipid tail saturation degree, (ii) lipid head group type, and (iii) membrane cholesterol fraction. Using encapsulated molecular probes, we study the influence of the above membrane components on the ingress of RONS into the vesicles, and subsequent DNA damage. Our results indicate that all of the above membrane components can enhance or suppress RONS uptake, depending on their relative concentration within the membrane. Further, we show that higher RONS uptake into the vesicles does not always correlate with increased DNA damage, which is attributed to ROS reactivity and lifetime. The MD simulations indicate the multifactorial chemical and physical processes at play, including (i) lipid oxidation, (ii) lipid packing, and (iii) lipid rafts formation. The methods and findings presented here provide a platform of knowledge that could be leveraged in the development of therapies relying on the action of plasma, in which the cell membrane and oxidative stress response in cells is targeted.

Abstract: The need for carbon negative technologies led to the development of a wide array of novel CO₂conversion techniques. Most of them either rely on high temperatures or generate highly reactive O species, which can lead to the undesirable formation of NO_xand N₂O when the CO₂feeds contain N₂. Here, we show that, for plasma-based CO₂conversion, adding a hydrogen source, as a chemical oxygen scavenger, can suppress their formation,<italic>in situ</italic>. This allows the use of low-cost N₂containing (industrial and direct air capture) feeds, rather than expensive purified CO₂. To demonstrate this, we add CH₄to a dielectric barrier discharge plasma used for converting impure CO₂. We find that when adding a stoichiometric amount of CH₄, 82% less NO₂and 51% less NO are formed. An even higher reduction (96 and 63%) can be obtained when doubling this amount. However, in that case the excess radicals promote the formation of by-products, such as HCN, NH₃and CH₃OH. Thus, we believe that by using an appropriate amount of chemical scavengers, we can use impure CO₂feeds, which would bring us closer to ‘real world’ conditions and implementation.

Abstract: The peculiarities of sputtering processes at 0.5-5 keV Ne grazing ion bombardment of Si(001) and SiC(001) surfaces and their possible application for the surface modification have been studied by computer simulation. Sputtering yields in the primary knock-on recoil atoms regime versus the initial energy of incident ions (E(0) = 0.5-5 keV) and angle of incidence (psi = 0-30 degrees) counted from a target surface have been calculated. Comparative studies of layer-by-layer sputtering for Si(001) and SiC(001) surfaces versus the initial energy of incident ions as well as an effective sputtering and sputtering threshold are discussed.

Abstract: The phase transition processes induced by ultrashort, 100 fs pulsed laser irradiation of Au, Cu, and Ni are studied by means of a combined atomistic-continuum approach. A moderately low absorbed laser fluence range, from 200 to 600 J/m2 is considered to study phase transitions by means of a local and a nonlocal order parameter. At low laser fluences, the occurrence of layer-by-layer evaporation has been observed, which suggests a direct solid to vapor transition. The calculated amount of molten material remains very limited under the conditions studied, especially for Ni. Therefore, our results show that a kinetic equation that describes a direct solid to vapor transition might be the best approach to model laser-induced phase transitions by continuum models. Furthermore, the results provide more insight into the applicability of analytical superheating theories that were implemented in continuum models and help the understanding of nonequilibrium phase transitions.

Abstract: The physical processes in an Ar/O2 magnetron discharge used for the reactive sputter deposition of TiOx thin films were simulated with a 2d3v particle-in-cell/Monte Carlo collisions (PIC/MCC) model. The plasma species taken into account are electrons, Ar+ ions, fast Arf atoms, metastable Arm* atoms, Ti+ ions, Ti atoms, O+ ions, O2+ ions, O− ions and O atoms. This model accounts for plasmatarget interactions, such as secondary electron emission and target sputtering, and the effects of target poisoning. Furthermore, the deposition process is described by an analytical surface model. The influence of the O2/Ar gas ratio on the plasma potential and on the species densities and fluxes is investigated. Among others, it is shown that a higher O2 pressure causes the region of positive plasma potential and the O− density to be more spread, and the latter to decrease. On the other hand, the deposition rates of Ti and O are not much affected by the O2/Ar proportion. Indeed, the predicted stoichiometry of the deposited TiOx film approaches x=2 for nearly all the investigated O2/Ar proportions.

Abstract: The plasma behavior in a parallel-plate dielectric barrier discharge (DBD) is simulated by a two-dimensional particle-in-cell/Monte Carlo collision model, comparing for the first time an unpacked (empty) DBD with a packed bed DBD, i.e., a DBD filled with dielectric spheres in the gas gap. The calculations are performed in air, at atmospheric pressure. The discharge is powered by a pulse with a voltage amplitude of −20 kV. When comparing the packed and unpacked DBD reactors with the same dielectric barriers, it is clear that the presence of the dielectric packing leads to a transition in discharge behavior from a combination of negative streamers and unlimited surface streamers on the bottom dielectric surface to a combination of predominant positive streamers and limited surface discharges on the dielectric surfaces of the beads and plates. Furthermore, in the packed bed DBD, the electric field is locally enhanced inside the dielectric material, near the contact points between the beads and the plates, and therefore also in the plasma between the packing beads and between a bead and the dielectric wall, leading to values of $4\times {10}

Abstract: The present study analyses a specific laser ablation cell, the High Efficiency Aerosol Dispersion (HEAD) cell (see J. Pisonero et al., J. Anal. At. Spectrom., 2006, 21, 922931), by means of computational fluid dynamics (CFD) simulations. However, this cell consists of different modular parts, therefore, the results are probably of interest for the further development of other ablation cells. In the HEAD cell, the ablation spot is positioned below an orifice in the ceiling of the sample chamber. The particle transport through this orifice has been analysed for a ceiling height of 0.8 mm. The critical velocity for the onset of particle losses was found to be independent on the ejection angle at the crater spot. The deceleration of the particles can be described as the stopping in an effectively steady gas. Particle losses were negligible in this modular part of the cell at the evaluated laser ablation conditions. The transport efficiency through the Venturi chamber was investigated for different sample gas flow rates. In this case, it was found that small particles were predominantly lost at low flow rates, the large particles at higher flow rates. Making use of the simulation results, it was possible to design a modification of the HEAD cell that results in extremely short calculated washout times. The simulations yielded a signal of less than 10 ms, which was produced by more than 99% of the introduced mass.

Abstract: The purpose of this Brief Report is to point out the mistake in a formula for anisotropic electron scattering, previously published in Phys. Rev. A 41, 1112 (1990), which is widely used in Monte Carlo models of gas discharges. Anisotropic electron scattering is investigated based on the screened Coulomb potential between electrons and neutral atoms. The approach is also applied for electron scattering by nonpolar neutral molecules. Differential cross sections for electron scattering by Ar, N2, and CH4 are constructed on the basis of momentum and integrated cross sections. The formula derived in this paper is useful for Monte Carlo simulations of gas discharges.

Abstract: The reaction behavior of species that may affect the growth of ultrananocrystal line and nanocrystalline diamond ((U)NCD) films is investigated by means of molecular dynamics simulations. Impacts of CHx (x = 0 – 4), C2Hx (x=0-6), C3Hx (x=0-2), C4Hx (x = 0 – 2), H, and H-2 on clean and hydrogenated diamond (100)2 x 1 and (111) 1 x 1 surfaces at two different substrate temperatures are simulated. We find that the different bonding structures of the two surfaces cause different temperature effects on the sticking efficiency. These results predict a temperature-dependent ratio of diamond (100) and (111) growth. Furthermore, predictions of which are the most important hydrocarbon species for (U)NCD growth are made.

Abstract: The requirements for new technologies to serve as anticancer agents go far beyond their toxicity potential. Novel applications also need to be safe on a molecular and patient level. In a broader sense, this also relates to cancer metastasis and inflammation. In a previous study, the toxicity of an atmospheric pressure argon plasma jet in four human pancreatic cancer cell lines was confirmed and plasma treatment did not promote metastasis in vitro and in ovo. Here, these results are extended by additional types of analysis and new models to validate and define on a molecular level the changes related to metastatic processes in pancreatic cancer cells following plasma treatment in vitro and in ovo. In solid tumors that were grown on the chorion-allantois membrane of fertilized chicken eggs (TUM-CAM), plasma treatment induced modest to profound apoptosis in the tissues. This, however, was not associated with a change in the expression levels of adhesion molecules, as shown using immunofluorescence of ultrathin tissue sections. Culturing of the cells detached from these solid tumors for 6d revealed a similar or smaller total growth area and expression of ZEB1, a transcription factor associated with cancer metastasis, in the plasma-treated pancreatic cancer tissues. Analysis of in vitro and in ovo supernatants of 13 different cytokines and chemokines revealed cell line-specific effects of the plasma treatment but a noticeable increase of, e.g., growth-promoting interleukin 10 was not observed. Moreover, markers of epithelial-to-mesenchymal transition (EMT), a metastasis-promoting cellular program, were investigated. Plasma-treated pancreatic cancer cells did not present an EMT-profile. Finally, a realistic 3D tumor spheroid co-culture model with pancreatic stellate cells was employed, and the invasive properties in a gel-like cellular matrix were investigated. Tumor outgrowth and spread was similar or decreased in the plasma conditions. Altogether, these results provide valuable insights into the effect of plasma treatment on metastasis-related properties of cancer cells and did not suggest EMT-promoting effects of this novel cancer therapy.