No Arabic abstract
We propose a generalized gradient approximation (GGA) for the angle- and system-averaged exchange-correlation hole of a many-electron system. This hole, which satisfies known exact constraints, recovers the PBEsol (Perdew-Burke-Ernzerhof for solids) exchange-correlation energy functional, a GGA that accurately describes the equilibrium properties of densely packed solids and their surfaces. We find that our PBEsol exchange-correlation hole describes the wavevector analysis of the jellium exchange-correlation surface energy in agreement with a sophisticated time-dependent density-functional calculation (whose three-dimensional wavevector analysis we report here).
Successful modern generalized gradient approximations (GGA) are biased toward atomic energies. Restoration of the first-principles gradient expansion for the exchange energy over a wide range of density gradients eliminates this bias. We introduce PBEsol, a revised Perdew-Burke-Ernzerhof GGA that improves equilibrium properties for many densely-packed solids and their surfaces.
Successful modern generalized gradient approximations (GGAs) are biased toward atomic energies. Restoration of the first-principles gradient expansion for exchange over a wide range of density gradients eliminates this bias. We introduce PBEsol, a revised Perdew-Burke-Ernzerhof GGA that improves equilibrium properties of densely-packed solids and their surfaces.
We present a detailed study of the coupling-constant-averaged exchange-correlation hole density at a jellium surface, which we obtain in the random-phase approximation (RPA) of many-body theory. We report contour plots of the exchange-only and exchange-correlation hole densities, the integration of the exchange-correlation hole density over the surface plane, the on-top correlation hole, and the energy density. We find that the on-top correlation hole is accurately described by local and semilocal density-functional approximations. We also find that for electrons that are localized far outside the surface the main part of the corresponding exchange-correlation hole is localized at the image plane.
We develop numerical methods to calculate electron dynamics in crystalline solids in real-time time-dependent density functional theory employing exchange-correlation potentials which reproduce band gap energies of dielectrics; a meta generalized gradient approximation (meta-GGA) proposed by Tran and Blaha [Phys. Rev. Lett. 102, 226401 (2009)] (TBm-BJ) and a hybrid functional proposed by Heyd, Scuseria, and Ernzerhof [J. Chem. Phys. 118, 8207 (2003)] (HSE). In time evolution calculations employing the TB-mBJ potential, we have found it necessary to adopt a predictor-corrector step for stable time-evolution. Since energy functional is not known for the TB-mBJ potential, we propose a method to evaluate electronic excitation energy without referring to the energy functional. Calculations using the HSE hybrid functional is computationally expensive due to the nonlocal Fock-like term. We develop a computational method for the operation of the Fock-like term in Fourier space, for which we employ massively parallel computers equipped with graphic processing units. To demonstrate significances of utilizing potentials providing correct band gap energies, we compare electronic excitations induced by femtosecond laser pulses using the TB-mBJ, HSE, and a simple local density approximation (LDA). At low laser intensities, electronic excitations are found to be sensitive to the band gap energy: results using TB-mBJ and HSE are close to each other, while the excitation of the LDA calculation is more intensive than the others. At high laser intensities close to a damage threshold, we have found that electronic excitation energies are similar among the three cases.
In our previous study [Phys. Rev. B 86, 201104 (2012)] we introduced the so called quasi-non-uniform gradient-level exchange-correlation approximation (QNA) and demonstrated its strength in producing highly accurate equilibrium volumes for metals and their alloys within the density-functional theory. In this paper we extend the scheme to include the accuracy of bulk modulus as an additional figure of merit and show that this scheme is flexible enough as to allow the computation of accurate equilibrium volumes and bulk moduli at the same time. The power and feasibility of this scheme is demonstrated on NiAl and FeV binary alloys.