Abstract: Optical and electronic properties of black phosphorus strongly depend on the number of layers and type of stacking. Using first-principles calculations within the framework of density functional theory, we investigate the electronic properties of bilayer black phosphorus with an interlayer twist angle of 90 degrees. These calculations are complemented with a simple (k) over right arrow . (p) over right arrow model which is able to capture most of the low energy features and is valid for arbitrary twist angles. The electronic spectrum of 90 degrees twisted bilayer black phosphorus is found to be x-y isotropic in contrast to the monolayer. However x-y anisotropy, and a partial return to monolayer-like behavior, particularly in the valence band, can be induced by an external out-of-plane electric field. Moreover, the preferred hole effective mass can be rotated by 90 degrees simply by changing the direction of the applied electric field. In particular, a +0.4 (-0.4) V angstrom(1) out-of-plane electric field results in a similar to 60% increase in the hole effective mass along the y (x) axis and enhances the m(y)*/m(x)* (m(x)*/m(y)*) ratio as much as by a factor of 40. Our DFT and (k) over right arrow . (p) over right arrow simulations clearly indicate that the twist angle in combination with an appropriate gate voltage is a novel way to tune the electronic and optical properties of bilayer phosphorus and it gives us a new degree of freedom to engineer the properties of black phosphorus based devices.

Abstract: The self rotation of a graphene flake over graphite is controlled by the size, initial misalignment and temperature. Using both ab initio calculations and molecular dynamics simulations, we investigate annealing effects on the self rotation of a graphene flake on a graphene substrate. The energy barriers for rotation and drift of a graphene flake over graphene is found to be smaller than 25 meV/atom which is comparable to thermal energy. We found that small flakes (of about similar to 4 nm) are more sensitive to temperature and initial misorientation angles than larger one (beyond 10 nm). The initial stacking configuration of the flake is found to be important for its dynamics and time evolution of misalignment. Large flakes, which are initially in the AA-or AB-stacking state with small misorientation angle, rotate and end up in the AB-stacking configuration. However small flakes can they stay in an incommensurate state specially when the initial misorientation angle is larger than 2 degrees. Our results are in agreement with recent experiments.

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.

Abstract: Dirac plasmons in graphene hybridize with phonons of transition metal dichalcogenides (TMDs) when the materials are combined in so-called van der Waals heterostructures (vdWh), thus forming surface plasmon-phonon polaritons (SPPPs). The extend to which these modes are coupled depends on the TMD composition and structure, but also on the plasmons' properties. By performing realistic simulations that account for the contribution of each layer of the vdWh separately, we calculate how the strength of plasmon-phonon coupling depends on the number and composition of TMD layers, on the graphene Fermi energy and the specific phonon mode. From this, we present a semiclassical theory that is capable of capturing all relevant characteristics of the SPPPs. We find that it is possible to realize both strong and ultra-strong coupling regimes by tuning graphene's Fermi energy and changing TMD layer number.

Abstract: Topological semimetal phases in two-dimensional (2D) materials have gained widespread interest due to their potential applications in novel nanoscale devices. Despite the growing number of studies on 2D topological nodal lines (NLs), candidates with significant topological features that combine nontrivial topological semimetal phase with superconductivity are still rare. Herein, we predict Al2B2 and AlB4 monolayers as new 2D nonmagnetic Dirac nodal line semimetals with several novel features. Our extensive electronic structure calculations combined with analytical studies reveal that, in addition to multiple Dirac points, these 2D configurations host various highly dispersed NLs around the Fermi level, all of which are semimetal states protected by time-reversal and in-plane mirror symmetries. The most intriguing NL in Al2B2 encloses the K point and crosses the Fermi level, showing a considerable dispersion and thus providing a fresh playground to explore exotic properties in dispersive Dirac nodal lines. More strikingly, for the AlB4 monolayer, we provide the first evidence for a set of 2D nonmagnetic open type-II NLs coexisting with superconductivity at a rather high transition temperature. The coexistence of superconductivity and nontrivial band topology in AlB4 not only makes it a promising material to exhibit novel topological superconducting phases, but also a rather large energy dispersion of type-II nodal lines in this configuration may offer a platform for the realization of novel topological features in the 2D limit.

Abstract: Silicon-based materials are attractive for thermoelectric applications due to their thermal stability, chemical inertness, and natural abundance of silicon. Here, using a combination of first-principles and Boltzmann transport calculations we report the thermoelectric properties of the recently synthesized compound Mg3Si2Te6. Our analysis reveals that Mg3Si2Te6 is a direct bandgap semiconductor with a bandgap of 1.6 eV. The combination of heavy and light valence bands, along with a high valley degeneracy, results in a large power factor under p-type doping. We also find that Mg is weakly bonded both within and between the layers, leading to low phonon group velocities. The vibrations of the Mg atoms are localized and make a significant contribution to phonon-phonon scattering. This high anharmonicity, coupled with low phonon group velocity, results in a low lattice thermal conductivity of & kappa;(l) = 0.5 W m(-1) K-1 at room temperature, along the cross-plane direction. Combining excellent electronic transport properties and low & kappa;(l), p-type Mg3Si2Te6 achieves figure-of-merit (zT) values greater than 1 at temperatures above 600 K. Specifically, a zT of 2.0 is found at 900 K along the cross-plane direction. Our findings highlight the importance of structural complexity and chemical bonding in electronic and phonon transport, providing guiding insights for further design of Si-based thermoelectrics.

Abstract: Epitaxial graphene is grown on a non-polar n(+) 6H-SiC m-plane substrate and studied using atomic scale scanning tunneling microscopy. Multilayer graphene is found throughout the surface and exhibits rotational disorder. Moire patterns of different spatial periodicities are found, and we found that as the wavelength increases, so does the amplitude of the modulations. This relationship reveals information about the interplay between the energy required to bend graphene and the interaction energy, i.e. van der Waals energy, with the graphene layer below. Our experiments are supported by theoretical calculations which predict that the membrane topographical amplitude scales with the Moire pattern wavelength, L as L-1 + alpha L-2. (C) 2014 Elsevier Ltd. All rights reserved.

Abstract: Two-dimensional (2D) carbon nitride materials play an important role in energy-harvesting, energy-storage and environmental applications. Recently, a new carbon nitride, 2D polyaniline (C3N) was proposed [PNAS 113 (2016) 7414-7419]. Based on the structure model of this C3N monolayer, we propose two new carbon nitride monolayers, named dumbbell (DB) C4N-I and C4N-II. Using first-principles calculations, we systematically study the structure, stability, and band structure of these two materials. In contrast to other carbon nitride monolayers, the orbital hybridization of the C/N atoms in the DB C4N monolayers is sp(3). Remarkably, the band structures of the two DB C4N monolayers have a Dirac cone at the K point and their Fermi velocities (2.6/2.4 x 10(5) m/s) are comparable to that of graphene. This makes them promising materials for applications in high-speed electronic devices. Using a tight-binding model, we explain the origin of the Dirac cone. (C) 2017 Elsevier Ltd. All rights reserved.

Abstract: The flexoelectric and electronic properties of zig-zag graphene nanoribbons are explored under mechanical bending using state of the art first principles calculations. A linear dependence of the bending induced out of plane polarization on the applied strain gradient is found. The inferior flexoelectric properties of graphene nanoribbons can be improved by more than two orders of magnitude by hydrogen and fluorine functionalization (CH and CF nanoribbons). A large out of plane flexoelectric effect is predicted for CF nanoribbons. The origin of this enhancement lies in the electro-negativity difference between carbon and fluorine atoms, which breaks the out of plane charge symmetry even for a small strain gradient. The flexoelectric effect can be further improved by co-functionalization with hydrogen and fluorine (CHF Janus-type nanoribbon), where a spontaneous out of plane dipole moment is formed even for flat nanoribbons. We also find that bending can control the charge localization of valence band maxima and therefore enables the tuning of the hole effective masses and band gaps. These results present an important advance towards the understanding of flexoelectric and electronic properties of hydrogen and fluorine functionalized graphene nanoribbons, which can have important implications for flexible electronic applications. (C) 2020 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: Water permeation between stacked layers of hBN sheets forming 2D nanochannels is investigated using large-scale ab initio-quality molecular dynamics simulations. A high-dimensional neural network potential trained on density-functional theory calculations is employed. We simulate water in van der Waals nanocapillaries and study the impact of nanometric confinement on the structure and dynamics of water using both equilibrium and nonequilibrium methods. At an interlayer distance of 10.2 A confinement induces a first-order phase transition resulting in a well-defined AA-stacked bilayer of hexagonal ice. In contrast, for h < 9 A, the 2D water monolayer consists of a mixture of different locally ordered patterns of squares, pentagons, and hexagons. We found a significant change in the transport properties of confined water, particularly for monolayer water where the water-solid friction coefficient decreases to half and the diffusion coefficient increases by a factor of 4 as compared to bulk water. Accordingly, the slip-velocity is found to increase under confinement and we found that the overall permeation is dominated by monolayer water adjacent to the hBN membranes at extreme confinements. We conclude that monolayer water in addition to bilayer ice has a major contribution to water transport through 2D nanochannels.

Abstract: Motivated by the recent synthesis of two-dimensional alpha-Fe2O3 (Balan et al 2018 Nat. Nanotechnol. 13 602), we analyze the structural, vibrational, electronic and magnetic properties of single- and few-layer alpha-Fe2O3 compared to bulk, by ab initio and Monte-Carlo simulations. We reveal how monolayer alpha-Fe2O3 (hematene) can be distinguished from the few-layer structures, and how they all differ from bulk through observable Raman spectra. The optical spectra exhibit gradual shift of the prominent peak to higher energy, as well as additional features at lower energy when alpha-Fe2O3 is thinned down to a monolayer. Both optical and electronic properties have strong spin asymmetry, meaning that lower-energy optical and electronic activities are allowed for the single-spin state. Finally, our considerations of magnetic properties reveal that 2D hematite has anti-ferromagnetic ground state for all thicknesses, but the critical temperature for Morin transition increases with decreasing sample thickness. On all accounts, the link to available experimental data is made, and further measurements are prompted.

Abstract: A twin boundary in one of the layers of a twisted van der Waals heterostructure separates regions with near opposite inter-layer twist angles. In a MoS₂/WSe₂bilayer, the regions with<inline-formula><tex-math><?CDATA $Rh^h$?></tex-math><math overflow=“scroll”><msubsup><mi>R</mi><mi>h</mi><mi>h</mi></msubsup></math><inline-graphic href=“tdmac529dieqn1.gif” type=“simple” /></inline-formula>and<inline-formula><tex-math><?CDATA $Rh^X$?></tex-math><math overflow=“scroll”><msubsup><mi>R</mi><mi>h</mi><mi>X</mi></msubsup></math><inline-graphic href=“tdmac529dieqn2.gif” type=“simple” /></inline-formula>stacking registry that defined the sub-lattices of the moiré honeycomb pattern would be mirror-reflected across such a twist boundary. In that case, we demonstrate that topologically protected chiral moiré exciton states are confined at the twist boundary. These are one-dimensional and uni-directional excitons with opposite velocities for excitons composed by electronic states with opposite valley/spin character, enabling intrinsic, guided, and far reaching valley-polarized exciton currents.

Abstract: Dirac plasmons in a two-dimensional (2D) crystal are strongly affected by the dielectric properties of the environment, due to interaction of their electric field lines with the surrounding medium. Using graphene as a 2D reservoir of free carriers, one can engineer a material configuration that provides an anisotropic environment to the plasmons. In this work, we discuss the physical properties of Dirac plasmons in graphene surrounded by an arbitrary anisotropic dielectric and exemplify how h-BN-based heterostructures can be designed to bear the required anisotropic characteristics. We calculate how dielec-tric anisotropy impacts the spatial propagation of the plasmons and find that an anisotropy-induced plasmon mode emerges, together with a damping pathway, that stem from the out-of-plane off-diagonal elements in the dielectric tensor. Furthermore, we find that one can create hyperbolic plasmons by inher-iting the dielectric hyperbolicity of the designed material environment. Strong control over plasmon propagation patterns can be realized in a similar manner. Finally, we show that in this way one can also control the polarization of the light-matter excitations that constitute the plasmon. Taken together, our results promote the design of the dielectric environment as an effective path to tailor the plasmonic response of graphene on the nanoscopic level.

Abstract: Using density-functional theory in combination with the nonequilibrium Green's function formalism, we study the effect of iodide/chloride and iodide/bromide mixing on the electronic transport in lead based organometallic perovskite CH3NH3PbI3, which is known to be an effective tool to tune the electronic and optical properties of such materials. We found that depending on the level and position of the halide mixing, the electronic transport can be increased by more than a factor of 4 for a given voltage biasing. The largest current is observed for small concentration of bromide substitutions located at the equatorial sites. However, full halide substitution has a negative effect on the transport properties of this material: the current drops by an order of magnitude for both CH3NH3PbCl3 and CH3NH3PbBr3 samples. (C) 2015 Elsevier B.V. All rights reserved.