Abstract: Two-dimensional (2D) honeycomb ferromagnets, such as monolayer chromium trihalides, are predicted to behave as topological magnon insulators, characterized by an insulating bulk and topologically protected edge states, giving rise to a thermal magnon Hall effect. Here we report the behavior of the topological magnons in monolayer CrI3 under external stimuli, including biaxial and uniaxial strain, electric gating, as well as in-plane and out-of-plane magnetic field, revealing that one can thereby tailor the magnetic states as well as the size and the topology of the magnonic bandgap. These findings broaden the perspective of using 2D magnetic materials to design topological magnonic devices.

Abstract: A mode converter for electron vortex beams is described. Numerical simulations, confirmed by experiment, show that the converter transforms a vortex beam with a topological charge m = +/- 1 into beams closely resembling Hermite-Gaussian HG(10) and HG(01) modes. The converter can be used as a mode discriminator or filter for electron vortex beams. Combining the converter with a phase plate turns a plane wave into modes with topological charge m = +/- 1. This combination serves as a generator of electron vortex beams of high brilliance.

Abstract: The dynamics of the pinned vortex, antivortex and interstitial vortex have been studied in superconducting/magnetic hybrids consisting of arrays of Co/Pd multilayer nanodots embedded in Nb films. The magnetic nanodots show out-of-plane magnetization at the remanent state. This magnetic state allows for superconducting vortex lattices of different types in an applied homogeneous magnetic field. We experimentally and theoretically show three such lattices: (i) a lattice containing only antivortices; (ii) a vortex lattice entirely pinned on the dots; and (iii) a vortex lattice with pinned and interstitial vortices. Between the flux creep (low vortex velocity) and the free flux flow (high vortex velocity) regimes the interaction between the magnetic array and the vortex lattice governs the vortex dynamics, which in turn enables distinguishing experimentally the type of vortex lattice which governs the dissipation. We show that the vortex lattice with interstitial vortices has the highest onset velocity where the lattice becomes ordered, whereas the pinned vortex lattice has the smallest onset velocity. Further, for this system, we directly estimate that the external force needed to depin vortices is 60% larger than the one needed to depin antivortices; therefore we are able to decouple the antivortex-vortex motion.

Abstract: Electron-energy-loss spectroscopy (EELS) is used to investigate single layers of LaAlO3 grown on SrTiO3 having an n-type interface as well as multilayers of LaAlO3 and SrTiO3 in which both n- and p-type interfaces occur. Only minor changes in Ti valence at the n-type interface are observed. This finding seems to contradict earlier experiments for other SrTiO3/LaAlO3 systems where large deviations in Ti valency were assumed to be responsible for the conductivity of these interfaces. Ab initio calculations have been carried out in order to interpret our EELS results. Using the concept of Bader charges, it is demonstrated that the so-called polar discontinuity is mainly resolved by lattice distortions and to a far lesser extent by changes in valency for both single layer and multilayer geometries.

Abstract: Methane needle-plane discharge has practical application prospect and scientific research significance since methane conversion heavy oil hydrogenation is formed by coupling methane needle-plane discharge with heavy oil hydrogenation, which can achieve high-efficient heavy oil hydrogenation and increase the yields of high value-added light olefins. In this paper, a two-dimensional fluid model is built up for numerically simulating the methane needle-plane discharge plasma at atmospheric pressure. Spatial and axial distributions of electric intensity, electron temperature and particle densities are obtained. Reaction yields are summarized and crucial pathways to produce various kinds of charged and neutral particles are found out. Simulation results indicate that axial evolutions of CH3+ and CH4+ densities, electric intensity and electron temperature are similar and closely related. The CH5+ and C2H5+ densities first increase and then decrease along the axial direction. The CH3 and H densities have nearly identical spatial and axial distributions. Particle density distributions of CH2, C2H4 and C2H5 are obviously different in the area near the cathode but comparatively resemblant in the positive column region. The CH3+ and CH4+ are produced by electron impact ionizations between electrons and CH4. The CH5+ and C2H5+ are respectively generated by molecular impact dissociations between CH3+ and CH4 and between CH4+ and CH4. Electron impact decomposition between electrons and CH4 is a dominated reaction to produce CH3, CH2, CH and H. The reactions between CH2 and CH4 and between electrons and C2H4 are critical pathways to produce C2H4 and C2H2, respectively. In addition, the yields of electron impact decomposition reactions between electrons and CH4 and reactions between CH2 and CH4 account for 52.15% and 47.85% of total yields of H-2 respectively.

Abstract: Coupled shallow impurity states in a freestanding semiconductor nanowire and in a semiconductor nanowire surrounded by a metallic gate are studied within the effective-mass approximation. Bonding and antibonding states are found due to the coupling of the two impurities, and their energy converges with increasing distance di between the two impurities. The dependences of the binding energy on the wire radius R, the distance di between the two impurities, and the impurity radial position in the nanowire are examined.

Abstract: In this paper we investigate the basic physics of charge carriers (electrons) leaking out of the inversion layer of a metal-oxide-semiconductor capacitor with a biased gate. In particular, we treat the gate leakage current as resulting from two combined processes: (1) the time-dependent decay of electron wave packets representing the inversion-layer charge and (2) the local generation of new electrons replacing those that have leaked away. As a result, the gate current simply emerges as the ratio of the total charge in the inversion layer to the tunneling lifetime. The latter is extracted from the quantum dynamics of the decaying wave packets, while the generation rate is incorporated as a phenomenological source term in the continuity equation. Not only do the gate currents calculated with this model agree very well with experiment, the model also provides an onset to solve the paradox of the current-free bound states representing the resonances of the Schrödinger equation that governs the fully coupled metal-oxide-semiconductor system.

Abstract: Using the tight-binding approach, we study the influence of a perpendicular magnetic field on the energy levels of hexagonal, triangular, and circular bilayer graphene (BLG) quantum dots (QDs) with zigzag and armchair edges. We obtain the energy levels for AB (Bernal)-stacked BLG QDs in both the absence and the presence of a perpendicular electric field (i.e., biased BLG QDs). We find different regions in the spectrum of biased QDs with respect to the crossing point between the lowest-electron and -hole Landau levels of a biased BLG sheet. Those different regions correspond to electron states that are localized at the center, edge, or corner of the BLG QD. Quantum Hall corner states are found to be absent in circular BLG QDs. The spatial symmetry of the carrier density distribution is related to the symmetry of the confinement potential, the position of zigzag edges, and the presence or absence of interlayer inversion symmetry.

Abstract: By making use of circularly polarized light and electrostatic gating, monolayer molybdenum disulfide (ML – MoS2) can form a platform supporting multiple types of charge carriers. They can be discriminated by their spin, valley index, or whether they are electrons or holes. We investigate the collective properties of those charge carriers and are able to identify distinct plasmon modes. We analyze the corresponding dispersion relation, lifetime, and oscillator strength, and calculate the phase relation between the oscillations in the different components of the plasmon modes. All platforms in ML-MoS2 support a long-wavelength root q plasmon branch at zero kelvins. In addition to this, for an n-component system, n-1 distinct plasmon modes appear as acoustic modes with linear dispersion in the long-wavelength limit. These modes correspond to out-of-phase oscillations in the different fermion liquids and have, although being damped, a relatively long lifetime. Additionally, we also find distinct modes at large wave vectors that are more strongly damped by intraband processes.

Abstract: When slow ions incident on a surface are neutralized, the excess potential energy is passed on to an electron inside the surface, leading to emission of secondary electrons. The microscopic description of this process, as
well as the calculation of the secondary electron yield, is a challenging problem due to its complexity as well
as its sensitivity to surface properties. One of the first quantitative descriptions was articulated in the 1950s by
Hagstrum, who based his calculation on a parametrization of the density of states of the material. In this paper, we
present a model for calculating the secondary electron yield, derived from Hagstrum’s initial approach. We use
first-principles density functional theory calculations to acquire the necessary input and introduce the concept of
electron cascades to Hagstrum’s model in order to improve the calculated spectra, as well as remove its reliance
on fitting parameters. We apply our model to He+ and Ne+ ions incident on Ge(111) and Si(111) and obtain
yield spectra that match closely to the experimental results of Hagstrum.

Abstract: Using the tight-binding (TB) approximation with inclusion of the spin-orbit interaction, we predict a topological phase transition in the electronic band structure of phosphorene in the presence of axial strains. We derive a low-energy TB Hamiltonian that includes the spin-orbit interaction for bulk phosphorene. Applying a compressive biaxial in-plane strain and perpendicular tensile strain in ranges where the structure is still stable leads to a topological phase transition. We also examine the influence of strain on zigzag phosphorene nanoribbons (zPNRs) and the formation of the corresponding protected edge states when the system is in the topological phase. For zPNRs up to a width of 100 nm the energy gap is at least three orders of magnitude larger than the thermal energy at room temperature.

Abstract: Recent experimental studies showed that formation of a two-dimensional crystal structure of nitrogenated holey graphene (NHG) is possible. Similar to graphene, NHGs have an atomically thin and strong crystal structure. Using first-principles calculations, we investigate the structural, phononic, and thermal properties of monolayer NHG crystal. Our charge analysis reveals that the charged holey sites of NHG provide a reactive ground for further functionalization by adatoms or molecules. We also found that similar to graphene, the NHG structure has quite high-frequency phonon modes and the presence of nitrogen atoms leads to the emergence of additional vibrational modes. Our phonon analysis reveals the presence of three characteristic Raman-active modes of NHG. Furthermore, the analysis of constant-volume heat capacity showed that the NHG structure has a linear temperature dependence in the low-temperature region. The strong lattice structure and unique thermal properties of the NHG crystal structure are desirable in nanoscale device applications.

Abstract: We present a unified theory of the phonon dispersions and elastic properties of graphene, graphite, and graphene multilayer systems. Starting from a fifth-nearest-neighbor force-constant model derived from full in-plane phonon dispersions of graphite [Mohr et al., Phys. Rev. B 76, 035439 (2007)], we use Born's long-wave method to calculate the tension and bending coefficients of graphene. Extending the model by interplanar interactions, we study the phonon dispersions and the elastic constants of graphite, and the phonon spectra of graphene multilayers. We find that the inner displacement terms due to sublattice shifts between inequivalent C atoms are quantitatively important in determining the elastomechanical properties of graphene and of graphite. The overall agreement between theory and experiment is very satisfactory. We investigate the evolution from graphene to graphite by studying the increase in the rigid plane optical mode as a function of the number of layers N. At N=10 the graphite value B2g1127 cm−1 is attained within a few percent.

Abstract: Motivated by a recent experiment that reported the successful synthesis of hexagonal (h) AlN [Tsipas et al., Appl. Phys. Lett. 103, 251605 (2013)], we investigate structural, electronic, and vibrational properties of bulk, bilayer, and monolayer structures of h-AlN by using first-principles calculations. We show that the hexagonal phase of the bulk h-AlN is a stable direct-band-gap semiconductor. The calculated phonon spectrum displays a rigid-layer shear mode at 274 cm(-1) and an E-g mode at 703 cm(-1), which are observable by Raman measurements. In addition, single-layer h-AlN is an indirect-band-gap semiconductor with a nonmagnetic ground state. For the bilayer structure, AA'-type stacking is found to be the most favorable one, and interlayer interaction is strong. While N-layered h-AlN is an indirect-band-gap semiconductor for N = 1 – 9, we predict that thicker structures (N >= 10) have a direct band gap at the Gamma point. The number-of-layer-dependent band-gap transitions in h-AlN is interesting in that it is significantly different from the indirect-to-direct crossover obtained in the transition-metal dichalcogenides.

Abstract: Using first-principles calculations, we study the dependence of the electronic and vibrational properties of multilayered PbI2 crystals on the number of layers and focus on the electronic-band structure and the Raman spectrum. Electronic-band structure calculations reveal that the direct or indirect semiconducting behavior of PbI2 is strongly influenced by the number of layers. We find that at 3L thickness there is a direct-to-indirect band gap transition (from bulk-to-monolayer). It is shown that in the Raman spectrum two prominent peaks, A(1g) and E-g, exhibit phonon hardening with an increasing number of layers due to the interlayer van der Waals interaction. Moreover, the Raman activity of the A(1g) mode significantly increases with an increasing number of layers due to the enhanced out-of-plane dielectric constant in the few-layer case. We further characterize rigid-layer vibrations of low-frequency interlayer shear (C) and breathing (LB) modes in few-layer PbI2. A reduced monatomic (linear) chain model (LCM) provides a fairly accurate picture of the number of layers dependence of the low-frequency modes and it is shown also to be a powerful tool to study the interlayer coupling strength in layered PbI2.