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

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.