No Arabic abstract
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.
A systemically theoretical study has been presented to explored the crystal structures and electronic characteristics of polycyclic aromatic hydrocarbons (PAHs), such as solid phenanthrene, picene, 1,2;8,9-dibenzopentacene, and 7-phenacenes, since these PAHs exhibited the superconductivity when potassium doping into. For tripotassium-doped phenanthrene and picene, we demonstrate the K atomic positions to fit the experimental lattice parameters, and analyze the distinction between the stablest configuration and the fitted experimental one. Based on the first-principles calculations, for the first time, we predict the possible crystal configurations of pristine and tripotassium-doped 1,2;8,9-dibenzopentacene and 7-phenacenes, respectively. For these four PAHs, the electronic structures after doping are investigated in details. The results show that the electronic characters near the Fermi level are high sensitive to structure. Because of the change of the benzene rings arrangement, the 1,2;8,9-dibenzopentacene exhibits visibly different band structures from other three PAHs. In these metallic PAHs, two bands cross the Fermi level which results in the complicated multiband feature of Fermi surfaces. Fascinatingly, we find that the electronic states of potassium contribute to the Fermi surfaces especially for K-3$d$ electrons, which improves a way to understand this superconductivity. As a result, we suggest that the rigid-band picture is invalidated due to the hybridization between K atoms and PAH molecules as well as the rearrangement and distortion of PAH molecules.
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.
A new approach is developed to calculate temperature dependent Seebeck coefficient of heavily doped systems by using Boltzmann transport theory and electron density of states (DOS) obtained from density functional calculations. This approach is applied to heavily doped La:STO with DOS from tetrahedral method and Fermi energy using Fermi integrals. The calculated Seebeck coefficient agrees with the experimental data nearly quantitatively, which proved the accuracy of this approach. The influence of the Fermi energy and asymmetry of DOS on the Seebeck coefficient is analyzed.
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.
Using the first-principles density-functional theory plan-wave pseudopotential method, we investigate the structure and magnetism in 25% Mn substitutive and interstitial doped monoclinic, tetragonal and cubic ZrO2 systematically. Our studies show that the introduction of Mn impurities into ZrO2 not only stabilizes the high temperature phase, but also endows ZrO2 with magnetism. Based on the simple crystal field theory (CFT), we discuss the origination of magnetism in Mn doped ZrO2. Moreover, we discuss the effect of electron donor on magnetic semiconductors, and the possibility as electronic structure modulator.