Abstract: In this tutorial review paper, we illustrate how computer modeling can contribute to a better insight in inductively coupled plasma-mass spectrometry (ICP-MS). We start with a brief overview on previous efforts, studying the fundamentals of the ICP and ICP-MS, with main focus on previous modeling activities. Subsequently, we explain in detail the model that we developed in previous years, and we show typical calculation results, illustrating the plasma characteristics, gas flow patterns and the sample transport, evaporation and ionization. We also present the effect of various experimental parameters, such as operating conditions, geometrical aspects and sample characteristics, to illustrate how modeling can help to elucidate the optimal conditions for improved analytical performance.

Abstract: In this work we evaluate the chemical kinetics of dry reforming of methane in warm plasmas (1000–4000 K) using modelling with a newly developed chemistry set, for a broad range of parameters (temperature, power density and CO2/CH4 ratio). We compare the model against thermodynamic equilibrium concentrations, serving as validation of the thermal chemical kinetics. Our model reveals that plasma-specific reactions (i.e., electron impact collisions) accelerate the kinetics compared to thermal conversion, rather than altering the overall kinetics pathways and intermediate products, for gas temperatures below 2000 K. For higher temperatures, the kinetics are dominated by heavy species collisions and are strictly thermal, with negligible influence of the electrons and ions on the overall kinetics. When studying the effects of different gas mixtures on the kinetics, we identify important intermediate species, side reactions and side products. The use of excess CO2 leads to H2O formation, at the expense of H2 formation, and the CO2 conversion itself is limited, only approaching full conversion near 4000 K. In contrast, full conversion of both reactants is only kinetically limited for mixtures with excess CH4, which also gives rise to the formation of C2H2, alongside syngas. Within the given parameter space, our model predicts the 30/70 ratio of CO2/CH4 to be the most optimal for syngas formation with a H2/CO ratio of 2.

Abstract: In this work we study energy transport in a gliding arc discharge with two diverging flat
electrodes in argon gas at atmospheric pressure. The discharge is ignited at the shortest electrode
gap and it is pushed downstream by a forced gas flow. The current values considered are
relatively low and therefore a non-equilibrium plasma is produced. We consider two cases, i.e.
with high and low discharge current (28 mA and 2.8mA), and a constant gas flow of 10 lmin −1 ,
with a significant turbulent component to the velocity. The study presents an analysis of the
various energy transport mechanisms responsible for the redistribution of Joule heating to the
plasma species and the moving background gas. The objective of this work is to provide a
general understanding of the role of the different energy transport mechanisms in arc formation
and sustainment, which can be used to improve existing or new discharge designs. The work is
based on a three-dimensional numerical model, combining a fluid plasma model, the shear stress
transport Reynolds averaged Navier–Stokes turbulent gas flow model, and a model for gas
thermal balance. The obtained results show that at higher current the discharge is constricted
within a thin plasma column several hundred kelvin above room temperature, while in the low-
current discharge the combination of intense convective cooling and low Joule heating prevents
discharge contraction and the plasma column evolves to a static non-moving diffusive plasma,
continuously cooled by the flowing gas. As a result, the energy transport in the two cases is
determined by different mechanisms. At higher current and a constricted plasma column, the
plasma column is cooled mainly by turbulent transport, while at low current and an unconstricted
plasma, the major cooling mechanism is energy transport due to non-turbulent gas convection. In
general, the study also demonstrates the importance of turbulent energy transport in
redistributing the Joule heating in the arc and its significant role in arc cooling and the formation
of the gas temperature profile. In general, the turbulent energy transport lowers the average gas
temperature in the arc, thus allowing additional control of thermal non-equilibrium in the
discharge.

Abstract: In this work, a hybrid model is used to investigate the effect of different gas ratios on the Si etching and polymer film deposition characteristics in an Ar/CF4 inductively coupled plasma. The influence of the surface processes on the bulk plasma properties is studied, and also the spatial characteristics of important gas phase and etched species. The densities of F and CF2 decrease when the surface module is included in the simulations, due to the species consumption caused by etching and polymer deposition. The influence of the surface processes on the bulk plasma depends on the Ar/CF4 gas ratio. The deposited polymer becomes thicker at high CF4 content because of more abundant CFx radicals. As a result of the competition between the polymer thickness and the F flux, the etch rate first increases and then decreases upon increasing the CF4 content. The electron properties, more specifically the electron density profile, affect the Si etch characteristics substantially by determining the radical density and flux profiles. In fact, the radial profile of the etch rate is more uniform at low CF4 content since the electron density has a smooth distribution. At high CF4 content, the etch rate is less uniform with a minimum halfway along the wafer radius, because the electron density distribution is more localized. Therefore, our calculations predict that it is better to work at relatively high Ar/CF4 gas ratios, in order to obtain high etch rate and good profile uniformity for etch applications. This, in fact, corresponds to the typical experimental etch conditions in Ar/CF4 gas mixtures as found in the literature, where Ar is typically present at a much higher concentration than CF4.

Abstract: In this work, the role of vibrationally excited HBr (HBr(vib)) is computationally investigated for a HBr/He inductively coupled plasma applied for Si etching. It is found that at least 50% of all dissociations of HBr occur through HBr(vib). This additional dissociation pathway through HBr(vib) makes the plasma significantly more atomic. It also results in a slightly higher electron temperature (i.e. about 0.2 eV higher compared to simulation results where HBr(vib) is not included), as well as a higher gas temperature (i.e. about 50 K higher than without including HBr(vib)), due to the enhanced Franck–Condon heating through HBr(vib) dissociation,
at the conditions investigated. Most importantly, the calculated etch rate with HBr(vib) included in the model is a factor 3 higher than in the case without HBr(vib), due to the higher fluxes of etching species (i.e. H and Br), while the chemical composition of the wafer surface shows no significant difference. Our calculations clearly show the importance of including HBr(vib) for accurate modeling of HBr-containing plasmas.

Abstract: In this work, we study the abatement of dilute trichloroethylene (TCE) in air with a negative direct current corona discharge. A numerical model is used to theoretically investigate the underlying plasma chemistry for the removal of TCE, and a reaction pathway for the abatement of TCE is proposed. The Cl atom, mainly produced by dissociation of COCl, is one of the controlling species in the TCE destruction chemistry and contributes to the production of chlorine containing by-products. The effect of humidity on the removal efficiency is studied and a good agreement is found between experiments and the model for both dry (5 % relative humidity (RH)) and humid air (50 % RH). An increase of the relative humidity from 5 % to 50 % has a negative effect on the removal efficiency, decreasing by ±15 % in humid air. The main loss reactions for TCE are with ClO·, O· and CHCl2. Finally, the by-products and energy cost of TCE abatement are discussed.

Abstract: In this work, we use density functional theory calculations to study the combined effect of external electric fields, surface morphology, and surface charge on CO2 activation over Cu(111), Cu(211), Cu(110), and Cu(001) surfaces. We observe that the binding energy of the CO2 molecule on Cu surfaces increases significantly upon increasing the applied electric field strength. In addition, rougher surfaces respond more effectively to the presence of the external electric field toward facilitating the formation of a carbonate-like CO2 structure and the transformation of the most stable adsorption mode from physisorption to chemisorption. The presence of surface charges further strengthens the electric field effect and consequently causes an improved bending of the CO2 molecule and C−O bond length elongation. On the other hand, a net charge in the absence of an externally applied electric field shows only a marginal effect on CO2 binding. The chemisorbed CO2 is more stable and further activated when the effects of an external electric field, rough surface, and surface charge are combined. These results can help to elucidate the underlying factors that control CO2 activation in heterogeneous and plasma catalysis, as well as in electrochemical processes.

Abstract: Increasing the selectivity of cancer treatments is attractive, as it has the potential to reduce side-effects of therapy. Cold atmospheric plasma (CAP) is a novel cancer treatment that disrupts the intracellular oxidative balance. Several reports claim CAP treatment to be selective, but retrospective analysis of these studies revealed discrepancies in several biological factors and culturing methods. Before CAP can be conclusively stated as a selective cancer treatment, the importance of these factors must be investigated. In this study, we evaluated the influence of the cell type, cancer type, and cell culture medium on direct and indirect CAP treatment. Comparison of cancerous cells with their non-cancerous counterparts was performed under standardized conditions to determine selectivity of treatment. Analysis of seven human cell lines (cancerous: A549, U87, A375, and Malme-3M; non-cancerous: BEAS-2B, HA, and HEMa) and five different cell culture media (DMEM, RPMI1640, AM, BEGM, and DCBM) revealed that the tested parameters strongly influence indirect CAP treatment, while direct treatment was less affected. Taken together, the results of our study demonstrate that cell type, cancer type, and culturing medium must be taken into account before selectivity of CAP treatment can be claimed and overlooking these parameters can easily result in inaccurate conclusions of selectivity.

Abstract: Industrial production of HCN from NH3 and CH4 not only uses precious Pt or Pt−Rh catalysts but also requires extremely high temperatures (∼1600 K). From an energetic, operational, and safety perspective, a drastic decrease in temperature is highly desirable. Here, we report ammonia reforming of methane for the production of HCN and H2 at 673 K by the combination of CH4/NH3 plasma and a supported Cu/silicalite-1 catalyst. 30% CH4 conversion has been achieved with 79% HCN selectivity. Catalyst characterization and plasma diagnostics reveal that the excellent reaction performance is attributed to metallic Cu active sites. In addition, we propose a possible reaction pathway, viz. E-R reactions with N, NH, NH2, and CH radicals produced in the plasma, for the production of HCN, based on density functional theory calculations. Importantly, the Cu/silicalite-1 catalyst costs less than 5% of the commercial Pt mesh catalyst.

Abstract: Intercellular communication plays a crucial role in cancer, as well as other diseases, such as inflammation, tissue degeneration, and neurological disorders. One of the proteins responsible for this, are connexins (Cxs), which come together to form a hemichannel. When two hemichannels of opposite cells interact with each other, they form a gap junction (GJ) channel, connecting the intracellular space of these cells. They allow the passage of ions, reactive oxygen and nitrogen species (RONS), and signaling molecules from the interior of one cell to another cell, thus playing an essential role in cell growth, differentiation, and homeostasis. The importance of GJs for disease induction and therapy development is becoming more appreciated, especially in the context of oncology. Studies have shown that one of the mechanisms to control the formation and disruption of GJs is mediated by lipid oxidation pathways, but the underlying mechanisms are not well understood. In this study, we performed atomistic molecular dynamics simulations to evaluate how lipid oxidation influences the channel properties of Cx26 hemichannels, such as channel gating and permeability. Our results demonstrate that the Cx26 hemichannel is more compact in the presence of oxidized lipids, decreasing its pore diameter at the extracellular side and increasing it at the amino terminus domains, respectively. The permeability of the Cx26 hemichannel for water and RONS molecules is higher in the presence of oxidized lipids. The latter may facilitate the intracellular accumulation of RONS, possibly increasing oxidative stress in cells. A better understanding of this process will help to enhance the efficacy of oxidative stress-based cancer treatments.

Abstract: It is generally known that antidiabetic activity is associated with an increased level of glucose uptake in adipocytes and skeletal muscle cells. However, the role of exogenous reactive oxygen and nitrogen species (RONS) in muscle development and more importantly in glucose uptake is largely unknown. We investigate the effect of RONS generated by cold atmospheric plasma (CAP) in glucose uptake. We show that the glucose uptake is significantly enhanced in differentiated L6 skeletal muscle cells after CAP treatment. We also observe a significant increase of the intracellular Ca++ and ROS level, without causing toxicity. One of the possible reasons for an elevated level of glucose uptake as well as intracellular ROS and Ca++ ions is probably the increased oxidative stress leading to glucose transport.

Abstract: It is well established that lipid aldehydes (LAs) are able to increase the permeability of cell membranes and induce their rupture. However, it is not yet clear how LAs are distributed in phase-separated membranes (PSMs), which are responsible for the transport of selected molecules and intracellular signaling. Thus, we investigate here the distribution of LAs in a PSM by coarse-grained molecular dynamics simulations. Our results reveal that LAs derived from mono-unsaturated lipids tend to accumulate at the interface between the liquid-ordered/liquiddisordered domains, whereas those derived from poly-unsaturated lipids remain in the liquid-disordered domain. These results are important for understanding the effects caused by oxidized lipids in membrane structure, properties and organization.

Abstract: Laser-induced plasmas are transient systems rapidly aging in few nanoseconds of evolution. Time-of-flight spectrometry allowed studying initial plasma characteristics based on frozen translational degrees of freedom, hence overcoming intrinsic limitations of optical spectroscopy. Experimental ion velocity distributions were reconstructed as developed during the longitudinal plasma expansion. The obtained onset plasma temperatures are in the range of similar to 18-45 eV depending on the ablated metals. Also the ion angular spreads were found to be a function of ablated metal, e.g., the narrowest for Fe, the broadest for Al, due to different collisional coupling in the plasma population. (c) 2007 American Institute of Physics.

Abstract: Live-cell imaging with fluorescence microscopy is a powerful tool, especially in cancer research, widely-used for capturing dynamic cellular processes over time. However, light-induced toxicity (phototoxicity) can be incurred from this method, via disruption of intracellular redox balance and an overload of reactive oxygen species (ROS). This can introduce confounding effects in an experiment, especially in the context of evaluating and screening novel therapies. Here, we aimed to unravel whether phototoxicity can impact cellular homeostasis and response to non-thermal plasma (NTP), a therapeutic strategy which specifically targets the intracellular redox balance. We demonstrate that cells incorporated with a fluorescent reporter for live-cell imaging have increased sensitivity to NTP, when exposed to ambient light or fluorescence excitation, likely through altered proliferation rates and baseline intracellular ROS levels. These changes became even more pronounced the longer the cells stayed in culture. Therefore, our results have important implications for research implementing this analysis technique and are particularly important for designing experiments and evaluating redox-based therapies like NTP.

Abstract: Low-temperature plasma chemical kinetic models are particularly important to the plasma community. These models typically require dozens of inputs, especially rate coefficients. The latter are not always precisely known and it is not surprising that the error on the rate coefficient data can propagate to the model output. In this paper, we present a model that uses N = 400 different combinations of rate coefficients based on the uncertainty attributed to each rate coefficient, giving a good estimation of the uncertainty on the model output due to the rate coefficients. We demonstrate that the uncertainty varies a lot with the conditions and the type of output. Relatively low uncertainties (about 15%) are found for electron density and temperature, while the uncertainty can reach more than an order of magnitude for the population of the vibrational levels in some cases and it can rise up to 100% for the CO2 conversion. The reactions that are mostly responsible for the largest uncertainties are identified. We show that the conditions of pressure, gas temperature and power density have a great effect on the uncertainty and on which reactions lead to this uncertainty. In all the cases tested here, while the absolute values may suffer from large uncertainties, the trends observed in previous modeling work are still valid. Finally, in accordance with the work of Turner, a number of ‘good practices’ is recommended.