Efficient non-equidistant FFT approach to the measurement of single- and two-particle quantities in continuous time Quantum Monte Carlo methods

Continuous time cluster solvers allow us to measure single- and two-particle Greens functions in the Matsubara frequency domain with unprecedented accuracy. Currently, the usage of the two-particle functions is limited due to a lack of an efficient measurement method that can deal with the random times of the vertices. In this paper, we show how the Non-equidistant Fast Fourier Transform (NFFT) algorithm can be modified in order to obtain a very efficient measurement algorithm. For the single particle case, we propose a delayed-NFFT (d-NFFT) scheme, which reduces the arithmetical operations from O(N log(N)) in NFFT to O(N), a huge improvement compared to the standard O(N2) of the Non-equidistant Discrete Fourier Transform (NDFT), currently used in most continuous time cluster solvers. For the two-particle case, we discuss how the NFFT can be applied to measure the two-particle Greens functions and how to exploit its properties to further optimize the NFFT. We then apply these algorithms to the half-filled 2D Hubbard model at U/t = 8 in order to study the anti-ferromagnetic transition. In particular, we confirm the logarithmic decay of the Neel-temperatures versus cluster-sizes in accordance with the Mermin-Wagner theorem.

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

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

[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]  Matthias Troyer,et al.  Continuous-time solver for quantum impurity models. , 2005, Physical review letters.

[5]  Stefan Kunis,et al.  Using NFFT 3---A Software Library for Various Nonequispaced Fast Fourier Transforms , 2009, TOMS.

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

[7]  M. Troyer,et al.  Continuous-time Monte Carlo methods for quantum impurity models , 2010, 1012.4474.

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

[9]  Philipp Werner,et al.  Hybridization expansion impurity solver: General formulation and application to Kondo lattice and two-orbital models , 2006 .

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

[11]  Steven G. Johnson,et al.  The Design and Implementation of FFTW3 , 2005, Proceedings of the IEEE.

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

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