Abstract: One of the most challenging and recurring problems when modeling plasmas is the lack of data on the key atomic and molecular reactions that drive plasma processes. Even when there are data for some reactions, complete and validated datasets of chemistries are rarely available. This hinders research on plasma processes and curbs development of industrial applications. The QDB project aims to address this problem by providing a platform for provision, exchange, and validation of chemistry datasets. A new data model developed for QDB is presented. QDB collates published data on both electron scattering and heavy-particle reactions. These data are formed into reaction sets, which are then validated against experimental data where possible. This process produces both complete chemistry sets and identifies key reactions that are currently unreported in the literature. Gaps in the datasets can be filled using established theoretical methods. Initial validated chemistry sets for SF6/CF4/O2 and SF6/CF4/N2/H2 are presented as examples.

Abstract: The formation process of a microdischarge (MD) in both μm- and nm-sized catalyst pores is simulated by a two-dimensional particle-in-cell/Monte Carlo collision model. A parallel-plate dielectric barrier discharge configuration in filamentary mode is considered in ambient air. The discharge is powered by a high voltage pulse. Our calculations reveal that a streamer can penetrate into the surface features of a porous catalyst and MDs can be formed inside both μm- and nm-sized pores, yielding ionization inside the pore. For the μm-sized pores, the ionization mainly occurs inside the pore, while for the nm-sized pores the ionization is strongest near and inside the pore. Thus, enhanced discharges near and inside the mesoporous catalyst are observed. Indeed, the maximum values of the electric field, ionization rate and electron density occur near and inside the pore. The maximum electric field and electron density inside the pore first increase when the pore size rises from 4 nm to 10 nm, and then they decrease for the 100 nm pore, due to
a more pronounced surface discharge for the smaller pores. However, the ionization rate is highest for the 100 nm pore due to the largest effective ionization region.

Abstract: Packed bed plasma reactors (PBPRs) are gaining increasing interest for use in environmental applications, such as greenhouse gas conversion into value-added chemicals or renewable fuels and volatile pollutant removal (e.g. NOx, VOC, K), as they enhance the conversion and energy efficiency of the process compared to a non-packed reactor. However, the plasma behaviour in a PBPR is not well understood. In this paper we demonstrate, by means of a fluid model, that the discharge behaviour changes considerably when changing the size of the packing beads and their dielectric constant, while keeping the interelectrode spacing constant. At low dielectric constant, the plasma is spread out over the full discharge gap, showing significant density in the voids as well as in the connecting void channels. The electric current profile shows a strong peak during each half cycle. When the dielectric constant increases, the plasma becomes localised in the voids, with a current profile consisting of many smaller peaks during each half cycle. For large bead sizes, the shift from full gap discharge to localised discharges takes place at a higher dielectric constant than for smaller beads. Furthermore, smaller beads or beads with a lower dielectric constant require a higher breakdown voltage to cause plasma formation.

Abstract: A database-driven approach combined with ab initio density functional theory (DFT) simulations is used to identify and simulate alternative ferroelectric materials beyond Hf(Zr)O-2. The database-driven screening method identifies a class of wurtzite ferroelectric materials. DFT simulations of wurtzite magnesium chalcogenides, including MgS, MgSe, and MgTe, show their potential to achieve improved ferroelectric (FE) stability, simple atomistic unit cell structure, and large FE polarization. Strain engineering can effectively modulate the FE switching barrier height for facilitating FE switching. The effect of the piezoelectric property on the FE switching barrier heights is also examined.

Abstract: Plasma catalysis continues to gain increasing scientific interest, both in established fields like toxic waste abatement and emerging fields like greenhouse gas conversion into value-added chemicals. Attention is typically focused on the obtained conversion process selectivity, rates and energy efficiency. Much less attention is usually paid to the underlying mechanistic aspects of the processes that occur. In this contribution, we critically examine a number of fundamentally important nanoscale thermodynamic aspects of plasma catalysis, which are very relevant to these processes but so far have been overlooked or insufficiently covered in the plasma catalysis literature.

Abstract: A packed bed dielectric barrier discharge plasma reactor is computationally studied with a fluid model. Two different complementary axisymmetric 2D geometries are used to mimic the intrinsic 3D problem. It is found that a packing enhances the electric field strength and electron temperature at the contact points of the dielectric material due to polarization of the beads by the applied potential. As a result, these contact points prove to be of direct importance to initiate the plasma. At low applied potential, the discharge stays at the contact points, and shows the properties of a Townsend discharge. When a high enough potential is applied, the plasma will be able to travel through the gaps in between the beads from wall to wall, forming a kind of glow discharge. Therefore, the inclusion of a so-called ‘channel of voids’ is indispensable in any type of packed bed modelling.

Abstract: In this experimental study, a flowing dielectric barrier discharge operating at atmospheric pressure is used for the splitting of CO2 into O2 and CO. The influence of the applied frequency and plasma power on the microdischarge properties is investigated to understand their role on the CO2 conversion. Electrical measurements are carried out to explain the conversion trends and to characterize the microdischarges through their number, their lifetime,
their intensity and the induced electrical charge. Their influence on the gas and electrode temperatures is also evidenced through optical emission spectroscopy and infrared imaging. It is shown that, in our configuration, the conversion depends mostly on the charge delivered in the plasma and not on the effective plasma voltage when the applied power is modified. Similarly, at constant total current, a better conversion is observed at low frequencies, where a less filamentary discharge regime with a higher effective plasma voltage than that at a higher
frequency is obtained.

Abstract: Dielectric barrier discharges (DBDs) are commonly used to generate cold plasmas at
atmospheric pressure. Whatever their configuration (tubular or planar), the presence of a dielectric barrier is mandatory to prevent too much charge build up in the plasma and the formation of a thermal arc. In this article, the role of the barrier thickness (2.0, 2.4 and 2.8 mm) and of the kind of dielectric material (alumina, mullite, pyrex, quartz) is investigated on the filamentary behavior in the plasma and on the CO2 conversion in a tubular flowing DBD, by means of mass spectrometry measurements correlated with electrical characterization and IR imaging. Increasing the barrier thickness decreases the capacitance, while preserving the electrical charge. As a result, the voltage over the dielectric increases and a larger number of microdischarges is generated, which enhances the CO2 conversion. Furthermore, changing the dielectric material of the barrier, while keeping the same geometry and dimensions, also affects the CO2 conversion. The highest CO2 conversion and energy efficiency are obtained for quartz and alumina, thus not following the trend of the relative permittivity. From the
electrical characterization, we clearly demonstrate that the most important parameters are the somewhat higher effective plasma voltage (yielding a somewhat higher electric field and electron energy in the plasma) for quartz, as well as the higher plasma current (and thus larger electron density) and the larger number of microdischarge filaments (mainly for alumina, but also for quartz). The latter could be correlated to the higher surface roughness for alumina and to the higher voltage over the dielectric for quartz.

Abstract: CO2 conversion into value-added products has gained significant interest over the few last years, as the greenhouse gas concentrations constantly increase due to anthropogenic activities. Here we report on experiments for CO2 conversion by means of a cold atmospheric plasma using a cylindrical flowing dielectric barrier discharge (DBD) reactor. A detailed comparison of this DBD ignited in a so-called burst mode (i.e. where an AC voltage is applied during a limited amount of time) and pure AC mode is carried out to evaluate their effect on the conversion of CO2 as well as on the energy efficiency. Decreasing the duty cycle in the burst mode from 100% (i.e. corresponding to pure AC mode) to 40% leads to a rise in the
conversion from 16–26% and to a rise in the energy efficiency from 15 to 23%. Based on a detailed electrical analysis, we show that the conversion correlates with the features of the microfilaments. Moreover, the root-mean-square voltage in the burst mode remains constant as a function of the process time for the duty cycles <70%, while a higher duty cycle or the usual pure AC mode leads to a clear voltage decay by more than 500 V, over approximately 90 s, before reaching a steady state regime. The higher plasma voltage in the burst mode yields a higher electric field. This causes the increasing the electron energy, and therefore their
involvement in the CO2 dissociation process, which is an additional explanation for the higher CO2 conversion and energy efficiency in the burst mode.

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: In this paper, we show that the apparent delocalization of the conduction band reported from first-principles simulations for the high-mobility amorphous oxide semiconductor InGaZnO4 (a-IGZO) is an artifact induced by the periodic conditions imposed to the model. Given a sufficiently large unit-cell dimension (over 40 angstrom), the conduction band becomes localized. Such a model size is up to four times the size of commonly used models for the study of a-IGZO. This finding challenges the analyses done so far on the nature of the defects and on the interpretation of numerous electrical measurements. In particular, we re-interpret the meaning of the computed effective mass reported so far in literature. Our finding also applies to materials such as SiZnSnO, ZnSnO, InZnSnO, In2O3 or InAlZnO4 whose models have been reported to display a fully delocalized conduction band in the amorphous phase.

Abstract: In this study we report on a 2D fluid model of a gliding arc discharge in argon. Despite the 3D nature of the discharge, 2D models are found to be capable of providing very useful information about the operation of the discharge. We employ two modelsan axisymmetric and a Cartesian one. We show that for the considered experiment and the conditions of a low current arc (around 30 mA) in argon, there is no significant heating of the cathode surface and the discharge is sustained by field electron emission from the cathode accompanied by the formation of a cathode spot. The obtained discharge power and voltage are relatively sensitive to the surface properties and particularly to the surface roughness, causing effectively an amplification of the normal electric field. The arc body and anode region are not influenced by this and depend mainly on the current value. The gliding of the arc is modelled by means of a 2D Cartesian model. The arcelectrode contact points are analysed and the gliding mechanism along the electrode surface is discussed. Following experimental observations, the cathode spot is simulated as jumping from one point to another. A complete arc cycle is modelled from initial ignition to arc decay. The results show that there is no interaction between the successive gliding arcs.

Abstract: The chemical complexity of non-equilibrium plasmas poses a challenge for plasma modeling because of the computational load. This paper presents a dimension reduction method for such chemically complex plasmas based on principal component analysis (PCA). PCA is used to identify a low-dimensional manifold in chemical state space that is described by a small number of parameters: the principal components. Reduction is obtained since continuity equations only need to be solved for these principal components and not for all the species. Application of the presented method to a CO2 plasma model including state-to-state vibrational kinetics of CO2 and CO demonstrates the potential of the PCA method for dimension reduction. A manifold described by only two principal components is able to predict the CO2 to CO conversion at varying ionization degrees very accurately.

Abstract: An inductively coupled plasma (ICP) is analyzed by means of experiments and numerical simulation. Important plasma properties are analyzed, namely, the effective temperature inside the central channel and the mean flow velocity inside the plasma. Furthermore, the effect of torches with different injector diameters is studied by the model. The temperature inside the central channel is determined from the end-on collected line-to-background ratio in dependence of the injector gas flow rates. Within the limits of 3% deviation, the results of the simulation and the experiments are in good agreement in the range of flow rates relevant for the analysis of relatively large droplets, i.e., 50 μm. The deviation increases for higher gas flow rates but stays below 6% for all flow rates studied. The velocity of the gas inside the coil region was determined by side-on analyte emission measurements with single monodisperse droplet introduction and by the analysis of the injector gas path lines in the simulation. In the downstream region significantly higher velocities were found than in the upstream region in both the simulation and the experiment. The quantitative values show good agreement in the downstream region. In the upstream region, deviations were found in the absolute values which can be attributed to the flow conditions in that region and because the methods used for velocity determination are not fully consistent. Eddy structures are found in the simulated flow lines. These affect strongly the way taken by the path lines of the injector gas and they can explain the very long analytical signals found in the experiments at low flow rates. Simulations were performed for different injector diameters in order to find conditions where good analyte transport and optimum signals can be expected. The results clearly show the existence of a transition flow rate which marks the lower limit for effective analyte transport conditions through the plasma. A rule-of-thumb equation was extracted from the results from which the transition flow rate can be estimated for different injector diameters and different injector gas compositions.

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.