Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:cond-mat/0306077 · Physical review. B, Condensed matter

Short-time-evolved wave functions for solving quantum many-body problems

Orion Ciftja, Siu A. Chin

Published Oct 6, 2003• 20 citations
arXivHTMLPDF

Abstract

The exact ground state of a strongly interacting quantum many-body system can be obtained by evolving a trial state with finite overlap with the ground state to infinite imaginary time. In many cases, since the convergence is exponential, the system converges essentially to the exact ground state in a relatively short time. Thus a short-time evolved wave function can be an excellent approximation to the exact ground state. Such a short-time-evolved wave function can be obtained by factorizing, or splitting, the evolution operator to high order. However, for the imaginary time Schr\"odinger equation, which contains an irreversible diffusion kernel, all coefficients, or time steps, must be positive. (Negative time steps would require evolving the diffusion process backward in time, which is impossible.) Heretofore, only second-order factorization schemes can have all positive coefficients, but without further iterations, these cannot be used to evolve the system long enough to be close to the exact ground state. In this work, we use a newly discovered fourth-order positive factorization scheme which requires knowing both the potential and its gradient. We show that the resulting fourth-order wave function alone, without further iterations, gives an excellent description of strongly interacting quantum systems such as liquid ${}^{4}\mathrm{He},$ comparable to the best variational results in the literature. This suggests that such a fourth-order wave function can be used to study the ground state of diverse quantum many-body systems, including Bose-Einstein condensates and Fermi systems.

DOI
10.1103/physrevb.68.134510

Related papers

Fourth-order diffusion Monte Carlo algorithms for solving quantum many-body problems

Harald Forbert, Siu A. Chin

By decomposing the important sampled imaginary time Schr\"odinger evolution operator to fourth order with positive coefficients, we derived a number of distinct fourth-order diffusion Monte Carlo algorithms. These sophisticated algorithms require higher derivatives of the drift velocity and local energy and are more complicated to program. However, they allowed very large time steps to be used, converged faster with lesser correlations, and virtually eliminated the step size error. We demonstrated the effectiveness of these quartic algorithms by solving for the ground-state energy and radial density distribution of bulk liquid helium.

Physical review. B, Condensed matter•Mar 21, 2001•36 citations

Modified splitting methods for Gross-Pitaevskii systems modelling Bose-Einstein condensates: Time evolution and ground state computation

Mechthild Thalhammer, Gregor Thalhammer

The year 2025 marks the 100 and 30 years anniversaries of the discovery of Bose--Einstein condensation and its successful experimental realisation. Inspired by these important research achievements, a conceptually simple approach is proposed to facilitate reliable and efficient numerical simulations. The structure of the underlying systems of coupled Gross--Pitaevskii equations suggests the use of optimised high-order operator splitting methods for dynamical evolution and ground state computation. A second-order barrier, however, prevents the applicability of standard operator splitting methods for both, time evolution as well as imaginary time propagation. An innovative alternative approach accomplishes the design of novel modified operator splitting methods that remain stable under moderate smallness assumptions on the time increments. The core idea is to incorporate commutators of the defining differential and nonlinear multiplication operators, since this permits to fulfill the basic stability requirement of positive method coefficients. Further improvements with respect to convergence at the targeted precision arise from automatic adjustments of the time stepsizes by an inexpensive local error control. The presented numerical experiments confirm the favourable performance of a specific fourth-order modified operator splitting method. Amongst others, it is demonstrated that the excellent mass and energy conservation in long-term evolutions, intrinsic attributes of geometric numerical integrators for Hamiltonian systems, is maintained for a sensible variation of the time stepsizes. Moreover, the benefits of adaptive higher-order approximations in ground state computations are illustrated.

arXiv:2601.19838•Jan 27, 2026•0 citations

A FOURTH ORDER DIFFUSION MONTE CARLO ALGORITHM FOR SOLVING QUANTUM MANY-BODY PROBLEMS

Harald Forbert, Siu A. Chin

We derive a fourth-order diffusion Monte Carlo algorithm for solving quantum many-body problems. The method uses a factorization of the imaginary time propagator in terms of the usual local energy E and Langevin operators L as well as an additional pseudo-potential consisting of the double commutator [EL, [L, EL]]. A new factorization of the propagator of the Fokker-Planck equation enables us to implement the Langevin algorithm to the necessary fourth order. We achieve this by the addition of correction terms to the drift steps and the use of a position-dependent Gaussian random walk. We show that in the case of bulk liquid helium the systematic step size errors are indeed fourth order over a wide range of step sizes.

arXiv:nucl-th/0009068•International Journal of Modern Physics B•May 10, 2001•8 citations

Exact two-body expansion of the many-particle wave function

David A. Mazziotti

Progress toward the solution of the strongly correlated electron problem has been stymied by the exponential complexity of the wave function. Previous work established an exact two-body exponential product expansion for the ground-state wave function. By developing a reduced density-matrix analog of Dalgarno-Lewis perturbation theory, we prove here that (i) the two-body exponential product expansion is rapidly and globally convergent with each operator representing an order of a renormalized perturbation theory, (ii) the energy of the expansion converges quadratically near the solution, and (iii) the expansion is exact for both ground and excited states. The two-body expansion offers a reduced parametrization of the many-particle wave function as well as the two-particle reduced density matrix with potential applications on both conventional and quantum computers for the study of strongly correlated quantum systems. We demonstrate the result with the exact solution of the contracted Schr\"odinger equation for the molecular chains ${\mathrm{H}}_{4}$ and ${\mathrm{H}}_{5}$.

arXiv:2010.02191•Physical Review A•Sep 21, 2020•19 citations

Bang-bang algorithms for quantum many-body ground states: A tensor network exploration

Ruoshui Wang, Timothy H. Hsieh, Guifré Vidal

We use matrix product techniques to investigate the performance of two algorithms for obtaining the ground state of a quantum many-body Hamiltonian $H={H}_{A}+{H}_{B}$ in infinite systems. The first algorithm is a generalization of the quantum approximate optimization algorithm and uses a quantum computer to evolve an initial product state into an approximation of the ground state of $H$ by alternating between ${H}_{A}$ and ${H}_{B}$. We show for the one-dimensional (1D) quantum Ising model that the accuracy in representing a gapped ground state improves exponentially with the number of alternations. The second algorithm is the variational imaginary time ansatz, which uses a classical computer to simulate the ground state via alternating imaginary time steps with ${H}_{A}$ and ${H}_{B}$. We find for the 1D quantum Ising model that an accurate approximation to the ground state is obtained with a total imaginary time $\ensuremath{\tau}$ that grows only logarithmically with the inverse energy gap $1/\mathrm{\ensuremath{\Delta}}$ of $H$. This is much faster than imaginary time evolution by $H$, which would require $\ensuremath{\tau}\ensuremath{\sim}1/\mathrm{\ensuremath{\Delta}}$.

arXiv:2208.00271•Physical review. B./Physical review. B•Nov 17, 2022•5 citations

Exact factorization ansatz of the ground state many-body wave function: A case study of quartic potential problem

Zhiyuan Yin, Yi Hua Deng, Xinyao Wang +1 more

We propose an explicit factorized formula for the ground state many-body wave function of a model quartic potential problem in two and three dimensions, using our recently developed method for solving Schrödinger equations. Our novel formula has a clear real-space structure: it is exactly factorized into three parts, including a non-integer pre-exponential power factor, dominant decaying terms on the exponent, and a modulator function, which is of minor importance. This result is a generalization of our previously obtained exact formula for the 1-body wave function. As in the 1-body case, here we show that for many-body ground wave functions that cannot achieve variable separation, our method is still far more efficient than the conventional basis expansion method both in representing the wave function and in calculating the energy. The exact factorization ansatz proposed in the present case study shall provide important insight into the general structure of many-body wave functions.

The Journal of Chemical Physics•Aug 28, 2025•1 citation
The exact ground state of a strongly interacting quantum many-body system can be obtained by evolving a trial state with finite overlap with the ground state to infinite imaginary time. In many cases, since the convergence is exponential, the system converges essentially to the exact ground state in a relatively short time. Thus a short-time evolved wave function can be an excellent approximation to the exact ground state. Such a short-time-evolved wave function can be obtained by factorizing, or splitting, the evolution operator to high order. However, for the imaginary time Schr"odinger equation, which contains an irreversible diffusion kernel, all coefficients, or time steps, must be positive. (Negative time steps would require evolving the diffusion process backward in time, which is impossible.) Heretofore, only second-order factorization schemes can have all positive coefficients, but without further iterations, these cannot be used to evolve the system long enough to be close to the exact ground state. In this work, we use a newly discovered fourth-order positive factorization scheme which requires knowing both the potential and its gradient. We show that the resulting fourth-order wave function alone, without further iterations, gives an excellent description of strongly interacting quantum systems such as liquid 4He,{}^{4}\mathrm{He},4He, comparable to the best variational results in the literature. This suggests that such a fourth-order wave function can be used to study the ground state of diverse quantum many-body systems, including Bose-Einstein condensates and Fermi systems.
By decomposing the important sampled imaginary time Schr\"odinger evolution operator to fourth order with positive coefficients, we derived a number of distinct fourth-order diffusion Monte Carlo algorithms. These sophisticated algorithms require higher derivatives of the drift velocity and local energy and are more complicated to program. However, they allowed very large time steps to be used, converged faster with lesser correlations, and virtually eliminated the step size error. We demonstrated the effectiveness of these quartic algorithms by solving for the ground-state energy and radial density distribution of bulk liquid helium.
The year 2025 marks the 100 and 30 years anniversaries of the discovery of Bose--Einstein condensation and its successful experimental realisation. Inspired by these important research achievements, a conceptually simple approach is proposed to facilitate reliable and efficient numerical simulations. The structure of the underlying systems of coupled Gross--Pitaevskii equations suggests the use of optimised high-order operator splitting methods for dynamical evolution and ground state computation. A second-order barrier, however, prevents the applicability of standard operator splitting methods for both, time evolution as well as imaginary time propagation. An innovative alternative approach accomplishes the design of novel modified operator splitting methods that remain stable under moderate smallness assumptions on the time increments. The core idea is to incorporate commutators of the defining differential and nonlinear multiplication operators, since this permits to fulfill the basic stability requirement of positive method coefficients. Further improvements with respect to convergence at the targeted precision arise from automatic adjustments of the time stepsizes by an inexpensive local error control. The presented numerical experiments confirm the favourable performance of a specific fourth-order modified operator splitting method. Amongst others, it is demonstrated that the excellent mass and energy conservation in long-term evolutions, intrinsic attributes of geometric numerical integrators for Hamiltonian systems, is maintained for a sensible variation of the time stepsizes. Moreover, the benefits of adaptive higher-order approximations in ground state computations are illustrated.
We derive a fourth-order diffusion Monte Carlo algorithm for solving quantum many-body problems. The method uses a factorization of the imaginary time propagator in terms of the usual local energy E and Langevin operators L as well as an additional pseudo-potential consisting of the double commutator [EL, [L, EL]]. A new factorization of the propagator of the Fokker-Planck equation enables us to implement the Langevin algorithm to the necessary fourth order. We achieve this by the addition of correction terms to the drift steps and the use of a position-dependent Gaussian random walk. We show that in the case of bulk liquid helium the systematic step size errors are indeed fourth order over a wide range of step sizes.
Progress toward the solution of the strongly correlated electron problem has been stymied by the exponential complexity of the wave function. Previous work established an exact two-body exponential product expansion for the ground-state wave function. By developing a reduced density-matrix analog of Dalgarno-Lewis perturbation theory, we prove here that (i) the two-body exponential product expansion is rapidly and globally convergent with each operator representing an order of a renormalized perturbation theory, (ii) the energy of the expansion converges quadratically near the solution, and (iii) the expansion is exact for both ground and excited states. The two-body expansion offers a reduced parametrization of the many-particle wave function as well as the two-particle reduced density matrix with potential applications on both conventional and quantum computers for the study of strongly correlated quantum systems. We demonstrate the result with the exact solution of the contracted Schr\"odinger equation for the molecular chains H4{\mathrm{H}}_{4}H4​ and H5{\mathrm{H}}_{5}H5​.
We use matrix product techniques to investigate the performance of two algorithms for obtaining the ground state of a quantum many-body Hamiltonian H=HA+HBH={H}_{A}+{H}_{B}H=HA​+HB​ in infinite systems. The first algorithm is a generalization of the quantum approximate optimization algorithm and uses a quantum computer to evolve an initial product state into an approximation of the ground state of HHH by alternating between HA{H}_{A}HA​ and HB{H}_{B}HB​. We show for the one-dimensional (1D) quantum Ising model that the accuracy in representing a gapped ground state improves exponentially with the number of alternations. The second algorithm is the variational imaginary time ansatz, which uses a classical computer to simulate the ground state via alternating imaginary time steps with HA{H}_{A}HA​ and HB{H}_{B}HB​. We find for the 1D quantum Ising model that an accurate approximation to the ground state is obtained with a total imaginary time $\ensuremath{\tau}$ that grows only logarithmically with the inverse energy gap $1/\mathrm{\ensuremath{\Delta}}$ of HHH. This is much faster than imaginary time evolution by HHH, which would require $\ensuremath{\tau}\ensuremath{\sim}1/\mathrm{\ensuremath{\Delta}}$.
We propose an explicit factorized formula for the ground state many-body wave function of a model quartic potential problem in two and three dimensions, using our recently developed method for solving Schrödinger equations. Our novel formula has a clear real-space structure: it is exactly factorized into three parts, including a non-integer pre-exponential power factor, dominant decaying terms on the exponent, and a modulator function, which is of minor importance. This result is a generalization of our previously obtained exact formula for the 1-body wave function. As in the 1-body case, here we show that for many-body ground wave functions that cannot achieve variable separation, our method is still far more efficient than the conventional basis expansion method both in representing the wave function and in calculating the energy. The exact factorization ansatz proposed in the present case study shall provide important insight into the general structure of many-body wave functions.