No Arabic abstract
Wannier tight-binding models are effective models constructed from first-principles calculations. As such, they bridge a gap between the accuracy of first-principles calculations and the computational simplicity of effective models. In this work, we extend the existing methodology of creating Wannier tight-binding models from first-principles calculations by introducing the symmetrization post-processing step, which enables the production of Wannier-like models that respect the symmetries of the considered crystal. Furthermore, we implement automatic workflows, which allow for producing a large number of tight-binding models for large classes of chemically and structurally similar compounds, or materials subject to external influence such as strain. As a particular illustration, these workflows are applied to strained III-V semiconductor materials. These results can be used for further study of topological phase transitions in III-V quantum wells.
An empirical multiorbital (spd) tight binding (TB) model including magnetism and spin-orbit coupling is applied to calculations of magnetic anisotropy energy (MAE) in CoPt L1_0 structure. A realistic Slater-Koster parametrisation for single-element transition metals is adapted for the ordered binary alloy. Spin magnetic moment and density of states are calculated using a full-potential linearized augmented plane-wave (LAPW) ab initio method and our TB code with different variants of the interatomic parameters. Detailed mutual comparison of this data allows for determination of a subset of the compound TB parameters tuning of which improves the agreement of the TB and LAPW results. MAE calculated as a function of band filling using the refined parameters is in broad agreement with ab initio data for all valence states and in quantitative agreement with ab initio and experimental data for the natural band filling. Our work provides a practical basis for further studies of relativistic magnetotransport anisotropies by means of local Greens function formalism which is directly compatible with our TB approach.
The topological order of single-crystal Bi and its surface states on the (111) surface are studied in detail based on empirical tight-binding (TB) calculations. New TB parameters are presented that are used to calculate the surface states of semi-infinite single-crystal Bi(111), which agree with the experimental angle-resolved photoelectron spectroscopy results. The influence of the crystal lattice distortion is surveyed and a topological phase transition is found that is driven by in-plane expansion. In contrast with the semi-infinite system, the surface-state dispersions on finite-thickness slabs are non-trivial irrespective of the bulk topological order. The role of the interaction between the top and bottom surfaces in the slab is systematically studied, and it is revealed that a very thick slab is required to properly obtain the bulk topological order of Bi from the (111) surface state: above 150 biatomic layers in this case.
We study the different ways of introducing light-matter interaction in first-principle tight-binding (TB) models. The standard way of describing optical properties is the velocity gauge, defined by linear coupling to the vector potential. In finite systems a transformation to represent the electromagnetic radiation by the electric field instead is possible, albeit subtleties arise in periodic systems. The resulting dipole gauge is a multi-orbital generalization of Peierls substitution. In this work, we investigate accuracy of both pathways, with particular emphasis on gauge invariance, for TB models constructed from maximally localized Wannier functions. Focusing on paradigmatic two-dimensional materials, we construct first-principle models and calculate the response to electromagnetic fields in linear response and for strong excitations. Benchmarks against fully converged first-principle calculations allow for ascertaining the accuracy of the TB models. We find that the dipole gauge provides a more accurate description than the velocity gauge in all cases. The main deficiency of the velocity gauge is an imperfect cancellation of paramagnetic and diamagnetic current. Formulating a corresponding sum rule however provides a way to explicitly enforce this cancellation. This procedure corrects the TB models in the velocity gauge, yielding excellent agreement with dipole gauge and thus gauge invariance.
Maximally localized Wannier functions are localized orthogonal functions that can accurately represent given Bloch eigenstates of a periodic system at a low computational cost, thanks to the small size of each orbital. Tight-binding models based on the maximally localized Wannier functions obtained from different systems are often combined to construct tight-binding models for large systems such as a semi-infinite surface. However, the corresponding maximally localized Wannier functions in the overlapping region of different systems are not identical, and this discrepancy can introduce serious artifacts to the combined tight-binding model. Here, we propose two methods to seamlessly stitch two different tight-binding models that share some basis functions in common. First, we introduce a simple and efficient method: (i) finding the best matching maximally localized Wannier function pairs in the overlapping region belonging to the two tight-binding models, (ii) rotating the spin orientations of the two corresponding Wannier functions to make them parallel to each other, and (iii) making their overall phases equal. Second, we propose a more accurate and generally applicable method based on the iterative minimization of the difference between the Hamiltonian matrix elements in the overlapping region. We demonstrate our methods by applying them to the surfaces of diamond, GeTe, Bi$_2$Se$_3$, and TaAs. Our methods can be readily used to construct reliable tight-binding models for surfaces, interfaces, and defects.
{it Ab initio} investigations of the full static dielectric response and Born effective charge of BN nanotubes (BN-NTs) have been performed for the first time using finite electric field method. It is found that the ionic contribution to the static dielectric response of BN-NTs is substantial and also that a pronounced chirality-dependent oscillation is superimposed on the otherwise linear relation between the longitudinal electric polarizability and the tube diameter ($D$), as for a thin dielectric cylinderical shell. In contrast, the transverse dielectric response of the BN-NTs resemble the behavior of a thin (non-ideal) conducting cylindrical shell of a diameter of $D+4$AA$ $, with a screening factor of 2 for the inner electric field. The medium principal component $Z_y^*$ of the Born effective charge corresponding to the transverse atomic displacement tangential to the BN-NT surface, has a pronounced $D$-dependence (but independent of chirality), while the large longitudinal component $Z_z^*$ exhibits a clear chirality dependence (but nearly $D$-independent), suggesting a powerful way to characterize the diameter and chirality of a BN-NT.