We extent the recently developed Multi-Layer Multi-Configuration Time-Dependent Hartree method for Bosons (ML-MCTDHB) for simulating the correlated quantum dynamics of bosonic mixtures to the fermionic sector and establish a unifying approach for the investigation of the correlated quantum dynamics of mixture of indistinguishable particles, be it fermions or bosons.Relying on a multi-layer wave-function expansion, the resulting Multi-Layer Multi-Configuration Time-Dependent Hartree method for Mixtures (ML-MCTDHX) can be adapted to efficiently resolve system-specific intra-and inter-species correlations. The versatility and efficiency of ML-MCTDHX is demonstrated by applying it to the problem of colliding few-atom mixtures of both Bose-Fermi and Fermi-Fermi type. Thereby, we elucidate the role of correlations in the transmission and reflection properties of the collisional events. In particular, we present examples where the reflection (transmission) at the other atomic species is a correlation-dominated effect, i.e. it is suppressed in the mean-field approximation.

We develop the multi-layer multi-configuration time-dependent Hartree method for bosons (ML-MCTDHB), a variational numerically exact ab initio method for studying the quantum dynamics and stationary properties of general bosonic systems. ML-MCTDHB takes advantage of the permutation symmetry of identical bosons, which allows for investigations of the quantum dynamics from few to many-body systems. Moreover, the multi-layer feature enables ML-MCTDHB to describe mixed bosonic systems consisting of arbitrary many species. Multi-dimensional as well as mixed-dimensional systems can be accurately and efficiently simulated via the multi-layer expansion scheme. We provide a detailed account of the underlying theory and the corresponding implementation. We also demonstrate the superior performance by applying the method to the tunneling dynamics of bosonic ensembles in a one-dimensional double well potential, where a single-species bosonic ensemble of various correlation strengths and a weakly interacting two-species bosonic ensemble are considered.

We develop and apply the multi-layer multi-configuration timedependent Hartree method for bosons, which represents an ab initio method for investigating the non-equilibrium quantum dynamics of multi-species bosonic systems. Its multi-layer feature allows for tailoring the wave function ansatz to describe intra-and inter-species correlations accurately and efficiently. To demonstrate the beneficial scaling and efficiency of the method, we explored the correlated tunneling dynamics of two species with repulsive intra-and inter-species interactions, to which a third species with vanishing intra-species interaction was weakly coupled. The population imbalances of the first two species can feature a temporal equilibration and their time evolution significantly depends on the coupling to the third species. Bosons of the first and second species exhibit a bunching tendency, whose strength can be influenced by their coupling to the third species.