Numerical Realization of Bethe Rapidities in cold quenched systems by Feynman-Kac path integral method

TL;DR

This paper employs Feynman-Kac path integral Monte Carlo to numerically analyze the non-equilibrium dynamics of a 1D Bose gas after a quench, providing an efficient alternative to Bethe Ansatz for solving the Lieb-Liniger model.

Contribution

It introduces a Monte Carlo method to compute Bethe rapidities in quenched 1D Bose gases, demonstrating improved efficiency over traditional techniques.

Findings

Monte Carlo efficiently computes density distributions.

Density remains self-similar after changing box length.

Method captures non-equilibrium dynamics effectively.

Abstract

We apply Quantum Monte Carlo technique to analyze the non equlibrium state of a trapped 1d Bose gas just after the quenching of the confining potential. As a matter of fact we solve the time dependent Schroedinger equation for the system of one-dimensional bosons interacting via delta potential in an infinite square well (namely Lieb-Liniger model) using Feynman-Kac path integral Monte Carlo technique. These 1d systems are extremely interesting and worth investigating in the context of non-equilibrium dynamics of interacting many body systems. Even though the systems can be realized experimentally and are exactly solvable by Bethe Ansatz, the diffusion Monte Carlo is proven to be more efficient in most circumstances than other mean value techniques as the numerical method can incorporate the finite interaction very easily. Using N particle ground state wavefunction for one-dimensional hard core bosons in a harmonic trap, we develop an algorithm to calculate density. We also observe the change in the density distribution by changing the length of the hard wall box. After an increase in the box length, we still get a self-similar density distribution.