Simulation of Liquids with the Tight-Binding Density-Functional Approach and Improved Atomic Charges

Theoretical description of liquids, especially liquid water, is an ongoing subject with important implications in various domains such as homogeneous catalysis; solvation of molecular, ionic, and biomolecular species; and reactivity. Various formalisms exist to describe liquids, each one displaying its own balance between accuracy and computational cost that defines its range of applications. The present article revisits the ability of the density-functional-based tight-binding (SCC-DFTB) approach to model liquids by focusing on liquid water and liquid benzene under ambient conditions. To do so, we benchmark a recent correction for the SCC-DFTB atomic charges that allows for a drastic improvement of the pair radial distribution functions of liquid water as compared to both experimental data and density-functional theory results performed in the generalized-gradient approximation. We also report the coupling of the deMonNano and i-PI codes to perform path-integral molecular dynamics. This allows us to rationalize the impact of nuclear quantum effects on the SCC-DFTB description of liquid water. This study evidences the rather good ability of SCC-DFTB to describe liquid water and liquid benzene. As the first example of application, we also present results for a benzene molecule solvated in water with the perspectives of further studies devoted to solvent/water interfaces.