Taking a quantum leap in time to solution for simulations of high-Tc superconductors

We present a new quantum cluster algorithm to simulate models of high-Tc superconductors. This algorithm extends current methods with continuous lattice self-energies, thereby removing artificial long-range correlations. This cures the fermionic sign problem in the underlying quantum Monte Carlo solver for large clusters and realistic values of the Coulomb interaction in the entire temperature range of interest. We find that the new algorithm improves time-to-solution by nine orders of magnitude compared to current, state of the art quantum cluster simulations. An efficient implementation is given, which ports to multi-core as well as hybrid CPU-GPU systems. Running on 18,600 nodes on ORNL's Titan supercomputer enables us to compute a converged value of Tc/t = 0.053±0.0014 for a 28 site cluster in the 2D Hubbard model with U/t = 7 at 10% hole doping. Typical simulations on Titan sustain between 9.2 and 15.4 petaflops (double precision measured over full run), depending on configuration and parameters used.

[1]  L. Lucy.  An iterative technique for the rectification of observed distributions , 1974 .

[2]  M. Norman,et al.  The Challenge of Unconventional Superconductivity , 2011, Science.

[3]  Eduardo F. D'Azevedo,et al.  New algorithm to enable 400+ TFlop/s sustained performance in simulations of disorder effects in high- T c superconductors , 2008, HiPC 2008.

[4]  O. Jepsen,et al.  LDA energy bands, low-energy hamiltonians, t′, t″, t⊥ (k), and J⊥ , 1995 .

[5]  Antonios Gonis,et al.  Green functions for ordered and disordered systems , 1992 .

[6]  White,et al.  Sign problem in the numerical simulation of many-electron systems. , 1990, Physical review. B, Condensed matter.

[7]  John M. Levesque,et al.  New algorithm to enable 400+ TFlop/s sustained performance in simulations of disorder effects in high-Tc superconductors , 2008, 2008 SC - International Conference for High Performance Computing, Networking, Storage and Analysis.

[8]  A. B. Migdal,et al.  INTERACTION BETWEEN ELECTRONS AND THE LATTICE VIBRATIONS IN A NORMAL METAL , 1958 .

[9]  Peter Staar.  Dynamical cluster approximation with continuous lattice self-energy , 2013 .

[10]  T. Pruschke,et al.  Quantum cluster theories , 2004, cond-mat/0404055.

[11]  Michael S. Summers,et al.  Fast update algorithm for the quantum Monte Carlo simulation of the Hubbard model , 2009 .

[12]  O. Jepsen,et al.  Calculations of Hubbard U from first-principles , 2006 .

[13]  India,et al.  Dynamical cluster approximation: Nonlocal dynamics of correlated electron systems , 2000 .

[14]  O. K. Andersen,et al.  Muffin-tin orbitals of arbitrary order , 2000 .

[15]  Mila.  Parameters of a Hubbard Hamiltonian to describe superconducting Cu oxides. , 1988, Physical review. B, Condensed matter.

[16]  I. A. Nekrasov,et al.  Full orbital calculation scheme for materials with strongly correlated electrons , 2004, cond-mat/0407359.

[17]  Thomas C. Schulthess,et al.  Efficient non-equidistant FFT approach to the measurement of single- and two-particle quantities in continuous time Quantum Monte Carlo methods , 2012 .

[18]  Martin,et al.  Cuprate parameters from numerical Wannier functions. , 1990, Physical review. B, Condensed matter.

[19]  Mark Jarrell,et al.  Quantum Monte Carlo algorithm for nonlocal corrections to the dynamical mean-field approximation , 2001 .

[20]  C. Marianetti,et al.  Electronic structure calculations with dynamical mean-field theory , 2005, cond-mat/0511085.

[21]  G. M. Éliashberg,et al.  Interactions between electrons and lattice vibrations in a superconductor , 1960 .

[22]  Emanuel Gull,et al.  Submatrix updates for the continuous-time auxiliary-field algorithm , 2010, 1010.3690.

[23]  T. Maier,et al.  Structure of the pairing interaction in the two-dimensional Hubbard model. , 2005, Physical review letters.

[24]  Physikalisches Institut,et al.  Fictive impurity models: An alternative formulation of the cluster dynamical mean-field method , 2003 .

[25]  Fye,et al.  Monte Carlo method for magnetic impurities in metals. , 1986, Physical review letters.

[26]  P R C Kent,et al.  Systematic study of d-wave superconductivity in the 2D repulsive Hubbard model. , 2005, Physical review letters.

[27]  K. Cheng.  Theory of Superconductivity , 1948, Nature.

[28]  N. Marzari,et al.  Maximally-localized Wannier Functions: Theory and Applications , 2011, 1112.5411.

[29]  Thomas C. Schulthess,et al.  DCA$^+$: Incorporating self-consistently a continuous momentum self-energy in the Dynamical Cluster Approximation , 2013 .

[30]  Emanuel Gull,et al.  Continuous-time auxiliary-field Monte Carlo for quantum impurity models , 2008, 0802.3222.

[31]  William H. Richardson,et al.  Bayesian-Based Iterative Method of Image Restoration , 1972 .

[32]  Thomas C. Schulthess,et al.  Toward First Principles Electronic Structure Simulations of Excited States and Strong Correlations in Nano- and Materials Science , 2010, 2010 ACM/IEEE International Conference for High Performance Computing, Networking, Storage and Analysis.