Abstract: A gliding arc (GA) plasma has great potential for CO2 conversion into value-added chemicals, because of its high energy efficiency. To improve the application, a 2D/3D fluid model is needed to investigate the CO2 conversion mechanisms in the actual discharge geometry. Therefore, the complex CO2 chemical kinetics description must be reduced due to the huge computational cost associated with 2D/3D models. This paper presents a chemistry reduction method for CO2 plasmas, based on the so-called directed relation graph method. Depending on the defined threshold values, some marginal species are identified. By means of a sensitivity analysis, we can further reduce the chemistry set by removing one by one the marginal species. Based on the socalled flux-sensitivity coupling, we obtain a reduced CO2 kinetics model, consisting of 36 or 15 species (depending on whether the 21 asymmetric mode vibrational states of CO2 are explicitly included or lumped into one group), which is applied to a GA discharge. The results are compared with those predicted with the full chemistry set, and very good agreement is reached. Moreover, the range of validity of the reduced CO2 chemistry set is checked, telling us that this reduced set is suitable for low power GA discharges. Finally, the time and spatial evolution of the CO2 plasma characteristics are presented, based on a 2D model with the reduced kinetics.

Abstract: This study focuses on a magnetically stabilized gliding arc (MGA) plasma. Two fully coupled flow-plasma models (in 3D and 2D) are presented. The 3D model is applied to compare the arc dynamics of the MGA with a traditional gas-driven gliding arc. The 2D model is used for a detailed parametric study on the effect of the external magnetic field. The results show that the relative velocity between the plasma and feed gas is generated due to the Lorentz force, which can increase the plasma-treated gas fraction. The magnetic field also helps to decrease the gas temperature by enhancing heat transfer and to increase the electron number density. This work shows the potential of an external magnetic field to control the gliding arc behavior, for enhanced gas conversion at low gas flow rates.

Abstract: We propose that low-strain Si/O superlattices can be constructed by connecting reconstructed Si{001} surfaces by Si-O-Si bridges. Ab initio calculations show that our models are energetically more favorable than all the models proposed so far. The part of our Si/O superlattice model is experimentally accessible just by oxidizing a Si( 001) substrate. To complete our Si/O superlattice model, we propose a three-step method. We also explore the potential of our Si/O superlattice models for new materials used in future Si electronics. We find that the location of the channel where the carriers travel can be controlled between the interfaces and the Si layers by the insertion of O atoms into the Si-Si dimers. By revealing the origins of the interface electron and hole states, we find that similar interface states should be easily achieved for Si slabs and Si substrates. Interestingly, the interface electrons and holes have small effective masses in the direction parallel to the channel and large effective masses in the direction normal to the channel, which makes the Si/O superlattices attractive to be used for channel materials. We also find that the valley splitting of Si is enhanced by the formation of the Si/O/Si interfaces, which is ideal for developing Si-based qubits. Our findings open new perspectives to design and control the electronic properties of Si.

Abstract: Monte Carlo (MC) methods have a long-standing history as partners of molecular dynamics (MD) to simulate the evolution of materials at the atomic scale. Among these techniques, the uniform-acceptance force-bias Monte Carlo (UFMC) method [ G. Dereli Mol. Simul. 8 351 (1992)] has recently attracted attention [ M. Timonova et al. Phys. Rev. B 81 144107 (2010)] thanks to its apparent capacity of being able to simulate physical processes in a reduced number of iterations compared to classical MD methods. The origin of this efficiency remains, however, unclear. In this work we derive a UFMC method starting from basic thermodynamic principles, which leads to an intuitive and unambiguous formalism. The approach includes a statistically relevant time step per Monte Carlo iteration, showing a significant speed-up compared to MD simulations. This time-stamped force-bias Monte Carlo (tfMC) formalism is tested on both simple one-dimensional and three-dimensional systems. Both test-cases give excellent results in agreement with analytical solutions and literature reports. The inclusion of a time scale, the simplicity of the method, and the enhancement of the time step compared to classical MD methods make this method very appealing for studying the dynamics of many-particle systems.

Abstract: Recently, core-shell Ag-Co, Ag-Cu, and “onionlike” Cu-Co equilibrium configurations were predicted in the case of isolated face centered cubic (fcc) bimetallic clusters, and three shell onionlike configurations were predicted in the case of ternary metallic clusters with spherical and truncated octahedral morphologies. In the present paper, immiscible binary CuCo and ternary AgCuCo clusters with icosahedral structures are studied as functions of their size and composition. Clusters studied are formed by 13, 55, 147, 309, and 561 atoms corresponding to the five smallest possible closed shell icosahedral structures. An embedded atom model potential is used to describe their cohesion. Equilibrium configurations are investigated by means of Metropolis Monte Carlo free energy minimization in the (NPT) canonical ensemble. Most simulations are achieved at 10 and 300 K. The effect of temperature on segregation ordering is systematically investigated. Selected cases are used to identify the effect of size and composition on melting. In contrast with fcc clusters, homogeneous onionlike configurations of binary clusters are not predicted. When it is allowed by the composition, a complete outer shell is formed by Cu in binary Cu-Co clusters and by Ag in ternary Ag-Cu-Co clusters. Depending on temperature, Co may precipitate into decahedral groups under the Cu vertices of the icosahedra in binary clusters, while the Co-Cu configuration in ternary clusters drastically depends on the Ag coating. Despite the multicomponent character of the clusters and the immiscibility of the species forming them, for most compositions and sizes, equilibrium structures remain close to perfectly icosahedral at 10 K as well as at 300 K.

Abstract: Amorphous semiconductors are usually characterized by a low charge carrier mobility, essentially related to their lack of long-range order. The development of such material with higher charge carrier mobility is hence challenging. Part of the issue comes from the difficulty encountered by first-principles simulations to evaluate concepts such as the electron effective mass for disordered systems since the absence of periodicity induced by the disorder precludes the use of common concepts derived from condensed matter physics. In this paper, we propose a methodology based on first-principles simulations that partially solves this problem, by quantifying the degree of delocalization of a wave function and of the connectivity between the atomic sites within this electronic state. We validate the robustness of the proposed formalism on crystalline and molecular systems and extend the insights gained to disordered/amorphous InGaZnO4 and Si. We also explore the properties of p-type oxide semiconductor candidates recently reported to have a low effective mass in their crystalline phases [G. Hautier et al., Nat. Commun. 4, 2292 (2013)]. Although in their amorphous phase none of the candidates present a valence band with delocalization properties matching those found in the conduction band of amorphous InGaZnO4, three of the seven analyzed materials show some potential. The most promising candidate, K2Sn2O3, is expected to possess in its amorphous phase a slightly higher hole mobility than the electron mobility in amorphous silicon.

Abstract: Nanosecond pulsed discharge plasma shows a high degree of non-equilibrium, and exhibits relatively high conversions in the dry reforming of methane. To further improve the application, a good insight of the underlying mechanisms is desired. We developed a chemical kinetics model to explore the underlying plasma chemistry in nanosecond pulsed discharge. We compared the calculated conversions and product selectivities with experimental results, and found reasonable agreement in a wide range of specific energy input. Hence, the chemical kinetics model is able to provide insight in the underlying plasma chemistry. The modeling results predict that the most important dissociation reaction of CO₂and CH₄is electron impact dissociation. C₂H₂is the most abundant hydrocarbon product, and it is mainly formed upon reaction of two CH₂radicals. Furthermore, the vibrational excitation levels of CO₂contribute for 85% to the total dissociation of CO₂.

Abstract: Vibrational excitation represents an efficient channel to drive the dissociation of CO₂in a non-thermal plasma. Its viability is investigated in low-pressure pulsed discharges, with the intention of selectively exciting the asymmetric stretching mode, leading to stepwise excitation up to the dissociation limit of the molecule. Gas heating is crucial for the attainability of this process, since the efficiency of vibration–translation (V–T) relaxation strongly depends on temperature, creating a feedback mechanism that can ultimately thermalize the discharge. Indeed, recent experiments demonstrated that the timeframe of V–T non-equilibrium is limited to a few milliseconds at ca. 6 mbar, and shrinks to the<italic>μ</italic>s-scale at 100 mbar. With the aim of backtracking the origin of gas heating in pure CO₂plasma, we perform a kinetic study to describe the energy transfers under typical non-thermal plasma conditions. The validation of our kinetic scheme with pulsed glow discharge experiments enables to depict the gas heating dynamics. In particular, we pinpoint the role of vibration–vibration–translation relaxation in redistributing the energy from asymmetric to symmetric levels of CO₂, and the importance of collisional quenching of CO₂electronic states in triggering the heating feedback mechanism in the sub-millisecond scale. This latter finding represents a novelty for the modelling of low-pressure pulsed discharges and we suggest that more attention should be paid to it in future studies. Additionally, O atoms convert vibrational energy into heat, speeding up the feedback loop. The efficiency of these heating pathways, even at relatively low gas temperature and pressure, underpins the lifetime of V–T non-equilibrium and suggests a redefinition of the optimal conditions to exploit the ‘ladder-climbing’ mechanism in CO₂discharges.

Abstract: Density functional theory is widely used to gain insights into molecule−metal surface reaction networks, which is important for a better understanding of catalysis. However, it is well-known that generalized gradient approximation (GGA)
density functionals (DFs), most often used for the study of reaction networks, struggle to correctly describe both gas-phase molecules and metal surfaces. Also, GGA DFs typically underestimate reaction barriers due to an underestimation of the selfinteraction energy. Screened hybrid GGA DFs have been shown to reduce this problem but are currently intractable for wide usage. In this work, we use a more affordable meta-GGA (mGGA) DF in combination with a nonlocal correlation DF for the first time to study and gain new insights into a catalytically important surface
reaction network, namely, CO2 hydrogenation on Cu. We show that the mGGA DF used, namely, rMS-RPBEl-rVV10, outperforms typical GGA DFs by providing similar or better predictions for metals and molecules, as well as molecule−metal surface adsorption
and activation energies. Hence, it is a better choice for constructing molecule−metal surface reaction networks.

Abstract: We report a plasma‐assisted CO₂hydrogenation to CH₃OH over Fe₂O₃/γ‐Al₂O₃catalysts, achieving 12% CO₂conversion and 58% CH₃OH selectivity at a temperature of nearly 80°C atm pressure. We investigated the effect of various supports and loadings of the Fe‐based catalysts, as well as optimized reaction conditions. We characterized catalysts by X‐ray powder diffraction (XRD), hydrogen temperature programmed reduction (H₂‐TPR), CO₂and CO temperature programmed desorption (CO₂/CO‐TPD), high‐resolution transmission electron microscopy (HRTEM), scanning transmission electron microscopy (STEM), x‐ray photoelectron spectroscopy (XPS), Mössbauer, and Fourier transform infrared<bold>(</bold>FTIR). The XPS results show that the enhanced CO₂conversion and CH₃OH selectivity are attributed to the chemisorbed oxygen species on Fe₂O₃/γ‐Al₂O₃. Furthermore, the diffuse reflectance infrared Fourier transform spectroscopy (DRIFTs) and TPD results illustrate that the catalysts with stronger CO₂adsorption capacity exhibit a higher reaction performance.<italic>In situ</italic>DRIFTS gain insight into the specific reaction pathways in the CO₂/H₂plasma. This study reveals the role of chemisorbed oxygen species as a key intermediate, and inspires to design highly efficient catalysts and expand the catalytic systems for CO₂hydrogenation to CH₃OH.

Abstract: Plasma is gaining interest for CH4 conversion into higher hydrocarbons and H2. However, the performance in terms of conversion and selectivity toward different hydrocarbons is different for different plasma types, and the underlying mechanisms are not yet fully understood. Therefore, we study here these mechanisms in different plasma sources, by means of a chemical kinetics model. The model is first validated by comparing the calculated conversions and hydrocarbon/H2 selectivities with experimental results in these different plasma types and over a wide range of specific energy input (SEI) values. Our model predicts that vibrational−translational nonequilibrium is negligible in all CH4 plasmas investigated, and instead, thermal conversion is important. Higher gas temperatures also lead to a more selective production of unsaturated hydrocarbons (mainly C2H2) due to neutral dissociation of CH4 and subsequent dehydrogenation processes, while three-body recombination reactions into saturated hydrocarbons (mainly C2H6, but also higher hydrocarbons) are dominant in low temperature plasmas.

Abstract: Plasma technology is gaining increasing interest for CO2 conversion. The gas temperature in (and after) the plasma reactor largely affects the performance. Therefore, we examine the effect of cooling/quenching, during and after the plasma, on the CO2 conversion and energy efficiency, for typical “warm” plasmas, by means of chemical kinetics modeling. For plasmas at low specific energy input (SEI ∼ 0.5 eV/molecule), it is best to quench at the plasma end, while for high-SEI plasmas (SEI ∼ 4 eV/molecule), quenching at maximum conversion is better. For low-SEI plasmas, quenching can even increase the conversion beyond the dissociation in the plasma, known as superideal quenching. To better understand the effects of quenching at different plasma conditions, we study the dissociation and recombination rates, as well as the vibrational distribution functions (VDFs) of CO2, CO, and O2. When a high vibrational−translational (VT) nonequilibrium exists at the moment of quenching, the dissociation and recombination reaction rates both increase. Depending on the conversion degree at the moment of quenching, this can lead to a net increase or decrease of CO2 conversion. In general, however, and certainly for equilibrium plasmas at high temperature, quenching after the plasma helps prevent recombination reactions and clearly enhances the final CO2 conversion. We also investigate the effect of different quenching cooling rates on the CO2 conversion and energy efficiency. Finally, we compare plasma-based conversion to purely thermal conversion. For warm plasmas with typical temperatures of 3000−4000 K, the conversion is roughly thermal.

Abstract: We developed a microkinetic model to reveal the effects of plasma-generated radicals, intermediates, and vibrationally excited species on the catalytic hydrogenation of CO2 to CH3OH on a Cu(111) surface. As a benchmark, we first present the mechanisms of thermal catalytic CH3OH formation. Our model predicts that the reverse water-gas shift reaction followed by CO hydrogenation, together with the formate path, mainly contribute to CH3OH formation in thermal catalysis. Adding plasma-generated radicals and intermediates results in a higher CH3OH turnover frequency (TOF) by six to seven orders of magnitude, showing the potential of plasma-catalytic CO2 hydrogenation into CH3OH, in accordance with the literature. In addition, CO2 vibrational excitation further increases the CH3OH TOF, but the effect is limited due to relatively low vibrational temperatures under typical plasma catalysis conditions. The predicted increase in CH3OH formation by plasma catalysis is mainly attributed to the increased importance of the formate path. In addition, the conversion of plasma-generated CO to HCO* and subsequent HCOO* or H2CO* formation contribute to CH3OH formation. Both pathways bypass the HCOO* formation from CO2, which is the main bottleneck in the process. Hence, our model points toward the important role of CO, but also O, OH, and H radicals, as they influence the reactions that consume CO2 and CO. In addition, our model reveals that the H pressure should not be smaller than ca. half of the O pressure in the plasma as this would cause O* poisoning, which would result in very small product TOFs. Thus, plasma conditions should be targeted with a high CO and H content as this is favorable for CH3OH formation, while the O content should be minimized.

Abstract: In this paper we study the cold atmospheric pressure plasma jet, called kinpen, operating in Ar with different admixture fractions up to 1% pure , and + . Moreover, the device is operating with a gas curtain of dry air. The absolute net production rates of the biologically active ozone () and nitrogen dioxide () species are measured in the far effluent by quantum cascade laser absorption spectroscopy in the mid-infrared. Additionally, a zero-dimensional semi-empirical reaction kinetics model is used to calculate the net production rates of these reactive molecules, which are compared to the experimental data. The latter model is applied throughout the entire plasma jet, starting already within the device itself. Very good qualitative and even quantitative agreement between the calculated and measured data is demonstrated. The numerical model thus yields very useful information about the chemical pathways of both the and the generation. It is shown that the production of these species can be manipulated by up to one order of magnitude by varying the amount of admixture or the admixture type, since this affects the electron kinetics significantly at these low concentration levels.