Publications

Recently, it has been shown that under pressure, unexpected and counterintuitive chemical compounds become stable. Laser shock experiments (A. Rode, unpublished) on alumina (Al 2 O 3 ) have shown non-equilibrium decomposition of alumina with the formation of free Al and a mysterious transparent phase. Inspired by these observations, we have explored the possibility of the formation of new chemical compounds in the system Al-O. Using the variable-composition structure prediction algorithm USPEX, in addition to the well-known Al 2 O 3 , we have found two extraordinary compounds Al 4 O 7 and AlO 2 to be thermodynamically stable in the pressure ranges 330-443 GPa and above 332 GPa, respectively. Both of these compounds at the same time contain oxide O 2− and peroxide O 2 2− ions and both are insulating. Peroxo-groups are responsible for gap states, which significantly reduce the electronic band gap of both Al 4 O 7 and AlO 2 .

When wurtzite ZnO is sliced perpendicular to the (0001) axis, two different polar surfaces, the (0001)-Zn and (000$\overline{1}$)-O terminated surfaces, are formed. In a simple ionic picture, both surfaces are electrostatically unstable due to a diverging electrostatic energy. Although the ionic picture is an oversimplification, the surfaces adopt a modified surface structure to compensate for the polarity. In close collaboration with experiment, a hexagonal honeycomblike reconstruction has been suggested [J. V. Lauritsen et al., ACS Nano 5, 5987 (2011)]. The remarkable observation is that the (000$\overline{1}$)-O surface behaves very differently than the (0001)-Zn surface. Here, we present a detailed density functional theory investigation of the (000$\overline{1}$)-O surface, including a systematic investigation of H and Zn coverage as well as an investigation of various surface reconstructions. The difference between the two polar surfaces is explained by the different bonding preferences of Zn and O atoms: as a $d$ element, Zn atoms are more flexible in their bond formation than O atoms.

We demonstrate that natural orbitals allow for reducing the computational cost of wave function based correlated calculations, especially for atoms and molecules in a large box, when a plane wave basis set under periodic boundary conditions is used. The employed natural orbitals are evaluated on the level of second-order Møller–Plesset perturbation theory (MP2), which requires a computational effort that scales as O(N5), where N is a measure of the system size. Moreover, we find that a simple approximation reducing the scaling to O(N4) yields orbitals that allow for a similar reduction of the number of virtual orbitals. The MP2 natural orbitals are applied to coupled-cluster singles and doubles (CCSD) as well as full configuration interaction Quantum Monte Carlo calculations of the H2 molecule to test our implementation. Finally, the atomization energies of the LiH molecule and solid are calculated on the level of MP2 and CCSD.

Optical conductivity of an interacting polaron gas is calculated within an extended random phase approximation which takes into account mixing of collective excitations of the electron gas with LO phonons. This mixing is important for the optical response of strongly polar crystals where the static dielectric constant is rather high: strontium titanate is the case. The present calculation sheds light on unexplained features of experimentally observed optical conductivity spectra in $n$-doped SrTiO$_{3}$. These features appear to be due to dynamic screening of the electron-electron interaction by polar optical phonons and hence do not require additional mechanisms for the explanation.

Electron-phonon physics has become an integral part of today's materials science world. In this work, the authors present $a\phantom{\rule{0}{0ex}}b$ $i\phantom{\rule{0}{0ex}}n\phantom{\rule{0}{0ex}}i\phantom{\rule{0}{0ex}}t\phantom{\rule{0}{0ex}}i\phantom{\rule{0}{0ex}}o$ results for the zero-point renormalization of the band gap of various semiconductors and insulators. They use a new implementation of the VASP code that utilizes the sophisticated projector augmented wave method to calculate the electron-phonon perturbation. Two complementary ways of treating the all-electron information are discussed in this context.

No abstract is provided for this article.

A Comment on the Letter by B. Arnaud et al., [Phys. Rev. Lett. 96, 026402 (2006)]. The authors of the Letter offer a Reply.Received 8 October 2007DOI:https://doi.org/10.1103/PhysRevLett.100.189701©2008 American Physical Society

Several proton-disordered crystalline ice structures are known to proton order at sufficiently low temperatures, provided that the right preparation procedure is used. For cubic ice, ice Ic, however, no proton ordering has been observed so far. Here, we subject ice Ic to an experimental protocol similar to that used to proton order hexagonal ice. In situ FT-IR spectroscopy carried out during this procedure reveals that the librational band of the spectrum narrows and acquires a structure that is observed neither in proton-disordered ice Ic nor in ice XI, the proton-ordered variant of hexagonal ice. On the basis of vibrational spectra computed for ice Ic and four of its proton-ordered variants using classical molecular dynamics and ab initio simulations, we conclude that the features of our experimental spectra are due to partial proton ordering, providing the first evidence of proton ordering in cubic ice. We further find that the proton-ordered structure with the lowest energy is ferroelectric, while the structure with the second lowest energy is weakly ferroelectric. Both structures fit the experimental spectral similarly well such that no unique assignment of proton order is possible based on our results.

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.

We present quasiparticle (QP) energies from fully self-consistent $GW$ (sc$GW$) calculations for a set of prototypical semiconductors and insulators within the framework of the projector-augmented wave methodology. To obtain converged results, both finite basis-set corrections and $k$-point corrections are included, and a simple procedure is suggested to deal with the singularity of the Coulomb kernel in the long-wavelength limit, the so called head correction. It is shown that the inclusion of the head corrections in the sc$GW$ calculations is critical to obtain accurate QP energies with a reasonable $k$-point set. We first validate our implementation by presenting detailed results for the selected case of diamond, and then we discuss the converged QP energies, in particular the band gaps, for the entire set of gapped compounds and compare them to single-shot $G_0W_0$, QP self-consistent $GW$, and previously available sc$GW$ results as well as experimental results.