No Arabic abstract
The self-consistent charge density functional tight-binding (DFTB) theory is a useful tool for realizing the electronic structures of large molecular complex systems. In this study, we analyze the electronic structure of C61, formed by fullerene C60 with a carbon adatom, using the fully localized limit and pseudo self-interaction correction methods of DFTB to adjust the Hubbard U parameter (DFTB+U). The results show that both the methods used to adjust U can significantly reduce the molecular orbital energy of occupied states localized on the defect carbon atom and improve the gap between highest occupied molecular orbital(HOMO) and lowest unoccupied molecular orbital(LUMO) of C61. This work will provide a methodological reference point for future DFTB calculations of the electronic structures of carbon materials.
The time-dependent density functional based tight-binding (TD-DFTB) approach is generalized to account for fractional occupations. In addition, an on-site correction leads to marked qualitative and quantitative improvements over the original method. Especially, the known failure of TD-DFTB for the description of sigma -> pi* and n -> pi* excitations is overcome. Benchmark calculations on a large set of organic molecules also indicate a better description of triplet states. The accuracy of the revised TD-DFTB method is found to be similar to first principles TD-DFT calculations at a highly reduced computational cost. As a side issue, we also discuss the generalization of the TD-DFTB method to spin-polarized systems. In contrast to an earlier study [Trani et al., JCTC 7 3304 (2011)], we obtain a formalism that is fully consistent with the use of local exchange-correlation functionals in the ground state DFTB method.
We obtain parameters for non-orthogonal and orthogonal TB models from two-atomic molecules for all combinations of elements of period 1 to 6 and group 3 to 18 of the periodic table. The TB bond parameters for 1711 homoatomic and heteroatomic dimers show clear chemical trends. In particular, using our parameters we compare to the rectangular d-band model, the reduced sp TB model as well as canonical TB models for sp- and d-valent systems which have long been used to gain qualitative insight into the interatomic bond. The transferability of our dimer-based TB bond parameters to bulk systems is discussed exemplarily for the bulk ground-state structures of Mo and Si. Our dimer-based TB bond parameters provide a well-defined and promising starting point for developing refined TB parameterizations and for making the insight of TB available for guiding materials design across the periodic table.
The self-consistent evaluation of Hubbard parameters using linear-response theory is crucial for quantitatively predictive calculations based on Hubbard-corrected density-functional theory. Here, we extend a recently-introduced approach based on density-functional perturbation theory (DFPT) for the calculation of the on-site Hubbard $U$ to also compute the inter-site Hubbard $V$. DFPT allows to reduce significantly computational costs, improve numerical accuracy, and fully automate the calculation of the Hubbard parameters by recasting the linear response of a localized perturbation into an array of monochromatic perturbations that can be calculated in the primitive cell. In addition, here we generalize the entire formalism from norm-conserving to ultrasoft and projector-augmented wave formulations, and to metallic ground states. After benchmarking DFPT against the conventional real-space Hubbard linear response in a supercell, we demonstrate the effectiveness of the present extended Hubbard formulation in determining the equilibrium crystal structure of Li$_x$MnPO$_4$ (x=0,1) and the subtle energetics of Li intercalation.
We present ab initio results for the electron-phonon interaction of the Gamma-point phonons in the tetragonal high-temperature phase of La2CuO4. Eigenfrequencies and eigenvectors for the symmetry-allowed phonon modes are calculated with the full-potential augmented plane wave+local orbitals method using the frozen phonon approach. It is found that the Gamma-point phonons with the strongest electron-phonon interaction are the A{2u} modes with 236 cm^{-1}, 131 cm^{-1} and 476 cm^{-1}. To take effect of strong electron on-site interaction into account we use generalized tight-binding method that results in the interaction of phonons with Hubbard fermions forming quasiparticles band structure. Finally, the matrix elements of Hubbard fermion-phonon interaction and their reduction due to strong electron correlation are obtained.
Density-functional theory (DFT) has revolutionized computational prediction of atomic-scale properties from first principles in physics, chemistry and materials science. Continuing development of new methods is necessary for accurate predictions of new classes of materials and properties, and for connecting to nano- and mesoscale properties using coarse-grained theories. JDFTx is a fully-featured open-source electronic DFT software designed specifically to facilitate rapid development of new theories, models and algorithms. Using an algebraic formulation as an abstraction layer, compact C++11 code automatically performs well on diverse hardware including GPUs. This code hosts the development of joint density-functional theory (JDFT) that combines electronic DFT with classical DFT and continuum models of liquids for first-principles calculations of solvated and electrochemical systems. In addition, the modular nature of the code makes it easy to extend and interface with, facilitating the development of multi-scale toolkits that connect to ab initio calculations, e.g. photo-excited carrier dynamics combining electron and phonon calculations with electromagnetic simulations.