The cohesive energy of superheavy element copernicium determined from accurate relativistic coupled-cluster theory.

The cohesive energy of bulk copernicium is accurately determined using the incremental method within a relativistic coupled-cluster approach. For the lowest energy structure of hexagonal close-packed (hcp) symmetry, we obtain a cohesive energy of -36.3 kJ mol-1 (inclusion of uncertainties leads to a lower bound of -39.6 kJ mol-1), in excellent agreement with the experimentally estimated sublimation enthalpy of -38 kJ mol-1 [R. Eichler et al., Angew. Chem. Int. Ed., 2008, 47, 3262]. At the coupled-cluster singles, doubles and perturbative triples level of theory, we find that the hcp structure is energetically quasi-degenerate with both face-centred and body-centred cubic structures. These results provide a basis for testing various density-functionals, of which the PBEsol functional yields a cohesive energy of -34.1 kJ mol-1 in good agreement with our coupled-cluster value.

Embedding for bulk systems using localized atomic orbitals.

We extend projection-based embedding techniques to bulk systems to treat point defects in semiconductors and insulators. To avoid non-additive kinetic energy contributions, we construct the density partition using orthogonal subsets of orbitals. We have implemented our approach in the popular Vienna ab initio simulation package software package. We demonstrate its power for defect structures in silicon and polaron formation in titania, two challenging cases for conventional Kohn-Sham density functional theory.

NMR shieldings from density functional perturbation theory: GIPAW versus all-electron calculations.

We present a benchmark of the density functional linear response calculation of NMR shieldings within the gauge-including projector-augmented-wave method against all-electron augmented-plane-wave+local-orbital and uncontracted Gaussian basis set results for NMR shieldings in molecular and solid state systems. In general, excellent agreement between the aforementioned methods is obtained. Scalar relativistic effects are shown to be quite large for nuclei in molecules in the deshielded limit. The small component makes up a substantial part of the relativistic corrections.

The random phase approximation applied to ice.

Standard density functionals without van der Waals interactions yield an unsatisfactory description of ice phases, specifically, high density phases occurring under pressure are too unstable compared to the common low density phase Ih observed at ambient conditions. Although the description is improved by using functionals that include van der Waals interactions, the errors in relative volumes remain sizable. Here we assess the random phase approximation (RPA) for the correlation energy and compare our results to experimental data as well as diffusion Monte Carlo data for ice. The RPA yields a very balanced description for all considered phases, approaching the accuracy of diffusion Monte Carlo in relative energies and volumes. This opens a route towards a concise description of molecular water phases on surfaces and in cavities.

Accurate surface and adsorption energies from many-body perturbation theory.

Kohn-Sham density functional theory is the workhorse computational method in materials and surface science. Unfortunately, most semilocal density functionals predict surfaces to be more stable than they are experimentally. Naively, we would expect that consequently adsorption energies on surfaces are too small as well, but the contrary is often found: chemisorption energies are usually overestimated. Modifying the functional improves either the adsorption energy or the surface energy but always worsens the other aspect. This suggests that semilocal density functionals possess a fundamental flaw that is difficult to cure, and alternative methods are urgently needed. Here we show that a computationally fairly efficient many-electron approach, the random phase approximation to the correlation energy, resolves this dilemma and yields at the same time excellent lattice constants, surface energies and adsorption energies for carbon monoxide and benzene on transition-metal surfaces.

Cohesive properties and asymptotics of the dispersion interaction in graphite by the random phase approximation.

The structural properties of graphite, such as the interlayer equilibrium distance, the elastic constant, and the net layer binding energy, are obtained using the adiabatic-connection fluctuation-dissipation theorem in the random phase approximation. Excellent agreement is found with the available experimental data; however, our computed binding energy of 48 meV per atom is somewhat smaller than the one obtained by quantum Monte Carlo methods. The asymptotic behavior of the interlayer dispersion interaction, previously derived from analytic approximations, is explicitly demonstrated to follow a d-3 behavior at very large distances.

Second-order Møller-Plesset perturbation theory applied to extended systems. I. Within the projector-augmented-wave formalism using a plane wave basis set.

We present an implementation of the canonical formulation of second-order Møller-Plesset (MP2) perturbation theory within the projector-augmented-wave method under periodic boundary conditions using a plane wave basis set. To demonstrate the accuracy of our approach we show that our result for the atomization energy of a LiH molecule at the Hartree-Fock+MP2 level is in excellent agreement with well converged Gaussian-type-orbital calculations. To establish the feasibility of employing MP2 perturbation theory in its canonical form to systems that are periodic in three dimensions we calculated the cohesive energy of bulk LiH.

Hybrid functionals applied to extended systems.

We present an overview of the description of structural, thermochemical, and electronic properties of extended systems using several well known hybrid Hartree-Fock/density-functional-theory functionals (PBE0, HSE03, and B3LYP). In addition we address a few aspects of the evaluation of the Hartree-Fock exchange interactions in reciprocal space, relevant to all methods that employ a plane wave basis set and periodic boundary conditions.

Accurate band gaps and dielectric properties from one-electron theories (abstract only).

For semiconductor modeling, a major shortcoming of density functional theory is that the predicted band gaps are usually significantly too small. It is generally argued that this shortcoming is related to the fact that density functional theory is a ground state theory, and as a result, one is not allowed to associate the one-electron energies with the energies of quasi-particles. Although this fundamental objection is certainly correct, the modeling of the positioning of donor and acceptor levels in semiconductors faces serious limitations with present density functionals. Several solutions to this problem have been suggested. A particular attractive and fairly simple one is the inclusion of a small fraction of the non-local exchange in the Hamiltonian (hybrid functionals). This approach leads to sensible band gaps for most semiconductors, but fails for ionic solids. A more reliable approach is via many-electron Green's function techniques, which have made tremendous advances in recent years. Here GW calculations in various flavors are presented for small gap and large gap systems, comprising typical semiconductors (Si, SiC, GaAs, GaN, ZnO, ZnS, CdS and AlP), small gap semiconductors (PbS, PbSe, PbTe), insulators (C, BN, MgO, LiF) and noble gas solids (Ar, Ne). The general finding is that single-shot G(0)W(0) calculations based on wavefunctions obtained from conventional density functional theory yield too small band gaps, whereas G(0)W(0) calculations following hybrid functional calculations tend to overestimate the band gaps by roughly the same amount. This is at first sight astonishing, since the hybrid functionals yield very good band gaps themselves. The contradiction is resolved by showing that the inclusion of the attractive electron-hole interactions (excitonic effects) is required to obtain good static and dynamic dielectric functions using hybrid functionals. The corrections are usually incorporated in GW calculations using 'vertex corrections', and in fact inclusion of these vertex corrections rectifies the predicted band gaps. Finally, in order to remove the dependence on the initial wavefunctions, self-consistent GW calculations are presented, again including an approximate treatment of vertex corrections. The results are in excellent agreement with experiment, with a few per cent deviation for all materials considered. We conclude that predictive band gap engineering is now possible with the theoretical description approaching experimental accuracy.

Tuning the balance between dispersion and entropy to design temperature-responsive flexible metal-organic frameworks.

Temperature-responsive flexibility in metal-organic frameworks (MOFs) appeals to the imagination. The ability to transform upon thermal stimuli while retaining a given crystalline topology is desired for specialized sensors and actuators. However, rational design of such shape-memory nanopores is hampered by a lack of knowledge on the nanoscopic interactions governing the observed behavior. Using the prototypical MIL-53(Al) as a starting point, we show that the phase transformation between a narrow-pore and large-pore phase is determined by a delicate balance between dispersion stabilization at low temperatures and entropic effects at higher ones. We present an accurate theoretical framework that allows designing breathing thermo-responsive MOFs, based on many-electron data for the dispersion interactions and density-functional theory entropy contributions. Within an isoreticular series of materials, MIL-53(Al), MIL-53(Al)-FA, DUT-4, DUT-5 and MIL-53(Ga), only MIL-53(Al) and MIL-53(Ga) are proven to switch phases within a realistic temperature range.

Charge self-consistent many-body corrections using optimized projected localized orbitals.

In order for methods combining ab initio density-functional theory and many-body techniques to become routinely used, a flexible, fast, and easy-to-use implementation is crucial. We present an implementation of a general charge self-consistent scheme based on projected localized orbitals in the projector augmented wave framework in the Vienna Ab Initio Simulation Package. We give a detailed description on how the projectors are optimally chosen and how the total energy is calculated. We benchmark our implementation in combination with dynamical mean-field theory: first we study the charge-transfer insulator NiO using a Hartree-Fock approach to solve the many-body Hamiltonian. We address the advantages of the optimized against non-optimized projectors and furthermore find that charge self-consistency decreases the dependence of the spectral function-especially the gap-on the double counting. Second, using continuous-time quantum Monte Carlo we study a monolayer of SrVO3, where strong orbital polarization occurs due to the reduced dimensionality. Using total-energy calculation for structure determination, we find that electronic correlations have a non-negligible influence on the position of the apical oxygens, and therefore on the thickness of the single SrVO3 layer.

Mapping of the catalytic site of CHO-t-PA and the t-PA variant BM 06.022 by synthetic inhibitors and substrates.

BM 06.022 is a t-PA deletion variant that is produced as inactive inclusion bodies in Escherichia coli and transformed into the native form by an in vitro refolding process. Until now, no X-ray and NMR structures of BM 06.022 were available. Therefore a detailed kinetic analysis of the hydrolysis of peptide substrates and of the inhibition by several benzamidine-derived inhibitors was carried out in order to assess that the active site region of the protease domain of BM 06.022 is correctly structured in comparison with t-PA. Our data reveal that the single-chain as well as the two-chain form of BM 06.022 and native t-PA are similar in catalytic and in inhibitor binding properties. This indicates that the active site and the highly complex rearrangement of t-PA upon cleavage of the Arg275-Ile276 bond are maintained in BM 06.022.

Maximally localized Wannier functions in LaMnO3 within PBE + U, hybrid functionals and partially self-consistent GW: an efficient route to construct ab initio tight-binding parameters for eg perovskites.

Using the newly developed VASP2WANNIER90 interface we have constructed maximally localized Wannier functions (MLWFs) for the e(g) states of the prototypical Jahn-Teller magnetic perovskite LaMnO(3) at different levels of approximation for the exchange-correlation kernel. These include conventional density functional theory (DFT) with and without the additional on-site Hubbard U term, hybrid DFT and partially self-consistent GW. By suitably mapping the MLWFs onto an effective e(g) tight-binding (TB) Hamiltonian we have computed a complete set of TB parameters which should serve as guidance for more elaborate treatments of correlation effects in effective Hamiltonian-based approaches. The method-dependent changes of the calculated TB parameters and their interplay with the electron-electron (el-el) interaction term are discussed and interpreted. We discuss two alternative model parameterizations: one in which the effects of the el-el interaction are implicitly incorporated in the otherwise 'noninteracting' TB parameters and a second where we include an explicit mean-field el-el interaction term in the TB Hamiltonian. Both models yield a set of tabulated TB parameters which provide the band dispersion in excellent agreement with the underlying ab initio and MLWF bands.

First-principles study of Cu2ZnSnS4 and the related band offsets for photovoltaic applications.

First-principles calculations of the band offsets between Cu(2)ZnSnS(4) (CZTS) and XS (X = Cd, Zn) are performed. While the interface dipole contribution for the band offsets is calculated using the Perdew-Burke-Ernzerhof functional, the Heyd-Scuseria-Ernzerhof hybrid functional is employed to introduce the quasiparticle corrections to the band offsets. The calculated conduction band offset between CZTS and CdS is 0.2 eV, validating CdS for the buffer layer of the CZTS solar cell. The small conduction band offset stems from the band gap narrowing of CdS under the interface strain caused by the lattice misfit with CZTS. A large valence band offset over 0.9 eV between CZTS and ZnS indicates that precipitated ZnS is regarded as an inactive insulator phase in CZTS absorbers.