# Phantom vortices: hidden angular momentum in ultracold dilute Bose-Einstein condensates

###### Abstract

Vortices are essential to angular momentum in quantum systems such as
ultracold atomic gases.
The existence of quantized vorticity in bosonic systems stimulated the development of the Gross-Pitaevskii mean-field approximation.
However, the true dynamics of angular momentum in finite, interacting many-body systems like trapped Bose-Einstein condensates is enriched by the emergence of quantum correlations whose description demands more elaborate methods.
Herein we theoretically investigate the full many-body dynamics of the acquisition of angular momentum by a gas of ultracold bosons in two dimensions using a standard rotation procedure.
We demonstrate the existence of a novel mode of quantized vorticity, which we term the phantom vortex that, contrary to the conventional mean-field vortex, can be detected as a topological
defect of spatial coherence, but not of the density.
We describe previously unknown many-body mechanisms of vortex nucleation and show that angular momentum is hidden in phantom vortex modes which so far seem to have evaded experimental detection.

###### pacs:

03.75Lm, 03.75.Kk, 03.75.Nt, 05.30.Jp## Introduction

Quantized vortices are perhaps the most interesting way angular momentum is known to manifest in quantum many-body systems Pethick (); Ketterle2001 (). They appear in a variety of experiments with systems including atomic Bose-Einstein condensates (BECs) Vortex_Exp (); Vortex_Exp2 (); Ketterle2001 (); Kali (), and exciton-polariton condensates Baas2008 (). Typically, a quantum vortex is characterized by a density node and a phase discontinuity. Other known signatures of angular momentum, such as center-of-mass rotation, surface waves, and quadrupole modes are also features of the density profile Dalibard2000 (); Pethick (). However, it is well known that there exist aspects of many-body dynamics, like fragmentation James1982 (); Spekkens1999 (), that may not be visible in the density glauber (); Kaspar2008 (). Coreless vortices, i.e., vortices invisible in the density, are known to exist in mean-field spinor condensates Lovegrove2014 (); Kasamatsu2003 () and have been recently found in a multicomponent condensate Liu2013 (). Since the core of a vortex in one (spin-)component is filled with the other (spin-)component, coreless vortices are not visible in the density. However, they have not yet been seen in single-component condensates with no spin. Note that experiments studying vortex nucleation in rotating BECs have thus far observed a critical frequency below which no vortices are detected, see for instance Ketterle2001 (); Hodby (); Feder (). As we shall show in the following, however, vortices without cores are present in single-component BECs even below the critical rotation velocity. These coreless vortices, to date, have evaded experimental detection as they cannot be seen from the density.

Traditionally, the dynamics of BECs are aptly treated using the time-dependent Gross-Pitaevskii (GP) equation Dalfovo1999 () which assumes that the many-body wavefunction is coherent and condensed for all time. Naturally, this assumption neglects the existence of quantum correlations and precludes the possibility of fragmentation, i.e. macroscopic occupation of more than one single particle state Onsager1956 (). In other computational investigations, the rotation frequency was modulated at fixed anisotropy which also leads to instabilities and formation of vortex lattices, within the GP theory Tsubota2002 (); Suominen2003 (); Parker2005 (), the Bogoliubov approach Lobo2004 () and beyond Wright2008 (). More recently, much effort has been devoted to exploring the role of fragmentation in stationary states of bosonic systems Spekkens1999 (); Cederbaum2003 (); Streltsov2006 (); Mueller2006 (); Lewenstein2009 (); Fischer2009 (); Fischer2014 () including those with spin degrees of freedom Qi2013 (); Song2014 () using methods such as best mean-field theory or general variational many-body approaches. In particular, it has been shown that there exist energy eigenstates with definite angular momentum that are fragmented for repulsive interactions Mueller2006 () as well as in two dimensions for attractive interaction Marios2014 (). Furthermore, in the study of adiabatic vortex nucleation in few-body systems Lewenstein2009 (), it was shown that fragmentation is unavoidable. Generally, nontrivial dynamics involve many eigenstates, each of which may be fragmented. Hence, it stands to reason that a true many-body method is necessary to describe the dynamics of a two-dimensional condensate as it acquires angular momentum.

An appropriate many-body method is the multiconfigurational time-dependent Hartree for bosons (MCTDHB) MCTDHB () method, which has been shown to accurately capture the transition from coherence to fragmentation with great success in a variety of systems Alexej2007 (); Axel2012 (); streltsova2014 (). MCTDHB provides highly accurate results even for problems with time-dependent traps and time-dependent interparticle interactions Lode2012 (); AxelThesis ().

Here we investigate the dynamics of interacting bosons in two dimensions by numerically solving the time-dependent many-body Schrödinger equation for a standard rotation scenario using the recursive implementation of MCTDHB in the MCTDH-X software package RMCTDHB (). Our work is among the first that explore the many-body dynamics of this type of system beyond mean-field theory MariosJLTP (). The system that we describe here is similar to that simulated in Ref. Suominen2003 () using the time-dependent GP equation, and analyzed experimentally in Ref. Madison2000 (). There, regions in the parameter space were established for which vortices are observed. Even though the GP treatment for the parameters chosen herein predicts no vortex nucleation (Supplementary Information), the many-body analysis exhibits rich vortex dynamics that are invisible in the density. Our results demonstrate that angular momentum in fragmented condensates manifests itself in a new type of coreless vortex, which we name the phantom vortex due to its elusive nature. It is important to stress that we found phantom vortices for a wide range of particle numbers, interaction strength, interaction range, trap anisotropy, and rotation frequencies, see (Supplementary Information). An analysis of the energy eigenstates of a six-boson system in the corotating frame has shown signatures of phantom vortices (see Fig. 3 in Ref. Lewenstein2009 ()). Yet, their importance in the dynamical acquisition of angular momentum by many-body systems, even in the parameter regime where no vortices nucleate in the density, has been overlooked.

## System Description

We begin our analysis, in silico, by first computing the many-body ground state of interacting bosons in an isotropic harmonic trap. Subsequently, we gently transfer angular momentum into the system by elliptically deforming the harmonic trap while rotating the axis of anisotropy at a fixed frequency. The anisotropy parameter is varied in a piecewise-linear fashion; first, it is ramped from zero to a maximum value, then held constant over a fixed time interval, and finally ramped back to zero. The trap is then kept isotropic for the remainder of the simulation. See Fig. 1a and Supplementary Information for details.

The parameters chosen here correspond to rotating, weakly interacting, zero-temperature atoms, confined by a pancake-like trap. The particle density and interactions are similar to the small (i.e. ) Rb condensate reported in Schmied2016 (). The subcritical rotation frequency chosen herein is not fast enough to nucleate quantized vortices in the density. A GP () simulation with the same parameters only absorbs of angular momentum per particle (see Fig 1.d).

The dynamics of an interacting gas of bosons are governed by the time-dependent many-body Schrödinger equation

(1) |

We solve Eq. 1 using the MCTDHB ansatz, written in second quantization as

(2) |

That is, the bosons dynamically populate all configurations of particles in time-dependent variationally optimized single-particle states, called orbitals (Supplementary Information). The case where and is analyzed in the main text of this work whereas other are discussed in the Supplementary Information. In our setting the particle density at rest () is close to a Gaussian-like density distribution because the bosons are weakly interacting. The temperature is absolute zero throughout the dynamics and the external frequency of rotation of the mild anisotropy is at of the trapping frequency.

The many-body Hamiltonian in dimensionless units reads:
^{†}^{†}footnotetext: To convert to dimensionful units, the Hamiltonian is multiplied by . With the mass Kg corresponding to , and length scale m, we get a transverse trapping frequency of and a time scale of ms. With this choice of scale, corresponds to , and the linear extent of the simulation is m. The interaction parameter is related to the scattering length and transverse confinement by . For , . Using this scattering length and gives kHz for . So the transverse aspect ratio is and the maximum energy per particle, kHz is merely . Thus, we can truly treat the dynamics as two-dimensional.

(3) |

We use as a two-body potential a normalized Gaussian of width and interaction strength . This choice of a short-range interaction is motivated by the similarity of the physics for zero- and short-ranged interaction potentials in Ref. rosti2013 (). See Supplemental Material for a direct comparison. The one-body Hamiltonian is the sum of the kinetic energy and the trapping potential,

(4) |

where

(5) |

and . The anisotropy parameter varies in time so that it is ramped up from to over , held fixed at for , ramped down to zero over , and then held fixed at for the rest of simulation until as shown in Fig. 1a. The resulting healing length, is comparable to the oscillator length (1 in dimensionless units) so the particles are indeed weakly interacting.

The primary object of our analysis is the one-body reduced density matrix (RDM), which is defined as the partial trace of the -body density Kaspar2008 (); glauber ():

(6) |

with normalized to . The trace operation can be understood as eliminating knowledge of the many-body effects to obtain a single-particle operator, the RDM. The diagonal, , of the RDM is the particle density, , or simply the density. The RDM is written in its eigenbasis as

(7) |

The eigenvalues, , are ordered in decreasing magnitude and called natural occupations and the eigenfunctions, , are called the natural orbitals or synonymously, fragments. Note that the density, , is equal to the sum of the squared amplitudes of the natural orbitals, weighted with their natural occupations.

A state is considered coherent when only one natural orbital has significant occupation. In this case, the full information of the many-body wavefunction is contained in a single one-particle state and thus, the Gross-Pitaevskii mean-field approximation is valid. Otherwise, multiple orbitals have macroscopic occupation and the state is said to be fragmented because the density is a sum of multiple single-particle functions. We emphasize that, from the view of the RDM, the MCTDHB approach generalizes the mean-field approximation by allowing a dynamical transition from coherence to fragmentation. This means that, contrary to the mean-field case, the natural occupations are allowed to have non-integer values and vary in time.

Importantly, MCTDHB boils down to the GP mean-field for the case of , because then the ansatz for the method, Eq. (2), becomes identical to the ansatz of the GP mean-field. This allows us to straightforwardly compare the predictions of the two approaches.

Another way of describing coherence is through the first-order correlation function, , defined in terms of the RDM as

(8) |

A coherent state at time must satisfy , while any other value indicates fragmentation glauber (). Intuitively, measures how well the many-body wavefunction is described by a single particle state.

The off-diagonal, , terms in and the orbitals are complex-valued. Thus it is instructive to plot their magnitudes and phases separately. We hence define the phase of to be

(9) |

and the orbital phases to be

(10) |

For the sake of completeness, we likewise define the phase of the many-body wavefunction, , as

(11) |

## Simulation Results

We have found that there may exist many vortices in the natural orbitals despite there being no visible vortices in the density (Figs. 2 and 3), even with angular momentum per particle (Fig. 1). We name the vortices which exist in the natural orbitals phantom vortices. These phantom vortices persist on long time scales compared to the period of the harmonic trap. Phantom vortices that exist for long times near the center of the trap either nucleate on existing topological defects present in the initial natural orbitals or nucleate via a transfer of vorticity between natural orbitals (Video S1). Since the simulation starts with a condensed state and both mechanisms of phantom vortex nucleation rely on the dynamics and occupation of several fragments, they are inherently many-body phenomena that cannot be described by mean-field methods.

We now give a chronological description of the evolution of particles with orbitals starting in the ground state at to (see also Video S1 in Supplementary Information). At , resembled 1s, 2p, 2p, and 2s orbitals, respectively. The system began with occupation in and thus was almost entirely condensed (Fig. 1b). As the anisotropy was ramped up, deformed smoothly and rotated with the potential. By , due to a reordering of the natural occupations, and switched labels. At this time, , , and already showed faint hints of vortices at the edge of their density similar to those seen at later times in Figs. 2 and 3. These vortices are the fragmented counterparts of ghost vortices reported in Tsubota2002 (). Ghost vortices are distinct from phantom vortices. Ghost vortices are phase defects outside the bulk of the condensate density. Although they do not contribute significantly to the energy or angular momentum, they may – through interference effects – be responsible for surface waves Tsubota2002 () which are also observed in the present work. In contrast, phantom vortices do contribute significantly to the angular momentum and energy despite being invisible in the density of the system. By , many ghost vortices were established in and . Since the state was still more than condensed, the ghost vortices in manifested in the density and would thus have been detectable in high-fidelity absorption imaging. A one-dimensional cut along a core of a ghost vortex in showed that the outer density maximum was less than of the density maximum at the center of the cloud. By this time, the two initial lobes in had spread out and closed off the angular node into an elliptic shape. This was the first sign of a true phantom vortex in the bulk of the cloud (Fig. 2c). We term this first mechanism of phantom vortex nucleation node mutation, since a node in the fragment deforms and mutates into a phantom vortex. Around , ghost vortices in fused with a distorted angular node to split the two lobes into four with an “I” shaped node (see Video S1 at ). This node mutated to nucleate three persistent phantom vortices which merged into the charge- vortex shown in Fig. 3e,j.

By , the system was significantly fragmented, (Fig. 1b). Fragmentation obscured the vortex-induced density nodes of from being visible in because density from other orbitals filled the vortex cores present in . We comment that although larger particle numbers delay the onset of fragmentation, we have observed phantom vortices in simulations with up to (Supplementary Information). From to , each orbital had a complicated vortex structure while the density was largely featureless (Figs. 2a and 1a). This marks an important aspect of the dynamics: there were many phantom vortices in each orbital, but remarkably no vortices were detectable by directly observing the density (Figs. 2 and 3). In this time interval, the orbital angular momenta, , attained their maximal values and started to decay to their equilibrium values despite the maximal anisotropy.

After , the trap was still maximally anisotropic, yet the system energy and total angular momentum became saturated (Fig. 1d). By , the vortex structures of and reached their steady state. had a single phantom vortex near the center of the cloud which was nucleated from the single node of its initial shape. This node-mutated phantom vortex persisted for the length of the simulation. At , had three phantom vortices in a linear arrangement about the center of the cloud, which persisted for the duration of the anisotropy. This triplet fused into a triply charged phantom vortex at the center of the cloud at when the trap became symmetric again (Fig. 3e,j). Near , and swapped labels due to occupation reordering (see and in Fig. 1b), therefore the steady state of the phantom vortex structure can be seen in Fig. 3b,e,g,j. From on, the labeling of the orbitals remained fixed because the occupation numbers did not change order anymore (cf. Fig. 1b).

At , there were two corotating vortices near the center of the cloud in and none present in . This phantom vortex pair nucleated at the edge of the cloud and then moved towards the center, resembling mean-field vortex nucleation Suominen2003 (); Tsubota2002 (). However, in our treatment, this pair was transient. The phantom vortices then transferred from to , marking a second mechanism unique to phantom vortex nucleation: slow orbital-orbital vorticity transfer, which we now describe in detail. By , there were two prominent ghost vortices in and the phantom vortex pair in had returned from the center to the edge of this fragment. Gradually, the vortices in disappearred by exiting the edge of the orbital while the prominent phantom vortex pair in entered the orbital bulk. By , the phantom vortex nucleation in was complete, marked by the intervortex separation being as small as previously in at , before the transfer (Video S1). This direct interaction between and appears as a strong correlation in the orbital angular momentum, and , as seen from to in Fig. 1c. The transfer of vorticity between fragments occurred on a time scale much slower () than the trap rotation period ().

As the anisotropy was ramped down, the phantom vortex pair in approached the trap center and reached a minimal separation of by when the trap became symmetric. This was the only instance where we observed a phantom vortex nucleate at the edge of a fragment and persist in its bulk (Fig. 3d,i). Simultaneously, during the ramp-down, three singly charged phantom vortices in coalesced into a dynamically stable charge- phantom vortex. This contrasts with the known spontaneous decay of charge- vortices to single vortices Shin2004 () and can only occur for phantom vortices for which the surrounding density of the other fragments may appartently stabilize the higher charge.

Although phantom vortices are not detectable directly from the density, they are strikingly pronounced in the one-body correlation function, (Fig. 4). We have two key observations about the coherence of phantom vortices. First, we fix near the core of a phantom vortex; remarkably, is close to for all near and inside the cores of other phantom vortices. Second, we fix away from the cores of all phantom vortices and again observe that is close to for all near and inside the cores of phantom vortices (Fig. 4, Videos S2 and S3). Even phantom vortices in the same fragment are found to be incoherent. We conclude, therefore, that phantom vortices are incoherent both with respect to each other and the remaining bulk density. This implies that they are distinct objects in the fragmented condensate. Furthermore, phantom vortices are observable since is measurable via interference experiments, e.g. Ref. Hofferberth2007 ().

## Discussion

We have found that fragmentation increases – together with the energy and total angular momentum of the system – even for a rotation that is not fast enough to nucleate vortices in the density . Vortices, as seen in the laboratory (see for instance Hodby (); Fetter2009 (); Hall2010 ()), seem not emerge in our scenario. The GP mean-field description of the same system shows no vortices; however the absorbed angular momentum in that case is roughly three times less than the one predicted by the many-body theory. The angular momentum that the vortex-free GP state posseses is all due to surface excitations and deviations from a symmetric state. The many-body state, however, additionally contains angular momentum in the phantom vortices.

It is a natural follow-up question to ask whether or not the experimentally observed vortices are coherent or fragmented objects. If they are fragmented, phantom vortex cores across all relevant orbitals must coincide. We were able to show that if the coincident phantom vortices are of the same charge, there must be additional phantom vortices, such that the fragments can maintain orthonormality (see extended discussion in Supplementary Information). Furthermore, if all phantom vortices are coincident, they must be of different charge, see for instance and in Fig. 3b,e and Supplementary Information. In either case, and in direct contradiction to the mean-field result, the angular momentum per particle must be greater than unity in order to nucleate a fragmented vortex in . Perhaps some experimentally detected vortices are in fact phantom vortices rather than mean-field vortices, as suggested in Kaspar2016 (). More work is needed in that direction, in order to explore the existence and implication of fragmentation in supercritically rotated gases and the experimental visibility of phantom vortices in single-shots.

In summary, we have observed rich vortex dynamics within individual orbitals of a fragmented single component BEC. Since these vortex dynamics cannot be observed in the density, we termed the vortices in the fragments phantom vortices. We have identified two mechanisms of phantom vortex nucleation, namely node mutation and slow orbital-orbital vortex transfer, which have no mean-field analogue. In node mutation, phantom vortices nucleate on preexisting topological defects in a fragment, whereas in slow orbital-orbital transfer, vortices are transferred between fragments. Phantom vortices are clearly visible in the correlation function. A detailed analysis of the correlation function shows that phantom vortices are completely incoherent both with each other and the bulk density between other phantom vortices. Phantom vortices are thus distinct quantum objects that are experimentally observable, for instance, via interference experiments.

###### Acknowledgements.

We thank Alexej I. Streltsov, Ofir E. Alon, and Tomos Wells for helpful discussions. Computation time on the Cray clusters Hermit, Hornet, and Hazel Hen at the High Performance Computing Center Stuttgart (HLRS), and financial support by the DAAD-RISE program, the Swiss SNF, the NCCR Quantum Science and Technology, DFG, and FAPESP are gratefully acknowledged. We also acknowledge the hospitality of Vanderlei Bagnato and the CEPOF at the IFSC-USP, where part of the work was completed.## References

- (1) Abo-Shaeer, J.R., Raman, C., Vogels, J.M., Ketterle, W. Observation of vortex lattices in Bose-Einstein condensates. Science 292, 476 (2001).
- (2) Pethick, C., Smith, H. Bose-Einstein Condensation in Dilute Gases.
- (3) Donadello, S. et al. Observation of Solitonic Vortices in Bose-Einstein Condensates. Phys. Rev. Lett. 113, 065302 (2014).
- (4) Anderson, B.P. et al. Watching Dark Solitons Decay into Vortex Rings in a Bose-Einstein Condensate. Phys. Rev. Lett. 86, 2926 (2001).
- (5) Wilson, K.E., Newman, Z.L., Lowney, J.D., Anderson, B.P. In situ imaging of vortices in Bose-Einstein condensates. Phys. Rev. A 91, 023621 (2015).
- (6) Lagoudakis, K.G. et al. Quantized vortices in an exciton-polariton condensate. Nature Phys. 4, 706 (2008).
- (7) Chevy, F., Madison, K.W., Dalibard, J. Measurement of the angular momentum of a rotating Bose-Einstein condensate. Phys. Rev. Lett. 85, 2223 (2000).
- (8) Spekkens, R.W., Sipe, J.E. Spatial fragmentation of a Bose-Einstein condensate in a double-well potential. Phys. Rev. A 59, 3868 (1999).
- (9) Nozières, P., James, D.S. Particle vs. pair condensation in attractive Bose liquids. J. Phys. France 43, 1133 (1982).
- (10) Glauber, R.J. Quantum Theory of Optical Coherence (Wiley-VCH Verlag GmbH & Co. KGaA 2007).
- (11) Sakmann, K., Streltsov, A.I., Alon, O.E., Cederbaum, L.S. Reduced density matrices and coherence of trapped interacting bosons. Phys. Rev. A 78, 023615 (2008).
- (12) Kasamatsu K, Tsubota, M., Ueda, M. Vortex Phase Diagram in Rotating Two-Component Bose-Einstein Condensates. Phys. Rev. Lett. 91, 150406 (2003).
- (13) Lovegrove, J., Borgh, M.O., Ruostekoski, J. Energetic stability of coreless vortices in spin-1 Bose-Einstein condensates with conserved magnetization. Phys. Rev. Lett. 112, 075301 (2014).
- (14) Liu, C.F., Fan, H., Gou, S.C., Liu, W.M. Crystallized and amorphous vortices in rotating atomic-molecular Bose-Einstein condensates. Sci. Rep. 4, 4224 (2013).
- (15) David L. Feder, Anatoly A. Svidzinsky, Alexander L. Fetter and Charles W. Clark, Anomalous Modes Drive Vortex Dynamics in Confined Bose-Einstein Condensates. Phys. Let. Rev. 86, 564 (2001).
- (16) Hodby, E., Hechenblaikner, G., Hopkins, S.A., Maragò, O.M., and Foot, C.J., Vortex Nucleation in Bose-Einstein Condensates in an Oblate, Purely Magnetic Potential. Phys. Rev. Lett. 88, 010405 (2002).
- (17) Dalfovo, F., Giorgini, S., Pitaevskii, L.P., Stringari, S. Theory of Bose-Einstein condensation in trapped gases. Rev. Mod. Phys. 71,463 (1999).
- (18) Penrose, O., Onsager, L. Bose-Einstein condensation and liquid helium. Phys. Rev. 104,576 (1956).
- (19) Cederbaum, L.S., Streltsov, A.I. Best mean-field for condensates. Phys. Lett. A 318,564 (2003).
- (20) Streltsov, A.I., Alon, O.E., Cederbaum, L.S. General variational many-body theory with complete self-consistency for trapped bosonic systems. Phys. Rev. A 73,063626 (2006).
- (21) Mueller, E.J., Ho, T.L., Ueda, M., Baym, G. Fragmentation of Bose-Einstein condensates. Phys. Rev. A 74,033612 (2006).
- (22) Dagnino, D., Barberán, N., Lewenstein, M., Dalibard, J. Vortex nucleation as a case study of symmetry breaking in quantum systems. Nature Phys. 5,431 (2009).
- (23) Bader, P., Fischer, U.R. Fragmented many-body ground states for scalar bosons in a single trap. Phys. Rev. Lett. 103,060402 (2009).
- (24) Kang, M.K., Fischer, U.R. Revealing single-trap condensate fragmentation by measuring density-density correlations after time of flight. Phys. Rev. Lett. 113,140404 (2014).
- (25) Zhou, Q., Cui, X. Fate of a Bose-Einstein condensate in the presence of spin-orbit coupling. Phys. Rev. Lett. 110,140407 (2013).
- (26) Song, S.W., Zhang, Y.C., Zhao, H., Wang, X., Liu, W.M. Fragmentation of spin-orbit-coupled spinor Bose-Einstein condensates. Phys. Rev. A 89,063613 (2014).
- (27) Tsatsos, M.C. Attractive Bose gas in two dimensions: An analytical study of its fragmentation and collapse. Phys. Rev. A 89,043604 (2014).
- (28) Alon, O.E., Streltsov, A.I., Cederbaum, L.S. Multiconfigurational time-dependent Hartree method for bosons: Many-body dynamics of bosonic systems. Phys. Rev. A 77,033613 (2008).
- (29) Streltsov, A.I., Alon, O.E., Cederbaum, L.S. Role of excited states in the splitting of a trapped interacting Bose-Einstein condensate by a time-dependent barrier. Phys. Rev. Lett. 99,030402 (2007).
- (30) Lode, A.U.J., Streltsov, A.I., Sakmann, K., Alon, O.E., Cederbaum, L.S. How an interacting many-body system tunnels through a potential barrier to open space. Proc. Natl. Acad. Sci. 109,13521 (2012).
- (31) Streltsova, O.I., Alon, O.E., Cederbaum, L.S., Streltsov A.I. Generic regimes of quantum many-body dynamics of trapped bosonic systems with strong repulsive interactions. Phys. Rev. A (R) 89,061602 (2014).
- (32) Lode, A.U.J., Sakmann, K., Alon, O.E., Cederbaum, L.S., Streltsov, A.I. Numerically exact quantum dynamics of bosons with time-dependent interactions of harmonic type. Phys. Rev. A 86,063606 (2012).
- (33) Lode, A.U.J. Tunneling Dynamics in Open Ultracold Bosonic Systems (Springer (2015).
- (34) Lode, A.U.J., Tsatsos M-C, Fasshauer E, MCTDH-X: The multiconfigurational time-dependent Hartree for indistinguishable particles software, version 0.99. http://ultracold.org; http://mctdhx.org; http://schroedinger.org, http://mctdh.bf (2016).
- (35) Lundh, E., Martikainen, J.P., Suominen, K.A. Vortex nucleation in Bose-Einstein condensates in time-dependent traps. Phys. Rev. A 67, 063604 (2003).
- (36) Madison, K.W., Chevy, F., Wohlleben, W., Dalibard, J. Vortex formation in a stirred Bose-Einstein condensate. Phys. Rev. Lett. 84,806 (2000).
- (37) Tsatsos, M.C. & Lode, A.U.J. Resonances and Dynamical Fragmentation in a Stirred Bose-Einstein Condensate, J. L. Temp. Phys. 181, 171 (2015)
- (38) Parker, N.G. & Adams, C.S. Emergence and Decay of Turbulence in Stirred Atomic Bose-Einstein Condensates. Phys. Rev. Lett. 95, 145301 (2005).
- (39) Tsubota, M., Kasamatsu, K., Ueda, M. Vortex lattice formation in a rotating Bose-Einstein condensate. Phys. Rev. A 65,023603 (2002).
- (40) Lobo, C., Sinatra, A. & Castin, Y. Vortex Lattice Formation in Bose-Einstein Condensates. Phys. Rev. Lett. 92, 020403 (2004).
- (41) Wright, T.M., Ballagh, R.J., Bradley, A.S., Blakie, P.B. & Gardiner, C.W. Dynamical thermalization and vortex formation in stirred two-dimensional Bose-Einstein condensates. Phys. Rev. A 78, 063601 (2008).
- (42) Schmied, R. et al. Bell correlations in a Bose-Einstein condensate. Science 352,441 (2016).
- (43) Doganov, R.A., Klaiman, S., Alon, O.E., Streltsov, A.I., Cederbaum, L.S. Two trapped particles interacting by a finite-range two-body potential in two spatial dimensions. Phys. Rev. A 87,033631 (2013).
- (44) Shin, Y., et al. Dynamical Instability of a Doubly Quantized Vortex in a Bose-Einstein Condensate. Phys. Rev. Lett. 93, 160406 (2004).
- (45) Hofferberth, S., Lesanovsky, I., Fischer, B., Schumm, T., Schmiedmayer, J. Non-equilibrium coherence dynamics in one-dimensional Bose gases. Nature 449,324 (2007).
- (46) Fetter, A.L. Rotating trapped Bose-Einstein condensates. Rev. Mod. Phys. 81,647 (2009).
- (47) Freilich, D.V., Bianchi, D.M., Kaufman, A.M., Langin, T.K., & Hall, D.S. Real-time dynamics of single vortex lines and vortex dipoles in a Bose-Einstein condensate. Science 329,1182 (2010).
- (48) Sakmann, K. & Kasevich, M. Single-shot simulations of dynamic quantum many-body systems. Nature Physics 12,451-454 (2016).

## Author Contribution Statement

S. Weiner ran the simulations analyzed herein, generated the plots figures and videos, and contributed to writing the manuscript. M. C. Tsatsos helped write parts of the manuscript and performed complementary calculations. A. U. J. Lode and L. S. Cederbaum conceived the idea for this project and A. U. J. Lode contributed to writing the manuscript throughout. All authors interpreted the results, reviewed and corrected the manuscript.

## Additional information

All authors declare no competing interests.