Abstract
We present an accurate parametrization of density functional theory augmented with an empirical correction term to describe properly pi-stacking interactions in nucleic acids. The approach is based on the popular Perdew-Burke-Ernzerhof (PBE), Becke-Perdew (BP), and hybrid Becke-Lee-Yang-Parr (B3LYP) density functionals augmented by a classical London C(6)R(-6) dispersion term. The novelty of our implementation lies in the accurate tuning of the empirical parameters, included in the [formula: see text] damping function, to reproduce high-level post Hartree-Fock calculations. In particular, we present sets of parameters and the needed code to correct the PBE, BP, and B3LYP results from the Turbomole and ADF packages in connection with basis sets of double and triple zeta quality. The developed approach is validated by comparison with the JSCH-2005 benchmark and with best quality stacking energies reported in the literature for the stacking of H-bonded nucleic acids base pairs.