Abstract: Given the ecological problems associated to the CO2 emissions of fossil fuels, plasma technology has gained
interest for conversion of CO2 into value-added products. Microwave plasmas operating at atmospheric pressure
have proven to be especially interesting, due to the high gas temperatures inside the reactor (i.e. up to 6000 K)
allowing for efficient thermal dissociation of CO2 into CO and O2. However, the performance of these high
temperature plasmas is limited by recombination of CO back into CO2 once the gas cools down in the afterglow.
In this work, we computationally investigated several quenching nozzles, developed and experimentally tested
by Hecimovic et al., [1] for their ability to quickly cool the gas after the plasma, thereby quenching the CO
recombination reactions. Using a 3D computational fluid dynamics model and a quasi-1D chemical kinetics
model, we reveal that a reactor without nozzle lacks gas mixing between hot gas in the center and cold gas near
the reactor walls. Especially at low flow rates, where there is an inherent lack of convective cooling due to the
low gas flow velocity, the temperature in the afterglow remains high (between 2000 and 3000 K) for a relatively
long time (in the 0.1 s range). As shown by our quasi-1D chemical kinetics model, this results in a important loss
of CO due to recombination reactions. Attaching a nozzle in the effluent of the reactor induces fast gas quenching
right after the plasma. Indeed, it introduces (i) more convective cooling by forcing cool gas near the walls to mix
with hot gas in the center of the reactor, as well as (ii) more conductive cooling through the water-cooled walls of
the nozzle. Our model shows that gas quenching and the suppression of recombination reactions have more
impact at low flow rates, where recombination is the most limiting factor in the conversion process.

Abstract: Accurate barriers for rate controlling elementary reactions on metal surfaces are key to understanding, controlling, and predicting the rate of heterogeneously catalyzed processes. While barrier heights for gas phase reactions have been extensively benchmarked, dissociative chemisorption barriers for the reactions of molecules on metal surfaces have received much less attention. The first database called SBH10 and containing 10 entries was recently constructed based on the specific reaction parameter approach to density functional theory (SRP-DFT) and experimental results. We have now constructed a new and improved database (SBH17) containing 17 entries based on SRP-DFT and experiments. For this new SBH17 benchmark study, we have tested three algorithms (high, medium, and light) for calculating barrier heights for dissociative chemisorption on metals, which we have named for the amount of computational effort involved in their use. We test the performance of 14 density functionals at the GGA, GGA+vdW-DF, and meta-GGA rungs. Our results show that, in contrast with the previous SBH10 study where the BEEF-vdW-DF2 functional seemed to be most accurate, the workhorse functional PBE and the MS2 density functional are the most accurate of the GGA and meta-GGA functionals tested. Of the GGA+vdW functionals tested, the SRP32-vdW-DF1 functional is the most accurate. Additionally, we found that the medium algorithm is accurate enough for assessing the performance of the density functionals tested, while it avoids geometry optimizations of minimum barrier geometries for each density functional tested. The medium algorithm does require metal lattice constants and interlayer distances that are optimized separately for each functional. While these are avoided in the light algorithm, this algorithm is found not to give a reliable description of functional performance. The combination of relative ease of use and demonstrated reliability of the medium algorithm will likely pave the way for incorporation of the SBH17 database in larger databases used for testing new density functionals and electronic structure methods.

Abstract: A crucial property for implants is their biocompatibility. To ensure biocompatibility, thin coatings of hydroxyapatite (HA) are deposited on the actual implant. In this study, we investigate the effects of the addition of silicate anions to the structure of hydroxyapatite coatings on their adhesion strength via a scratch test and ab initio calculations. We find that both the grain size and adhesion strength decrease with the increase in the silicon content in the HA coating (SiHA). The increase in the silicon content to 1.2 % in the HA coating leads to a decrease in the average crystallite size from 28 to 21 nm, and in the case of 4.6 %, it leads to the formation of an amorphous or nanocrystalline film. The decreases in the grain and crystallite sizes lead to peeling and destruction of the coating from the titanium substrate at lower loads. Further, our ab initio simulations demonstrate an increased number of molecular bonds at the amorphous SiHA-TiO2 interface. However, the experimental results revealed that the structure and grain size have more pronounced effects on the adhesion strength of the coatings. In conclusion, based on the results of the ab initio simulations and the experimental results, we suggest that the presence of Si in the form of silicate ions in the HA coating has a significant impact on the structure, grain size, and number of molecular bonds at the interface and on the adhesion strength of the SiHA coating to the titanium substrate.

Abstract: Laser ablation is used in a growing number of applications in various areas including medicine, archaeology, chemistry, environmental and materials sciences. In this work the heat transfer and phase change phenomena during nanosecond laser ablation of a copper (Cu) target in a helium (He) background gas at atmospheric pressure are presented. An enthalpy model is outlined, which accounts for heating, melting, and vaporization of the target. As far as we know, this is the first model that connects the thermodynamics and underlying kinetics of this challenging phase change problem in a selfconsistent way.

Abstract: Cancer immunotherapy is a promising strategy that engages the patient's immune system to kill cancer cells selectively while sparing normal tissue. Treatment of macrophages with a nanosecond-pulsed dielectric barrier discharge directly enhanced their cytotoxic activity against tumor cells but not normal cells. These results underscore the clinical potential of plasma for cancer immunotherapy.

Abstract: Studies using xenograft mouse models have shown that plasma applied to the skin overlying tumors results in tumor shrinkage. Plasma is considered a nonpenetrating treatment; however, these studies demonstrate plasma effects that occur beyond the postulated depth of physical penetration of plasma components. The present study examines the propagation of plasma effects through a tissue model using three-dimensional, cell-laden extracellular matrices (ECMs). These ECMs are used as barriers against direct plasma penetration. By placing them onto a monolayer of target cancer cells to create an in-vitro analog to in-vivo studies, we distinguished between cellular effects from direct plasma exposure and cellular effects due to cell-to-cell signaling stimulated by plasma. We show that nanosecond-pulsed dielectric barrier discharge plasma treatment applied atop an acellular barrier impedes the externalization of calreticulin (CRT) in the target cells. In contrast, when a barrier is populated with cells, CRT externalization is restored. Thus, we demonstrate that plasma components stimulate signaling among cells embedded in the barrier to transfer plasma effects to the target cells.

Abstract: The evaluation of nucleation rates from molecular dynamics trajectories is hampered by the slow nucleation time scale and impact of finite size effects. Here, we show that accurate nucleation rates can be obtained in a very general fashion relying only on the free energy barrier, transition state theory, and a simple dynamical correction for diffusive recrossing. In this setup, the time scale problem is overcome by using enhanced sampling methods, in casu metadynamics, whereas the impact of finite size effects can be naturally circumvented by reconstructing the free energy surface from an appropriate ensemble. Approximations from classical nucleation theory are avoided. We demonstrate the accuracy of the approach by calculating macroscopic rates of droplet nucleation from argon vapor, spanning 16 orders of magnitude and in excellent agreement with literature results, all from simulations of very small (512 atom) systems.

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: Packed-bed DBD (PB-DBD) plasmas hold promise for effective degradation of greenhouse gases like SF6. In this work, we conducted a combined experimental and theoretical study to investigate the effect of the packing surface structure and the plasma surface discharge on the SF6 degradation in a gamma-Al2O3 packing DBD system. Experimental results show that both the hydration effect of the surface (upon moisture) and the presence of excessive reactive gases in the plasma can significantly reduce the SF6 degradation, but they hardly change the discharge behavior. DFT results show that the pre-adsorption of species such as H, OH, H2O and O-2 can occupy the active sites (Al-III site) which negatively impacts the SF6 adsorption. H2O molecules pre-adsorbed at neighboring sites can promote the activation of SF6 molecules and lower the reaction barrier for the S-F bond-breaking process. Surface-induced charges and local external electric fields caused by the plasma can both improve the SF6 adsorption and enhance the elongation of the S-F bonds. Our results indicate that both the surface structure of the packing material and the plasma surface discharge are crucial for SF6 degradation performance, and the packing beads should be kept dry during the degradation. This work helps to understand the underlying mechanisms of SF6 degradation in a PB-DBD system.

Abstract: Plasma catalysis has emerged as a promising approach for driving thermodynamically unfavorable chemical
reactions. Nevertheless, comprehending the mechanisms involved remains a challenge, leading to uncertainty
about whether the optimal catalyst in plasma catalysis aligns with that in thermal catalysis. In this research, we
explore this question by studying monometallic catalysts (Fe, Co, Ni and Mo) and bimetallic catalysts (Fe-Co, Mo-
Co, Fe-Ni and Mo-Ni) in both thermal catalytic and plasma catalytic NH3 decomposition. Our findings reveal that
the Fe-Co bimetallic catalyst exhibits the highest activity in thermal catalysis, the Fe-Ni bimetallic catalyst
outperforms others in plasma catalysis, indicating a discrepancy between the optimal catalysts for the two
catalytic modes in NH3 decomposition. Comprehensive catalyst characterization, kinetic analysis, temperature
program surface reaction experiments and plasma diagnosis are employed to discuss the key factors influencing
NH3 decomposition performance.

Abstract: Carbon nanotubes (CNTs) have been successfully used in biomedicine, including cancer therapy, due to their unique physico-chemical properties. Because pristine CNTs exhibit hydrophobic behaviour, they can have a cytotoxic effect on cells, which limits their practical use in biomedicine. The toxicity of CNTs can be reduced by adding water-soluble functional radicals to their surface, i.e. by increasing their hydrophilicity. Another possibility for increasing the hydrophilicity of CNTs is probably filling them with endohedral metal atoms, which has not yet been studied. Thus, in this study, we use computer simulations to investigate the combined effect of endohedral nickel atoms and functional groups on the hydrophilicity of CNTs. Our simulation results show that the introduction of endohedral nickel atoms into CNTs increases their binding energy with functional groups. We also find that the addition of functional groups to the surface of CNT, along with filling it with endohedral nickel atoms, leads to an increase in the dipole moment of the CNT as well as its interaction energy with water, thereby increasing the hydrophilicity of the CNT and, consequently, its solubility in water. This, in turn, can lead to a decrease in CNT toxicity.

Abstract: The overexpression of voltage dependent anion channels (VDACs), particularly VDAC1, in cancer cells compared to normal cells, plays a crucial role in cancer cell metabolism, apoptosis regulation, and energy homeostasis. In this study, we used molecular dynamics (MD) simulations to investigate the effect of a low level of VDAC1 oxidation (induced e.g., by cold atmospheric plasma (CAP)) on the pyruvate (Pyr) uptake by VDAC1. Inhibiting Pyr uptake through VDAC1 can suppress cancer cell proliferation. Our primary target was to study the translocation of Pyr across the native and oxidized forms of hVDAC1, the human VDAC1. Specifically, we employed MD simulations to analyze the hVDAC1 structure by modifying certain cysteine residues to cysteic acids and methionine residues to methionine sulfoxides, which allowed us to investigate the effect of oxidation. Our results showed that the free energy barrier for Pyr translocation through the native and oxidized channel was approximately 4.3 +/- 0.7 kJ mol-1 and 10.8 +/- 1.8 kJ mol-1, respectively. An increase in barrier results in a decrease in rate of Pyr permeation through the oxidized channel. Thus, our results indicate that low levels of CAP oxidation reduce Pyr translocation, resulting in decreased cancer cell proliferation. Therefore, low levels of oxidation are likely sufficient to treat cancer cells given the inhibition of Pyr uptake.

Abstract: In this letter, a dipole sheet method is proposed to theoretically study the adsorption and decomposition of C5F10O over-discharged Cu (111) and Al (111) surfaces. A synergistic effect of external electric fields and surface excess charges shows up for jointly promoting the adsorption of C5F10O, accompanied by the enhancement of C-F bond elongation and charge transfer process. The decomposition of C5F10O is facilitated in the discharged region and the initial decomposition is found most likely to occur via the cleavage of the C-F single bond. The results indicate that the decomposition of C5F10O over the metal electrode surfaces is much accelerated when discharge faults occur and free F atoms could be generated from C5F10O before its carbon chain breakage. These findings help to elucidate the underlying decomposition tendency of C5F10O in discharged systems and provide a practical method for evaluating and designing new insulation gases.

Abstract: The unique physical properties of carbyne, a novel carbon nanostructure, have attracted considerable interest in modern nanotechnology. While carbyne synthesis has been accomplished successfully using diverse techniques, the underlying mechanisms governing the carbon monoxide-dependent catalytic synthesis of endohedral carbyne remain poorly understood. In this simulation-based study, we investigate the synthesis of endohedral carbyne from carbon and carbon monoxide radicals in the presence of a nickel catalyst inside double-walled carbon nanotubes with a (5,5)@(10,10) structure. The outcome of our investigation demonstrates that the incorporation of the carbon atom within the Ni-n@(5,5)@(10,10) model system initiates the formation of an elongated carbon chain. In contrast, upon the introduction of carbon monoxide radicals, the growth of the carbyne chain is inhibited as a result of the oxidation of endohedral nickel clusters by oxygen atoms after the initial steps of nucleation. Our findings align with prior theoretical, simulation, and experimental investigations, reinforcing their consistency and providing valuable insights into the synthesis of carbyne-based nanodevices that hold promising potential for future advancements in nanotechnology.

Abstract: The development of compact, low power, charge-neutral propulsion sources is of significant recent interest due to the rising application of micro-scale satellite platforms. Among such sources, radio frequency (rf) electrothermal microthrusters present an attractive option due to their scalability, reliability, and tunable control of power coupling to the propellant. For micropropulsion applications, where available power is limited, it is of particular importance to understand how electrical power can be transferred to the propellant efficiently, a process that is underpinned by the plasma sheath dynamics. In this work, two-dimensional fluid/Monte Carlo simulations are employed to investigate the effects of applied voltage frequency on the electron, ion, and neutral heating in an rf capacitively coupled plasma microthruster operating in argon. Variations in the electron and argon ion densities and power deposition, and their consequent effect on neutral-gas heating, are investigated with relation to the phase-averaged and phase-resolved sheath dynamics for rf voltage frequencies of 6-108 MHz at 450 V. Driving voltage frequencies above 40.68 MHz exhibit enhanced volumetric ionization from bulk electrons at the expense of the ion heating efficiency. Lower driving voltage frequencies below 13.56 MHz exhibit more efficient ionization due to secondary electrons and an increasing fraction of rf power deposition into ions. Thermal efficiencies are improved by a factor of 2.5 at 6 MHz as compared to the more traditional 13.56 MHz, indicating a favorable operating regime for low power satellite applications.