BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Ofir E. Alon

Ofir E. Alon appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.quant-gasquant-phphysics.atom-phcond-mat.other

Catalog footprint

What is connected

29works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

29 published item(s)

preprint2021arXiv

Longitudinal and transversal resonant tunneling of interacting bosons in a two-dimensional Josephson junction: Mean-field and many-body dynamics

We unravel the out-of-equilibrium quantum dynamics of a few interacting bosonic clouds in a two-dimensional asymmetric double-well potential at the resonant tunneling scenario. At the single-particle level of resonant tunneling, particles tunnel under the barrier from, typically, the ground-state in the left well to an excited state in the right well, i.e., states of different shapes and properties are coupled when their one-particle energies coincide. In two spatial dimensions, two types of resonant tunneling processes are possible, to which we refer to as longitudinal and transversal resonant tunneling. Longitudinal resonant tunneling implies that the state in the right well is longitudinally-excited with respect to the state in the left well, whereas transversal resonant tunneling implies that the former is transversely-excited with respect to the latter. We show that interaction between bosons makes resonant tunneling phenomena in two spatial dimensions profoundly rich, and analyze these phenomena in terms of the loss of coherence of the junction and development of fragmentation, and coupling between transverse and longitudinal degrees-of-freedom and excitations. To this end, a detailed analysis of the tunneling dynamics is performed by exploring the time evolution of a few physical quantities, namely, the survival probability, occupation numbers of the reduced one-particle density matrix, and the many-particle position, momentum, and angular-momentum variances. In general, we display the impact of the transversal and longitudinal degrees-of-freedom in the many-boson tunneling dynamics at the resonant tunneling scenarios.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Many-body effects in the excitations and dynamics of trapped Bose-Einstein condensates

This review explores the dynamics and the low-energy excitation spectra of Bose-Einstein condensates (BECs) of interacting bosons in external potential traps putting particular emphasis on the emerging many-body effects beyond mean-field descriptions. To do so, methods have to be used that, in principle, can provide numerically exact results for both the dynamics and the excitation spectra in a systematic manner. Numerically exact results for the dynamics are presented employing the well-established multicongurational time-dependent Hartree for bosons (MCTDHB) method. The respective excitation spectra are calculated utilizing the more recently introduced linear-response theory atop it (LR-MCTDHB). The latter theory gives rise to an, in general, non-hermitian eigenvalue problem. The theory and its newly developed implementation are described in detail and benchmarked towards the exactly-solvable harmonic-interaction model. Several applications to BECs in one- and two-dimensional potential traps are discussed. With respect to dynamics, it is shown that both the out-of-equilibrium tunneling dynamics and the dynamics of trapped vortices are of many-body nature. Furthermore, many-body effects in the excitation spectra are presented for BECs in different trap geometries. It is demonstrated that even for essentially-condensed systems, the spectrum of the lowest-in-energy excitations computed at the many-body level can differ substantially from the standard mean-field description. In general, it is shown that bosons carrying angular momentum are more sensitive to many-body effects than bosons without. These effects are present in both the dynamics and the excitation spectrum.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Impact of the transverse direction on the many-body tunneling dynamics in a two-dimensional bosonic Josephson junction

Tunneling in a many-body system appears as one of the novel implications of quantum physics, in which particles move in space under an otherwise classically-forbidden potential barrier. Here, we theoretically describe the quantum dynamics of the tunneling phenomenon of a few intricate bosonic clouds in a closed system of a two-dimensional symmetric double-well potential. We examine how the inclusion of the transverse direction, orthogonal to the junction of the double-well, can intervene in the tunneling dynamics of bosonic clouds. We use a well-known many-body numerical method, called the multiconfigurational time-dependent Hartree for bosons (MCTDHB) method. MCTDHB allows one to obtain accurately the time-dependent many-particle wavefunction of the bosons which in principle entails all the information of interest about the system under investigation. We analyze the tunneling dynamics by preparing the initial state of the bosonic clouds in the left well of the double-well either as the ground, longitudinally or transversely excited, or a vortex state. We unravel the detailed mechanism of the tunneling process by analyzing the evolution in time of the survival probability, depletion and fragmentation, and the many-particle position, momentum, and angular-momentum expectation values and their variances. As a general rule, all objects lose coherence while tunneling through the barrier and the states which include transverse excitations do so faster. Implications are briefly discussed.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Many-body tunneling dynamics of Bose-Einstein condensates and vortex states in two spatial dimensions

In this work, we study the out-of-equilibrium many-body tunneling dynamics of a Bose-Einstein condensate in a two-dimensional radial double well. We investigate the impact of interparticle repulsion and compare the influence of angular momentum on the many-body tunneling dynamics. Accurate many-body dynamics are obtained by solving the full many-body Schrödinger equation. We demonstrate that macroscopic vortex states of definite total angular momentum indeed tunnel and that, even in the regime of weak repulsions, a many-body treatment is necessary to capture the correct tunneling dynamics. As a general rule, many-body effects set in at weaker interactions when the tunneling system carries angular momentum.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The uncertainty product of an out-of-equilibrium many-particle system

In the present work we show, analytically and numerically, that the variance of many-particle operators and their uncertainty product for an out-of-equilibrium Bose-Einstein condensate (BEC) can deviate from the outcome of the time-dependent Gross-Pitaevskii dynamics, even in the limit of infinite number of particles and at constant interaction parameter when the system becomes 100% condensed. We demonstrate our finding on the dynamics of the center-of-mass position--momentum uncertainty product of a freely expanding as well as of a trapped BEC. This time-dependent many-body phenomenon is explained by the existence of time-dependent correlations which manifest themselves in the system's reduced two-body density matrix used to evaluate the uncertainty product. Our work demonstrates that one has to use a many-body propagation theory to describe an out-of-equilibrium BEC, even in the infinite particle limit.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Variance as a sensitive probe of correlations enduring the infinite particle limit

Bose-Einstein condensates made of ultracold trapped bosonic atoms have become a central venue in which interacting many-body quantum systems are studied. The ground state of a trapped Bose-Einstein condensate has been proven to be 100% condensed in the limit of infinite particle number and constant interaction parameter [Lieb and Seiringer, Phys. Rev. Lett. {\bf 88}, 170409 (2002)]. The meaning of this result is that properties of the condensate, noticeably its energy and density, converge to those obtained by minimizing the Gross-Pitaevskii energy functional. This naturally raises the question whether correlations are of any importance in this limit. Here, we demonstrate both analytically and numerically that even in the infinite particle limit many-body correlations can lead to a substantial modification of the \textit{variance} of any operator compared to that expected from the Gross-Pitaevskii result. The strong deviation of the variance stems from its explicit dependence on terms of the reduced two-body density matrix which otherwise do not contribute to the energy and density in this limit. This makes the variance a sensitive probe of many-body correlations even when the energy and density of the system have already converged to the Gross-Pitaevskii result. We use the center-of-mass position operator to exemplify this persistence of correlations. Implications of this many-body effect are discussed.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Many-body excitation spectra of trapped bosons with general interaction by linear response

The linear-response theory of the multiconfigurational time-dependent Hartree for bosons method for computing many-body excitations of trapped Bose-Einstein condensates [Phys. Rev. A {\bf 88}, 023606 (2013)] is implemented for systems with general interparticle interaction. Illustrative numerical examples for repulsive and attractive bosons are provided. The many-body linear-response theory identifies the excitations not unraveled within Bogoliubov--de Gennes equations. The theory is herewith benchmarked against the exactly-solvable one-dimensional harmonic-interaction model. As a complementary result, we represent the theory in a compact block-diagonal form, opening up thereby an avenue for treating larger systems.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Spatially partitioned many-body vortices

A vortex in Bose-Einstein condensates is a localized object which looks much like a tiny tornado storm. It is well described by mean-field theory. In the present work we go beyond the current paradigm and introduce many-body vortices. These are made of {\it spatially-partitioned} clouds, carry definite total angular momentum, and are fragmented rather than condensed objects which can only be described beyond mean-field theory. A phase diagram based on a mean-field model assists in predicting the parameters where many-body vortices occur. Implications are briefly discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Universality of Fragmentation in the Schrödinger Dynamics of Bosonic Josephson Junctions

The many-body Schrödinger dynamics of a one-dimensional bosonic Josephson junction is investigated for up to ten thousand bosons and long times. The initial states are fully condensed and the interaction strength is weak. We report on a universal fragmentation dynamics on the many-body level: systems consisting of different numbers of particles fragment to the same value at constant mean-field interaction strength. The phenomenon manifests itself in observables such as the correlation functions of the system. We explain this universal fragmentation dynamics analytically based on the Bose-Hubbard model. We thereby show that the extent to which many-body effects become important at later times depends crucially on the initial state. Even for arbitrarily large particle numbers and arbitrarily weak interaction strength the dynamics is many-body in nature and the fragmentation universal. There is no weakly interacting limit where the Gross-Piatevskii mean-field is valid for long times.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Generic regimes of quantum many-body dynamics of trapped bosonic systems with strong repulsive interactions

Two generically different but universal dynamical quantum many-body behaviors are discovered by probing the stability of trapped fragmented bosonic systems with strong repulsive finite/long range inter-particle interactions. We use different time-dependent processes to destabilize the systems -- a sudden displacement of the trap is accompanied by a sudden quench of the strength of the inter-particle repulsion. A rather moderate non-violent evolution of the density in the first "topology-preserved" scenario is contrasted with a highly-non-equilibrium dynamics characterizing an explosive changes of the density profiles in the second scenario. The many-body physics behind is identified and interpreted in terms of self-induced time-dependent barriers governing the respective under- and over-a-barrier dynamical evolutions. The universality of the discovered scenarios is explicitly confirmed in 1D, 2D and 3D many-body computations in (a)symmetric traps and repulsive finite/long range inter-particle interaction potentials of different shapes. Implications are briefly discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Unified view on linear response of interacting identical and distinguishable particles from multiconfigurational time-dependent Hartree methods

A unified view on linear response of interacting systems utilizing multicongurational time-dependent Hartree methods is presented. The cases of one-particle and two-particle response operators for identical particles and up to all-system response operators for distinguishable degrees-of-freedom are considered. The working equations for systems of identical bosons (LR-MCTDHB) and identical fermions (LR-MCTDHF), as well for systems of distinguishable particles (LR-MCTDH) are explicitly derived. These linear-response theories provide numerically-exact excitation energies and system's properties, when numerical convergence is achieved in the calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Dynamics and symmetries of a repulsively bound atom pair in an infinite optical lattice

We investigate the dynamics of two bosons trapped in an infinite one-dimensional optical lattice potential within the framework of the Bose-Hubbard model and derive an exact expression for the wavefunction at finite time. As initial condition we chose localized atoms that are separated by a distance of $d$ lattice sites and carry a center of mass quasi-momentum. An initially localized pair ($d=0$) is found to be more stable as quantified by the pair probability (probability to find two atoms at the same lattice site) when the interaction and/or the center of mass quasi-momentum is increased. For initially separated atoms ($d \neq 0$) there exists an optimal interaction strength for pair formation. Simple expressions for the wavefunction, the pair probability and the optimal interaction strength for pair formation are computed in the limit of infinite time. Whereas the time-dependent wavefunction differs for values of the interaction strength that differ only by the sign, important observables like the density and the pair probability do not. With a symmetry analysis this behavior is shown to extend to the $N$-particle level and to fermionic systems. Our results provide a complementary understanding of the recently observed [Winkler \textit{et al.}, Nature (London) \textbf{441}, 853 (2006)] dynamical stability of atom pairs in a repulsively interacting lattice gas.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Excitation spectra of fragmented condensates by linear response: General theory and application to a condensate in a double-well potential

Linear response of simple (i.e., condensed) Bose-Einstein condensates is known to lead to the Bogoliubov- de Gennes equations. Here, we derive linear response for fragmented Bose-Einstein condensates, i.e., for the case where the many-body wave function is not a product of one, but of several single-particle states (orbitals). Our approach is based on the number-conserving variational time-dependent mean field theory [O. E. Alon, A. I. Streltsov, and L. S. Cederbaum, Phys. Lett. A 362, 453 (2007)], which describes the time evolution of best-mean field states. Correspondingly, we call our linear response theory for fragmented states LR-BMF. In the derivation it follows naturally that excitations are orthogonal to the ground-state orbitals. As applications excitation spectra of Bose-Einstein condensates in double-well potentials are calculated. Both symmetric and asymmetric double-wells are studied for several interaction strengths and barrier heights. The cases of condensed and two-fold fragmented ground states are compared. Interestingly, even in such situations where the response frequencies of the two cases are computed to be close to each other, which is the situation for the excitations well below the barrier, striking differences in the density response in momentum space are found. For excitations with an energy of the order of the barrier height, both the energies and the density response of condensed and fragmented systems are very different. In fragmented systems there is a class of "swapped" excitations where an atom is transfered to the neighboring well. The mechanism of its origin is discussed. In asymmetric wells, the response of a fragmented system is purely local (i.e., finite in either one or the other well) with different frequencies for the left and right fragments.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Two trapped particles interacting by a finite-ranged two-body potential in two spatial dimensions

We examine the problem of two particles confined in an isotropic harmonic trap, which interact via a finite-ranged Gaussian-shaped potential in two spatial dimensions. We derive an approximative transcendental equation for the energy and study the resulting spectrum as a function of the interparticle interaction strength. Both the attractive and repulsive systems are analyzed. We study the impact of the potential's range on the ground-state energy. Complementary, we also explicitly verify by a variational treatment that in the zero-range limit the positive delta potential in two dimensions only reproduces the non-interacting results, if the Hilbert space in not truncated. Finally, we establish and discuss the connection between our finite-range treatment and regularized zero-range results from the literature.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Number fluctuations of cold spatially split bosonic objects

We investigate the number fluctuations of spatially split many-boson systems employing a theorem about the maximally and minimally attainable variances of an observable. The number fluctuations of many-boson systems are given for different numbers of lattice sites and both mean-field and many-body wave functions. It is shown which states maximize the particle number fluctuations, both in lattices and double-wells. The fragmentation of the states is discussed, and it is shown that the number fluctuations of some fragmented states are identical to those of fully condensed states.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Optimal time-dependent lattice models for nonequilibrium dynamics

Lattice models are abundant in theoretical and condensed-matter physics. Generally, lattice models contain time-independent hopping and interaction parameters that are derived from the Wannier functions of the noninteracting problem. Here, we present a new concept based on time-dependent Wannier functions and the variational principle that leads to optimal time-dependent lattice models. As an application, we use the Bose-Hubbard model with time-dependent Wannier functions to study a quench scenario involving higher bands. We find a separation of times scales in the dynamics and show that under some circumstances the multi-band nonequilibrium dynamics of a quantum system can be obtained essentially at the cost of a single-band model.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Fragmented Many-Body states of definite angular momentum and stability of attractive 3D Condensates

A three dimensional attractive Bose-Einstein Condensate (BEC) is expected to collapse, when the number of the particles $N$ in the ground state or the interaction strength $λ_0$ exceeds a critical value. We study systems of different particle numbers and interaction strength and find that even if the overall ground state is collapsed there is a plethora of fragmented excited states that are still in the metastable region. Utilizing the \emph{configuration interaction} expansion we determine the spectrum of the ground (`yrast') and excited many-body states with definite total angular momentum quantum numbers $0\leqslant L\leqslant N$ and $-L\leqslant M_L\leqslant L$, and we find and examine states that survive the collapse. This opens up the possibility of realizing a metastable system with overcritical numbers of bosons in a ground state with angular momentum $L\neq0$. The multi-orbital mean-field theory predictions about the existence of fragmented metastable states with overcritical numbers of bosons are verified and elucidated at the many-body level. The descriptions of the total angular momentum within the mean-field and the many-body approaches are compared.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A novel general mapping for bosonic and fermionic operators in Fock space

In this paper we provide a novel and general way to construct the result of the action of any bosonic or fermionic operator represented in second quantized form on a state vector, without resorting to the matrix representation of operators and even to its elements. The new approach is based on our proposal to compactly enumerate the configurations (i.e., determinants for fermions, permanents for bosons) which are the elements of the state vector. This extremely simplifies the calculation of the action of an operator on a state vector. The computations of statical properties and of the evolution dynamics of a system become much more efficient and applications to systems made of more particles become feasible. Explicit formulations are given for spin-polarized fermionic systems and spinless bosonic systems, as well as to general (two-component) fermionic systems, two-component bosonic systems, and mixtures thereof.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Accurate multi-boson long-time dynamics in triple-well periodic traps

To solve the many-boson Schrödinger equation we utilize the Multiconfigurational time-dependent Hartree method for bosons (MCTDHB). To be able to attack larger systems and/or to propagate the solution for longer times, we implement a parallel version of the MCTDHB method thereby realizing the recently proposed [Streltsov {\it et al.} arXiv:0910.2577v1] novel idea how to construct efficiently the result of the action of the Hamiltonian on a bosonic state vector. We study the real-space dynamics of repulsive bosonic systems made of N=12, 51 and 3003 bosons in triple-well periodic potentials. The ground state of this system is three-fold fragmented. By suddenly strongly distorting the trap potential, the system performs complex many-body quantum dynamics. At long times it reveals a tendency to an oscillatory behavior around a threefold fragmented state. These oscillations are strongly suppressed and damped by quantum depletions. In spite of the richness of the observed dynamics, the three time-adaptive orbitals of MCTDHB(M=3) are capable to describe the many-boson quantum dynamics of the system for short and intermediate times. For longer times, however, more self-consistent time-adaptive orbitals are needed to correctly describe the non-equilibrium many-body physics. The convergence of the MCTDHB($M$) method with the number $M$ of self-consistent time-dependent orbitals used is demonstrated.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Quantum dynamics of attractive versus repulsive bosonic Josephson junctions: Bose-Hubbard and full-Hamiltonian results

The quantum dynamics of one-dimensional bosonic Josephson junctions with attractive and repulsive interparticle interactions is studied using the Bose-Hubbard model and by numerically-exact computations of the full many-body Hamiltonian. A symmetry present in the Bose-Hubbard Hamiltonian dictates an equivalence between the evolution in time of attractive and repulsive Josephson junctions with attractive and repulsive interactions of equal magnitude. The full many-body Hamiltonian does not possess this symmetry and consequently the dynamics of the attractive and repulsive junctions are different.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Pathway from condensation via fragmentation to fermionization of cold bosonic systems

For small scattering lengths, cold bosonic atoms form a condensate the density profile of which is smooth. With increasing scattering length, the density {\it gradually} acquires more and more oscillations. Finally, the number of oscillations equals the number of bosons and the system becomes {\it fermionized}. On this pathway from condensation to fermionization intriguing phenomena occur, depending on the shape of the trap. These include macroscopic fragmentation and {\it coexistence} of condensed and fermionized parts that are separated in space.

0 citations0 reviewsNo code linkedOpen access