- Physical Review Journals
All JournalsPhysics Magazine

Physical Review Research

Reuse & Permissions

It is not necessary to obtain permission to reuse this article or its components as it is available under the terms of the Creative Commons Attribution 4.0 International license. This license permits unrestricted use, distribution, and reproduction in any medium, provided attribution to the author(s) and the published article's title, journal citation, and DOI are maintained. Please note that some figures may have been included with permission from other third parties. It is your responsibility to obtain the proper permission from the rights holder directly for these figures.

Export citation

Export citation

Choose format for download:

Download Citation
  • Letter
  • Open Access

Denoising of imaginary time response functions with Hankel projections

Yang Yu1,*, Alexander F. Kemper2,†, Chao Yang3,‡, and Emanuel Gull1,§

  • *Contact author: umyangyu@umich.edu
  • †Contact author: akemper@ncsu.edu
  • ‡Contact author: cyang@lbl.gov
  • §Contact author: egull@umich.edu
PDF

Phys. Rev. Research 6, L032042 – Published 23 August, 2024

DOI: https://doi.org/10.1103/PhysRevResearch.6.L032042

Abstract

Imaginary-time response functions of finite-temperature quantum systems are often obtained with methods that exhibit stochastic or systematic errors. Reducing these errors comes at a large computational cost—in quantum Monte Carlo simulations, the reduction of noise by a factor of two incurs a simulation cost of a factor of four. In this paper, we relate certain imaginary-time response functions to an inner product on the space of linear operators on Fock space. We then show that data with noise typically does not respect the positive definiteness of its associated Gramian. The Gramian has the structure of a Hankel matrix. As a method for denoising noisy data, we introduce an alternating projection algorithm that finds the closest positive definite Hankel matrix consistent with noisy data. We test our methodology at the example of fermion Green's functions for continuous-time quantum Monte Carlo data and show remarkable improvements of the error, reducing noise by a factor of up to 20 in practical examples. We argue that Hankel projections should be used whenever finite-temperature imaginary-time data of response functions with errors is analyzed, be it in the context of quantum Monte Carlo, quantum computing, or in approximate semianalytic methodologies.

View figure in article

Physics Subject Headings (PhySH)

Article Text

Supplemental Material

References (55)

  1. A. Altland and B. D. Simons, Condensed Matter Field Theory (Cambridge University Press, Cambridge, 2010).
  2. R. Blankenbecler, D. J. Scalapino, and R. L. Sugar, Monte Carlo calculations of coupled boson-fermion systems. I, Phys. Rev. D 24, 2278 (1981).
  3. J. E. Hirsch and R. M. Fye, Monte Carlo method for magnetic impurities in metals, Phys. Rev. Lett. 56, 2521 (1986).
  4. N. V. Prokof'ev and B. V. Svistunov, Polaron problem by diagrammatic quantum Monte Carlo, Phys. Rev. Lett. 81, 2514 (1998).
  5. A. N. Rubtsov, V. V. Savkin, and A. I. Lichtenstein, Continuous-time quantum Monte Carlo method for fermions, Phys. Rev. B 72, 035122 (2005).
  6. P. Werner, A. Comanac, L. de' Medici, M. Troyer, and A. J. Millis, Continuous-time solver for quantum impurity models, Phys. Rev. Lett. 97, 076405 (2006).
  7. N. Prokof'ev and B. Svistunov, Bold diagrammatic Monte Carlo technique: When the sign problem is welcome, Phys. Rev. Lett. 99, 250201 (2007).
  8. A. Bazavov, D. Toussaint, C. Bernard, J. Laiho, C. DeTar, L. Levkova, M. B. Oktay, S. Gottlieb, U. M. Heller, J. E. Hetrick, P. B. Mackenzie, R. Sugar, and R. S. Van de Water, Nonperturbative QCD simulations with 2+1 Flavors of improved staggered quarks, Rev. Mod. Phys. 82, 1349 (2010).
  9. E. Gull, A. J. Millis, A. I. Lichtenstein, A. N. Rubtsov, M. Troyer, and P. Werner, Continuous-time Monte Carlo methods for quantum impurity models, Rev. Mod. Phys. 83, 349 (2011).
  10. G. Cohen, E. Gull, D. R. Reichman, and A. J. Millis, Taming the dynamical sign problem in real-time evolution of quantum many-body problems, Phys. Rev. Lett. 115, 266802 (2015).
  11. D. Neuhauser, R. Baer, and D. Zgid, Stochastic self-consistent second-order green's function method for correlation energies of large electronic systems, J. Chem. Theory Comput. 13, 5396 (2017).
  12. R. Rossi, Determinant diagrammatic Monte Carlo algorithm in the thermodynamic limit, Phys. Rev. Lett. 119, 045701 (2017).
  13. B. Winograd, E. Gull, and D. Zgid, Error propagation in the fully self-consistent stochastic second-order green's function method, arXiv:2009.10048.
  14. M. Qin, T. Schäfer, S. Andergassen, P. Corboz, and E. Gull, The Hubbard model: A computational perspective, Annu. Rev. Condens. Matter Phys. 13, 275 (2022).
  15. J. Gubernatis, N. Kawashima, and P. Werner, Quantum Monte Carlo Methods: Algorithms for Lattice Models (Cambridge University Press, Cambridge, 2016).
  16. H. J. Rothe, Lattice Gauge Theories: An Introduction (World Scientific Publishing Company, Singapore, 2012).
  17. G. D. Mahan, Many Particle Physics (Springer, New York, 2010).
  18. L. Mühlbacher and E. Rabani, Real-time path integral approach to nonequilibrium many-body quantum systems, Phys. Rev. Lett. 100, 176403 (2008).
  19. P. Werner, T. Oka, and A. J. Millis, Diagrammatic Monte Carlo simulation of nonequilibrium systems, Phys. Rev. B 79, 035320 (2009).
  20. A. Erpenbeck, E. Gull, and G. Cohen, Quantum Monte Carlo method in the steady state, Phys. Rev. Lett. 130, 186301 (2023).
  21. M. Jarrell and J. Gubernatis, Bayesian inference and the analytic continuation of imaginary-time quantum Monte Carlo data, Phys. Rep. 269, 133 (1996).
  22. J. Fei, C.-N. Yeh, and E. Gull, Nevanlinna analytical continuation, Phys. Rev. Lett. 126, 056402 (2021).
  23. L. Zhang and E. Gull, Minimal pole representation and controlled analytic continuation of Matsubara response functions, Phys. Rev. B 110, 035154 (2024).
  24. C. Charles, E. J. Gustafson, E. Hardt, F. Herren, N. Hogan, H. Lamm, S. Starecheski, R. S. Van de Water, and M. L. Wagman, Simulating z2 lattice gauge theory on a quantum computer, Phys. Rev. E 109, 015307 (2024).
  25. A. F. Kemper, C. Yang, and E. Gull, Denoising and extension of response functions in the time domain, Phys. Rev. Lett. 132, 160403 (2024).
  26. M. J. Hyrkäs, D. Karlsson, and R. van Leeuwen, Cutting rules and positivity in finite temperature many-body theory, J. Phys. A: Math. Theor. 55, 335301 (2022).
  27. L. Boehnke, H. Hafermann, M. Ferrero, F. Lechermann, and O. Parcollet, Orthogonal polynomial representation of imaginary-time green's functions, Phys. Rev. B 84, 075145 (2011).
  28. N. J. Higham, Computing a nearest symmetric positive semidefinite matrix, Linear Algebra Appl. 103, 103 (1988).
  29. M. G. Eberle and M. C. Maciel, Finding the closest toeplitz matrix, Comput. Appl. Math. 22, 1 (2003).
  30. J. P. Boyle and R. L. Dykstra, A method for finding projections onto the intersection of convex sets in Hilbert spaces, in Advances in Order Restricted Statistical Inference, edited by R. Dykstra, T. Robertson, and F. T. Wright (Springer, New York, 1986), pp. 28–47.
  31. See https://github.com/CQMP/HankelProjection.git for the C++ implementation of the algorithm, documentation, and examples.
  32. G. Kotliar, S. Y. Savrasov, K. Haule, V. S. Oudovenko, O. Parcollet, and C. A. Marianetti, Electronic structure calculations with dynamical mean-field theory, Rev. Mod. Phys. 78, 865 (2006).
  33. P. Werner and A. J. Millis, Hybridization expansion impurity solver: General formulation and application to Kondo lattice and two-orbital models, Phys. Rev. B 74, 155107 (2006).
  34. B. Bauer, L. D. Carr, H. G. Evertz, A. Feiguin, J. Freire, S. Fuchs, L. Gamper, J. Gukelberger, E. Gull, S. Guertler, A. Hehn, R. Igarashi, S. V. Isakov, D. Koop, P. N. Ma, P. Mates, H. Matsuo, O. Parcollet, G. Pawłowski, J. D. Picon et al., The ALPS project release 2.0: open source software for strongly correlated systems, J. Stat. Mech. (2011) P05001.
  35. P. Seth, I. Krivenko, M. Ferrero, and O. Parcollet, TRIQS/CTHYB: A continuous-time quantum Monte Carlo hybridisation expansion solver for quantum impurity problems, Comput. Phys. Commun. 200, 274 (2016).
  36. H. Shinaoka, E. Gull, and P. Werner, Continuous-time hybridization expansion quantum impurity solver for multi-orbital systems with complex hybridizations, Comput. Phys. Commun. 215, 128 (2017).
  37. M. Wallerberger, A. Hausoel, P. Gunacker, A. Kowalski, N. Parragh, F. Goth, K. Held, and G. Sangiovanni, w2dynamics: Local one- and two-particle quantities from dynamical mean field theory, Comput. Phys. Commun. 235, 388 (2019).
  38. W. Metzner and D. Vollhardt, Correlated lattice fermions in d=∞ dimensions, Phys. Rev. Lett. 62, 324 (1989).
  39. E. Müller-Hartmann, Correlated fermions on a lattice in high dimensions, Z. Phys. B 74, 507 (1989).
  40. A. Georges, G. Kotliar, W. Krauth, and M. J. Rozenberg, Dynamical mean-field theory of strongly correlated fermion systems and the limit of infinite dimensions, Rev. Mod. Phys. 68, 13 (1996).
  41. N. Trivedi and M. Randeria, Deviations from Fermi-liquid behavior above Tc in 2D short coherence length superconductors, Phys. Rev. Lett. 75, 312 (1995).
  42. T. A. Maier, M. Jarrell, T. Prushke, and M. Hettler, Quantum cluster theories, Rev. Mod. Phys. 77, 1027 (2005).
  43. E. Gull, P. Werner, O. Parcollet, and M. Troyer, Continuous-time auxiliary-field Monte Carlo for quantum impurity models, Europhys. Lett. 82, 57003 (2008).
  44. G. Pick, Über die beschränkungen analytischer funktionen durch vorgegebene funktionswerte, Math. Ann. 78, 270 (1917).
  45. J. Fei, C.-N. Yeh, D. Zgid, and E. Gull, Analytical continuation of matrix-valued functions: Carathéodory formalism, Phys. Rev. B 104, 165111 (2021).
  46. L. Ying, Analytic continuation from limited noisy Matsubara data, J. Comput. Phys. 469, 111549 (2022).
  47. Z. Huang, E. Gull, and L. Lin, Robust analytic continuation of Green's functions via projection, pole estimation, and semidefinite relaxation, Phys. Rev. B 107, 075151 (2023).
  48. L. Pollet, Recent developments in quantum Monte Carlo simulations with applications for cold gases, Rep. Prog. Phys. 75, 094501 (2012).
  49. G. Stefanucci, Y. Pavlyukh, A.-M. Uimonen, and R. van Leeuwen, Diagrammatic expansion for positive spectral functions beyond GW: Application to vertex corrections in the electron gas, Phys. Rev. B 90, 115134 (2014).
  50. Y. Núñez Fernández, M. Jeannin, P. T. Dumitrescu, T. Kloss, J. Kaye, O. Parcollet, and X. Waintal, Learning Feynman diagrams with tensor trains, Phys. Rev. X 12, 041018 (2022).
  51. A. Erpenbeck, W.-T. Lin, T. Blommel, L. Zhang, S. Iskakov, L. Bernheimer, Y. Núñez-Fernández, G. Cohen, O. Parcollet, X. Waintal, and E. Gull, Tensor train continuous time solver for quantum impurity models, Phys. Rev. B 107, 245135 (2023).
  52. See Supplemental Material at http://link.aps.org/supplemental/10.1103/PhysRevResearch.6.L032042 for the data of the figures.
  53. M. H. Hettler, A. N. Tahvildar-Zadeh, M. Jarrell, T. Pruschke, and H. R. Krishnamurthy, Nonlocal dynamical correlations of strongly interacting electron systems, Phys. Rev. B 58, R7475 (1998).
  54. M. Jarrell, T. Maier, C. Huscroft, and S. Moukouri, Quantum Monte Carlo algorithm for nonlocal corrections to the dynamical mean-field approximation, Phys. Rev. B 64, 195130 (2001).
  55. S. Li and E. Gull, Magnetic and charge susceptibilities in the half-filled triangular lattice Hubbard model, Phys. Rev. Res. 2, 013295 (2020).

Outline

Information

Phys. Rev. Research 6, L032042– Published 23 August, 2024
  • Received 19 March 2024
  • Accepted 30 July 2024
Crossmark - Check for updates button
Creative Commons Logo

Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article's title, journal citation, and DOI.

Published by the American Physical Society

Sign In to Your Journals Account

Filter

Filter

Article Lookup

Enter a citation