# Hawking Radiation from Ultrashort Laser Pulse Filaments

In 1974 Hawking predicted that the space-time curvature at the event horizon of a black hole is sufficient to excite photons out of the vacuum and induce a continuous flux, referred to as Hawking radiation [1,2]. In a simplified description of the process, vacuum fluctuation pairs close to the horizon are split so that the inner photon falls in and the outside photon escapes away from the black hole. As the outgoing photon cannot return to the vacuum, it necessarily becomes a real entity, gaining energy at the expense of the black hole. It was soon realized that the essential ingredient of Hawking radiation was not the astrophysical black hole itself but rather the space-time curvature associated with the event horizon [3–7]. There are a wealth of physical systems that may exhibit event horizons ranging from flowing water or Bose-Einstein-Condensates to a moving refractive index perturbation (RIP) in a dielectric medium [3,4,8–10]. In a few words, referring to the case of optical pulses in a dielectric medium proposed by Philbin *et al.* [9], a laser pulse with large intensity, $I$, propagating in a nonlinear Kerr medium will excite a RIP given by $\delta n={n}_{2}I$ where ${n}_{2}$ is the so-called nonlinear Kerr index [11]. Light experiences an increase in the local refractive index as it approaches the RIP and is thus slowed down. By choosing appropriate conditions (frequency of the light and velocity of the RIP) it is possible to bring the light waves to a standstill in the reference frame comoving with the RIP, thus forming a so-called white hole event horizon, i.e., a point beyond which light is unable to penetrate. A similar mechanism may be observed with water waves or with any kind of waves in a flowing medium, so that the formation of an analogue event horizon is a rather universal phenomenon that may be studied in accessible laboratory conditions. What remains to be established, is whether Hawking radiation is actually emitted in the presence of an event horizon of any kind, be it analogue or astrophysical.

Recently, an alternative approach was proposed for generating controllable RIPs, namely, ultrashort laser pulse filamentation [10].

Ultrashort laser pulse filaments are intense laser pulses in a transparent Kerr medium (i.e., with a third order optical nonlinearity) characterized by a high-intensity spike that propagates apparently without diffraction over distances much longer than the Rayleigh length associated with the spike dimensions and have been proposed for many applications ranging from white light generation to the control of atmospheric conditions [12–16]. Filaments may either occur spontaneously when a powerful Gaussian shaped beam is loosely focused into the Kerr medium [12,13] or, alternatively, they may be induced by preshaping the laser pulse into a Bessel beam [17,18].

In this Letter we use ultrashort laser pulse filaments to create a traveling RIP in a transparent dielectric medium (fused silica glass) and we report experimental evidence of photon emission that on the one hand bears the characteristics of Hawking radiation and on the other is distinguishable and thus separate from other known photon emission mechanisms. We therefore interpret the observed photon emission as an indication of Hawking radiation induced by the analogue event horizon.

Neglecting any dependence on the transverse coordinates and dispersion it is possible to give a complete description of the event horizon associated with the RIP and calculate a blackbody temperature of the emitted photons in the laboratory reference frame [9,19]. However, we point out that the dielectric medium in which the RIP is created will always be dominated by optical dispersion and this in turn implies that in real settings the spectrum will not be that of a perfect blackbody and that in any case only a limited spectral portion of the full spectrum will be observable. To show this last point we describe the RIP as a perturbation induced by the laser pulse on top of a uniform, dispersive background refractive index ${n}_{0}$, i.e. $n(z,t,\omega )={n}_{0}(\omega )+\delta nf(z-vt)$, where $\omega $ is the optical frequency, $f(z-vt)$ is a function bounded by 0 and 1, that describes the shape of the laser pulse. In the reference frame comoving at velocity $v$ with the RIP, the event horizon in a 2D geometry is defined by $c/v=n$ which admits solutions only for RIP velocities satisfying the inequality [19]

Therefore, Eq. (1) predicts an emission spectrum with well-defined boundaries and it is precisely this feature of the spectral emission that we consider to be peculiar to analogue Hawking radiation in the present settings.

We note that two horizons may be associated with each RIP: the leading edge is the analogue of a black hole horizon while the trailing edge has time reversed features and is the analogue of a white hole horizon [9,10]. The experimental layout is shown in Fig. 1. The laser pulses are provided by a regeneratively amplified, 10 Hz repetition rate Nd:glass laser. The pulse duration is 1 ps and maximum energy is 6 mJ. The Bessel pulse filament with a cone angle $\theta =7\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{deg}$ is generated by a 20° fused silica axicon (conical lens), $F$, placed directly in front of the 2 cm fused silica Kerr sample where the RIP is generated. The input energy is varied in the $100\u20131200\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mu \mathrm{J}$ range.

Radiation from the filament is then collected at 90° with respect to the laser pulse propagation axis using a lens that images the filament on to the input slit of an imaging spectrometer. The spectrum is then recorded with a 16 bit, cooled CCD camera. This arrangement, in particular, detection at 90°, was chosen in order to strongly suppress or eliminate any spurious effects. More specifically: (i) Cerenkov-like radiation, i.e., in this context, radiation from a superluminal perturbation: spontaneous emission from a strictly superluminal perturbation has been considered in detail in Ref. [20]. The most relevant difference with respect to the present measurements is that Cerenkov-like emission occurs with no upper bound in the spectral emission window, in contrast to the limited Hawing spectrum given by Eq. (1). (ii) Four wave mixing (FWM) and self phase modulation (SPM) will not occur at 90° due to phase-matching constraints which imply that any newly generated frequencies will be generated at small angles with respect to the propagation axis [21]. In any case, we directly verified that for the relatively large $\theta =7\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{deg}$ Bessel cone angle used in our experiments, no FWM, SPM, or spectral broadening was observed at any angle, even in the forward direction. (iii) Rayleigh scattering will occur at 90° only for vertically polarized light and the scattering process will maintain the polarization state. In our experiments we used horizontally polarized light and in any case, in virtue of point (i), there is no input or generated light at the frequencies relevant for the present experiment. (iv) Fluorescence is certainly the main problem in these experiments. It bears many features in common with Hawking radiation yet, it may still be clearly distinguished from the latter. Figure 2(a) shows the overall spectrum measured at 90° and integrated over 30 laser pulse shots. In the figure, $R$, $F1$, and $F2$ indicate the laser pulse induced spontaneous Raman, nonbridging oxygen hole center (NBOHC) and oxygen deficient center (ODC) fluorescences, respectively. These fluorescence peaks are well documented features of fused silica [22,23]. This therefore allows us to fit the peaks with Gaussian functions (in the frequency domain) and subsequently subtract out the fluorescence signals, thus leading to greatly improved contrast and cleaner spectra.

Nevertheless, looking at the fluorescence spectrum it is clear that, given the possibility to tune the Hawking emission window by tuning the group velocity of the laser pulse, the most advantageous situation is that in which the window is located between 800 and 900 nm where there is no known fluorescence emission and moreover the CCD response is maximum. In fused silica, ${n}_{2}\sim 3\times {10}^{-16}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{cm}}^{2}/\mathrm{W}$ and typical intensities of $I\sim {10}^{12}\u2013{10}^{13}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{W}/{\mathrm{cm}}^{2}$ may be obtained, thus leading to a $\delta n\sim {10}^{-2}\u2013{10}^{-4}$. Therefore, if we take, for example, a Gaussian pulse centered at 1055 nm, with $\delta n=0.001$ and group velocity, ${v}_{G}=d\omega /dk$, determined solely by material dispersion, then according to Eq. (1) we would expect emission in the region between 500 and 510 nm. This would clearly fall in between the $F1$ and $F2$ fluorescence peaks and leads to noisy results even after subtraction of the fluorescence signal. By using Bessel pulses we were able to solve this problem. Figure 2(b) shows the Hawking spectral window predicted from Eq. (1), where the dashed line shows the value $c/{v}_{B}$ with $\theta =7\xb0$ and the Bessel pulse velocity is ${v}_{B}={v}_{G}/\mathrm{cos}\theta $. This line intersects the refractive index curves $n$ and $n+\delta n$ with $\delta n=0.001$ delimiting a window (shaded area) which now lies in the desired region, i.e., between 800 and 875 nm. Figure 3(a) shows the resulting spectra, integrated over 3600 laser shots. The black line shows a reference spectrum obtained from a Gaussian pulse, that clearly shows the absence of any signal notwithstanding the same peak intensity of the Bessel pulses and thus underlining the absence of any fluorescence or other possible emission signals when the RIP $v$ does not satisfy Eq. (1). The other four curves show the emitted spectra for four different input energies of the Bessel filament, as indicated in the figure. A clear photon emission is registered in the wavelength window predicted by Eq. (1). Moreover we verified that the emitted radiation was unpolarized (data not shown) thus further supporting the interpretation of a spontaneous emission. In addition, the measurements clearly show that the bandwidth of the emission is increasing with the input energy. The Bessel pulse intensity evolution along the propagation direction $z$ may be estimated analytically from the input energy using the formula $I(z)=2\pi {I}_{0}kz{tan}^{2}(\theta )\mathrm{exp}[-({z}^{2}{tan}^{2}\theta )/{w}_{0}^{2}]$ [24] where $k=\omega /c$ and ${I}_{0}$ and $\sqrt{2}{w}_{0}$ are the input Gaussian peak intensity and radius at $1/{e}^{2}$, respectively. By fitting the measured spectra with Gaussian functions [dashed curves in Fig. 3(a)] we may therefore estimate the bandwidth as a function of input energy and Bessel intensity. Moreover, using the fused silica dispersion relation shown in Fig. 2(b) it is possible to map the bandwidth into values of $\delta n$ so that we finally obtain Fig. 3(b) that shows the bandwidth and the $\delta n$ as a function of input energy and Bessel pulse peak intensity (at $z=1\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{cm}$ where measurements were performed). There is a clear linear dependence which is in qualitative agreement with the fact that, according to our interpretation, the emission bandwidth is predicted to depend on $\delta n$ which in turn is a linear function $\delta n={n}_{2}I$ of the pulse intensity. The slope of the linear fit gives ${n}_{2}=2.8\pm 0.5\times {10}^{-16}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{cm}}^{2}/\mathrm{W}$ which is in good agreement with the tabulated value of ${n}_{2}\sim 3\times {10}^{-16}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{cm}}^{2}/\mathrm{W}$ [12,25]. Therefore this shows that there is also an agreement at the quantitative level between the measurements and the model based on Hawking-like radiation emission.

In Fig. 4 we present additional data regarding Hawking emission from a spontaneous filament obtained by replacing the axicon with a 20 cm focal length lens and thus loosely focusing a $50\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mu \mathrm{J}$ Gaussian pulse into the fused silica sample. The spontaneous nonlinear Kerr dynamics lead to the formation of a filament [12] that is very similar to the Bessel filament in the sense that it is characterized by a very high-intensity peak that propagates over long distances, thus creating a strong RIP in the medium that may be expected to excite Hawking radiation in a similar fashion to that illustrated above. We underline some important differences with respect to the Bessel filament: the spontaneous filament is characterized by an intense pulse that travels slower than the input Gaussian pulse, i.e., $v<{v}_{G}$ (see [26] for a detailed description). Moreover the pulse accelerates during propagation so that a certain range of velocity values is covered during propagation, with low velocity at the beginning of the filament and higher velocity at the end. This has the effect of broadening the emission window, which is predicted from Eq. (1) to be between 270 and 450 nm. The light gray curve in Fig. 4(a) shows the spectrum measured when collecting light from the whole filament. As can be seen, the measured spectral bandwidth agrees well with the prediction. We also note that no emission was observed in this region with the Bessel filament, indicating that this is not a fluorescence signal. Moreover, by shifting and closing the input slit of the imaging spectrometer, we could collect spectra of the ending and beginning sections of the filament, that are shown as dark-filled curves in Figs. 4(a) and 4(b), respectively. Remarkably, different sections of the filament emit only at selected portions of the overall spectrum. We underline that, accounting for the different pulse velocities in the different filament sections, this striking behavior is in quantitative agreement with the predictions of Eq. (1).

Finally, we note that Eq. (1) and subsequent analysis, refers to the existence of phase velocity horizons. However, other analogue systems rely on the existence of group velocity horizons, i.e., points at which $v={v}_{G}$ (see, e.g., [27,28]) so that there is still an open question regarding the relative role of the two different types of horizon. In the specific case of the Bessel pulse used in our experiments, a group horizon does not even exist; i.e., there are no frequencies that satisfy ${v}_{G}(\omega )={v}_{B}$. Our measurements therefore indicate that, at least in the RIP setting, the phase velocity horizon alone may lead to photon emission in the predicted spectral window.

In summary, the same physics that underlie black hole evaporation in the form of Hawking radiation may be found and studied in other, more accessible systems. Our measurements highlight spontaneous emission of Hawking radiation from an analogue event horizon generated by an “evaporating” refractive index perturbation and suggest a path towards the experimental study of phenomena traditionally relegated to the areas of quantum gravity and astrophysics.

## References

- S. Hawking, Nature (London)
**248**, 30 (1974). - S. Hawking, Commun. Math. Phys.
**43**, 199 (1975). - W. G. Unruh, Phys. Rev. Lett.
**46**, 1351 (1981). - C. Barceló, S. Liberati, and M. Visser, Living Rev. Relativity
**8**, 12 (2005). - C. Barceló, S. Liberati, S. Sonego, and M. Visser, Phys. Rev. Lett.
**97**, 171301 (2006). - W. G. Unruh and R. Schützhold, Phys. Rev. D
**71**, 024028 (2005). - M. Visser, Phys. Rev. Lett.
**80**, 3436 (1998). - M. Novello, M. Visser, and G. E. Volovik,
*Artificial Black Holes*(World Scientific, Singapore, 2002). - T. Philbin
*et al.*, Science**319**, 1367 (2008). - D. Faccio
*et al.*, Europhys. Lett.**89**, 34004 (2010). - R. Boyd,
*Nonlinear Optics*(Academic Press, New York, 2008), 3rd ed. - A. Couairon
*et al.*, Phys. Rep.**441**, 47 (2007). - L. Bergé
*et al.*, Rep. Prog. Phys.**70**, 1633 (2007). - A. Gaeta, Science
**301**, 54 (2003). - J. Kasparian
*et al.*, Science**301**, 61 (2003). - P. Rohwetter
*et al.*, Nat. Photon.**4**, 451 (2010). - P. Polesana
*et al.*, Phys. Rev. Lett.**99**, 223902 (2007). - P. Polesana
*et al.*, Phys. Rev. E**73**, 056612 (2006). - F. Belgiorno
*et al.*, arXiv:1003.4150. - F. Belgiorno
*et al.*, Phys. Rev. Lett.**104**, 140403 (2010). - R. Gadonas
*et al.*, Opt. Commun.**196**, 309 (2001). - A. Zoubir
*et al.*, Phys. Rev. B**73**, 224117 (2006). - L. Skuja, J. Non-Cryst. Solids
**239**, 16 (1998). - A. T. Friberg
*et al.*, J. Opt. Soc. Am. A**13**, 743 (1996). - R. DeSalvo
*et al.*, IEEE J. Quantum Electron.**32**, 1324 (1996). - D. Faccio
*et al.*, Phys. Rev. A**78**, 033825 (2008). - G. Rousseaux
*et al.*, New J. Phys.**10**, 053015 (2008). - S. Weinfurtner
*et al.*, arXiv:1008.1911.