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 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: The McMillan-Ginzburg-Landau (MGL) model for charge density waves (CDW) is employed in a systematic phenomenological study of the different phases that have been probed in recent experiments involving transition metal dichalcogenides. We implemented an efficient imaginary time evolution method to solve the MGL equations, which enabled us to investigate the role of different coupling parameters on the CDW patterns and to perform calculations with different energy functionals that lead to several experimentally observed singularities in the CDW phase profiles. In particular, by choosing the appropriate energy functionals, we were able to obtain phases that go beyond the well-known periodic phase slips (discommensurations), exhibiting also topological defects (i.e., vortex-antivortex pairs), domain walls where the CDW order parameter is suppressed, and even CDW with broken rotational symmetry. Finally, we briefly discuss the effect of these different CDW phases on the profile and critical temperature of the competing superconducting state.

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: We study the effect of mechanical deformations on the elec- tronic properties of hexagonal flakes of bilayer graphene. The behavior of electrons induced by triaxial strain can be de- scribed by an effective pseudo-magnetic field which is homo- geneous in the center of the flake. We find that in-plane strain, applied to both layers equally, can break the layer symmetry leading to different behavior in the top and bottom layers of graphene. At low energy, just one of the layers feels
the pseudo-magnetic field: the zero-energy pseudo-Landau level is missing in the second layer, thus creating a gap be- tween the lowest non-zero levels. While the layer asymmetry is most significant at zero energy, interaction with the edges of the flake extends the effect to higher pseudo-Landau lev- els. The behavior of the top and bottom layers may be re- versed by rotating the triaxial strain by 60°.

Abstract: The effect of strain in graphene is usually modeled by a pseudomagnetic vector potential which is, however, derived in the limit of small strain. In realistic cases deviations are expected in view of graphene's very high strain tolerance, which can be up to 25%. Here we investigate the pseudomagnetic field generated by a Gaussian bump and we show that it exhibits significant differences with numerical tight-binding results. Furthermore, we calculate the electronic states in the strained region for a hexagon shaped flake with armchair edges. We find that the sixfold symmetry of the wave functions inside the Gaussian bump is directly related to the different effects of strain along the fundamental directions of graphene: zigzag and armchair. Low energy electrons are strongly confined in the armchair directions and are localized on the carbon atoms of a single sublattice.

Abstract: Electrons in graphene, in addition to their spin, have two pseudospin degrees of freedom: sublattice and valley pseudospin. Valleytronics uses the valley degree of freedom as a carrier of information similarly to the way spintronics uses electron spin. We show how a double-barrier structure consisting of electric and vector potentials can be used to filter massive Dirac electrons based on their valley index. We study the resonant transmission through a finite number of barriers and we obtain the energy spectrum of a superlattice consisting of electric and vector potentials. When a mass term is included, the energy bands and energy gaps at the K and K′ points are different and they can be tuned by changing the potential.