We show that numerical atomic orbital basis sets that are variationally optimized for specific compounds are biased and not fully transferable to other compounds. The situation improves when the basis sets are optimized for several compounds and a compromise is made. We demonstrate this for basis sets for compounds containing H, O, N, Al, Cl, and Cu atoms, using the SIESTA code. Reaction energies are compared to benchmark calculations with the ADF code. Our double-zeta basis sets give a root mean square error of 0.2 eV. With triple-zeta basis sets, we obtain a root mean square error of only 0.1 eV. Finally, for the case of mixed phase systems, we propose that optimizing basis sets for gaseous molecules and then using them also for bulk compounds, is better than the other way around.

• Received 27 May 2011

DOI:https://doi.org/10.1103/PhysRevB.85.035108

This article is available under the terms of the Creative Commons Attribution 3.0 License. Further distribution of this work must maintain attribution to the author(s) and the published article’s title, journal citation, and DOI.

©2012 American Physical Society