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:1104.2627 · Physical Review A

Simulating strongly correlated multiparticle systems in a truncated Hilbert space

Thomas Ernst, D. W. Hallwood, Jake Gulliksen, Hans‐Dieter Meyer, Joachim Brand

Published Aug 16, 2011• 24 citations
arXivHTMLPDF

Abstract

Representing a strongly interacting multiparticle wave function in a finite product basis leads to errors. Simple rescaling of the contact interaction can preserve the low-lying energy spectrum and long-wavelength structure of wave functions in one-dimensional systems and thus correct for the basis set truncation error. The analytic form of the rescaling is found for a two-particle system where the rescaling is exact. A detailed comparison between finite Hilbert space calculations and exact results for up to five particles show that rescaling can significantly improve the accuracy of numerical calculations in various external potentials. In addition to ground-state energies, the low-lying excitation spectrum, density profile, and correlation functions are studied. The results give a promising outlook for numerical simulations of trapped ultracold atoms.

DOI
10.1103/physreva.84.023623

Related papers

Are smooth pseudopotentials a good choice for representing short-range interactions?

Péter Jeszenszki, Ali Alavi, Joachim Brand

When seeking a numerical representation of a quantum-mechanical multiparticle problem it is tempting to replace a singular short-range interaction by a smooth finite-range pseudopotential. Finite basis set expansions, e.g., in Fock space, are then guaranteed to converge exponentially. The need to faithfully represent the artificial length scale of the pseudopotential, however, places a costly burden on the basis set. Here we discuss scaling relations for the required size of the basis set and demonstrate the basis set convergence on the example of a two-dimensional system of few fermions with short-range $s$-wave interactions in a harmonic trapping potential. In particular we show that the number of harmonic oscillator basis functions needed to reach a regime of exponential convergence for a Gaussian pseudopotential scales with the fourth power of the pseudopotential length scale, which can be improved to quadratic scaling when the basis functions are rescaled appropriately. Numerical examples for three fermions with up to a few hundred single-particle basis functions are presented and implications for the feasibility of accurate numerical multiparticle simulations of interacting ultracold-atom systems are discussed.

arXiv:1812.06521•Physical Review A•Mar 11, 2019•6 citations

On the renormalization of contact interactions for the configuration-interaction method in two-dimensions

Massimo Rontani, G Eriksson, S Åberg +1 more

Abstract The configuration interaction (CI) method for calculating the exact eigenstates of a quantum-mechanical few-body system is problematic when applied to particles interacting through contact forces. In dimensions higher than one the approach fails due to the pathology of the Dirac δ -potential, making it impossible to reach convergence by gradually increasing the size of the Hilbert space. However, this problem may be cured in a rather simple manner by renormalizing the strength of the contact potential when diagonalizing in a truncated Hilbert space. One hereby relies on the comparison of the CI results to the two-body ground-state energy obtained by the exact solution of the Schrödinger equation for a regularized contact interaction. We discuss here a scheme that provides cutoff-independent few-body physical observables. The method is applied to a few-body system of ultracold atoms confined by a two-dimensional harmonic oscillator.

arXiv:0810.4305•Journal of Physics B Atomic Molecular and Optical Physics•Feb 14, 2017•20 citations

Configuration interaction approach to the few-body problem in a two-dimensional harmonic trap with contact interaction

Massimo Rontani, Sarah Reimann

The configuration interaction (CI) method for calculating the exact eigenstates of a quantum-mechanical few-body system is problematic when applied to contact interactions between the particles. In two and three dimensions, the approach fails due to the pathology of the Dirac delta-potential, making it impossible to reach convergence by gradually increasing the size of the Hilbert space. However, for practical applications this problem may be cured in a rather simple manner, by renormalizing the STRENGTH of the contact potential, which must be diagonalized in a TRUNCATED Hilbert space. The procedure relies on the comparison of CI energies and wave functions with those obtained by the exact solution of the two-body Schrodinger equation for the regularized contact interaction. The rather simple scheme, while keeping the numerical procedures still elementary, nevertheless provides both cut-off-independent few-body physical observables and an estimate of the error of the CI calculation.

•Oct 23, 2008•3 citations

Path-Integral Renormalization Group Method for Numerical Study on Ground States of Strongly Correlated Electronic Systems

Tsuyoshi Kashima, Masatoshi Imada

A new efficient numerical algorithm for interacting fermion systems is proposed and examined in detail. The ground state is expressed approximately by a linear combination of numerically chosen basis states in a truncated Hilbert space. Two procedures lead to a better approximation. The first is a numerical renormalization, which optimizes the chosen basis and projects onto the ground state within the fixed dimension, $L$, of the Hilbert space. The second is an increase of the dimension of the truncated Hilbert space, which enables the linear combination to converge to a better approximation. The extrapolation $L\to\infty$ after the convergence removes the approximation error systematically. This algorithm does not suffer from the negative sign problem and can be applied to systems in any spatial dimension and arbitrary lattice structure. The efficiency is tested and the implementation explained for two-dimensional Hubbard models where Slater determinants are employed as chosen basis. Our results with less than 400 chosen basis indicate good accuracy within the errorbar of the best available results as those of the quantum Monte Carlo for energy and other physical quantities.

arXiv:cond-mat/0104140•Journal of the Physical Society of Japan•Aug 15, 2001•89 citations

Systematic interpolatory ansatz for one-dimensional polaron systems

Lindgren, E. J., R. E. Barfknecht, Nikolaj Thomas Zinner

We explore a new variational principle for studying one-dimensional quantum systems in a trapping potential. We focus on the Fermi polaron problem, where a single distinguishable impurity interacts through a contact potential with a background of identical fermions. We can accurately describe this system at arbitrary finite repulsion by constructing a truncated basis containing states at both the limits of zero and infinite repulsion. We show how to construct this basis and how to obtain energies, density matrices and correlation functions, and provide results both for a harmonic well and a double well for various particle numbers. The results are compared both with matrix product states methods and with the analytical result for two particles in a harmonic well.

arXiv:1908.02849•Aug 7, 2019•0 citations

Importance-Truncated No-Core Shell Model for Fermionic Many-Body Systems

Helena Spies

The exact solution of quantum mechanical many-body problems is only possible for few particles. Therefore, numerical methods were developed in the fields of quantum physics and quantum chemistry for larger particle numbers. Configuration Interaction (CI) methods or the No-Core Shell Model (NCSM) allow ab initio calculations for light and intermediate-mass nuclei, without resorting to phenomenology. An extension of the NCSM is the Importance-Truncated No-Core Shell Model, which uses an a priori selection of the most important basis states. The importance truncation was first developed and applied in quantum chemistry in the 1970s and latter successfully applied to models of light and intermediate mass nuclei. Other numerical methods for calculations for ultra-cold fermionic many-body systems are the Fixed-Node Diffusion Monte Carlo method (FN-DMC) and the stochastic variational approach with Correlated Gaussian basis functions (CG). There are also such method as the Coupled-Cluster method, Green’s Function Monte Carlo (GFMC) method, et cetera, used for calculation of many-body systems. In this thesis, we adopt the IT-NCSM for the calculation of ultra-cold Fermi gases at unitarity. Ultracold gases are dilute, strongly correlated systems, in which the average interparticle distance is much larger than the range of the interaction. Therefore, the detailed radial dependence of the potential is not resolved, and the potential can be replaced by an effective contact interaction. At low energy, s-wave scattering dominates and the interaction can be described by the s-wave scattering length. If the scattering length is small and negative, Cooper-pairs are formed in the Bardeen-Cooper-Schrieffer (BCS) regime. If the scattering length is small and positive, these Cooper-pairs become strongly bound molecules in a Bose-Einstein-Condensate (BEC). In between (for large scattering lengths) is the unitary limit with universal properties. Calculations of the energy spectra (ground-state and first excited-state) have so far only been performed for up to five particles using CI or NCSM methods or for up to six particles using the CG method. Calculations with larger particle numbers have only been performed with Monte Carlo methods and only for the ground state of up to 30 particles. We extend ab initio calculations of the energy spectra of ultra-cold Fermi gases at unitarity for up to 20 particle using the IT-NCSM. For our calculations we use different interactions: an effective interaction introduced by Alhassid, Bertsch and Fang and an interaction constructed using an effective field theory (EFT) approach. Furthermore, we use a Gauss-shaped potential as it is also used for the calculations of ultra-cold Fermi gases at unitarity using CG or FN-DMC methods. Although more effort must be invested to make the Gauss-shaped potential suitable for IT-NCSM calculations of ultra-cold Fermi gases, calculations with the other two interactions yield ground-state energies which agree excellently with the results obtained using the FN-DMC method. The IT-NCSM extends the range of the NCSM and permits the calculation of energy spectra (ground-state and excited-state energies) of fermionic systems with large particle numbers for which previously only the ground states could be calculated using Monte Carlo methods.

•Jan 1, 2017•0 citations
Representing a strongly interacting multiparticle wave function in a finite product basis leads to errors. Simple rescaling of the contact interaction can preserve the low-lying energy spectrum and long-wavelength structure of wave functions in one-dimensional systems and thus correct for the basis set truncation error. The analytic form of the rescaling is found for a two-particle system where the rescaling is exact. A detailed comparison between finite Hilbert space calculations and exact results for up to five particles show that rescaling can significantly improve the accuracy of numerical calculations in various external potentials. In addition to ground-state energies, the low-lying excitation spectrum, density profile, and correlation functions are studied. The results give a promising outlook for numerical simulations of trapped ultracold atoms.
When seeking a numerical representation of a quantum-mechanical multiparticle problem it is tempting to replace a singular short-range interaction by a smooth finite-range pseudopotential. Finite basis set expansions, e.g., in Fock space, are then guaranteed to converge exponentially. The need to faithfully represent the artificial length scale of the pseudopotential, however, places a costly burden on the basis set. Here we discuss scaling relations for the required size of the basis set and demonstrate the basis set convergence on the example of a two-dimensional system of few fermions with short-range sss-wave interactions in a harmonic trapping potential. In particular we show that the number of harmonic oscillator basis functions needed to reach a regime of exponential convergence for a Gaussian pseudopotential scales with the fourth power of the pseudopotential length scale, which can be improved to quadratic scaling when the basis functions are rescaled appropriately. Numerical examples for three fermions with up to a few hundred single-particle basis functions are presented and implications for the feasibility of accurate numerical multiparticle simulations of interacting ultracold-atom systems are discussed.
Abstract The configuration interaction (CI) method for calculating the exact eigenstates of a quantum-mechanical few-body system is problematic when applied to particles interacting through contact forces. In dimensions higher than one the approach fails due to the pathology of the Dirac δ -potential, making it impossible to reach convergence by gradually increasing the size of the Hilbert space. However, this problem may be cured in a rather simple manner by renormalizing the strength of the contact potential when diagonalizing in a truncated Hilbert space. One hereby relies on the comparison of the CI results to the two-body ground-state energy obtained by the exact solution of the Schrödinger equation for a regularized contact interaction. We discuss here a scheme that provides cutoff-independent few-body physical observables. The method is applied to a few-body system of ultracold atoms confined by a two-dimensional harmonic oscillator.
The configuration interaction (CI) method for calculating the exact eigenstates of a quantum-mechanical few-body system is problematic when applied to contact interactions between the particles. In two and three dimensions, the approach fails due to the pathology of the Dirac delta-potential, making it impossible to reach convergence by gradually increasing the size of the Hilbert space. However, for practical applications this problem may be cured in a rather simple manner, by renormalizing the STRENGTH of the contact potential, which must be diagonalized in a TRUNCATED Hilbert space. The procedure relies on the comparison of CI energies and wave functions with those obtained by the exact solution of the two-body Schrodinger equation for the regularized contact interaction. The rather simple scheme, while keeping the numerical procedures still elementary, nevertheless provides both cut-off-independent few-body physical observables and an estimate of the error of the CI calculation.
A new efficient numerical algorithm for interacting fermion systems is proposed and examined in detail. The ground state is expressed approximately by a linear combination of numerically chosen basis states in a truncated Hilbert space. Two procedures lead to a better approximation. The first is a numerical renormalization, which optimizes the chosen basis and projects onto the ground state within the fixed dimension, LLL, of the Hilbert space. The second is an increase of the dimension of the truncated Hilbert space, which enables the linear combination to converge to a better approximation. The extrapolation L→∞L\to\inftyL→∞ after the convergence removes the approximation error systematically. This algorithm does not suffer from the negative sign problem and can be applied to systems in any spatial dimension and arbitrary lattice structure. The efficiency is tested and the implementation explained for two-dimensional Hubbard models where Slater determinants are employed as chosen basis. Our results with less than 400 chosen basis indicate good accuracy within the errorbar of the best available results as those of the quantum Monte Carlo for energy and other physical quantities.
We explore a new variational principle for studying one-dimensional quantum systems in a trapping potential. We focus on the Fermi polaron problem, where a single distinguishable impurity interacts through a contact potential with a background of identical fermions. We can accurately describe this system at arbitrary finite repulsion by constructing a truncated basis containing states at both the limits of zero and infinite repulsion. We show how to construct this basis and how to obtain energies, density matrices and correlation functions, and provide results both for a harmonic well and a double well for various particle numbers. The results are compared both with matrix product states methods and with the analytical result for two particles in a harmonic well.
The exact solution of quantum mechanical many-body problems is only possible for few particles. Therefore, numerical methods were developed in the fields of quantum physics and quantum chemistry for larger particle numbers. Configuration Interaction (CI) methods or the No-Core Shell Model (NCSM) allow ab initio calculations for light and intermediate-mass nuclei, without resorting to phenomenology. An extension of the NCSM is the Importance-Truncated No-Core Shell Model, which uses an a priori selection of the most important basis states. The importance truncation was first developed and applied in quantum chemistry in the 1970s and latter successfully applied to models of light and intermediate mass nuclei. Other numerical methods for calculations for ultra-cold fermionic many-body systems are the Fixed-Node Diffusion Monte Carlo method (FN-DMC) and the stochastic variational approach with Correlated Gaussian basis functions (CG). There are also such method as the Coupled-Cluster method, Green’s Function Monte Carlo (GFMC) method, et cetera, used for calculation of many-body systems. In this thesis, we adopt the IT-NCSM for the calculation of ultra-cold Fermi gases at unitarity. Ultracold gases are dilute, strongly correlated systems, in which the average interparticle distance is much larger than the range of the interaction. Therefore, the detailed radial dependence of the potential is not resolved, and the potential can be replaced by an effective contact interaction. At low energy, s-wave scattering dominates and the interaction can be described by the s-wave scattering length. If the scattering length is small and negative, Cooper-pairs are formed in the Bardeen-Cooper-Schrieffer (BCS) regime. If the scattering length is small and positive, these Cooper-pairs become strongly bound molecules in a Bose-Einstein-Condensate (BEC). In between (for large scattering lengths) is the unitary limit with universal properties. Calculations of the energy spectra (ground-state and first excited-state) have so far only been performed for up to five particles using CI or NCSM methods or for up to six particles using the CG method. Calculations with larger particle numbers have only been performed with Monte Carlo methods and only for the ground state of up to 30 particles. We extend ab initio calculations of the energy spectra of ultra-cold Fermi gases at unitarity for up to 20 particle using the IT-NCSM. For our calculations we use different interactions: an effective interaction introduced by Alhassid, Bertsch and Fang and an interaction constructed using an effective field theory (EFT) approach. Furthermore, we use a Gauss-shaped potential as it is also used for the calculations of ultra-cold Fermi gases at unitarity using CG or FN-DMC methods. Although more effort must be invested to make the Gauss-shaped potential suitable for IT-NCSM calculations of ultra-cold Fermi gases, calculations with the other two interactions yield ground-state energies which agree excellently with the results obtained using the FN-DMC method. The IT-NCSM extends the range of the NCSM and permits the calculation of energy spectra (ground-state and excited-state energies) of fermionic systems with large particle numbers for which previously only the ground states could be calculated using Monte Carlo methods.