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: Low-temperature plasma chemical kinetic models are particularly important to the plasma community. These models typically require dozens of inputs, especially rate coefficients. The latter are not always precisely known and it is not surprising that the error on the rate coefficient data can propagate to the model output. In this paper, we present a model that uses N = 400 different combinations of rate coefficients based on the uncertainty attributed to each rate coefficient, giving a good estimation of the uncertainty on the model output due to the rate coefficients. We demonstrate that the uncertainty varies a lot with the conditions and the type of output. Relatively low uncertainties (about 15%) are found for electron density and temperature, while the uncertainty can reach more than an order of magnitude for the population of the vibrational levels in some cases and it can rise up to 100% for the CO2 conversion. The reactions that are mostly responsible for the largest uncertainties are identified. We show that the conditions of pressure, gas temperature and power density have a great effect on the uncertainty and on which reactions lead to this uncertainty. In all the cases tested here, while the absolute values may suffer from large uncertainties, the trends observed in previous modeling work are still valid. Finally, in accordance with the work of Turner, a number of ‘good practices’ is recommended.

Abstract: Understanding the nature and effect of the multitude of plasma–surface interactions in plasma catalysis is a crucial requirement for further process development and improvement. A particularly intriguing and rather unique property of a plasma-catalytic setup is the ability of the plasma to modify the electronic structure, and hence chemical properties, of the catalyst through charging, i.e. the absorption of excess electrons. In this work, we develop a quantum chemical model based on density functional theory to study excess negative surface charges in a heterogeneous catalyst exposed to a plasma. This method is specifically applied to investigate plasma-catalytic CO2 activation on supported M/Al2O3 (M=Ti, Ni, Cu) single atom catalysts. We find that (1) the presence of a negative surface charge dramatically improves the reductive power of the catalyst, strongly promoting the splitting of CO2 to CO and oxygen, and (2) the relative activity of the investigated transition metals is also changed upon charging, suggesting that controlled surface charging is a powerful additional parameter to tune catalyst activity and selectivity. These results strongly point to plasma-induced surface charging of the catalyst as an important factor contributing to the plasma-catalyst synergistic effects frequently reported for plasma catalysis.

Abstract: Plasma catalysis is gaining increasing interest, but the underlying mechanisms are far from understood. Different catalyst materials will have different chemical effects, but in addition, they might also have different dielectric constants, which will affect surface charging, and thus the plasma behavior. In this work, we demonstrate that surface charging plays an important role in the streamer propagation and discharge enhancement inside catalyst pores, and in the plasma distribution along the dielectric surface, and this role greatly depends on the dielectric constant of the material. For εr50, surface charging causes the plasma to spread along the dielectric surface and inside the pores, leading to deeper plasma streamer penetration, while for εr>50 or for metallic coatings, the discharge is more localized, due to very weak surface charging. In addition, at εr=50, the significant surface charge density near the pore entrance causes a large potential drop at the sharp pore edges, which induces a strong electric field and results in most pronounced plasma enhancement near the pore entrance.

Abstract: In this work we study energy transport in a gliding arc discharge with two diverging flat
electrodes in argon gas at atmospheric pressure. The discharge is ignited at the shortest electrode
gap and it is pushed downstream by a forced gas flow. The current values considered are
relatively low and therefore a non-equilibrium plasma is produced. We consider two cases, i.e.
with high and low discharge current (28 mA and 2.8mA), and a constant gas flow of 10 lmin −1 ,
with a significant turbulent component to the velocity. The study presents an analysis of the
various energy transport mechanisms responsible for the redistribution of Joule heating to the
plasma species and the moving background gas. The objective of this work is to provide a
general understanding of the role of the different energy transport mechanisms in arc formation
and sustainment, which can be used to improve existing or new discharge designs. The work is
based on a three-dimensional numerical model, combining a fluid plasma model, the shear stress
transport Reynolds averaged Navier–Stokes turbulent gas flow model, and a model for gas
thermal balance. The obtained results show that at higher current the discharge is constricted
within a thin plasma column several hundred kelvin above room temperature, while in the low-
current discharge the combination of intense convective cooling and low Joule heating prevents
discharge contraction and the plasma column evolves to a static non-moving diffusive plasma,
continuously cooled by the flowing gas. As a result, the energy transport in the two cases is
determined by different mechanisms. At higher current and a constricted plasma column, the
plasma column is cooled mainly by turbulent transport, while at low current and an unconstricted
plasma, the major cooling mechanism is energy transport due to non-turbulent gas convection. In
general, the study also demonstrates the importance of turbulent energy transport in
redistributing the Joule heating in the arc and its significant role in arc cooling and the formation
of the gas temperature profile. In general, the turbulent energy transport lowers the average gas
temperature in the arc, thus allowing additional control of thermal non-equilibrium in the
discharge.

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: We studied the evolution of an isolated pulsed plasma in a vortex flow stabilised microwave (MW) discharge in N2 at 25 mbar via the combination of 0D kinetics modelling, iCCD imaging and laser scattering diagnostics. Quenching of electronically excited N2 results in fast gas heating and the onset of a thermal-ionisation instability, contracting the discharge volume. The onset of a thermal-ionisation instability driven by vibrational excitation pathways is found to facilitate significantly higher N2 conversion (i.e. dissociation to atomic N2 ) compared to pre-instability conditions, emphasizing the potential utility of this dynamic in future fixation applications. The instability onset is found to be instigated by super-elastic heating of the electron energy distribution tail via vibrationally excited N2 . Radial contraction of the discharge to the skin depth is found to occur post instability, while the axial elongation is found to be temporarily contracted during the thermal instability onset. An increase in power reflection during the thermal instability onset eventually limits the destabilising effects of exothermic electronically excited N2 quenching. Translational and vibrational temperature reach a quasi-non-equilibrium after the discharge contraction, with translational temperatures reaching ∼1200 K at the pulse end, while vibrational temperatures are found in near equilibrium with the electron energy (1 eV, or ∼11 600 K). This first description of the importance of electronically excited N2 quenching in thermal instabilities gives an additional fundamental understanding of N2 plasma behaviour in pulsed MW context, and thereby brings the eventual implementation of this novel N2 fixation method one step closer.

Abstract: Nitrogen oxyanions and oxyacids are important agents in atmospheric chemistry and medical biology. Although their chemical behavior in solution is relatively well understood, they may behave very differently at the water/air interface of atmospheric aerosols or at the membrane/water interface of cells. Here, we developed a fully classical model for molecular dynamics simulations of NO3−, NO2−, HNO3, and HNO2 in the framework of the GROMOS 53A6 and 54A7 force field versions. The model successfully accounted for the poorly structured solvation shell and ion pairing tendency of NO3−. Accurate pure-liquid properties and hydration free energies were obtained for the oxyacids. Simulations at the water/air interface showed a local enrichment of HNO3 and depletion of NO3−. The effect was discussed in light of earlier spectroscopic data and ab initio calculations, suggesting that HNO3 behaves as a weaker acid at the surface of water. Our model will hopefully allow for efficient and accurate simulations of nitrogen oxyanions and oxyacids in solution and at microheterogeneous interface environments.

Abstract: In the present work, the effect of various embedded atom impurities on tuning electronic and magnetic properties of C3N4 and C4N3 nanosheets have been studied using first-principles calculations. Our calculations show that C3N4 is a semiconductor and it exhibits extraordinary electronic properties such as dilute-magnetic semiconductor (with H, F, Cl, Be, V, Fe and Co); metal (with N, P, Mg and Ca), half-metal (with Li, Na, K, Al, Sc, Cr, Mn, and Cu) and semiconductor (with O, S, B, C, Si, Ti, Ni and Zn) with the band gaps in the range of 0.3-2.0 eV depending on the species of embedded atom. The calculated electronic properties reveal that C4N3 is a half-metal and it retains half-metallic character with embedded H, O, S, F, B, N, P, Be, Mg, Al, Sc, V, Fe, Ni and Zn atoms. The substitution of Cl, C, Cr and Mn atoms create ferromagnetic-metal character in the C4N3 nanosheet, embedded Co and Cu atoms exhibit a dilute-magnetic semiconductor nature, and embedded Ti atoms result in the system becoming a semiconductor. Therefore, our results reveal the fact that the band gap and magnetism can be modified or induced by various atom impurities, thus, offering effective possibilities to tune the electronic and magnetic properties of C3N4 and C4N3 nanosheets.

Abstract: Nonmetal doping is an effective approach to modify the electronic band structure and enhance the photocatalytic performance of bismuth oxyhalides. Using density functional theory, we systematically examine the fundamental properties of single-layer BiOBr doped with boron (B) and phosphorus (P) atoms. The stability of the doped models is investigated based on the formation energies, where the substitutional doping is found to be energetically more stable under O-rich conditions than under Bi-rich ones. The results showed that substitutional doping of P atoms reduced the bandgap of pristine BiOBr to a greater extent than that of boron substitution. The calculation of the effective masses reveals that B doping can render the electrons and holes of pristine BiOBr lighter and heavier, respectively, resulting in a slower recombination rate of photoexcited electron-hole pairs. Based on the results of HOMO-LUMO calculations, the introduction of B atoms tends to increase the number of photocatalytically active sites. The top of the valence band and the conduction band bottom of the B doped BiOBr monolayer match well with the water redox potentials in an acidic environment. The absorption spectra propose that B(P) doping causes a red-shift. Overall, the results predict that nonmetal-doped BiOBr monolayers have a reduced bandgap, a slow recombination rate, more catalytically active sites, enhanced optical absorption edges, and reduced work functions, which will contribute to superior photocatalytic performance.

Abstract: In this work, we investigated systematically the electronic and optical properties of B doped C3N monolayers as well as B and N doped graphene bilayers (BN-Gr@2L). We found that the doping of B atoms leads to an enlarged band gap of the C3N monolayer and when the dopant concentration reaches 12.5%, an indirect-to-direct band gap switching occurs. In addition, with co-doping of B and N atoms on the graphene monolayer in the hexagonal configuration, an electronic transition from semi-metal to semiconductor occurs. Our optical results for B-C3N show a broad absorption spectrum in a wide visible range starting from 400 nm to 1000 nm with strong absorption intensity, making it a suitable candidate for nanoelectronic and optoelectronic applications. Interestingly, a transition from semi-metal to semiconductor emerges in the graphene monolayer with doping of B and N atoms. Furthermore, our results demonstrate that the in-plane strain and out-of-plane strain (pressure) can modulate the band gap of the BN-Gr@2L. The controllable electronic properties and optical features of the doped graphene bilayer by strain engineering may facilitate their practical performance for various applications in future.

Abstract: Topochemical reduction of La0.9Ba0.1MnO3 with NaH at 225 degrees C yields the brownmillerite phase La0.9Ba0.1MnO2.5. However, reduction with CaH2 at 435 degrees C results in the formation of La0.9Ba0.1Mn0.96O2.43 via the deintercalation of both oxide anions and manganese cations from the parent perovskite phase. Electron and neutron diffraction data reveal La0.9Ba0.1Mn0.96O2.43 adopts a complex noncentrosymmetric structure, described in space group I23, confirmed by SHG measurements. Low-temperature neutron diffraction data reveal La0.9Ba0.1Mn0.96O2.43 adopts an ordered magnetic structure in which all the nearest neighbor interactions are antiferromagnetic. However, the presence of ordered manganese cation-vacancies results in a net ferrimagnetic structure with net saturated moment of 0.157(2) mu B per manganese center.

Abstract: We conduct ballistic transport calculations to study the electronic transport properties of diazatriptycene molecule which can be self-assembled on metallic surfaces with uniform coverage and upright orientation of the functional head group. Due to its structural asymmetry, the molecule shows a clear current rectification, where the level of the rectification depends on the nature of the head group. For example, current rectification can be increased by more than a factor of 2 by anchoring the molecules to the electrode by CN functional group or introducing insulating CH2 group between the thiol end group and the adjacent phenyl ring. Our findings show the possibility of creating self-assembled monolayer of DAT molecules with controlled electronic transport properties through functionalization of the head group.

Abstract: Two-dimensional Transition Metal Dichalcogenides (TMDs) possessing extraordinary physical properties at reduced dimensionality have attracted interest due to their promise in electronic and optical device applications. However, TMD monolayers can show a broad range of different properties depending on their crystal phase; for example, H phases are usually semiconductors, while the T phases are metallic. Thus, controlling phase transitions has become critical for device applications. In this study, the energetically low-lying crystal structures of pristine and Janus TMDs are investigated by using ab initio Nudged Elastic Band and molecular dynamics simulations to provide a general explanation for their phase stability and transition properties. Across all materials investigated, the T phase is found to be the least stable and the H phase is the most stable except for WTe2, while the T' and T '' phases change places according to the TMD material. The transition energy barriers are found to be large enough to hint that even the higher energy phases are unlikely to undergo a phase transition to a more stable phase if they can be achieved except for the least stable T phase, which has zero barrier towards the T ' phase. Indeed, in molecular dynamics simulations the thermodynamically least stable T phase transformed into the T ' phase spontaneously while in general no other phase transition was observed up to 2100 K for the other three phases. Thus, the examined T ', T '' and H phases were shown to be mostly stable and do not readily transform into another phase. Furthermore, so-called mixed phase calculations considered in our study explain the experimentally observed lateral hybrid structures and point out that the coexistence of different phases is strongly stable against phase transitions. Indeed, stable complex structures such as metal-semiconductor-metal architectures, which have immense potential to be used in future device applications, are also possible based on our investigation.

Abstract: We report a high-throughput first-principles characterization of the structural, mechanical, electronic, and vibrational properties of tetragonal single-layer transition metal Xenes (t-TMXs). Our calculations revealed 22 dynamically, mechanically and chemically stable structures among the 96 possible free-standing layers present in the t-TMX family. As a fingerprint for their structural identification, we identified four characteristic Raman active phonon modes, namely three in-plane and one out-of-plane optical branches, with various intensities and frequencies depending on the material in question. Spin-polarized electronic calculations demonstrated that anti-ferromagnetic (AFM) metals, ferromagnetic (FM) metals, AFM semiconductors, and non-magnetic semiconductor materials exist within this family, evidencing the potential of t-TMXs for further use in multifunctional heterostructures.