ABSTRACT
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.
- See, for example, M. Rasetti, ed., The Hubbard Model: Recent Results, World Scientific, Singapore, 1991; A. Montorsi, ed., The Hubbard Model, World Scientific, Singapore, 1992.Google Scholar
- M. Machida, S. Yamada, Y. Ohashi, and H. Matsumoto, Novel superfluidity in a trapped gas of Fermi atoms with repulsive interaction loaded on an optical lattice, Phys. Rev. Lett., 93, 200402, 2004.Google ScholarCross Ref
- M. Rigol, A. Muramatsu, G. G. Batrouni, and R. T. Scalettar, Local quantum criticality in confined fermions on optical lattices, Phys. Rev. Lett., 91, 130403, 2003.Google ScholarCross Ref
- For example, see, E. Dagotto, Correlated electrons in high-temperature superconductors, Rev. Mod. Phys., 66, 763, 1994.Google ScholarCross Ref
- S. Yamada, T. Imamura, M. Machida, 16.447 TFlops and 159-Billion-dimensional Exact-diagonalization for Trapped Fermion-Hubbard Model on the Earth Simulator, Proc. of SC2005, 2005. http://sc05.supercomputing.org/schedule/pdf/pap188.pdf Google ScholarDigital Library
- The Earth Simulator Center, http://www.es.jamstec.go.jp/esc/eng/index.htmlGoogle Scholar
- S. Shingu, H. Takahara, H. Fuchigami, M. Yamada, Y. Tsuda, W. Ohfuchi, Y. Sasaki, K. Kobayashi, T. Hagiwara, S. Habata, M. Yokokawa, H. Itoh and K. Otsuka, A 26.58 Tflops Global Atmospheric Simulation with the Spectral Transform Method on the Earth Simulator, Proc. of SC2002, 2002. http://sc-2002.org/paperpdfs/pap.pap331.pdf Google ScholarDigital Library
- H. Sakagami, H. Murai, Y. Seo and M. Yokokawa, 14.9 TFLOPS Three-dimensional Fluid Simulation for Fusion Science with HPF on the Earth Simulator, Proc. of SC2002, 2002. http://sc-2002.org/paperpdfs/pap.pap147.pdf Google ScholarDigital Library
- M. Yokokawa, K. Itakura, A. Uno, T. Ishihara and Y. Kaneda, 16.4 Tflops Direct Numerical Simulation of Turbulence by Fourier Spectral Method on the Earth Simulator, Proc. of SC2002, 2002. http://sc-2002.org/paperpdfs/pap.pap273.pdf Google ScholarDigital Library
- D. Komatitsch, S. Tsuboi, C, Ji and J. Tromp, A 14.6 billion degrees of freedom, 5 teraflops, 2.5 terabyte earthquake simulation on the Earth Simulator, Proc. of SC2003, 2003. http://www.sc-conference.org/sc2003/paperpdfs/pap124.pdf Google ScholarDigital Library
- A. Kageyama, M. Kameyama, S. Fujihara, M. Yoshida, M. Hyodo, and Y. Tsuda, A 15.2 TFlops Simulation of Geodynamo on the Earth Simulator, Proc. of SC2004, 2004. http://www.sc-conference.org/sc2004/schedule/pdfs/pap234.pdf Google ScholarDigital Library
- A. V. Knyazev, Preconditioned eigensolvers - An oxymoron?, Electronic Transactions on Numerical analysis, Vol. 7, 104--123, 1998.Google Scholar
- A. V. Knyazev, Toward the optimal eigensolver: Locally optimal block preconditioned conjugate gradient method, SIAM J. Sci. Comput., 23, 517--541, 2001. Google ScholarDigital Library
- H. Uehara, M. Tamura, K. Itakura and M. Yokokawa, MPI Performance Evaluation on the Earth Simulator (in Japanese), Transactions on High Performance Computing Systems, 44(SIG 1 (HPS 6)), 2003, 24--34.Google Scholar
- J. E. Hoffman, E. W. Hudson, K. M. Lang, V. Madhavan, H. Eisaki, S. Uchida, and J. C. Davis, A four unit cell periodic pattern of Quasi-particle states surrounding vertex cores in Bi2Sr2CaCuO8+Δ, SCIENCE, Vol. 295, 2002Google Scholar
- M Machida, and T. Koyama, Friedel oscillation in charge profile and position dependent screening around a superconducting vortex core, Phys. Rev. Lett., Vol. 90, 077003, 2003.Google ScholarCross Ref
- G. Levy, M. Kugler, A. A. Manuel, and Ø. Fischer, Fourfold structure of vertex-cores states in Bi2Sr2CaCuO8+Δ, Phys. Rev. Lett., Vol. 95, 257005, 2005.Google ScholarCross Ref
- M Machida, and T. Koyama, Structure of a quantized vortex near the BCS-BEC crossover in an atomic fermi gas, Phys. Rev. Lett., Vol. 94, 140401, 2005.Google ScholarCross Ref
- ScaLAPACK, http://www.netlib.org/scalapack/scalapack_home.htmlGoogle Scholar
- K. Naono, Y. Yamamoto, M. Igai, and H. Hirayama: "High performance implementation of tridiagonalization on the SR8000", In Proc. of the High Performance Computing in Asia-Pacific Region (HPC-ASIA2000), Beijing, pp.206--219, 2000.Google ScholarCross Ref
- C. Bischof, X. Sun, and B. Lang, "Parallel tridiagonalization through two-step band reduction", In Proc. of Scalable High Performance Computing Conference, pp.23--27. IEEE, 1994.Google ScholarCross Ref
- Y. Yamamoto, Performance and Accuracy of Algorithms for Computing the Eigenvalues of Real Symmetric Matrices on Cache-based Multiprocessors (in Japanese), Transactions on High Performance Computing Systems, 46(SIG 3 (ACS8)), 2005, 81--91.Google Scholar
- C. Bischof, C. V. Loan, "The WY Representation for Products of Householder Matrices", SIAM J. Sci. Stat. Comput. Vol.8, No.1, 1987. Google ScholarDigital Library
- P. Bientinest, I. S. Dhillon, R. V. Geijn: "A Parallel Eigensolver for Dense Symmetric Matrices Based on Multiple Relatively Robust Representations", SIAM J. Sci. Comput. Vol.27, No.1, pp.43--66, 2005. Google ScholarDigital Library
- P. Alpatov, et al.: "PLAPACK: Parallel Linear Algebra Package," in Proc. of the SIAM Parallel Processing Conference, 1997 http://www.cs.utexas.edu/~plapack/Google Scholar
- E. Breitmoser, A. G. Sunderland, "A performance study of the PLAPACK and ScaLAPACK Eigensolvers on HPCx for the standard problem", Technical Report from the HPCx Consortium, 2004. http://www.hpcx.ac.uk/research/hpc/HPCxTR0406.pdfGoogle Scholar
Index Terms
- High-performance computing for exact numerical approaches to quantum many-body problems on the earth simulator
Recommendations
Quantum information processing with quantum zeno many-body dynamics
We show how the quantum Zeno effect can be exploited to control quantum many-bodydynamics for quantum information and computation purposes. In particular, we considera one dimensional array of three level systems interacting via a nearest-...
Comments