Extremely Scalable Algorithm for 108-atom Quantum Material Simulation on the Full System of the K Computer

An extremely scalable linear-algebraic algorithm was developed for quantum material simulation (electronic state calculation) with 108 atoms or 100-nm-scale materials. The mathematical foundation is generalized shifted linear equations ((zB — A)x = b), instead of conventional generalized eigenvalue equations. The method has a highly parallelizable mathematical structure. The benchmark shows an extreme strong scaling and a qualified time-to-solution on the full system of the K computer. The method was demonstrated in a real material research for ultra-flexible (organic) devices, key devices of next-generation Internet-of-Things (IoT) products. The present paper shows that an innovative scalable algorithm for a real research can appear by the co-design among application, algorithm and architecture.

[1]  Richard Dronskowski,et al.  Crystal orbital Hamilton populations (COHP): energy-resolved visualization of chemical bonding in solids based on density-functional calculations , 1993 .

[2]  Kohn,et al.  Density functional and density matrix method scaling linearly with the number of atoms. , 1996, Physical review letters.

[3]  D. Sánchez-Portal,et al.  The SIESTA method for ab initio order-N materials simulation , 2001, cond-mat/0111138.

[4]  Andreas Frommer,et al.  BiCGStab(ℓ) for Families of Shifted Linear Systems , 2003, Computing.

[5]  Chris-Kriton Skylaris,et al.  Introducing ONETEP: linear-scaling density functional simulations on parallel computers. , 2005, The Journal of chemical physics.

[6]  Taisuke Ozaki.  O(N) Krylov-subspace method for large-scale ab initio electronic structure calculations , 2006 .

[7]  R. Takayama,et al.  Linear algebraic calculation of the Green’s function for large-scale electronic structure theory , 2006 .

[8]  M. J. Gillan,et al.  Order-N first-principles calculations with the conquest code , 2007, Comput. Phys. Commun..

[9]  Tomohiro Sogabe,et al.  ON A WEIGHTED QUASI-RESIDUAL MINIMIZATION STRATEGY FOR SOLVING COMPLEX SYMMETRIC SHIFTED LINEAR SYSTEMS , 2008 .

[10]  Tomohiro Sogabe,et al.  Shifted Conjugate-Orthogonal-Conjugate-Gradient Method and Its Application to Double Orbital Extended Hubbard Model(Condensed matter: electronic structure and electrical, magnetic, and optical properties) , 2008 .

[11]  T. Fujiwara,et al.  On a weighted quasi-residual minimization strategy of the QMR method for solving complex symmetric shifted linear systems , 2009, 0902.2614.

[12]  Tetsuya Sakurai,et al.  Filter diagonalization of shell-model calculations , 2010, 1008.0434.

[13]  S. Louie,et al.  GW method with the self-consistent Sternheimer equation , 2009, 0912.3087.

[14]  S. Yamamoto,et al.  Efficient and accurate linear algebraic methods for large-scale electronic structure calculations with nonorthogonal atomic orbitals , 2011, 1101.4768.

[15]  Taisuke Boku,et al.  First-principles calculations of electron states of a silicon nanowire with 100,000 atoms on the K computer , 2011, 2011 International Conference for High Performance Computing, Networking, Storage and Analysis (SC).

[16]  Tomohiro Sogabe,et al.  Solution of generalized shifted linear systems with complex symmetric matrices , 2012, J. Comput. Phys..

[17]  Kenichi Miura,et al.  The design of ultra scalable MPI collective communication on the K computer , 2012, Computer Science - Research and Development.

[18]  T. Fujiwara,et al.  An order-N electronic structure theory with generalized eigenvalue equations and its application to a ten-million-atom system , 2012, Journal of physics. Condensed matter : an Institute of Physics journal.

[19]  Tatsunori Tanaka,et al.  Ten-million-atom electronic structure calculations on the K computer with a massively parallel order-N theory , 2012, 1210.1531.

[20]  Yohei Akiyama,et al.  Novel linear algebraic theory and one-hundred-million-atom electronic structure calculation on the K computer , 2013, 1307.3218.

[21]  A Marek,et al.  The ELPA library: scalable parallel eigenvalue solutions for electronic structure theory and computational science , 2014, Journal of physics. Condensed matter : an Institute of Physics journal.

[22]  Tomohiro Sogabe,et al.  Novel linear algebraic theory and one-hundred-million-atom quantum material simulations on the K computer , 2014, CSE 2014.

[23]  Tomoya Ono,et al.  Numerical solver for first-principles transport calculation based on real-space finite-difference method. , 2015, Physical review. E, Statistical, nonlinear, and soft matter physics.

[24]  Hiroto Imachi,et al.  One-hundred-nm-scale electronic structure and transport calculations of organic polymers on the K computer , 2016, 1603.09616.

[25]  Hiroto Imachi,et al.  Hybrid Numerical Solvers for Massively Parallel Eigenvalue Computations and Their Benchmark with Electronic Structure Calculations , 2015, J. Inf. Process..

[26]  M. Kaltenbrunner,et al.  Ultraflexible organic photonic skin , 2016, Science Advances.