Transferable density functional tight binding for carbon, hydrogen, nitrogen, and oxygen: Application to shock compression

A new parameterization for density functional tight binding (DFTB) theory, lanl31, has been developed for molecules containing carbon, hydrogen, nitrogen, and oxygen. Optimal values for the Hubbard Us, on-site energies, and the radial dependences of the bond integrals and repulsive potentials were determined by numerical optimization using simulated annealing to a modest database of ab initio-calculated atomization energies and interatomic forces. The transferability of the optimized DFTB parameterization has been assessed using the CHNO subset of the QM-9 database [R. Ramakrishnan et al., Sci. Data 1, 140022 (2014)]. These analyses showed that the errors in the atomization energies and interatomic forces predicted by our model are small and in the vicinity of the differences between density functional theory calculations with different basis sets and exchange-correlation functionals. Good correlations between the molecular dipole moments and HOMO-LUMO gaps predicted by lanl31 and the QM-9 data set are al...