Scalar radiation from a radially infalling source into a Schwarzschild black hole in the framework of quantum field theory

We investigate the radiation to infinity of a massless scalar field from a source falling radially towards a Schwarzschild black hole using the framework of the quantum field theory at tree level. When the source falls from infinity, the monopole radiation is dominant for low initial velocities. Higher multipoles become dominant at high initial velocities. It is found that, as in the electromagnetic and gravitational cases, at high initial velocities the energy spectrum for each multipole with l≥1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$l \ge 1$$\end{document} approximately is constant up to the fundamental quasinormal frequency and then drops to zero. We also investigate the case where the source falls from rest at a finite distance from the black hole. It is found that the monopole and dipole contributions in this case are dominant. This case needs to be carefully distinguished from the unphysical process where the source abruptly appears at rest and starts falling, which would result in radiation of an infinite amount of energy. We also investigate the radiation of a massless scalar field to the horizon of the black hole, finding some features similar to the gravitational case.


Introduction
Black holes (BHs) stand out as the most relevant and simple objects described by General Relativity. BHs are trapped regions, even to light, due to their extremely intense gravitational field. The boundary of no return from which light cannot escape, the BH event horizon, is determined as a function of only three parameters associated with the BH, i.e. mass, angular momentum and electric charge [1]. Particles falling into BHs emit radiation which carries information about the a e-mail: laoliveira@ufpa.br b e-mail: crispino@ufpa.br c e-mail: atsushi.higuchi@york.ac.uk event horizon to infinity, as a "fingerprint" of the BHs [2][3][4][5]. Thus, in principle, evidence for the existence of an event horizon and therefore for the existence of BHs can be obtained by analyzing the radiation emitted from a source falling into BHs. The scientific literature about the dynamics of a test particle falling into BHs has developed significantly in the early 1970s. The existing results for the problem of radiation emission from a particle falling radially into BHs were obtained using the formalism of the Classical Field Theory (CFT) [2][3][4][5][6][7][8][9][10][11][12][13][14][15][16][17]. The investigation of this kind of problems from the viewpoint of the Quantum Field Theory (QFT) has not been carried out. The formalism of QFT applied to the problem of radiation emission has been used for the determination of the radiation emission by sources and charges rotating around a BH, known as synchrotron radiation [18][19][20][21][22][23][24][25]. Furthermore, QFT has been used to investigate radiation emission from an uniformly accelerated source in flat spacetimes [26,27] and also to investigate the interaction of sources with Hawking radiation [28][29][30][31][32][33][34]. In this paper, using QFT at tree level, we investigate in detail the properties of the radiation emission due to the radial infall of a particle source of a massless scalar field into a Schwarzschild BH. One of the main advantages in computing the emitted energy using the framework of QFT at level tree is that this approach makes the extension to the radiative quantum corrections more straightforward. We note in passing that the change in the geometry due to the Hawking radiation for an astrophysical black hole is extremely small [35] and, as a result, that it is legitimate to use the eternal black hole in our calculations.
The remainder of this paper is organized as follows. In Sect. 2 we briefly review the general formalism used in this paper. In Sect. 3 we describe the radial infall of a source into a Schwarzschild BH according to General Relativity. In Sect. 4 we obtain expressions for the emitted energy spectrum and the total emitted energy for a quantum scalar field minimally coupled to the scalar source. In Sect. 5 we obtain the zerofrequency limit of the emitted energy spectra, using approximate analytic solutions. In Sect. 6 we analyze and discuss our numerical results for the radiation emission obtained from the viewpoint of QFT. We summarize some features of our results in Sect. 7. We use natural units with G = c =h = 1, unless otherwise stated.

General formalism of QFT in Schwarzschild spacetime
The total Lagrangian density with a classical source j (x μ ), minimally coupled to a massless and chargeless scalar field Φ (x μ ), can be written by where g ≡ det(g μν ) is the determinant of the metric g μν . The line element ds 2 = g μν dx μ dx ν of the Schwarzschild spacetime can be written where f (r ) = 1 − 2M/r . Note that f (r = r h ) = 0, with r h ≡ 2M being the position of the event horizon of the Schwarzschild BH, and f (r → ∞) = 1. This spacetime is asymptotically flat. The scalar fieldΦ (x μ ) can be expanded in terms of a complete set of positive-and negative-frequency modes, u ω l m and u * ω l m , aŝ with ω > 0, where " * " denotes complex conjugation and a ω l m and its Hermitian adjointâ † ω l m are, respectively, the annihilation and the creation operators [36]. These operators satisfy the following non-vanishing commutation relations: Since the Schwarzschild spacetime is spherically symmetric, the positive-frequency modes u ω l m (x μ ), can be written by which satisfy the Klein-Gordon equation where C ω l m is a normalization constant, Y l m (θ, φ) are the spherical harmonics [37] and √ −g = r 2 sin θ [obtained from the line element (2)]. The Klein-Gordon inner product for the mode functions is defined as follows: with dS μ = r 2 sin θ dr dθ dφ δ μ 0 / f (r ), where S is a constanttime hypersurface, which is a Cauchy surface [38]. The commutation relations (4) imply that the modes u ω l m are normalized as follows: The conditions (8) determine the normalization constant C ω l m in Eq. (5). Substituting Eq. (5) into the Klein-Gordon equation (6), we obtain the following ordinary differential equation for ψ ω l (x): where V l (r ) is the effective potential, given by The Regge-Wheeler coordinate r * is defined by dr * /dr ≡ f (r ) −1 , which for the Schwarzschild BH can be explicitly written as Equation (10) admits two independent sets of solutions which can be represented by the modes ψ up ω l (r * ), purely incoming from the past horizon H − , and the modes ψ in ω l (r * ), purely incoming from the past null infinity J − . The solutions ψ up ω l (r * ) and ψ in ω l (r * ) satisfy, respectively, the following boundary conditions at the event horizon (r * → −∞) and at spatial infinity (r * → ∞): and The functions χ hor and χ inf are defined to be of the form exp(−iωr * ) and exp(iωr * ), respectively, at leading order in 1/r * . In our numerical computations we write these functions near the horizon and near spatial infinity as where the coefficients a j and b j are obtained from Eq. (10), with the choice a 0 = 1 and b 0 = 1. We let j max = 10 in our computation. The coefficients A up/in ω l and B up/in ω l are determined by matching the boundary conditions (13) and (14) with the numerical solution ψ ω l (r * ), obtained from Eq. (10), at the event horizon r * → −∞ and at spatial infinity r * → ∞, respectively. Note that the coefficients A The reflection and transmission coefficients satisfy the following relation: From the boundary conditions (13) and (14) and Eq. (8) we find that the normalization constant C ω l m can be written We note that, because of the complex conjugation in Eq. (5) of ψ ω l (r ), the modes u ω l m (x μ ) are the modes purely ingoing into the future horizon H + [with ψ up * ω l (r )] or purely outgoing to the future null infinity J + [with ψ in * ω l (r )]. These are the modes that are used for computing the radiation into the horizon and to null infinity in the next section.

Radial infall of a source into a Schwarzschild black hole
We consider a source falling radially into a Schwarzschild BH. The source has a zero angular momentum as a result. Without loss of generality we let the source fall along the z-axis. The stress-energy tensor for a point source can be written by where τ is the source's proper time.
By setting j (x μ ) = T ν ν (corresponding to a scalar source) we find the following expression for a massive source: where (r, θ, φ) = (r s , θ s , φ s ) refers to the spatial coordinates of the source at given time t in spherical polar coordinates, q is a coupling constant between the source and a massless and chargeless scalar fieldΦ [18,32], and v t is the contravariant t-component of the four-velocity of the source. The factor 1/v t makes the source boost invariant along its trajectory. The four-velocity v μ of a source infalling radially is given by the following expression: where E is the source's conserved energy divided by its rest mass [1]. If the source has initial position r = r 0 and velocity v 0 in the ingoing radial direction (at t = 0), then [38] Using Eqs. (2) and (23) it is possible to find an expression for the modulus of the velocity of the source (falling radially) at position r in the static frame [39], namely In Fig. 1 we exhibit plots of U r for selected values of v 0 and r 0 . We write the radial coordinate of the source as r s ≡ r s (t s ) with t s being the time coordinate in Eq. (2) associated with the trajectory of the radial infall along its geodesic. The function t s (r s ) is obtained from the relation This formula follows from Eq. (23). Using the properties of the Dirac delta function and Eqs. (23) and (26) Eq. (22) can be rewritten as follows: where t s (r ) is the inverse function of r = r s (t s ). The constants θ s will be set to 0 (and then φ s will be ambiguous and can be set to any value).

Total emitted energy
Now the computation of the total emitted energy (for each multipole number l and azimuthal number m) by a source j (x μ ) minimally coupled to a scalar fieldΦ (x μ ) in a background of a spherically symmetric spacetime using QFT at tree level can be done. The starting point is the following expression [18,19]: where the labels "inf" and "hor" correspond to the energy radiated, respectively, to infinity and to the event horizon. A up/in ω l m are the emission amplitudes at tree level given by corresponding to a transition between the vacuum state and one scalar particle state. By recalling that with |0 being the Boulware vacuum [36,40], and |ω l m = a † ω l m |0 , we find, using the commutation relations (4), Then, substituting Eq. (3) into Eq. (29) and using Eqs. (30) and (31), we find If we had chosen the Unruh or Hartle-Hawking vacuum [41,42] there would be absorption and stimulated emission of the scalar particles. The rates of these two processes would be exactly the same and, as a result, the net emission would be the same as in the Boulware vacuum (see, e.g. [18,43]). It is well known that the Boulware vacuum is unphysical because the expectation value of the stress-energy tensor is singular at the past and future horizons for this state [44]. Our results can be regarded to be about the net emission from the scalar source in the Unruh vacuum, which is more physical.
As stated before, the motion is along the z-axis without loss of generality because of the spherical symmetry of Schwarzschild spacetime. Thus we consider θ s = 0 and φ s = 0. (The value of φ s is arbitrary once we have θ s = 0. However, we choose this value for definiteness.) As a consequence only the mode m = 0 contributes to the emission amplitude [9,17]. Because of this fact the only spherical harmonics which will be associated to non-vanishing amplitudes are Y l 0 (θ, φ) with From now on we will omit the azimuthal quantum number m from A up/in ω l m and E hor/inf l m for the reason stated above. By substituting Eqs. (5), (20), (27) and (33) into Eq. (32) the following expression for the emission amplitude is found: Using Eq. (26), Eq. (34) can be rewritten by integrating by parts as follows: where Equation (34), or equivalently (35), in fact represents the emission amplitude from a scalar source that suddenly appears at t s = 0 and starts falling. To obtain the amplitude from a scalar source that is static until t s = 0 and then starts falling, we first convert the r s integral in Eq. (34) for the amplitude back to the t s integral using Eq. (26) and find Since r s (t s ) = r 0 for t s ≤ 0 we can readily evaluate the integral over (−∞, 0] by changing exp(iωt s ) to exp(iωt s + t s ) with > 0 and letting → 0. We integrate by parts over [0, ∞) and find that the boundary term is canceled by the integral over (−∞, 0]. Thus the result turns out to be Eq. (35) with B = 0.

The zero-frequency limit
In this section we obtain the zero-frequency limit (ZFL) of the spectra of the energy emitted to infinity (Eq. (28)), using approximate analytic solutions of Eq. (10) in order to check the results obtained for the energy spectra considering the full numerical solution.
Since there is only one parameter M with dimensions in Schwarzschild spacetime the low-frequency limit is the limit where ω M −1 . Since the energy spectra E hor/inf ω l are dimensionless, they are functions of Mω for r 0 = ∞. Hence, their ω → 0 limit for r 0 = ∞ is achieved by letting M → 0 [45]. (Note that r 0 introduces another parameter with dimensions if it is finite.) In this limit we may replace the potential V l (r ) in Eq. (11) by its leading term for a large r , and we let Equation (10) in this approximation is the wave equation in Minkowski spacetime with the following familiar solutions: where j l (ωr ) are the spherical Bessel functions [37], up to a phase factor. Thus, in the limit M → 0 Eq. (34) for r 0 = ∞ becomes 1 Some authors define E hor/inf ω l to be 1/2 times ours because the total energy emitted is obtained by integrating E hor/inf over ω from −∞ to ∞ in their case.
where Q l (z) is the Legendre function of the second kind with the branch cut [−1, 1]. Hence for r 0 = ∞ we find the ω → 0 limit of the spectra of the energy emitted to infinity, E inf Now, if r is held fixed, then (see, e.g. Sec. VI of Ref. [32]) For r 0 finite this limit can be used for all r ≤ r 0 in Eq. (35) (with B = 0). Thus, we find Note that this formula is not valid for r 0 = ∞ because the limits ω → 0 and r 0 → ∞ do not commute. The zerofrequency limit of (A up ω l ) −1 ψ up ω l (r ), relevant to the radiation emitted to the horizon, is also known (see, e.g. Ref. [32]). It can be used to find the zero-frequency limit of the spectra of energy emitted to the horizon, E hor ω l , for finite r 0 in a similar manner with the following result:

Numerical results
In Sect. 4 we described how to find the energy spectra of massless scalar radiation from a source freely falling radially. In this section numerical evaluation of the energy spectra is presented. As for the source's motion, the following two distinct cases are being considered: (i) the source coming from r = ∞ with a non-vanishing initial velocity v 0 and (ii) the source released from rest at a certain position r = r 0 . The results for the spectra for the energy emitted to the horizon, E hor ω l , and to infinity, E inf ω l , for the radial infall of a source starting from r = ∞ for selected values of initial velocity v 0 are shown in Fig. 2. These results were obtained numerically from Eqs. (34) and (38)  (plots on the left) and E inf ω l (plots on the right) as functions of ω for selected values of the multipole number, obtained from Eqs. (34) and (38) for the radial infall of a source from spatial infinity (r 0 → ∞), with the initial velocity v 0 = 0.5 (plots at the top), v 0 = 0.75 (plots in the middle) and v 0 = 0.95 (plots at the bottom) same results by using Eq. (35), instead of Eq. (34), and setting B = 0, because the boundary term (36) vanishes.] We note that the spectra of energy emitted to the horizon, E hor ω l , (plots on the left in Fig. 2) starts from zero (at ω = 0), reaches a maximum and then slowly decreases to zero (for high values of ω). Note that the energy emitted to the horizon does not decrease as the multipole number l increases unlike the energy emitted to infinity. This behavior of the spectra E hor ω l is similar to that of the corresponding spectrum for the gravitational radiation [5]. It reflects the fact that the source has infinite self energy due to the Coulomb-like potential and that the region of large energy density passes through the horizon as the source approaches it. The spectrum of energy emitted to infinity, E inf ω l , (plots on the right in Fig. 2) starts from a non-vanishing finite value (at ω = 0) and goes to zero for high frequencies. The contribution of higher multipoles decreases rapidly with increasing l as in the gravitational case [4]. These spectra were studied by Brito [45] and our results are in agreement with his. It is interesting that for high ini-tial velocities the spectrum for each l ≥ 1 is approximately constant and drops to zero around the fundamental quasinormal frequency. [These frequencies ω qn are given as ω qn r h = 0.215 (l = 0), 0.586 (l = 1), 0.967 (l = 2), 1.351 (l = 3) and 1.733 (l = 4) (see, e.g. Ref. [46]).] The spectra for the gravitational and electromagnetic radiation behave in a similar manner [7,13]. (plots on the left) and E inf ω l (plots on the right) as a function of ω for selected values of the multipole number, obtained from Eqs. (34) and (38), considering a source infalling from rest (v 0 = 0) at positions r 0 = 5 r h (plots at the top), r 0 = 2.5 r h (plots in the middle) and r 0 = 1.1 r h (plots at the bottom) The results for the emitted energy spectra E hor ω l and E inf ω l for the radial infall of a source starting from rest for selected values of position r 0 obtained numerically from Eqs. (35) (with B = 0) and (38) are shown in Fig. 3 as functions of ω and for selected values of the multipole number l. We note that the emission to infinity is dominated by lower multipoles (l = 0, 1) while there is a substantial contribution from the modes with l ≥ 2 to the emission to the horizon, reflecting the region of high energy density surrounding the source passing through the horizon.
Next we show the same spectra using Eq. (34) instead of Eq. (35) (with B = 0). As we stated before, the emission in this case is from a source that emerges suddenly at r = r 0 and starts falling. The results are shown in Fig. 4. We note that Numerical estimates of the emitted energy radiated to infinity, E inf l , as a function of v 0 for multipole numbers l = 0, 1, 2, 3, 4, obtained from Eq. (28) the spectra of the emitted energy both to the horizon and to infinity, E hor ω l and E inf ω l , vanish for both ω = 0 and ω → ∞, (i) oscillating between these limits for mid-to-large values of r 0 (e.g., r 0 = 5 r h ) or (ii) behaving with a Gaussian-like profile, for small values of r 0 (e.g., r 0 = 2.5 r h and r 0 = 1.1 r h ).
In the plots of Fig. 5 we show the energy emitted to infinity E inf l , obtained from Eq. (28), as a function of the initial velocity v 0 . The monopole (l = 0) emission is dominant for low initial velocities while higher multipoles (l ≥ 1) become significant for high initial velocities. We note that for l = 0 as the initial velocity v 0 increases, the value of E inf 0 decreases. It is interesting that this behavior for the multipole number l = 0 is the opposite to that for the multipole numbers l ≥ 1. In the plots of Fig. 6 we show the emitted energy E inf l , obtained from Eq. (28), as a function of the position r 0 . We note that, as the position r 0 gets closer to the BH event horizon, the emitted energy E inf l goes to zero as expected. The plot on the right uses Eq. (35) with B = 0. It can be Table 1 Spectra of the energy radiated to infinity in the ZFL for multipole numbers l = 0, 1, 2, considering the source released from infinity with non-vanishing initial velocity. We compare the results obtained numerically by using Eq. (10) with the ones obtained analytically in Eq. (43)  seen that the emission is mainly with l = 0 and l = 1 and that it increases as a function of r 0 . The plots on the left use Eq. (34). The spectrum does not decrease as a function of l. This is an ultraviolet effect arising from a sudden emergence of a source at r = r 0 .
In Tables 1, 2 and 3, we compare the zero-frequency limit (ZFL) of the energy spectra obtained numerically with the corresponding analytic results. In Table 1 we compare the ZFL of the spectrum of the energy radiated to infinity from a source infalling from r = ∞ obtained analytically in Eq. (43) with the ones obtained using the numerical solution of Eq. (10) in order to check the numerical results. Similarly, in Table 2 we compare the ZFL for the spectrum of the energy radiated to infinity for the source infalling radially from rest from a certain position r 0 , given by Eq. (45), with the ones obtained using the numerical solution of Eq. (10).
In Table 3 we compare the ZFL for the spectrum of the energy radiated to the event horizon for the source infalling radially from rest at a certain position r 0 , given by Eq. (46), with the ones obtained using the numerical solution of  Eq. (10). In all cases the ZFL of the numerical results agree very well with the analytic expressions.

Summary
In this paper we studied the radiation emission of massless scalar field from the radial infall of a source into a Schwarzschild BH, using the formalism of QFT at tree level. We numerically computed the spectra of the emitted energy, E hor ω l (energy radiated to the event horizon) and E inf ω l (energy radiated to infinity), in two distinct situations related to the initial condition of the radial infall of a source, namely: (i) the source starting with non-vanishing velocity v 0 from spatial infinity and (ii) the source infalling from rest at a certain finite position r 0 .
Some aspects of the case in which the source comes from infinity with a certain non-vanishing velocity v 0 are summarized here. For all multipole numbers l the spectra of the energy emitted to infinity are nonzero in the low-frequency limit. For high initial velocities the spectrum is approximately constant until the frequency is around the fundamental quasinormal frequency and then rapidly goes to zero. For low initial velocities the monopole (l = 0) radiation and the dipole (l = 1) (to a lesser extent) are dominant while higher multipoles are significant at higher initial velocities. Interestingly, as the initial velocity v 0 is increased, the emitted energy for the monopole radiation (l = 0) decreases whereas that for the higher multipoles increases. The spectrum of the energy emitted to the event horizon starts from zero increases to a maximum and then decreases very slowly, reflecting the fact that, as the source falls toward the horizon, the Coulomb-like energy of the source passes through the black hole.
Next, some properties of the energy spectra when the source starts falling from rest at a certain distance from the black hole are being discussed. The monopole spectrum of the energy emitted to infinity starts from a nonzero value, while for l ≥ 1 the spectra start at zero. These spectra all rise to a maximum and then decrease to zero. As for the spectrum of energy emitted to the event horizon, the emission is mainly with l = 0 (monopole) and with l = 1 (dipole). The emission to the horizon has more contribution from higher multipoles. It is also remarked that a naïve calculation would lead to a source appearing abruptly and then starting to fall. A boundary term needs to be subtracted in order to calculate the emission from a source at rest and then starting to fall.
Finally, it is noted that the behavior of the spectra of the emitted energy, E inf ω l , and the total emitted energy E inf l to infinity for multipole numbers l ≥ 1 are similar to the electromagnetic and gravitational cases (see e.g. Refs. [7,8,10,13]).