Monte Carlo Ground State Energy for Trapped Boson Systems

M. E. Grypeos

Jing Jin, Tang Yi

The diffusion Monte Carlo method is applied to study the ground-state properties of charged bosons in one dimension confined in a harmonic double-well trap. The particles interact repulsively through a Coulombic 1/ r potential. Numerical results show that the well separation has significant influence on the ground-state properties of the system. When the interaction of the system is weak, ground-state energy decreases with the increasing well separation and has a minimal value. If the well separation increases continually, the ground-state energy increases and approaches to a constant gradually. This effect will be abatable in the strong interacting system. In addition, by calculating the density of the systems for different interaction strengths with various well separations, we find that the density increases abnormally when the well separation is large at the centre of the system.

Augusto González, B. Partoens, Algirdas Matulis +1 more

The Pad\'e approximant technique and the variational Monte Carlo method are applied to determine the ground-state energy of a finite number of charged bosons in two dimensions confined by a parabolic trap. The particles interact repulsively through a Coulombic $1/r$ potential. Analytic expressions for the ground-state energy are obtained. The convergence of the Pad\'e sequence and comparison with the Monte Carlo results show that the error of the Pad\'e estimate is less than 4% at any boson density and is exact in the extreme situations of very dilute and high density.

Hans A. De Raedt

Naoki Kawashima, 川島 直輝

Zheng Rong-Jie, Jing Jin, Yi Ping Tang

The ground-state properties of a system with a small number of interacting bosons over a wide range of densities are investigated. The system is confined in a two-dimensional isotropic harmonic trap, where the interaction between bosons is treated as a hard-core potential. By using variational Monte Carlo method, we diagonalize the one-body density matrix of the system to obtain the ground-state energy, condensate wavefunction and the condensate fraction. We find that in the dilute limit the depletion of central condensate in the 2D system is larger than in a 3D system for the same interaction strength; however as the density increases, the depletion at the centre of 2D trap will be equal to or even lower than that at the centre of 3D trap, which is in agreement with the anticipated in Thomas–Fermi approximation. In addition, in the 2D system the total condensate depletion is still larger than in a 3D system for the same scattering length.