Abstract: Atomistic simulations can in principle provide an unbiased description of all mechanisms, intermediates, and products of complex chemical processes. However, due to the severe time scale limitation of conventional simulation techniques, unrealistically high simulation temperatures are usually applied, which are a poor approximation of most practically relevant low-temperature applications. In this work, we demonstrate the direct observation at the atomic scale of the pyrolysis and oxidation of n-dodecane at temperatures as low as 700 K through the use of a novel simulation technique, collective variable-driven hyperdynamics (CVHD). A simulated timescale of up to 39 seconds is reached. Product compositions and dominant mechanisms are found to be strongly temperature-dependent, and are consistent with experiments and kinetic models. These simulations provide a first atomic-level look at the full dynamics of the complicated fuel combustion process at industrially relevant temperatures and time scales, unattainable by conventional molecular dynamics simulations.

Abstract: van der Waals (vdW) forces play an important role in the adsorption of molecules on the surface of solids. However, the choice of the most suitable vdW functional for different systems is an essential problem which must be addressed for different systems. The lack of a systematic study on the performance of the vdW functionals in the bulk and adsorption properties of metal-oxides motivated us to examine different vdW approaches and compute the bulk and molecular adsorption properties of alpha-Cr2O3, alpha-Fe2O3, and alpha-Al2O3. For the bulk properties, we compared our results for the heat of formation, cohesive energy, lattice parameters and bond distances between the different vdW functionals and available experimental data. Next we studied the adsorption of benzene and CH3 molecules on top of different oxide surfaces. We employed different approximations to exchange and correlation within DFT, namely, the Perdew-Burke-Ernzerhof (PBE) GGA, (PBE)+U, and vdW density functionals [ DFT(vdW-DF/DF2/optPBE/optB86b/optB88)+U] as well as DFT-D2/D3(+U) methods of Grimme for the bulk calculations and optB86b-vdW(+U) and DFT-D2(+U) for the adsorption energy calculations. Our results highlight the importance of vdW interactions not only in the adsorption of molecules, but importantly also for the bulk properties. Although the vdW contribution in the adsorption of CH3 (as a chemisorption interaction) is less important compared to the adsorption of benzene (as a physisorption interaction), this contribution is not negligible. Also adsorption of benzene on ferryl/chromyl terminated surfaces shows an important chemisorption contribution in which the vdW interactions become less significant.

Abstract: To investigate the impact of both the CoFe ferromagnetic layer thickness and the capping paramagnetic layer on the tunnel magnetoresistance (TMR), we performed first-principles simulations on epitaxial magnetic tunnel junctions contacted with either CoFe or Ta paramagnetic capping layers. We observed a strong oscillation of the TMR amplitude with respect to the thickness of the ferromagnetic layer. The TMR is found to be amplified whenever the MgO spin tunnel barrier is thickened. Quantization of the electronic structure of the ferromagnetic layers is found to be at the origin of this oscillatory behavior. Metals such as Ta contacting the magnetic layer are found to enhance the amplitude of the oscillations due to the occurrence of an interface dipole. The latter drives the band alignment and tunes the nature of the spin channels that are active during the tunneling process. Subsequently, the regular transmission spin channels are modulated in the magnetic tunnel junction stack and other complex ones are being activated.

Abstract: Experiments show that the etch rate of Si with SF6 inductively coupled plasma (ICP) is significantly influenced by the absolute gas flow rate in the range of 50–600 sccm, with a maximum at around 200 sccm. Therefore, we numerically investigate the effects of the gas flow rate on the bulk plasma properties and on the etch rate, to obtain more insight in the underlying reasons of this effect. A hybrid Monte Carlo—fluid model is applied to simulate an SF6 ICP. It is found that the etch rate is influenced by two simultaneous effects: (i) the residence time of the gas and (ii) the temperature profile of the plasma in the ICP volume, resulting indeed in a maximum etch rate at 200 sccm.

Abstract: Plasma–liquid interactions represent a growing interdisciplinary area of research involving plasma science, fluid dynamics, heat and mass transfer, photolysis, multiphase chemistry and aerosol science. This review provides an assessment of the state-of-the-art of this multidisciplinary area and identifies the key research challenges. The developments in diagnostics, modeling and further extensions of cross section and reaction rate databases that are necessary to address these challenges are discussed. The review focusses on nonequilibrium plasmas.

Abstract: Electrical breakdown by the application of an electric field occurs more easily in hot gases than in cold gases because of the extra electron-species interactions that occur as a result of dissociation, ionization and excitation at higher temperature. This paper discusses some overlooked physics and clarifies inaccuracies in the evaluation of the effective ionization coefficients and the critical reduced breakdown electric field of CO2 at elevated temperature, considering the influence of excited states and ion kinetics. The critical reduced breakdown electric field is obtained by balancing electron generation and loss mechanisms using the electron energy distribution function (EEDF) derived from the Boltzmann transport equation under the two-term approximation. The equilibrium compositions of the hot gas mixtures are determined based on Gibbs free energy minimization considering the ground states as well as vibrationally and electronically excited states as independent species, which follow a Boltzmann distribution with a fixed excitation temperature. The interaction cross sections between electrons and the excited species, not reported previously, are properly taken into account. Furthermore, the ion kinetics, including electron–ion recombination, associative electron detachment, charge transfer and ion conversion into stable negative ion clusters, are also considered. Our results indicate that the excited species lead to a greater population of high-energy electrons at higher gas temperature and this affects the Townsend rate coefficients (i.e. of electron impact ionization and attachment), but the critical reduced breakdown electric field strength of CO2 is only affected when also properly accounting for the ion kinetics. Indeed, the latter greatly influences the effective ionization coefficients and hence the critical reduced breakdown electric field at temperatures above 1500 K. The rapid increase of the dissociative electron attachment cross-section of molecular oxygen with rising vibrational quantum number leads to a larger electron loss rate and this enhances the critical reduced breakdown electric field strength in the temperature range where the concentration of molecular oxygen is relatively high. The results obtained in this work show reasonable agreement with experimental results from literature, and are important for the evaluation of the dielectric strength of CO2 in a highly reactive environment at elevated temperature.

Abstract: The idea of stacking multiple monolayers of different two-dimensional materials has become a global pursuit. In this work, a silicene armchair nanoribbon of width W and van der Waals-bonded to different transition-metal dichalcogenide (TMD) bilayer substrates MoX2 and WX2, where X = S, Se, Te is considered. The orbital resolved electronic structure and ballistic transport properties of these systems are simulated by employing van der Waals-corrected density functional theory and nonequilibrium Green's functions. We find that the lattice mismatch with the underlying substrate determines the electronic structure, correlated with the silicene buckling distortion and ultimately with the contact resistance of the two-terminal system. The smallest lattice mismatch, obtained with the MoTe2 substrate, results in the silicene ribbon properties coming close to those of a freestanding one. With the TMD bilayer acting as a dielectric layer, the electronic structure is tunable from a direct to an indirect semiconducting layer, and subsequently to a metallic electronic dispersion layer, with a moderate applied perpendicular electric field.

Abstract: Plasma catalysis is gaining increasing interest for various environmental applications, but the crucial question is whether plasma can be created inside catalyst pores and under which conditions. In practice, various catalytic support materials are used, with various dielectric constants. We investigate here the influence of the dielectric constant on the plasma properties inside catalyst pores and in the sheath in front of the pores, for various pore sizes. The calculations are performed by a two-dimensional fluid model for an atmospheric pressure dielectric barrier discharge in helium. The electron impact ionization rate, electron temperature, electron and ion density, as well as the potential distribution and surface charge density, are analyzed for a better understanding of the discharge behavior inside catalyst pores. The results indicate that, in a 100 μm pore, the electron impact ionization in the pore, which is characteristic for the plasma generation inside the pore, is greatly enhanced for dielectric constants below 300. Smaller pore sizes only yield enhanced ionization for smaller dielectric constants, i.e., up to εr = 200, 150, and 50 for pore sizes of 50, 30, and 10 μm. Thus, the most common catalyst supports, i.e., Al2O3 and SiO2, which have dielectric constants around εr = 8−11 and 4.2, respectively, should allow more easily that microdischarges can be formed inside catalyst pores, even for smaller pore sizes. On the other hand, ferroelectric materials with dielectric constants above 300 never seem to yield plasma enhancement inside catalyst pores, not even for 100 μm pore sizes. Furthermore, it is clear that the dielectric constant of the material has a large effect on the extent of plasma enhancement inside the catalyst pores, especially in the range between εr = 4 and εr = 200. The obtained results are explained in detail based on the surface charge density at the pore walls,
and the potential distribution and electron temperature inside and above the pores. The results obtained with this model are
important for plasma catalysis, as the production plasma species in catalyst pores might affect the catalyst properties, and thus
improve the applications of plasma catalysis.

Abstract: In recent years, the ability of cold atmospheric pressure plasmas (CAPS) to selectively induce cell death in cancer cells has been widely established. This selectivity has been assigned to the reactive oxygen and nitrogen species (RONS) created in CAPs. To provide new insights in the search for an explanation
for the observed selectivity, we calculate the transfer free energy of multiple ROS across membranes containing a varying amount of cholesterol. The cholesterol fraction is investigated as a selectivity parameter because membranes of cancer cells are known to contain lower fractions of cholesterol compared to healthy cells. We find that cholesterol has a significant effect on the permeation of
reactive species across a membrane. Indeed, depending on the specific reactive species, an increasing cholesterol fraction can lead to (i) an increase of the transfer free energy barrier height and width, (ii) the formation of a local free energy minimum in the center of the membrane and (iii) the creation of extra free energy barriers due to the bulky sterol rings. In the context of plasma oncology, these observations suggest that the increased ingress of RONS in cancer cells can be explained by the decreased cholesterol fraction of their cell membrane.

Abstract: Plasma-surface interactions are in general highly complex due to the interplay of many concurrent processes. Molecular dynamics simulations provide insight in some of these processes, subject to the accessible time and length scales, and the availability of suitable force fields. In this introductory tutorial-style review, we aim to describe the current capabilities and limitations of molecular dynamics simulations in this field, restricting ourselves to low-temperature nonthermal plasmas. Attention is paid to the simulation of the various fundamental processes occurring, including sputtering, etching, implantation, and deposition, as well as to what extent the basic plasma components can be accounted for, including ground state and excited species, electric fields, ions, photons, and electrons. A number of examples is provided, giving an bird’s eye overview of the current state of the field.

Abstract: In this work, the electrostatic control in metal-oxide-semiconductor field-effect transistors based on MoS2 is studied, with respect to the number of MoS2 layers in the channel and to the equivalent oxide thickness of the gate dielectric, using first-principles calculations combined with a quantum transport formalism. Our simulations show that a compromise exists between the drive current and the electrostatic control on the channel. When increasing the number of MoS2 layers, a degradation of the device performances in terms of subthreshold swing and OFF currents arises due to the screening of the MoS2 layers constituting the transistor channel. Published by AIP Publishing.

Abstract: CO2 conversion into value-added chemicals is gaining increasing interest in recent years, and a gliding arc plasma has great potential for this purpose, because of its high energy efficiency. In this study, a chemical reaction kinetics model is presented to study the CO2 splitting in a gliding arc discharge. The calculated
conversion and energy efficiency are in good agreement with experimental data in a range of different operating conditions. Therefore, this reaction kinetics model can be used to elucidate the dominant chemical reactions contributing to CO2 destruction and formation. Based on this reaction pathway analysis, the restricting factors for CO2 conversion are figured out, i.e., the reverse reactions and the small treated gas fraction. This allows us to propose some solutions in order to improve the CO2 conversion, such as decreasing the gas temperature, by using a high frequency discharge, or increasing the power
density, by using a micro-scale gliding arc reactor, or by removing the reverse reactions, which could be realized in practice by adding possible scavengers for O atoms, such as CH4. Finally, we compare our results with other types of plasmas in terms of conversion and energy efficiency, and the results illustrate that gliding arc discharges are indeed quite promising for CO2 conversion, certainly when keeping in mind the possible solutions for further performance improvement.

Abstract: The discharge properties of the plasma bulk flow near the surface of heated compound-materials strongly affects the kinetic layer parameters modeled and manifested in the Knudsen layer. This paper extends the widely used two-layer kinetic ablation model to the ablation controlled non-equilibrium discharge due to the fact that the local thermodynamic equilibrium (LTE) approximation is often violated as a result of the interaction between the plasma and solid walls. Modifications to the governing set of equations, to account for this effect, are derived and presented by assuming that the temperature of the electrons deviates from that of the heavy particles. The ablation characteristics of one typical material, polytetrafluoroethylene (PTFE) are calculated with this improved model. The internal degrees of freedom as well as the average particle mass and specific heat ratio of the polyatomic vapor, which strongly depends on the temperature, pressure and plasma non-equilibrium degree and plays a crucial role in the accurate determination of the ablation behavior by this model, are also taken into account. Our assessment showed the significance of including such modifications related to the non-equilibrium effect in the study of vaporization of heated compound materials in ablation controlled arcs. Additionally, a two-temperature magneto-hydrodynamic (MHD) model accounting for the thermal non-equilibrium occurring near the wall surface is developed and applied into an ablation-dominated discharge for an electro-thermal chemical launch device. Special attention is paid to the interaction between the non-equilibrium plasma and the solid propellant surface. Both the mass exchange process caused by the wall ablation and plasma species deposition as well as the associated momentum and energy exchange processes are taken into account. A detailed comparison of the results of the non-equilibrium model with those of an equilibrium model is presented. The non-equilibrium results show a non-equilibrium region near the plasma-wall interaction region and this indicates the need for the consideration of the influence of the possible departure from LTE in the plasma bulk on the determination of ablation rate.

Abstract: Here, we present routes to increase CO2 conversion into CO using an atmospheric pressure dielectric-barrier discharge. The change in conversion as a function of simple plasma parameters, such as power, flow rate, but also frequency, on-and-off power pulse, thickness and the chemical nature of the dielectric, wall and gas temperature, are described. By means of an in-depth electrical characterization of the discharge (effective plasma voltage, dielectric voltage, plasma current, number and lifetime of the microdischarges), combined with infrared analysis of the walls of the reactor, optical emission spectroscopy for the gas temperature, and mass spectrometry for the CO2 conversion, we propose a global interpretation of the effect of all the experimental parameters on the conversion and efficiency of the reaction.