Abstract: Using the four-band continuum model we derive a general expression for the infinite-mass boundary condition in bilayer graphene. Applying this new boundary condition we analytically calculate the confined states and the corresponding wave functions in a bilayer graphene quantum dot in the absence and presence of a perpendicular magnetic field. Our results for the energy spectrum show an energy gap between the electron and hole states at small magnetic fields. Furthermore the electron (e) and hole (h) energy levels corresponding to the K and K' valleys exhibit the E-K(e(h)) (m) = E-K'(e(h)) (m) symmetry, where m is the angular momentum quantum number. (C) 2014 Elsevier Ltd. All rights reserved.

Abstract: We present the Tight-Binding Studio (TB Studio) software package that calculates the different parameters of a tight-binding Hamiltonian from a set of Bloch energy bands obtained from first principle theories such as density functional theory, Hartree-Fock calculations or semi-empirical band-structure theory. This will be helpful for scientists who are interested in studying electronic and optical properties of structures using Green's function theory within the tight-binding approximation. TB Studio is a cross-platform application written in C++ with a graphical user interface design that is user-friendly and easy to work with. This software is powered by Linear Algebra Package C interface library for solving the eigenvalue problems and the standard high performance OpenGL graphic library for real time plotting. TB Studio and its examples together with the tutorials are available for download from tight-binding.com. Program summary Program Title: Tight-Binding Studio Program Files doi:http://dx.doi.org/10.17632/j6x5mwzm2d.1 Licensing provisions: LGPL Programming language: C++ External routines: BLAS, LAPACK, LAPACKE, wxWidgets, OpenGL, MathGL Nature of problem: Obtaining Tight-Binding Hamiltonian from a set of Bloch energy bands obtained from first-principles calculations. Solution method: Starting from the simplified LCAO method, a tight-binding model in the two-center approximation is constructed. The Slater and Koster (SK) approach is used to calculate the parameters of the TB Hamiltonian. By using non-linear fitting approaches the optimal values of the SK parameters are obtained such that the TB energy eigenvalues are as close as possible to those from first-principles calculations. We obtain the expression for the Hamiltonian and the overlap matrix elements between the different orbitals of the different atoms in an orthogonal or non-orthogonal basis set. (C) 2020 Elsevier B.V. All rights reserved.

Abstract: Dibismuth dioxychalcogenides, Bi2O2Ch (Ch = S, Se, Te), are a promising class of materials for next-generation electronics and thermoelectrics due to their ultrahigh carrier mobility and excellent air stability. An interesting member of this family is Bi2O2S, which has a stereochemically active 6s2 lone pair of Bi3+ cations, heterogeneous bonding, and a high mass contrast between its constituent elements. In the present study, we have used first-principles calculations in combination with Boltzmann transport theory to systematically investigate the effect of hydrostatic pressure on lattice dynamics and phonon transport properties of Bi2O2S. We found that the ambient Pnmn phase has a low average lattice thermal conductivity (kappa l) of 1.71 W/(m K) at 300 K. We also predicted that Bi2O2S undergoes a structural phase transition from a low-symmetry (Pnmn) to a high-symmetry (I4/mmm) structure at around 4 GPa due to centering of Bi3+ cations with pressure. Upon compression, the lone pair activity of Bi3+ cations is suppressed, which increases kappa l by almost 3 times to 4.92 W/ (m K) at 5 GPa for the I4/mmm phase. The computed phonon lifetimes and Gru''neisen parameters show that anharmonicity decreases with increasing pressure due to further suppression of the lone pair activity and strengthening of intra-and intermolecular interactions, leading to an average room-temperature kappa l of 12.82 W/(m K) at 20 GPa. Overall, this study provides a comprehensive understanding of the effect of hydrostatic pressure on the stereochemical activity of the lone pair of Bi3+ cations and its implications on the phonon transport properties of Bi2O2S.

Abstract: By means of density functional theory-based first-principle calculations, the structural, vibrational and electronic properties of single-layer Ge3N4 are investigated. Structural optimizations and phonon band dispersions reveal that single-layer ultrathin form of Ge3N4 possesses a dynamically stable buckled structure with large hexagonal holes. Predicted Raman spectrum of single-layer Ge3N4 indicates that the buckled holey structure of the material exhibits distinctive vibrational features. Electronic band dispersion calculations indicate the indirect band gap semiconducting nature of single-layer Ge3N4. It is also proposed that single-layer Ge3N4 forms type-II vertical heterostructures with various planar and puckered 2D materials except for single-layer GeSe which gives rise to a type-I band alignment. Moreover, the electronic properties of single-layer Ge3N4 are investigated under applied external in-plane strain. It is shown that while the indirect gap behavior of Ge3N4 is unchanged by the applied strain, the energy band gap increases (decreases) with tensile (compressive) strain.

Abstract: Molecular dynamics simulations are used to investigate the effect of an AFM tip when indenting graphene nano bubbles filled by a noble gas (i.e. He, Ne and Ar) up to the breaking point. The failure points resemble those of viral shells as described by the Foppl-von Karman (FvK) dimensionless number defined in the context of elasticity theory of thin shells. At room temperature, He gas inside the bubbles is found to be in the liquid state while Ne and Ar atoms are in the solid state although the pressure inside the nano bubble is below the melting pressure of the bulk. The trapped gases are under higher hydrostatic pressure at low temperatures than at room temperature.

Abstract: In this study, reactivity of single-layer GaAs against Ge atoms is studied by means of ab initio density functional theory calculations. Firstly, it is shown that Ge atoms interact quite strongly with the GaAs layer which allows the formation of Ge islands while it hinders the growth of detached germanene monolayers. It is also predicted that adsorption of Ge atoms on GaAs single-layer lead to formation of two novel stable single-layer crystal structures, namely 1H-GaGeAs and 1H(A)-GaGeAs. Both the total energy optimizations and the calculated vibrational spectra indicate the dynamical stability of both single layer structures. Moreover, although both structures crystallize in 1H phase, 1H-GaGeAs and 1H(A)-GaGeAs exhibit distinctive vibrational features in their Raman spectra which is quite important for distinguishing the structures. In contrast to the semiconducting nature of single-layer GaAs, both polytypes of GaGeAs exhibit metallic behavior confirmed by the electronic band dispersions. Furthermore, the linear-elastic constants, in-plane stiffness and Poisson ratio, reveal the ultrasoft nature of the GaAs and GaGeAs structures and the rigidity of GaAs is found to be slightly enhanced via Ge adsorption. With their stable, ultra-thin and metallic properties, predicted single-layer GaGeAs structures can be promising candidates for nanoscale electronic and mechanical applications.

Abstract: Motivated by the recent synthesis of two-dimensional VI3 [Kong et al. Adv. Mater. 31, 1808074 (2019)], we investigate the effect of V doping on the magnetic and electronic properties of monolayer VI3 by means of first-principles calculations. The dynamically stable semiconducting ferromagnetic (FM) and antiferromagnetic (AFM) phases of monolayer VI3 are found to display distinctive vibrational features that the magnetic state can be distinguished by Raman spectroscopy. In order to clarify the effect of experimentally observed excessive V atoms, the magnetic and electronic properties of the V-doped VI3 structures are analyzed. Our findings indicate that partially doped VI3 structures display FM ground state while the fully-doped structure exhibits AFM ground state. The fully-doped monolayer VI3 is found to be a semiconductor with a relatively larger band gap than its pristine structure. In addition, strain-dependent electronic and magnetic properties of fully- and partially-doped VI3 structures reveal that pristine monolayer displays a FM-to-AFM phase transition with robust semiconducting nature for 5% of compressive strain, while fully-doped monolayer VI3 structure possesses AFM-to-FM semiconducting transition at tensile strains larger than 4%. In contrast, the partially-doped VI3 monolayers are found to display robust FM ground state under biaxial strain. Its dopant and strain tunable electronic and magnetic nature makes monolayer VI3 a promising material for applications in nanoscale spintronic devices.

Abstract: MXenes are an emerging class of two-dimensional materials, which in their thinnest limit consist of a monolayer of carbon or nitrogen (X) sandwiched between two transition metal (M) layers. We have systematically searched for superconductivity among MXenes for a range of transition metal elements, based on a full first-principles characterization in combination with the Eliashberg formalism. Thus, we identified six superconducting MXenes: three carbides (Mo2C, W2C and Sc2C) and three nitrides (Mo2N, W2N and Ta2N). The highest critical temperature of similar to 16 K is found in Mo2N, for which a successful synthesis method has been established [Urbankowskiet al.,Nanoscale, 2017,9, 17722-17730]. Moreover, W2N presents a novel case of competing superconducting and charge density wave phases.

Abstract: We explore the flatness of conduction and valence bands of interlayer excitons in MoS2/WSe2 van der Waals heterobilayers, tuned by interlayer twist angle, pressure, and external electric field. We employ an efficient continuum model where the moire pattern from lattice mismatch and/or twisting is represented by an equivalent mesoscopic periodic potential. We demonstrate that the mismatch moire potential is too weak to produce significant flattening. Moreover, we draw attention to the fact that the quasi-particle effective masses around the Gamma-point and the band flattening are reduced with twisting. As an alternative approach, we show (i) that reducing the interlayer distance by uniform vertical pressure can significantly increase the effective mass of the moire hole, and (ii) that the moire depth and its band flattening effects are strongly enhanced by accessible electric gating fields perpendicular to the heterobilayer, with resulting electron and hole effective masses increased by more than an order of magnitude – leading to record-flat bands. These findings impose boundaries on the commonly generalized benefits of moire twistronics, while also revealing alternative feasible routes to achieve truly flat electron and hole bands to carry us to strongly correlated excitonic phenomena on demand.

Abstract: We detail the effects of Cl and S functionalization on the superconducting properties of layered (bulk) and monolayer niobium carbide (Nb2C) MXene crystals, based on first-principles calculations combined with Eliashberg theory. For bulk layered Nb2CCl2, the calculated superconducting transition temperature (T-c) is in very good agreement with the recently measured value of 6 K. We show that T-c is enhanced to 10 K for monolayer Nb2CCl2, due to an increase in the density of states at the Fermi level, and the corresponding electron-phonon coupling. We further demonstrate feasible gate- and strain-induced enhancements of T-c for both bulk-layered and monolayer Nb2CCl2 crystals, resulting in T-c values of around 38 K. In the S-functionalized Nb2CCl2 crystals, our calculations reveal the importance of phonon softening in understanding their superconducting properties. Finally, we predict that Nb3C2S2 in bulk-layered and monolayer forms is also superconducting, with a T-c of around 28 K. Considering that Nb2C is not superconducting in pristine form, our findings promote functionalization as a pathway towards robust superconductivity in MXenes.

Abstract: We present extensive calculations of the optical and plasmonic properties of a graphene sheet carrying a dc current. By calculating analytically the density-density response function of current-carrying states at finite temperature, we demonstrate that an applied dc current modifies the Pauli blocking mechanism and that absorption acquires a birefringent character with respect to the angle between the in-plane light polarization and current flow. Employing the random phase approximation at finite temperature, we show that graphene plasmons display a degree of non-reciprocity and collimation that can be tuned with the applied current. We discuss the possibility to measure these effects.

Abstract: We study the electronic properties of commensurate faulted bilayer graphene by diagonalizing the one-particle Hamiltonian of the bilayer system in a complete basis of Bloch states of the individual graphene layers. Our novel approach is very general and can be easily extended to any commensurate graphene-based heterostructure. Here, we consider three cases: (i) twisted bilayer graphene, (ii) bilayer graphene where triaxial stress is applied to one layer and (iii) bilayer graphene where uniaxial stress is applied to one layer. We show that the resulting superstructures can be divided into distinct classes, depending on the twist angle or the magnitude of the induced strain. The different classes are distinguished from each other by the interlayer coupling mechanism, resulting in fundamentally different low-energy physics. For the cases of triaxial and uniaxial stress, the individual graphene layers tend to decouple and we find significant charge transfer between the layers. In addition, this piezoelectric effect can be tuned by applying a perpendicular electric field. Finally, we show how our approach can be generalized to multilayer systems.

Abstract: <script type='text/javascript'>document.write(unpmarked('Quantum electrodynamics predicts that heavy atoms (Z \u003E Z(c) approximate to 170) will undergo the process of atomic collapse where electrons sink into the positron continuum and a new family of so-called collapsing states emerges. The relativistic electrons in graphene exhibit the same physics but at a much lower critical charge (Z(c) approximate to 1) which has made it possible to confirm this phenomenon experimentally. However, there exist conflicting predictions on the effect of a magnetic field on atomic collapse. These theoretical predictions are based on the continuum Dirac-Weyl equation, which does not have an exact analytical solution for the interplay of a supercritical Coulomb potential and the magnetic field. Approximative solutions have been proposed, but because the two effects compete on similar energy scales, the theoretical treatment varies depending on the regime which is being considered. These limitations are overcome here by starting from a tight-binding approach and computing exact numerical results. By avoiding special limit cases, we found a smooth evolution between the different regimes. We predict that the atomic collapse effect persists even after the magnetic field is activated and that the critical charge remains unchanged. We show that the atomic collapse regime is characterized: (1) by a series of Landau level anticrossings and (2) by the absence of root B scaling of the Landau levels with regard to magnetic field strength.'));

Abstract: Using the non-equilibrium Green's function formalism, we study carrier transport through imperfect two-dimensional (2D) topological insulator (TI) ribbons. In particular, we investigate the effect of vacancy defects on the carrier transport in 2D TI ribbons with hexagonal lattice structure. To account for the random distribution of the vacancy defects, we present a statistical study of varying defect densities by stochastically sampling different defect configurations. We demonstrate that the topological edge states of TI ribbons are fairly robust against a high concentration (up to 2%) of defects. At very high defect densities, we observe an increased inter-edge interaction, mediated by the localisation of the edge states within the bulk region. This effect causes significant back-scattering of the, otherwise protected, edge-states at very high defect concentrations (>2%), resulting in a loss of conduction through the TI ribbon. We discuss how this coherent vacancy scattering can be used to our advantage for the development of TI-based transistors. We find that there is an optimal concentration of vacancies yielding an ON-OFF current ratio of up to two orders of magnitude. Finally, we investigate the importance of spin-orbit coupling on the robustness of the edge states in the TI ribbon and show that increased spin-orbit coupling could further increase the ON-OFF ratio.

Abstract: Atomic collapse is a phenomenon inherent to relativistic quantum mechanics where electron states dive in the positron continuum for highly charged nuclei. This phenomenon was recently observed in graphene. Here we investigate a novel collapse phenomenon when multiple sub- and supercritical charges of equal strength are put close together as in a molecule. We construct a phase diagram which consists of three distinct regions: (1) subcritical, (2) frustrated atomic collapse, and (3) molecular collapse. We show that the single impurity atomic collapse resonances rearrange themselves to form molecular collapse resonances which exhibit a distinct bonding, anti-bonding and non-bonding character. Here we limit ourselves to systems consisting of two and three charges. We show that by tuning the distance between the charges and their strength a high degree of control over the molecular collapse resonances can be achieved.