Source author record

Axel U. J. Lode

Axel U. J. Lode appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

21works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

21 published item(s)

preprint2022arXiv

Pauli crystal melting in shaken optical traps

Pauli crystals are ordered geometric structures that emerge in trapped noninteracting fermionic systems due to their underlying Pauli repulsion. The deformation of Pauli crystals - often called melting - has been recently observed in experiments, but the mechanism that leads to it remains unclear. We address this question by studying the melting dynamics of N=6 fermions as a function of periodic driving and experimental imperfections in the trap (anisotropy and anharmonicity) by employing a combination of numerical simulations and Floquet theory. Surprisingly, we reveal that the melting of Pauli crystals is not simply a direct consequence of an increase in system energy, but is instead related to the trap geometry and the population of the Floquet modes. We show that the melting is absent in traps without imperfections and triggered only by a sufficiently large shaking amplitude in traps with imperfections.

preprint2021arXiv

Fermionization via cavity-assisted infinite-range interactions

We study a one-dimensional array of bosons with infinite-range interactions mediated by a laser-driven dissipative optical cavity. The cavity-mediated infinite-range interactions open up a new pathway to fermionization, hitherto only known for dipolar bosons due to their long-range interactions. In parameter ranges attainable in state-of-the-art experiments, we systematically compare observables for bosons and fermions with infinite-range interactions. At large enough laser pump power, many observables, including density distributions in real and momentum space, correlation functions, eigenvalues of the one-body density matrix, and superradiance order parameter, become identical for bosons and fermions. We map out the emergence of this cavity-induced fermionization as a function of pump power and contact interactions. We discover that cavity-mediated interactions can compensate a reduction by several orders of magnitude in the strength of the contact interactions needed to trigger fermionization.

preprint2021arXiv

Many-body density and coherence of trapped cold bosons

Many-body densities and correlation functions are of paramount importance for understanding quantum many-body physics. Here, we present a method to compute them; our approach is general and based on the action of bosonic or fermionic annihilation field operators on the many-body wavefunction. We analyze $N = 6$ quasi-one-dimensional harmonically-trapped bosons with weak to strong contact interaction strength up to the Tonks-Girardeau limit with infinite repulsion using the MultiConfigurational Time-Dependent Hartree method for indistinguishable particles (MCTDH-X). We compare our MCTDH-X solutions to the analytical ones in the infinite repulsion regime as well as to the so-called correlated pair wavefunction approach and find a good agreement. Since numerical approximations are not bound to the cases where analytical solutions are known, we thus demonstrate a general method to investigate high-order reduced density matrices and correlation functions in systems for which analytical solutions are unknown. We trace the build-up of correlation features in the crossover from weak interactions to the Tonks-Girardeau limit and find that the higher-order correlation functions and densities resemble those in the Tonks-Girardeau limit for way smaller interactions than anticipated from just the one-body density.

preprint2020arXiv

Detecting One-Dimensional Dipolar Bosonic Crystal Orders via Full Distribution Functions

We explore the groundstates of a few dipolar bosons in optical lattices with incommensurate filling. The competition of kinetic, potential, and interaction energies leads to the emergence of a variety of crystal state orders with characteristic one- and two-body densities. We probe the transitions between these orders and construct the emergent state diagram as a function of the dipolar interaction strength and the lattice depth. We show that the crystal state orders can be observed using the full distribution functions of the particle number extracted from simulated single-shot images.

preprint2020arXiv

MCTDH-X: The multiconfigurational time-dependent Hartree method for indistinguishable particles software

We introduce and describe the multiconfigurational time-depenent Hartree for indistinguishable particles (MCTDH-X) software. This powerful tool allows the investigation of ground state properties and dynamics of interacting quantum many-body systems in different spatial dimensions. The MCTDH-X software is a set of programs and scripts to compute, analyze, and visualize solutions for the time-dependent and time-independent many-body Schrödinger equation for indistinguishable quantum particles. As the MCTDH-X software represents a general solver for the Schrödinger equation, it is applicable to a wide range of problems in the fields of atomic, optical, molecular physics as well as condensed matter systems. In particular, it can be used to study light-matter interactions, correlated dynamics of electrons, as well as some aspects related to quantum information and computing. The MCTDH-X software solves a set of non-linear coupled working equations based on the application of the variational principle to the Schrödinger equation. These equations are obtained by using an ansatz for the many-body wavefunction that is a time-dependent expansion in a set of time-dependent many-body basis states. The time-dependence of the basis set enables MCTDH-X to deal with quantum dynamics at a superior accuracy as compared to, for instance, exact diagonalization approaches. Herein, we give an introduction to the MCTDH-X software via an easy-to-follow tutorial with a focus on accessibility. We use the double well to illustrate the fermionization of bosonic particles, the crystallization of fermionic particles, characteristics of the superfluid and Mott-insulator quantum phases in Hubbard models, and even dynamical quantum phase transitions. Our tutorial guides the potential user to apply the MCTDH-X software also to more complex systems.

preprint2020arXiv

Multiconfigurational time-dependent Hartree approaches for indistinguishable particles

In this Colloquium, the wavefunction-based Multiconfigurational Time-Dependent Hartree approaches to the dynamics of indistinguishable particles (MCTDH-F for Fermions and MCTDH-B for Bosons) are reviewed. MCTDH-B and MCTDH-F or, together, MCTDH-X are methods for describing correlated quantum systems of identical particles by solving the time-dependent Schrödinger equation from first principles. MCTDH-X is used to accurately model the dynamics of real-world quantum many-body systems in atomic, molecular, and optical physics. The key feature of these approaches is the time-dependence and optimization of the single-particle states employed for the construction of a many-body basis set, which yields nonlinear working equations. We briefly describe the historical developments that have lead to the formulation of the MCTDH-X methods and motivate the necessity for wavefunction-based approaches. We sketch the derivation of the unified MCTDH-F and MCTDH-B equations of motion for complete and also specific restricted configuration spaces. The strengths and limitations of the MCTDH-X approach are assessed via benchmarks against an exactly solvable model and via convergence checks. We highlight some applications to instructive and experimentally-realized quantum many-body systems: the dynamics of atoms in Bose-Einstein condensates in magneto-optical and optical traps and of electrons in atoms and molecules. We discuss the current development and frontiers in the field of MCTDH-X: theories and numerical methods for indistinguishable particles, for mixtures of multiple species of indistinguishable particles, the inclusion of nuclear motion for the nonadiabatic dynamics of atomic and molecular systems, as well as the multilayer and second-quantized-representation approaches, and the orbital-adaptive time-dependent coupled-cluster theory are discussed.

preprint2020arXiv

Pathway to chaos through hierarchical superfluidity in a cavity-BEC system

We explore the role of atomic correlations in a harmonically trapped Bose-Einstein condensate coupled to a dissipative cavity, where both the atoms and the cavity are blue detuned from the external pumping laser. Using a genuine many-body approach that goes beyond mean-field, we extract density distributions and many-body correlations to unveil a pathway to chaos at large pump power through a hierarchical self-organization of the atoms, where the atoms transition from a single-well optical lattice to a double-well optical lattice. Correlated states of the atoms emerge and are characterized by local superfluid correlations in phases which are globally superfluid or Mott insulating. Local superfluid-Mott transitions are precluded by a dynamical instability to chaos which occurs via quasiperiodic attractors. Our results explain the mechanism behind the dynamical instabilities observed in experiments.

preprint2020arXiv

Spectral Structure and Many-Body Dynamics of Ultracold Bosons in a Double-Well

We examine the spectral structure and many-body dynamics of two and three repulsively interacting bosons trapped in a one-dimensional double-well, for variable barrier height, inter-particle interaction strength, and initial conditions. By exact diagonalization of the many-particle Hamiltonian, we specifically explore the dynamical behaviour of the particles launched either at the single particle ground state or saddle point energy, in a time-independent potential. We complement these results by a characterisation of the cross-over from diabatic to quasi-adiabatic evolution under finite-time switching of the potential barrier, via the associated time-evolution of a single particle's von Neumann entropy. This is achieved with the help of the multiconfigurational time-dependent Hartree method for indistinguishable particles (\textsc{Mctdh-x}) -- which also allows us to extrapolate our results for increasing particle numbers.

preprint2016arXiv

Dynamics of Hubbard Hamiltonians with the multiconfigurational time-dependent Hartree method for indistinguishable particles

We apply the multiconfigurational time-dependent Hartree method for indistinguishable particles (MCTDH-X) to systems of bosons or fermions in lattices described by Hubbard type Hamiltonians with long-range or short-range interparticle interactions. The wavefunction is expanded in a variationally optimized time-dependent many-body basis generated by a set of effective creation operators that are related to the original particle creation operators by a time-dependent unitary transform. We use the time-dependent variational principle for the coefficients of this transform as well as the expansion coefficients of the wavefunction in the time-dependent many-body basis as variational parameters to derive equations of motion. The convergence of MCTDH-X is shown by comparing its results to the exact diagonalization of one-, two-, and three-dimensional lattices filled with bosons with contact interactions. We use MCTDH-X to study the buildup of correlations in the long-time splitting dynamics of a Bose-Einstein condensate loaded into a large two-dimensional lattice subject to a barrier that is ramped up in the center. We find that the system is split into two parts with emergent time-dependent correlations that depend on the ramping time -- for most barrier-raising-times the system becomes two-fold fragmented, but for some of the very fast ramps, the system shows revivals of coherence.

preprint2016arXiv

The multiconfigurational time-dependent Hartree for fermions: Implementation, exactness and application to few-fermion tunneling to open space

We report on an implementation of the multiconfigurational time-dependent Hartree method (MCTDH) for spin-polarized fermions (MCTDHF). Our approach is based on a mapping for opera- tors in Fock space that allows a compact and efficient application of the Hamiltonian and solution of the MCTDHF equations of motion. Our implementation extends, builds on and exploits the recursive implementation of MCTDH for bosons (R-MCTDHB) package. Together with R-MCTDHB, the present implementation of MCTDHF forms the MCTDH-X package. We benchmark the accuracy of the algorithm with the harmonic interaction model and a time-dependent generalization thereof. These models consider parabolically trapped particles that interact through a harmonic interaction potential. We demonstrate, that MCTDHF is capable of solving the time-dependent many-fermion Schrödinger equation to an in principle arbitrary degree of precision and can hence yield numerically exact results even in the case of Hamiltonians with time-dependent one-body and two-body potentials. As an application we study the problem of two initially parabolically confined and charged fermions tunneling through a barrier to open space. We demonstrate the validity of a model proposed previously for the many-body tunneling to open space of bosonic particles with contact interactions [Proc. Natl. Acad. Sci. USA 109, 13521-13525 (2012)]. The many-fermion tunneling can be built up from sequentially happening single-fermion tunneling processes. The characteristic momenta of these processes are determined by the chemical potentials of trapped subsystems of smaller particle numbers: the escaped fermions convert the different chemical potentials into kinetic energy. Using the two-body correlation function, we present a detailed picture of the sequentiality of the process and are able to tell tunneling from over-the-barrier escape.

preprint2016arXiv

The multiconfigurational time-dependent Hartree method for bosons with internal degrees of freedom: Theory and composite fragmentation of multi-component Bose-Einstein condensates

In this paper the multiconfigurational time-dependent Hartree for bosons method (MCTDHB) is derived for the case of $N$ identical bosons with internal degrees of freedom. The theory for bosons with internal degrees of freedom constitutes a generalization of the MCTDHB method that substantially enriches the many-body physics that can be described. We demonstrate that the numerically exact solution of the time-dependent many-body Schrödinger equation for interacting bosonic particles with internal degrees of freedom is now feasible. We report on the MCTDHB equations of motion for bosons with internal degrees of freedom and their implementation for a general many-body Hamiltonian with one-body and two-body terms that, both, may depend on the internal states of the considered particles. To demonstrate the capabilities of the theory and its software implementation integrated in the MCTDH-X software, we apply MCTDHB to the emergence of fragmentation of parabolically trapped bosons with two internal states: we study the groundstate of $N=100$ parabolically confined bosons as a function of the splitting between the state-dependent minima of the two parabolic potentials. To quantify the coherence of the system we compute its normalized one-body correlation function. We find that the coherence within each internal state of the atoms is maintained, while it is lost between the different internal states. This is a hallmark of a new kind of fragmentation which is absent in bosons without internal structure. We term the emergent phenomenon "composite fragmentation".

preprint2015arXiv

Condensate fragmentation as a sensitive measure of the quantum many-body behavior of bosons with long-range interactions

The occupation of more than one single-particle state and hence the emergence of fragmentation is a many-body phenomenon universal to systems of spatially confined interacting bosons. In the present study, we investigate the effect of the range of the interparticle interactions on the fragmentation degree of one- and two-dimensional systems. We solve the full many-body Schrödinger equation of the system using the recursive implementation of the multiconfigurational time-dependent Hartree for bosons method, R-MCTDHB. The dependence of the degree of fragmentation on dimensionality, particle number, areal or line density and interaction strength is assessed. It is found that for contact interactions, the fragmentation is essentially density independent in two dimensions. However, fragmentation increasingly depends on density the more long-ranged the interactions become. The degree of fragmentation is increasing, keeping the particle number $N$ fixed, when the density is decreasing as expected in one spatial dimension. We demonstrate that this remains, nontrivially, true also for long-range interactions in two spatial dimensions. We, finally, find that within our fully self-consistent approach, the fragmentation degree, to a good approximation, decreases universally as $N^{-1/2}$ when only $N$ is varied.

preprint2014arXiv

Breaking the resilience of a two-dimensional Bose-Einstein condensate to fragmentation

A two-dimensional Bose-Einstein condensate (BEC) split by a radial potential barrier is investigated. We determine on an accurate many-body level the system's ground-state phase diagram as well as a time-dependent phase diagram of the splitting process. Whereas the ground state is condensed for a wide range of parameters, the time-dependent splitting process leads to substantial fragmentation. We demonstrate for the first time the dynamical fragmentation of a BEC despite its ground state being condensed. The results are analyzed by a mean-field model and suggest that a large manifold of low-lying fragmented excited states can significantly impact the dynamics of trapped two-dimensional BECs.

preprint2014arXiv

Vortex nucleation through fragmentation in a stirred resonant Bose-Einstein condensate

Superfluids are distinguished from ordinary fluids by the quantized manner the rotation is manifested in them. Precisely, quantized vortices are known to appear in the bulk of a superfluid subject to external rotation. In this work we study a trapped ultracold Bose gas of $N=101$ atoms in two spatial dimensions that is stirred by a rotating beam. We use the multiconfigurational Hartree method for bosons, that extends the mainstream mean-field theory, to calculate the dynamics of the gas in real time. As the gas is rotated the wavefunction of the system changes symmetry and topology. We see a series of resonant rotations as the stirring frequency is increased. Fragmentation accompanies the resonances and change of symmetry of the wavefunction of the gas. We conclude that fragmentation of the gas appears hand-in-hand with resonant absorption of energy and angular momentum from the external agent of rotation.

preprint2013arXiv

Controlling the Velocities and Number of Emitted Particles in the Tunneling to Open Space Dynamics

A scheme to control the many-boson tunneling process to open space is derived and demonstrated. The number of ejected particles and their velocities can be controlled by two parameters, the threshold of the potential and the interparticle interaction. Since these parameters are fully under experimental control, this is also the case for the number of ejected particles and their emission spectrum. The process of tunneling to open space can hence be used, for example, for the quantum simulation of complicated tunneling ionization processes and atom lasers. To understand the many-body tunneling process, a generalization of the model introduced in [Proc. Natl. Acad. Sci. USA, 109, 13521 (2012)] for tunneling in the absence of a threshold is put forward and proven to apply for systems with a non-zero threshold value. It is demonstrated that the model is applicable for general interparticle interaction strengths, particle numbers and threshold values. The model constructs the many-body process from single-particle emission processes. The rates and emission momenta of the single-particle processes are determined by the chemical potentials and energy differences to the threshold value of the potential for systems with different particle numbers. The chemical potentials and these energy differences depend on the interparticle interaction. Both the number of confined particles and their rate of emission thus allow for a control by the manipulation of the interparticle interaction and the threshold. Numerically exact results for two, three and one hundred bosons are shown and discussed. The devised control scheme for the many-body tunneling process performs very well for the dynamics of the momentum density, the correlations, the coherence and of the final state, i.e., the number of particles that remain confined in the potential.

preprint2013arXiv

Elastic scattering of a Bose-Einstein condensate at a potential landscape

We investigate the elastic scattering of Bose-Einstein condensates at shallow periodic and disorder potentials. We show that the collective scattering of the macroscopic quantum object couples to internal degrees of freedom of the Bose-Einstein condensate such that the Bose-Einstein condensate gets depleted. As a precursor for the excitation of the Bose-Einstein condensate we observe wave chaos within a mean-field theory.

preprint2013arXiv

Excitation spectra of many-body systems by linear response: General theory and applications to trapped condensates

We derive a general linear-response many-body theory capable of computing excitation spectra of trapped interacting bosonic systems, e.g., depleted and fragmented Bose-Einstein condensates (BECs). To obtain the linear-response equations we linearize the multiconfigurational time-dependent Hartree for bosons (MCTDHB) method, which provides a self-consistent description of many-boson systems in terms of orbitals and a state vector (configurations), and is in principle numerically-exact. The derived linear-response many-body theory, which we term LR-MCTDHB, is applicable to systems with interaction potentials of general form. From the numerical implementation of the LR-MCTDHB equations and solution of the underlying eigenvalue problem, we obtain excitations beyond available theories of excitation spectra, such as the Bogoliubov-de Gennes (BdG) equations. The derived theory is first applied to study BECs in a one-dimensional harmonic potential. The LR-MCTDHB method contains the BdG excitations and, also, predicts a plethora of additional many-body excitations which are out of the realm of standard linear response. In particular, our theory describes the exact energy of the higher harmonic of the first (dipole) excitation not contained in the BdG theory. We next study a BEC in a very shallow one-dimensional double-well potential. We find with LR-MCTDHB low-lying excitations which are not accounted for by BdG, even though the BEC has only little fragmentation and, hence, the BdG theory is expected to be valid. The convergence of the LR-MCTDHB theory is assessed by systematically comparing the excitation spectra computed at several different levels of theory.

preprint2012arXiv

Exact quantum dynamics of bosons with finite-range time-dependent interactions of harmonic type

The exactly solvable quantum many-particle model with harmonic one- and two-particle interaction terms is extended to include time-dependency. We show that when the external trap potential and finite-range interparticle interaction have a time-dependency the exact solutions of the corresponding time-dependent many-boson Schrödinger equation are still available. We use these exact solutions to benchmark the recently developed multiconfigurational time-dependent Hartree method for bosons (MCTDHB) [Phys. Rev. Lett. {\bf 99}, 030402 (2007), Phys. Rev. A {\bf 77}, 033613 (2008)]. In particular, we benchmark the MCTDHB method for: (i) the ground state; (ii) the breathing many-body dynamics activated by a quench scenario where the interparticle interaction strength is suddenly turned on to a finite value; (iii) the non-equilibrium dynamic for driven scenarios where both the trap- and interparticle-interaction potentials are {\it time-dependent}. Excellent convergence of the ground state and dynamics is demonstrated. The great relevance of the self-consistency and time-adaptivity, which are the intrinsic features of the MCTDHB method, is demonstrated by contrasting the MCTDHB predictions and those obtained within the standard full configuration interaction method spanning the Fock space of the same size, but utilizing as one-particle basis set the fixed-shape eigenstates of the one-particle potential. Connections of the model's results to ultra-cold Bose-Einstein condensed systems are addressed.

preprint2012arXiv

Wave chaos as signature for depletion of a Bose-Einstein condensate

We study the expansion of repulsively interacting Bose-Einstein condensates (BECs) in shallow one-dimensional potentials. We show for these systems that the onset of wave chaos in the Gross-Pitaevskii equation (GPE), i.e. the onset of exponential separation in Hilbert space of two nearby condensate wave functions, can be used as indication for the onset of depletion of the BEC and the occupation of excited modes within a many-body description. Comparison between the multiconfigurational time-dependent Hartree for bosons (MCTDHB) method and the GPE reveals a close correspondence between the many-body effect of depletion and the mean-field effect of wave chaos for a wide range of single-particle external potentials. In the regime of wave chaos the GPE fails to account for the fine-scale quantum fluctuations because many-body effects beyond the validity of the GPE are non-negligible. Surprisingly, despite the failure of the GPE to account for the depletion, coarse grained expectation values of the single-particle density such as the overall width of the atomic cloud agree very well with the many-body simulations. The time dependent depletion of the condensate could be investigated experimentally, e.g., via decay of coherence of the expanding atom cloud.

preprint2011arXiv

Recursive formulation of the multiconfigurational time-dependent Hartree method for fermions, bosons and mixtures thereof in terms of one-body density operators

The multiconfigurational time-dependent Hartree method (MCTDH) [Chem. Phys. Lett. {\bf 165}, 73 (1990); J. Chem. Phys. {\bf 97}, 3199 (1992)] is celebrating nowadays entering its third decade of tackling numerically-exactly a broad range of correlated multi-dimensional non-equilibrium quantum dynamical systems. Taking in recent years particles' statistics explicitly into account, within the MCTDH for fermions (MCTDHF) and for bosons (MCTDHB), has opened up further opportunities to treat larger systems of interacting identical particles, primarily in laser-atom and cold-atom physics. With the increase of experimental capabilities to simultaneously trap mixtures of two, three, and possibly even multiple kinds of interacting composite identical particles together, we set up the stage in the present work and specify the MCTDH method for such cases. Explicitly, the MCTDH method for systems with three kinds of identical particles interacting via all combinations of two- and three-body forces is presented, and the resulting equations-of-motion are briefly discussed. All four possible mixtures of fermions and bosons are presented in a unified manner. Particular attention is paid to represent the coefficients' part of the equations-of-motion in a compact recursive form in terms of one-body density operators only. The recursion utilizes the recently proposed Combinadic-based mapping for fermionic and bosonic operators in Fock space [Phys. Rev. A {\bf 81}, 022124 (2010)] and successfully applied and implemented within MCTDHB. Our work sheds new light on the representation of the coefficients' part in MCTDHF and MCTDHB without resorting to the matrix elements of the many-body Hamiltonian with respect to the time-dependent configurations. It suggests a recipe for efficient implementation of the schemes derived here for mixtures which is suitable for parallelization.

preprint2010arXiv

Mechanism of Tunneling in Interacting Open Ultracold Few-Boson Systems

We investigate the mechanism in the tunneling dynamics of open ultracold few-boson systems by numerically solving the time-dependent few-boson Schrödinger equation exactly. By starting from a weakly repulsive, initially coherent two-boson system we demonstrate that the decay dynamics incorporate fragmentation. The wavefunction of the tunneling state exhibits a pronounced dynamically-stable pattern which we explain by an analytical model. By studying more bosons and stronger interactions we arrive to the conclusion that the decay by tunneling is not a coherent process and exhibits a wealth of phenomena depending on the interaction between the particles.