Abstract: The effect of the number of stacking layers and the type of stacking on the electronic and optical properties of bilayer and trilayer black phosphorus are investigated by using first-principles calculations within the framework of density functional theory. We find that inclusion of many-body effects (i.e., electron-electron and electron-hole interactions) modifies strongly both the electronic and optical properties of black phosphorus. While trilayer black phosphorus with a particular stacking type is found to be a metal by using semilocal functionals, it is predicted to have an electronic band gap of 0.82 eV when many-body effects are taken into account within the G(0)W(0) scheme. Though different stacking types result in similar energetics, the size of the band gap and the optical response of bilayer and trilayer phosphorene are very sensitive to the number of layers and the stacking type. Regardless of the number of layers and the type of stacking, bilayer and trilayer black phosphorus are direct band gap semiconductors whose band gaps vary within a range of 0.3 eV. Stacking arrangements that are different from the ground state structure in both bilayer and trilayer black phosphorus exhibit significant modified valence bands along the zigzag direction and result in larger hole effective masses. The optical gap of bilayer (trilayer) black phosphorus varies by 0.4 (0.6) eV when changing the stacking type. The calculated binding energy of the bound exciton hardly changes with the type of stacking and is found to be 0.44 (0.30) eV for bilayer (trilayer) phosphorous.

Abstract: The effect of vacancy defects on the structural properties and the thermal stability of free standing silicene – a buckled structure of hexagonally arranged silicon atoms – is studied using reactive molecular dynamics simulations. Pristine silicene is found to be stable up to 1500 K, above which the system transits to a three-dimensional amorphous configuration. Vacancy defects result in local structural changes in the system and considerably reduce the thermal stability of silicene: depending on the size of the vacancy defect, the critical temperature decreases by more than 30%. However, the system is still found to be stable well above room temperature within our simulation time of 500 ps. We found that the, stability of silicene can be increased by saturating the dangling bonds at the defect edges by foreign atoms (e.g., hydrogen).

Abstract: The effects due to intersubband coupling and screening on the ionized impurity scattering are studied for a quasi-two-dimensional electron system in delta-doped semiconductors. We found that intersubband coupling plays an essential role in describing the screening properties and the effect of ionized impurity scattering on the mobility in a multisubband system. At the onset of the occupation of a higher subband, the screening due to the intersubband coupling leads to a reduction of the small angle scattering rate in the lower subband. We showed that such an effect is significant in a delta-doped quantum well and results in a pronounced increase of the quantum mobility at the onset of the occupation of a higher subband. (C) 1996 American Institute of Physics.

Abstract: The effects of interface optical-phonon and confined slab LO-phonon modes on the polaron cyclotron-resonance frequency are investigated for a GaAs/AlAs quantum well. Using degenerate second-order perturbation theory, the polaron Landau levels are calculated and the polaron resonant region is investigated. In order to know the relative importance of the different resonant frequencies we present a full calculation of the magneto-optical absorption spectrum. At a fixed magnetic field we found four different peaks in the absorption spectrum. The relative oscillator strength of the different peaks changes with increasing magnetic field. For comparative purposes, the polaron Landau levels and cyclotron mass are also calculated using only the bulk LO-phonon modes. The influence of the finiteness of the confinement potential is investigated. We found that the interface-phonon modes influence the magnetopolaron resonance considerably near the optical-phonon frequencies for narrow wells. In the limit of zero magnetic field we recover our previous results and in the case of an infinite-barrier quantum well we are able to recover the results for a two- and three-dimensional system.

Abstract: The elastic constant C-11 and piezoelectric stress constant e(1),(11) of two-dimensional (2D) dielectric materials comprising h-BN, 2H-MoS2, and other transition-metal dichalcogenides and dioxides are calculated using lattice dynamical theory. The results are compared with corresponding quantities obtained with ab initio calculations. We identify the difference between clamped-ion and relaxed-ion contributions with the dependence on inner strains which are due to the relative displacements of the ions in the unit cell. Lattice dynamics allows us to express the inner-strain contributions in terms of microscopic quantities such as effective ionic charges and optoacoustical couplings, which allows us to clarify differences in the piezoelectric behavior between h-BN and MoS2. Trends in the different microscopic quantities as functions of atomic composition are discussed.

Abstract: The electrical transport properties of a MoGe thin film with a honeycomb array of nanoscale holes are investigated. The critical current of the system shows nonmatching anomalies as a function of applied magnetic field, enabling us to distinguish between multiquanta vortices trapped in the holes and interstitial vortices located between the holes. The number of vortices trapped in each hole is found to be larger than the saturation number predicted for an isolated hole and shows a nonlinear field dependence, leading to the caging effect as predicted from the Ginzburg-Landau (GL) theory. Our experimental results are supplemented by numerical simulations based on the GL theory.

Abstract: The electron and hole spectrum in single and bilayer graphene is derived from known one-dimensional models, and the relation between the spectrum and symmetry of the lattice is shown.

Abstract: The electron states in axially symmetric quantum wires are computed by means of the effective-mass Schrodinger equation, which is written in cylindrical coordinates phi, rho, and z. We show that a direct discretization of the Schrodinger equation by central finite differences leads to a non-symmetric Hamiltonian matrix. Because diagonalization of such matrices is more complex it is advantageous to transform it in a symmetric form. This can be done by the Liouville-like transformation proposed by Rizea et al. (2008), which replaces the wave function psi(rho) with the function F(rho) = psi(rho)root rho and transforms the Hamiltonian accordingly. Even though a symmetric Hamiltonian matrix is produced by this procedure, the computed wave functions are found to be inaccurate near the origin, and the accuracy of the energy levels is not very high. In order to improve on this, we devised a finite-difference scheme which discretizes the Schrodinger equation in the first step, and then applies the Liouville-like transformation to the difference equation. Such a procedure gives a symmetric Hamiltonian matrix, resulting in an accuracy comparable to the one obtained with the finite element method. The superior efficiency of the new finite-difference scheme (FDM) is demonstrated for a few p-dependent one-dimensional potentials which are usually employed to model the electron states in free-standing and core shell quantum wires. The new scheme is compared with the other FDM schemes for solving the effective-mass Schrodinger equation, and is found to deliver energy levels with much smaller numerical error for all the analyzed potentials. It also gives more accurate results than the scheme of Rizea et al., except for the ground state of an infinite rectangular potential in freestanding quantum wires. Moreover, the PT symmetry is invoked to explain similarities and differences between the considered FDM schemes. (C) 2015 Elsevier B.V. All rights reserved.

Abstract: The electron transport properties in delta-doped semiconductor systems-are studied. The subband electronic structure of the delta-doped system is obtained by solving the coupled Schrodinger and Poisson equations. The screening of the quasi-two-dimensional electron gas is taken into account for the ionized impurity scattering through the matrix dielectric function within the random-phase approximation. The quantum and transport mobilities are calculated numerically as a function of the total electron density and the width of the doped layer at zero temperature. The intersubband scattering and the effect of empty subbands above the Fermi level on the electron mobilities are investigated. The calculated mobilities are in reasonable agreement with the available experimental results.

Abstract: The electronic and optical properties of zincblende nanowires are investigated in the presence of a uniform magnetic field directed along the [001] growth direction within the k . p method. We focus our numerical study on core-shell nanowires consisting of the III-V materials GaAs, AlxGa1-xAs and AlyGa1-y/0.51In0.49P. Nanowires with electrons confined in the core exhibit a Fock-Darwin-like spectrum, whereas nanowires with electrons confined in the shell show Aharonov-Bohm oscillations. Thus, by properly choosing the core and the shell materials of the nanowire, the optical properties in a magnetic field can be tuned in very different ways.

Abstract: The electronic and transport properties of a three-dimensional (3D) Dirac system are investigated theoretically, which is motivated by recent experimental measurements on quantum and transport mobilities in the 3D Dirac semimetal Na3Bi by J. Xiong et al. [Science 350, 413 (2015); Europhys. Lett. 114, 27002 (2016)]. The electron Hamiltonian is taken from a simplified k center dot p approach. From the obtained electronic band structure and the Fermi energy, we explain why the anomalous effect induced by the chiral anomaly and the Berry curvature in the energy band can be observed experimentally in magnetotransport coefficients in both low-and high-density samples. Moreover, the quantum and transport mobilities are calculated on the basis of the momentum-balance equation derived from a semiclassical Boltzmann equation with the electron-impurity interaction. The quantum and transport mobilities obtained from this study agree both qualitatively and quantitatively with those measured experimentally. We also examine the electron mobilities along different crystal directions in Na3Bi and find them largely anisotropic. The theoretical findings from this work can be helpful in gaining an in-depth understanding of the experimental results and of the basic electronic and transport properties of newly developed 3D Dirac systems.

Abstract: The electronic properties of a triaxially strained hexagonal graphene flake with either armchair or zigzag edges are investigated using molecular dynamics simulations and tight-binding calculations. We found that (i) the pseudomagnetic field in strained graphene flakes is not uniform neither in the center nor at the edge of zigzag terminated flakes, (ii) the pseudomagnetic field is almost zero in the center of armchair terminated flakes but increases dramatically near the edges, (iii) the pseudomagnetic field increases linearly with strain, for strains lower than 15% but increases nonlinearly beyond it, (iv) the local density of states in the center of the zigzag hexagon exhibits pseudo-Landau levels with broken sublattice symmetry in the zeroth pseudo-Landau level, and in addition there is a shift in the Dirac cone due to strain induced scalar potentials, and (v) there is size effect in pseudomagnetic field. This study provides a realistic model of the electronic properties of inhomogeneously strained graphene where the relaxation of the atomic positions is correctly included together with strain induced modifications of the hopping terms up to next-nearest neighbors.

Abstract: The electronic properties of graphene nano-flakes (GNFs) with different edge passivation are investigated by using density functional theory. Passivation with F and H atoms is considered: C-Nc X-Nx (X = F or H). We studied GNFs with 10 < N-c < 56 and limit ourselves to the lowest energy configurations. We found that: (i) the energy difference Delta between the highest occupied molecular orbital and the lowest unoccupied molecular orbital decreases with N-c, (ii) topological defects (pentagon and heptagon) break the symmetry of the GNFs and enhance the electric polarization, (iii) the mutual interaction of bilayer GNFs can be understood by dipole-dipole interaction which were found sensitive to the relative orientation of the GNFs, (iv) the permanent dipoles depend on the edge terminated atom, while the energy gap is independent of it, and (v) the presence of heptagon and pentagon defects in the GNFs results in the largest difference between the energy of the spin-up and spin-down electrons which is larger for the H-passivated GNFs as compared to F-passivated GNFs. Our study shows clearly the effect of geometry, size, termination, and bilayer on the electronic properties of small GNFs. This study reveals important features of graphene nano-flakes which can be detected using Raman spectroscopy. (C) 2014 AIP Publishing LLC.

Abstract: The electronic properties of the two-dimensional material silicene are strongly influenced by the application of a perpendicular electric field E-z and of an exchange field M due to adatoms positioned on the surface or a ferromagnetic substrate. Within the random phase approximation, we investigate how electron-electron interactions are affected by these fields and present analytical and numerical results for the dispersion of plasmons, their lifetime, and their oscillator strength. We find that the combination of the fields E-z and M brings a spin and valley texture to the particle-hole excitation spectrum and allows the formation of spin-and valley-polarized plasmons. When the Fermi level lies in the gap of one spin in one valley, the intraband region of the corresponding spectrum disappears. For zero E-z and finite M the spin symmetry is broken and spin polarization is possible. The lifetime and oscillator strength of the plasmons are shown to depend strongly on the number of spin and valley type electrons that form the electron-hole pairs.

Abstract: The electronic spectrum of three coupled graphene layers (graphene trilayers) is investigated in the presence of an external magnetic field. We obtain analytical expressions for the Landau level spectrum for both the ABA and ABC type of stacking, which exhibit very different dependence on the magnetic field. We show that layer asymmetry and an external gate voltage can strongly influence the properties of the system.