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

Storm Weiner, Marios Tsatsos, Lorenz Cederbaum, Axel Lode

► To cite this version:

Storm Weiner, Marios Tsatsos, Lorenz Cederbaum, Axel Lode. Phantom vortices: hidden angular momentum in ultracold dilute Bose-Einstein condensates. Scientific Reports, 2017, 7, pp.40122. 10.1038/srep40122 . hal-02369838

HAL Id: hal-02369838

https://hal.science/hal-02369838v1

Submitted on 19 Nov 2019

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
SCIENTIFIC REPORTS

<figure><img src="image_2.png" /></figure>

OPEN

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

Received: 14 June 2016

Accepted: 02 December 2016

Published: 16 January 2017

Storm E. Weiner<fnref n="1" />, Marios C. Tsatsos², Lorenz S. Cederbaum³ & Axel U. J. Lode⁴

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. 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 vortices modes which so far seem to have evaded experimental detection. This phenomenon is likely important in the formation of the Abrikosov lattice and the onset of turbulence in superfluids.

Quantized vortices are perhaps the most interesting way angular momentum is known to manifest in quantum many-body systems<fnref n="1" />,². They appear in a variety of experiments with systems including atomic Bose-Einstein condensates (BECs)<fnref n="1" />,³⁻⁵, and exciton-polariton condensates⁶. 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²,⁷. However, it is well known that there exist aspects of many-body dynamics, like fragmentation⁸,⁹, that may not be visible in the density¹⁰,¹¹. Coreless vortices, i.e., vortices invisible in the density, are known to exist in mean-field spinor condensates¹²,¹³ and have been recently found in a multicomponent condensate¹⁴. 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, coreless vortices 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 refs 1,15,16. 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, seem to 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¹⁷ 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 hyphenate state¹⁸. 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¹⁹⁻²¹, the Bogoliubov approach²² and beyond²³. More recently, much effort has been devoted to exploring the role of fragmentation in stationary states of bosonic systems⁸,²⁴⁻²⁹ including those with spin degrees of freedom³⁰⁻³² 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²⁶ as well as in two dimensions for attractive interaction³³. In the study of adiabatic vortex nucleation in few-body systems²⁷, it was shown that fragmentation is unavoidable. Generally, nontrivial dynamics involve many eigenstates, each of which may be fragmented. Hence, it stands to

<footnote n="1">¹Department of Physics, University of California at Berkeley, CA, USA. ²Instituto de Física de São Carlos, Universidade de São Paulo, São Carlos, São Paulo, Brazil. ³Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, Germany. ⁴Department of Physics, University of Basel, Switzerland. Correspondence and requests for materials should be addressed to S.E.W. (email: StormWeiner@Berkeley.edu)</footnote>
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)³⁴ method, which has been shown to accurately capture the transition from coherence to fragmentation with great success in a variety of systems³⁵⁻³⁷. MCTDHB provides highly accurate results even for problems with time-dependent traps and time-dependent interparticle interactions³⁸,³⁹.

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⁴⁰. Our work is among the first that explore the many-body dynamics of this type of system beyond mean-field theory⁴¹. The system that we describe here is similar to that simulated in ref. 19 using the time-dependent GP equation, and analyzed experimentally in ref. 42. 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. 27). 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. $N < 500$) ⁸⁷Rb condensate reported in ref. 43. The subcritical rotation frequency chosen herein is not fast enough to nucleate quantized vortices in a mean-field condensate. A GP ($M = 1$) simulation with the same parameters only absorbs 0.5ħ of angular momentum per particle (see Fig. 1d).

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

$$i \partial_t |\Psi\rangle = \hat{H} |\Psi\rangle. \qquad (1)$$

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

$$|\Psi\rangle = \sum_{n_1 + \dots + n_M = N} C_{n_1, \dots, n_M}(t) |n_1, \dots, n_M; t\rangle. \qquad (2)$$

That is, the bosons dynamically populate all configurations of N particles in M time-dependent variationally optimized single-particle states, called orbitals (Supplementary Information). The case where $M = 4$ and $N = 100$ is analyzed in the main text of this work whereas other N, M are discussed in the Supplementary Information. In our setting the particle density at rest ($t = 0$) 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 78% of the trapping frequency.

The many-body Hamiltonian in dimensionless units reads:

$$\hat{H}(\mathbf{r}_1, \dots, \mathbf{r}_N; t) = \sum_{i=1}^{N} \hat{h}_i(\mathbf{r}_i; t) + \lambda_0 \sum_{i<j}^{N} \hat{W}(\mathbf{r}_i - \mathbf{r}_j). \qquad (3)$$

We use a normalized gaussian with width $\sigma = 0.25$ as a two-body potential $\hat{W}$ with a two-body interaction strength $\lambda_0(N - 1) = 17.1$. This choice of a short-range interaction is motivated by the similarity of the physics for zero- and short-ranged interaction potentials shown in ref. 44. See Supplemental Material for a direct comparison. The one-body Hamiltonian $\hat{h}_i(\mathbf{r}_i; t)$ is the sum of the kinetic energy $\hat{T}_i = -\frac{1}{2} \partial_{\mathbf{r}_i}^2$ and the trapping potential,

$$\hat{V}(\mathbf{r}, t) = \frac{1}{2}(x(t)^2 + y(t)^2) + \frac{1}{2}\eta(t)(x(t)^2 - y(t)^2), \qquad (4)$$

where

$$\begin{pmatrix} x(t) \\ y(t) \end{pmatrix} = \begin{pmatrix} \cos(\omega t) & \sin(\omega t) \\ -\sin(\omega t) & \cos(\omega t) \end{pmatrix} \begin{pmatrix} x \\ y \end{pmatrix} \qquad (5)$$

and $\omega = 0.78$. The anisotropy parameter varies in time so that it is ramped up from 0 to $\eta_{max} = 0.1$ over $t_r = 80$, held fixed at $\eta_{max}$ for $t_f = 220$, ramped down to zero over $t_r$, and then held fixed at 0 for the rest of a simulation until
<figure><img src="image_3.png" /><figcaption>**Figure 1.** (a) Sketch of the ramping procedure of the anisotropy parameter $\eta(t)$ with plots of the density $\rho(\mathbf{r})$ at representative times in each part. (b) The onset of fragmentation occurs around $t=80$, which corresponds to the time of maximum anisotropy. By the end of the time of maximal anisotropy, the state is completely fragmented, $\rho_1^{(NO)} \approx 40\%$. (c) Orbital angular momenta reach their maximum values early in the period of maximum anisotropy ($t=80-150$), but evolve to their equilibrium values before the anisotropy is removed. Discontinuities in $(L_z)_{ii}$ arise from occupation reordering [see panel (b)]. Fluctuations in $(L_z)_{22}$ and $(L_z)_{33}$ from $t=250$ to $t=500$ are strongly correlated due to vorticity transfer between $\phi_2$ and $\phi_3$. (d) Comparison of energy and angular momentum curves for $M=4$ (solid lines) and $M=1$ (GP, dotted lines). Both quantities agree between the two cases until around $t=80$ when fragmentation becomes significant and the GP ansatz breaks down. Energy and angular momentum oscillate about their maximal values by $t=200$, while the anisotropy is still maximal. The angular momentum stabilizes at $L_z \approx 1.25$ per particle ($M=4$) and $L_z < 0.5$ per particle ($M=1$). The strong correlation between energy and $L_z$ indicates that the perturbation strictly excites angular momentum modes in both simulations. To guide the eye, on the bottom plot is marked with a horizontal line at 1 on the vertical axis. All quantities shown are dimensionless.</figcaption></figure>

$t=500$ as shown in Fig. 1a. The resulting healing length, $\xi = \frac{1}{\sqrt{2\lambda_0(N-1)}} = 0.17$ 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 N-body density$^{10,11}$:

$$ \rho^{(1)}(\mathbf{r}_1|\mathbf{r}'_1; t) = N \int \Psi(\mathbf{r}_1, \mathbf{r}_2, \dots, \mathbf{r}_N; t) \times \Psi^*(\mathbf{r}'_1, \mathbf{r}'_2, \dots, \mathbf{r}'_N; t) d\mathbf{r}_2 \dots d\mathbf{r}_N, \quad (6) $$

with $\Psi$ normalized to 1. The trace operation can be understood as eliminating knowledge of the many-body effects to obtain a single-particle operator, the RDM. The diagonal, $\mathbf{r} = \mathbf{r}'$, of the RDM is the particle density, $\rho(\mathbf{r}; t) = \rho^{(1)}(\mathbf{r}|\mathbf{r}; t)$, or simply the density. The RDM is written in its eigenbasis as

$$ \rho^{(1)}(\mathbf{r}|\mathbf{r}'; t) = \sum_{i=1}^{M} \rho_i^{(NO)}(t) \phi_i^*(\mathbf{r}'; t) \phi_i(\mathbf{r}; t) \quad (7) $$

The eigenvalues, $\rho_i^{(NO)}(t)$, are ordered in decreasing magnitude and called natural occupations and the eigenfunctions, $\phi_i(\mathbf{r}; t)$, are called the natural orbitals or synonymously, fragments. Note that the density, $\rho(\mathbf{r}; t) = \rho^{(1)}(\mathbf{r}|\mathbf{r}; t)$, is equal to the sum of the squared amplitudes of the natural orbitals, weighted with their natural occupations.

This treatment bears some resemblance to the two-fluid model of Helium below its lambda point, where the total density has a superfluid component and thermal component interacting through a mutual viscosity. The MCTDHB treatment can likewise be thought of as an M-fluid model, where there are M interpenetrating
<figure><img src="image_4.png" /><figcaption>**Figure 2.** (a) The density $\rho(\mathbf{r})$ shows some density minima due to phantom vortices in $\phi_1$, but no true density node. (b–e) The natural orbital densities, $|\phi_i(\mathbf{r})|^2$ ($i=1, 2, 3, 4$), are plotted. Many phantom vortices are present in each orbital. (f) The phase, $S_g(\mathbf{r}|0)$, of $g^{(1)}(\mathbf{r}|0)$ is plotted. Note, $S_g$ is not the many-body phase, $S_{MB}(\mathbf{r}_1, \dots, \mathbf{r}_N; t)$, which is too complicated to visualize. Since $\phi_1$ carries most of the particles at this time, $S_g$ bears strong resemblance to $S_1$. (g–j) The phases, $S_\phi$ of the natural orbitals $\phi_i$ ($i=1, 2, 3, 4$) are plotted. Each $2\pi$ phase discontinuity marks a phantom vortex core. The central phantom vortex in $\phi_2$ (c,d) is mutated from its initial angular node, and persisted for the duration of simulation. The two central phantom vortices in $\phi_3$ (d,i) were mutated from its initial angular nodes and are transient. The three centermost phantom vortices in $\phi_4$ (e,j) were mutated from an “I” shaped node and persisted for the length of the simulation. All other phantom vortices nucleated at the edge of their orbital density and are transient. All panels are plotted at $t=115.0$, when $\rho_1^{(NO)} = 82.0\%$, $\rho_2^{(NO)} = 11.8\%$, $\rho_3^{(NO)} = 4.1\%$, and $\rho_4^{(NO)} = 2.1\%$. See complementary Video S1 (Supplementary Information). All quantities shown are dimensionless.</figcaption></figure>

superfluid components coupled through the MCTDHB equations of motion [39], but no thermal component. In this treatment, there are hence $\$M\$ order parameters, one for each fragment.

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 $M=1$, 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 and fragmentation is through the first-order correlation function, $g^{(1)}(\mathbf{r}|\mathbf{r}'; t)$, defined in terms of the RDM as

$$g^{(1)}(\mathbf{r}|\mathbf{r}'; t) = \frac{\rho^{(1)}(\mathbf{r}|\mathbf{r}'; t)}{\sqrt{\rho^{(1)}(\mathbf{r}|\mathbf{r}; t)\rho^{(1)}(\mathbf{r}'|\mathbf{r}'; t)}} \quad (8)$$

A coherent state at time $\tau$ must satisfy $|g^{(1)}(\mathbf{r}|\mathbf{r}'; \tau)|=1$, while any other value indicates fragmentation$^{10}$. Intuitively, $|g^{(1)}|$ measures how well the many-body wavefunction is described by a hyphenate state.

The off-diagonal, $\mathbf{r} \neq \mathbf{r}'$, terms in $g^{(1)}$ and the orbitals $\phi_i$ are complex-valued. Thus it is instructive to plot their magnitudes and phases separately. We hence define the phase of $g^{(1)}(\mathbf{r}|\mathbf{r}'; t)$ to be

$$S_g(\mathbf{r}|\mathbf{r}'; t) = \arg[g^{(1)}(\mathbf{r}|\mathbf{r}'; t)], \quad (9)$$

and the orbital phases to be

$$S_i(\mathbf{r}; t) = \arg[\phi_i(\mathbf{r}; t)]. \quad (10)$$

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

$$S_{MB}(\mathbf{r}_1, \dots, \mathbf{r}_N; t) = \arg[\Psi(\mathbf{r}_1, \dots, \mathbf{r}_N; t)]. \quad (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 >$1h$ (Fig. 1). We name the vortices which
<figure><img src="image_5.png" /><figcaption>**Figure 3.** (a) The density $\rho(\mathbf{r})$ shows a density minimum at the origin which is 75% the maximum density, but no true node. (b–e) The natural orbital densities, $|\phi_1(\mathbf{r})|^2$ ($i$ = 1, 2, 3, 4), are plotted. (f) The phase, $S_g(\mathbf{r}|0)$, of $g^{(1)}(\mathbf{r}|0)$ is plotted. Note, $S_g$ is not the many-body phase, $S_{MB}(\mathbf{r}_1, \dots, \mathbf{r}_N; t)$, which is too complicated to visualize. There are singularities in $S_g$ that are associated with phantom vortices. Since $\phi_1$ and $\phi_4$ have a node at $\mathbf{r}' = 0$, they do not contribute to $S_g$. (g–j) The phases, $S_p$, of the natural orbitals $\phi_i$ ($i$ = 1, 2, 3, 4) are plotted. Each $2\pi$ phase discontinuity marks a phantom vortex core. The depicted phantom vortex configuration is stable for more than 100 trap periods. The phantom vortices in $\phi_1$ (b,g) and $\phi_4$ (e,j) nucleated along preexisting topological defects through node mutation whereas the pair of phantom vortices in $\phi_3$ (d,i) was nucleated by slow orbital-orbital vortex transfer as ghost vortices which then moved to the center. The nucleation of the phantom vortex pair was coupled to the destruction of a pair in $\phi_2$. (c,h) See the coupled oscillation of $(L_z)_{22}$ and $(L_z)_{33}$ [(Fig. 1c) and Video S1 (Supplementary Information)]. All panels are plotted at $t = 450.0$, when $\rho_1^{(NO)} = 40.8\%$, $\rho_2^{(NO)} = 25.0\%$, $\rho_3^{(NO)} = 20.8\%$, and $\rho_4^{(NO)} = 13.5\%$. See complementary Video S1 (Supplementary Information). All quantities shown are dimensionless.</figcaption></figure>

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 $N=100$ particles with $M=4$ orbitals starting in the ground state at $t=0$ to $t=500$ (see also Video S1 in Supplementary Information). At $t=0$, $\phi_1, \dots, \phi_4$ resem- bled 1 s, $2p_x$, $2p_y$, and $2s$ orbitals, respectively. The system began with $\rho_1 > 99.7\%$ occupation in $\phi_1$ and thus was almost entirely condensed (Fig. 1b). As the anisotropy was ramped up, $\phi_1, \dots, \phi_4$ deformed smoothly and rotated with the potential. By $t=29.6$, due to a reordering of the natural occupations, $\phi_3$ and $\phi_4$ switched labels. At this time, $\phi_1, \phi_3$, and $\phi_4$ already showed faint hints of vortices at the edge of their density similar to those seen at later times in and Video S1. These vortices are the fragmented counterparts of *ghost vortices* reported in ref. 21. 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²¹ 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 $t=45$, many ghost vortices were established in $\phi_1$ and $\phi_4$. Since the state was still more than 99% condensed, the ghost vortices in $\phi_1$ manifested in the density and would thus have been detectable in high-fidelity density measurement. A one-dimensional cut along a core of a ghost vortex in $\rho(\mathbf{r})$ showed that the outer density maximum was less than 1% of the density maximum at the center of the cloud. By this time, the two initial lobes in $\phi_2$ 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 and Video S1. 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 $t=70$, ghost vortices in $\phi_4$ fused with a distorted angular node to split the two lobes into four with an “I” shaped node (see Video S1 at $t=99.4$). This node mutated to nucleate three persistent phantom vortices which merged into the charge-3 vortex shown in Fig. 3e,j.

By $t=90$, the system was significantly fragmented, $\rho_1^{(NO)} \approx 90\%$ (Fig. 1b). Fragmentation obscured the vortex-induced density nodes of $\phi_1$ from being visible in $\rho(\mathbf{r})$ because density from other orbitals filled the vortex cores present in $\phi_1$. We comment that although larger particle numbers delay the onset of fragmentation, we have observed phantom vortices in simulations with up to $N=10^4$ (Supplementary Information). From $t=90$ to $t=140$, each orbital had a complicated vortex structure while the density $\rho(\mathbf{r})$ was largely featureless (Figs 1a and 2a). 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 and Video S1. In this time interval, the orbital
<figure><img src="image_6.png" /><figcaption>**Figure 4.** To visualize the four-dimensional single-particle correlation function, $g^{(1)}(\mathbf{r}|\mathbf{r}', t=219.4)$, we fix a reference point at $\mathbf{r}' = (0, 0)$ (a,c) and $\mathbf{r}' = (1.25, 0)$ (b,d). The function is complex, so we plot the magnitudes (a,b) and phases (c,d) separately. In all panels, we mark the reference point, $\mathbf{r}'$, with a white circle. In panels (a) and (c) [(b) and (d)], the $\mathbf{r}'$ is colocated with a phantom vortex core in $\phi_2$ [$\phi_1$]. In both cases, phantom vortex cores in all orbitals appear as spots of almost complete incoherence, $|g^{(1)}| \approx 0$, while the core colocated with the reference point has a $|g^{(1)}| \approx 1$ (full coherence). $\rho_1^{(\text{NO})} = 35.3\%$, $\rho_2^{(\text{NO})} = 35.2\%$, $\rho_3^{(\text{NO})} = 17.5\%$, and $\rho_4^{(\text{NO})} = 12.0\%$. See complementary Videos S2 and S3 for visualizations of $g^{(1)}$ at $t=115$ and $t=450$ with different reference points $\mathbf{r}'$ (Supplementary Information). All quantities shown are dimensionless.</figcaption></figure>

angular momenta, $(L_z)_{ii} \equiv \langle \phi_i | L_z | \phi_i \rangle$, attained their maximal values and started to decay to their equilibrium values despite the maximal anisotropy.

After $t=150$, the trap was still maximally anisotropic, yet the system energy and total angular momentum became saturated (Fig. 1d). By $t=144$, the vortex structures of $\phi_2$ and $\phi_4$ reached their steady state. $\phi_2$ 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 $t=150$, $\phi_4$ 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 $t=380$ when the trap became symmetric again (Fig. 3e,j). Near $t=220$, $\phi_2$ and $\phi_1$ swapped labels due to occupation reordering (see $\rho_1^{(\text{NO})}$ and $\rho_2^{(\text{NO})}$ in Fig. 1b), therefore the steady state of the phantom vortex structure can be seen in Fig. 3b,e,g,j. From $t=220$ on, the labeling of the orbitals remained fixed because the occupation numbers did not change order anymore (cf. Fig. 1b).

At $t=220$, there were two corotating vortices near the center of the cloud in $\phi_2$ and none present in $\phi_3$. This phantom vortex pair nucleated at the edge of the cloud and then moved towards the center, resembling mean-field vortex nucleation$^{19,21}$. However, in our treatment, this pair was transient. The phantom vortices then transferred from $\phi_2$ to $\phi_3$, marking a second mechanism unique to phantom vortex nucleation: slow orbital-orbital vorticity transfer, which we now describe in detail. By $t=275$, there were two prominent ghost vortices in $\phi_3$ and the phantom vortex pair in $\phi_2$ had returned from the center to the edge of this fragment. Gradually, the vortices in $\phi_2$ disappeared by exiting the edge of the orbital while the prominent phantom vortex pair in $\phi_3$ entered the orbital bulk. By $t=340$, the phantom vortex nucleation in $\phi_3$ was complete, marked by the intervortex separation being as small as previously in $\phi_2$ at $t=220$, before the transfer (Video S1). This direct interaction between $\phi_2$ and $\phi_3$
appears as a strong correlation in the orbital angular momentum, ($L_z$)$_{22}$ and ($L_z$)$_{33}$, as seen from $t=250$ to $t=500$ in Fig. 1c. The transfer of vorticity between fragments occurred on a time scale much slower ($\tau \approx 50–100$) than the trap rotation period ($\tau \approx 8$).

As the anisotropy was ramped down, the phantom vortex pair in $\phi_3$ approached the trap center and reached a minimal separation of $\approx 1.5$ by $t=380$ 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 $\phi_4$ coalesced into a dynamically stable charge-3 phantom vortex. This contrasts with the known spontaneous decay of charge-n vortices to n single vortices$^{45}$ 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, $|g^{(1)}(\mathbf{r}|\mathbf{r}'; t)|$ (Fig. 4). We have two key observations about the coherence of phantom vortices. First, we fix $\mathbf{r}'$ near the core of a phantom vortex; remarkably, $|g^{(1)}(\mathbf{r}|\mathbf{r}'; t)|$ is close to 0 for all $\mathbf{r}$ near and inside the cores of other phantom vortices. Second, we fix $\mathbf{r}'$ away from the cores of all phantom vortices and again observe that $|g^{(1)}(\mathbf{r}|\mathbf{r}'; t)|$ is close to 0 for all $\mathbf{r}$ 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 $g^{(1)}(\mathbf{r}|\mathbf{r}'; t)$ is measurable via interference experiments, e.g. ref. 46.

## 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 $\rho(\mathbf{r}; t)$. Vortices, as seen in the laboratory (see for instance refs 16,47,48), 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 possesses 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 $\phi_1$ and $\phi_4$ 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 $\rho(\mathbf{r}; t)$. Perhaps some experimentally detected vortices are in fact phantom vortices rather than mean-field vortices, as also suggested in ref. 49. 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.

## 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* (2008).

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. Feder, D. L., Svidzinsky, A. A., Fetter, A. L. & Clark, C. W. 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. & 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. Lundh, E., Martikainen, J. P. & Suominen, K. A. Vortex nucleation in Bose-Einstein condensates in time-dependent traps. *Phys. Rev. A* **67**, 063604 (2003).

20. Parker, N. G. & Adams, C. S. Emergence and Decay of Turbulence in Stirred Atomic Bose-Einstein Condensates. *Phys. Rev. Lett.* **95**, 145301 (2005).

21. Tsubota, M., Kasamatsu, K. & Ueda, M. Vortex lattice formation in a rotating Bose-Einstein condensate. *Phys. Rev. A* **65**, 023603 (2002).

22. Lobo, C., Sinatra, A. & Castin, Y. Vortex Lattice Formation in Bose-Einstein Condensates. *Phys. Rev. Lett.* **92**, 020403 (2004).

23. 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).

24. Cederbaum, L. S. & Streltsov, A. I. Best mean-field for condensates. *Phys. Lett. A* **318**, 564 (2003).

25. 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).

26. Mueller, E. J., Ho, T. L., Ueda, M. & Baym, G. Fragmentation of Bose-Einstein condensates. *Phys. Rev. A* **74**, 033612 (2006).

27. 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).

28. Bader, P. & Fischer, U. R. Fragmented many-body ground states for scalar bosons in a single trap. *Phys. Rev. Lett.* **103**, 060402 (2009).

29. 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).

30. Zhou, Q. & Cui, X. Fate of a Bose-Einstein condensate in the presence of spin-orbit coupling. *Phys. Rev. Lett.* **110**, 140407 (2013).

31. 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).

32. Lode, A. U. J. Multiconfigurational time-dependent Hartree method for bosons with internal degrees of freedom: Theory and composite fragmentation of multicomponent Bose-Einstein condensates. *Phys. Rev. A* **89**, 063601 (2016).

33. Tsatsos, M. C. Attractive Bose gas in two dimensions: An analytical study of its fragmentation and collapse. *Phys. Rev. A* **93**, 043604 (2014).

34. 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).

35. 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).

36. 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).

37. 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).

38. 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).

39. Lode, A. U. J. * Tunneling Dynamics in Open Ultracold Bosonic Systems* (Springer, 2015).

40. 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).

41. Tsatsos, M. C. & Lode, A. U. J. Resonances and Dynamical Fragmentation in a Stirred Bose-Einstein Condensate. *J. L. Temp. Phys.* **181**, 171 (2015).

42. Madison, K. W., Chevy, F., Wohlleben, W. & Dalibard, J. Vortex formation in a stirred Bose-Einstein condensate. *Phys. Rev. Lett.* **84**, 806 (2000).

43. Schmied, R. et al. Bell correlations in a Bose-Einstein condensate. *Science* **352**, 441 (2016).

44. 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).

45. Shin, Y. et al. Dynamical Instability of a Doubly Quantized Vortex in a Bose-Einstein Condensate. *Phys. Rev. Lett.* **93**, 160406 (2004).

46. Hofferberth, S., Lesanovsky, I., Fischer, B., Schumm, T. & Schmiedmayer, J. Non-equilibrium coherence dynamics in one-dimensional Bose gases. *Nature* **449**, 324 (2007).

47. Fetter, A. L. Rotating trapped Bose-Einstein condensates. *Rev. Mod. Phys.* **81**, 647 (2009).

48. 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).

49. Sakmann, K. & Kasevich, M. Single-shot simulations of dynamic quantum many-body systems. *Nature Physics* **12**, 451–454 (2016).

## 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.

## Author Contributions

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

Supplementary information accompanies this paper at http://www.nature.com/srep

Competing financial interests: The authors declare no competing financial interests.

How to cite this article: Weiner, S. E. et al. Phantom vortices: hidden angular momentum in ultracold dilute Bose-Einstein condensates. *Sci. Rep.* **7**, 40122; doi: 10.1038/srep40122 (2017).

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
<figure><img src="image_7.png" /></figure>

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

© The Author(s) 2017