Analogue gravity simulation of superpositions of spacetimes

Taking the principles of quantum mechanics as they stand and applying them to gravity, leads to the conclusion that one might be able to generate superpositions of spacetimes, at least formally. We analyze such a possibility from an analogue gravity perspective. We present an analogue toy model consisting of a Bose-Einstein condensate in a double-well potential and identify the states that could potentially be interpreted as superposition of effective spacetimes. These states are unstable and the source of instability from a microscopic point of view can be related to the absence of a well-defined causal structure in the effective geometric description. We explore the consequences of these instabilities and argue that they resonate with Penrose's ideas about the decay that superpositions of states with sufficiently different gravitational fields associated should experience.


I. INTRODUCTION
The search for a theory of quantum gravity has been one of the main driving forces in theoretical physics over the last century. Many approaches exist toward building a successful theory of quantum gravity, among them string theory [1,2] and canonical quantum gravity [3] in its modern formulation in terms of loop quantization [4] are the leading approaches; although there are many other approaches as Asymptotic Safety [5,6], Causal Dynamical Triangulations [7], or Causal Sets [8]. In any case, it is fair to say that none of them is regarded as being completely satisfactory. The laws of quantum physics are maintained in all of these theories and one simply applies different quantization schemes to the geometric degrees of freedom, or regards them as emergent and applies the standard quantum rules to other microscopic degrees of freedom. However, one remarkable property shared by many of these approaches is that it is possible to generate superpositions of almost classical states, i.e. states that we can almost describe as approximate smooth geometries at low energies. This very same notion of superposing semiclassical spacetimes has also a counterpart with the electromagnetic field, where even Schrödinger-cat-like states can be generated [9].
Our aim in this work is to dig into the physical meaning of superposing spacetimes, trying to simulate them within an analogue gravity framework [10]. The analogue gravity program intends to shed new light into gravitational physics by reproducing some of its behaviours in laboratory systems. Currently, analogue gravity is just able to reproduce kinematical properties of gravitational fields, not dynamical properties. In that sense, since our understanding of kinematical aspects of gravity is based on the theory of General Relativity, in which spacetime is described as a curved geometry [11], we try to to simulate the propagation of signals on top of effectively curved geometries. Among the analogue gravity systems at our disposal, fluid systems and, in particular, quantum fluids based on Bose-Einstein Condensates (BECs) as their substratum, have proved to be very useful both at the theoretical and at the experimental level [12,13]. This analogy relies on the long-wavelength regime of acoustic perturbations: such perturbations behave as if they were a scalar field moving in an effectively curved spacetime.
Whereas the standard analogy for a classical spacetime is well understood (see [14] for a systematic presentation) our aim here is to explore whether such analogy can be extended to the case of having superpositions of semiclassical spacetimes. The effective spacetime metric depends on the characteristics of the macroscopic wave function of a BEC which is a fully quantum mechanical object. Therefore, within this framework it appears in principle possible to engineer situations such that there exist acoustic waves perceiving a superposition of two different spacetimes, i.e. two different causalities. Whether this can be done or not is the main inquiry that we address in this work. For that purpose, we study a toy model consisting in a BEC in a double-well potential and look for the states that resemble a superposition of spacetimes.
Given the insatisfactory situation regarding our current theories of quantum gravity, there has been a relatively small but steady trend pointing out that maybe one should not quantize gravity but gravitizate quantum mechanics, in Penrose's words [15]. One of the striking phenomena that we face when assuming that the laws of quantum physics are valid at all scales and for arbitrary objects is that we observe no phenomenology associated with quantum superpositions of macroscopic states although they are allowed within the conceptual framework of quantum mechanics. The standard explanation relies on a process known as environmental decoherence [16][17][18]. However, Penrose has heuristically argued that gravity could be at the heart of some more fundamental source of decoherence [15,[19][20][21]. As we will see, our results in the analogue gravity model strongly resonate with Penrose's ideas, although some differences also appear. Essentially we find that whenever some aspects of our analogue quantum gravity involve pushing the causal behaviour of General Relativity to its limits, either by superposing causalities as we do here, or by engineering chronologically pathological spacetimes as in [22], we find that there exist kinematical limitations toward building such states. However, the dynamical mechanism leading to their destabilization when we try to build them is highly dependent on the particular dynamics of the analogue model.
Here is an outline of the article. In Section. II we review the standard approach toward analogue gravity in which single spacetimes are simulated. Section III constitutes the core of the article. In subsection III A we will review the main properties of the toy model that we will consider for bosons in a double well, the Bose-Hubbard model with two sites. In subsection III B, we will explain which are the main properties that this toy model needs to fullfill so that its sound wave excitations admit an effective description in terms of a Klein-Gordon equation propagating on top of a curved geometry. In sub-section III C, we will explore the stability of the ground states of the Bose-Hubbard model depending on the possible values of its parameters. In subsection III D, we will study the implications of this stability analysis and discuss the constraints that we have to impose to the model in order to describe an effective spacetime. Furthermore, we will explore the possibility of building an effective superposition of spacetimes in this model. Finally, in subsection III E, we will conclude by extending our discussion to the possibility of creating superposition of similar geometries in the same location but with different causalities. In Section. IV we will discuss the similarities and differences between our findings in the analogue gravity system and Penrose's ideas on the superpositions of spacetimes. We will finish the paper by summarizing the main conclusions that can be drawn up from this work in Section V.
Notation and conventions. We work in units in which = c = 1, unless explicitly stated. We will use the signature (−, +, +, +) for the spacetime metric. Greek indices (µ, ν, ...) will run from 0 to 3, representing spacetime indices, whereas latin indices (i, j...) will run from 1 to 3 and represent spatial indices. Einstein's summation convention is used throughout the work unless otherwise stated.

II. ANALOGUE GRAVITY: THE STANDARD PICTURE
In general terms, for every barotropic and inviscid 1 fluid whose flow is irrotational, it is possible to arrange the velocity potential φ describing sound waves on top of the fluid background to satisfy a massless Klein-Gordon equation [14] 1 where we have introduced the metric g µν , its determinant g = detg µν and its inverse g µν , defined in terms of the physical parameters as with c the local speed of sound, v i the velocity of the fluid, ρ its density, and where we are writing everything in the (t, x) laboratory coordinates. This acoustic metric provides the sound-cones at each point of the spacetime. One of the main properties of this acoustic geometry is that it is causally well-behaved. This was noticed in [14], where it was shown that this metric automatically inherits the stable causality property [11,23] For a careful analysis on this point and a discussion on the possibility of building causally pathological acoustic geometries, see [22]. This will be the departure point of our discussion. Actually, up to this point, we have not made any comments on the nature of the fluid under consideration. It can be either an ordinary classical fluid or a quantum fluid. We will focus here on a particular kind of quantum fluid, a Bose-Einstein condensate, which is experimentally feasible to simulate spacetimes of interest [13,24]. Actually, among the whole plethora of analogue systems that we have [10], BECs are one of the best analogues since they are very clean and insensitive to external noise.
From the perspective of this work the crucial characteristics of using a BEC as analogue system is that it is an essentially quantum system: it is close to the quantum vacuum state and in this sense very far from classicallity. The fluid analogy with a BEC starts from describing many-body system of bosons through a quantum field Ψ, whose dynamics is encoded in with g the leading order interaction among bosons. At sufficiently low temperatures, these systems of weakly interacting bosons can develop condensation in the sense that a macroscopic number of particles populates the same state. To describe this condensation, it is customary to break the field into a macroscopic condensate and a fluctuation as Ψ = ψ + δψ, with the condensation expressed as Ψ = ψ. Using a Madelung representation and considering a long wavelength approximation one can describe the system as quantum acoustic perturbations moving in a curved background spacetime of the form (2). The velocity of the background fluid is v = ∇θ/m associated with the phase θ of ψ: The density ρ 0 = mn c and c s represents the speed of the phonons in the medium As mentioned above, for this picture to be consistent, we must require that the wavelengths of the sound excitations are large. In fact, they have to be larger than the so-called healing length of the condensate ξ = 1/(mc s ), i.e. λ ≫ 2πξ. This healing length is the characteristic length scale at which the hydrodynamic approximation breaks down and one is able to unveil the acoustic description of the fluid begins to break down.
Given this setup that we have introduced in which we have a quantum fluid, it is possible to think about generating superpositions of macroscopic states of the BEC. Using the language of analogue models in principle one could, for example, engineer a BEC to be in a superposition of two rather different acoustic black hole spacetimes: two black holes with different surface gravities. Another option is to engineer the BEC to be in two equivalent spacetimes, for example a single acoustic black hole, centered around different laboratory points. In a quantum gravity language, this could be taken also to be the geometries of two different coexistent universes. In the following we are going to discuss this last example by analyzing a BEC in a double-well and how the different possible ground states of the system depend on the values of the different parameters. Thus, our inquiry is the following: what does this superposition mean from an analogue gravity point of view? Is it possible to interpret this superposition as an effective superposition of spacetimes? And what is the fate of acoustic excitations propagating on top of it? The next section is deserved to analyze a toy model in which this questions can be sharply formulated.

A. Condensates in a double well
Let us assume that we have scalar bosons in a symmetric double-well potential where we will label the wells by i = 1, 2. We will make the simplifying assumption for the moment that there is only one relevant state in each well and particles within a well have a contact interaction controlled by the parameter U which can be either repulsive (U > 0) or attractive (U < 0). Also, we will assume that we have a term describing tunneling between the wells controlled by the parameter t which is the amplitude of probability for tunneling between both wells (not to be confused with the time t introduced in the previous section). The bigger the parameter t, the higher the potential barrier will be. The Hamiltonian describing these particular interactions is the Bose-Hubbard Hamiltonian [25] with just "two sites" in the lattice where a † i , a i represent the usual creation and annihilation operators that annihilate and create particles in the well i, with n i = a † i a i being the number of particles in the well i. Furthermore, we will work under the assumption of having a fixed number of particles N = n 1 + n 2 . For the non-interacting Hamiltonian, i.e. Eq. (7) with only the tunneling term, the solutions are explicit and can be found via an ordinary Bogolyubov transformation. In the regime in which the tunneling is highly suppressed, just the interacting term survives. It is diagonal in the number basis and the solutions are also explicit since they simply correspond to the number basis. Notice that we are omitting the hats to denote the quantum operators from now on, to avoid a cumbersome notation.
We will now describe the main properties of the states which will be interesting for our purposes: the ground state of the free Hamiltonian and the ground states for the Hamiltonian with t = 0, i.e., the ones with no hopping between the two wells, either with attractive U < 0 or repulsive interactions U > 0. They will serve as approximate ground states for the regimes of hopping domination t ≫ |U |N and interaction domination t ≪ |U |N , respectively. The intermediate ground states between both regimes are not so useful for our purposes and can be found numerically. However, understanding the evolution of the system from one regime to the other as we adiabatically vary the parameters sheds some light onto the behaviour of the model and it can be well captured by the ansatz provided in [26].
Prior to discussing the ground states on the extreme regimes above, we need to introduce the concept of the single-particle reduced density matrix ρ sp . Such density matrix is defined as where ρ represents the density matrix of the system. According to Leggett's nomenclature [27], a condensate will be fragmented if its single-particle density matrix has two macroscopic eigenvalues, i.e., two eigenstates with eigenvalue of order O(N ). Intuitively, this would correspond to a situation in which condensation occurs simultaneously in two different single-particle states. Armed with this tool, let us now analyze the three ground states of interest: the free case, the repulsive interacting case and the attractive interacting case. The corresponding states are pictorially represented in Fig. 1 Free condensate: The free case can be explicitly solved via an ordinary Bogolyubov transformation, described by the following SO(2) rotation within the space Thus, we find that for the single-particle subspace the eigenstates of the Hamiltonian are the symmetric and antisymmetric states with energies −t and t, respectively. Thus, the ground state will be given by the fully symmetric state because they do not interact among themselves. |0 represents the Fock vacuum, i.e., the state annihilated by a 1 and a 2 . The C stands for the coherent state, since it corresponds to N bosons coherently delocalized among the two wells with equal probability of finding them in The opacity in the figures above represents how localized the state is. The coherent state is a completely delocalized state with a non-vanishing projection onto all of the basis number states, the Fock state corresponds to two localized places at which the state exhibits a peak and the cat state corresponds to a delocalized state corresponding to all the particles being in one of the wells or the other.
one of the wells. The single particle density matrix for this state reads This state has a single macroscopic eigenvalue N . According to Leggett's classification [27], it corresponds to a non-fragmented condensate since we just have one macroscopic eigenvalue of the single-particle density matrix. Physically, it corresponds to a condensation on the single-particle state which is approximately a superposition of the two Gaussians peaked around the center of each well in position space, describing the ground state of each well. Since the ground state |C is a linear combination of number states |n 1 , n 2 = a †n1 1 a †n2 2 |0 / √ n 1 !n 2 !, the number of particles in each well has enormous fluctuations. We can compute them by writing the coherent state |C in the number basis. For even N (we consider such case for simplicity), we have where |ℓ ≡ N 2 + ℓ, N 2 − ℓ and The fluctuations in the number of particles on each well are given by [26] with i = 1, 2. Let us move on to discuss the strongly interacting case, i.e. t = 0, U = 0. For this purpose, taking into account that the total number of particles is conserved n 1 + n 2 = N , we can rewrite the Hamiltonian as which parametrizes the interaction by the difference in the number of particles on each well. Strong repulsive interactions: We will begin with the repulsive interactions (U > 0). Such case clearly favors the minimum difference in the number of particles between the two wells. Thus, for even N , the ground state for the repulsive interactions is clearly where the F stands for the Fock state. Its single particle density matrix is which corresponds to a fragmented condensate since it has a macroscopic eigenvalue N/2 with multiplicity 2. It corresponds to two uncorrelated condensates (notice that the state is a product state) having half the particles, each one located at one of the two wells. Such a result for the ground state in this limit is quite intuitive since having one more particle than the half of them on one of the two wells would imply paying a penalty in energy. The fluctuations in the number of particles on each well vanish identically Thus, in this case, we have a fragmented condensate with a well-defined number of particles on each well. Strong attractive interactions: For attractive interactions U < 0, the situation is the opposite: the favored states are those containing a huge difference in the number of particles on each well. This means that the ground state is the subspace spanned by the vectors {|N, 0 , |0, N }. However, if we begin with the coherent state |C , the ground state of the free theory, and we turn on adiabatically the attractive interactions, the state that we will reach once t becomes negligible will be a concrete superposition of both vectors (further details of how this state is reached will be provided at the end of the section). Actually, it corresponds to a Schrödinger-cat like state It is also fragmented since its single-particle denstiy matrix has also two eigenvalues being equal to ρ sp (F ), the one that we had for the Fock state. Particle number fluctuations on each of the wells allow us to distinguish between both states though: in this case they do not vanish but Thus, we have seen that although the Fock and cat states are fragmented according to the standard definition of Leggett (their one-particle density matrices have two macroscopic eigenvalues), they physically correspond to very different states. While the Fock state has a definite number of particles on each well, the cat state has an indefinite number of particles on each well; the latter corresponds to having the bosons delocalized. Until now, we have presented the three states that will be relevant for our purposes. They are the coherent state |C , which is the ground state of the free theory t = 0, U = 0; the Fock state |F , which is the ground state of the strongly repulsive interacting theory t = 0, U > 0; and the cat state |cat , which is the ground state of the strongly attractive interacting theory t = 0, U < 0. The coherent state corresponds to a nonfragmented condensation while the Fock and cat states correspond to fragmented condensates. With respect to fluctuations in the number of particles on each well, the Fock state has no fluctuations on the number of particles on each well, while the coherent and cat states have enormous fluctuations on the number of particles on each well. These are the relevant features of these three states which will play a role in building our analogue models. We now provide additional details of how the model behaves for intermediate regimes.
Intermediate interactions: We will analyze what happens as we turn on the interactions adiabatically, i.e., we begin with the free Hamiltonian and assume we slowly turn on the interactions in such a way that the ground state of the free theory accommodates to the ground state of the interacting theory. Although numerically it is possible to obtain explicitly the wave function of the ground state for the interacting theory, as we advanced, we will take advantage of the ansatz introduced in [26] and use their families of states as analytic states capturing the main properties of the ground states and being close to them.
These ansätze for the interacting theory will be based on the ground state wave function for the non-interacting theory, the U = 0 case of the Hamiltonian (7). The starting point will be to write the ground state in terms of the number basis and write the Schrödinger equation for this system H |Ψ = E |Ψ , which gives the following expression with The problem is equivalent then to a one-dimensional tight-binding model in a harmonic potential with non-uniform tunneling matrix elements [26]. The nonuniformity is such that wave functions Ψ ℓ with large amplitudes near ℓ ∼ 0 always have less energy than the ones that spread around different values of ℓ. Actually, in the free case, the wave function was a narrow Gaussian centered at ℓ = 0. With this expression for the eigenvalue problem, let us start by considering the repulsive case (U > 0). These interactions make the coherent state (13) squeeze in an even narrower distribution. The family of states introduced in [26] for capturing this evolution from the coherent to the Fock state is As σ 2 varies from √ N to small values, the initial coherent state |C starts looking much more like the Fock state. Actually we can obtain a relation between σ and the value of U by taking a continuum limit on (23), see [26] for details. In this limit, the equation reduces to that of a harmonic oscillator potential and we can obtain the value of σ(U ) since in that case the problem is exactly solvable: σ −2 = (2/N )(1 + U N/t) 1/2 . The single particle density matrix for these states reads It has eigenvalues N 2 (1±e −1/σ 2 ) and the fluctuation in the number of particles on each well is (∆n i ) 2 = σ 2 /2. As the gaussian width increases, the number of particles on each well vary from (N, 0) to (N/2, N/2) and the number fluctuations (∆n i ) 2 varies from N to 0.
Let us discuss now the case of attractive interactions U < 0. When the interaction is attractive, states having a huge difference in the number of particles on each well, or, in other words, a huge amount of particles on a single well, are favored. Thus, the effect of slowly turning on an atractive interaction is to split the Gaussian of the noninteracting state (13) into a symmetric distribution with two peaks, a process which is captured by the family of states being 2a is the separation between the peaks, σ ′ their width and K is a normalization factor. As a varies from 0 to N/2 and σ ′ reduces at the same time from 1/ √ N to 0, the coherent state of the free theory evolves to the cat state (which is reached in the U/t → −∞ limit).

B. Effective spacetimes from BEC's
We discussed above that BECs are good analogue systems. Let us dig in the properties that our toy model needs to fullfill in order for the excitations propagating on top of it to obey a Klein-Gordon equation on a curved geometry. In order to build an analogue gravity model, a crucial necessary condition is that the system can exhibit in some regime causality properties. For the Bose-Einstein condensate, this means that it should be able to accommodate the propagation of sound-like excitations within the system, which correspond to some of the collective excitations of the system. A fundamental ingredient for the propagation of this collective excitations is the existence of an effective repulsion among the constituents of the BEC. In addition to this condition, we require the additional condition that the wavelengths of the excitations are much greater than the healing length of the condensate, which roughly speaking is the length at which inhomogeneities of the atomic density composing the condensate are smoothened out, as we discussed above.
The coherent state resulting from taking U = 0 or the cat state that results from taking U < 0 will not give rise to sound excitations in principle since the condensate is in a regime in which no repulsive interactions are present. Hence, it seems that it is not possible that any quantum sound emerges. Let us focus on the cat state, since it is the state that would give rise to an analogue of a superposition of spacetimes, as long as it admitted sound excitations. One could build initially a Schrödinger-cattype state from Eq. (19) displaying no sound-like excitations initially. For example, one could adiabatically switch on the attractive interactions (U < 0) beginning with the U = 0 case through Feshbach resonances. Once the system settles in the appropriate cat state, one could turn the interaction towards the repulsive regime in such a way that sound excitations emerge. We will examine the validity of this picture in the next section through a stability analysis.

C. Stability of the ground states
We saw that even though one-particle density matrices of the Fock and cat states were the same, they described completely different states: while the Schrödinger-cat state had enormous number fluctuations, the Fock state had zero number fluctuations; that is, higher order correlation functions are needed to characterize these different fragmented states.
Moreover, these huge number fluctuations in the Schrödinger-cat state are telling us something about the stability of the system: it is obvious that such a delocalized system must tend to be unstable under local perturbations which act independently on each of the wells. In fact, let us consider a generic interaction modelled by the Hamiltonian where ǫ is a small parameter controlling the perturbation, a 1 is the usual annihilation operator for bosons in the first well and c † 1 is the creation operator for a generic excited state of the first well and, consequently, representing a local perturbation of the first well. Under this small perturbation, the Fock state (16) is robust since it simply changes to another Fock state |F ′ under the action of this Hamiltonian operator with K ′ an irrelevant normalization constant, which shows the robustness of Fock states under standard local perturbations. However, the cat state does not display such robustness. Actually, the cat state tends to collapse to a localized condensate once we take them into account with K ′′ another irrelevant normalization constant, since this state corresponds to a localized state on the first well. It corresponds to a number state with N − 1 particles on the ground state and one particle in the excited state described by the operator c † 1 . The conclusion we extract from this analysis is that the system tends to avoid being in delocalized macroscopic superpositions since they are unstable under local perturbations that do not coherently hit the system on the two places. We have just considered interactions of the type (28) which will be contained in a generic condensed matter Hamiltonian. This can be seen by noting that, in general, we can work in the second quantized formalism in which the fundamental object we consider is the field ψ(x), whose expansion in a concrete basis of one-particle states reads The previous simplification was equivalent to saying that we just focused on the mode f 1 which was an excited state of the first well. In practice we will have a tower of excited states above them that can be regarded as perturbations of the form (28).
D. Attempting to generate a superposition of two effective spacetimes Now that we have described the three main states that can appear in a Bose-Hubbard model, and their stability properties, it is time to discuss them at the light of their ability to provide superpositions of effective spacetimes.
Our aim here is to consider a model whose substratum is an intrinsic quantum system, admitting an effective geometrical description in some regime, and analyze the implications and viability of producing a state which would be understood within the effective geometrical description as a putative superposition or mixture of spacetimes. For that purpose, we first need to identify the states that would correspond to superpositions of spacetimes. We have two natural candidates for this purpose: the coherent and the cat states. From an effective spacetime perspective, the Fock state, being a product state, would correspond to an incoherent mixture of two spacetimes. This means that the propagation of sound on each of the two wells would occur independently of what is happening on the other well.
The coherent state can be taken as representing a system in a macroscopically delocalized state. However, we have seen that this state appears at the cost of eliminating the local interactions between the composing bosons. Having negligible local interactions, this state lacks the possibility of producing a rich causal physics. For instance, for U = 0 the state develops no acoustic excitations, it does not really serves as an analogue model of gravity. However, given our analysis of the previous section in which we see that there is a smooth interpolation among the states with U = 0 and U > 0, we can turn on a small repulsion without breaking the qualitative features of the state and allowing sound-like excitations. However, this delocalized state would have a wave-function with a profile varying too fast spatially to allow for the hydrodynamic approximation to be accurate. So we see that the coherent state does not seem suitable to describe superpositions of spacetimes.
We are left with the cat state as the only one that we can identify, at least momentarily, with a superposition of two different spacetimes, since it is a superposition of two states representing an effectively curved geometry on each of the wells. Actually, it seems to be the naive state that one would also build in a quantum gravity theory as a putative superposition of spacetimes.
Hence, the following question appears now: if we construct a condensate in the cat state, i.e., what we would claim to be a would-be "superposition of spacetimes" in the analogue system, what happens with the poten-tial acoustic perturbations propagating on top of such a spacetime: is there an effective spacetime resulting from this superposition? The idea is that there should be a limitation in describing the acoustic perturbations propagating on top of such a superposition as acoustic perturbations on a new spacetime. Otherwise, this would mean that the superposition of analogue spacetimes would be also a spacetime, contrary to what one would expect to happen even with real gravitational systems. Indeed, such limitation appears in the form of a huge instability of the system under consideration. As we discussed in Subsection III C, the cat state tends to decay to a Fock-like state under small perturbations. The propagation of sound-like excitations itself would be enough to destabilize such superposition. This translates into a decay of the putative superpositon of spacetimes into a mixture of two spacetimes, i.e., two effective spacetimes which behave independently from each other on each of the wells.
Apart from the instability under perturbations of this state which we have already discussed, an additional impediment appears. If we consider the evolution of such a system according to the Gross-Pitaevskii equation which describes generic fluid condensates, the interference is such that the healing length of this condensate is much bigger than the healing length of either of the condensates. Consequently, the effective spacetime picture breaks at much bigger wavelengths and we conclude that the hydrodynamic picture is, at best, hard to maintain for these superpositions. This manifests the impossibility of finding an analogue of a superposition of spacetimes.
Actually, this impossibility recalls Penrose's ideas on the possibility of superposing spacetimes [15,20,21]. On the one hand, Penrose argues that the superposition of spacetimes is an ill-posed concept in gravity, since it leads to a causal structure which is not well-defined. Since already in gravity there are hints that these superpositions are ill-posed, we could have expected in advance a fundamental limitation to find it in an analogue version. On the other hand, Penrose relates this ill-definitness of the causality of superposed spacetimes with the mechanisms of quantum state reduction. We will come back to Penrose's idea in Section IV, where we will contrast it with our findings in this analogue setup.

E. Same localization, different causality
We will now discuss the viability of generating a superposition of spacetimes, but now we will consider a different situation. Instead of superposing two spacetimes with the same causal properties but based around different locations, we will now consider superposing two spacetimes at the same location with different causal properties. Roughly speaking, instead of attempting to superpose two geometries generated by two lumps of matter spatially translated, we want to superpose two lumps of matter with different shapes and densities located around the same region. The double well Bose-Hubbard model does not allow us to consider this situation and it is hard to design a toy model that allows to consider it. However, we will try to analyze in general terms whether this is at all possible or not.
As the simplest situation, one could think of a homogeneous BEC at rest but with a superposition of sound velocities. We recall that the sound velocity in a condensate, c 2 s = gn c /m, is controlled by the coupling constant (from the |ψ| 4 term of the Gross-Pitaevskii equation) g, the effective mass of the bosons m and their number density n c . The numbers g and m are parameters of the system which are in principle not subject to quantum rules. The masses are renormalized due to the interactions and, in principle, it is even possible to achieve negative masses with suitable microscopic structure of the system [28]. On the other hand, the coupling constant can be controlled by an external magnetic field by using Feshbach resonances [29]. Then, we can ask, whether we can engineer this magnetic field to be in a quantum superposition. The magnetic field we speak about is a macroscopic entity so what we are doing with this inquiry is just translating the problem of generating a macroscopic superposition from one place to another. Instead, it is typically assumed that the values of the parameter of a system are phenomenological values to which one could in principle associate infinitely precise values. As the causality of the system depends directly on these values one could think that they imprint a classical flavor into the notion of causality. Hence, this tuning of magnetic fields does not seem to be a good way of achieving such a superposition.
The other quantity entering the definition of the sound velocity is the density of condensed particles n c . This quantity can have fluctuations if we deal with a system allowing flows of particles (grand-canonical ensembles), but in a condensed matter system the number of particles is completely fixed: the creation of new particles is forbidden by a huge energy gap. If one were hypothetically able to generate a BEC in an ultrarelativistic regime, the possibility of having a superposition of states with different number of particles would exist. But again, would it be possible to have a stable superposition of two mutually non-interacting gases each having a different number of particles? The possibility does not sound realistic, the very notion of macroscopic (thermodynamic) behaviour would try to produce states with a peaked distribution for the number of particles. Actually, the standard ensemble in terms of which condensed matter systems are described is the Grand Canonical Ensemble, which allows not only for energy fluxes but also particle density fluxes between the system and the environment. In such scenarios, thermal homogeneity is needed to reach thermodynamical equilibrium although it is not sufficient. Chemical homogeneity is also required, which means having equal chemical potentials for the bath and the system. If the number of particles is not highly peaked, it is impossible to reach a stationary state. Hence, it seems that this is also a no-go pathway, since an almost stationary situation in which the condensate is equilibrated is required for the hydrodynamic picture to be sensible at all.
The only element of the causality in a BEC remaining is the flow velocity of the fluid. The flow velocity that appears in the effective analogue gravity metric in a BEC corresponds to a macroscopic occupation of a particular phase structure for the mono-particle wave functions. Once again, the natural question is whether it is possible to have a stable superposition of condensates with different phases or not. The situation would be parallel to our discussion of the stability properties of the cat and Fock states, respectively. Depending on whether we perform a superposition of the two states with all the particles on each phase state or a more democratic Focklike state in which half of the particles are on each of the phase states, we would produce an unstable state or a stable state, respectively. The additional difference is the external potential which, for the formation of standard condensates, always favors one of the two phase states. Thus, although under the absence of external potentials it would be possible to generate a stable Fock-like state, in a real systems the particles would always tend to be projected onto one of the two fixed phases states.
Thus, this qualitative discussion suggests also the inability to perform superpositions of different causalities, i.e. to have a BEC sitting on a single location but giving rise to different causal structures for the propagation of sound-like waves. Nonetheless, further study is required to fully understand and characterize whether there is a universal constraint on the possibility of simulating such superposition of analogue geometries in arbitrary analogue systems.
Penrose argues that superposing spacetimes is an illdefined concept because the resulting state does not incorporate a clear definition of time translation. In turn, he argues that this uncertainty in the definition of timetranslation should lead to an uncertainty in the definition of energy. Furthermore, this should entail the decay of FIG. 2. We pictorially represent the two configurations that we are want to superpose. They correspond to the gravitational field associated with a certain lump of matter (a star, for instance), which is located in two different places for each of the states.
such states into mixed states where correlations among the superposed spacetimes are lost. His ideas are not at the level of a precise dynamical theory and hence we do not attempt to make a precise analogy here: we simply discuss its similarities and differences with our analogue toy model.
In the case in which we have a superposition of configurations based on different locations, for example, two spacetimes corresponding to a lump of matter (for example a star) located at two different positions. If the star is not too compact, the effective spacetimes can be welldescribed within the Newtonian limit. For this setup, which is pictorially depicted in Fig. 2, Penrose is able to quantify the incompatibility of the superpositions. He defines a quantity ∆ which is basically the gravitational self-energy of the difference between the mass distributions of each of the two locations of the lump of mass: where ρ, ρ ′ are the mass densities of the two lumps. This quantity would control the lifetime τ of the superposition for these simple examples of cat-like states for the matter distribution.
For a careful analysis of the physical impact of this mechanism of gravitationally induced decoherence and a comparison with the more standard enviromental decoherence [16][17][18], see [60]. Moreover, recently it has been argued that natural models with τ = η /∆ and η of order 1 are falsified experimentally 2 : the relaxation times these models imply appear to be much shorter than the ones inferred in the experiment [59]. However, it is important to recall that the indirect results of this experiment crucially rely on assuming a random diffusive behaviour for the density matrix of the system. For us it is not clear to what extent this needs to be a compulsory assumption.
The mechanism causing decoherence in our model, although similar in spirit to Penrose's proposal differs from it in several points. In our toy model the instability under local perturbations that the cat state suffers comes from the huge fluctuations in the particle number density that the system displays. Recalling that the speed of sound at a point in the condensate c s is proportional to the local number density of particles n c , huge fluctuations in the number of particles would induce huge fluctuations of the underlying causal structure for the propagation of sound waves in the effective geometry. Such fluctuations arise because the cat state is highly delocalized. However, Penrose argues that the source of instability for the superpositions of spacetimes arises due to a purely geometric reason: it arises from the inability to make sense of the notion of stationarity for such states. The absence of a well-defined notion of time, i.e. the absence of a well-defined parameter with respect to which we can talk about evolution when considering these kind of superpositions makes impossible to regard the system as having some kind of causal structure. Although some fuzziness of the causal structure when entering the quantum regime is to be expected, attempting to superpose two radically different causalities would drive us completely into this regime where the causal structure itself dilutes and becomes a meaningless concept. Thus, the reason for the decay of a delocalized superposition seems to be the huge instability these states develop due to being far from having a well-defined causal structure.
All mechanisms that try to account for the reduction of state in quantum mechanics need to introduce parameters controlling how macroscopic a system is. Standard environmentally induced decoherence generically uses the number of particles in the system N as a measure of how macroscopic the system is. Penrose's model measures how complex and how macroscopic a system is in terms of the gravitational field it can create. In general terms, we can expect that for ordinary macroscopic matter both measures display similar behaviours regarding the decoherence time-scales, but conceptual differences appear. If one were able to monitor the whole set of particles composing a system, it would be possible in principle to distinguish among both mechanisms [60].
An important difference of our model with respect to Penrose's idea is that in BECs we can make sense of the superposition of two Minkoswki spacetimes with different speeds of sound. In Penrose's geometrical language these spacetimes are one and the same. If causality were an emergent property, much like the sound speed of collective excitations, then it would have additional properties beyond those associated with the effective geometrical description. Our model suggest that the inability of superposing different causalities is a robust idea that could survive in different approaches to quantum gravity. Although in a completely different setup, similar resonant conclusions are reached in Euclidean Dynamical Triangulations and Causal Dynamical Triangulations. In the former, one sums over arbitrary Euclidean manifolds and the result is dominated by configurations that do not resemble at all a smooth, almost flat on suitable scales, spacetime that we live in. On the other hand, in the latter one just sums over geometries that lead to a welldefined causal structure, hence constraining the geometries one sums over. In this case, it is possible to find almost macroscopic spacetimes [7].
Something quite similar happens when one tries to use analogue models for constructing effective spacetimes with causal pathological behaviours, for instance, with chronological horizons [22]. In the general relativistic description, the presence of causal pathologies appears to be tied up just to the presence of exotic matter: the semiclassical regime does not uncover further limitations. However, playing with analogue models suggests that one might find further fundamental limitations, like the existence of an underlying fundamental causality, making the difficulty of violating causality principles also a robust idea.
One additional difference between the two mechanisms is that, strictly speaking, in our formalism the loss of coherence is only effective as the underlying system is a Nbody quantum mechanical system. On the contrary, in Penrose's description, the gravitational decoherence appears as something fundamental. In a sense, our proposal seems to be closer to the environmental decoherence proposal than Penrose's. However, we rather think that our analysis highlights that the line separating an effective decoherence from a fundamental decoherence is blurred. Our toy model suggests that even taking a purely quantum theory as the high-energy theory, any practical and effective physical theory constructed by macroscopic observers would contain quantum elements but also some elements that should not be treated quantum mechanically: there are elements in the effective theory that are collective excitations of already quantized degrees of freedom and they must not be quantized again. The difference with standard environmental decoherence is that the environmental degrees of freedom in this case are not so controllable from the effective theory point of view.
As a final comment, let us discuss the interplay of our model and the problems that quantum mechanics faces, since they were Penrose's original motivation to propose his ideas. Quantum mechanics, as described by just a dynamical wavefunction, faces two related problems. The first one is how to accommodate the apparent absence of quantum superpositions of macroscopic systems. As the formalism per se does not suggest a division between microscopic and macroscopic systems, the same rules for superposition should apply to both realms. The second problem, is how to account for the specific although probabilistic results of actual acts of measurement in the laboratories, i.e., how to account for the apparent intrinsically statistical result of individual acts of measurement. Without the additional Born rule, the unitary evolution of quantum mechanics does not even describes what is a measurement nor what is the source of non-determinism. As a disclaimer, let us recall here that a hidden-variable take on quantum mechanics such as that of Bohm provides a solution to these problems at the cost of introduc-ing some explicit non-locality in its formulation [61,62]. As explained in the introduction, our discussion in this paper is focused on interpretations of quantum mechanics giving some reality just to the wavefunction.
Our model, like Penrose's proposal, directly accounts for the first of these problems: how a superposition ends up becoming a mixture (a state lacking quantum coherences). For the second problem, it does not offer a complete solution since it would require making an explicit model of measurement. Actually, this is a longstanding problem that has accumulated a vast literature [63]. Nonetheless, we feel that our analogue toy model adds some interesting new perspectives to this second conundrum.
The philosophy underlying the analogue gravity program and in a more complete sense, emergent gravity scenarios (i.e. hypothetical scenarios in which Einstein equations are also emergent [64]) is that everything that we perceive is constituted by low-energy collective excitations of some fundamental microscopic degrees of freedom. For instance, the spacetime itself and the quantum fields describing the Standard Model of Particle Physics would be part of such collective excitations. In that sense, although the microscopic theory may not display a lorentzian causality, the effective low-energy regime does. Within this logic, observers would also be "emergent" in the sense of being composed of excitations of these effective fields. The notion of observer that we are using is that of a conscious being, understood not literally as a human being but as any primitive form of consciousness or even a recording device. In general any sufficiently complex system relying on causality does the job. In fact a crucial prerequisite for consciousness to appear is causality itself 3 . Now, recall that the main conclusion that we have drawn up from our analysis is that causality is a special macroscopic property in the sense that it is not amenable to be in quantum superposition. Then, an immediate consequence of this tight relation between causality and consciousness seems to be that different states of consciousness will also tend to rapidly decohere.
Our construction in this paper, takes standard nonrelativistic quantum mechanics, call it NRQM1 as the microscopic fundamental theory with no special commitment to any of its interpretations. Then, inner observers living in the system could make another non-relativistic quantum theory (NRQM2) to describe the phenomena they perceive whenever the regimes they are probing are non-relativistic. NRQM2 should of course become a Lorentz-invariant quantum field theory at high energies from the internal perspective but still low energies from the perspective of laboratory observers, that would still use NRQM1. Now, consider for example that NRQM1 as a high-energy theory happens to provide only practical decoherence and nothing similar to a collapse to one of the reduced states. Even in this case the description of NRQM2 will tend to incorporate in its rules a collapse postulate. This will be so because conscious beings will never have fuzzy experiences associated with being in cat-like states, i.e. they will never perceived macroscopic quantum coherences between states with very different causalities. If on the contrary the formulation of NRQM1 incorporated a mechanism for the selection of one and only one alternative under individual acts of measurement (by incorporating for instance a stochastic rule, for exampleà la Ghirardi-Rimini-Weber [36]), the NRQM2 would have the same form as before. Thus, the phenomenology of NRQM2 is independent of the measurement theory taken in NRQM1. In this way our example illustrates how difficult it is to learn things about a microscopic fundamental theory if we only have access to the collective excitations belonging to its low-energy sector.
This is closely related to Wigner's discussion on the role of consciousness in our physical description [31]. Essentially, he advocated that sufficiently macroscopic or complex processes (not just those involving conscious beings as complex as human beings) will in practice decohere, marking a separation between the classical and quantum realms. Notice that here, complexity and macroscopic are used in a non-technical sense referring to an imprecise characterization of such magnitudes. The difference of our approach is that it highlights the decoherent character of causality as a more primitive notion than consciousness.

V. SUMMARY AND DISCUSSION
Our aim in this work has been to present an analogue model of gravity where it is possible to analyze the properties of the states corresponding to would-be superposition of spacetimes. Although they are logically possible in theories of quantum gravity, we find difficulties when we try to engineer them in analogue gravity set-ups.
We have studied a toy model of a Bose-Einstein condensate in a double well which in a suitable regime gives rise to what we could call a superposition of spacetimes. This toy model shows three phases: a coherent phase, when interactions are negligible; the fragmented uncorrelated phase, whose representative in the regime in which the hopping dominates is the Fock state; and the fragmented correlated phase, whose representative when the interactions dominate is the so-called cat state. In principle, the coherent state cannot lead to causal behaviour of signals propagating on top of it. Furthermore, even if we turn on some repulsive interactions for the system to display such behaviour, the resulting state does not resemble at all what we would expect that a superposition of spacetimes looks like. The Fock state lacks quantum coherences and we regard it as a mixture of spacetimes, instead of a quantum superposition. It is stable under local perturbations and hence offers no difficulty either from the condensed matter perspective or from the point of view of its analogue gravity interpretation. The cat state, which we would interpret as an analogue of the superposition of spacetimes, is highly unstable under local perturbations, that tend to drive it to a Fock-like state. This is exactly the same phenomenology advocated by Penrose. Superposed spacetimes should rapidly decay into mixtures of spacetimes. The rate at which this decay occurs is proportional to a measure of how macroscopic the system is. Penrose measure is the magnitude of the gravitational field that the system is able to generate. In our model, this measure is given by the number of condensed particles and the ratio between the hopping parameter and the interactions (t/U ). Penrose proposal is based on the fact that a superposition of two spacetimes is not a well-posed concept, since it does not exhibit any sensible causal structure. This happens even if one superposes stationary causal structures. In our model, the fluctuations in the number of particles on each well is what quantifies the instability of the BEC under local perturbations. We recall that the sound speed is proportional to the local number density of bosons, and, thus, huge fluctuations in the number density are to be interpreted as a huge undefiniteness of the effective causal structure of the system. In that sense, we find that our analogue model illustrates the difficulty in pushing the causal structure of analogue systems to its limits: systems avoid causal pathologies. This is similar to our previous findings that analogue systems avoid giving rise to chronologically illbehaved spacetimes [22]. To put it explictly again, we have found that analogue systems tend to avoid configurations exhibiting superpositions of analogue spacetimes (or causalities); this means, configurations that do not have a well-defined causality.
Thus, in tune with Penrose's ideas, the lesson we seem to be learning is that, although macroscopic superpositions might be relevant for very high energies (this situation is of course logically possible), there are arguments to suspect that they might be qualitatively irrelevant at low energies at which gravity could be described by an effective geometry.

ACKNOWLEDGMENTS
GGM would like to thank Alfredo Luis for very useful discussions.
CB and GGM acknowledge financial support from the State Agency for Research of the Spanish MCIU through the "Center of Excellence Severo Ochoa" award to the Instituto de Astrofísica de Andalucía (SEV-2017-0709). GGM is funded by the Spanish Government fellowship FPU20/01684.