Publications

Order by: publication date number of citations

• Triangular Gross-Pitaevskii breathers and Damski-Chandrasekhar shock waves

  M. Olshanii, D. Deshommes, J. Torrents, M. Gonchenko, V. Dunjko, G. E. Astrakharchik

  SciPost Phys. 10, 114 (2021) · published 26 May 2021 | Toggle abstract · pdf

  The recently proposed map [arXiv:2011.01415] between the hydrodynamic equations governing the two-dimensional triangular cold-bosonic breathers [Phys. Rev. X 9, 021035 (2019)] and the high-density zero-temperature triangular free-fermionic clouds, both trapped harmonically, perfectly explains the former phenomenon but leaves uninterpreted the nature of the initial ($t=0$) singularity. This singularity is a density discontinuity that leads, in the bosonic case, to an infinite force at the cloud edge. The map itself becomes invalid at times $t<0$. A similar singularity appears at $t = T/4$, where $T$ is the period of the harmonic trap, with the Fermi-Bose map becoming invalid at $t > T/4$. Here, we first map---using the scale invariance of the problem---the trapped motion to an untrapped one. Then we show that in the new representation, the solution [arXiv:2011.01415] becomes, along a ray in the direction normal to one of the three edges of the initial cloud, a freely propagating one-dimensional shock wave of a class proposed by Damski in [Phys.~Rev.~A 69, 043610 (2004)]. There, for a broad class of initial conditions, the one-dimensional hydrodynamic equations can be mapped to the inviscid Burgers' equation, which is equivalent to a nonlinear transport equation. More specifically, under the Damski map, the $t=0$ singularity of the original problem becomes, verbatim, the initial condition for the wave catastrophe solution found by Chandrasekhar in 1943 [Ballistic Research Laboratory Report No. 423 (1943)]. At $t=T/8$, our interpretation ceases to exist: at this instance, all three effectively one-dimensional shock waves emanating from each of the three sides of the initial triangle collide at the origin, and the 2D-1D correspondence between the solution of [arXiv:2011.01415] and the Damski-Chandrasekhar shock wave becomes invalid.

  6 citations

• Collective oscillations of a two-component Fermi gas on the repulsive branch

  Tomasz Karpiuk, Piotr T. Grochowski, Mirosław Brewczyk, Kazimierz Rzążewski

  SciPost Phys. 8, 066 (2020) · published 21 April 2020 | Toggle abstract · pdf

  We calculate frequencies of collective oscillations of two-component Fermi gas that is kept on the repulsive branch of its energy spectrum. Not only is a paramagnetic phase explored, but also a ferromagnetically separated one. Both in-, and out-of-phase perturbations are investigated, showing contributions from various gas excitations. Additionally, we compare results coming from both time-dependent Hartree-Fock and density-functional approaches.

  6 citations

• Invariant approach to the driven Jaynes-Cummings model

  Ivan Alejandro Bocanegra-Garay, L. Hernández-Sánchez, Irán Ramos-Prieto, Francisco Soto-Eguibar, Héctor Manuel Moya-Cessa

  SciPost Phys. 16, 007 (2024) · published 8 January 2024 | Toggle abstract · pdf

  We investigate the dynamics of the driven Jaynes-Cummings model, where a two-level atom interacts with a quantized field and both, atom and field, are driven by an external classical field. Via an invariant approach, we are able to transform the corresponding Hamiltonian into the one of the standard Jaynes-Cummings model. Subsequently, the exact analytical solution of the Schrödinger equation for the driven system is obtained and employed to analyze some of its dynamical variables.

  6 citations

• On computing non-equilibrium dynamics following a quench

  Neil J. Robinson, Albertus J. J. M. de Klerk, Jean-Sébastien Caux

  SciPost Phys. 11, 104 (2021) · published 9 December 2021 | Toggle abstract · pdf

  Computing the non-equilibrium dynamics that follows a quantum quench is difficult, even in exactly solvable models. Results are often predicated on the ability to compute overlaps between the initial state and eigenstates of the Hamiltonian that governs time evolution. Except for a handful of known cases, it is generically not possible to find these overlaps analytically. Here we develop a numerical approach to preferentially generate the states with high overlaps for a quantum quench starting from the ground state or an excited state of an initial Hamiltonian. We use these preferentially generated states, in combination with a "high overlap states truncation scheme" and a modification of the numerical renormalization group, to compute non-equilibrium dynamics following a quench in the Lieb-Liniger model. The method is non-perturbative, works for reasonable numbers of particles, and applies to both continuum and lattice systems. It can also be easily extended to more complicated scenarios, including those with integrability breaking.

  6 citations

• Phases and dynamics of quantum droplets in the crossover to two-dimensions

  Jose Carlos Pelayo, George Bougas, Thomás Fogarty, Thomas Busch, Simeon I. Mistakidis

  SciPost Phys. 18, 129 (2025) · published 17 April 2025 | Toggle abstract · pdf

  We explore the ground states and dynamics of ultracold atomic droplets in the crossover region from three to two dimensions by solving the two-dimensional and the quasi two-dimensional extended Gross-Pitaevskii equations numerically and with a variational approach. By systematically comparing the droplet properties, we determine the validity regions of the pure two-dimensional description, and therefore the dominance of the logarithmic nonlinear coupling, as a function of the sign of the averaged mean-field interactions and the size of the transverse confinement. One of our main findings is that droplets become substantially extended upon transitioning from negative-to-positive averaged mean-field interactions. This is accompanied by a significant reduction of their binding energies which are approximately inversely proportional to the square of their size. To explore fundamental dynamical properties in the crossover region, we study interaction quenches and show that the droplets perform a periodic breathing motion for modest quench strengths, while larger quench amplitudes lead to continuous expansion exhibiting density ring structures. We also showcase that it is possible to form complex bulk and surface density patterns in anisotropic geometries following the quench. Since we are working with realistic parameters, our results can directly facilitate future experimental realizations.

  5 citations

• Many-body quantum dynamics with MCTDH-X

  Paolo Molignini, Sunayana Dutta, Elke Fasshauer

  SciPost Phys. Lect. Notes 94 (2025) · published 26 March 2025 | Toggle abstract · pdf

  The lecture notes on "Many-body Quantum Dynamics with MCTDH-X", adapted from the 2023 Heidelberg MCTDH Summer School, provide an in-depth exploration of the Multiconfigurational Time-Dependent Hartree approach for indistinguishable particles. They serve as a comprehensive guide for understanding and utilizing the MCTDH-X software for both bosonic and fermionic systems. The tutorial begins with an introduction to the MCTDH-X software, highlighting its capability to handle various quantum systems, including those with internal degrees of freedom and long-range interactions. The theoretical foundation is then laid out on how to solve the time-dependent and time-independent Schrödinger equations for many-body systems. The workflow section provides practical instructions on setting up and executing simulations using MCTDH-X. Detailed benchmarks against exact solutions are presented, showcasing the accuracy and reliability of the software in ground-state and dynamic simulations. The notes then delve into the dynamics of quantum systems, covering relaxation processes, time evolution, and the analysis of propagation for both bosonic and fermionic particles. The discussion includes the interpretation of various physical quantities such as energy, density distributions, and orbital occupations. Advanced features of MCTDH-X are also explored in the last section, including the calculation of correlation functions and the creation of visualizations through video tutorials. The notes conclude with a Linux/UNIX command cheat sheet, facilitating ease of use for users operating the software on different systems. Overall, these lecture notes provide a valuable resource for researchers and students in the field of quantum dynamics, offering both theoretical insights and practical guidance on the use of MCTDH-X for studying complex many-body systems.

  5 citations

• Multimode-polariton superradiance via Floquet engineering

  Christian Høj Johansen, Johannes Lang, Andrea Morales, Alexander Baumgärtner, Tobias Donner, Francesco Piazza

  SciPost Phys. 12, 094 (2022) · published 16 March 2022 | Toggle abstract · pdf

  We consider an ensemble of ultracold bosonic atoms within a near-planar cavity, driven by a far detuned laser whose phase is modulated at a frequency comparable to the transverse cavity mode spacing. We show that a strong, dispersive atom-photon coupling can be reached for many transverse cavity modes at once. The resulting Floquet polaritons involve a superposition of a set of cavity modes with a density excitation of the atomic cloud. The mutual interactions between these modes lead to distinct avoided crossings between the polaritons. Increasing the laser drive intensity, a low-lying multimode Floquet polariton softens and eventually becomes undamped, corresponding to the transition to a superradiant, self-organized phase. We demonstrate the stability of the stationary state for a broad range of parameters.

  5 citations

• Analysis of the buildup of spatiotemporal correlations and their bounds outside of the light cone

  Nils O. Abeling, Lorenzo Cevolani, Stefan Kehrein

  SciPost Phys. 5, 052 (2018) · published 26 November 2018 | Toggle abstract · pdf

  In non-relativistic quantum theories the Lieb-Robinson bound defines an effective light cone with exponentially small tails outside of it. In this work we use it to derive a bound for the time evolution of the correlation function of two local disjoint observables if the initial state has a power-law decay. We show that the exponent of the power-law of the bound is identical to the initial (equilibrium) decay. We explicitly verify this result by studying the full dynamics of the susceptibilities and correlations in the exactly solvable Luttinger model after a sudden quench from the non-interacting to the interacting model.

  5 citations

• Time-dependent variational Monte Carlo study of the dynamic response of bosons in an optical lattice

  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.

  5 citations

• Quantum phases of hardcore bosons with repulsive dipolar density-density interactions on two-dimensional lattices

  J. A. Koziol, G. Morigi, K. P. Schmidt

  SciPost Phys. 17, 111 (2024) · published 10 October 2024 | Toggle abstract · pdf

  We analyse the ground-state quantum phase diagram of hardcore Bosons interacting with repulsive dipolar potentials. The bosons dynamics is described by the extended-Bose-Hubbard Hamiltonian on a two-dimensional lattice. The ground state results from the interplay between the lattice geometry and the long-range interactions, which we account for by means of a classical spin mean-field approach limited by the size of the considered unit cells. This extended classical spin mean-field theory accounts for the long-range density-density interaction without truncation. We consider three different lattice geometries: square, honeycomb, and triangular. In the limit of zero hopping the ground state is always a devil's staircase of solid (gapped) phases. Such crystalline phases with broken translational symmetry are robust with respect to finite hopping amplitudes. At intermediate hopping amplitudes, these gapped phases melt, giving rise to various lattice supersolid phases, which can have exotic features with multiple sublattice densities. At sufficiently large hoppings the ground state is a superfluid. The stability of phases predicted by our approach is gauged by comparison to the known quantum phase diagrams of the Bose-Hubbard model with nearest-neighbour interactions as well as quantum Monte Carlo simulations for the dipolar case on the square and triangular lattice. Our results are of immediate relevance for experimental realisations of self-organised crystalline ordering patterns in analogue quantum simulators, e.g., with ultracold dipolar atoms in an optical lattice.

  5 citations