Abstract: In this paper, simulations are performed to gain a better insight into the properties of a Cl2/Ar plasma, with and without O2, during plasma etching of Si. Both plasma and surface properties are calculated in a self-consistent manner. Special attention is paid to the behavior of etch products coming from the wafer or the walls, and how the chamber walls can affect the plasma and the resulting etch process. Two modeling cases are considered. In the first case, the reactor walls are defined as clean (Al2O3), whereas in the second case a SiO2 coating is introduced on the reactor walls before the etching process, so that oxygen will be sputtered from the walls and introduced into the plasma. For this reason, a detailed reaction set is presented for a Cl2/O2/Ar plasma containing etched species, as well as an extensive reaction set for surface processes, including physical and chemical sputtering, chemical etching and deposition processes. Density and flux profiles of various species are presented for a better understanding of the bulk plasma during the etching process. Detailed information is also given on the composition of the surfaces at various locations of the reactor, on the etch products in the plasma and on the surface loss probabilities of the plasma species at the walls, with different compositions. It is found that in the clean chamber, walls are mostly chlorinated (Al2Cl3), with a thin layer of etch products residing on the wall. In the coated chamber, an oxy-chloride layer is grown on the walls for a few nanometers during the etching process. The Cl atom wall loss probability is found to decrease significantly in the coated chamber, hence increasing the etch rate. SiCl2, SiCl4 and SiCl3 are found to be the main etch products in the plasma, with the fraction of SiCl2 being always slightly higher. The simulation results compare well with experimental data available from the literature.

Abstract: Dielectric barrier discharges (DBDs) are commonly used to generate cold plasmas at
atmospheric pressure. Whatever their configuration (tubular or planar), the presence of a dielectric barrier is mandatory to prevent too much charge build up in the plasma and the formation of a thermal arc. In this article, the role of the barrier thickness (2.0, 2.4 and 2.8 mm) and of the kind of dielectric material (alumina, mullite, pyrex, quartz) is investigated on the filamentary behavior in the plasma and on the CO2 conversion in a tubular flowing DBD, by means of mass spectrometry measurements correlated with electrical characterization and IR imaging. Increasing the barrier thickness decreases the capacitance, while preserving the electrical charge. As a result, the voltage over the dielectric increases and a larger number of microdischarges is generated, which enhances the CO2 conversion. Furthermore, changing the dielectric material of the barrier, while keeping the same geometry and dimensions, also affects the CO2 conversion. The highest CO2 conversion and energy efficiency are obtained for quartz and alumina, thus not following the trend of the relative permittivity. From the
electrical characterization, we clearly demonstrate that the most important parameters are the somewhat higher effective plasma voltage (yielding a somewhat higher electric field and electron energy in the plasma) for quartz, as well as the higher plasma current (and thus larger electron density) and the larger number of microdischarge filaments (mainly for alumina, but also for quartz). The latter could be correlated to the higher surface roughness for alumina and to the higher voltage over the dielectric for quartz.

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: Although CO2 conversion by plasma technology is gaining increasing interest, the
underlying mechanisms for an energy-efficient process are still far from understood. In this work, a reduced non-equilibrium CO2 plasma chemistry set, based on level lumping of the vibrational levels, is proposed and the reliability of this level-lumping method is tested by a self-consistent zero-dimensional code. A severe reduction of the number of equations to be solved is achieved, which is crucial to be able to model non-equilibrium CO2 plasmas by 2-dimensional models. Typical conditions of pressure and power used in a microwave plasma for CO2 conversion are investigated. Several different sets, using different numbers of lumped groups, are considered. The lumped models with 1, 2 or 3 groups are able to reproduce the gas temperature, electron density and electron temperature profiles, as calculated by the full model treating all individual excited levels, in the entire pressure range investigated. Furthermore, a 3-groups model is also able to reproduce the shape of the vibrational distribution function (VDF) and gives the most reliable prediction of the CO2 conversion. A strong influence of the vibrational excitation on the plasma characteristics is observed. Finally, the limitations of the lumped-levels method are discussed.

Abstract: We show that a domain wall separating single-layer graphene and AA-stacked bilayer graphene (AA-BLG) can be used to generate highly collimated electron beams which can be steered by a magnetic field. Two distinct configurations are studied, namely, locally delaminated AA-BLG and terminated AA-BLG whose terminal edge types are assumed to be either zigzag or armchair. We investigate the electron scattering using semiclassical dynamics and verify the results independently with wave-packet dynamics simulations. We find that the proposed system supports two distinct types of collimated beams that correspond to the lower and upper cones in AA-BLG. Our computational results also reveal that collimation is robust against the number of layers connected to AA-BLG and terminal edges.

Abstract: Amorphous semiconductors are usually characterized by a low charge carrier mobility, essentially related to their lack of long-range order. The development of such material with higher charge carrier mobility is hence challenging. Part of the issue comes from the difficulty encountered by first-principles simulations to evaluate concepts such as the electron effective mass for disordered systems since the absence of periodicity induced by the disorder precludes the use of common concepts derived from condensed matter physics. In this paper, we propose a methodology based on first-principles simulations that partially solves this problem, by quantifying the degree of delocalization of a wave function and of the connectivity between the atomic sites within this electronic state. We validate the robustness of the proposed formalism on crystalline and molecular systems and extend the insights gained to disordered/amorphous InGaZnO4 and Si. We also explore the properties of p-type oxide semiconductor candidates recently reported to have a low effective mass in their crystalline phases [G. Hautier et al., Nat. Commun. 4, 2292 (2013)]. Although in their amorphous phase none of the candidates present a valence band with delocalization properties matching those found in the conduction band of amorphous InGaZnO4, three of the seven analyzed materials show some potential. The most promising candidate, K2Sn2O3, is expected to possess in its amorphous phase a slightly higher hole mobility than the electron mobility in amorphous silicon.

Abstract: <script type='text/javascript'>document.write(unpmarked('We show that delaminations in bilayer graphene (BLG) with electrostatically induced interlayer symmetry can provide one with ballistic channels for electrons with energies inside the electrostatically induced BLG gap. These channels are formed by a combination of valley-polarized evanescent states propagating along the delamination edges (which persist in the presence of a strong magnetic field) and standing waves bouncing between them inside the delaminated region (in a strong magnetic field, these transform into Landau levels in the monolayers). For inverted stackings in BLGs on the left and right of the delamination (AB-2ML-BA or BA-2ML-AB, where 2ML indicates two decoupled monolayers of graphene), the lowest-energy ballistic channels are gapless, have linear dispersion, and appear to be weakly topologically protected. When BLG stackings on both sides of the delamination are the same (AB-2ML-AB or BA-2ML-BA), the lowest-energy ballistic channels are gapped, with a gap epsilon(g) scaling as epsilon(g) alpha W-1 with delamination width and epsilon(g) alpha delta(-1) with the on-layer energy difference in the delaminated part of the structure. Depending on the width, delaminations may also support several \u0022higher-energy\u0022 waveguide modes. Our results are based on both the analytical study of the wave matching of Dirac states and tight-binding model calculations, and we analyze in detail the dependence of the delamination spectrum on the electrostatic conditions in the structure, such as the vertical displacement field.'));

Abstract: Vibrational thermal properties of CuZn2InTe4, AgZn2InTe4, and Cu2CdSnTe4, derived from binary II-VI zinc-blendes, are reported based on first-principles calculations. While the chalcogenide atoms in these materials have the same lattice positions, the cation atom arrangements vary, resulting in different crystal symmetries and subsequent properties. The compositional differences have important effects on the vibrational thermal characteristics of the studied materials, which demonstrate that low-frequency optical phonons hybridize with acoustic phonons and lead to enhanced phonon-phonon scattering and low lattice thermal conductivities. The phonon density of states, mode Gruneisen parameters, and phonon scattering rates are also calculated, enabling deeper insight into the microscopic thermal conduction processes in these materials. Compositional variations drive differences among the three materials considered here; nonetheless, their structural similarities and generally low thermal conductivities (0.5-4 W/mK at room temperature) suggest that other similar II-VI zinc-blende derived materials will also exhibit similarly low values, as also corroborated by experimental data. This, combined with the versatility in designing a variety of motifs on the overall structure, makes quaternary chalcogenides interesting for thermal management and energy conversion applications that require low thermal conductivity.

Abstract: The magnetic-field dependence of the energy spectrum, wave function, binding energy, and oscillator strength of exciton states confined in a circular graphene quantum dot (CGQD) is obtained within the configuration interaction method. We predict that (i) excitonic effects are very significant in the CGQD as a consequence of a combination of geometric confinement, magnetic confinement, and reduced screening; (ii) two types of excitons (intravalley and intervalley) are present in the CGQD because of the valley degree of freedom in graphene; (iii) the intravalley and intervalley exciton states display different magnetic-field dependencies due to the different electron-hole symmetries of the single-particle energy spectra; (iv) with increasing magnetic field, the exciton ground state in the CGQD undergoes an intravalley to intervalley transition accompanied by a change of angular momentum; (v) the exciton binding energy does not increase monotonically with the magnetic field due to the competition between geometric and magnetic confinements; and (vi) the optical transitions of the intervalley and intravalley excitons can be tuned by the magnetic field, and valley-dependent excitonic transitions can be realized in a CGQD.

Abstract: Ripples in pristine freestanding graphene naturally orient themselves in an array that is alternately curved-up and curved-down; maintaining an average height of zero. Using scanning tunneling microscopy (STM) to apply a local force, the graphene sheet will reversibly rise and fall in height until the height reaches 60%-70% of its maximum at which point a sudden, permanent jump occurs. We successfully model the ripples as a spin-half Ising magnetic system, where the height of the graphene plays the role of the spin. The permanent jump in height, controlled by the tunneling current, is found to be equivalent to an antiferromagnetic-to-ferromagnetic phase transition. The thermal load underneath the STM tip alters the local tension and is identified as the responsible mechanism for the phase transition. Four universal critical exponents are measured from our STM data, and the model provides insight into the statistical role of graphene's unusual negative thermal expansion coefficient.

Abstract: By performing density functional theory based first-principles calculations, the electronic, vibrational, elastic, and piezoelectric properties of two dynamically stable crystal phases of monolayer Janus MoSTe, namely 1H-MoSTe and 1T'-MoSTe, are investigated. Vibrational frequency analysis reveals that the other possible crystal structure, 1T-MoSTe, of this Janus monolayer does not exhibit dynamical stability. The 1H-MoSTe phase is found to be an indirect band-gap semiconductor while 1T'-MoSTe is predicted as small-gap semiconductor. Notably, in contrast to the direct band-gap nature of monolayers 1H-MoS2 and 1H-MoTe2, 1H-MoSTe is found to be an indirect gap semiconductor driven by the induced surface strains on each side of the structure. The calculated Raman spectrum of each structure shows unique character enabling us to clearly distinguish the stable crystal phases via Raman measurements. The systematic piezoelectric stress and strain coefficient analysis reveals that out-of-plane piezoelectricity appears in 1H-MoSTe and the noncentral symmetric 1T'-MoSTe has large piezoelectric coefficients. Static total-energy calculations show clearly that the formation of 1T'-MoSTe is feasible by using 1T'-MoTe2 as a basis monolayer. Therefore, we propose that the Janus MoSTe structure can be fabricated in two dynamically stable phases which possess unique electronic, dynamical, and piezoelectric properties.