Back-reaction in canonical analogue black holes

We study the back-reaction associated with Hawking evaporation of an acoustic canonical analogue black hole in a Bose-Einstein condensate. We show that the emission of Hawking radiation induces a local back-reaction on the condensate, perturbing it in the near-horizon region, and a global back-reaction in the density distribution of the atoms. We discuss how these results produce useful insights into the process of black hole evaporation and its compatibility with a unitary evolution.


Introduction
Hawking's prediction of black hole radiation in 1974 [1,2] has been a milestone in our understanding of gravitation, by linking it to quantum field theory, thermodynamics and information theory, with far reaching consequences till nowadays.
The presence of a horizon has two important consequences. First of all it determines the causal separation of a region of spacetime, containing information that is no longer accessible by observers outside. Moreover, it induces the phenomenon of emission of black-body radiation which has the effect of extracting mass from the black hole leading to evaporation, and therefore to the apparent loss of the information contained inside that region.
The production of radiation can be predicted through the tools of quantum field theory in curved spacetime [3], but the description of the evaporation process requires a consistent formulation of semiclassical gravity, in which the degrees of freedom of geometry interact with the quantum degrees of freedom of matter fields. From a geometric point of view, this process entails the shrinking of the radius of the horizon and therefore the raise of the radiation temperature.
Furthermore, Hawking radiation introduces the so called transplanckian problem [3], i.e. the apparent necessity for Hawking's derivation of strong assumptions about the structure of spacetime and of the quantum field theory beyond the Planck scale. The theory of general relativity is not valid in the high energy regime that is required to fully define how geometry is related to the degrees of freedom of quantum fields and how it is influenced by high-energy wave packets. Settling these issues seems then to require a knowledge of the ultraviolet completion of the theory, so rising the question about the robustness of Hawking radiation with respect to the details of the latter.
Unfortunately, we presently do not have the ability of settling these issues experimentally in gravitational systems, as measuring Hawking radiation is basically impossible for astrophysical objects. However, Unruh's realisation that hydrodynamical systems can be used for simulating black hole spacetimes [4,5], paved the way for devising more general table top experiments in which the propagation of some perturbations on a given background can simulate classical and/or quantum field theory on a curved spacetime. This realisation marked the dawn of the so called Analogue Gravity research field [6].
In particular, in these analogue models the correspondence of their mathematical equations with those of black holes, opens the possibility to reproduce in a lab the analogue of Hawking radiation in the presence of a horizon, and so to start addressing long standing issues like the aforementioned ones. Indeed, in these systems the robustness of Hawking radiation with respect to a well understood ultraviolet completion of the effective spacetime and quantum field theory was soon argued for on theoretical grounds [7,8,9], further stimulating the attempt to achieve a direct observation.
Bose-Einstein condensates have been among the most successful analogue models due to their intrinsic quantum nature, simplicity and experimental realizability [10,11,12,6]. In these systems, Bogoliubov quasi-particles are quantum excitations that propagate following the acoustic metric induced by the wavefunction of the condensate, and the atomic structure provides a natural regularization of the theory at short distances. As such, they have provided an ideal setting for the study of analogue Hawking radiation and cosmological particle production [13,14,15,16,17,18,19,20,21,22,23,24]. Even more remarkably, strong experimental evidence for Hawking radiation in BEC analogue systems was recently presented [25,26,27,28,29], largely confirming the theoretical expectations (at least in some suitable regime of the analogue black hole dynamics) and quite convincingly settling the issue of the robustness against the ultraviolet completion of spacetime.
These results have of course stimulated new studies, also with regard to the interplay of the microscopical (the analogue quantum gravitational) dynamics and the emergent phenomenology of spacetime and quantum fields. In particular, it was argued by the present authors that the other big issue related to Hawking radiation, the information loss problem -the apparent loss of unitary evolution implied by this phenomenoncan be understood using again the analogue gravity insight [30]. Indeed, in a isolated BEC system unitarity must be preserved and it is so once the evaporation process is considered at the microscopical level, tracking down the generation of correlations between the atoms underlying the emergent spacetime and the Hawking quanta.
Although in analogue gravity the back-reaction of the quantum fields on the acoustic geometry would not follow the Einstein equations, the observation of evaporation in analogue black holes could provide information over the interplay of quantum fields and classical geometry within a semiclassical scheme in a broad class of quantum gravity scenarios. If one sees this analogue gravity system as a toy model for emergent gravity, it could then give a physical intuition of how Hawking radiation can be pictured as a feature emerging from an underlying full quantum theory.
So, seeking this kind of theoretical insights, as well as the possibility to test this understanding of the back-reaction in experimental realisations, attracted an increasing interest on the nature of the back-reaction in analogue Hawking radiation [31,30,32]. For example, in recent experiments with surface waves on a draining vortex water flow [32], it has been shown that the back-reaction is indeed observable, and it is possible to measure the exchange of energy and angular momentum between the background flow and the waves incident on the horizon. In Bose-Einstein condensates, a next step in this direction should be achievable as it should be possible to observe the back-reaction of the Hawking Bogoliubov quasi-particles on the acoustic geometryi.e. on the substratum of condensate atoms -and interpret it through the lens of the underlying atomic theory.
With this aim in mind in this paper we focus our attention on the back-reaction of the Hawking radiation on a canonical analogue black hole within a BEC analogue system characterised by a λφ 4 interaction [33,34,35]. The canonical analogue black hole is for our study a remarkable geometry because it is the spherically symmetric stationary solution of the Gross-Pitaevskii equation with homogeneous atom number density. The spherical symmetry allows to reduce the problem to a time-radius problem that presents a horizon at the radius where the velocity of the stationary ingoing superfluid equals the speed of sound.
In what follows we first describe the acoustic metric in this system and how the solutions of the Klein-Gordon equation behave at the horizon, and we calculate the Hawking radiation. These are the basis for showing how the back-reaction exerted by the Hawking radiation can be studied and how it affects the geometry of the black hole. We then propose to study a regime in which the Hawking radiation would lead to a sizeable evaporation of the analogue black hole. Finally, we show how this evaporation corresponds, in the evolution of the system, to the depletion of atoms from the condensate to the non-condensed part.

Analogue gravity in Bose-Einstein condensates
We start by giving a brief overview of analogue gravity with Bose-Einstein condensates. This fixes the notation and provides results used in the following Sections.
Consider the local λφ 4 theory for a non-relativistic complex scalar field ( ≡ 1) in 3+1 dimensions, with the usual bosonic commutation relations. The dynamics of the field is described by the Schrödinger field equation Here and in the following, unless explicitly stated otherwise, we assume all the operators in a given correlation function to be evaluated at same time and same position.
In presence of Bose-Einstein condensation we can describe the field by splitting it in mean-field and quantum fluctuations as [δφ (x) , δφ (y)] = 0 .
At leading order the mean-field is described by the solution of the Gross-Pitaevskii equation [36] i∂ The linearized quantum fluctuation follows the Bogoliubov-de Gennes equations, which are coupled to the mean-field equation. When the mean-field is described by the Gross-Pitaevskii equation, the linearized quantum fluctuation can be solved separately after the solution of the Gross-Pitaevskii equation is found: If instead, we include the back-reaction of the quantum fluctuation in the dynamics of the condensate wavefunction, the equations for the mean-field and the quantum fluctuation cannot be separated. The modified Gross-Pitaevskii equation includes the anomalous terms n and m [37] and the Bogoliubov-de Gennes equations are consequently modified, thus depending on this better approximation of the mean-field instead of the previous expression These expressions will be used in section 5.
With the usual Madelung representation introducing number and phase variables, the real and imaginary part of the Gross-Pitaevskii equation can be reorganised in the two quantum Euler equations for the condensate We can also apply the Madelung representation to the quantum fluctuations and rewrite the Bogoliubov-de Gennes equations obtaining two coupled dynamical equations for θ 1 and ρ 1 from which we see that these two real quantum fields are one the conjugate of the other. If we consider the quantum pressure to be negligible we obtain Eq. (22), which provides easily workable linearized expressions: when the spatial derivatives of ρ 1 are negligible, ρ 1 is proportional to the flux derivative of the field θ 1 . The Bogoliubov equations can be combined into a second order partial differential equation for the real scalar field θ 1 , which appears to be a massless scalar field propagating in a curved spacetime. The resulting equation takes the form of a Klein-Gordon equation for an acoustic metric, analogue to that of a curved spacetime. One gets From this equation we can directly extract the elements of the inverse of the metric tensor. We do so rewriting the coupling λ in terms of the speed of sound c s = λ ρ 0 /m of the condensate: which correspond to the metric tensor of components This acoustic metric is analogous to the metric of a 3+1 spacetime in Painlevé-Gullstrand coordinates up to the conformal factor ρ 0 /c s . The role of the velocity V is played by the superfluid velocity In the case of spherical symmetry with ∂ r θ 0 < 0 we are in the case of ingoing Painlevé-Gullstrand coordinates, which allow to study horizon penetrating trajectories and the behaviour outside a black hole horizon. The solutions for the null radial geodesics are one regular ingoing solution and one outgoing solution that peels off the horizon (at the radius where |V | = 1) In the following we will make use of the useful notation since it is most convenient to describe the behaviour at the horizon.

Canonical analogue black holes
Among the acoustic metrics simulated with Bose-Einstein condensates, it is convenient for our investigation to consider the canonical analogue black hole. A canonical analogue black hole could be obtained as the stationary solution of the Gross-Pitaevskii equation with spherical symmetry and with homogeneous atom number density. Due to the spherical symmetry, it presents a horizon at the radius where the velocity of the ingoing flow is equal to the speed of sound. The stationarity of the background metric, of the mean-field, means that the deviations from this configuration should be understood as the effect of the phenomenology previously neglected: the subleading effect of the quasi-particles and their back-reaction. The homogeneity of the atom density ρ 0 means that the conformal factor ρ 0 /c s is a number, and can trivially be transformed away.
Taking a constant and homogeneous number density ρ 0 , and assuming θ 0 to be a function of the radius only, from the Eqs. (13)- (14) we obtain the stationary spherically symmetric solution of the Gross-Pitaevskii equation. The superfluid velocity and the external potential setting the corresponding velocity profile are respectively with an ingoing flux of atoms crossing every closed surface containing the origin at a rate of 4πc s r 2 H ρ 0 atoms per unit time, which must be pumped in the system and move towards the singularity in the origin. The radius of the horizon r H is set by the external potential V ext .
Omitting the conformal factor, the analogue metric implied by Eq. (36) is therefore The singularity in the origin poses a problem that would be unavoidable in devising an experiment: one should either remove the atoms reaching the origin and put them back in the system at large radius [38], or should consider only the region in which the hydrodynamical limit holds with a good approximation, at a distance from the origin and for a period of time such that it is not affected by atoms which flowed towards the origin. 1 These are however practical concerns that we shall assume dealt with in our theoretical study. In addition to them, we observe that the possible removal of atoms from the system would affect the atomic correlation functions and introduce statistical mixture, but the dynamics of the quasi-particles in the near-horizon region, on which we focus in the rest of the paper, would not be sensibly affected.

Klein-Gordon equation and field modes
The massless scalar field living in our analogue spacetime can be described in terms of its modes, found by solving the Klein-Gordon equation. For a stationary spherically symmetric spacetime in Painlevé-Gullstrand coordinates, the equation to solve is where Ω is a compact notation for the solid angle and L 2 is the usual angular momentum operator We can find a complete set of solutions looking for those which preserve the symmetries of the system [34,35], i.e. the eigenfunctions of both the time-derivative and the angular momentum operator where the functions Y lm (Ω) above are the usual spherical harmonics. So Eq. (39) becomes We have written Eq. (45) reorganising the various contribution to have a clearer view of the problem. The operator h∂ r can be understood as the directional derivative with respect to the tortoise radial coordinate. The term h∂ r h/r is a potential which induces mode mixing. This second order differential equation will have two linearly independent solutions for each set of eigenvalues, an outgoing mode and an ingoing mode. Splitting the solutions in norm and phase we get two coupled real differential equations Different integrating constants C set the sign of the derivative of the phase and give the different solutions.

Near-horizon behaviour
We are interested in the solutions in proximity of the horizon, for a small distance δr = r − r H , and such that n ωl (r) is regular. The horizon-crossing solutions will have a regular behaviour in phase, and the horizon-tangent solutions will be the remaining linearly independent solutions. In presence of a horizon, h vanishes at r H and so do all the terms depending on it in the Klein-Gordon equation. We can therefore set the phase-defining constant C of the solutions with regular norm by requiring it to satisfy the Klein-Gordon equation at the horizon Considering the phase of the function f ωl (r) we get, in the specific case of the canonical analogue black hole, When the constant C takes positive values, the above expression describes solutions that peel off from the horizon and are relevant for the production of Hawking quanta, while for negative values of C it describes horizon-crossing ingoing solutions. From Eq. (53) it follows that the phase of each mode with positive C diverges logarithmically, and on the two sides of the horizon we have two independent solutions. In the outer region these are outgoing modes, while in the inner region they peel off the horizon towards the singularity. Considering the radial part of the modes, we have defined them either for r > r H or for 0 < r < r H . These functions can be analytically continued in the complex plane only assuming the presence of a branch cut.
In conclusion, the various modes are the horizon crossing modes f HC ωlm , the outgoing modes in the outer region f ext ωlm and their counterparts in the inner region f int ωlm .

Far-field behaviour
At large radii the equation Eq. (49) is such that h ≈ 1 and the mode-mixing potential h∂ r h/r is negligible, and the modes can be approximated with a linear combination of Bessel functions of the first and of the second kind. At infinity even the angular momentum potential can be neglected and the properly normalized solutions take the form where ingoing and outgoing modes are mixed at infinity, as expressed by the phase Θ ωl due to the potential encountered along the radial propagation. In this region, apart from the overall phase, which is given by the limit lim r→∞ r √ 1−h h dr, and by the conditions at the horizon, the radial part of the outgoing and the ingoing modes are one the conjugate of the other, with

Mode decomposition of the scalar field
In conclusion one obtains the mode decomposition of the scalar field, as a superposition of modes with fixed angular momentum and frequency Depending on the radius, only one set of modes between f ext and f int appears in the definition of the field. The field operators follow the canonical commutation relations and the ladder operators relative to the modes are defined to follow the usual commutation relations in the quasi-particle Fock space and all the other commutation relations between ladder operators vanish identically.

Hawking radiation
We discuss now the creation of Hawking quanta. The presence of the horizon induces the creation of Hawking radiation, and the velocity profile defines the temperature of the canonical analogue black hole. Then we study the properties of the initial state, the Unruh vacuum state. Until now we have focused on the region outside the horizon, where for every positive ω the modes f ext ωlm and f HC ωlm -with positive norm with respect to the Klein-Gordon internal product -must have the time dependence e −iωt usual for stationary spacetimes, while in the inner region the time dependence is inverted for the modes with positive norm f int ωlm , because of the change of sign of the function h (r). Moreover, each of the modes peeling off from the horizon has a phase that diverges logarithmically, as can be deduced from Eq. (53), implying that on the two sides of the horizon the two independent sets of modes f ext ωlm and f int ωlm can be analytically continued and put in superposition with eachother.
The operation of mixing destruction operators of the modes defined in the outer region with creation operators of the modes defined in the inner region is described by a Bogoliubov transformation.
In particular we are interested in the mixing between modes at the horizon with those associated to a vacuum state for static observers at past null infinity.
Near the horizon, at the two sides, the modes are (dropping the eigenvalues from the notation) Linear combinations of these modes define new solutions, in particular we are interested in the analytic continuations of the radial part on the real r-axis. One solution can be extended in the upper half-plane -with a branch cut in the lower, and one in the lower half-plane -with a branch cut in the upper.
Due to the logarithmic term in the phase of the functions, passing from one side of the horizon to the other, the analytic continuations gain a phase of ±πωr H /2c s , and therefore we can write The coefficients α and β of the Bogoliubov transformation can be found imposing that the field can be rewritten in terms of these new modes and imposing the canonical commutation relations for the new ladder operators: Therefore it is found that, apart from irrelevant overall phases, the Bogoliubov transformation is given by The Unruh vacuum |∅ is defined such that so Eqs. (72)-(73) imply that an observer outside the horizon will observe a thermal distribution of radiation coming from the black hole horizon. Indeed, through the Bogoliubov transformation, the occupation of the modes f ext is given by This occupation number corresponds to the thermal Bose-Einstein distribution associated to the temperature where κ is the surface gravity of the canonical analogue black hole, which is half the radial derivative of h evaluated at the horizon

Back-reaction
When the Gross-Pitaevskii equation is modified by introducing the terms of anomalous mass and density, it becomes possible to study the back-reaction that the analogous scalar field induces on the condensate wavefunction.
The study of the back-reaction in analogue systems enables to address various issues, such as the problem of the evaporation of black holes and in general the dynamics of analogue horizon. To study the dynamics of the horizon it is necessary to express the anomalous mass and density as correlation functions of the analogous scalar field and its conjugate. While the anomalous density is strictly real, the anomalous mass contains a real and an imaginary part.
It is particularly convenient to study these quantities in dimensionless terms: We observe that in Eq. (83), in the hydrodynamic limit, the imaginary part of the anomalous density is the derivative along the flow of the θ 1 θ 1 correlation function, the vacuum polarization. Note that in Eq. (81) we have regularized the definition of anomalous density, thus eliminating a δ 3 (0) term. Formally this is done defining the correlation function as the limit for vanishing distance of the 2-point correlation function.
To describe the mean-field dynamics with the inclusion of the anomalous terms we modify the Gross-Pitaveskii equation as in Eq. (9) and apply the Madelung representation. 2 We study separately the square norm of the wavefunction -the atom density of the condensate -and the phase: The equation for the dynamics of the atom density function is affected by the imaginary part of the anomalous mass, while the equation for the dynamics of the phase of the condensate wavefunction is modified by a combination of anomalous density and real part of the anomalous mass.
We study the regime in which these classical perturbations to the condensate wavefunction are small deviations from the stationary solution of the Gross-Pitaevkii equation that do not break the hydrodynamic approximation.
Therefore, we neglect the terms that are suppressed by the healing length scale ξ = 1/ λ ρ m, because they are only relevant at a high energy scale where the interactions between the individual atoms become dominant in determining the phenomenology. In Eq. (84) and Eq. (85) we neglect the contributions to the anomalous density and to the real part of the anomalous mass which goes as ρ 1 ρ 1 . Since the definition of the conjugate field ρ 1 depends on the healing length scale as the anomalous density in Eq. (81) and the real part of the anomalous mass in Eq. (82) can be approximated to the contributions of the correlation function θ 1 θ 1 : Under these assumptions, from Eq. (84) and Eq. (85) we proceed linearizing the equations for the mean-field for small perturbations -of the same order of magnitude of the anomalous terms -of the solution of the Gross-Pitaevskii equation, getting We can see that the classical fluctuation of the phase of the condensate δθ evolves like a Klein-Gordon mode, with conjugate momentum δρ ρ 0 + θ 1 θ 1 .
These equations allow to study the dynamics of the classical perturbation, and are differential equations that can be solved after the source term θ 1 θ 1 and the initial conditions for δρ and δθ are provided. Unless broken by the initial conditions, the evolution of the classical perturbation will preserve the spherical symmetry and the time translation symmetry.

Dynamics of the horizon
The production of Hawking radiation emitted by a black hole formed by gravitational collapse is a process that transfers energy towards the future null infinity. The outgoing flux of Hawking quanta induces a loss of mass of the black hole, whose apparent horizon gradually shrinks, until its eventual evaporation. In this non-stationary process, the gravitational collapse triggers the production of Hawking radiation quanta, consequently exerting a back-reaction on the classical metric tensor.
In the case of the canonical analogue black hole we can study the dynamics of the horizon assuming the condensate to be initially in the stationary condensate configuration described in section 3. Then, turning on the production of quasi-particles of the field θ 1 , we observe a back-reaction on the wavefunction.
In general, having spherical symmetry, the horizon radius r H is defined as the radius at which the radial velocity of the condensate equals the speed of sound. In the unperturbed case we have If the condensate is in a state which is a small perturbation of the stationary solution of the Gross-Pitaevskii equation, the density and the phase of the condensate are modified as described in Eqs. (88)-(89). Consequently the velocity of the condensate and the speed of sound change and, under this variation, Eq. (92) is satisfied at a different radius r H + δr H , where r H is the unperturbed horizon radius: The denominator in Eq. (95) is different for each spherically symmetric system and is proportional to the surface gravity of the horizon, i.e. to the temperature of the black hole in Eq. (78) The radius of the perturbed horizon is smaller or larger than r H depending on the sign of the variation δ (V 2 /c 2 s ).
For the canonical analogue black hole we substitute the expression of V 0 of Eq. (36) in Eq. (95) and we get Deriving the variation of the horizon radius with respect to time we observe which quantities determine the dynamics of δr H For a stationary spacetime the second term vanishes. We consider the case of the canonical analogue black hole substituting from Eq. (90) and Eq. (91) the expressions for δV and δρ : The source term on the RHS is a quantity that depends on the quantum state of the field θ 1 -through the correlation function θ 1 θ 1 -and on the state of the classical perturbation. It does not have definite sign: different initial conditions for the classical perturbation determine different regimes of the black hole, i.e expansion or contraction. This is expected, as Eq. (99) applies to every possible regime of the black hole. Nonetheless, such different initial conditions for the classical perturbation do not affect, the contribution from the quantum field θ 1 , which remains always the same. In the hierarchy of the equations, the analogue scalar field affects the perturbation of the analogue metric, not viceversa. A self-consistent semiclassical approach would require including the perturbation of the wavefunction in the Bogoliubov-de Gennes equations, explicitly or by iteratively recalculating θ 1 θ 1 and the pair ( δρ , δθ ).
The effect of the analogue Hawking radiation is understood when θ 1 θ 1 is the contribution leading the phenomenology. We therefore assume the perturbation terms δθ and δρ and their radial derivatives to be negligible. In the analogy between the canonical analogue black hole and the gravitational black hole, this corresponds to assuming the initial perturbation of the metric to be negligible, and then to excite it through the presence of Hawking radiation. This reflects the evolution of a gravitational black hole, where the free falling matter induces the formation of the horizon, which triggers the production of Hawking quanta slowly driving the spacetime far from the stationary configuration.
Therefore we consider Eq. (99) and assume the perturbation of the wavefunction to be null and the corresponding terms to vanish. Moreover, we can assume the time derivative of the Hawking radiation term to be negligible. The Hawking quanta are produced in pairs of equal frequency, and in the correlation functions their time dependent phases cancel each other out. The term ∂ t θ 1 θ 1 is therefore suppressed.
In this regime, we are left with The quantity θ 1 θ 1 is of paramount interest in the investigation of the back-reaction that quantum fields exert on curved geometries. Together with the stress-energy tensor, the vacuum polarization is a quantity that in literature is studied in association with the production of Hawking quanta in gravitational systems [39,40].
As described in [41,42] this function is obtained through the coincidence limit of the Green function of the Klein-Gordon operator, i.e. the limit for x ′ → x and t ′ → t of G (t, x; t ′ , x ′ ). Different boundary conditions determine different expressions for θ 1 θ 1 : the occupation numbers of the modes of the field, which are associated at the horizon with the Hawking radiation and with the horizon-crossing quanta, determine the behaviour of this function. In particular, the occupation numbers are zero for the Boulware vacuum; are thermal for the outgoing modes and null for the ingoing modes for the Unruh vacuum (as seen in Eqs. (74)-(75)-(76)); and are all thermal for the Hartle-Hawking vacuum. To obtain the expression for θ 1 θ 1 in the various cases one must subtract different renormalization counterterms. Renormalizing is strictly necessary as the Green function generally includes divergent terms. In the coincidence limit already in the 3+1 flat Minkowski spacetime it is found that The most commonly used renormalization scheme is based on the point splitting regularization method [43,44,41] which removes the divergence by splitting the point in which the Green functions is evaluated in two nearby by points characterised by their geodesic distance, so regularizing the vacuum polarization as measured along geodesic trajectories. For Hadamard states [45] the resulting divergent structure in the coincidence limit is universal in curved spacetimes: one gets the above mentioned divergent terms together with other logarithmically divergent terms typical of the Hadamard structure. The universality of such ultraviolet divergences (of the same functional form of those obtainable in Minkowski spacetime), allows to safely discard them. However, other irregular behaviours may nonetheless arise from the peculiarities of the curved geometry and the vacuum state. In particular, in presence of a horizon, the Boulware vacuum gives a divergent vacuum polarization due to the vanishing of the time-time element of the metric h. As argued by Candelas, in proximity of the horizon the Green function is divergent as G = −1/hǫ 2 [42]. Also, let us mention that while the above mentioned regularization schemes have been mostly applied in Ricci flat spacetimes, in non-Ricci flat spacetimes like ours, they will generically include an extra contribution, which however we can expect to provide at the horizon at most a prefactor of order unity, which can be therefore neglected for our considerations.
When we consider the Bogoliubov transformations, the Green function depends on the occupation numbers of the quasi-particle states. For states different from the Boulware vacuum, in presence of non-null occupation numbers, the vacuum polarization is given by the limit where the first contribution is fixed for every state.
To appropriately regularize the contribution of the Hawking radiation, it is therefore necessary to subtract the quadratic divergence that defines the Boulware vacuum and consider the occupation numbers n ω . The Unruh vacuum is associated with the occupation numbers of Eq. (77), for which the integral expression in Eq. (102) gives This quantity is divergent at the horizon, but is regularized by subtraction of the second term in the series expansion of the geodesic distance c 2 s 6h h ′ 2 2 , which removes the divergence due to h in the limit r → r H . Reintroducing the proper numerical factors, we therefore obtain that in proximity of the horizon the leading effect of the Hawking radiation induces the perturbative expression This approximation does not hold for every radius, e.g. it includes not integrable divergences in the origin, but is such that it is consistent with the expected behaviour of the Rindler spacetime -the near-horizon region [46,47]. In Eq. (104) we retain the term of order δr, as we are interested in the radial derivative of this expression.
Note that the adimensionality of the field is a feature that differs from the usual notation, as in 3 + 1 quantum field theory the Klein-Gordon scalar field has generally the dimension of an inverse length. In Eq. (104) we introduced the length scale λ/c s , which is the proper quantity needed to replicate the usual relations [θ (x) , π (y)] = −iδ 3 (x, y) .
Considering the definition of ρ 1 given in Eq. (22), a straightforward calculation gives the length scale required. Therefore we obtain the estimate for the time derivative of the horizon radius, and from the expression of the surface gravity we can also make an estimate of its rate of These derivatives are very small: while typically the atom separation -equal to ρ −1/3 0 -is of an order comparable to the healing lenght, the radius of the horizon can safely be assumed being much larger than both.
Given that the rate of change of the surface gravity is very small, the system can be assumed to evolve adiabatically along the evolution: the hydrodynamic approximation is broken before the rate of change of the curvature in Eq. (111) becomes comparable to 1.
We can therefore take Eq. (109) to provide a rough approximate description of the evolution of the black hole, promoting r H on the RHS to be the dynamical horizon radius. We obtaining that the expected lifetime of the black hole is long and proportional to the inverse of the interaction coupling λ Let us note that the cubic dependence on r H in Eq. (113) resembles that of a Schwarzschild black hole, and follows from the proportionality between δṙ H and r −2 H . To make this prediction it is necessary to assume the regime in which the Hawking radiation, i.e. the contribution of the Bogoliubov quasi-particles obtained from the solution of the Bogoliubov-de Gennes equations, is dominant over the pertubation obtained by backreaction in the modified Gross-Pitaevskii equation.
We briefly point out that another regime which could be of interest is that of an analogue black hole in equilibrium with the analogue scalar field, where the back-reaction is not neglected but the classical perturbation is stationary. This would be the analogue of the Hartle-Hawking vacuum state for our system. In this case a solution can be found without changing the speed profile, i.e. with δV ≈ 0, but with a perturbation of the number density δρ such that it counterbalances θ 1 θ 1 : In this case the horizon for the stationary configuration is found at a radius larger than r H , but this solution will be driven out of equilibrium by the terms of order ξ 2 and the non-linarities previously neglected. In this case the density-density correlation function differs from the unperturbed case by a negative quantity on the diagonal (when the density operators are evaluated at the same position): If in a realization of the system the evolution is kept under control allowing only adiabatic transformations between near equilibrium configurations, it would be reasonable to assume these initial conditions for the system; but this equilibrium configuaration would not be analogous to the Unruh vacuum.

Dynamics of the number of atoms in the condensate
From Eq. (90) we observe that the dynamics of the perturbation of the condensate affects not only the position of the horizon, but also the global properties of the system, such as the number of atoms in the condensate state. Let us consider the dynamics of δρ . The integral of this quantity in space gives the rate at which the atoms leave the condensate and move to the excited part. It is a global process which can already be described at the level of the Bogoliubov-de Gennes equations, considering the derivative with respect to time of the total number of atoms in the condensate 1-particle state, and assuming the stationarity of the unperturbed condensate wavefunction: The same derivative can be expressed as the derivative of the integral of the number density, i.e. of the squared norm of the condensate wavefunction. Since the unperturbed configuration is assumed stationary, the only contribution is made by the classical perturbation, i.e. by the back-reaction of the analogue quantum field. From Eq. (90) we get We observe that the time derivative of the number of atoms in the condensatein the entire occupied volume -depends on three factors: the expectation number of condensate atoms in the region within the horizon the Hawking temperature of the black hole, as defined in Eq. (78), that goes with the inverse of the horizon radius (plus subleading corrections describing the imperfection in the adiabaticity of the evolution); the value that the vacuum polarization θ 1 θ 1 assumes at the origin, which is effectively a dimensionless structure factor depending on the velocity profile. The vacuum polarization requires a proper renormalization and it is generally difficult to extend it down to the singularity in the origin, not only when studying the canonical analogue black hole but also, e.g. for the Schwarzschild black hole [48]. Anyhow, in the integral form as in Eq. (120) it is clear that it is the radial derivative of the vacuum polarization, on both sides of the horizon, that determines locally the contribution to the depletion.

Discussion and conclusions
We have studied the canonical analogue black hole as a particular realization of analogue gravity with Bose-Einstein condensates.
The linearized quantum fluctuation over the background -the mean-field described by the Gross-Pitaevskii equation -propagates as a massless scalar quantum field in the curved spacetime of a black hole. In the limit of negligible quantum pressure the Bogoliubov-de Gennes equations can be reorganised to show how the quantum fluctuation of the phase is governed by an acoustic metric, defined by the condensate wavefunction, that presents a horizon. The calculation of the Hawking radiation produced in this acoustic geometry follows closely that of a massless scalar field near a gravitational black hole, and we have set the conditions for the modes of the field to replicate the Unruh vacuum, the state expected in a black hole spacetime formed through gravitational collapse.
Being interested in the back-reaction that the quasi-particles exert on the condensate, we have modified the Gross-Pitaevskii equation to include the anomalous terms, and thus observed how the acoustic metric is modified by the quasi-particle dynamics. The study of the quasi-particle back-reaction pushes analogue gravity towards an understanding of semiclassical gravity, where the back-reaction of the quantum fields is expected to be included in Einstein's equations.
We have provided expressions for the back-reaction in analogue gravity with Bose-Einstein condensates, in an approximation based on the suppression of the terms depending on the healing length scale -subleading with respect to scale set by the Hawking temperature -that leads to an equation that focuses on the vacuum polarization of the analogue scalar field.
From this general result, we have specialised to the canonical analogue black hole and have obtained an expression for the dynamics of the horizon and for how it is affected by the state of the quantum field. In particular, we have argued that in a regime in which the Hawking radiation is the leading contribution in the dynamics of the acoustic geometry (that would otherwise be stationary) it leads to a shrinking of the horizon radius.
This process of evaporation induces the depletion of atoms from the condensate to the excited part, and the rate of depletion is driven by the radial derivative of the vacuum polarization. The exchange of atoms and therefore the exchange of information between the two sectors of the system happens not only at the horizon, but also in the black hole region enclosed by it. In particular we have obtained an expression which ties together the rate of depletion to the number of atoms in the region within the the horizon, the temperature of the black hole and the value of the vacuum polarization near the singularity. These are all results in agreement with the conclusions reported in [30] with regards to the nature of a possible resolution of the information loss problem associated to black hole evaporation. What is peculiar of the presented calculation is the role played by the vacuum polarization in the inner region. In the gravitational case of the Schwarzschild black hole there is neither access to the information in the inner region nor knowledge of an underlying quantum gravity structure from which the classical geometry would emerge. Instead, in the canonical analogue black whole there are both these features. Not only in an experimental realisation one has access to the whole space involved in the dynamics, but it is now clear that the role played by the vacuum polarization of the analogue field, the Bogoliubov quasi-particles, is fundamental in understanding the exchange of information between the condensate and the excited part of the system.
We can then conclude that a deeper understanding of the dynamics of acoustic horizons in BEC based analogue gravity, and in particular of the associated phenomenology of depletion, could give -in spite of the different form of the gravitodynamic equations -precious insights towards a deeper understanding of the role of back-reaction in semiclassical gravity in the entanglement between quantum matter and quantum spacetime degrees of freedom. An entanglement which might be key for understanding the compatibility Hawking radiation with unitarity and hence for understanding black holes as thermodynamic objects.