Abstract: Plasma catalysis is promising for various environmental, energy and chemical synthesis applications, but the underlying mechanisms are far from understood. Modeling can help to obtain a better insight in these mechanisms. Some burning questions relate to the plasma behavior inside packed bed reactors and whether plasma can penetrate into catalyst pores. In this paper, we try to provide answers to these questions, by means of both fluid modeling and particle-in-cell/Monte Carlo collision simulations. We present a short overview of recent findings obtained in our group by means of modeling, i.e., the enhanced electric field near the contact points and the streamer propagation through the packing in packed bed reactors, as well as the plasma behavior in catalyst pores, to determine the minimum pore size in which plasma streamers can penetrate.

Abstract: The electrical asymmetry effect is realized by applying multiple frequency power sources
(13.56 MHz and 27.12 MHz) to a capacitively biased substrate electrode in a specific inductively
coupled plasma reactor. On the one hand, by adjusting the phase angle θ between the multiple
frequency power sources, an almost linear self-bias develops on the substrate electrode, and
consequently the ion energy can be well modulated, while the ion flux stays constant within a
large range of θ. On the other hand, the plasma density and ion flux can be significantly
modulated by tuning the inductive power supply, while only inducing a small change in the self-
bias. Independent control of self-bias/ion energy and ion flux can thus be realized in this specific
inductively coupled plasma reactor.

Abstract: We investigated the mode transition from volume to surface discharge in a packed bed dielectric barrier discharge reactor by a two-dimensional particle-in-cell/Monte Carlo collision method. The calculations are performed at atmospheric pressure for various driving voltages and for gas mixtures with different N2 and O2 compositions. Our results reveal that both a change of the driving voltage and gas mixture can induce mode transition. Upon increasing voltage, a mode transition from hybrid (volume+surface) discharge to pure surface discharge occurs, because the charged species can escape much more easily to the beads and charge the bead surface due to the strong electric field at high driving voltage. This significant surface charging will further enhance the tangential component of the electric field along the dielectric bead surface, yielding surface ionization waves (SIWs). The SIWs will give rise to a high concentration of reactive species on the surface, and thus possibly enhance the surface activity of the beads, which might be of interest for plasma catalysis. Indeed, electron impact excitation and ionization mainly take place near the bead surface. In addition, the propagation speed of SIWs becomes faster with increasing N2 content in the gas mixture, and slower with increasing O2 content, due to the loss of electrons by attachment to O2
molecules. Indeed, the negative O-2 ion density produced by electron impact attachment is much higher than the electron and positive O+2 ion density. The different ionization rates between N2 and O2 gases will create different amounts of electrons and ions on the dielectric bead surface, which might also have effects in plasma catalysis.

Abstract: Although plasma catalysis is gaining increasing interest for various environmental applications, the underlying mechanisms are still far from understood. For instance, it is not yet clear whether and how plasma streamers can propagate in catalyst pores, and what is the minimum pore size to make this happen. As this is crucial information to ensure good plasma-catalyst interaction, we study here the mechanism of plasma streamer propagation in a catalyst pore, by means of a twodimensional particle-in-cell/Monte Carlo collision model, for various pore diameters in the nm range to μm-range. The so-called Debye length is an important criterion for plasma penetration into catalyst pores, i.e. a plasma streamer can penetrate into pores when their diameter is larger than the Debye length. The Debye length is typically in the order of a few 100 nm up to 1 μm at the conditions under study, depending on electron density and temperature in the plasma streamer. For pores in the range of ∼50 nm, plasma can thus only penetrate to some extent and at
very short times, i.e. at the beginning of a micro-discharge, before the actual plasma streamer reaches the catalyst surface and a sheath is formed in front of the surface. We can make plasma streamers penetrate into smaller pores (down to ca. 500 nm at the conditions under study) by increasing the applied voltage, which yields a higher plasma density, and thus reduces the Debye length. Our simulations also reveal that the plasma streamers induce surface charging of the catalyst pore sidewalls, causing discharge enhancement inside the pore, depending on pore diameter and depth.

Abstract: A controllable ferroelastic switching in ferroelectric/multiferroic oxides is highly desirable due to the non-volatile strain and possible coupling between lattice and other order parameter in heterostructures. However, a substrate clamping usually inhibits their elastic deformation in thin films without micro/nano-patterned structure so that the integration of the non-volatile strain with thin film devices is challenging. Here, we report that reversible in-plane elastic switching with a non-volatile strain of approximately 0.4% can be achieved in layered-perovskite Bi2WO6 thin films, where the ferroelectric polarization rotates by 90 degrees within four in-plane preferred orientations. Phase-field simulation indicates that the energy barrier of ferroelastic switching in orthorhombic Bi2WO6 film is ten times lower than the one in PbTiO3 films, revealing the origin of the switching with negligible substrate constraint. The reversible control of the in-plane strain in this layered-perovskite thin film demonstrates a new pathway to integrate mechanical deformation with nanoscale electronic and/or magnetoelectronic applications.