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.
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 (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 re
vised Perdew-Burke-Ernzerhof GGA that improves equilibrium properties of densely-packed solids and their surfaces.
We apply the coupled dynamics of time-dependent density functional theory and Maxwell equations to the interaction of intense laser pulses with crystalline silicon. As a function of electromagnetic field intensity, we see several regions in the respo
nse. At the lowest intensities, the pulse is reflected and transmitted in accord with the dielectric response, and the characteristics of the energy deposition is consistent with two-photon absorption. The absorption process begins to deviate from that at laser intensities ~ 10^13 W/cm^2, where the energy deposited is of the order of 1 eV per atom. Changes in the reflectivity are seen as a function of intensity. When it passes a threshold of about 3 times 1012 W/cm2, there is a small decrease. At higher intensities, above 2 times 10^13 W/cm^2, the reflectivity increases strongly. This behavior can be understood qualitatively in a model treating the excited electron-hole pairs as a plasma.
We consider several aspects of high-order harmonic generation in solids: the effects of elastic and inelastic scattering; varying pulse characteristics; and inclusion of material-specific parameters through a realistic band structure. We reproduce ma
ny observed characteristics of high harmonic generation experiments in solids including the formation of only odd harmonics in inversion-symmetric materials, and the nonlinear formation of high harmonics with increasing field. We find that the harmonic spectra are fairly robust against elastic and inelastic scattering. Furthermore, we find that the pulse characteristics play an important role in determining the harmonic spectra.
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 gra
dient 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.