Abstract: Plasma-treated liquids have great potential for biomedical applications. However, insight into the underlying mechanisms and the exact chemistry is still scarce. In this study, we present the combination of a 0D chemical kinetics and a 2D fluid dynamics model to investigate the plasma treatment of a buffered water solution with the kINPen (R) plasma jet. Using this model, we calculated the gas and liquid flow profiles and the transport and chemistry of all species in the gas and the liquid phase. Moreover, we evaluated the stability of the reactive oxygen and nitrogen species after plasma treatment. We found that of all species, only H2O2, HNO2/NO2-, and HNO3/NO3- are stable in the buffered solution after plasma treatment. This is because both their production and loss processes in the liquid phase are dependent on short-lived radicals (e.g. OH, NO, and NO2). Apart from some discrepancy in the absolute values of the concentrations, which can be explained by the model, all general trends and observations in our model are in qualitative agreement with experimental data and literature.

Abstract: The application of coreshell SiSiO2 nanowires as nanoelectronic devices strongly depends on their structure, which is difficult to tune precisely. In this work, we investigate the formation of the coreshell nanowires at the atomic scale, by reactive molecular dynamics simulations. The occurrence of two temperature-dependent oxidation mechanisms of ultra-small diameter Si-NWs is demonstrated. We found that control over the Si-core radius and the SiOx (x ≤ 2) oxide shell is possible by tuning the growth temperature and the initial Si-NW diameter. Two different structures were obtained, i.e., ultrathin SiO2 silica nanowires at high temperature and Si core|ultrathin SiO2 silica nanowires at low temperature. The transition temperature is found to linearly decrease with the nanowire curvature. Finally, the interfacial stress is found to be responsible for self-limiting oxidation, depending on both the initial Si-NW radius and the oxide growth temperature. These novel insights allow us to gain control over the exact morphology and structure of the wires, as is needed for their application in nanoelectronics.

Abstract: In this review paper, we present several examples of reactive molecular dynamics simulations, which contribute to a better understanding of the underlying mechanisms in plasma medicine on the atomic scale. This includes the interaction of important reactive oxygen plasma species with the outer cell wall of both gram-positive and gram-negative bacteria, and with lipids present in human skin. Moreover, as most biomolecules are surrounded by a liquid biofilm, the behavior of these plasma species in a liquid (water) layer is presented as well. Finally, a perspective for future atomic scale modeling studies is given, in the field of plasma medicine in general, and for cancer treatment in particular.

Abstract: Cold atmospheric plasma in contact with solutions has many applications, but its chemistry contains many unknowns such as the undescribed reactions with solutes. By combining experiments and modelling, we report the first direct demonstration of the reaction of chloride with oxygen atoms in aqueous solutions exposed to cold plasma.

Abstract: Molecular dynamics simulations using the Brenner potential have been performed to investigate reaction mechanisms of various hydrocarbon radicals with low kinetic energies on amorphous hydrogenated carbon (a-C:H) surfaces and to simulate thin a-C:H film growth. Experimental data from an expanding thermal plasma setup were used as input for the simulations. The hydrocarbon reaction mechanisms were studied both during growth of the films and on a set of surface sites specific for a-C:H surfaces. Thin film growth was studied using experimentally detected growth species. It is found that the reaction mechanisms and sticking coefficients are dependent on the specific surface sites, and the structural properties of the growth radicals. Furthermore, it is found that thin a-C:H films can be densified using an additional H-flux towards the substrate.

Abstract: The modelling of a gliding arc discharge (GAD) is studied by means of the quasineutral (QN) plasma modelling approach. The model is first evaluated for reliability and proper description of a gliding arc discharge at atmospheric pressure, by comparing with a more elaborate non-quasineutral (NQN) plasma model in two different geometries – a 2D axisymmetric and a Cartesian geometry. The NQN model is considered as a reference, since it provides a continuous self-consistent plasma description, including the near electrode regions. In general, the results of the QN model agree very well with those obtained from the NQN model. The small differences between both models are attributed to the approximations in the derivation of the QN model. The use of the QN model provides a substantial reduction of the computation time compared to the NQN model, which is crucial for the development of more complex models in three dimensions or with complicated chemistries. The latter is illustrated for (i) a reverse vortex flow(RVF) GAD in argon, and (ii) a GAD in CO2. The RVF discharge is modelled in three dimensions and the effect of the turbulent heat transport on the plasma and gas characteristics is
discussed. The GAD model in CO2 is in a 1D geometry with axial symmetry and provides results for the time evolution of the electron, gas and vibrational temperature of CO2, as well as for the molar fractions of the different species.

Abstract: Boron neutron capture therapy (BNCT) is an extensively studied radiotherapeutic strategy for cancer treatment. BNCT is based on irradiation of malignant tumour cells with neutrons after uptake of a B-10 containing molecule. Alpha particles, locally produced by neutron irradiation kill the cancer cells. Important for ongoing research regarding cellular uptake and cytotoxicity of a large variety of B-10 containing molecules is the accurate determination of boron concentrations in cell cultures. In this work, the sample preparation for quantitative inductively coupled plasma mass spectrometry (ICP-MS) analysis on cell cultures was optimized. By making use of acid digestion combined with UV digestion, low detection limits (0.4 mu g L-1) and full recoveries of boron could be achieved while measurements were free of spectral and non-spectral interferences. Finally, cell-associated boron in the form of 4-borono-l-phenylalanine (l-BPA) in vascular endothelial cell cultures could be determined with ICP-MS as (1.26 +/- 0.10) x 10(9) boron atoms per cell. The developed method can prove its importance for further BNCT research and elemental analysis of cell cultures.

Abstract: In this letter we compare the effect of a radio-frequency sine, a low frequency sine, a rectangular and a pulsed dc voltage profile on the calculated electron production and power consumption in the dielectric barrier discharge. We also demonstrate using calculated potential distribution profiles of high time and space resolution how the pulsed dc discharge generates a secondary discharge pulse by deactivating the power supply.

Abstract: Although plasma catalysis is gaining increasing interest for various environmental applications, the underlying mechanisms are still far from understood. For instance, it is not yet clear whether and how plasma streamers can propagate in catalyst pores, and what is the minimum pore size to make this happen. As this is crucial information to ensure good plasma-catalyst interaction, we study here the mechanism of plasma streamer propagation in a catalyst pore, by means of a twodimensional particle-in-cell/Monte Carlo collision model, for various pore diameters in the nm range to μm-range. The so-called Debye length is an important criterion for plasma penetration into catalyst pores, i.e. a plasma streamer can penetrate into pores when their diameter is larger than the Debye length. The Debye length is typically in the order of a few 100 nm up to 1 μm at the conditions under study, depending on electron density and temperature in the plasma streamer. For pores in the range of ∼50 nm, plasma can thus only penetrate to some extent and at
very short times, i.e. at the beginning of a micro-discharge, before the actual plasma streamer reaches the catalyst surface and a sheath is formed in front of the surface. We can make plasma streamers penetrate into smaller pores (down to ca. 500 nm at the conditions under study) by increasing the applied voltage, which yields a higher plasma density, and thus reduces the Debye length. Our simulations also reveal that the plasma streamers induce surface charging of the catalyst pore sidewalls, causing discharge enhancement inside the pore, depending on pore diameter and depth.

Abstract: Plasma catalysis, a promising technology for conversion of CO2 into value-added chemicals near room temperature, is gaining increasing interest. A dielectric barrier discharge (DBD) plasma has attracted attention due to its simple design and operation at near ambient conditions, ease to implement catalysts in the plasma zone and upscaling ability to industrial applications. To improve its main drawbacks, being relatively low conversion and energy efficiency, a packing material is used in the plasma discharge zone of the reactor, sometimes decorated by a catalytic material. Nevertheless, the extent to which different properties of the packing material influence plasma performance is still largely unexplored and unknown. In this study, the particular effect of synthesis induced differences in the morphology of a TiO2 shell covering a SiO2 core packing material on the plasma conversion of CO2 is studied. TiO2 has been successfully deposited around 1.6–1.8 mm sized SiO2 spheres by means of spray coating, starting from aqueous citratoperoxotitanate(IV) precursors. Parameters such as concentration of the Ti(IV) precursor solutions and addition of a binder were found to affect the shells’ properties and surface morphology and to have a major impact on the CO2 conversion in a packed bed DBD plasma reactor. Core-shell SiO2@TiO2 obtained from 0.25 M citratoperoxotitante(IV) precursors with the addition of a LUDOX binder showed the highest CO2 conversion 37.7% (at a space time of 70 s corresponding to an energy efficiency of 2%) and the highest energy efficiency of 4.8% (at a space time of 2.5 s corresponding to a conversion of 3%).

Abstract: Mesoporous carbon, which can be replicated from mesoporous silica and whose surface is hydrophobic, can be an ideal catalyst for the esterification of fatty acids. Here we report an easy and low cost way to prepare sulfonic acid group-functionalized mesoporous carbon. A sample of calcined mesoporous silica SBA-15 was added to an aqueous sucrose solution followed by drying and calcination at different temperatures. In contrast to existing procedures, the obtained hybrid Si/C material was then first sulfonated in H2SO4, before the final removal of the silica template in order to stabilize the porous structure towards the liquid phase sulfonation treatment. Thus the silicacarbon composites, instead of the mesoporous carbon, were successfully sulfonated to introduce SO3H groups, while keeping the ordered mesoporous structure intact. The influence of carbonization temperature was investigated, suggesting an optimum temperature of 873 K. The SO3H group-functionalized mesoporous carbon, denoted as CMK-3-873-SO3H, was characterized by means of XRD, N2 physisorption, SEM, FT-IR, elemental analysis and TEM. It followed that a uniform mesoporous carbon was obtained with an average pore size of 3.89 nm, a specific surface of 807 m2/g and a SO3H group loading of 0.39 meq/g of dry material. Compared with other solid acid catalysts, the resulting material shows enhanced activity in the acid-catalyzed esterification of oleic acid with methanol, and can be used repeatedly. The increased catalytic performance is attributed to the hydrophobic surface and larger pore size of the new catalyst. It can effectively accommodate long chain fatty acids and reject formed water, making the active sites easily accessible.

Abstract: Mcl1 is a primary member of the Bcl–2 family—anti–apoptotic proteins (AAP)—that is overexpressed in several cancer pathologies. The apoptotic regulation is mediated through the binding of pro-apoptotic peptides (PAPs) (e.g., Bak and Bid) at the canonical hydrophobic binding groove (CBG) of Mcl1. Although all PAPs form amphipathic α-helices, their amino acid sequences vary to different degree. This sequence variation exhibits a central role in the binding partner selectivity towards different AAPs. Thus, constructing a novel peptide or small organic molecule with the ability to mimic the natural regulatory process of PAP is essential to inhibit various AAPs. Previously reported experimental binding free energies (BFEs) were utilized in the current investigation aimed to understand the mechanistic basis of different PAPs targeted to mMcl1. Molecular dynamics (MD) simulations used to estimate BFEs between mMcl1—PAP complexes using Molecular Mechanics-Generalized Born Solvent Accessible (MMGBSA) approach with multiple parameters. Predicted BFE values showed an excellent agreement with the experiment (R2 = 0.92). The van–der Waals (ΔGvdw) and electrostatic (ΔGele) energy terms found to be the main energy components that drive heterodimerization of mMcl1—PAP complexes. Finally, the dynamic network analysis predicted the allosteric signal transmission pathway involves more favorable energy contributing residues. In total, the results obtained from the current investigation may provide valuable insights for the synthesis of a novel peptide or small organic inhibitor targeting Mcl1.

Abstract: We investigate discharge contraction in a microwave plasma at sub-atmospheric pressure, operating in CO2 and CO2/CH4 mixtures. The rise of the electron number density with plasma contraction intensifies the gas heating in the core of the plasma. This, in turn, initiates fast core-periphery transport and defines the rate of thermal chemistry over plasma chemistry. In this context, power concentration describes the overall mechanism including plasma contraction and chemical kinetics. In a complex chemistry such as dry reforming of methane, transport of reactive species is essential to define the performance of the reactor and achieve the desired outputs. Thus, we couple experimental observations and thermodynamic calculations for model validation and understanding of reactor performance. Adding CH4 alters the thermodynamic properties of the mixture, especially the reactive component of the heat conductivity. The increase in reactive heat conductivity increases the pressure at which plasma contraction occurs, because higher rates of gas heating are required to reach the same temperature. In addition, we suggest that the predominance of heat conduction over convection is a key condition to observe the effect of heat conductivity on gas temperature.

Abstract: Plasma catalysis is emerging for plasma-assisted gas conversion
processes. However, the underlying mechanisms of plasma catalysis are poorly
understood. In this work, we present a 1D heterogeneous catalysis model with axial
dispersion (i.e., accounting for back-mixing and molecular diffusion of fluid elements in
the process stream in the axial direction), for plasma-catalytic NO production from
N2/O2 mixtures. We investigate the concentration and reaction rates of each species
formed as a function of time and position across the catalyst, in order to determine the
underlying mechanisms. To obtain insights into how the performance of the process
can be further improved, we also study how changes in the postplasma gas flow
composition entering the catalyst bed and in the operation conditions of the catalytic
stage affect the performance of NO production.

Abstract: Transforming CO2 into value-added chemicals is crucial to realizing a carbon–neutral economy, and plasma-based conversion, a Power-2-X technology, offers a promising route to realizing an efficient and scalable process. This paper investigates the effects of post-plasma placement of a converging–diverging nozzle in a vortex-stabilized 2.45 GHz CO2 microwave plasma reactor to increase energy efficiency and conversion. The CDN leads to a 21 % relative increase in energy efficiency (31 %) and CO2 conversion (13 %) at high flow rates and near-atmospheric conditions. The most significant performance improvement was seen at low flow rates and sub-atmospheric pressure (300 mbar), where energy efficiency was 23 % and conversion was 28 %, a 71 % relative increase over conditions without the CDN. Using CFD simulations, we found that the CDN produces a change in the flow geometry, leading to a confined temperature profile at the height of the plasma, and forced extraction of CO to the post-CDN region.