Gaussian basis sets for accurate calculations on molecular systems in gas and condensed phases

Abstract

(2007). Gaussian basis sets for accurate calculations on molecular systems in gas and condensed phases. The Journal of Chemical Physics, 127 (11) Gaussian basis sets for accurate calculations on molecular systems in gas and condensed phases AbstractWe present a library of Gaussian basis sets that has been specifically optimized to perform accurate molecular calculations based on density functional theory. It targets a wide range of chemical environments, including the gas phase, interfaces, and the condensed phase. These generally contracted basis sets, which include diffuse primitives, are obtained minimizing a linear combination of the total energy and the condition number of the overlap matrix for a set of molecules with respect to the exponents and contraction coefficients of the full basis. Typically, for a given accuracy in the total energy, significantly fewer basis functions are needed in this scheme than in the usual split valence scheme, leading to a speedup for systems where the computational cost is dominated by diagonalization. More importantly, binding energies of hydrogen bonded complexes are of similar quality as the ones obtained with augmented basis sets, i.e., have a small (down to 0.2 kcal/mol) basis set superposition error, and the monomers have dipoles within 0.1 D of the basis set limit. However, contrary to typical augmented basis sets, there are no near linear dependencies in the basis, so that the overlap matrix is always well conditioned, also, in the condensed phase. The basis can therefore be used in first principles molecular dynamics simulations and is well suited for linear scaling calculations. We present a library of Gaussian basis sets that has been specifically optimized to perform accurate molecular calculations based on density functional theory. It targets a wide range of chemical environments, including the gas phase, interfaces, and the condensed phase. These generally contracted basis sets, which include diffuse primitives, are obtained minimizing a linear combination of the total energy and the condition number of the overlap matrix for a set of molecules with respect to the exponents and contraction coefficients of the full basis. Typically, for a given accuracy in the total energy, significantly fewer basis functions are needed in this scheme than in the usual split valence scheme, leading to a speedup for systems where the computational cost is dominated by diagonalization. More importantly, binding energies of hydrogen bonded complexes are of similar quality as the ones obtained with augmented basis sets, i.e., have a small ͑down to 0.2 kcal/ mol͒ basis set superposition error, and the monomers have dipoles within 0.1 D of the basis set limit. However, contrary to typical augmented basis sets, there are no near linear dependencies in the basis, so that the overlap matrix is always well conditioned, also, in the condensed phase. The basis can therefore be used in first principles molecular dynamics simulations and is well suited for linear scaling calculations. Gaus...

Wetting of water on hexagonal boron nitride@Rh(111): a QM/MM model based on atomic charges derived for nano-structured substrates

“…40,41 Single-ζ valence plus polarization (SZVP) basis sets are used to describe the valence electrons of Rh, whereas basis sets of at least double-ζ quality (DZVP) are used for the other elements. 42 Norm-conserving pseudo-potentials [43][44][45] are employed to approximate the interactions between valence electrons and atomic cores. The exchange and correlation potential is given by the revised Perdew-Burke-Ernzerhof (revPBE) 46 functional.…”

Microscopic properties of liquid water from combined ab initio molecular dynamics and energy decomposition studies

“…In the linearly dependent limit, where the basis functions on one fragment can exactly mimic functions on another fragment, this ceases to be the case 34 . This is not an issue for the AO basis sets used routinely in quantum chemistry and AIMD simulations 55 .…”

“…Efforts towards a full inear scaling algorithm are underway [118][119][120] . Here, the orbitals were represented by an accurate triple-ζ basis set with two set of polarization functions (TZV2P) 55 , while a density cutoff of 320 Ry was used for the charge density.…”

“…The MOs in our ALMO EDA calculations were represented by a triple-ζ Gaussian basis set with two sets of polarization functions (TZV2P) optimized specifically for molecular systems 55 . A very high density cutoff of 1000 Ry was used to describe the electron density.…”

Ligand/cluster/support catalytic complexes in heterogeneous ultrananocatalysis: NO oxidation on Ag3/MgO(100)

“…22 We choose pseudopotentials derived by Goedecker, Teter and Hutter (GTH) to describe the core electrons 23 and DZVP basis sets 24 to represent the DFT Kohn-Sham orbitals. The cut-off for the auxiliary plane wave representation of the density is 300 Ry.…”