Portal de Periódicos da CAPES

Lazy skip-lists: An algorithm for fast hybridization-expansion quantum Monte Carlo

2014; American Physical Society; Volume: 90; Issue: 7 Linguagem: Inglês

10.1103/physrevb.90.075149

1550-235X

P. Sémon, Chuck-Hou Yee, Kristjan Haule, A.–M. S. Tremblay,

Quantum many-body systems

The solution of a generalized impurity model lies at the heart of electronic structure calculations with dynamical mean-field theory (DMFT). In the strongly-correlated regime, the method of choice for solving the impurity model is the hybridization expansion continuous time quantum Monte Carlo(CT-HYB). Enhancements to the CT-HYB algorithm are critical for bringing new physical regimes within reach of current computational power. Taking advantage of the fact that the bottleneck in the algorithm is a product of hundreds of matrices, we present optimizations based on the introduction and combination of two concepts of more general applicability: a) skip lists and b) fast rejection of proposed configurations based on matrix bounds. Considering two very different test cases with $d$ electrons, we find speedups of $\sim 25$ up to $\sim 500$ compared to the direct evaluation of the matrix product. Even larger speedups are likely with $f$ electron systems and with clusters of correlated atoms.