Eigenvector Continuation with Subspace Learning.

A common challenge faced in quantum physics is finding the extremal eigenvalues and eigenvectors of a Hamiltonian matrix in a vector space so large that linear algebra operations on general vectors are not possible. There are numerous efficient methods developed for this task, but they generally fail when some control parameter in the Hamiltonian matrix exceeds some threshold value. In this Letter we present a new technique called eigenvector continuation that can extend the reach of these methods. The key insight is that while an eigenvector resides in a linear space with enormous dimensions, the eigenvector trajectory generated by smooth changes of the Hamiltonian matrix is well approximated by a very low-dimensional manifold. We prove this statement using analytic function theory and propose an algorithm to solve for the extremal eigenvectors. We benchmark the method using several examples from quantum many-body theory.

[1]  D. Sorensen Numerical methods for large eigenvalue problems , 2002, Acta Numerica.

[2]  Willy Govaerts,et al.  Numerical computation of bifurcations in large equilibrium systems in matlab , 2014, J. Comput. Appl. Math..

[3]  Andrew G. Glen,et al.  APPL , 2001 .

[4]  Dean Lee,et al.  Nuclear lattice simulations using symmetry-sign extrapolation , 2015, 1502.06787.

[5]  Dean Lee,et al.  Nucleon-deuteron scattering using the adiabatic projection method , 2016, 1603.02333.

[6]  Dean Lee,et al.  Ab initio alpha–alpha scattering , 2015, Nature.

[7]  F. Iachello Analytic description of critical point nuclei in a spherical-axially deformed shape phase transition. , 2001, Physical review letters.

[8]  Ning Li,et al.  Ab initio Calculations of the Isotopic Dependence of Nuclear Clustering. , 2017, Physical review letters.

[9]  Morten Hjorth-Jensen,et al.  An advanced course in computational nuclear physics , 2017 .

[10]  H. Sebastian Seung,et al.  The Manifold Ways of Perception , 2000, Science.

[11]  Fred Wubs,et al.  Numerical Bifurcation Methods and their Application to Fluid Dynamics: Analysis beyond Simulation , 2012 .

[12]  Dean Lee Lattice simulations for few- and many-body systems , 2008, 0804.3501.

[13]  Ning Li,et al.  Nuclear Binding Near a Quantum Phase Transition. , 2016, Physical review letters.

[14]  G. C. Knollman,et al.  Quantum Cell Model for Bosons , 1963 .

[15]  James Demmel,et al.  Continuation of Invariant Subspaces in Large Bifurcation Problems , 2008, SIAM J. Sci. Comput..

[16]  Sandro Stringari,et al.  Theory of ultracold atomic Fermi gases , 2007, 0706.3360.

[17]  C. Lanczos An iteration method for the solution of the eigenvalue problem of linear differential and integral operators , 1950 .

[18]  Daniel Kressner,et al.  Subspace Acceleration for Large-Scale Parameter-Dependent Hermitian Eigenproblems , 2015, SIAM J. Matrix Anal. Appl..

[19]  Y. Saad Numerical Methods for Large Eigenvalue Problems , 2011 .

[20]  Yousef Saad Analysis of Subspace Iteration for Eigenvalue Problems with Evolving Matrices , 2016, SIAM J. Matrix Anal. Appl..