Large-Scale Quantum-Mechanical Molecular Dynamics Simulations Using Density-Functional Tight-Binding Combined with the Fragment Molecular Orbital Method.

The fully analytic gradient is developed for density-functional tight-binding (DFTB) combined with the fragment molecular orbital (FMO) method (FMO-DFTB). The response terms arising from the coupling of the electronic state to the embedding potential are derived, and the gradient accuracy is demonstrated on water clusters and a polypeptide. The radial distribution functions (RDFs) obtained with FMO-DFTB are found to be similar to those from conventional DFTB, while the computational cost is greatly reduced; for 256 water molecules one molecular dynamics (MD) step takes 73.26 and 0.68 s with full DFTB and FMO-DFTB, respectively, showing a speed-up factor of 108. FMO-DFTB/MD is applied to 100 ps MD simulations of liquid hydrogen halides and is found to reproduce experimental RDFs reasonably well.