No Arabic abstract
For a previously published study of the titanium hcp (alpha) to omega (omega) transformation, a tight-binding model was developed for titanium that accurately reproduces the structural energies and electron eigenvalues from all-electron density-functional calculations. We use a fitting method that matches the correctly symmetrized wavefuctions of the tight-binding model to those of the density-functional calculations at high symmetry points. The structural energies, elastic constants, phonon spectra, and point-defect energies predicted by our tight-binding model agree with density-functional calculations and experiment. In addition, a modification to the functional form is implemented to overcome the collapse problem of tight-binding, necessary for phase transformation studies and molecular dynamics simulations. The accuracy, transferability and efficiency of the model makes it particularly well suited to understanding structural transformations in titanium.
We extend a tight-binding method to include the effects of spin-orbit coupling, and apply it to the study of the electronic properties of the actinide elements Th, U, and Pu. These tight-binding parameters are determined for the fcc crystal structure using the equivalent equilibrium volumes. In terms of the single particle energies and the electronic density of states, the overall quality of the tight-binding representation is excellent and of the same quality as without spin-orbit coupling. The values of the optimized tight-binding spin-orbit coupling parameters are comparable to those determined from purely atomic calculations.
A method for incorporating electromagnetic fields into empirical tight-binding theory is derived from the principle of local gauge symmetry. Gauge invariance is shown to be incompatible with empirical tight-binding theory unless a representation exists in which the coordinate operator is diagonal. The present approach takes this basis as fundamental and uses group theory to construct symmetrized linear combinations of discrete coordinate eigenkets. This produces orthogonal atomic-like orbitals that may be used as a tight-binding basis. The coordinate matrix in the latter basis includes intra-atomic matrix elements between different orbitals on the same atom. Lattice gauge theory is then used to define discrete electromagnetic fields and their interaction with electrons. Local gauge symmetry is shown to impose strong restrictions limiting the range of the Hamiltonian in the coordinate basis. The theory is applied to the semiconductors Ge and Si, for which it is shown that a basis of 15 orbitals per atom provides a satisfactory description of the valence bands and the lowest conduction bands. Calculations of the dielectric function demonstrate that this model yields an accurate joint density of states, but underestimates the oscillator strength by about 20% in comparison to a nonlocal empirical pseudopotential calculation.
New interatomic potentials describing defects, plasticity and high temperature phase transitions for Ti are presented. Fitting the martensitic hcp-bcc phase transformation temperature requires an efficient and accurate method to determine it. We apply a molecular dynamics (MD) method based on determination of the melting temperature of competing solid phases, and Gibbs-Helmholtz integration, and a lattice-switch Monte Carlo method (LSMC): these agree on the hcp-bcc transformation temperatures to within 2 K. We were able to develop embedded atom potentials which give a good fit to either low or high temperature data, but not both. The first developed potential (Ti1) reproduces the hcp-bcc transformation and melting temperatures and is suitable for the simulation of phase transitions and bcc Ti. Two other potentials (Ti2 and Ti3) correctly describe defect properties, and can be used to simulate plasticity or radiation damage in hcp Ti. The fact that a single EAM potential cannot describe both low and high temperature phases may be attributed to neglect of electronic degrees of freedom, notably bcc has a much higher electronic entropy. A temperature-dependent potential obtained from the combination of potentials Ti1 and Ti2 may be used to simulate Ti properties at any temperature.
By deriving a tight-binding model, we demonstrate a mechanism of forming a nodal line of Dirac points in a single-component molecular conductor [Pt(dmtd)$_2$] [Zhou {it et al.}, Chem. Commun. {bfseries 55}, 3327 (2019)], consisting of HOMO and LUMO. The nodal line is obtained as the intersection of two surfaces, where one corresponds to the HOMO-LUMO band crossing and another is vanishing of the HOMO-LUMO couplings due to their different symmetries. The latter property is essential for the Dirac electron in molecular conductors. The nature of the open nodal line is discussed in terms of the parity of the wavefunctions at eight TRIMs (time reversal invariant momenta).
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.