Abstract: When the charge Z of an atom exceeds the critical value of 170, it will undergo a process called atomic collapse which triggers the spontaneous creation of electron-positron pairs. The high charge requirements have prevented the observation of this phenomenon with real atomic nuclei. However, thanks to the relativistic nature of the carriers in graphene, the same physics is accessible at a much lower scale. The atomic collapse analogue in graphene is realized using artificial nuclei which can be created via the deposition of impurities on the surface of graphene or using charged vacancies. These supercritically charged artificial nuclei trap electrons in a sequence of quasi-bound states which can be observed experimentally as resonances in the local density of states.

Abstract: We study ballistic electron transport through a finite chain of quantum circular rings in the presence of spin-orbit interaction (SOI) of strength alpha. The transmission and reflection coefficients for a single ring, obtained analytical lylead to the conductance for a chain of rings as a function of alpha and of the wave vector k of the incident electron. Due to destructive spin interferences the chain can be totaly opaque for certain ranges of k the width of which depends on the value of alpha. A periodic modulation of a widens up the gaps considerably and produces a nearly binary conductance output.

Abstract: Using improved Wigner-Brillouin perturbation theory we study resonant electron-phonon interaction in a semiconductor quantum dot. We predict pinning of the excited energy levels to the ground state level plus one optical phonon as a function of the strength of the confinement potential. This effect should be observable through optical spectroscopic measurements.

Abstract: Using a numerical tight-binding approach based on the Chebyshev–Bogoliubov–de Gennes method we describe Josephson junctions made of multilayer graphene contacted by top superconducting gates. Both Bernal (ABA) and rhombohedral (ABC) stacking are considered and we find that the type of stacking has a strong effect on the proximity effect and the supercurrent flow. For both cases the pair amplitude shows a polarization between dimer and nondimer atoms, being more pronounced for rhombohedral stacking. Even though the proximity effect in nondimer sites is enhanced when compared to single-layer graphene, we find that the supercurrent is suppressed. The spatial distribution of the supercurrent shows that for Bernal stacking the current flows only in the topmost layers while for rhombohedral stacking the current flows throughout the whole structure.

Abstract: Using highly efficient GPU-based simulations of the tight-binding Bogoliubov-de Gennes equations we solve self-consistently for the pair correlation in rhombohedral (ABC) and Bernal (ABA) multilayer graphene by considering a finite intrinsic s-wave pairing potential. We find that the two different stacking configurations have opposite bulk/surface behavior for the order parameter. Surface superconductivity is robust for ABC stacked multilayer graphene even at very low pairing potentials for which the bulk order parameter vanishes, in agreement with a recent analytical approach. In contrast, for Bernal stacked multilayer graphene, we find that the order parameter is always suppressed at the surface and that there exists a critical value for the pairing potential below which no superconducting order is achieved. We considered different doping scenarios and find that homogeneous doping strongly suppresses surface superconductivity while nonhomogeneous field-induced doping has a much weaker effect on the superconducting order parameter. For multilayer structures with hybrid stacking (ABC and ABA) we find that when the thickness of each region is small (few layers), high-temperature surface superconductivity survives throughout the bulk due to the proximity effect between ABC/ABA interfaces where the order parameter is enhanced. DOI: 10.1103/PhysRevB.87.134509

Abstract: Starting from the simplified linear combination of atomic orbitals method in combination with first-principles calculations, we construct a tight-binding (TB) model in the two-centre approximation for borophene and hydrogenated borophene (borophane). The Slater and Koster approach is applied to calculate the TB Hamiltonian of these systems. We obtain expressions for the Hamiltonian and overlap matrix elements between different orbitals for the different atoms and present the SK coefficients in a nonorthogonal basis set. An anisotropic Dirac cone is found in the band structure of borophane. We derive a Dirac low-energy Hamiltonian and compare the Fermi velocities with that of graphene.

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: Using the simplified linear combination of atomic orbitals (LCAO) method in combination with ab initio calculations, we construct a tight-binding (TB) model for two different crystal structures of monolayer gallium: a(100)- and b(010)-Gallenene. The analytical expression for the Hamiltonian and numerical results for the overlap matrix elements between different orbitals of the Ga atoms and for the Slater and Koster (SK) integrals are obtained. We find that the compaction of different structures affects significantly the formation of the orbitals. The results for a(100)-Gallenene can be very well explained with an orthogonal basis set, while for b(010)-Gallenene we have to assume a non-orthogonal basis set in order to construct the TB model. Moreover, the transmission properties of nanoribbons of both monolayers oriented along the AC and ZZ directions are also investigated and it is shown that both AC- and ZZ-b(010)-Gallenene nanoribbons exhibit semiconducting behavior with zero transmission while those of a(100)-Gallenene nanoribbons are metallic.

Abstract: Using an extension of the Gutzwiller approximation for an inhomogeneous system, we study the two-band Hubbard model with unequal band widths for a slab geometry. The aim is to investigate the mutual effect of individual bands on the spatial distribution of quasi-particle weight and charge density, especially near the surface of the slab. The main effect of the difference in band width is the presence of two different length scales corresponding to the quasi-particle profile of each band. This is enhanced in the vicinity of the critical interaction of the narrow band where an orbitally selective Mott transition occurs and a surface dead layer forms for the narrow band. For the doped case, two different regimes of charge transfer between the surface and the bulk of the slab are revealed. The charge transfer from surface/ center to center/ surface depends on both the doping level and the average relative charge accumulated in each band. Such effects could also be of importance when describing the accumulation of charges at the interface between structures made of multi-band strongly correlated materials.

Abstract: The influence of triaxial in-plane strain on the electronic properties of a hexagonal boron-nitride sheet is investigated using density functional theory. Different from graphene, the triaxial strain localizes the molecular orbitals of the boron-nitride flake in its center depending on the direction of the applied strain. The proposed technique for localizing the molecular orbitals that are close to the Fermi level in the center of boron nitride flakes can be used to actualize engineered nanosensors, for instance, to selectively detect gas molecules. We show that the central part of the strained flake adsorbs polar molecules more strongly as compared with an unstrained sheet.

Abstract: Recent experiments showed that nonuniform strain can be produced by depositing graphene over pillars. We employed atomistic calculations to study the nonuniform strain and the induced pseudomagnetic field in graphene on top of nanopillars. By decreasing the distance between the nanopillars a complex distribution for the pseudomagnetic field can be generated. Furthermore, we performed tight-binding calculations of the local density of states (LDOS) by using the relaxed graphene configuration obtained from atomistic calculations. We find that the quasiparticle LDOS are strongly modified near the pillars, both at low energies showing sublattice polarization and at high energies showing shifts of the van Hove singularity. Our study shows that changing the specific pattern of the nanopillars allows us to create a desired shape of the pseudomagnetic field profile while the LDOS maps provide an input for experimental verification by scanning tunneling microscopy.

Abstract: We investigate the stability of circular monolayer graphene subjected to a radial load using non-equilibrium molecular dynamics simulations. When monolayer graphene is radially stressed, after some small circular strain (~0.4%) it buckles and bends into a new bowl-like shape. Young's modulus is calculated from the linear relation between stress and strain before the buckling threshold, which is in agreement with experimental results. The prediction of elasticity theory for the buckling threshold of a radially stressed plate is presented and its results are compared to the one of our atomistic simulation. The Jarzynski equality is used to estimate the difference between the free energy of the non-compressed states and the buckled states. From a calculation of the free energy we obtain the optimum radius for which the system feels the minimum boundary stress.