Linear scaling 3D fragment method for large-scale electronic structure calculations [electronic resource].
- Berkeley, Calif. : Lawrence Berkeley National Laboratory, 2008. and Oak Ridge, Tenn. : Distributed by the Office of Scientific and Technical Information, U.S. Dept. of Energy.
- Additional Creators:
- Lawrence Berkeley National Laboratory and United States. Department of Energy. Office of Scientific and Technical Information
- Restrictions on Access:
- Free-to-read Unrestricted online access
- We present a new linearly scaling three-dimensional fragment (LS3DF) method for large scale ab initio electronic structure calculations. LS3DF is based on a divide-and-conquer approach, which incorporates a novel patching scheme that effectively cancels out the artificial boundary effects due to the subdivision of the system. As a consequence, the LS3DF program yields essentially the same results as direct density functional theory (DFT) calculations. The fragments of the LS3DF algorithm can be calculated separately with different groups of processors. This leads to almost perfect parallelization on tens of thousands of processors. After code optimization, we were able to achieve 35.1 Tflop/s, which is 39% of the theoretical speed on 17,280 Cray XT4 processor cores. Our 13,824-atom ZnTeO alloy calculation runs 400 times faster than a direct DFT calculation, even presuming that the direct DFT calculation can scale well up to 17,280 processor cores. These results demonstrate the applicability of the LS3DF method to material simulations, the advantage of using linearly scaling algorithms over conventional O(N³) methods, and the potential for petascale computation using the LS3DF method.
- Published through SciTech Connect., 07/11/2008., "lbnl-603e", Supercomputing 2008., Wang, Lin-Wang; Strohmaier, Erich; Shan, HongZhang; Bailey, David; Meza, Juan; Lee, Byounghak; Zhao, Zhengji., and Computational Research Division
- Funding Information:
View MARC record | catkey: 14654167