Advances in High Energy Physics

Volume 2014 (2014), Article ID 713574, 8 pages

http://dx.doi.org/10.1155/2014/713574

## Amplifying the Hawking Signal in BECs

^{1}Dipartimento di Fisica dell’Università di Bologna and INFN Sezione di Bologna, Via Irnerio 46, 40126 Bologna, Italy^{2}Museo Storico della Fisica e Centro Studi e Ricerche “Enrico Fermi”, Piazza del Viminale 1, 00184 Roma, Italy^{3}Dipartimento di Fisica dell’Università di Bologna, Via Irnerio 46, 40126 Bologna, Italy^{4}Departamento de Física Teòrica and IFIC, Universidad de Valencia-CSIC, C. Dr. Moliner 50, 46100 Burjassot, Spain

Received 7 January 2014; Accepted 6 March 2014; Published 27 April 2014

Academic Editor: Piero Nicolini

Copyright © 2014 Roberto Balbinot and Alessandro Fabbri. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The publication of this article was funded by SCOAP^{3}.

#### Abstract

We consider simple models of Bose-Einstein condensates to study analog pair-creation effects, namely, the Hawking effect from acoustic black holes and the dynamical Casimir effect in rapidly time-dependent backgrounds. We also focus on a proposal by Cornell to amplify the Hawking signal in density-density correlators by reducing the atoms’ interactions shortly before measurements are made.

#### 1. Introduction

Analogue models in condensed matter systems are nowadays an active field of investigation, not only on the theoretical side but, more importantly, also on the experimental one. The underlying idea is to reproduce in a condensed matter context peculiar and interesting quantum effects predicted by Quantum Field Theory in curved space, whose experimental verification in the gravitational context appears at the moment by far out of reach.

Many efforts are devoted to find the most famous of these effects, namely, the thermal emission by black holes predicted by Hawking in 1974 [1]. Among the condensed matter systems under examination, Bose-Einstein condensates appear as the most promising setting to achieve this goal [2–13]. The major problem one has to face experimentally is the correct identification of the signal corresponding to the analogue of Hawking radiation, namely, a thermal emission of phonons as a consequence of a sonic horizon formation, since it can be covered by other competing effects, like large thermal fluctuations.

A major breakthrough to overcome this problem came in 2008, when it was predicted that, as a consequence of Hawking radiation being a genuine pair creation process, a characteristic peak in the density correlation function of the condensate should appear for points situated on opposite sides with respect to the horizon [14]. This is the “smoking gun” of the Hawking effect. Soon after this proposal, Eric Cornell at the first meeting on “experimental Hawking radiation” held in Valencia in 2009 suggested that one can amplify this characteristic signal by reducing the interaction coupling among the atoms of the BEC shortly before measuring the density correlations [15].

Here we review in a simple pedagogical way, using toy models, how the analogous of Hawking radiation occurs in a supersonic flowing BEC and how the corresponding characteristic peak in the correlation function can be amplified according to Cornell’s suggestion. It should be stressed that nowadays correlation functions measurements are becoming the basic experimental tool to investigate Hawking-like radiation in condensed matter systems.

#### 2. BECs: The Gravitational Analogy and Hawking Radiation

A Bose gas in the dilute gas approximation is described by a field operator with equal-time commutator (see, e.g., [16]) satisfying the time-dependent Schrödinger equation where is the mass of the atoms, is the external potential, and is the nonlinear atom-atom interaction coupling constant. At sufficiently low temperatures a large fraction of the atoms condense into a common ground state which is described, in the mean field approach, by a -number field .

To consider linear fluctuations around this classical macroscopic condensate, one writes the bosonic field operator as where is a small (quantum) perturbation. and satisfy, respectively, Gross-Pitaevski where is the number density, and Bogoliubov-de Gennes equations with being the speed of sound.

Contact with the gravitational analogy (see, e.g., [13]) is achieved in the (long wavelength) hydrodynamic approximation, more easily realised by considering the density-phase representation for the Bose operator and the splitting , in which represent the linear (quantum) density and phase fluctuations, respectively. In terms of and we have Provided that the condensate density and velocity vary on length scales much bigger than the healing length (the fundamental length scale of the condensate), the BdG equation reduces to the continuity and Euler equations for and and these can be combined to give a second-order differential equation for which is mathematically equivalent to a Klein-Gordon (KG) equation where is the covariant KG operator from the acoustic metric For a flow which presents a transition from a subsonic () flow to a supersonic one ( in some region) the metric (8) describes an acoustic black hole, with horizon located at the surface where . The same analysis performed by Hawking in the gravitational case can be repeated step by step, leading to the prediction [17] that acoustic black holes will emit a thermal flux of phonons at the temperature where , with being the normal to the horizon, is the horizon’s surface gravity.

#### 3. The Model

To simplify the mathematics involved in the process, we will consider a 1D configuration (in 1D one should more correctly speak of quasi-condensation [18, 19]) in which and are constant and where the only nontrivial quantity is the speed of sound . As explained in [20], this can be achieved by varying the coupling constant (and therefore ) and the external potential but keeping the sum constant. In this way, the plane-wave function , where is the condensate velocity and , where is the chemical potential of the gas, is a solution of (4) everywhere. Note that such a stationary configuration is difficult to reach experimentally; nevertheless it gives results similar to those obtained by more realistic configurations [21].

The fluctuation operator is expanded in the usual form in terms of positive and negative norm modes as where and are quasi-particle’s annihilation and creation operators. From (5) and its Hermitean conjugate, we see that the modes and satisfy the coupled differential equations The normalizations are fixed, via integration of the equal-time commutator obtained from (1), namely, by

In order to get simple analytical expressions, in the following we will consider simple models with step-like discontinuities in the speed of sound and impose the appropriate boundary conditions for the modes that are solutions to (11). For more general profiles a numerical analysis can be performed; see, for example, [22–25].

##### 3.1. Acoustic Black Holes and the Hawking Effect

A simple analytical model of an acoustic black hole [26] can be obtained by gluing two semi-infinite stationary and homogeneous 1D condensates, one subsonic () and the other supersonic (), along a spatial discontinuity at (see [27], to which we refer for more detailed explanations throughout this paragraph, and the references therein): . We take ; that is, the flow is from right to left, and . We denote the modes solutions in each homogeneous region and corresponding to the fields and as so that (11) become while the normalization condition (13) gives The combination of the two equations (15) gives the Bogoliubov dispersion relation for a one-dimensional Bose liquid flowing at constant velocity containing the positive and negative norm branches which, for the subsonic and supersonic regions, are given, respectively, in Figures 1 and 2. Moreover, inserting the relation between and from (15) into (16) we find the mode normalizations where are the roots of the quartic equation (17) at fixed .

In the subsonic case (17) admits two real and two complex solutions. Regarding the real solutions, Figure 1, we call and the ones corresponding to, respectively, negative and positive group velocity (the other two complex conjugated solutions correspond to, respectively, spatially decaying and growing modes). In the supersonic case (see Figure 2) we see that, for the most interesting regime (), there are now four real solutions, corresponding to four propagating modes: (present also in the hydrodynamical limit ) and , two of which ( and ) belong to the negative norm branch.

To find modes evolution for all one needs to write down the general solutions for () in the left supersonic () and the right subsonic () regions (we restrict to the most interesting case ) (the expansions for are the same up to the replacement ) and impose, from (11), the matching conditions where indicates the variation across the jump at , allowing us to write down the relations between left and right amplitudes through a scattering matrix in the form This allows us to construct explicitly the decomposition of the field operator in terms of the “in” and “out” bases. The “in” basis is constructed with modes propagating from the asymptotic regions () towards the discontinuity (), while the “out” basis is constructed with modes propagating away from the discontinuity to . Looking at Figures 1 and 2, we see that unit amplitude modes defined on the left moving and right moving momenta define the ingoing scattering states, while unit amplitude modes defined on the right moving and left moving momenta define the outgoing scattering states. One can then write down the “in” decomposition in terms of the “in” scattering states or, equivalently, on the basis of the “out” scattering ones. Note that since belongs to the negative norm branch the corresponding “in” mode is multiplied by a creation operator (the same thing happens, in the “out” decomposition, for ). Using (21) one can construct the S-matrix relating and modes which is not trivial since it mixes positive and negative norm modes. As a consequence, the Bogoliubov transformation between “in” and “out” creation and annihilation operators is also not trivial because it mixes creation and annihilation operators. This has the crucial consequence that the “in” and “out” Hilbert spaces are not unitary related; in particular the corresponding vacua are different; that is, . The physical consequence is that if we prepare the system in the vacuum state, so that there are no incoming phonons at , we will have, at late times, outgoing quanta on both sides of the horizon: the vacuum has spontaneously emitted phonons, mainly in the channel (Hawking quanta) and (partners). The analytical calculations show that the number of emitted Hawking quanta [26] and partners follow an approximate (low-frequency) thermal spectrum [21], the proportionality factor allowing the identification of a Hawking temperature () in this idealised setting. We can understand the mechanism by which Hawking radiation is emitted by looking, using (6), at the equal-time density-density correlator whose main contribution in the and sector comes from the term The existence of the peak at was first pointed out in [14] in the hydrodynamical approximation using QFT in curved space techniques. The physical picture that emerges is that Hawking quanta and partners are continuously created in pairs from the horizon at each time , propagate on opposite directions at speeds and , and after time are located at and related as in (28). The existence of the Hawking peak was nicely confirmed by numerical “ab initio” simulations with more realistic configurations performed in [20]. How correlation measurements can reveal the quantum nature of Hawking radiation is discussed in [28–30].

##### 3.2. Analog Dynamical Casimir Effect

A distinct type of pair-creation takes place in time-dependent backgrounds, with one important example being quantum particle creation in cosmology. We will study the analogue of these phenomena in BECs with another simple model in which the speed of sound has a step-like discontinuity in time at separating two “initial” and “final” infinite homogeneous condensates: ; see for more details [31, 32]. The modes solutions in the initial () and final () regions are now of the type for which (11) become while the normalization condition (13) yields giving Here, corresponds to the two real solutions to (17), which is quadratic in at fixed . These read where corresponds to the positive norm branch and corresponds to the negative norm one. To find modes evolution for all one first writes down the general solutions in the initial and final regions for () (for we have the same expansion with replaced by ) and impose the matching conditions, from (11), where now indicates the variation across the discontinuity at . They allow the final amplitudes to be related to the initial ones through the matrix : where One can easily construct “in” and “fin” decompositions for the field by considering initial (final) unit amplitude positive norm modes () () modes for all using (36): The modes are related through a nontrivial Bogoliubov transformation mixing positive and negative norm modes where and, consequently, also the relation between “in” and “fin” annihilation and creation operators will mix annihilation and creation operators, implying again that the two decompositions are inequivalent and in particular the two vacuum states and are different. The physical consequence is that the step-like discontinuity at will induce particle creation, the features of which can be understood by looking at the time-dependent terms of the one-time density-density correlator which in the hydrodynamical limit reads At and everywhere in space correlated pairs of particles with opposite momentum are created out of the vacuum state, with velocities (left moving) and (right moving). At time such particles are separated by a distance which is indeed the correlation displayed in (41). This effect was recently observed in [33] by considering homogeneous condensates with trapping potential rapidly varying in time, where correlation functions in velocity/momentum space were measured.

##### 3.3. Amplification of the Hawking Signal in Density Correlators

It has been argued by Cornell [15] that a way to amplify the Hawking signal in density-density correlators is to reduce the interactions shortly before measuring the density correlations. Since , reducing means that the speed of sound is also reduced. We will model this situation by matching our idealised acoustic black hole configuration of Section 3.1 with a final infinite homogeneous condensate characterised by a small sound velocity (); see Figure 3. To study this situation, in which a spatial step-like discontinuity in at is combined with a temporal step-like discontinuity at some , we will use the tools introduced in the previous two subsections. We will calculate the density-density correlator in the “in” vacuum by expanding the density operator in the “in” decomposition to get By expressing the “in” modes in terms of the “out” modes and in the absence of the temporal step-like discontinuity (say ) the Hawking signal is given by In the presence of the temporal step-like discontinuity we need to evolve the relevant modes () and (), at the same value of , across the discontinuity at . Going to the basis and considering small () we have where , , and with , .

It is useful to rewrite (46) and (47) in terms of to compare with the standard result without the temporal step-like discontinuity The standard result is, from (45), a stationary peak at weighted by the -independent factor (see (24) and (25)) The effect of the temporal step-like discontinuity at (see (48)) is to modify this signal into a main signal located at with strength and three smaller signals located at = − , − = − , and − = − with strengths given by (52) in which is substituted, respectively, by , , and .

We see immediately that we loose the stationarity of the Hawking signal (49) and that the main signal is multiplied by with respect to the standard result. This results indeed in an amplification (i.e., the above term is 1) when . Being the Hawking peak expected to be of order for realistic experimental settings [20] (where was considered), we obtain an amplification factor for .

#### 4. Conclusions

In this paper we have briefly reviewed the analysis of the analog Hawking effect and of the analog dynamical Casimir effect by considering simple analytical models of Bose-Einstein condensates in which the speed of sound has step-like discontinuities. We focussed on the study of the density-density correlators which show, in the former case, the existence of a characteristic stationary Hawking quanta-partner peak located at (28) and, in the latter, that of a time-dependent feature (42). Following a suggestion by Cornell, we combined these two analyses to construct a model in which the atoms’ interactions are rapidly lowered (and so the speed of sound) before the correlations are measured. This results in an amplification of the main Hawking peak, now time-dependent and located at (51), by the factor (53) that could be useful in the experimental search.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgment

The authors thank I. Carusotto for useful discussions.

#### References

- S. W. Hawking, “Particle creation by black holes,”
*Communications in Mathematical Physics*, vol. 43, no. 3, pp. 199–220, 1975. View at Publisher · View at Google Scholar · View at MathSciNet - L. J. Garay, J. R. Anglin, J. I. Cirac, and P. Zoller, “Sonic black holes in dilute Bose-Einstein condensates,”
*Physical Review A*, vol. 63, no. 2, Article ID 023611, 2001. View at Publisher · View at Google Scholar · View at Scopus - L. J. Garay, J. R. Anglin, J. I. Cirac, and P. Zoller, “Sonic analog of gravitational black holes in Bose-Einstein condensates,”
*Physical Review Letters*, vol. 85, no. 22, pp. 4643–4647, 2000. View at Publisher · View at Google Scholar · View at Scopus - O. Lahav, A. Itah, A. Blumkin et al., “Realization of a sonic black hole analog in a Bose-Einstein condensate,”
*Physical Review Letters*, vol. 105, no. 24, Article ID 240401, 2010. View at Publisher · View at Google Scholar · View at Scopus - I. Shammass, S. Rinott, A. Berkovitz, R. Schley, and J. Steinhauer, “Phonon dispersion relation of an atomic Bose-Einstein condensate,”
*Physical Review Letters*, vol. 109, Article ID 195301, 2012. View at Publisher · View at Google Scholar - R. Schley, A. Berkovitz, S. Rinott, I. Shammass, A. Blumkin, and J. Steinhauer, “Planck distribution of phonons in a Bose-Einstein condensate,”
*Physical Review Letters*, vol. 111, Article ID 055301, 2013. View at Publisher · View at Google Scholar - C. Barceló, S. Liberati, and M. Visser, “Towards the observation of Hawking radiation in Bose-Einstein condensates,”
*International Journal of Modern Physics A*, vol. 18, no. 21, pp. 3735–3745, 2003. View at Publisher · View at Google Scholar · View at Scopus - C. Barceló, S. Liberati, and M. Visser, “Probing semiclassical analog gravity in Bose-Einstein condensates with widely tunable interactions,”
*Physical Review A*, vol. 68, Article ID 053613, 2003. View at Publisher · View at Google Scholar - S. Giovanazzi, C. Farrell, T. Kiss, and U. Leonhardt, “Conditions for one-dimensional supersonic flow of quantum gases,”
*Physical Review A*, vol. 70, no. 6, Article ID 063602, 2004. View at Publisher · View at Google Scholar · View at Scopus - R. Schützhold, “Detection scheme for acoustic quantum radiation in Bose-Einstein condensates,”
*Physical Review Letters*, vol. 97, Article ID 190405, 2006. View at Publisher · View at Google Scholar - S. Wüster and C. M. Savage, “Limits to the analog Hawking temperature in a Bose-Einstein condensate,”
*Physical Review A*, vol. 76, Article ID 013608, 2007. View at Publisher · View at Google Scholar - Y. Kurita and T. Morinari, “Formation of a sonic horizon in isotropically expanding Bose-Einstein condensates,”
*Physical Review A*, vol. 76, no. 5, Article ID 053603, 2007. View at Publisher · View at Google Scholar · View at Scopus - C. Barcelo, S. Liberati, and M. Visser, “Analogue gravity,”
*Living Reviews in Relativity*, vol. 8, article 12, 2005. View at Publisher · View at Google Scholar - R. Balbinot, A. Fabbri, S. Fagnocchi, A. Recati, and I. Carusotto, “Nonlocal density correlations as a signature of Hawking radiation from acoustic black holes,”
*Physical Review A*, vol. 78, no. 2, Article ID 021603, 2008. View at Publisher · View at Google Scholar · View at Scopus - E. Cornell, “Towards the observation of Hawking radiation in condensed matter systems,” Talk given at the workshop, Valencia, Spain, 2009, http://www.uv.es/workshopEHR.
- L. Pitaevskii and S. Stringari,
*Bose-Einstein Condensation*, vol. 116 of*International Series of Monographs on Physics*, The Clarendon Press, Oxford, UK, 2003. View at MathSciNet - W. G. Unruh, “Experimental black-hole evaporation?”
*Physical Review Letters*, vol. 46, no. 21, pp. 1351–1353, 1981. View at Publisher · View at Google Scholar · View at Scopus - Y. Castin, “Simple theoretical tools for low dimension Bose gases,”
*Journal de Physique IV France*, vol. 116, pp. 89–132, 2004. View at Publisher · View at Google Scholar - C. Mora and Y. Castin, “Extension of Bogoliubov theory to quasicondensates,”
*Physical Review A*, vol. 67, Article ID 053615, 2003. View at Publisher · View at Google Scholar - I. Carusotto, S. Fagnocchi, A. Recati, R. Balbino, and A. Fabbri, “Numerical observation of Hawking radiation from acoustic black holes in atomic Bose-Einstein condensates,”
*New Journal of Physics*, vol. 10, Article ID 103001, 2008. View at Publisher · View at Google Scholar · View at Scopus - P.-É. Larré, A. Recati, I. Carusotto, and N. Pavloff, “Quantum fluctuations around black hole horizons in Bose-Einstein condensates,”
*Physical Review A*, vol. 85, no. 1, Article ID 013621, 2012. View at Publisher · View at Google Scholar · View at Scopus - J. Macher and R. Parentani, “Black-hole radiation in Bose-Einstein condensates,”
*Physical Review A*, vol. 80, no. 4, Article ID 043601, 2009. View at Publisher · View at Google Scholar · View at Scopus - S. Finazzi and R. Parentani, “Spectral properties of acoustic black hole radiation: broadening the horizon,”
*Physical Review D*, vol. 83, no. 8, Article ID 084010, 2011. View at Publisher · View at Google Scholar · View at Scopus - S. Finazzi and R. Parentani, “Hawking radiation in dispersive theories, the two regimes,”
*Physical Review D*, vol. 85, Article ID 124027, 2012. View at Publisher · View at Google Scholar - S. Bar-Ad, R. Schilling, and V. Fleurov, “Nonlocality and fluctuations near the optical analog of a sonic horizon,”
*Physical Review A*, vol. 87, Article ID 013802, 2013. View at Publisher · View at Google Scholar - A. Recati, N. Pavloff, and I. Carusotto, “Bogoliubov theory of acoustic Hawking radiation in Bose-Einstein condensates,”
*Physical Review A*, vol. 80, no. 4, Article ID 043603, 2009. View at Publisher · View at Google Scholar · View at Scopus - R. Balbinot, I. Carusotto, A. Fabbri, C. Mayoral, and A. Recati, “Understanding hawking radiation from simple models of atomic Bose-Einstein condensates,” in
*Analogue Gravity Phenomenology*, vol. 870 of*Lecture Notes in Physics*, pp. 181–219, 2013. View at Google Scholar - K. V. Kheruntsyan, J.-C. Jaskula, P. Deuar et al., “Violation of the Cauchy-Schwarz inequality with matter waves,”
*Physical Review Letters*, vol. 108, Article ID 260401, 2012. View at Publisher · View at Google Scholar - J. R. M. de Nova, F. Sols, and I. Zapata, “Violation of Cauchy-Schwarz inequalities by spontaneous Hawking radiation in resonant boson structures,” http://arxiv.org/abs/1211.1761.
- S. Finazzi and I. Carusotto, “Entangled phonons in atomic Bose-Einstein condensates,” http://arxiv.org/abs/1309.3414.
- I. Carusotto, R. Balbinot, A. Fabbri, and A. Recati, “Density correlations and analog dynamical Casimir emission of Bogoliubov phonons in modulated atomic Bose-Einstein condensates,”
*European Physical Journal D*, vol. 56, no. 3, pp. 391–404, 2010. View at Publisher · View at Google Scholar · View at Scopus - C. Mayoral, A. Fabbri, and M. Rinaldi, “Steplike discontinuities in Bose-Einstein condensates and Hawking radiation: dispersion effects,”
*Physical Review D*, vol. 83, Article ID 124047, 2011. View at Publisher · View at Google Scholar - J. C. Jaskula, G. B. Partridge, M. Bonneau et al., “Acoustic analog to the dynamical casimir effect in a Bose-Einstein condensate,”
*Physical Review Letters*, vol. 109, Article ID 220401, 2012. View at Publisher · View at Google Scholar