- 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

Trimmed sampling algorithm for the noisy generalized eigenvalue problem

Caleb Hicks and Dean Lee

  • Facility for Rare Isotope Beams and Department of Physics and Astronomy, Michigan State University, East Lansing, Michigan 48824, USA

PDF

Phys. Rev. Research 5, L022001 – Published 3 April, 2023

DOI: https://doi.org/10.1103/PhysRevResearch.5.L022001

Abstract

Solving the generalized eigenvalue problem is a useful method for finding energy eigenstates of large quantum systems. It uses projection onto a set of basis states which are typically not orthogonal. One needs to invert a matrix whose entries are inner products of the basis states, and the process is unfortunately susceptible to even small errors. The problem is especially bad when matrix elements are evaluated using stochastic methods and have significant error bars. In this work, we introduce the trimmed sampling algorithm in order to solve this problem. Using the framework of Bayesian inference, we sample prior probability distributions determined by uncertainty estimates of the various matrix elements and likelihood functions composed of physics-informed constraints. The result is a probability distribution for the eigenvectors and observables which automatically comes with a reliable estimate of the error and performs far better than standard regularization methods. The method should have immediate use for a wide range of applications involving classical and quantum computing calculations of large quantum systems.

View figure in article

Physics Subject Headings (PhySH)

Article Text

Supplemental Material

References (37)

  1. D. L. Hill and J. A. Wheeler, Nuclear constitution and the interpretation of fission phenomena, Phys. Rev. 89, 1102 (1953).
  2. J. J. Griffin and J. A. Wheeler, Collective Motions in Nuclei by the Method of Generator Coordinates, Phys. Rev. 108, 311 (1957).
  3. C. Wa Wong, Generator-coordinate methods in nuclear physics, Phys. Rep. 15, 283 (1975).
  4. P. Magierski, P. H. Heenen, and W. Nazarewicz, Generator-coordinate method study of hexadecapole correlations in superdeformed Hg194, Phys. Rev. C 51, R2880 (1995).
  5. K. Varga and Y. Suzuki, Precise solution of few-body problems with the stochastic variational method on a correlated Gaussian basis, Phys. Rev. C 52, 2885 (1995).
  6. D. Blume and K. M. Daily, Universal relations for a trapped four-fermion system with arbitrary s-wave scattering length, Phys. Rev. A 80, 053626 (2009).
  7. T. Togashi, Y. Tsunoda, T. Otsuka, N. Shimizu, and M. Honma, Novel Shape Evolution in Sn Isotopes from Magic Numbers 50 to 82, Phys. Rev. Lett. 121, 062501 (2018).
  8. D. Frame, R. He, I. Ipsen, D. Lee, D. Lee, and E. Rrapaj, Eigenvector Continuation with Subspace Learning, Phys. Rev. Lett. 121, 032501 (2018).
  9. P. Demol, T. Duguet, A. Ekström, M. Frosini, K. Hebeler, S. König, D. Lee, A. Schwenk, V. Somà, and A. Tichai, Improved many-body expansions from eigenvector continuation, Phys. Rev. C 101, 041302(R) (2020).
  10. S. König, A. Ekström, K. Hebeler, D. Lee, and A. Schwenk, Eigenvector continuation as an efficient and accurate emulator for uncertainty quantification, Phys. Lett. B 810, 135814 (2020).
  11. A. Ekström and G. Hagen, Global Sensitivity Analysis of Bulk Properties of an Atomic Nucleus, Phys. Rev. Lett. 123, 252501 (2019).
  12. R. J. Furnstahl, A. J. Garcia, P. J. Millican, and X. Zhang, Efficient emulators for scattering using eigenvector continuation, Phys. Lett. B 809, 135719 (2020).
  13. A. Sarkar and D. Lee, Convergence of Eigenvector Continuation, Phys. Rev. Lett. 126, 032501 (2021).
  14. A. Quarteroni, A. Manzoni, and F. Negri, Reduced Basis Methods for Partial Differential Equations: An Introduction (Springer, New York, 2015), Vol. 92.
  15. E. Bonilla, P. Giuliani, K. Godbey, and D. Lee, Training and projecting: A reduced basis method emulator for many-body physics, Phys. Rev. C 106, 054322 (2022).
  16. J. A. Melendez, C. Drischler, R. J. Furnstahl, A. J. Garcia, and X. Zhang, Model reduction methods for nuclear emulators, J. Phys. G 49, 102001 (2022).
  17. L. Liu, G. Moir, M. Peardon, S. M. Ryan, C. E. Thomas, P. Vilaseca, J. J. Dudek, R. G. Edwards, B. Joo, and D. G. Richards (Hadron Spectrum), Excited and exotic charmonium spectroscopy from lattice QCD, J. High Energy Phys. 07 (2012) 126.
  18. R. G. Edwards, N. Mathur, D. G. Richards, and S. J. Wallace (Hadron Spectrum), Flavor structure of the excited baryon spectra from lattice QCD, Phys. Rev. D 87, 054506 (2013).
  19. S. Elhatisari, Adiabatic projection method with euclidean time subspace projection, Europhys. J. A 55, 144(2019).
  20. S. Shen, T. A. Lähde, D. Lee, and Ulf-G. Meißner, Wigner SU(4) symmetry, clustering, and the spectrum of C12, Eur. Phys. J. A 57, 276 (2021).
  21. S. Shen, T. A. Lähde, D. Lee, and Ulf-G. Meißner, Emergent geometry and duality in the carbon nucleus, arXiv:2202.13596.
  22. A. Francis, A. A. Agrawal, J. H. Howard, E. Kökcü, and A. F. Kemper, Subspace diagonalization on quantum computers using eigenvector continuation, arXiv:2209.10571.
  23. A. Peruzzo, J. McClean, P. Shadbolt, M.-H. Yung, X.-Q. Zhou, P. J. Love, A. Aspuru-Guzik, and J. L. O'Brien, A variational eigenvalue solver on a photonic quantum processor, Nat. Commun. 5, 4213 (2014).
  24. E. F. Dumitrescu, A. J. McCaskey, G. Hagen, G. R. Jansen, T. D. Morris, T. Papenbrock, R. C. Pooser, D. J. Dean, and P. Lougovski, Cloud Quantum Computing of an Atomic Nucleus, Phys. Rev. Lett. 120, 210501 (2018).
  25. Z. Qian, J. Watkins, G. Given, J. Bonitati, K. Choi, and D. Lee, Demonstration of the rodeo algorithm on a quantum computer, arXiv:2110.07747.
  26. M. Bee-Lindgren, Z. Qian, M. DeCross, N. C. Brown, C. N. Gilbreth, J. Watkins, X. Zhang, and D. Lee, Rodeo algorithm with controlled reversal gates, arXiv:2208.13557.
  27. A. N. Tikhonov, On the stability of inverse problems, Dokl. Akad. Nauk SSSR 39, 195 (1943).
  28. M. I. Alheety and B. M. Golam Kibria, Choosing Ridge parameters in the linear regression model with AR(1) error: A comparative simulation study, Intl. J. Stat. Econ. 7, 10 (2011).
  29. A. Gelman, J.B. Carlin, H.S. Stern, and D.B. Rubin, Bayesian Data Analysis (Chapman and Hall/CRC, Boca Raton, FL, 1995).
  30. D. R. Phillips et al., Get on the BAND Wagon: a Bayesian framework for quantifying model uncertainties in nuclear dynamics, J. Phys. G: Nucl. Part. Phys. 48, 072001 (2021).
  31. W. K. Hastings, Monte Carlo sampling methods using Markov chains and their applications, Biometrika 57, 97 (1970).
  32. A. F. M. Smith and A. E. Gelfand, Bayesian statistics without tears: A sampling-resampling perspective, Amer. Stat. 46, 84 (1992).
  33. W. Jiang and C. Forssén, Bayesian probability updates using sampling/importance resampling: Applications in nuclear theory, Front. Phys. 10, 1058809 (2022).
  34. See Supplemental Material at http://link.aps.org/supplemental/10.1103/PhysRevResearch.5.L022001 for descriptions of the Bose-Hubbard model, additional benchmark calculations, trimmed sampling error estimates, and sensitivity studies.
  35. K. Choi, D. Lee, J. Bonitati, Z. Qian, and J. Watkins, Rodeo Algorithm for Quantum Computing, Phys. Rev. Lett. 127, 040505 (2021).
  36. D. Lee, Lattice simulations for few- and many-body systems, Prog. Part. Nucl. Phys. 63, 117 (2009).
  37. T. A. Lähde and Ulf-G. Meißner, Nuclear Lattice Effective Field Theory: An introduction (Springer, Cham, Switzerland, 2019), Vol. 957.

Outline

Information

Phys. Rev. Research 5, L022001– Published 3 April, 2023
  • Received 14 October 2022
  • Accepted 15 March 2023
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