Source author record

Alexej I. Streltsov

Alexej I. Streltsov 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

28works
4topics
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

28 published item(s)

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.

preprint2016arXiv

Many-body excitations and de-excitations in trapped ultracold bosonic clouds

We employ the MultiConfiguraional Time-Dependent Hartree for Bosons (MCTDHB) method to study excited states of interacting Bose-Einstein condensates confined by harmonic and double-well trap potentials. Two approaches to access excitations, a static and a dynamic one, have been studied and contrasted. In static simulations the low-lying excitations have been computed by utilizing the LR-MCTDHB method - a linear response theory constructed on-top of a static MCTDHB solution. Complimentary, we propose two dynamic protocols that address excitations by propagating the MCTDHB wave-function. In particular, we investigate dipole-like oscillations induced by shifting the origin of the confining potential and breathing-like excitations by quenching frequency of a parabolic part of the trap. To contrast static predictions and dynamic results we have computed time-evolutions and their Fourier transforms of several local and non-local observables. Namely, we study evolution of the $\left< x(t) \right>$, its variance $\operatorname{Var}(x(t))$, and of a local density computed at a selected position. We found out that the variance is the most sensitive and informative quantity - along with excitations it contains information about the de-excitations even in a linear regime of the induced dynamics. The dynamic protocols are found to access the many-body excitations predicted by the static LR-MCTDHB approach.

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.

preprint2015arXiv

Quantum Speed Limit and Optimal Control of Many-Boson Dynamics

We extend the concept of quantum speed limit -- the minimal time needed to perform a driven evolution -- to complex interacting many-body systems. We investigate a prototypical many-body system, a bosonic Josephson junction, at increasing levels of complexity: (a) within the two-mode approximation {corresponding to} a nonlinear two-level system, (b) at the mean-field level by solving the nonlinear Gross-Pitaevskii equation in a double well potential, and (c) at an exact many-body level by solving the time-dependent many-body Schrödinger equation. We propose a control protocol to transfer atoms from the ground state of a well to the ground state of the neighbouring well. Furthermore, we show that the detrimental effects of the inter-particle repulsion can be eliminated by means of a compensating control pulse, yielding, quite surprisingly, an enhancement of the transfer speed because of the particle interaction -- in contrast to the self-trapping scenario. Finally, we perform numerical optimisations of both the nonlinear and the (exact) many-body quantum dynamics in order to further enhance the transfer efficiency close to the quantum speed limit.

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.

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

Interferometry with correlated matter-waves

Matter-wave interferometry of ultra-cold atoms with attractive interactions is studied at the full many-body level. First, we study how a coherent light-pulse applied to an initially-condensed solitonic system splits it into two matter-waves. The split system looses its coherence and develops correlations with time, and inevitably becomes fragmented due to inter-particle attractions. Next, we show that by re-colliding the sub-clouds constituting the split density together, along with a simultaneous application of the same laser-pulse, one creates three matter-waves propagating with different momenta. We demonstrate that the number of atoms in the sub-cloud with zero-momentum is directly proportional to the degree of fragmentation in the system. This interferometric-based protocol to discriminate, probe, and measure the fragmentation is general and can be applied to ultra-cold systems with attractive, repulsive, short- and long-range interactions.

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.

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.

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.

preprint2013arXiv

Quantum systems of ultra-cold bosons with customized inter-particle interactions

Recent progress in cooling and trapping of polarized clouds of chromium $^{52}Cr$, dysprosium $^{164}Dy$ and erbium $^{168}Er$ opens a road-map to quantum systems where shapes of inter-particle interactions can be customized. The main purpose of this work is to get a deeper insight on a role which the overall shape of the inter-particle interaction plays in a context of trapped ultra-cold bosons. We show that strong inter-particle repulsion inevitably leads to multi-hump fragmentation of the ground state. The fragmentation phenomenon is universal -- it takes place in traps of different dimensionality and topologies and for very broad classes of repulsive inter-particle potentials. The physics behind is identified and explained.

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.

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.

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

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.

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.

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

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.

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.

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

Death of soliton trains in attractive Bose-Einstein condensates

Experiments on ultra-cold attractive Bose-Einstein Condensates (BECs) have demonstrated that at low dimensions atomic clouds can form localized objects, propagating for long times without significant changes in their shapes and attributed to bright matter-wave solitons, which are coherent objects. We consider the dynamics of bright soliton trains from the perspective of many-boson physics. The fate of matter-wave soliton trains is actually to quickly loose their coherence and become macroscopically fragmented BECs. The death of the coherent matter-wave soliton trains gives birth to fragmented objects, whose quantum properties and experimental signatures differ substantially from what is currently assumed.

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.

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.

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.

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.

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.