Efficient Parallel Linear Scaling Construction of the Density Matrix for Born鈥揙ppenheimer Molecular Dynamics

详细信息    查看全文

• 作者：S. M. Mniszewski ; M. J. Cawkwell ; M. E. Wall ; J. Mohd-Yusof ; N. Bock ; T. C. Germann ; A. M. N. Niklasson
• 刊名：Journal of Chemical Theory and Computation
• 出版年：2015
• 出版时间：October 13, 2015
• 年：2015
• 卷：11
• 期：10
• 页码：4644-4654
• 全文大小：822K
• ISSN：1549-9626

文摘

We present an algorithm for the calculation of the density matrix that for insulators scales linearly with system size and parallelizes efficiently on multicore, shared memory platforms with small and controllable numerical errors. The algorithm is based on an implementation of the second-order spectral projection (SP2) algorithm [Niklasson, A. M. N. Phys. Rev. B 2002, 66, 155115] in sparse matrix algebra with the ELLPACK-R data format. We illustrate the performance of the algorithm within self-consistent tight binding theory by total energy calculations of gas phase poly(ethylene) molecules and periodic liquid water systems containing up to 15,000 atoms on up to 16 CPU cores. We consider algorithm-specific performance aspects, such as local vs nonlocal memory access and the degree of matrix sparsity. Comparisons to sparse matrix algebra implementations using off-the-shelf libraries on multicore CPUs, graphics processing units (GPUs), and the Intel many integrated core (MIC) architecture are also presented. The accuracy and stability of the algorithm are illustrated with long duration Born鈥揙ppenheimer molecular dynamics simulations of 1000 water molecules and a 303 atom Trp cage protein solvated by 2682 water molecules.