Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userSelf-consistent Hubbard parameters from density-functional perturbation theory in the ultrasoft and projector-augmented wave formulations
129
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
We present in full detail a newly developed formalism enabling density functional perturbation theory (DFPT) calculations from a DFT+$U$ ground state. The implementation includes ultrasoft pseudopotentials and is valid for both insulating and metallic systems. It aims at fully exploiting the versatility of DFPT combined with the low-cost DFT+$U$ functional. This allows to avoid computationally intensive frozen-phonon calculations when DFT+$U$ is used to eliminate the residual electronic self-interaction from approximate functionals and to capture the localization of valence electrons e.g. on $d$ or $f$ states. In this way, the effects of electronic localization (possibly due to correlations) are consistently taken into account in the calculation of specific phonon modes, Born effective charges, dielectric tensors and in quantities requiring well converged sums over many phonon frequencies, as phonon density of states and free energies. The new computational tool is applied to two representative systems, namely CoO, a prototypical transition metal monoxide and LiCoO$_2$, a material employed for the cathode of Li-ion batteries. The results show the effectiveness of our formalism to capture in a quantitatively reliable way the vibrational properties of systems with localized valence electrons.
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.
We present a computational scheme for orbital-free density functional theory (OFDFT) that simultaneously provides access to all-electron values and preserves the OFDFT linear scaling as a function of the system size. Using the projector augmented-wave method (PAW) in combination with real-space methods we overcome some obstacles faced by other available implementation schemes. Specifically, the advantages of using the PAW method are two fold. First, PAW reproduces all-electron values offering freedom in adjusting the convergence parameters and the atomic setups allow tuning the numerical accuracy per element. Second, PAW can provide a solution to some of the convergence problems exhibited in other OFDFT implementations based on Kohn-Sham codes. Using PAW and real-space methods, our orbital-free results agree with the reference all-electron values with a mean absolute error of 10~meV and the number of iterations required by the self-consistent cycle is comparable to the KS method. The comparison of all-electron and pseudopotential bulk modulus and lattice constant reveal an enormous difference, demonstrating that in order to assess the performance of OFDFT functionals it is necessary to use implementations that obtain all-electron values. The proposed combination of methods is the most promising route currently available. We finally show that a parametrized kinetic energy functional can give lattice constants and bulk moduli comparable in accuracy to those obtained by the KS PBE method, exemplified with the case of diamond.
We present a method of calculating crystal field coefficients of rare-earth/transition-metal (RE-TM) magnets within density-functional theory (DFT). The principal idea of the method is to calculate the crystal field potential of the yttrium analogue (Y-analogue) of the RE-TM magnet, i.e. the material where the lanthanide elements have been substituted with yttrium. The advantage of dealing with Y-analogues is that the methodological and conceptual difficulties associated with treating the highly-localized 4f electrons in DFT are avoided, whilst the nominal valence electronic structure principally responsible for the crystal field is preserved. In order to correctly describe the crystal field potential in the core region of the atoms we use the projector-augmented wave formalism of DFT, which allows the reconstruction of the full charge density and electrostatic potential. The Y-analogue crystal field potentials are combined with radial 4f charge densities obtained in self-interaction-corrected calculations on the lanthanides to obtain crystal field coefficients. We demonstrate our method on a test set of 10 materials comprising 9 RE-TM magnets and elemental Tb. We show that the calculated easy directions of magnetization agree with experimental observations, including a correct description of the anisotropy within the basal plane of Tb and NdCo$_5$. We further show that the Y-analogue calculations generally agree quantitatively with previous calculations using the open-core approximation to treat the 4f electrons, and argue that our simple approach may be useful for large-scale computational screening of new magnetic materials.
The knowledge of effective masses is a key ingredient to analyze numerous properties of semiconductors, like carrier mobilities, (magneto-)transport properties, or band extrema characteristics yielding carrier densities and density of states. Currently, these masses are usually calculated using finite-difference estimation of density functional theory (DFT) electronic band curvatures. However, finite differences require an additional convergence study and are prone to numerical noise. Moreover, the concept of effective mass breaks down at degenerate band extrema. We assess the former limitation by developing a method that allows to obtain the Hessian of DFT bands directly, using density functional perturbation theory (DFPT). Then, we solve the latter issue by adapting the concept of `transport equivalent effective mass to the $vec{k} cdot hat{vec{p}}$ framework. The numerical noise inherent to finite-difference methods is thus eliminated, along with the associated convergence study. The resulting method is therefore more general, more robust and simpler to use, which makes it especially appropriate for high-throughput computing. After validating the developed techniques, we apply them to the study of silicon, graphane, and arsenic. The formalism is implemented into the ABINIT software and supports the norm-conserving pseudopotential approach, the projector augmented-wave method, and the inclusion of spin-orbit coupling. The derived expressions also apply to the ultrasoft pseudopotential method.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
