In this work we derive closed expressions for the head of the frequency-dependent microscopic polarizability matrix in the projector-augmented wave ͑PAW͒ methodology. Contrary to previous applications, the longitudinal expression is utilized, resulting in dielectric properties that are largely independent of the applied potentials. The improved accuracy of the present approach is demonstrated by comparing the longitudinal and transversal expressions of the polarizability matrix for a number of cubic semiconductors and one insulator, i.e., Si, SiC, AlP, GaAs, and diamond ͑C͒, respectively. The methodology is readily extendable to more complicated nonlocal Hamiltonians or to the calculation of the macroscopic dielectric matrix including local field effects in the random phase or density functional approximation, which is demonstrated for the previously mentioned model systems. Furthermore, density functional perturbation theory is extended to the PAW method, and the respective results are compared to those obtained by summation over the conduction band states.

Hybrid Fock exchange/density functional theory functionals have shown to be very successful in describing a wide range of molecular properties. For periodic systems, however, the long-range nature of the Fock exchange interaction and the resultant large computational requirements present a major drawback. This is especially true for metallic systems, which require a dense Brillouin zone sampling. Recently, a new hybrid functional [HSE03, J. Heyd, G. E. Scuseria, and M. Ernzerhof, J. Chem. Phys. 118, 8207 (2003)] that addresses this problem within the context of methods that evaluate the Fock exchange in real space was introduced. We discuss the advantages the HSE03 functional brings to methods that rely on a reciprocal space description of the Fock exchange interaction, e.g., all methods that use plane wave basis sets. Furthermore, we present a detailed comparison of the performance of the HSE03 and PBE0 functionals for a set of archetypical solid state systems by calculating lattice parameters, bulk moduli, heats of formation, and band gaps. The results indicate that the hybrid functionals indeed often improve the description of these properties, but in several cases the results are not yet on par with standard gradient corrected functionals. This concerns in particular metallic systems for which the bandwidth and exchange splitting are seriously overestimated.

No abstract

We present converged ab initio calculations of the optical absorption spectra of single-layer, double-layer, and bulk MoS 2 . Both the quasiparticle-energy calculations (on the level of the GW approximation ) and the calculation of the absorption spectra (on the level of the Bethe-Salpeter equation) explicitly include spin-orbit coupling, using the full spinorial Kohn-Sham wave functions as input. Without excitonic effects, the absorption spectra would have the form of a step function, corresponding to the joint density of states of a parabolic band dispersion in two dimensions. This profile is deformed by a pronounced bound excitonic peak below the continuum onset. The peak is split by spin-orbit interaction in the case of single-layer and (mostly) by interlayer interaction in the case of double-layer and bulk MoS 2 . The resulting absorption spectra are thus very similar in the three cases, but the interpretation of the spectra is different. Differences in the spectra can be seen in the shape of the absorption spectra at 3 eV where the spectra of the single and double layers are dominated by a strongly bound exciton.

No abstract