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/9707107 · Physical Review A

Exact diagonalization of the Hamiltonian for trapped interacting bosons in lower dimensions

Tor Haugset, Hårek Haugerud

Published May 1, 1998• 60 citations
arXivHTMLPDF

Abstract

We consider systems of a small number of interacting bosons confined to harmonic potentials in one and two dimensions. By exact numerical diagonalization of the many-body Hamiltonian we determine the low-lying excitation energies and the ground-state energy and density profile. We discuss the dependence of these quantities on both interaction strength $g$ and particle number $N$. The ground-state properties are compared to the predictions of the Gross-Pitaevskii equation, and the agreement is surprisingly good even for relatively low particle numbers. We also calculate the specific heat based on the obtained energy spectra.

DOI
10.1103/physreva.57.3809

Related papers

Ground state properties of trapped boson system with finite-range Gaussian repulsion: exact diagonalization study

Mohd. Imran, M. A. H. Ahsan

Abstract We use exact diagonalization to study an interacting system of N spinless bosons with finite-range Gaussian repulsion, confined in a quasi-two-dimensional harmonic trap with and without an introduced rotation. The diagonalization of the Hamiltonian matrix using Davidson algorithm in subspaces of quantized total angular momentum L z is carried out to obtain the N -body lowest eigenenergy and eigenstate. To bring out the effect of quantum (Bose) statistics and consequent phase stiffness (rigidity) of the variationally obtained many-body wavefunction on various physical quantities, our study spans from few-body ( N = 2) to many-body ( N = 16) systems. Further, to examine the finite-range effect of the repulsive Gaussian potential on many-body ground state properties of the Bose-condensate, we obtain the lowest eigenstate, the critical angular velocity of single vortex state and the quantum correlation (measured) in terms of von Neumann entanglement entropy and degree of condensation. It is found that for small values of the range (measured by the parameter σ ) of Gaussian potential, the ground state energy increases for few-boson (2 ⩽ N ⩽ 8) systems but decreases for many-boson ( N > 8) systems. On the other hand for relatively large values of the range of Gaussian potential, the ground state energy exhibits a monotonic decrease, regardless of the number of bosons N . For a given N , there is found an optimal value of the range of Gaussian potential for which the first vortex (with L z = N ) nucleates at a lower value of the rotational angular velocity Ω c1 compared to the zero-range ( δ -function) potential. Further, we observe that the inter-particle interaction and the introduced rotation are competing effects with latter being dominant over the former. With increase in the range of Gaussian potential, the value of von Neumann entropy decreases and the degree of condensation increases implying an enhanced quantum correlation and phase rigidity.

arXiv:1602.03811•Journal of Physics B Atomic Molecular and Optical Physics•Mar 13, 2020•3 citations

Simple renormalization-group approximation of the ground-state properties of interacting bosonic systems

Roman Werpachowski, Jerzy Kijowski

We present a simple renormalization-group method of investigating ground-state properties of interacting bosonic systems. Our method reduces the number of particles in a system, which makes numerical calculations possible for large systems. It is conceptually simple and easy to implement, and allows investigation of properties unavailable through mean-field approximations, such as one- and two-particle reduced density matrices of the ground state. As an example, we model a weakly interacting one-dimensional Bose gas in a harmonic trap. Compared to the mean-field Gross-Pitaevskii approximation, our method provides a more accurate description of the ground-state one-particle density matrix. We have also obtained the Hall-Post lower bounds for the ground-state energy of the gas. All results have been obtained by straightforward numerical diagonalization of the Hamiltonian matrix.

arXiv:cond-mat/0510436•Physical Review A•Jun 21, 2007•0 citations

Finite number of trapped bosons interacting through the harmonic Calogero interaction in one dimension

Barnali Chakrabarti, Tapan Kumar Das

We consider a finite number of trapped bosons in one dimension (1D) interacting through the harmonic Calogero interaction. Using an exact analytic solution for the energy eigenvalues we calculate the condensate fraction, ground-state energy, and specific heat of the system. For weak interaction the behavior of thermodynamic properties is close to that of noninteracting trapped bosons.

Physical Review A•Jan 27, 2010•4 citations

Exact Diagonalization Study of Bose-Condensed Gas with Finite-Range Gaussian Interaction

Mohd. Imran, M. A. H. Ahsan

We investigate a system of $N$ spinless bosons confined in quasi-two-dimensional harmonic trap with repulsive two-body finite-range Gaussian interaction potential of large $s$-wave scattering length. Exact diagonalization of the Hamiltonian matrix is carried out to obtain the $N$-body ground state as well as low-lying excited states, using Davidson algorithm in beyond lowest-Landau-level approximation. We examine the finite-range effects of the interaction potential on the many-body ground state energy as also the degree of condensation of the Bose-condensed gas. The results obtained indicate that the finite-range Gaussian interaction potential enhances the degree of condensation compared to the zero-range interaction potential. We further analyze the effect of finite-range interaction potential on the breathing mode collective excitation. Our theoretical results may be relevant for experiments currently conducted on quasi-two-dimensional Bose gas with more realistic interaction potential.

arXiv:1511.03165•Advanced Science Letters•Sep 1, 2015•10 citations

Reduced density matrix approach to one-dimensional ultracold bosonic systems

Mitchell J. Knight, Harry M. Quiney, Andy M. Martin

The variational determination of the two-boson reduced density matrix is described for a one-dimensional system of N N harmonically trapped bosons interacting via contact interaction, where N N ranges from 2 2 to 10^4 10 4 . The ground-state energies are calculated, and compared to existing methods, including the analytic case (for N=2) N = 2 ) and mean-field approaches such as the one-dimensional Gross-Pitaevskii equation and its variations. Structural properties such as the density and correlation functions are also derived, including the behaviour of the correlation function when boson coordinates coincide. This collectively demonstrates the capacity of the reduced density matrix method to accurately calculate collective ground-state properties of bosonic systems comprising few to many bosons, including the cross-over region between these extremes, across a large range of interaction strengths.

arXiv:2503.15811•SciPost Physics•Sep 28, 2026•0 citations

Interacting ultracold few-boson systems

Pere Mujal

In this thesis, we study the physical properties of several ultracold few-boson systems depending on the interactions between their constituents. Nowadays, experimentally, it is possible to have great control with high precision over the geometry and the interactions between the particles, making them an excellent setup to test directly the principles of quantum mechanics. A very interesting point is to study the evolution of their properties with the number of particles. The theoretical study of these systems pretends to microscopically understand the current experimental results and give support to new experimental developments.\n\n\t\t\t\t The method that will be used is the exact diagonalization of the Hamiltonian of the system. As we will see, in spite of the attempts to improve it, the method is limited by the fact that, in practice, it is only useful to study few-particle systems. The method has several advantages. First of all, one has access to both the ground and the excited states. In second place, the method is variational and converges to the exact solution as long as the Hilbert space in which we diagonalize is enlarged.\n\t\t\t\t Moreover, since we have access to the states of the system, it is possible to calculate any observable quantity of interest.\n\n\t\t\t\t First, we will study a system of spinless bosons trapped in a two-dimensional harmonic potential. The effect of the trap is to keep the system bound. It will be seen how the presence of a repulsive interaction changes the energy spectrum and other properties of the system. For instance, the density profile, which is usually measurable, and also the two-body distribution function, which is intimately related to the existence of correlations.\n\n\t\t\t\t Afterwards, the focus will be on the particular case of having only two bosons in the system interacting through a strong repulsive force. Inspired by the one-dimensional case where the fermionization phenomenon takes place in the strongly-interacting limit, we will study whether in two dimensions there is a resembling reminiscent effect. In other words, we will analyze if there are properties of the two strongly-interacting bosons in two dimensions that are like the ones of two noninteracting fermions.\n\n\t\t\t\t After that, we will tackle the localization phenomenon in a one-dimensional system that is caused by an external speckle potential that introduces disorder in the system. We will show that the localization is a robust phenomenon against repulsive contact interactions.\n\n\t\t\t\t Finally, we will study the influence of the spin-orbit coupling in a system of bosons with two possible pseudospin components, associated, for instance, to two hyperfine levels, confined in a two-dimensional harmonic trap. We will present an exhaustive analysis of the combined effects of the interaction and the spin-orbit coupling in the spectrum and the properties of the system. In particular we show the existence of a crossover in the ground state of the system susceptible to be experimentally identified.

•Nov 22, 2019•6 citations
We consider systems of a small number of interacting bosons confined to harmonic potentials in one and two dimensions. By exact numerical diagonalization of the many-body Hamiltonian we determine the low-lying excitation energies and the ground-state energy and density profile. We discuss the dependence of these quantities on both interaction strength ggg and particle number NNN. The ground-state properties are compared to the predictions of the Gross-Pitaevskii equation, and the agreement is surprisingly good even for relatively low particle numbers. We also calculate the specific heat based on the obtained energy spectra.
Abstract We use exact diagonalization to study an interacting system of N spinless bosons with finite-range Gaussian repulsion, confined in a quasi-two-dimensional harmonic trap with and without an introduced rotation. The diagonalization of the Hamiltonian matrix using Davidson algorithm in subspaces of quantized total angular momentum L z is carried out to obtain the N -body lowest eigenenergy and eigenstate. To bring out the effect of quantum (Bose) statistics and consequent phase stiffness (rigidity) of the variationally obtained many-body wavefunction on various physical quantities, our study spans from few-body ( N = 2) to many-body ( N = 16) systems. Further, to examine the finite-range effect of the repulsive Gaussian potential on many-body ground state properties of the Bose-condensate, we obtain the lowest eigenstate, the critical angular velocity of single vortex state and the quantum correlation (measured) in terms of von Neumann entanglement entropy and degree of condensation. It is found that for small values of the range (measured by the parameter σ ) of Gaussian potential, the ground state energy increases for few-boson (2 ⩽ N ⩽ 8) systems but decreases for many-boson ( N > 8) systems. On the other hand for relatively large values of the range of Gaussian potential, the ground state energy exhibits a monotonic decrease, regardless of the number of bosons N . For a given N , there is found an optimal value of the range of Gaussian potential for which the first vortex (with L z = N ) nucleates at a lower value of the rotational angular velocity Ω c1 compared to the zero-range ( δ -function) potential. Further, we observe that the inter-particle interaction and the introduced rotation are competing effects with latter being dominant over the former. With increase in the range of Gaussian potential, the value of von Neumann entropy decreases and the degree of condensation increases implying an enhanced quantum correlation and phase rigidity.
We present a simple renormalization-group method of investigating ground-state properties of interacting bosonic systems. Our method reduces the number of particles in a system, which makes numerical calculations possible for large systems. It is conceptually simple and easy to implement, and allows investigation of properties unavailable through mean-field approximations, such as one- and two-particle reduced density matrices of the ground state. As an example, we model a weakly interacting one-dimensional Bose gas in a harmonic trap. Compared to the mean-field Gross-Pitaevskii approximation, our method provides a more accurate description of the ground-state one-particle density matrix. We have also obtained the Hall-Post lower bounds for the ground-state energy of the gas. All results have been obtained by straightforward numerical diagonalization of the Hamiltonian matrix.
We investigate a system of NNN spinless bosons confined in quasi-two-dimensional harmonic trap with repulsive two-body finite-range Gaussian interaction potential of large sss-wave scattering length. Exact diagonalization of the Hamiltonian matrix is carried out to obtain the NNN-body ground state as well as low-lying excited states, using Davidson algorithm in beyond lowest-Landau-level approximation. We examine the finite-range effects of the interaction potential on the many-body ground state energy as also the degree of condensation of the Bose-condensed gas. The results obtained indicate that the finite-range Gaussian interaction potential enhances the degree of condensation compared to the zero-range interaction potential. We further analyze the effect of finite-range interaction potential on the breathing mode collective excitation. Our theoretical results may be relevant for experiments currently conducted on quasi-two-dimensional Bose gas with more realistic interaction potential.
The variational determination of the two-boson reduced density matrix is described for a one-dimensional system of N N harmonically trapped bosons interacting via contact interaction, where N N ranges from 2 2 to 10^4 10 4 . The ground-state energies are calculated, and compared to existing methods, including the analytic case (for N=2) N = 2 ) and mean-field approaches such as the one-dimensional Gross-Pitaevskii equation and its variations. Structural properties such as the density and correlation functions are also derived, including the behaviour of the correlation function when boson coordinates coincide. This collectively demonstrates the capacity of the reduced density matrix method to accurately calculate collective ground-state properties of bosonic systems comprising few to many bosons, including the cross-over region between these extremes, across a large range of interaction strengths.
In this thesis, we study the physical properties of several ultracold few-boson systems depending on the interactions between their constituents. Nowadays, experimentally, it is possible to have great control with high precision over the geometry and the interactions between the particles, making them an excellent setup to test directly the principles of quantum mechanics. A very interesting point is to study the evolution of their properties with the number of particles. The theoretical study of these systems pretends to microscopically understand the current experimental results and give support to new experimental developments.\n\n\t\t\t\t The method that will be used is the exact diagonalization of the Hamiltonian of the system. As we will see, in spite of the attempts to improve it, the method is limited by the fact that, in practice, it is only useful to study few-particle systems. The method has several advantages. First of all, one has access to both the ground and the excited states. In second place, the method is variational and converges to the exact solution as long as the Hilbert space in which we diagonalize is enlarged.\n\t\t\t\t Moreover, since we have access to the states of the system, it is possible to calculate any observable quantity of interest.\n\n\t\t\t\t First, we will study a system of spinless bosons trapped in a two-dimensional harmonic potential. The effect of the trap is to keep the system bound. It will be seen how the presence of a repulsive interaction changes the energy spectrum and other properties of the system. For instance, the density profile, which is usually measurable, and also the two-body distribution function, which is intimately related to the existence of correlations.\n\n\t\t\t\t Afterwards, the focus will be on the particular case of having only two bosons in the system interacting through a strong repulsive force. Inspired by the one-dimensional case where the fermionization phenomenon takes place in the strongly-interacting limit, we will study whether in two dimensions there is a resembling reminiscent effect. In other words, we will analyze if there are properties of the two strongly-interacting bosons in two dimensions that are like the ones of two noninteracting fermions.\n\n\t\t\t\t After that, we will tackle the localization phenomenon in a one-dimensional system that is caused by an external speckle potential that introduces disorder in the system. We will show that the localization is a robust phenomenon against repulsive contact interactions.\n\n\t\t\t\t Finally, we will study the influence of the spin-orbit coupling in a system of bosons with two possible pseudospin components, associated, for instance, to two hyperfine levels, confined in a two-dimensional harmonic trap. We will present an exhaustive analysis of the combined effects of the interaction and the spin-orbit coupling in the spectrum and the properties of the system. In particular we show the existence of a crossover in the ground state of the system susceptible to be experimentally identified.