Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userDensity functional theory method for twisted geometries with application to torsional deformations in group-IV nanotubes
139
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present a real-space formulation and implementation of Kohn-Sham Density Functional Theory suited to twisted geometries, and apply it to the study of torsional deformations of X (X = C, Si, Ge, Sn) nanotubes. Our formulation is based on higher order finite difference discretization in helical coordinates, uses ab intio pseudopotentials, and naturally incorporates rotational (cyclic) and screw operation (i.e., helical) symmetries. We discuss several aspects of the computational method, including the form of the governing equations, details of the numerical implementation, as well as its convergence, accuracy and efficiency properties. The technique presented here is particularly well suited to the first principles simulation of quasi-one-dimensional structures and their deformations, and many systems of interest can be investigated using small simulation cells containing just a few atoms. We apply the method to systematically study the properties of single-wall zigzag and armchair group-IV nanotubes, as they undergo twisting. For the range of deformations considered, the mechanical behavior of the tubes is found to be largely consistent with isotropic linear elasticity, with the torsional stiffness varying as the cube of the nanotube radius. Furthermore, for a given tube radius, this quantity is seen to be highest for carbon nanotubes and the lowest for those of tin, while nanotubes of silicon and germanium have intermediate values close to each other. We also describe different aspects of the variation in electronic properties of the nanotubes as they are twisted. In particular, we find that akin to the well known behavior of armchair carbon nanotubes, armchair nanotubes of silicon, germanium and tin also exhibit bandgaps that vary periodically with imposed rate of twist, and that the periodicity of the variation scales in an inverse quadratic manner with the tube radius.
rate research
Read More
We propose a method to decompose the total energy of a supercell containing defects into contributions of individual atoms, using the energy density formalism within density functional theory. The spatial energy density is unique up to a gauge transformation, and we show that unique atomic energies can be calculated by integrating over Bader and charge-neutral volumes for each atom. Numerically, we implement the energy density method in the framework of the Vienna ab initio simulation package (VASP) for both norm-conserving and ultrasoft pseudopotentials and the projector augmented wave method, and use a weighted integration algorithm to integrate the volumes. The surface energies and point defect energies can be calculated by integrating the energy density over the surface region and the defect region, respectively. We compute energies for several surfaces and defects: the (110) surface energy of GaAs, the mono-vacancy formation energies of Si, the (100) surface energy of Au, and the interstitial formation energy of O in the hexagonal close-packed Ti crystal. The surface and defect energies calculated using our method agree with size-converged calculations of the difference between the total energies of the system with and without the defect. Moreover, the convergence of the defect energies with size can be found from a single calculation.
We propose a new method for accurately calculating electrical transport properties of a lightly-doped thermoelectric material from density functional theory (DFT) calculations, based on experimental data and density functional theory results for the corresponding undoped material. We employ this approach because hybrid DFT calculations are prohibitive for the large supercells required to model low dopant concentrations comparable to those achieved experimentally for high-performing thermoelectrics. Using zinc antimonide as our base material, we find that the electrical transport properties calculated with DFT and Boltzmann transport theory exhibit the same trends with changes in chemical potential as those computed with hybrid DFT, and propose a fitting algorithm that involves adjusting the computed Fermi energy so that the resulting Seebeck coefficient trends with temperature match experimental trends. We confirm the validity of this approach in reproducing the experimental trends in electrical conductivity and Seebeck coefficient versus temperature for Bi-doped $beta-$Zn$_4$Sb$_3$. We then screen various transition metal cation dopants, including copper and nickel, and find that a Cu dopant concentration of 2.56% in Zn$_{39}$Sb$_{30}$ exhibited a 14% increase in the thermoelectric power factor for temperatures between 300-400 K. We thus propose that transition metal dopants may significantly improve the thermoelectric performance of the host material, compared to heavy and/or rare-earth dopants.
We present a symmetry-adapted real-space formulation of Kohn-Sham density functional theory for cylindrical geometries and apply it to the study of large X (X=C, Si, Ge, Sn) nanotubes. Specifically, starting from the Kohn-Sham equations posed on all of space, we reduce the problem to the fundamental domain by incorporating cyclic and periodic symmetries present in the angular and axial directions of the cylinder, respectively. We develop a high-order finite-difference parallel implementation of this formulation, and verify its accuracy against established planewave and real-space codes. Using this implementation, we study the band structure and bending properties of X nanotubes and Xene sheets, respectively. Specifically, we first show that zigzag and armchair X nanotubes with radii in the range 1 to 5 nm are semiconducting. In particular, we find an inverse linear dependence of the bandgap with respect to the radius for all nanotubes, other than the armchair and zigzag type III carbon variants, for which we find an inverse quadratic dependence. Next, we exploit the connection between cyclic symmetry and uniform bending deformations to calculate the bending moduli of Xene sheets in both zigzag and armchair directions. We find Kirchhoff-Love type bending behavior for all sheets, with graphene and stanene possessing the largest and smallest moduli, respectively. In addition, other than graphene, the sheets demonstrate significant anisotropy, with larger bending moduli along the armchair direction. Finally, we demonstrate that the proposed approach has very good parallel scaling and is highly efficient, enabling ab initio simulations of unprecedented size for systems with a high degree of cyclic symmetry. In particular, we show that even micron-sized nanotubes can be simulated with modest computational effort.
The electronic structure of surfaces plays a key role in the properties of quantum devices. However, surfaces are also the most challenging to simulate and engineer. Here, we study the electronic structure of InAs(001), InAs(111), and InSb(110) surfaces using a combination of density functional theory (DFT) and angle-resolved photoemission spectroscopy (ARPES). We were able to perform large-scale first principles simulations and capture effects of different surface reconstructions by using DFT calculations with a machine-learned Hubbard U correction [npj Comput. Mater. 6, 180 (2020)]. To facilitate direct comparison with ARPES results, we implemented a bulk unfolding scheme by projecting the calculated band structure of a supercell surface slab model onto the bulk primitive cell. For all three surfaces, we find a good agreement between DFT calculations and ARPES. For InAs(001), the simulations clarify the effect of the surface reconstruction. Different reconstructions are found to produce distinctive surface states. For InAs(111) and InSb(110), the simulations help elucidate the effect of oxidation. Owing to larger charge transfer from As to O than from Sb to O, oxidation of InAs(111) leads to significant band bending and produces an electron pocket, whereas oxidation of InSb(110) does not. Our combined theoretical and experimental results may inform the design of quantum devices based on InAs and InSb semiconductors, e.g., topological qubits utilizing the Majorana zero modes.
The hyperfine interaction is an important probe for understanding the structure and symmetry of defects in a semiconductor. Density-functional theory has shown that it can provide useful first-principles predictions for both the hyperfine tensor and the hyperfine constants that arise from it. Recently there has been great interest in using group-IV impurity-vacancy color centers X$V^-$ (where X = Si, Ge, Sn, or Pb and $V$ is a carbon vacancy) for important applications in quantum computing and quantum information science. In this paper, we have calculated the hyperfine tensors for these X$V^-$ color centers using the HSE06 screened Hartree-Fock hybrid exchange-correlation functional with the inclusion of core electron spin polarization. We have compared our results to calculations which only use the PBE exchange-correlation functional without the inclusion of core electron spin polarization and we have found our results are in very good agreement with available experimental results. Finally, we have theoretically shown that these X$V^-$ color centers exhibit a Jahn-Teller distortion which explains the observed anisotropic distribution of the hyperfine constants among the neighboring $^{13}$C nuclear spins.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
