No Arabic abstract
We present a scheme for the calculation of linear optical properties by the all-electron full-potential linearized augmented planewave (LAPW) method. A summary of the theoretical background for the derivation of the dielectric tensor within the random-phase approximation is provided together with symmetry considerations and the relation between the optical constants. The momentum matrix elements are evaluated in detail for the LAPW basis, and the interband as well as the intraband contributions to the dielectric tensor are given. Results are presented for the metals aluminum and gold, where we crosscheck our results by sumrules. We find that the optical spectra can be extremely sensitive to the Brillouin zone sampling. For gold, the influence of relativistic effects on the dielectic function is investigated. It is shown that the scalar-relativistic effect is much more important than spin-orbit coupling. The interpretability of the Kohn-Sham eigenstates in terms of excited states is discussed.
In solid-state physics, energies of crystals are usually computed with a plane-wave discretization of Kohn-Sham equations. However the presence of Coulomb singularities requires the use of large plane-wave cut-offs to produce accurate numerical results. In this paper, an analysis of the plane-wave convergence of the eigenvalues of periodic linear Hamiltonians with Coulomb potentials using the variational projector-augmented wave (VPAW) method is presented. In the VPAW method, an invertible transformation is applied to the original eigenvalue problem, acting locally in balls centered at the singularities. In this setting, a generalized eigenvalue problem needs to be solved using plane-waves. We show that cusps of the eigenfunctions of the VPAW eigenvalue problem at the positions of the nuclei are significantly reduced. These eigenfunctions have however a higher-order derivative discontinuity at the spheres centered at the nuclei. By balancing both sources of error, we show that the VPAW method can drastically improve the plane-wave convergence of the eigenvalues with a minor additional computational cost. Numerical tests are provided confirming the efficiency of the method to treat Coulomb singularities.
We propose a simple linear scaling expression in reciprocal space for evaluating the ion--electron potential of crystalline solids. The expression replaces the long-range ion--electron potential with an equivalent localized charge distribution and corresponding boundary conditions on the unit cell. Given that no quadratic scaling structure factor is required---as used in traditional methods---the expression shows inherent linear behavior, and is well suited to simulating large-scale systems within orbital-free density functional theory. The scheme is implemented in the ATLAS software package and benchmarked by using a solid Mg bcc lattice containing tens of thousands of atoms in the unit cell. The test results show that the method can efficiently model large crystals with high computational accuracy.
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 many 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.
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 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 response. 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.