Consistent basis sets of triple-zeta valence with polarization quality for main group elements and transition metals from row one to three have been derived for periodic quantum-chemical solid-state calculations with the crystalline-orbital program CRYSTAL. They are based on the def2-TZVP basis sets developed for molecules by the Ahlrichs group. Orbital exponents and contraction coefficients have been modified and reoptimized, to provide robust and stable self-consistant field (SCF) convergence for a wide range of different compounds. We compare results on crystal structures, cohesive energies, and solid-state reaction enthalpies with the modified basis sets, denoted as pob-TZVP, with selected standard basis sets available from the CRYSTAL basis set database. The average deviation of calculated lattice parameters obtained with a selected density functional, the hybrid method PW1PW, from experimental reference is smaller with pob-TZVP than with standard basis sets, in particular for metallic systems. The effects of basis set expansion by diffuse and polarization functions were investigated for selected systems.

Revised versions of our published pob-TZVP Oliveira, D. V. and Bredow, T., J. Comput. Chem., 2013, 34 (6), 451-459.] and unpublished pob-DZVP basis sets, denoted as pob-TZVP-rev2 and pob-DZVP-rev2, have been derived for the elements H Br. It was observed that the pob basis sets suffer from the basis set superposition error (BSSE). In order to reduce this effect, we took into account the counterpoise energy of hydride dimers as an additional parameter in the basis set optimization. The overall performance, portability, and SCF stability of the resulting rev2 basis sets are significantly improved compared to the original pob basis sets.

We extend the previously developed geometrical correction for the inter- and intramolecular basis set superposition error (gCP) to periodic density functional theory (DFT) calculations. We report gCP results compared to those from the standard Boys-Bernardi counterpoise correction scheme and large basis set calculations. The applicability of the method to molecular crystals as the main target is tested for the benchmark set X23. It consists of 23 noncovalently bound crystals as introduced by Johnson et al. (J. Chem. Phys. 2012, 137, 054103) and refined by Tkatchenko et al. (J. Chem. Phys. 2013, 139, 024705). In order to accurately describe long-range electron correlation effects, we use the standard atom-pairwise dispersion correction scheme DFT-D3. We show that a combination of DFT energies with small atom-centered basis sets, the D3 dispersion correction, and the gCP correction can accurately describe van der Waals and hydrogen-bonded crystals. Mean absolute deviations of the X23 sublimation energies can be reduced by more than 70% and 80% for the standard functionals PBE and B3LYP, respectively, to small residual mean absolute deviations of about 2 kcal/mol (corresponding to 13% of the average sublimation energy). As a further test, we compute the interlayer interaction of graphite for varying distances and obtain a good equilibrium distance and interaction energy of 6.75 Å and -43.0 meV/atom at the PBE-D3-gCP/SVP level. We fit the gCP scheme for a recently developed pob-TZVP solid-state basis set and obtain reasonable results for the X23 benchmark set and the potential energy curve for water adsorption on a nickel (110) surface.

A recently developed empirical dispersion correction (Grimme et al., J. Chem. Phys. 2010, 132, 154104) to standard density functional theory (DFT-D3) is implemented in the plane-wave program package VASP. The DFT-D3 implementation is compared with an implementation of the earlier DFT-D2 version (Grimme, J. Comput. Chem. 2004, 25, 1463; Grimme, J. Comput. Chem. 2006, 27, 1787). Summation of empirical pair potential terms is performed over all atom pairs in the reference cell and over atoms in shells of neighboring cells until convergence of the dispersion energy is obtained. For DFT-D3, the definition of coordination numbers has to be modified with respect to the molecular version to ensure convergence. The effect of three-center terms as implemented in the original molecular DFT-D3 version is investigated. The empirical parameters are taken from the original DFT-D3 version where they had been optimized for a reference set of small molecules. As the coordination numbers of atoms in bulk and surfaces are much larger than in the reference compounds, this effect has to be discussed. The results of test calculations for bulk properties of metals, metal oxides, benzene, and graphite indicate that the original parameters are also suitable for solid-state systems. In particular, the interlayer distance in bulk graphite and lattice constants of molecular crystals is considerably improved over standard functionals. With the molecular standard parameters (Grimme et al., J. Chem. Phys. 2010, 132, 154104; Grimme, J. Comput. Chem. 2006, 27, 1787) a slight overbinding is observed for ionic oxides where dispersion should not contribute to the bond. For simple adsorbate systems, such as Xe atoms and benzene on Ag(111), the DFT-D implementations reproduce experimental results with a similar accuracy as more sophisticated approaches based on perturbation theory (Rohlfing and Bredow, Phys. Rev. Lett. 2008, 101, 266106).

The structure, properties and stability of 7 thermodynamically stable, meta-stable and high-pressure Al2O3 polymorphs and the structure and relative stability of 4 Al hydroxides were calculated with periodic hybrid density functional theory and compared with available experimental data.