Fragmentation and correlations in a rotating Bose–Einstein condensate undergoing breakup

The theoretical investigation of rotating Bose–Einstein condensates has mainly focused on the emergence of quantum vortex states and the condensed properties of such systems. In the present work, we concentrate on other facets by examining the impact of rotation on the ground state of weakly interacting bosons confined in anharmonic potentials computed both at the mean-field level and particularly at the many-body level of theory. For the many-body computations, we employ the well-established many-body method known as the multiconfigurational time-dependent Hartree method for bosons. We present how various degrees of fragmentation can be generated following the breakup of the ground state densities in anharmonic traps without ramping up a potential barrier for strong rotations. The breakup of the densities is found to be associated with the acquisition of angular momentum in the condensate due to the rotation. In addition to fragmentation, the presence of many-body correlations is examined by computing the variances of the many-particle position and momentum operators. For strong rotations, the many-body variances become smaller than their mean-field counterparts, and one even finds a scenario with opposite anisotropies of the mean-field and many-body variances. Further, it is observed that for higher discrete symmetric systems of order k, namely three-fold and four-fold symmetry, breakup to k sub-clouds and emergence of k-fold fragmentation take place. All in all, we provide a thorough many-body investigation of how and which correlations build up when a trapped Bose–Einstein condensate breaks up under rotation.

generated following the breakup of the ground state densities in anharmonic traps without ramping up a potential barrier for strong rotations.The breakup of the densities is found to be associated with the acquisition of angular momentum in the condensate due to the rotation.In addition to fragmentation, the presence of many-body correlations is examined by computing the variances of the many-particle position and momentum operators.For strong rotations, the many-body variances become smaller than their mean-field counterparts, and one even finds a scenario with opposite anisotropies of the mean-field and many-body variances.Further, it is observed that for higher discrete symmetric systems of order k, namely three-fold and four-fold symmetry, breakup to k sub-clouds and emergence of k-fold fragmentation take place.All in all, we provide a thorough many-body investigation of how and which correlations build up when a trapped Bose-Einstein condensate breaks up under rotation.

I. INTRODUCTION
The successful experimental realization of a rotating Bose-Einstein condensate (BEC) has paved the way to explore various rich physics of correlated quantum systems [1][2][3][4][5].Butts and Rokhsar [6] first evaluated the wave function of a rotating BEC using the lowest Landau level approximation with the help of a Gross-Pitaevskii functional.The rotating ultracold bosonic gases have led to the investigation of the occurrence of quantized vortices [1,3,6], vortex nucleation [7], emergence of quantum fluctuations [8,9], and presence of the fractional quantum Hall effect in weakly interacting quantum systems [10].
There have been studies where breaking up of fast rotating objects is widely detected in systems extending from astronomical objects, e.g., galaxies and supermassive rotating stars [11] to the quantum systems in nuclear physics [12].In atomic physics, the emergence of superfluid flow of a rotating quantum gas has been explored experimentally in Ref. [13], in an anharmonic potential.Theoretically, the emergence of breakup in a rotating dipolar condensate was investigated in three-dimension in Ref. [14] and finally, in a rotating pancakelike asymmetric quartic-quadratic potential in Ref. [15].A substantial volume of literature exists corresponding to the rotating BEC, investigating various rich quantum features by employing Gross-Pitaevskii mean-field approximation [16][17][18].However, there exists much less studies of rotating BEC in the many-body domain, see [19][20][21][22][23][24], that explore interesting many-body quantum features like the fragmentation of the condensate and correlations.
In the regime of "ultrafast rotation", when the rotation frequency comes closer to the trapping frequency in case of a harmonic trapping potential, the system would tend to escape as the centrifugal force would cancel the trapping force.This problem can be resolved either by introducing an anharmonic term into the confining potential or by adding anisotropy in the harmonic potential.Anisotropy makes the confining potential elongated.Ref. [49] investigated the appearance of vortex rows of BEC in an anisotropic potential under rotation using Gross-Pitaevskii formalism.Motivated also by these findings, in this work, we study a system of bosons subject to rotation confined in various anharmonic trap geometries with discrete rotational symmetry in two-dimension 2D, namely an elongated trap, a three-fold symmetric trap and finally four-fold symmetric trap.
There exist several ideas to induce mechanical rotation and coupling of the internal states in the condensate, such as, phase imprinting (using an electromagnetic field) as proposed by Williams and Holland in Ref. [50] and optical spoon stirring method where quantized vortices are observed in a stirred gaseous condensate of atomic rubidium [1].Here, we study the physics of a 2D weakly interacting ultracold bosons confined in different anharmonic potentials both at the mean-field level and particularly at the many-body level of theory.Specifically, we investigate how the system transformed from a fully-condensed state into a fragmented state followed by breaking up of the condensate density induced by the rotational motion.As variance is a sensitive probe that characterizes many-body correlations even for fully condensed systems [51,52], hence we also analyze the variances of many-particle operators of the fragmented state in the rotating frame.The qualitative difference between the mean-field and many-body variances is a useful tool to explain the nature of many-body correlations.The many-particle position variance characterizes to what extent a wave-packet spreads or narrows down, similarly the momentum variance is associated with the size of the wave-packet in momentum space.Hence, we emphasize on the emergence of correlations by investigating variances of many-particle operators like the position, momentum and angular momentum.Interestingly, the fluctuations present in the system are hardly observed in the angular momentum variance.Thus, we only present the position and momentum variances in the main text and give the detail discussion of the angular momentum variance in the supplemental material.The time-dependent Gross-Pitaevskii mean-field theory [53,54] is the most celebrated theoretical model to investigate the many-particle systems of ultracold bosonic atoms.However, this method is unable to study fragmentation and correlations owing to its building via the mean-field ansatz.In this paper, we employ a well established many-body numerical method named the multiconfigurational time-dependent Hartree method for bosons (MCTDHB) [55,56] to accurately solve the Schrödinger equation at the many-body level for ultracold atoms subject to a rotation.The MCTDHB method is a bosonic version of the MCTDH family of methods [57][58][59][60][61][62][63][64][65][66][67][68][69][70][71][72] which is able to self-consistently describe the physics involving the presence of many-body correlations.The main focus of the applications of MCTDHB has been the emergence of fragmentation of the condensate, where the one-body reduced density matrix has multiple significant eigenvalues.For numerical simulation of the results presented in this work, we use the MCTDH-X software [72][73][74][75].Finally, the supplemental material reports the benchmarking of MCTDHB for an exactly solvable many-body model under rotation and also presents the convergence of the many-body results of our present work.

II. SETUP AND THEORETICAL TOOLS
We consider a system of weakly interacting bosonic atoms in two spatial dimensions 2D confined in non-spherically symmetric trapping potentials in the rotating frame.The properties of these trapped bosons can be described by the (time-dependent) many-body Schrödinger equation.The Schrödinger equation dealing with a many-boson system is usually solved by employing the mean-field Gross-Pitaevskii approximation.However, the reduced density matrix involved in the Gross-Pitaevskii approximation has only a single eigenvalue and it involves a single basis state and thereby is unable to capture the many-body features such as fragmentation and correlations.
In MCTDHB, the (time-dependent) optimized one-body basis is used.Here the basis set and the expansion coefficients in the basis are optimized variationally [55,56].The MCTDHB is a numerically exact method [76] and can describe both coherent and fragmented condensates.MCTDHB includes the theory of Gross-Pitaevskii approximation as a special case when only a single one-body state is considered.

A. Hamiltonian
The general Hamiltonian of N interacting bosons is given as where the single-particle Hamiltonian ĥ is composed of the kinetic energy and the external potential energy, respectively.Here, the interaction of ultracold dilute bosonic gases is considered to be a finite range interaction and modelled by a Gaussian function [63,64,77], with σ = 0.25.
This avoids the regularization of the delta contact potential in 2D.The interaction strength λ 0 is scaled with the number of bosons N as Λ = λ 0 (N − 1), where Λ is the interaction parameter.One uses the interaction parameter to define the mean-field regime.In our study, we work in the units = m = 1 and all the quantities are dimensionless.We also consider three different trapping potentials V (r) that we shall discuss in the next section.The first setup is the elongated trap − that leads to breaking up of the ground state density into two clouds.After that we move to more complex traps, namely three-fold symmetric and four-fold symmetric traps to investigate and establish the generality of the results.Hence, our strategy is to first study an elongated trap, and then a three-fold symmetric trap, and finally, a four-fold symmetric trap to see what stays between the two to three and the three to four-fold symmetric traps.
In the rotating frame, the kinetic energy operator is modified and can be written as, here ω r is the rotation frequency and lz = xp y − ŷ px is the angular-momentum operator.
An alternative way to mimic the rotational effect in the condensate is by introducing a synthetic gauge field A(r) as Consider the following general form of the gauge field: Then, expansion of Equation ( 4) leads to, For the specific case b = −a, Equation (6) becomes Combining Equations ( 2), ( 3) and (7) we have ĥ(r) = T (r) + V (r), (8) where the modified confining potential is V (r) = V (r) + 1 2 a 2 r 2 and a = ω r corresponds to the rotation frequency of the condensate.

B. Many-body method
The MCTDHB method employs time-adaptive orbitals to represent the field operator as a sum of the M time-dependent single-particle states The ansatz of the MCTDHB wavefunction is The summation in Equation ( 10) runs over all N +M −1 N possible time-dependent configurations n = (n 1 , ..., n M ) with fixed particle number N = M i=1 n i .To derive the MCT-DHB equations, the time-dependent variational principle [78][79][80] is employed for the ansatz in Equation (10).Thus, resulting in two-coupled equations of motion -a set of linear equations for the coefficients {C n } and a set of non-linear equations for the orbitals {φ j (r); j = 1, ..., M }, see Refs.[55,56] for details and derivation of the equations of motion.
In the following work, the self-consistent ground state is achieved by relaxing the system via imaginary-time propagation and is hence determined by the variational principle.Thus, in the following sections, we omit the time-dependency from the various quantities and observables that are involved in the many-body simulations.

C. Quantities of interest
In this section, we define the quantities of interest, namely the one-body density, the eigenvalues of the one-body reduced density matrix (RDM), the expectation value of the angular momentum operator, and finally, the many-particle variances of the position, momentum, and angular momentum operators.

One-body reduced density matrix (RDM), one-body density, and natural occupations
The one-body RDM of the N-boson state |Ψ is a hermitian matrix and is defined as in its eigenbasis {φ q (r)}.The matrix elements ρ kq = Ψ| b † k bq |Ψ represent the one-body RDM using M orbitals corresponding to the creation (annihilation) operators b † k ( bq ).The diagonal of ρ (1) (r, r ) is referred to as the one-body density ρ(r) which is ρ(r) = ρ (1) (r, r = r).
The eigenvalues of the one-body RDM are obtained by the diagonalization of Equation (11) which corresponds to a unitary transformation of the orbitals φ q (r) to the natural orbitals φ Here, the eigenvalues n j are normalized as M j=1 n j = 1 and, without loss of generality, they are sorted in magnitude such that n 1 ≥ n 2 ≥ ... throughout this work.The eigenvalues n j are termed natural occupations and characterize the degree of condensation and fragmentation of the bosons.Thus, the system with one-body RDM consisting only a single macroscopicallycontributing eigenvalue n 1 is said to be condensed [31].When the one-body RDM has k macroscopically-occupied eigenvalues, the system is referred to as k-fold fragmented [34].

Angular momentum
For a 2D many-particle systems, there is only a single component of the angular momentum operator, i.e., LZ = Bosonic systems with angular momentum provide rich quantum features beyond mean-field theory, such as, phantom vortices [81], spatially partitioned vortices [51], and fragmentation [51,81,82].In the following studies, we investigate the expectation value of the angular momentum operator per particle 1 N Ψ| LZ |Ψ , for three different confining traps to see the effect of the rotation.In our system, we expect an intricate dynamics of the angular momentum acquisition under rotation and build up of correlations.

Many-particle variances
The variance of a many-particle observable Ô = N j=1 ô(r j ) can be written as [83] 1 Here, the expectation value of Ô = N j=1 ô(r j ) is dependent only on the one-body operator, whereas the expectation of Ô2 is a combination of one-and two-body operators Ô2 = N i=1 ô2 (r j ) + j<k 2ô(r j )ô(r k ).ρ jklm are the two-particle reduced density matrix elements, In the following work, to analyze the emergence of many-body correlations in the rotating condensate, we study the many-particle variances per particle of the position, momentum, and angular momentum operators [52].

III. MANY-BODY PHYSICS OF BOSE-EINSTEIN CONDENSATES BREAKING UP UNDER ROTATION: RESULTS AND DISCUSSION
We investigate the impact of rotation on the ground state of weakly interacting bosonic atoms by dividing the analysis into two main parts depending on the confining anharmonic potentials.First, the breakup and fragmentation processes in an elongated trap are investigated and then the breakup and fragmentation in more complicated traps of discrete spatial symmetry are explored.In the following sections, we analyze the effect of rotation on various static properties, namely the ground state energy, one-body density, natural occupations, expectation value of the angular momentum operator, and finally, the variances of the many-particle position, momentum, and angular momentum operators as a function of the rotation frequency ω r .For the numerical computations, the MCTDH-X implementation of the MCTDHB theory is employed [72][73][74][75].In our work, we consider N = 10 weakly interacting ultracold bosonic atoms interacting via a Gaussian interaction [63,64,77] with the interaction parameter chosen as Λ = 0.1 throughout the computations.The grid used to represent the Hamiltonian Equation ( 18

A. Breaking up to two clouds and pathway to two-fold fragmentation
Let us discuss the impact of rotation on the weakly interacting bosons trapped in an elongated confinement.The elongated trap is an anharmonic potential that elongates the condensate in the x-direction and can be written as For this trap, we consider M = 4 self-consistent orbitals to investigate the ground state properties of the condensate.We recomputed these results with M = 8 self-consistent orbitals to check the convergence of the system, see the elaborate discussion in the supplemental material.The rotation frequency range is ω r = [0, 2.0] in this trap.to 10 −3 for larger rotation frequencies ω r ≥ 1.6.The presence of these structures at the intermediate rotation might suggest that something interesting is happening at the manybody level.Therefore, we dig deeper than the energy of the system to see the many-body features.Now let us discuss the behaviour of the ground-state densities per particle ρ(r) N of the rotating condensate confined in the elongated trap.The densities computed at the manybody level of theory are shown for three different rotation frequencies ω r in Figure 2. In absence of rotation, the density displays a single cloud where all bosons accumulate in the center of the trap.This behaviour persists with the inclusion of rotation for slow rotation frequencies, e.g., at ω r = 0.5.Further increase in rotation induces a breaking of the density into two clouds.The breakup of the ground-state density confined in the potential given by Equation ( 33) corresponds to the scenario where the minimum of the potential is splitted into two parts and shifted with rotation by creating an effective double-well potential.It is also observed that the distance between the two densities increases with increase of ω r [Figure 2(b)-(c)].We have also computed the density per particle at the mean-field level and identical features are observed in the density profile.Hence, it is observed that the densities computed at the mean-field and many-body levels in the real space show identical pattern.
In the density profile, we show three specific frequencies that correspond to slow, fast, and faster rotations.
To understand whether the splitting of the density into two clouds [as shown in Figure 2(c)] for faster rotations leads to fragmentation of the condensate, we further discuss the dependence of the natural occupations n j N .In the elongated trap, four natural occupations corresponding to the four natural orbitals are employed and found to vary with the rotation frequency ω r (see Figure 3).The system remains fully condensed, i.e.,   As we know, the variance is a sensitive probe of correlations that allows one to study the quantum fluctuations present in a system [51,52].Thus, it would be interesting to investigate the variance of many-particle operators which signifies many-body correlations for the fragmented condensate in the rotating frame.Therefore, we further analyze the impact of rotation on the behaviour of the many-particle variances of position and momentum operators which are sensitive to rotation.
Figures 5(a)-(b) display the behaviour of the many-particle position variance per particle Ŷ along the x and y-directions respectively, as a function of the rotation frequency ω r in the elongated trap.The mean-field and many-body position variances along the x- , coincide till about ω r = 0.9, indicating the essential absence of correlations in the system [Figure 5 (a)].Further, the mean-field position variance monotonically increases with ω r which signifies spreading of the density as observed in Figure 2. Also, the manybody position variance slowly increases till ω r = 1.2.But, now, we observe a deep at ω r = 1.3 that corresponds to the loss of coherence in the natural occupation, see Figure 3.
After increasing and reaching a maximum at rotation frequency of about ω r ∼ 1.4, 1 N ∆ 2 X starts to decrease which incorporates the emergence of a small amount of depletion of the condensate [see Figure 5(a)].For a faster rotation, at ω r = 2, the position variance 1 N ∆ 2 X decreases significantly which goes hand in hand with the emergence of fragmentation of the condensate.Similar behaviour of the position variance along the x-direction is observed in a two-dimensional double well in [84], albeit without the resonant-like behaviour described above.The position variance per particle along the y-direction 1 N ∆ 2 Ŷ is almost frozen and varies slowly as shown in Figure 5(b).This might suggest that excitations along the tighter y-directions are practically not involved, at least as far as the position variance is considered.
The 1 N ∆ 2 Ŷ computed at the mean-field level matches that at the many-body level almost for all the rotation frequencies except of small difference for the intermediate ω r .As above, convergence of the results is detailed in the supplemental material.Further, we can describe the anisotropy of the variance by considering two facts, (1) by comparing two quantities, one along the x-direction and the other along the y-direction, whether they are similar or different.
(2) by comparing these quantities at the mean-field and many-body levels which demonstrates whether the anisotropy of the many-particle variances are alike or opposite.Hence, in the elongated trap and for small rotations, it is observed that Thus, the anisotropy of the many-particle position variance computed at the mean-field and many-body levels are alike at slow rotation.However, for fast rotation we find Hence, it indicates that the many-particle position variances display opposite anisotropy when computed at the mean-field and many-body levels.Now, we investigate the behaviour of the many-particle momentum variance per particle Inverse to the position variance, in the case of the momentum variance for slow rotation we find Hence, the anistropy of the momentum variances are alike when computed at the many-body and mean-field levels of theory for slow rotation.For a faster rotation however, we find This signifies the presence of opposite anisotropy of the momentum variance.This feature was not found before for static double well [84].Hence, this indicates one of the distinct features of the rotation, which makes both the many-particle position and momentum variances show opposite anisotropies with respect to the mean-field and many-body levels in the elongated trap.The rotation localized the position variance in the long direction and the momentum variance in the narrow direction and this is purely a many-body effect.
Let us briefly summarize the results so far.The ground state of a Bose-Einstein condensate in a two-dimensional elongated trap is analyzed in presence of rotation.Here, the ground state density splits into two clouds for fast rotation.It is fascinating to observe the effect of rotation on the condensate as it leads to a transition from condensed state to a fully two-fold fragmented state in a single well without ramping up a barrier.Thus, rotation can be used as a probe to manipulate various degrees of fragmentation which we will further discuss in the next section where higher-order fragmentations are explored.Further, we observe the presence of an opposite anisotropy both in the position and momentum variances by comparing the mean-field and many-body results for fast rotations.

B. Breaking up to several clouds and pathway to higher-order fragmentation
So far in the elongated trap, we observe breaking up of the ground state density into two clouds, followed by the emergence of two-fold fragmentation and presence of opposite anisotropy both in the position and momentum variances for faster rotations.Now, it would be fascinating to explore some more complex potentials with k-fold rotational symmetry, explicitly three-fold and four-fold symmetric potentials, to analyze the impact of rotation on generic properties.Would the density break into more than two clouds?Is higher-order fragmentation possible?How would the angular momentum enter into the condensate?
Finally, when correlations set in, who wins, the many-body or the mean-field variance?
Both for the position and momentum variances?
1. Bosons in a three-fold symmetric trap under rotation In this section, we analyze the impact of rotation on weakly interacting bosons confined in a three-fold symmetric trap of the form In this trap, we consider M = 3 self-consistent orbitals to compute the ground state properties of the system.We also computed the following results with M = 6 self-consistent orbitals for convergence and checked the consistency of our results (see the supplemental material).The rotation frequency range is ω r = [0, 2.5] for the following analysis.We computed the ground-state energy E/N , and found it to display similar pattern as that for the elongated trap.That is for slow rotation, E/N remains almost constant and then E/N drops gradually with increase in ω r .In addition, we found that the mean-field and many-body energies are practically identical for all ω r .The results are shown in the supplemental material, see Figure S1(b).
Figure (6) displays the ground-state densities per particle of a rotating condensate for three different rotation frequencies ω r , confined in the three-fold symmetric trap [Equation (34)] computed at the many-body level.In the absence of rotation, the one-body density depicts a single cloud where all the bosons accumulate in the center of the trap.However, faster rotations lead to breakup of the density profile into three clouds.We have also computed the ground state density at the mean-field level and it shows essentially identical features as the many-body level density in real space.
To further understand the breaking of the ground-state density, we discuss the behaviour of the natural occupations N as a function of the rotation frequency ω r as shown in Figure 7.It is evident from Figure (7) that the system remains fully condensed, i.e., with the inclusion of rotation till about ω r = 1.1.From about ω r = 1.2 onwards, the condensate starts to deplete with gradual decrease in the popula- tion of the first natural orbital followed by corresponding increase in the populations of the second and third natural orbitals, n 2 N ∼ n 3 N ∼ 10 −4 [see Figure S2(b) in the supplemental material].For a faster rotation, say ω r = 2, the system transits to a fragmented state with decrease in the population of the first natural orbital followed by macroscopic population of the second and third natural orbitals.Further increase in the rotation to ω r = 2.5 leads to a three-fold fragmented state having the natural occupations n 1 N ∼ 60% and n 2 N ∼ n 3 N ∼ 20% of the first, second, and third natural orbitals, respectively.It is observed that to achieve a fully three-fold fragmentation in this three-fold symmetric trap, we would need a much faster rotation in comparison with the elongated trap.Therefore, in case of a weakly interacting bosons confined in the three-fold symmetric trap given by Equation ( 34), switching on rotation leads to a transition from a fully condensed state to three-fold fragmentation.
Let us move to the behaviour of average angular-momentum per particle LZ /N , as a function of the rotation frequency ω r for the three-fold symmetric trap computed both at the mean-field and many-body levels.It is observed from Figure ( 8) that the mean-field and many-body angular momenta exactly match each other for all ω r .The angular momentum remains minimum till the rotation frequency of about ω r = 1.1.For about ω r = 1.2, the rotation produces a state where significant value of angular momentum generates with LZ /N > 1 at which the breakup of the density is observed.The angular momentum gradually increases with further increase in the rotation of the condensate.Even for strong rotation, the angular momentum computed at the mean-field and many-body levels coincide each other.Therefore, we can conclude that, at least for the ground state, the angular momentum and its variance [see Figure S6(b) in the supplemental material] are not good indicators for many-body effects.Nonetheless, we stress that the acquisition of angular momentum in the condensate follows the breakup of the density, the emergence of depletion, and the eventual fragmentation in the three-fold symmetric trap.Now we move to the discussion of the impact of rotation on the behaviour of the manyparticle variances of the position and momentum operators along the x-and y-directions, to further characterize the many-body properties of the rotating condensate undergoing breakup.
Figure 9(a) displays the behaviour of the many-particle position variance per particle Ŷ as a function of the rotation frequency ω r in the three-fold symmetric trap computed at the mean-field and many-body levels.It is observed that the mean-field and many-body monotonously depicting the spreading and finally the breakup of the density.However, the many-body 1 N ∆ 2 X, Ŷ first increase and then, after reaching a maximal value starts decreasing with further increase in ω r , because of the depletion and eventual fragmentation.One of the important features of the position variance is that the variance along the x-and y-direction exactly coincide each other for both the mean-field and many-body regimes, thereby indicating the fact that the three-fold rotational symmetry of the condensate is preserved even for fast rotations.Clearly, there is no anisotropy of the variance for all ω r in this trap.Now we move to a more complicated system, a four-fold symmetric trap, to show the stability of the ground-state properties found above for the three-fold symmetric trap.The potential of the four-fold symmetric trap is given by Here, we consider M = 4 self-consistent orbitals to obtain the ground-state properties of bosons under rotation in this trap.We also compute the results with M = 8 self-consistent orbitals to verify the numerical convergence, see the supplemental material.The range of rotation frequencies is taken to be ω r = [0, 2.0] for the following study.
We computed the ground state energy E/N , and found it to display similar pattern as that for the elongated and three-fold symmetric traps.That is for slow rotation, E/N remains almost constant and then E/N drops gradually with increase in ω r .In addition, we found that the mean-field and many-body energies are practically identical for all ω r .The results are shown in the supplemental material, see Figure S1(c).
Figure 10 shows the behaviour of the ground-state densities per particle of a rotating BEC confined in the four-fold symmetric trap for three different rotation frequencies ω r at For a deeper understanding of the many-boson density profile, the behaviour of the natural occupations N as a function of the rotation frequency ω r is analyzed in Figure (11).It is found that the system preserves the fully condensed state, with , till a rotation frequency of about ω r = 1.0.Further increase in ω r leads to slow depletion of the condensate.For a faster rotation, ω r = 2.0, four-fold fragmentation of the condensate with finite population of all the four natural orbitals is observed.
Let us discuss the behaviour of the average angular momentum per particle as a function of ω r , see Figure 12.The angular momentum remains minimal for ω r = 1.From about ω r = 1.1, the rotation produces a state where significant value of angular momentum generates with LZ /N > 1 at which the breakup of the density is observed.The angular momentum gradually increases with further increase in the rotation.Even for strong rotation, the angular momenta computed at the mean-field and many-body levels coincide each other.Finally, it can be concluded for the four-fold trap as well that, at least for the ground state, the angular momentum and its variance [see Figure S6(C) in the supplemental material] do not precisely signifies many-body effects.However, we can conclude that accumulation of angular momentum in the condensate, the breakup of the density, emergence of depletion, and eventual fragmentation in the four-fold symmetric trap are in sync similar to the elongated and three-fold symmetric traps.
Figures 13 displays the behaviour of the position and momentum variances computed at the mean-field and many-body levels as a function of rotation frequency ω r for fourfold symmetric trap along the x-and y-directions.The mean-field and many-body position Ŷ coincide till about ω r = 1.For faster rotations, from about ω r = 1.1 onward, the mean-field position variance deviates from the many-body position variance and the former is always larger than the latter.Further, the position variances along the x-and y-directions exactly coincide with each other both at the mean-field and many-body levels similar to the three-fold symmetric trap, see Figure 9 displays similar feature as that of the position variance, both at the mean-field and many-body levels.In particular, mean-field is larger than manybody from about ω r = 1.5.
Finally, we can conclude that the depletion, angular momentum, and the position and momentum variances follow hand in hand both in the three-fold and four-fold symmetric traps.The rotating interacting bosons acquire unique many-body properties while undergoing breakup in space.

IV. CONCLUDING REMARKS
In the present work, we have studied the impact of rotation on the ground state of weakly interacting bosonic atoms confined in two-dimensional anharmonic potentials, first in, an elongated trap and, subsequently, in three-fold and four-fold symmetric traps.Here, the multiconfigurational time-dependent Hartree method for bosons, which is particularly suitable to describe many-body properties, is employed to investigate the ground-state energy, density, the depletion and fragmentation, angular momentum, and finally, many-particle variances as a function of the rotation frequency to characterize the correlations present in the system.
In the elongated trap, it is observed that the ground-state density breaks up into two clouds with rotation.The splitting of the density is followed by the emergence of the twofold fragmentation.Interestingly, the ground state exhibits opposite anistoropy both for the many-particle position and momentum variances when computed at the many-body and mean-field levels.The rotation squeezes the position variance in the elongated direction and the momentum variance in the narrow direction, thereby producing unique correlations.
Finally, a synchronized pattern among the density breakup, eventual fragmentation, acquisition of angular momentum in the condensate, and many-particle variances is observed.
When the angular momentum sets in the breakup of density and the many-particle position and momentum variances start to increase.
For the three-fold and four-fold symmetric traps, the ground-state density eventually splits into three and four clouds, respectively, with the inclusion of the rotation.Side-byside, the rotation leads to transition from condensed to three-fold and four-fold fragmented condensates, respectively, at the many-body level of theory.We find that the depletion, the accumulation of angular momentum, and finally, the increase in the variances of position and momentum follow hand in hand.
The rotating frame in our work can be viewed as a specific case of synthetic gauge fields.
In future, continuation of this investigation includes the extension of the study in various synthetic gauge fields.It would be interesting to explore the ground-state breakup, eventual condensation or fragmentation, and finally, various correlations of the condensate in presence of synthetic gauge fields.It would be fascinating to discover many-body features beyond the capacity of a "simple" rotation.

SUPPLEMENTAL MATERIAL TO FRAGMENTATION AND CORRELATIONS IN A ROTATING BOSE-EINSTEIN CONDENSATE UNDERGOING BREAKUP
This Supplemental material provides further many-body analysis that supports our main results.Sec.S1 benchmarks the multiconfigurational time-dependent Hartree for bosons (MCTDHB) method [72,76,85] in a harmonic interaction model (HIM) under rotation that gives an exactly solvable many-body model.This method is employed for the numerical computations of the main results.Sec.S2 gives the convergence checks of the ground state energy per particle of our main results with respect to the number of self-consistent orbitals in different anharmonic traps, namely, elongated, three-fold symmetric, and four-fold symmetric traps.Sec.S3 illustrates the convergence of natural occupations in these confining traps with respect to the number of orbitals.Sec.S4 presents the many-particle variances discussed in the main text, along with their convergences with respect to the number of selfconsistent orbitals.This section is further divided into two subsections.The first subsection deals with convergence of the position and momentum variances.Convergence with respect to the number of self-consistent orbitals of the expectation value of the angular momenta and angular momentum variance is elaborated in the second subsection that supports our main results.

VI. BENCHMARKING OF MCTDHB IN THE ROTATING FRAME: AN EX-ACTLY SOLVABLE MANY-BODY MODEL
In this section, we present a numerical benchmark of the MCTDHB method in the rotating frame in two dimensions 2D.The many-body Hamiltonian in the rotating frame can be written as where the single-particle Hamiltonian is ĥ(r) = − 1 2 ∂ 2 ∂r 2 + V (r).We consider a harmonic interaction model (HIM) in 2D.In this HIM, the confining trapping potential V (r) is a translated 2D harmonic oscillator of the following form and the expectation values of the many-particle momentum operators are given as Similarly, the expectation value of the many-particle angular momentum operator reads The many-particle position variances per particle are obtained as, and the many-particle momentum variances per particle are the following It is evident that the many-particle variances of the position and momentum operators in the rotating frame boil down to those for the non-rotating case.
However, the many-particle angular momentum variance per particle is affected by the rotation and takes on the following form In Table I, we show the numerical convergence of the energy per particle E N .Exact analytical versus numerical MCTDHB ground state energies are obtained for three different rotation frequencies ω r .The numerical convergence of the many-body energy per particle is achieved by using M = 6 self-consistent orbitals.Thus, we conclude that for all ω r , the ground state energy of the system is converged.Table II shows the comparison of the exact analytical and numerical MCTDHB results of the many-particle position variances for various rotation frequencies ω r .The numerical convergence of MCTDHB results computed with M = 6 self-consistent orbitals are clearly evident from the table.
Table III and IV correspond to benchmarking of the many-particle momentum and angular momentum variances for various rotation frequencies computed with M = 6 selfconsistent orbitals.As the many-particle variances of the momentum and angular momentum are very sensitive quantities, it is generally more difficult to achieve numerical convergence of the analytical and numerical results.However, from these tables, it is evident that the MCTDHB results are converged with exact analytical results.The convergence of the many-body energies reported in the main text with M self-consistent orbitals is here verified by taking 2M self-consistent orbitals and comparing the respective results.
Figure 14 shows the behavior of the ground state energy E N as a function of the rotation frequency ω r for different numbers of self-consistent orbitals for three confining potentials given by Equations ( 33), (34), and (35).The inset in all the three panels of Figure 14 display the energy difference between the energies computed for two different orbital numbers M , with M = 4, 8 orbitals as evident from Figure 14(c).We conclude that the ground state energies are fully converged with respect to the number of orbitals in all the three confining potentials at all rotation frequencies.

VIII. CONVERGENCE OF THE NATURAL OCCUPATIONS
The behavior of the natural occupations n j /N for the confining potentials defined in Equations ( 33), (34), and ( 35    It is evident that the convergence with respect to the number orbitals of the many-particle variances of the position and momentum operators for all the three confining potentials is excellent.

B. Angular momentum properties
Before we discuss the angular momentum variance, we wish to show that the expectation value of the angular momentum operator discussed in the main text converges at the manybody level.Then, we proceed to the angular momentum variance, and augment the text by its analysis.Figure (18) proves the convergence of the expectation values of angular momentum operator LZ = N j=1 lz j per particle 1 N Ψ| LZ |Ψ , with respect to the number of orbitals M for the three confining potentials; the elongated trap [Equation (33)], the three-fold symmetric trap [Equation (34)] and the four-fold symmetric trap [Equation (35)].Finally, we discuss the convergence of the many-particle angular momentum variance  three confining potentials at all ω r .The insets of Figure (19) demonstrate that in the present rotating traps and interaction strengths the mean-field and many-body angular momentum variances essentially coincide.

FIG. 1 .
FIG. 1. Behaviour of the ground state energy per particle E/N as a function of the rotation frequency ω r in an elongated trap computed at the mean-field (MF) and many-body (MB) levels (with M=1 and M=4 self-consistent orbitals, respectively).The inset shows the energy difference between the mean-field and many-body energies.All quantities shown are dimensionless.

Figure 1 1 FIG. 2 .
Figure1shows the behaviour of the ground-state energy per particle E/N, in the rotating frame computed both at the mean-field and many-body levels.Initially, E/N remains almost constant for slow rotation.Then, E/N drops gradually with further increase in ω r as evident FIG. 3. Pathway from condensation to fragmentation in the rotating elongated trap.two-fold fragmentation is observed with increase in rotation.The two leading natural occupations n 1 /N and n 2 /N , are shown as a function of the rotation frequency ω r .M = 4 self-consistent orbitals are used.The third and fourth natural occupations satisfy n 3 /N , n 4 /N ≤ 10 −5 , see the supplemental material.All quantities shown are dimensionless.

n 3 N ∼ n 4 N 3 N and n 4 N
≤ 10 −5 with the inclusion of rotation till about ω r = 1.1 [see Figure S2(a) of the supplemental material, which displays the depletion as a function ω r on a log scale].As ω r increases further, the first natural occupation number n 1 /N , falls of gradually following an increase in population of the second natural occupation number n 2 /N .The other two natural occupation numbers remain almost the same as n 3 N ∼ n 4 N ≤ 10 −5 .For faster rotation at ω r = 2, the state becomes essentially fully two-fold fragmented with natural occupations of n 1 N ≈ n 2 N ≈ 50%.This signifies equally populated two leading natural orbitals, whereas n remain essentially unpopulated.For an intermediate frequency, ω r = 1.3, we observe a scenario where n 1 shows a deep followed by a peak in n 2 .To understand this feature, we zoom in at the intermediate points between the rotation frequencies ω r = 1.2 and ω r = 1.4 as shown in the inset of Figure 3.A smooth transition is found, from coherence via loss of coherence to build up of coherence.This transition might suggest the presence of a resonant-like behaviour of the interacting bosons in the elongated 2D trap at this specific rotation frequency.To intermediately summarize, for weakly interacting bosons confined in an elongated trap [Equation (33)], inclusion of rotation triggers a transition from a fully condensed state to a fully two-fold fragmented state with equally populated two leading natural orbitals.Thus, it can be concluded that the rotation can be used as a tool to manipulate fragmentation.

FIG. 4 .
FIG. 4. Expectation value of angular momentum operator LZ /N , computed at the mean-field (MF) and many-body (MB) levels [with M = 1 and M = 4 self-consistent orbitals, respectively] as a function of the rotation frequency ω r for the elongated trap.Actual data are points.The continuous curves are only to guide the eye.All quantities shown are dimensionless.

FIG. 5 .
FIG. 5. Dependence of the the many-particle position and momentum variances on the rotation in the elongated trap.Shown are (a) 1 N ∆ 2 X , (b) 1 N ∆ 2 Ŷ , (c) 1 N ∆ 2 PX , and (d) 1 N ∆ 2 PY as a function of ω r at the many-body level (MB) [M = 4 self-consistent orbitals] and at the mean-field level (MF) [M = 1 self-consistent orbital].All quantities shown are dimensionless.

1 N ∆ 2 PX 1 N ∆ 2 PY displays similar behaviour as the 1 N ∆ 2 X
, PY along the x-and y-directions as a function of rotation frequency ω r in the elongated trap [Figure 5(c)-(d)].It is observed that unlike the position variance, the mean-field and many-body momentum variances along the x-direction are almost similar, see Figure5(c).The momentum variances computed at the mean-field and many-body levels gradually increase from ω r = 1.2 onwards.This corroborates the narrowing of the density lobes along the x-direction in real space.The momentum variance in the y-direction displays a completely different picture, see Figure5(d).The mean-field and many-body variances match each other till ω r = 1.1.However, for a faster rotation, the momentum variance computed at the mean-field deviates from that computed at the many-body level in the y-direction.as clear from Figure5(a).Thus, excitations along the y-direction plays a role in the momentum space as a result of rotation.

FIG. 6 .
FIG. 6.The one-body densities per particle are shown for three different rotation frequencies ω r in the three-fold symmetric trap at the many-body level.M=3 self-consistent orbitals are used.The density computed at the mean-field (not shown) and many-body levels depict identical features for all ω r .All quantities shown are dimensionless.

FIG. 7 .
FIG. 7. Pathway from condensation to fragmentation in the rotating three-fold symmetric trap.three-fold fragmentation is observed with increase in rotation.The three leading natural occupations n 1 /N , n 2 /N and n 3 /N , are shown as a function of the rotation frequency ω r .M = 3 self-consistent orbitals are used.All quantities shown are dimensionless.

FIG. 8 .
FIG. 8. Expectation value of angular momentum operator LZ /N , computed at the mean-field (MF) and many-body (MB) levels [with M = 1 and M = 3 self-consistent orbitals, respectively] as a function of the rotation frequency ω r for the three-fold symmetric trap.All quantities shown are dimensionless.

1 N ∆ 2 X , ŶFIG. 9 .
FIG. 9. Dependence of the the many-particle position and momentum variances on the rotation in the three-fold symmetric trap.Shown are (a) 1 N ∆ 2 X and 1 N ∆ 2 Ŷ and (b) 1 N ∆ 2 PX and 1 N ∆ 2 PY as a function of ω r at the many-body level (MB) [M = 3 self-consistent orbitals] and at the mean-field level (MF) [M = 1 self-consistent orbital] All quantities shown are dimensionless.

Finally, Figure 9
(b)  shows the behaviour of the momentum variance per particle 1 N ∆ 2 PX , PY , computed at the mean-field and many-body levels as a function of ω r .The momentum variances remain small even with increase in rotation till about ω r = 1.1.From about ω r = 1.2 onwards, the momentum variances computed both at the mean-field and manybody levels slowly start to increase.Here, the mean-field momentum variance only increases monotonously.The many-body momentum variance increases, decreases, and again increases but remain much smaller than the mean-field momentum variance, indicates the presence of depletion and fragmentation.Similar to the position variance, the momentum variance along the x-and y-directions exactly coincide each other both at the mean-field and many-body levels, due to the rotational symmetry.

2 .
FIG. 10.The one-body densities per particle are shown for three different rotation frequencies ω r in the four-fold symmetric trap at the many-body level.M=4 self-consistent orbitals are used.The density computed at the mean-field (not shown) and many-body levels display identical features at all ω r .All quantities shown are dimensionless.

FIG. 11 .
FIG. 11.Rotation leads to four-fold fragmentation in the four-fold symmetric trap.M = 4 selfconsistent orbitals are used.The variation of four natural occupations n 1 /N , n 2 /N , n 3 /N and n 4 /N are shown as a function of rotation frequency ω r .All quantities computed are dimensionless.

FIG. 12 .
FIG. 12. Expectation value of angular momentum operator LZ /N , computed at the mean-field (MF) and many-body (MB) levels [with M = 1 and M = 4 self-consistent orbitals, respectively] as a function of the rotation frequency ω r for the four-fold symmetric trap.All quantities shown are dimensionless.
(a).Side by side, the behaviour of the momentum variance 1 N ∆ 2 PX , PY

FIG. 13 . 2 PX , 1 N ∆ 2 PY
FIG. 13.Dependence of the the many-particle position and momentum variances on the rotation in the four-fold symmetric trap.(a) depicts the position variances 1 N ∆ 2 X , 1 N ∆ 2 Ŷ and (b) displays momentum variances 1 N ∆ 2 PX , 1 N ∆ 2 PY at the many-body level (MB) [M = 4 self-consistent orbitals] and at the mean-field level (MF) [M = 1 self-consistent orbital].All quantities shown are dimensionless.

) 1 .
With a solvable many-boson model in the rotating frame, we can now proceed and benchmark MCTDHB in the rotating frame.In the following simulations, we work in the units = m = ω = 1.To compute the ground state of MCTDHB in the rotating frame, we consider N = 10 weakly interacting bosons with repulsive interaction parameter Λ = −0.The box extensions of the confining potential are chosen to be = [−12, 12) × [−12, 12) with 128 × 128 DVR functions to represent each of the orbitals.We translated the confining potential by L = 2 in Equation (19).
Benchmarking of the expectation values and variances of the many-particle position operators along the x-and y-directions with respect to the rotation frequency ω r computed with M = 6 self-consistent orbitals.All the quantities are shown in dimensionless units.

FIG. 14 .
FIG. 14. Many-body ground state energy per particle E N as a function of the rotation frequency ω r .(a) for the elongated trap M = 8 self-consistent orbitals are used to show the convergence of the M = 4 results, (b) for the three-fold symmetric trap M = 6 self-consistent orbitals are used to check the convergence of M = 3 results, and finally in (c) M = 8 self-consistent orbitals are used to check the convergence of M = 4 results of the four-fold symmetric trap.The inset shows the energy difference per particle as defined in Equation (36).The quantities shown are dimensionless.

Figure 14 (
Figure 14(a) corresponds to the elongated trap [Equation (33)].The convergence of energy is checked for two different orbitals M = 4, 8.We find that the two energy curves E N computed with M = 4, 8 orbitals fall on top of each other for all the frequencies.It is evident that our results converges with respect to the number of orbitals for elongated trap.Similar convergence is observed in case of the three-fold symmetric trap with M = 3, 6 orbitals in Figure 14(b) [see Equation (34)] and the four-fold symmetric trap [Equation (35)]

4 N , n 5 N , and n 6 N 4 N ∼ n 5 N ∼ n 6 N
Figures 15 for different orbital numbers M .It is observed that all the natural occupations are fully converged for the three confining potentials.The inset in Figure 15(b) shows the variation of the smallest natural occupations n 4 N , n 5 N , and n 6 N with the rotation frequency ω r for the three-fold symmetric trap.All the three natural occupations stays minimal with n 4 N ∼ n 5 N ∼ n 6 N ∼ 10 −6 for all ω r .Similarly, the inset of Figure 15(c) depicts the natural occupations n 5 N ∼ 10 −6 , n 6 N ∼ n 7 N ∼ n 8 N ∼ 10 −7 − 10 −6 for the four-fold symmetric trap.The convergence of the natural occupations is evident in all three confining traps.

A 2 PY
Figures 16(a)-(b) show the convergence of the many-particle position variances and Figures 16(c)-(d) correspond to the many-particle momentum variances for the elongated trap [Equation (33)] with M = 4, 8 orbitals.

FIG. 17 .
FIG. 17. Convergence of the many-particle position (a) and momentum variances (c) with M=3,6 self-consistent orbitals for the three-fold symmetric trap.Panels (b) and (d) correspond to manyparticle position and momentum variances with M=4,8 self-consistent orbitals for the four-fold symmetric trap.

FIG. 18 .
FIG. 18. Convergence of the expectation values of the angular momentum operator per particle.(a) for the elongated trap with M=4,8 self-consistent orbitals, (b) for the three-fold symmetric trap with M=3,6 self-consistent orbitals and (c) for the four-fold symmetric trap with M=4,8 self-consistent orbitals.

1 N ∆ 2 LZ
in all the three confining traps.It is evident from Figure(19) that the angular momentum variances are fully converged with respect to the number of orbitals for all the

FIG. 19 .
FIG. 19.Convergence of the angular momentum variance.(a) for the elongated trap with M=4,8 self-consistent orbitals, (b) for the three-fold symmetric trap with M=3,6 self-consistent orbitals and (c) for the four-fold symmetric trap with M=4,8 self-consistent orbitals.The insets of (a), (b), and (c) demonstrate the angular momentum variances at the mean-field and many-body levels that are essentially the same in our rotating traps.

TABLE I .
Benchmarking of the ground state many-body energy E/N with respect to the rotation frequency ω r .M = 6 self-consistent orbitals are used.The ground state energies shown are in dimensionless units.