Three pillars for achieving quantum mechanical molecular dynamics simulations of huge systems: Divide‐and‐conquer, density‐functional tight‐binding, and massively parallel computation. Issue 21 (18th June 2016)

• Record Type: Journal Article
• Authors: Nishizawa, Hiroaki; Nishimura, Yoshifumi; Kobayashi, Masato; Irle, Stephan; Nakai, Hiromi
• Is Part Of: Journal of computational chemistry. Volume 37:Issue 21(2016)
• Page Start: 1983
• Page End: 1992
• Publication Date: 2016-06-18
• Subjects: quantum mechanical molecular dynamics -- linear‐scaling divide‐and‐conquer method -- density‐functional tight‐binding method -- massively parallel computation; Chemistry -- Data processing -- Periodicals; 542.85
• DOI: 10.1002/jcc.24419
• Languages: English
• ISSNs: 0192-8651

Abstract

Abstract : The linear‐scaling divide‐and‐conquer (DC) quantum chemical methodology is applied to the density‐functional tight‐binding (DFTB) theory to develop a massively parallel program that achieves on‐the‐fly molecular reaction dynamics simulations of huge systems from scratch. The functions to perform large scale geometry optimization and molecular dynamics with DC‐DFTB potential energy surface are implemented to the program called DC‐DFTB‐K. A novel interpolation‐based algorithm is developed for parallelizing the determination of the Fermi level in the DC method. The performance of the DC‐DFTB‐K program is assessed using a laboratory computer and the K computer. Numerical tests show the high efficiency of the DC‐DFTB‐K program, a single‐point energy gradient calculation of a one‐million‐atom system is completed within 60 s using 7290 nodes of the K computer. © 2016 Wiley Periodicals, Inc. Abstract : The linear‐scaling divide‐and‐conquer (DC) quantum chemical methodology is applied to the density‐functional tight‐binding (DFTB) theory to develop a massively parallel program called DC‐DFTB‐K that can be routinely applied to on‐the‐fly molecular reaction dynamics simulations of large systems. Numerical tests based on calculations of water clusters in a cubic box show a single‐point energy gradient calculation of a one‐million‐atom system is completed within 60 s using 7290 nodes of the K computer.