High-performance computing for exact numerical approaches to quantum many-body problems on the earth simulator

In order to study quantum many-body problems, we develop two matrix diagonalization codes, which solve only the ground state and all quantum states, respectively. The target model in both codes is the Hubbard model with confinement potential which describes an atomic Fermi gas loaded on an optical lattice and partly High-Tc cuprate superconductor. For the former code, we obtain 18.692TFlops (57% of the peak) as the best performance on the Earth Simulator when calculating the ground state of 100-billion dimensional matrix. From these large-scale calculations, we find atomic-scale inhomogeneous superfluid state which is now a challenging subject for physicists. For the latter code, we succeed in solving the matrix whose dimension is 375,000 with locally 24.6TFlops (75% of the peak). The calculations reveal that a change from Schrodinger's cat to classical like one can be controlled by tuning the interaction. This is a marked contrast to the general concept.

[1]  Oystein Fischer,et al.  Four-fold structure of vortex core states in Bi2Sr2CaCu2O8 (Bi2212) , 2005 .

[2]  Mitsuo Yokokawa,et al.  16.4-Tflops Direct Numerical Simulation of Turbulence by a Fourier Spectral Method on the Earth Simulator , 2002, ACM/IEEE SC 2002 Conference (SC'02).

[3]  P. Alpatov,et al.  PLAPACK Parallel Linear Algebra Package Design Overview , 1997, ACM/IEEE SC 1997 Conference (SC'97).

[4]  Andrew Knyazev,et al.  Preconditioned Eigensolvers - an Oxymoron? , 1998 .

[5]  Masanori Kameyama,et al.  A 15.2 TFlops Simulation of Geodynamo on the Earth Simulator , 2004, Proceedings of the ACM/IEEE SC2004 Conference.

[6]  Andrew V. Knyazev,et al.  Toward the Optimal Preconditioned Eigensolver: Locally Optimal Block Preconditioned Conjugate Gradient Method , 2001, SIAM J. Sci. Comput..

[7]  M. Machida,et al.  Structure of a quantized vortex near the BCS-BEC crossover in an atomic Fermi gas. , 2005, Physical review letters.

[8]  Chen Ji,et al.  A 14.6 billion degrees of freedom, 5 teraflops, 2.5 terabyte earthquake simulation on the Earth Simulator , 2003, ACM/IEEE SC 2003 Conference (SC'03).

[9]  M. Machida,et al.  Novel superfluidity in a trapped gas of Fermi atoms with repulsive interaction loaded on an optical lattice. , 2004, Physical review letters.

[10]  Xiaobai Sun,et al.  Parallel tridiagonalization through two-step band reduction , 1994, Proceedings of IEEE Scalable High Performance Computing Conference.

[11]  Hiroshi Takahara,et al.  A 26.58 Tflops Global Atmospheric Simulation with the Spectral Transform Method on the Earth Simulator , 2002, ACM/IEEE SC 2002 Conference (SC'02).

[12]  M Machida,et al.  Friedel oscillation in charge profile and position dependent screening around a superconducting vortex core. , 2003, Physical review letters.

[13]  M. Rasetti.  The Hubbard model : recent results , 1991 .

[14]  Yoshiki Seo,et al.  14.9 TFLOPS Three-Dimensional Fluid Simulation for Fusion Science with HPF on the Earth Simulator , 2002, ACM/IEEE SC 2002 Conference (SC'02).

[15]  Elena Breitmoser,et al.  A performance study of the PLAPACK and ScaLAPACK Eigensolvers on HPCx for the standard problem , 2003 .

[16]  Ken Naono,et al.  High performance implementation of tridiagonalization on the SR8000 , 2000, Proceedings Fourth International Conference/Exhibition on High Performance Computing in the Asia-Pacific Region.

[17]  Toshiyuki Imamura,et al.  16.447 TFlops and 159-Billion-dimensional Exact-diagonalization for Trapped Fermion-Hubbard Model on the Earth Simulator , 2005, ACM/IEEE SC 2005 Conference (SC'05).

[18]  H. Eisaki,et al.  A Four Unit Cell Periodic Pattern of Quasi-Particle States Surrounding Vortex Cores in Bi2Sr2CaCu2O8+δ , 2002, Science.

[19]  Christian H. Bischof,et al.  The WY representation for products of householder matrices , 1985, PPSC.

[20]  R. Scalettar,et al.  Local quantum criticality in confined fermions on optical lattices. , 2003, Physical review letters.

[21]  Robert A. van de Geijn,et al.  PLAPACK: Parallel Linear Algebra Package , 1997, PPSC.

[22]  Robert A. van de Geijn,et al.  A Parallel Eigensolver for Dense Symmetric Matrices Based on Multiple Relatively Robust Representations , 2005, SIAM J. Sci. Comput..

[23]  E. Dagotto.  Correlated electrons in high-temperature superconductors , 1993, cond-mat/9311013.