Abstract: In this paper we analyse the gas phase production and loss pathways for several biomedically active species, i.e. N2(A), O, O3, O2(a), N, H, HO2, OH, NO, NO2, N2O5, H2O2, HNO2 and HNO3, in an argon plasma jet flowing into an open humid air atmosphere. For this purpose, we employ a zero-dimensional reaction kinetics model to mimic the typical experimental conditions by fitting several parameters to experimentally measured values. These include ambient air diffusion, the gas temperature profile and power deposition along the jet effluent. We focus in detail on how the pathways of the biomedically active species change as a function of the position in the effluent, i.e. inside the discharge device, active plasma jet effluent and afterglow region far from the nozzle. Moreover, we demonstrate how the reaction kinetics and species production are affected by different ambient air humidities, total deposited power into the plasma and gas temperature along the jet. It is shown that the dominant pathways can drastically change as a function of the distance from the nozzle exit or experimental conditions.

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: The aim of this work consists of the evaluation of atmospheric pressure dielectric barrier discharges for the conversion of greenhouse gases into useful compounds. Therefore, pure CO2 feed flows are administered to the discharge zone at varying discharge frequency, power input, gas temperature and feed flow rates, aiming at the formation of CO and O2. The discharge obtained in CO2 is characterized as a filamentary mode with a microdischarge zone in each half cycle of the applied voltage. It is shown that the most important parameter affecting the CO2-conversion levels is the gas flow rate. At low flow rates, both the conversion and the CO-yield are significantly higher. In addition, also an increase in the gas temperature and the power input give rise to higher conversion levels, although the effect on the CO-yield is limited. The optimum discharge frequency depends on the power input level and it cannot be unambiguously stated that higher frequencies give rise to increased conversion levels. A maximum CO2 conversion of 30% is achieved at a flow rate of 0.05 L min−1, a power density of 14.75 W cm−3 and a frequency of 60 kHz. The most energy efficient conversions are achieved at a flow rate of 0.2 L min−1, a power density of 11 W cm−3 and a discharge frequency of 30 kHz.

Abstract: We study vortex configurations in mesoscopic superconducting squares. Our theoretical approach is based on the analytical solution of the London equation using the Green's function method. The potential energy landscape found is then used in Langevin-type molecular-dynamics simulations to obtain stable vortex configurations. We show that the filling rules for vortices in squares with increasing applied magnetic field can be formulated, although in a different manner than in disks, in terms of the formation of vortex 'shells'. We discuss metastable states and the stability of the vortex configurations found with respect to variations of the material parameters and deformations of the shape of the sample.

Abstract: Nanostructured carbon materials show a tremendous variety in atomic structure, morphology, properties, and applications. As all properties are ultimately determined by the structure of the material, a thorough understanding of the growth mechanisms that give rise to the particular structure is critical. On many occasions, it has been shown that plasma enhanced growth can be strongly beneficial. This review will describe the authors current understanding of plasma enhanced growth of carbon nanotubes, the prototypical example of nanostructured carbon materials, as obtained from experiments, simulations, and modeling. Specific emphasis is put on where experiments and computational approaches correspond, and where they differ. Also, the current status on simulating PECVD growth of some other carbon nanomaterials is reviewed, including amorphous carbon, graphene, and metallofullerenes. Finally, computational challenges with respect to the simulation of PECVD growth are identified.

Abstract: A hitherto unknown composition is highlighted in the red and black inks preserved on ancient Egyptian papyri from the Roman period (circa 100 to 200 CE). Synchrotron-based macro-X-ray fluo-rescence (XRF) mapping brings to light the presence of iron (Fe) and lead (Pb) compounds in the majority of the red inks inscribed on 12 papyrus fragments from the Tebtunis temple library. The iron-based compounds in the inks can be assigned to ocher, notably due to the colocalization of Fe with aluminum, and the detection of hematite (Fe2O3) by micro-X-ray diffraction. Using the same techniques together with micro-Fourier transform infrared spectroscopy, Pb is shown to be associated with fatty acid phosphate, sulfate, chloride, and carboxylate ions. Moreover, microXRF maps reveal a peculiar distribution and colocalization of Pb, phosphorus (P), and sulfur (S), which are present at the micrometric scale resembling diffused “coffee rings” surrounding the ocher particles imbedded in the red letters, and at the submicrometric scale concentrated in the papyrus cell walls. A similar Pb, P, and S composition was found in three black inks, suggesting that the same lead components were employed in the manufacture of carbon-based inks. Bearing in mind that pigments such as red lead (Pb3O4) and lead white (hydrocerussite [Pb-3(CO3)(2)(OH)(2)] and/or cerussite [PbCO3]) were not detected, the results presented here suggest that the lead compound in the ink was used as a drier rather than as a pigment. Accordingly, the study calls for a reassessment of the composition of lead-based components in ancient Mediterranean pigments.

Abstract: Particulate matter (PM) deposition on urban green enables the collection of particulate pollution from a diversity of contexts, and insight into the physico-chemical profiles of PM is key for identifying main polluting sources. This study reports on the morphological and elemental characterization of PM2-10 deposited on ivy leaves from five different environments (forest, rural, roadside, train, industry) in the region of Antwerp, Belgium. Ca. 40,000 leaf-deposited particles were thoroughly investigated by particle-based analysis using scanning electron microscopy coupled with energy-dispersive X-ray spectroscopy (SEM/EDX) and their physico-chemical characteristics were explored for PM source apportionment purposes. The size distribution of all deposited particles was biased towards small-sized PM, with 32% of the particles smaller than 2.5 mu m (PM2.5) and median diameters of 2.80-3.09 mu m. The source type influenced both the particles' size and morphology (aspect ratio and shape), with roadside particles being overall the smallest in size and the most spherical. While forest and rural elemental profiles were associated with natural PM, the industry particles revealed the highest anthropogenic metal input. PM2-10 profiles for roadside and train sites were rather comparable and only distinguishable when evaluating the fine (2-2.5 mu m) and coarse (2.5-10 mu m) PM fractions separately, which enabled the identification of a larger contribution of combustion-derived particles (small, circular, Fe-enriched) at the roadside compared to the train. Random forest prediction model classified the source type correctly for 61-85% of the leaf-deposited PM. The still modest classification accuracy denotes the influence of regional background PM and demands for additional fingerprinting techniques to facilitate source apportionment. Nonetheless, the obtained results demonstrate the utility of leaf particle-based analysis to fingerprint and pinpoint source-specific PM, particularly when considering both the composition and size of leaf-deposited particles.

Abstract: In this work particle-in-cell/Monte Carlo collision simulations are performed to study the heating mechanism and plasma characteristics in direct current (dc) superposed radio-frequency (RF) capacitively coupled plasmas, operated both in single-frequency (SF) and dual-frequency (DF) regimes. An RF (60/2 MHz) source is applied on the bottom electrode to sustain the discharge, and a dc source is fixed on the top electrode. The heating mechanism appears to be very different in dc superposed SF and DF discharges. When only a single source of 60 MHz is applied, the plasma bulk region is reduced by the dc source, thus the ionization rate and hence the electron density decrease with rising dc voltage. However, when a DF source of 60 and 2 MHz is applied, the electron density can increase upon addition of a dc voltage, depending on the gap length and applied dc voltage. This is explained from the spatiotemporal ionization rates in the DF discharge. In fact, a completely different behavior is observed for the ionization rate in the two half-periods of the LF source. In the first LF half-period, the situation resembles the dc superposed SF discharge, and the reduced plasma bulk region due to the negative dc bias results in a very small effective discharge area and a low ionization rate. On the other hand, in the second half-period, the negative dc bias is to some extent counteracted by the LF voltage, and the sheath close to the dc electrode becomes particularly thin. Consequently, the amplitude of the high-frequency sheath oscillations at the top electrode is largely enhanced, while the LF sheath at the bottom electrode is in its expanding phase and can thus well confine the high-energy electrons. Therefore, the ionization rate increases considerably in this second LF half-period. Furthermore, in addition to the comparison between SF and DF discharges and the effect of gap length and dc voltage, the effect of secondary electrons is examined.

Abstract: The electron bounce resonance heating (BRH) in dual-frequency capacitively coupled plasmas operated in oxygen is studied by different experimental methods and a particle-in-cell/Monte Carlo collision (PIC/MCC) simulation, and compared with the electropositive argon discharge. In comparison with argon, the experimental results show that in an oxygen discharge the resonance peaks in positive-ion density and light intensity tend to occur at larger electrode gaps. Moreover, at electrode gaps L > 2.5 cm, the positive-ion (and electron) density and the light emission drop monotonically in the oxygen discharge upon increasing L, whereas they rise (after an initial drop) in the argon case. At resonance gap the electronegativity reaches its maximum due to the BRH. All these experimental observations are explained by PIC/MCC simulations, which show that in the oxygen discharge the bulk electric field becomes quite strong and is out of phase with the sheath field. Therefore, it retards the resonance electrons when traversing the bulk, resulting in a suppressed BRH. Both experiment and simulation results show that this effect becomes more pronounced at lower high-frequency power, when the discharge mode changes from electropositive to electronegative. In a pure oxygen discharge, the BRH is suppressed with increasing pressure and almost diminishes at 12 Pa. Finally, the driving frequency significantly affects the BRH, because it determines the phase relation between bulk electric field and sheath electric field.

Abstract: Recent advances in manufacturing arrays of artificial pinning sites, i.e., antidots, blind holes and magnetic dots, allowed an effective control of magnetic flux in superconductors. An array of magnetic bars deposited on top of a superconducting film was shown to display different pinning regimes depending on the direction of the in-plane magnetization of the bars. Changing the sign of their magnetization results in changes in the induced magnetic pinning potentials. By numerically solving the time-dependent Ginzburg-Landau equations in a superconducting film with periodic arrays of zigzag-arranged magnetic bars, we revealed various flux dynamics regimes. In particular, we demonstrate flux pinning and flux flow, depending on the direction of the magnetization of the magnetic bars. Remarkably, the revealed different flux-motion regimes are associated with different mechanisms of vortex-antivortex dynamics. For example, we found that for an 'antiparallel' configuration of magnetic bars this dynamics involves a repeating vortex-antivortex generation and annihilation. We show that the depinning transition and the onset of flux flow can be manipulated by the magnetization of the bars and the geometry of the array. This provides an effective control of the depinning critical current that can be useful for possible fluxonics applications.

Abstract: We investigate the coupling factor phi( mu) that quantifies the magnetic flux phi per magnetic moment mu of a point-like magnetic dipole that couples to a superconducting quantum interference device (SQUID). Representing the dipole by a tiny current-carrying (Amperian) loop, the reciprocity of mutual inductances of SQUID and Amperian loop provides an elegant way of calculating phi(mu)(r,e(mu)) vs. position r and orientation e(mu) of the dipole anywhere in space from the magnetic field B-J(r) produced by a supercurrent circulating in the SQUID loop. We use numerical simulations based on London and Ginzburg-Landau theory to calculate phi (mu) from the supercurrent density distributions in various superconducting loop geometries. We treat the far-field regime ( r greater than or similar to a= inner size of the SQUID loop) with the dipole placed on (oriented along) the symmetry axis of circular or square shaped loops. We compare expressions for phi (mu) from simple filamentary loop models with simulation results for loops with finite width w (outer size A > alpha), thickness d and London penetration depth lambda(L )and show that for thin ( d << alpha ) and narrow (w < alpha) loops the introduction of an effective loop size a(eff) in the filamentary loop-model expressions results in good agreement with simulations. For a dipole placed right in the center of the loop, simulations provide an expression phi(mu)(a,A,d,lambda(L)) that covers a wide parameter range. In the near-field regime (dipole centered at small distance z above one SQUID arm) only coupling to a single strip representing the SQUID arm has to be considered. For this case, we compare simulations with an analytical expression derived for a homogeneous current density distribution, which yields excellent agreement for lambda(L)>w,d . Moreover, we analyze the improvement of phi(mu) provided by the introduction of a narrow constriction in the SQUID arm below the magnetic dipole.

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: Non-equilibrium inductively coupled plasmas (ICPs) operating in hydrogen are of significant interest for applications including large-area materials processing. Increasing control of spatial gas heating, which drives the formation of neutral species density gradients and the rate of gas-temperature-dependent reactions, is critical. In this study, we use 2D fluid-kinetic simulations with the Hybrid Plasma Equipment Model to investigate the spatially resolved production of atomic hydrogen in a low-pressure planar ICP operating in pure hydrogen (10-20 Pa or 0.075-0.15 Torr, 300 W). The reaction set incorporates self-consistent calculation of the spatially resolved gas temperature and 14 vibrationally excited states. We find that the formation of neutral-gas density gradients, which result from spatially non-uniform electrical power deposition at constant pressure, can drive significant variations in the vibrational distribution function and density of atomic hydrogen when gas heating is spatially resolved. This highlights the significance of spatial gas heating on the production of reactive species in relatively high-power-density plasma processing sources.

Abstract: The possibility of manipulating vortex matter by using various artificial pinning arrays is of significant importance for possible applications in nano and micro fluxonics devices. By numerically solving the time-dependent GinzburgLandau equations, we study the vortexantivortex (vav) dynamics in a hybrid structure consisting of a finite-size superconductor with magnetic dipoles on top which generate vav pairs in the presence of an external current. The vav dynamics is analyzed for different arrangements and magnetic moments of the dipoles, as a function of angle α between the direction of the magnetic dipole and that of the Lorentz force produced by the applied current. The interplay of the attractive interaction between a vav pair and the Lorentz force leads either to the separation of (anti)vortices and their motion in opposite directions or to their annihilation. We found a critical angle αc, below which vortices and antivortices are repelled, while for larger angles they annihilate. In case of a single (few) magnetic dipole(s), this magnetic dipole induced vav guidance is influenced by the self-interaction of the vav pairs with their images in a finite-size sample, while for a periodic array of dipoles the guidance is determined by the interaction of a vav pair with other dipoles and vav pairs created by them. This effect is tunable through the external current and the magnetization and size of the magnetic dipoles.

Abstract: Using a combination of the phenomenological GinzburgLandau theory and micromagnetic simulations, we study properties of a superconducting film with an array of soft magnetic dots on top. An external in-plane magnetic field gradually drives the magnets from an out-of-plane or magnetic vortex state to an in-plane single-domain state, which changes spatially the distribution of the superconducting condensate. If induced by the magnets, the vortexantivortex molecules exhibit rich transitions as a function of the applied in-plane field. At the same time, we show how the magnetic dots act as very effective dynamic pinning centers for vortices in an applied perpendicular magnetic field.