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 paper, an O2 inductively coupled plasma used for plasma enhanced atomic layer deposition of Al2O3 thin films is investigated by means of modeling. This work intends to provide more information about basic plasma properties such as species densities and species fluxes to the substrate as a function of power and pressure, which might be hard to measure experimentally. For this purpose, a hybrid model developed by Kushner et al is applied to calculate the plasma characteristics in the reactor volume for different chamber pressures ranging from 1 to 10 mTorr and different coil powers ranging from 50 to 500 W. Density profiles of the various oxygen containing plasma species are reported as well as fluxes to the substrate under various operating conditions. Furthermore, different orientations of the substrate, which can be placed vertically or horizontally in the reactor, are taken into account. In addition, special attention is paid to the recombination process of atomic oxygen on the different reactor walls under the stated operating conditions. From this work it can be concluded that the plasma properties change significantly in different locations of the reactor. The plasma density near the cylindrical coil is high, while it is almost negligible in the neighborhood of the substrate. Ion and excited species fluxes to the substrate are found to be very low and negligible. Finally, the orientation of the substrate has a minor effect on the flux of O2, while it has a significant effect on the flux of O. In the horizontal configuration, the flux of atomic oxygen can be up to one order of magnitude lower than in the vertical configuration.

Abstract: This work examines the properties of a dielectric barrier discharge (DBD) reactor, built for CO2 decomposition, by means of electrical characterization, optical emission spectroscopy and gas chromatography. The discharge, formed in an electronegative gas (such as CO2, but also O2), exhibits clearly different electrical characteristics, depending on the surface conductivity of the reactor walls. An asymmetric current waveform is observed in the metaldielectric (MD) configuration, with sparse high-current pulses in the positive half-cycle (HC) and a more uniform regime in the negative HC. This indicates that the discharge is operating in two alternating regimes with rather different properties. At high CO2 conversion regimes, a conductive coating is deposited on the dielectric. This so-called coated MD configuration yields a symmetric current waveform, with current peaks in both the positive and negative HCs. In a double-dielectric (DD) configuration, the current waveform is also symmetric, but without current peaks in both the positive and negative HC. Finally, the DD configuration with conductive coating on the inner surface of the outer dielectric, i.e. so-called coated DD, yields again an asymmetric current waveform, with current peaks in the negative HC. These different electrical characteristics are related to the presence of the conductive coating on the dielectric wall of the reactor and can be explained by an increase of the local barrier capacitance available for charge transfer. The different discharge regimes affect the CO2 conversion, more specifically, the CO2 conversion is lowest in the clean DD configuration. It is somewhat higher in the coated DD configuration, and still higher in the MD configuration. The clean and coated MD configuration, however, gave similar CO2 conversion. These results indicate that the conductivity of the dielectric reactor walls can highly promote the development of the high-amplitude discharge current pulses and subsequently the CO2 conversion.

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: A two-dimensional fluid model, including the full set of Maxwell equations, has been developed and applied to investigate the effect of a phase shift between two power sources on the radial uniformity of several plasma characteristics in a hydrogen capacitively coupled plasma. This study was carried out at various frequencies in the range 13.56200 MHz. When the frequency is low, at 13.56 MHz, the plasma density is characterized by an off-axis peak when both power sources are in-phase (phgr = 0), and the best radial uniformity is obtained at phgr = π. This trend can be explained because the radial nonuniformity caused by the electrostatic edge effect can be effectively suppressed by the phase-shift effect at a phase difference equal to π. When the frequency rises to 60 MHz, the plasma density profiles shift smoothly from edge-peaked over uniform to centre-peaked as the phase difference increases, due to the pronounced standing-wave effect, and the best radial uniformity is reached at phgr = 0.3π. At a frequency of 100 MHz, a similar behaviour is observed, except that the maximum of the plasma density moves again towards the radial edge at the reverse-phase case (phgr = π), because of the dominant skin effect. When the frequency is 200 MHz, the bulk plasma density increases significantly with increasing phase-shift values, and a better uniformity is obtained at phgr = 0.4π. This is because the density in the centre increases faster than at the radial edge as the phase difference rises, due to the increasing power deposition Pz in the centre and the decreasing power density Pr at the radial edge. As the phase difference increases to π, the maximum near the radial edge becomes obvious again. This is because the skin effect has a predominant influence on the plasma density under this condition, resulting in a high density at the radial edge. Moreover, the axial ion flux increases monotonically with phase difference, and exhibits similar profiles to the plasma density. The calculation results illustrate that the radial uniformity of the various plasma characteristics is strongly dependent on the applied frequency and the phase shift between both power sources, which is important to realize, for controlling the uniformity of the plasma etch and deposition processes.

Abstract: Because of its unique properties, plasma technology has gained much prominence in the
microelectronics industry. Recently, environmental and energy applications of plasmas have gained a lot of attention. In this area, the focus is on converting CO 2 and reforming hydrocarbons, with the goal of developing an efficient single-step ‘gas-to-liquid’ (GTL) process. Here we show that applying tri-reforming principles to plasma—further called ‘plasma-based multi-reforming’—allows us to better control the plasma chemistry and thus the formed products. To demonstrate this, we used chemical kinetics calculations supported by experiments and reveal that better control of the plasma chemistry can be achieved by adding O 2 or H 2 O to a mixture containing CH 4 and CO 2 (diluted in N 2 ). Moreover, by adding O 2 and H 2 O simultaneously, we can tune the plasma chemistry even further, improving the conversions, thermal efficiency and methanol yield. Unlike thermocatalytic reforming, plasma-based reforming is capable of producing methanol in a single step; and compared with traditional plasma-based dry reforming, plasma-based multi-reforming increases the methanol yield by more than seven times and the thermal efficiency by 49%, as revealed by our model calculations. Thus, we believe that by using plasma-based multi-reforming, ‘gas-to-liquid’ conversion may be made efficient and scalable.

Abstract: The growth mechanism and chirality formation of a single-walled carbon nanotube (SWNT) on a surface-bound nickel nanocluster are investigated by hybrid reactive molecular dynamics/force-biased Monte Carlo simulations. The validity of the interatomic potential used, the so-called ReaxFF potential, for simulating catalytic SWNT growth is demonstrated. The SWNT growth process was found to be in agreement with previous studies and observed to proceed through a number of distinct steps, viz., the dissolution of carbon in the metallic particle, the surface segregation of carbon with the formation of aggregated carbon clusters on the surface, the formation of graphitic islands that grow into SWNT caps, and finally continued growth of the SWNT. Moreover, it is clearly illustrated in the present study that during the growth process, the carbon network is continuously restructured by a metal-mediated process, thereby healing many topological defects. It is also found that a cap can nucleate and disappear again, which was not observed in previous simulations. Encapsulation of the nanoparticle is observed to be prevented by the carbon network migrating as a whole over the cluster surface. Finally, for the first time, the chirality of the growing SWNT cap is observed to change from (11,0) over (9,3) to (7,7). It is demonstrated that this change in chirality is due to the metal-mediated restructuring process.

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: In this study, we investigated the diffusion of H-atoms to the subsurface and their further diffusion into the bulk of a Ni(111) crystal by means of density functional theory calculations in the context of thermal and plasma-assisted catalysis. The H-atoms at the surface can originate from the dissociative adsorption of H2 or CH4 molecules, determining the surface H-coverage. When a threshold H-coverage is passed, corresponding to 1.00 ML for the crystalline Ni(111) surface, the surface-bound H-atoms start to diffuse to the subsurface. A similar threshold coverage is observed for the interstitial H-coverage. Once the interstitial sites are filled up with a coverage above 1.00 ML of H, dissolution of interstitial H-atoms to the layer below the interstitial sites will be initiated. Hence, by applying a high pressure or inducing a reactive plasma and high temperature, increasing the H-flux to the surface, a large amount of hydrogen can diffuse in a crystalline metal like Ni and can be absorbed. The formation of metal hydride may modify the entire reaction kinetics of the system. Equivalently, the H-atoms in the bulk can easily go back to the surface and release a large amount of heat. In a plasma process, H-atoms are formed in the plasma, and therefore the energy barrier for dissociative adsorption is dismissed, thus allowing achievement of the threshold coverage without applying a high pressure as in a thermal process. As a result, depending on the crystal plane and type of metal, a large number of H-atoms can be dissolved (absorbed) in the metal catalyst, explaining the high efficiency of plasma-assisted catalytic reactions. Here, the mechanism of H-dissolution is established as a chemical identifier for the investigation of the reaction kinetics of a chemical process.

Abstract: This study reports on the possible effects of OH radical impact on the transmembrane domain 6 of P-glycoprotein, TM6, which plays a crucial role in drug binding in human cells. For the first time, we employ molecular dynamics (MD) simulations based on the self-consistent charge density functional tight binding (SCC-DFTB) method to elucidate the potential sites of fragmentation and mutation in this domain upon impact of OH radicals, and to obtain fundamental information about the underlying reaction mechanisms. Furthermore, we apply non-reactive MD simulations to investigate the long-term effect of this mutation, with possible implications for drug binding. Our simulations indicate that the interaction of OH radicals with TM6 might lead to the breaking of C-C and C-N peptide bonds, which eventually cause fragmentation of TM6. Moreover, according to our simulations, the OH radicals can yield mutation in the aromatic ring of phenylalanine in TM6, which in turn affects its structure. As TM6 plays an important role in the binding of a range of cytotoxic drugs with P-glycoprotein, any changes in its structure are likely to affect the response of the tumor cell in chemotherapy. This is crucial for cancer therapies based on reactive oxygen species, such as plasma treatment.

Abstract: The adsorption probability and reaction behavior of CHx plasma species on various nickel catalyst surfaces is investigated by means of reactive molecular dynamics (MD) simulations using the ReaxFF potential. Such catalysts are used in the reforming of hydrocarbons and in the growth of carbon nanotubes, and further insight in the underlying mechanisms of these processes is needed to increase their applicability. Single and consecutive impacts of CHx radicals (x={1,2,3}) were performed on four different Ni surfaces, at a temperature of 400 K. The adsorption probability is shown to be related to the number of free electrons, i.e. a higher number leads to more adsorptions, and the steric hindrance caused by the hydrogen atoms bonded to the impacting CHx species. Furthermore, some of the CH bonds break after adsorption, which generally leads to diffusion of the hydrogen atom over the surface. Additionally, these adsorbed H-atoms can be used in reactions to form new molecules, such as CH4 and C2Hx, although this is dependent on the precise morphology of the surface. New molecules are also formed by subtraction of H-atoms from adsorbed radicals, leading to occasional formation of H2 and C2Hx molecules.

Abstract: The adsorption, dissociation, and diffusion of CO2 on the anatase (001) surface was studied using DFT by means of the generalized gradient approximation using the Perdew−Burcke−Ernzerhof (PBE)-functional and applying corrections for long-range dispersion interactions. Different stable adsorption configurations were identified for the fully oxidized surface. The most stable adsorption configuration is the monodentated carbonate-like structure. Small energy barriers were identified for the conversion of a physisorbed to a chemisorbed configuration.
CO2 dissociation is found to be unfeasible on the stoichiometric surface. The introduction of oxygen vacancy defects gives rise to new highly stable adsorption configurations with a stronger activation of the C−O bonds. This leads to the possibility of exothermic dissociation of CO2 with barriers up to 22.2 kcal/mol,
corresponding to chemical lifetimes of less than 4 s at 300 K. These reactions cause a CO molecule to be formed, which will easily desorb, and the reduced surface to become oxidized. It is clear that oxygen vacancy defects play a key role in the catalytic activity of an anatase (001) surface. Oxygen vacancies play an important role in the dissociation of CO2 on the anatase (001) surface, and will play a significant role in complex problems, such as the catalytic conversion of CO2 to value-added chemicals.

Abstract: By means of chemical kinetics modeling, it is possible to elucidate the main dissociation mechanisms of CO2 in a gliding arc plasmatron (GAP). We obtain good agreement between the calculated and experimental conversions and energy efficiencies, indicating that the model can indeed be used to study the underlying mechanisms. The calculations predict that vibration-induced dissociation is the main dissociation mechanism of CO2, but it occurs mainly from the lowest vibrational levels because of fast thermalization of the vibrational distribution. Based on these findings, we propose ideas for improving the performance of the GAP, but testing of these ideas in the simulations reveals that they do not always lead to significant enhancement, because of other side effects, thus illustrating the complexity of the process. Nevertheless, the model allows more insight into the underlying mechanisms to be obtained and limitations to be identified.