Abstract: Plasma degradation is one of the most effective methods for the abatement of greenhouse gas sulfur hexafluoride
(SF6). To evaluate the potential of MgO as a catalyst in plasma degradation, we investigate the catalytic properties
of MgO on SF6 adsorption and activation by density functional theory (DFT) where the O-defective and
hydroxylated surfaces are considered as two typical plasma-generated surfaces. Our results show that perfect
MgO (001) and (111) surfaces cannot interact with SF6 and only physical adsorption happens. In case of Odefective
MgO surfaces, the O vacancy is the most stable adsorption site. SF6 undergoes a decomposition to SF5
and F over the O-defective MgO (001) surface and undergoes an elongation of the bottom S-F bond over the Odefective
(111) surface. Besides, SF6 shows a physically adsorption at the stepsite of the MgO (001) surface,
accompanied by small changes in its bond angle and length. Furthermore, SF6 is found to be physically and
chemically adsorbed over 0.5 and 1.0 ML (monolayer) H-covered O-terminated MgO (111) surfaces, respectively.
The SF6 molecule undergoes a self-decomposition on the 1.0 ML hydroxylated surface via a surface bonding
process. This study shows that defective and hydroxylated MgO surfaces have the surface capacities for SF6
activation, which shows that MgO has potential as packing material in SF6 waste treatment in packed-bed
plasmas.

Abstract: A combination of a gliding arc plasmatron (GAP) reactor and a newly designed tubular catalyst bed (N-bed) was applied to investigate the post-plasma catalytic (PPC) effect for dry reforming of methane (DRM). As comparison, a traditional plasma catalyst bed (T-bed) was also utilized. The post-plasma catalytic effect of a Ni-based mixed oxide (Ni/MO) catalyst with a thermal catalytic performance of 77% CO2 and 86% CH4 conversion at 700 ℃ was studied. Although applying the T-bed had little effect on plasma based CO2 and CH4 conversion, an increase in selectivity to H2 was obtained with a maximum value of 89% at a distance of 2 cm. However, even when only α-Al2O3 packing material was used in the N-bed configuration, compared to the plasma alone and the T-bed, an increase of the CO2 and CH4 conversion from 53% and 53% to 69% and 69% to 83% was achieved. Addition of the Ni/MO catalyst further enhanced the DRM reaction, resulting in conversions of 79% for CO2 and 91% for
CH4. Hence, although no insulation nor external heating was applied to the N-bed post plasma, it provides a slightly better conversion than the thermal catalytic performance with the same catalyst, while being fully electrically driven. In addition, an enhanced CO selectivity to 96% was obtained and the energy cost was reduced from ~ 6 kJ/L (plasma alone) to 4.3 kJ/L. To our knowledge, it is the first time that a post-plasma catalytic system achieves this excellent catalytic performance for DRM without extra external heating or insulation.

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: We perform first-principles calculations to investigate the preferred positions of B and P dopants, both neutral and in their preferred charge state, in Si and Si/SiO2 core-shell nanowires (NWs). In order to understand the observed trends in the formation energy, we isolate the different effects that determine these formation energies. By making the distinction between the unrelaxed and the relaxed formation energy, we separate the impact of the relaxation from that of the chemical environment. The unrelaxed formation energies are determined by three effects: (i) the effect of strain caused by size mismatch between the dopant and the host atoms, (ii) the local position of the band edges, and (iii) a screening effect. In the case of the SiNW (Si/SiO2 NW), these effects result in an increase of the formation energy away from the center (interface). The effect of relaxation depends on the relative size mismatch between the dopant and host atoms. A large size mismatch causes substantial relaxation that reduces the formation energy considerably, with the relaxation being more pronounced towards the edge of the wires. These effects explain the surface segregation of the B dopants in a SiNW, since the atomic relaxation induces a continuous drop of the formation energy towards the edge. However, for the P dopants, the formation energy starts to rise when moving from the center but drops to a minimum just next to the surface, indicating a different type of behavior. It also explains that the preferential location for B dopants in Si/SiO2 core-shell NWs is inside the oxide shell just next to the interface, whereas the P dopants prefer the positions next to the interface inside the Si core, which is in agreement with recent experiments. These preferred locations have an important impact on the electronic properties of these core-shell NWs. Our simulations indicate the possibility of hole gas formation when B segregates into the oxide shell.

Abstract: Plasma catalysis is gaining increasing interest for various environmental applications. Catalytic
material can be inserted in different shapes in the plasma, e.g., as pellets, (coated) beads, but also
as honeycomb monolith and 3DFD structures, also called ‘structured catalysts’, which have high
mass and heat transfer properties. In this work, we examine the streamer discharge propagation
and the interaction between plasma and catalysts, inside the channels of such structured catalysts,
by means of a two-dimensional particle-in-cell/Monte Carlo collision model. Our results reveal
that plasma streamers behave differently in various structured catalysts. In case of a honeycomb
structure, the streamers are limited to only one channel, with low or high plasma density when
the channels are parallel or perpendicular to the electrodes, respectively. In contrast, in case of a
3DFD structure, the streamers can distribute to different channels, causing discharge
enhancement due to surface charging on the dielectric walls of the structured catalyst, and
especially giving rise to a broader plasma distribution. The latter should be beneficial for plasma
catalysis applications, as it allows a larger catalyst surface area to be exposed to the plasma.

Abstract: Nowadays, lasers are commonly applied in spectrochemical analysis methods, for sampling, plasma formation or a combination of both. Despite the numerous investigations that have been performed on these applications, the underlying processes are still insufficiently understood. In order to fasten progress in the field and in honor of the lifework of professor Rick Russo, we here provide a brief overview of the fundamental mechanisms in lasermatter interaction as proposed in literature, and throw the spotlight on some aspects that have not received much attention yet. For an organized discussion, we choose laser ablation, laser desorption and the associated gaseous plasma formation as the central processes in this perspective article, based on a classification of the laserbased spectrochemical analysis techniques and the corresponding laser-matter interaction regimes. First, we put the looking glass over the excitation and thermalization mechanisms in the laser-irradiated condensed phase, for which we propose the so-called multi-plasma model. This novel model can be understood as an extension of the well-known two-temperature model, featuring multiple thermodynamic dimensions, each of which corresponds to a quasi-particle type. Next, the focus is placed on the mass transfer and ionization mechanisms, after which we shortly highlight the possible role of anisotropic and magnetic effects in the laser-excited material. We hope this perspective article motivates more fundamental research on laser-matter interaction, as a continuation of the lifework of Rick Russo.

Abstract: Piezoelectricity is considered to be one of the key functionalities in biomaterials to boost bone tissue regeneration, however, integrating biocompatibility, biodegradability and 3D structure with pronounced piezoresponse remains a material challenge. Herein, novel hybrid biocompatible 3D scaffolds based on biodegradable poly(3-hydroxybutyrate) (PHB) and reduced graphene oxide (rGO) flakes have been developed. Nanoscale insights revealed a more homogenous distribution and superior surface potential values of PHB fibers (33 +/- 29 mV) with increasing rGO content up to 1.0 wt% (314 +/- 31 mV). The maximum effective piezoresponse was detected at 0.7 wt% rGO content, demonstrating 2.5 and 1.7 times higher out-of-plane and in-plane values, respectively, than that for pure PHB fibers. The rGO addition led to enhanced zigzag chain formation between paired lamellae in PHB fibers. In contrast, a further increase in rGO content reduced the alpha-crystal size and prevented zigzag chain conformation. A corresponding model explaining structural and molecular changes caused by rGO addition in electrospun PHB fibers is proposed. In addition, finite element analysis revealed a negligible vertical piezoresponse compared to lateral piezoresponse in uniaxially oriented PHB fibers based on alpha-phase (P2(1)2(1)2(1) space group). Thus, the present study demonstrates promising results for the development of biodegradable hybrid 3D scaffolds with an enhanced piezoresponse for various tissue engineering applications.

Abstract: Biological membranes are under constant attack of free radicals, which may lead to lipid nitro-oxidation, pro­ ducing a complex mixture of nitro-oxidized lipids that are responsible for structural and dynamic changes on the membrane. Despite the latter, nitro-oxidized lipids are also associated with several inflammatory and neuro­ degenerative diseases, the underlying mechanisms of which remain elusive. We perform atomistic molecular dynamics simulations using several isomers of nitro-oxidized lipids to study their effect on the structure and permeability of the membrane, as well as the interaction between the mixture of these products in the phos­pholipid membrane environment. Our results show that the stereo- and positional isomers have a stronger effect on the properties of the membrane composed of oxidized lipids compared to that containing nitrated lipids. Nevertheless, nitrated lipids lead to three-fold increase in water permeability compared to oxidized lipids. In addition, we show that in a membrane consisting of combined nitro-oxidized lipid products, the presence of oxidized lipids protects the membrane from transient pores. Is well stablished that plasma application and photodynamic therapy produces a number of oxidative species used to kill cancer cells, through membrane damage induced by nitro-oxidative stress. This study is important to elucidate the mechanisms and the molecular level properties involving the reactive species produced during that cancer therapies.

Abstract: The limitation in time and length scale is a major issue of molecular dynamics (MD) simulation. Although several methods have been developed to extend the MD time scale, their performance usually deteriorates with increasing system size. Therefore, an acceleration method which is applicable to large systems is required to bridge the gap between the MD simulations and target phenomena. In this study, an accelerated MD method for large system is developed based on the collective variable-driven hyperdynamics (CVHD) method [K.M. Bal and E.C. Neyts, 2015]. The key idea is to run CVHD in parallel with rate control and accelerate multiple possible events simultaneously. Using this novel method, carbon diffusion in bcc-iron bicrystal with grain boundary is examined as an application for practical materials. Carbon atoms reaching at the grain boundary are trapped whereas carbon atoms in the bulk region diffuse randomly, and both dynamic regimes can be simultaneously accelerated with the parallel CVHD technique.

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 metallized film capacitors, there exists an air gap of about 0.2 μm between the films, with a pressure ranging generally from 130 atm. Because of the created potential difference between the two films, a microdischarge is formed in this gap. In this paper, we use an implicit particle-in-cell Monte Carlo collision simulation method to study the discharge properties in this direct-current microdischarge with 0.2 μm gap in a range of different voltages and pressures. The discharge process is significantly different from a conventional high pressure discharge. Indeed, the high electric field due to the small gap sustains the discharge by field emission. At low applied voltage (~15 V), only the electrons are generated by field emission, while both electrons and ions are generated as a stable glow discharge at medium applied voltage (~50 V). At still higher applied voltage (~100 V), the number of electrons and ions rapidly multiplies, the electric field reverses, and the discharge changes from a glow to an arc regime.

Abstract: Molecular dynamics simulations using improved TersoffBrenner potential parameters were performed to investigate Cl+ etching of a {2×1} reconstructed Si(100) surface. Steady-state Si etching accompanying the Cl coverage of the surface is observed. Furthermore, a steady-state chlorinated reaction layer is formed. The thickness of this reaction layer is found to increase with increasing energy. The stoichiometry of SiClx species in the reaction layer is found to be SiCl:SiCl2:SiCl3 = 1.0:0.14:0.008 at 50 eV. These results are in excellent agreement with available experimental data. While elemental Si products are created by physical sputtering, most SiClx (0<x<4) etch products are produced by chemical-enhanced physical sputtering.

Abstract: Superimposition of dual frequencies (DFs) is one of the methods used for controlling plasma distribution in an inductively coupled plasma (ICP) source. The effects of a superimposed DF on the argon plasma characteristics have been investigated using a two-dimensional self-consistent fluid model. When both currents are fixed at 6A, the plasma density drops with decrease in one of the source frequencies due to less efficient heating and the plasma uniformity improves significantly. Moreover, for ICP operated with superimposed DFs (i.e., 4.52MHz/13.56MHz and 2.26MHz/13.56MHz), the current source exhibits the same period as the low frequency (LF) component, and the plasma density is higher than that obtained at a single frequency (i.e., 4.52 and 2.26MHz) with the same total current of 12A. However, at superimposed current frequencies of 6.78MHz/13.56MHz, the plasma density is lower than that obtained at a single frequency of 6.78MHz due to the weaker negative azimuthal electric field between two positive maxima during one period of 6.78MHz. When the superimposed DF ICP operates at 2.26 and 13.56MHz, the rapid oscillations of the induced electric field become weaker during one period of 2.26MHz as the current ratio of 2.26MHz/13.56MHz rises from 24A/7 A to 30A/1 A, and the plasma density drops with the current ratio due to weakened electron heating. The uniformity of plasma increases due to sufficient diffusion under the low-density condition.

Abstract: In this letter, CMOS-compatible Ni/HfO2/TiN resistive random access memory stacks demonstrated attractive unipolar switching properties, showing >10(3) endurance and long retention at 150 degrees C. The Ni bottom electrode (BE) improved the switching yield over the NiSiPt BE. To better understand the unipolar forming mechanism, ab initio simulation and time of flight-secondary ion mass spectroscopy were utilized. Compared to the NiSiPt BE, Ni BE gives larger Ni diffusion in the HfO2 and lower formation enthalpy of Ni2+ species during electrical forming. Both the electrical and physical results supported a Ni-injection mechanism for the filament formation. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.3695078]

Abstract: Atomistic simulation methods for the quantification of free energies are in wide use. These methods operate by sampling the probability density of a system along a small set of suitable collective variables (CVs), which is, in turn, expressed in the form of a free energy surface (FES). This definition of the FES can capture the relative stability of metastable states but not that of the transition state because the barrier height is not invariant to the choice of CVs. Free energy barriers therefore cannot be consistently computed from the FES. Here, we present a simple approach to calculate the gauge correction necessary to eliminate this inconsistency. Using our procedure, the standard FES as well as its gauge-corrected counterpart can be obtained by reweighing the same simulated trajectory at little additional cost. We apply the method to a number of systems—a particle solvated in a Lennard-Jones fluid, a Diels–Alder reaction, and crystallization of liquid sodium—to demonstrate its ability to produce consistent free energy barriers that correctly capture the kinetics of chemical or physical transformations, and discuss the additional demands it puts on the chosen CVs. Because the FES can be converged at relatively short (sub-ns) time scales, a free energy-based description of reaction kinetics is a particularly attractive option to study chemical processes at more expensive quantum mechanical levels of theory.