Contributor info: Dr Robert Zillich

Mathias Gartner, David Miesbauer, Michael Kobler, Julia Freund, Giuseppe Carleo, Robert E. Zillich

SciPost Phys. 18, 123 (2025) · published 9 April 2025 | Toggle abstract · pdf

We present a new variational method to study the dynamics of a closed bosonic many-body system, the time-dependent hypernetted-chain Euler-Lagrange method, tHNC. Based on the Jastrow ansatz, it accounts for quantum fluctuations in a non-perturbative way. The tHNC method scales well with the number of dimensions, as demonstrated by our results on one-, two-, and three-dimensional systems. We apply the tHNC method to interaction quenches, i.e. sudden changes of the interaction strength, in homogeneous Bose gases. When the quench is strong enough that the final state has roton excitations (as found and predicted for dipolar and Rydberg-dressed Bose-Einstein condensates, respectively), the pair distribution function exhibits stable oscillations. For validation, we compare tHNC results with time-dependent variational Monte Carlo results in one and two dimensions.

Mathias Gartner, Ferran Mazzanti, Robert E. Zillich

SciPost Phys. 13, 025 (2022) · published 18 August 2022 | Toggle abstract · pdf

We study the dynamics of a one-dimensional Bose gas at unit filling in both shallow and deep optical lattices and obtain the dynamic structure factor ${S(k,\omega)}$ by monitoring the linear response to a weak probe pulse. We introduce a new procedure, based on the time-dependent variational Monte Carlo method (tVMC), which allows to evolve the system in real time, using as a variational model a Jastrow-Feenberg wave function that includes pair correlations. Comparison with exact diagonalization results of ${S(k,\omega)}$ obtained on a lattice in the Bose-Hubbard limit shows good agreement of the dispersion relation for sufficiently deep optical lattices, while for shallow lattices we observe the influence of higher Bloch bands. We also investigate non-linear response to strong pulses. From the power spectrum of the density fluctuations we obtain the excitation spectrum, albeit broadened, by higher harmonic generation after a strong pulse with a single low wave number. As a remarkable feature of our simulations we furthermore demonstrate that the full excitation spectrum can be retrieved from the power spectrum of the density fluctuations due to the stochastic noise inherent in any Monte Carlo method, without applying an actual perturbation.