Abstract: Silicene and germanene are the silicon and germanium counterparts of graphene, respectively. Recent experimental works have reported the growth of silicene on (1 1 1)Ag surfaces with different atomic configurations, depending on the growth temperature and surface coverage. We first theoretically study the structural and electronic properties of silicene on (1 1 1) Ag surfaces, focusing on the (4 x 4) silicene/Ag structure. Due to symmetry breaking in the silicene layer (nonequivalent number of top and bottom Si atoms), the corrugated silicene layer, with the Ag substrate removed, is predicted to be semiconducting, with a computed energy bandgap of about 0.3 eV. However, the hybridization between the Si 3p orbitals and the Ag 5s orbital in the silicene/(1 1 1)Ag slab model leads to an overall metallic system, with a distribution of local electronic density of states, which is related to the slightly disordered structure of the silicene layer on the (1 1 1)Ag surface. We next study the interaction of silicene and germanene with different hexagonal non-metallic substrates, namely ZnS and ZnSe. On reconstructed (0 0 0 1)ZnS or ZnSe surfaces, which should be more energetically stable for very thin layers, silicene and germanene are found to be semiconducting. Remarkably, the nature and magnitude of their energy bandgap can be controlled by an out-of-plane electric field, an important finding for the potential use of these materials in nanoelectronic devices. (C) 2013 Elsevier B. V. All rights reserved.

Abstract: The limitation in time and length scale is a major issue of molecular dynamics (MD) simulation. Although several methods have been developed to extend the MD time scale, their performance usually deteriorates with increasing system size. Therefore, an acceleration method which is applicable to large systems is required to bridge the gap between the MD simulations and target phenomena. In this study, an accelerated MD method for large system is developed based on the collective variable-driven hyperdynamics (CVHD) method [K.M. Bal and E.C. Neyts, 2015]. The key idea is to run CVHD in parallel with rate control and accelerate multiple possible events simultaneously. Using this novel method, carbon diffusion in bcc-iron bicrystal with grain boundary is examined as an application for practical materials. Carbon atoms reaching at the grain boundary are trapped whereas carbon atoms in the bulk region diffuse randomly, and both dynamic regimes can be simultaneously accelerated with the parallel CVHD technique.

Abstract: Oxocarbons, known for over two centuries, have recently revealed a long-awaited facet: two-dimensional crystalline structures. Employing an intelligent global optimization algorithm (IGOA) alongside densityfunctional calculations, we unearthed a quasi -flat oxocarbon (C 6 0 6 ), featuring an oxygen -decorated hole, and a novel 3D-borocarbon. Comparative analyses with recently synthesized isostructures, such as 2D -porous carbon nitride (C 6 N 6 ) and 2D -porous boroxine (B 6 0 6 ), highlight the unique attributes of these compounds. All structures share a common stoichiometry of X 6 Y 6 (which we call COF-66), where X = B, C, and Y = B, N, O (with X not equal Y), exhibiting a 2D -crystalline structure, except for borocarbon C 6 B 6 , which forms a 3D crystal. In our comprehensive study, we conducted a detailed exploration of the electronic structure of X 6 Y 6 compounds, scrutinizing their thermodynamic properties and systematically evaluating phonon stability criteria. With expansive surface areas, diverse pore sizes, biocompatibility, pi-conjugation, and distinctive photoelectric properties, these structures, belonging to the covalent organic framework (COF) family, present enticing prospects for fundamental research and hold potential for biosensing applications.

Abstract: The structural and transport properties of YBa2Cu3O7−x films grown by pulsed laser deposition with mixed 2.5 mol% Ba2YTaO6 (BYTO) and 2.5 mol% Ba2YNbO6 (BYNO) double-perovskite secondary phases are investigated in an extended film growth rate, R = 0.02–1.8 nm s−1. The effect of R on the film microstructure analyzed by TEM techniques shows an evolution from sparse and straight to denser, thinner and splayed continuous columns, with mixed BYNO + BYTO (BYNTO) composition, as R increases from 0.02 nm s−1 to 1.2 nm s−1. This microstructure results in very efficient flux pinning at 77 K, leading to a remarkable improvement in the critical current density (J c) behaviour, with the maximum pinning force density F p(Max) = 13.5 GN m−3 and the irreversibility field in excess of 11 T. In this range, the magnetic field values at which the F p is maximized varies from 1 T to 5 T, being related to the BYNTO columnar density. The film deposited when R = 0.3 nm s−1 exhibits the best performances over the whole temperature and magnetic field ranges, achieving F p(Max) = 900 GN m−3 at 10 K and 12 T. At higher rates, R > 1.2 nm s−1, BYNTO columns show a meandering nature and are prone to form short nanorods. In addition, in the YBCO film matrix a more disordered structure with a high density of short stacking faults is observed. From the analysis of the F p(H, T) curves it emerges that in films deposited at the high R limit, the vortex pinning is no longer dominated by BYNTO columnar defects, but by a new mechanism showing the typical temperature scaling law. Even though this microstructure produces a limited improvement at 77 K, it exhibits a strong J c improvement at lower temperature with F p = 700 GN m−3 at 10 K, 12 T and 900 GN m−3 at 4.2 K, 18 T.

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: Nanosecond pulsed discharge plasma shows a high degree of non-equilibrium, and exhibits relatively high conversions in the dry reforming of methane. To further improve the application, a good insight of the underlying mechanisms is desired. We developed a chemical kinetics model to explore the underlying plasma chemistry in nanosecond pulsed discharge. We compared the calculated conversions and product selectivities with experimental results, and found reasonable agreement in a wide range of specific energy input. Hence, the chemical kinetics model is able to provide insight in the underlying plasma chemistry. The modeling results predict that the most important dissociation reaction of CO₂and CH₄is electron impact dissociation. C₂H₂is the most abundant hydrocarbon product, and it is mainly formed upon reaction of two CH₂radicals. Furthermore, the vibrational excitation levels of CO₂contribute for 85% to the total dissociation of CO₂.

Abstract: Vibrational excitation represents an efficient channel to drive the dissociation of CO₂in a non-thermal plasma. Its viability is investigated in low-pressure pulsed discharges, with the intention of selectively exciting the asymmetric stretching mode, leading to stepwise excitation up to the dissociation limit of the molecule. Gas heating is crucial for the attainability of this process, since the efficiency of vibration–translation (V–T) relaxation strongly depends on temperature, creating a feedback mechanism that can ultimately thermalize the discharge. Indeed, recent experiments demonstrated that the timeframe of V–T non-equilibrium is limited to a few milliseconds at ca. 6 mbar, and shrinks to the<italic>μ</italic>s-scale at 100 mbar. With the aim of backtracking the origin of gas heating in pure CO₂plasma, we perform a kinetic study to describe the energy transfers under typical non-thermal plasma conditions. The validation of our kinetic scheme with pulsed glow discharge experiments enables to depict the gas heating dynamics. In particular, we pinpoint the role of vibration–vibration–translation relaxation in redistributing the energy from asymmetric to symmetric levels of CO₂, and the importance of collisional quenching of CO₂electronic states in triggering the heating feedback mechanism in the sub-millisecond scale. This latter finding represents a novelty for the modelling of low-pressure pulsed discharges and we suggest that more attention should be paid to it in future studies. Additionally, O atoms convert vibrational energy into heat, speeding up the feedback loop. The efficiency of these heating pathways, even at relatively low gas temperature and pressure, underpins the lifetime of V–T non-equilibrium and suggests a redefinition of the optimal conditions to exploit the ‘ladder-climbing’ mechanism in CO₂discharges.

Abstract: A gliding arc (GA) plasma has great potential for CO2 conversion into value-added chemicals, because of its high energy efficiency. To improve the application, a 2D/3D fluid model is needed to investigate the CO2 conversion mechanisms in the actual discharge geometry. Therefore, the complex CO2 chemical kinetics description must be reduced due to the huge computational cost associated with 2D/3D models. This paper presents a chemistry reduction method for CO2 plasmas, based on the so-called directed relation graph method. Depending on the defined threshold values, some marginal species are identified. By means of a sensitivity analysis, we can further reduce the chemistry set by removing one by one the marginal species. Based on the socalled flux-sensitivity coupling, we obtain a reduced CO2 kinetics model, consisting of 36 or 15 species (depending on whether the 21 asymmetric mode vibrational states of CO2 are explicitly included or lumped into one group), which is applied to a GA discharge. The results are compared with those predicted with the full chemistry set, and very good agreement is reached. Moreover, the range of validity of the reduced CO2 chemistry set is checked, telling us that this reduced set is suitable for low power GA discharges. Finally, the time and spatial evolution of the CO2 plasma characteristics are presented, based on a 2D model with the reduced kinetics.

Abstract: This study focuses on a magnetically stabilized gliding arc (MGA) plasma. Two fully coupled flow-plasma models (in 3D and 2D) are presented. The 3D model is applied to compare the arc dynamics of the MGA with a traditional gas-driven gliding arc. The 2D model is used for a detailed parametric study on the effect of the external magnetic field. The results show that the relative velocity between the plasma and feed gas is generated due to the Lorentz force, which can increase the plasma-treated gas fraction. The magnetic field also helps to decrease the gas temperature by enhancing heat transfer and to increase the electron number density. This work shows the potential of an external magnetic field to control the gliding arc behavior, for enhanced gas conversion at low gas flow rates.

Abstract: We investigate discharge contraction in a microwave plasma at sub-atmospheric pressure, operating in CO2 and CO2/CH4 mixtures. The rise of the electron number density with plasma contraction intensifies the gas heating in the core of the plasma. This, in turn, initiates fast core-periphery transport and defines the rate of thermal chemistry over plasma chemistry. In this context, power concentration describes the overall mechanism including plasma contraction and chemical kinetics. In a complex chemistry such as dry reforming of methane, transport of reactive species is essential to define the performance of the reactor and achieve the desired outputs. Thus, we couple experimental observations and thermodynamic calculations for model validation and understanding of reactor performance. Adding CH4 alters the thermodynamic properties of the mixture, especially the reactive component of the heat conductivity. The increase in reactive heat conductivity increases the pressure at which plasma contraction occurs, because higher rates of gas heating are required to reach the same temperature. In addition, we suggest that the predominance of heat conduction over convection is a key condition to observe the effect of heat conductivity on gas temperature.

Abstract: Since the onset of pattern transfer technologies for chip manufacturing, various strategies have been developed to circumvent or overcome aspect ratio dependent etching (ARDE). These methods have, however, their own limitations in terms of etch non-idealities, throughput or costs. Moreover, they have mainly been optimized for individual in-device features and die-scale patterns, while occasionally ending up with poor patterning of metrology marks, affecting the alignment and overlay in lithography. Obtaining a better understanding of the underlying mechanisms of ARDE and how to mitigate them therefore remains a relevant challenge to date, for both marks and advanced nodes. In this work, we accordingly assessed the neutral transport mechanisms in ARDE by means of experiments and multiscale modeling for SiO₂etching with CHF₃/Ar and CF₄/Ar plasmas. The experiments revealed a local maximum in the etch rate for an aspect ratio around unity, i.e. the simultaneous occurrence of regular and inverse reactive ion etching lag for a given etch condition. We were able to reproduce this ARDE trend in the simulations without taking into account charging effects and the polymer layer thickness, suggesting shadowing and diffuse reflection of neutrals as the primary underlying mechanisms. Subsequently, we explored four methods with the simulations to regulate ARDE, by varying the incident plasma species fluxes, the amount of polymer deposition, the ion energy and angular distribution and the initial hardmask sidewall angle, for which the latter was found to be promising in particular. Although our study focusses on feature dimensions characteristic to metrology marks and back-end-of-the-line integration, the obtained insights have a broader relevance, e.g. to the patterning of advanced nodes. Additionally, this work supports the insight that physisorption may be more important in plasma etching at room temperature than originally thought, in line with other recent studies, a topic on which we recommend further research.

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: Correction for ‘From the Birkeland–Eyde process towards energy-efficient plasma-based NO_{<italic>X</italic>}synthesis: a techno-economic analysis’ by Kevin H. R. Rouwenhorst<italic>et al.</italic>,<italic>Energy Environ. Sci.</italic>, 2021,<bold>14</bold>, 2520–2534, https://doi.org/10.1039/D0EE03763J.

Abstract: The response of electron systems to electrodynamic fields that change rapidly in space is endowed by unique features, including an exquisite spatial nonlocality. This can reveal much about the materials' electronic structure that is invisible in standard probes that use gradually varying fields. Here, we use graphene plasmons, propagating at extremely slow velocities close to the electron Fermi velocity, to probe the nonlocal response of the graphene electron liquid. The near-field imaging experiments reveal a parameter-free match with the full quantum description of the massless Dirac electron gas, which involves three types of nonlocal quantum effects: single-particle velocity matching, interaction-enhanced Fermi velocity, and interaction-reduced compressibility. Our experimental approach can determine the full spatiotemporal response of an electron system.