Path-integral derivation of black-hole radiance inside the de-Rham-Gabadadze-Tolley formulation of massive gravity

If we apply the path integral formulation in order to analyze the particle creation process of black-holes inside the non-linear formulation of massive gravity, it is possible to demonstrate that the effect of the extra-degrees of freedom is to deform the periodicity of the poles of the propagator in the complex $t$-plane. This might create the effect of extra-particle creation process at scales where the extra-degrees of freedom become relevant. For stationary solutions, depending on the values taken by the free parameters of the theory, the periodicity structure of the propagator reveal two effects. The first one is a shift on the positions of the pole of the propagator with respect to the GR case, affecting then the instant at which the particles are detected. The second one is the existence of branch points, affecting then the perception of particles. The branch point can be finite (including the zero order case) or infinite depending on the free-parameters of the theory.


I. INTRODUCTION
In [1,2], the first analysis of particle creation process of black-holes in dRGT massive gravity was performed. The analysis done in [1] used the Hawking method studied in [3], but in an implicit way. This means that the extra-degrees of freedom were not studied explicitly, but rather they appeared implicitly inside the definitions of extended coordinates (Stückelberg functions). This method just take the advantage of the fact that the extradegrees of freedom appear in a similar way as the coordinate transformations in general relativity (GR). This is just the well known Stückelberg trick applied at the non-linear level. In [2], the path integral formulation was explored. In such a case, the analysis of the periodicity of the poles of the propagator revealed the possibility of an extra component of radiation coming from the extra-degrees of freedom of the theory. This extra-component, is relevant at scales where the extra-degrees of freedom are important. If we take the origin of coordinates on the source, this happens for distances larger than the Vainshtein radius r V and only observers defining the time in an arbitrary direction with respect to the Stückelberg function (T 0 (r, t)) are able to detect it. On the other hand, observers defining the time in agreement with T 0 (r, t) will not be able to perceive the additional contribution to the Hawking result and all their measurements of temperature and particle creation process are reduced to the GR case. Although the result found in [2] is general, in such a case the author did not evaluate explicitly the new value of temperature. This value will depend on how the periodicity of the propagator is affected by the presence of the extra-degrees of freedom. The calculations performed by using the path integral method, are simplified if we work by defining the Stückelberg functions in agreement with the extended versions of the advanced and retarded coordinates as has been done in [1,2]. In the extended version of the coordinates, the Stückelberg function denoted by T 0 (r, t) appears in replacement of the time-coordinate. In this manuscript, I study the particle creation process of black-holes by using the path-integral formulation. I explore the analyticity properties of the propagator inside the non-linear formulation of massive gravity proposed by de-Rham, Gabadadze and Tolley (dRGT). The general result suggests that for Stückelberg functions given by T 0 (r, t) = St + A(r, t), with A(r, t) having some non-linear time-dependence, the periodicity properties of the poles of the propagator change and the effect of the extraparticle creation is evident. The dynamical metric in this case is however, non-stationary. If A(r, t) has a linear time-dependence, then the properties of the poles of the propagator are also affected and the stationary condition of the dynamical metric would be satisfied. These general situations are not explored in this manuscript. Here I concentrate the analysis on the black-hole solution found in [4] for the case β = (3/4)α 2 , with β and α representing two free parameters of the theory. This combination of parameters is special because it provides a solution with zero cosmological constant (Λ) and a non-zero graviton mass. Then it is easier to analyze the effect of the extra-degrees of freedom at large scales. Similar solution was found in [5]. The stationary condition of this solution, guarantees that the Stückelberg function is linear in time and given by T 0 (r, t) = St + A(r), with A(r) representing a spatial dependent function carrying the information of the extra-degrees of freedom. Here A(r) is time-independent. Then after analytical extension of the time-coordinate t, the effect of the function A(r), is to create a shift effect for the poles of the propagators if we compare them with respect to the GR case. In addition, if we express the extended version of the Kruskal coordinates as an explicit function of time t, then after the analytical extension of the time-coordinate, the periodicity of the poles of the propagator is given by 8πM/S 2 , with S depending on the free-parameters of the theory. The relation M/S 2 marks the possibility of having a branch point condition. The order of the branch point is given by 1/S 2 and it depends on the value taken by the parameters of the theory. For the case β = (3/4)α 2 , the order of the branch point becomes zero when α → −2. It is infinite for α = 0 and in general is finite for other values taken by the parameter α. If S = 1 (when α → ±∞), the periodicity of the propagator is the same as in the GR case, however the shift of the poles of the propagators produced by the function A(r) in this particular case, although does not affect the value of temperature with respect to the GR case (for the case S = 1), it affects the space-time locations where the observers detect particles. In other words, if an observer moving under the theory of GR travels through a line of constant r, θ, φ, then he will disagree with the instants (of time) at which the particles appear if he compare his results with the observers describing the physics in agreement with the dRGT theory of massive gravity. This is equivalent to a time-dilation effect reproduced by the extra-degrees of freedom at large scales. For the case of the branch points, when they become of order zero or infinite order, the effect of particle creation coming from the black-hole is completely suppressed for the observers describing this kind of solutions. The paper is organized as follows: In Section (II), I make a brief review of the Schwarzschild de-Sitter solution inside the non-linear formulation of massive gravity; this solution is discussed in detail in [4]. In Section (III), I make a brief review of the path integral formulation applied to the analysis of the particle creation process of black-holes. This analysis was originally proposed by Hartle and Hawking [6]. Section (IV) is dedicated to explore the non-linear version of the Stückelberg trick. The idea is to explain the distinction between coordinate transformations and the Stückelberg trick. From the perspective of GR, the Stückelberg trick looks like a dipheomorphism transformation. Inside dRGT massive gravity, this is not the case. In Section (V), I explore the analyticity properties of the propagator. This is the key section of the manuscript. In Section (VI), I explain the expected modifications for the black-hole radiance due to the extra-degrees of freedom of the theory for the case β = (3/4)α 2 . This section explores the connection between the rates of emission and absorption. Finally, in Section (VII), I conclude.

II. THE SCHWARZSCHILD DE-SITTER SOLUTION IN DRGT: UNITARY GAUGE
In [4], the S-dS solution can be written explicitly as: where: and f (Sr) = 1 − 2GM Sr − 1 3 Λ(Sr) 2 , with S = α α+1 being the scale factor which depends on the free-parameters of the theory. In this previous solution, all the degrees of freedom are inside the dynamical metric. The fiducial metric in this case is just the Minkowskian one given explicitly as: The solution (2) can be equivalently written as: where T 0 (r, t) corresponds to the Stückelberg function. This function contains the information of the extra degrees of freedom in agreement with the formulation introduced in [1,7]. The metric (4) is then diffeomorphism invariant [1,7]. The gravitational degrees of freedom of the Stückelberg fields appear inside T 0 (r, t) through the spatial (temporal) dependence of the function. For the family of solutions with two free-parameters and the Stückelberg function constrained as has been found in [4], the cosmological constant is given by [4]: This constant is zero when the two free-parameters of the theory α and β satisfy the relation β = (3/4)α 2 . Independent of the relation between the two free parameters of the theory, the Stückelberg function has to satisfy the constraint [4]: A global solution of this previous constraint, is given by the Finkelstein-type form [4]:

III. THE POSSIBLE PATHS FOR A PARTICLE COMING FROM THE BLACK-HOLE
The path integral method in GR, was introduced by Hartle and Hawking in [6]. Here I will use the same method in order to derive the results of the Hawking radiation in massive gravity and then analyze the possible contributions coming from the extra-degrees of freedom. In the Feynman path integral method, the amplitude K(x, x ) for a particle to propagate from a point x to another point x , is given by: where S(x, x ) is just the classical action for a particular path connecting x and x . The amplitude K(x, x ) is called propagator. In Fig. (1), we can observe the analytically extended Schwarzschild space. The shaded part of the diagram should be replaced by a gravitationally collapsing object. The path BCA, represents the process of a pair of particles being created at the point C. One of the particles propagates forward in time and reach the point A, where an observer with a detector is located. The other particle propagate backwards in time, reaching eventually the future singularity at B. There are two kind of observers in dRGT massive gravity. The first one defines the time in agreement with the Stückelberg function T 0 (r, t). This observer will not see any difference between GR and the dRGT formulation of massive gravity. This observer cannot detect any contribution coming from the extra-degrees of freedom of the theory. The second kind of observer, defines the time arbitrarily in agreement with t. In general, the directions of the Killing vectors defined in agreement with T 0 (r, t) and t, will disagree. When this happens, we have the possibility of extra-particle creation process due to the extra-degrees of freedom of the theory. In general metric solutions, the extra-particle creation process will emerge from the time-dependence (t) of the Stückelberg function. In the stationary case with A(r) being time-independent, which is the one explored in this manuscript, the Stückelberg function is linear with respect to t and there are two effects that will appear still at this level. The first one is the branch point condition reproduced by the factor S appearing in the metric components (2). The second effect is the shift of the poles of the propagator with respect to the GR case, even if S = 1. Later in this manuscript, I will explain in detail both effects. For the moment, I will define the path integral method and I will explain how can it be applied to the dRGT massive gravity case. Here I do not consider the paths coming from I − , because they correspond to propagation of incoming particles from the infinite past. Here I will assume that the infinite past does not contain any degree of freedom. The paths corresponding to the shaded region are ignored, since they correspond to the collapsing object itself. The propagation is then considered from the future singularity. If we want to use the stationary phase method in order to derive the particle creation process, we immediately find that it is impossible to find real stationary paths connecting the future singularity with a positive-frequency mode for a stationary exterior observer [6]. Then we have to analytically extend the path toward the past singularity, or equivalently, we have to translate the point B on the future singularity to complex values. This is equivalent to relate the rate of emission to the rate of absorption of the black-hole after thermodynamic equilibrium is reached. Mathematically, the connection is [6]: where N (E) and P (E) correspond to the probability emission and absorption respectively.
Here κ is the surface gravity of the black-hole, which in the standard case of GR becomes κ = 1/4GM . In massive gravity, some modification is expected for observers located at scales larger than r V if they define their corresponding time-like Killing vectors with respect to the usual notion of time coordinate. If we want to apply the path integral method for the analysis of black-holes in GR to the dRGT case, the Stückelberg function plays the role of the time-coordinate as has been explained in [2]. Then the time-coordinate t, will appear implicitly inside this function. This is equivalent to use the Stückelberg trick at the non-linear level explained in the coming section.

IV. THE NON-LINEAR STÜCKELBERG TRICK
The coming section will require to use the Stückelberg trick in order to reproduce the important results related to the Hawking radiation inside the non-linear formulation of massive gravity. At the non-linear level, the trick looks like the dipheomorphism transformations in GR. In unitary gauge (with a fiducial metric being Minkowski), it can be expressed as [7]: where Y α are the components of the Stückelberg functions. The previous equation looks like a standard gauge transformation in GR. It is however, the way of introducing redundant variables in order to restore the dipheomorphism invariance of the massive action in massive gravity. The redundant variables contain the information of the extra-gravitational degrees of freedom. Under the dipheomorphism transformations given by: the metric (10) is invariant. This is necessary for keeping the massive action invariant. Infinitesimally, the Stückelberg function expansion: provides the following result on the metric (10): The gauge transformations (11), after the infinitesimal expansion f (x) = x + ζ(x), become: The A µ -term corresponds to the Goldstone bosons that at the non-linear level carry the broken symmetry in massive gravity. It can be verified again that under the previous infinitesimal gauge transformations, eq. (13) provides the result: which demonstrates that the dynamical metric with the Stückelberg functions appearing explicitly is dipheomorphism invariant.

V. ANALYTICITY PROPERTIES OF THE PROPAGATOR
If we want to derive the analyticity properties of the propagator, the easiest way is by using the Stückelberg trick, already explained in the previous section. The trick permits us to define the Stückelberg function in terms of new objects which can be considered as extended coordinates. The basic idea is to replace the standard time-coordinate t by the use of the Stückelberg function T 0 (r, t) inside this new objects. We can then define the components of the Stückelberg functions as extended coordinates. Here we start with the basic differential definition of the propagator given by: with the appropriate boundary conditions. In dRGT massive gravity, the condition β = (3/4)α 2 is equivalent to the Schwarzschild solution (asymptotically flat), but with some background of Stückelberg fields. Inside the Schwarzschild geometry, we will consider the case where x is external to the black-hole and x is over the horizon as in the standard case. If the Vainshtein mechanism operates, as it is expected from these solutions, then the extra-degrees of freedom effects will become relevant after the Vainshtein scale and then the time coordinate orientation selected by the observer becomes important for the vacuum definition and as a consequence, for the definition of temperature. As has been explained before, for observers defining the time in agreement with T 0 (r, t), the perceived radiation will not change with respect to the GR case [1,2]. On the other hand, for observers defining the time in the standard way, the effect of the extra-degrees of freedom will appear and they will generate changes in the periodicity pattern for the poles of the propagator. We can define the Kruskal coordinates in dRGT by using the Stückelberg trick. They are defined as [1,2]: with: where the Stückelberg fields components are defined as: and T 0 (r, t) is the initial Stückelberg function. The objects (20) and (21) are the new Stückelberg functions. They can also be defined as extended coordinates inside the dRGT formulation. Note that the condition (19) is not affected by the presence of the extra-degrees of freedom of the theory (dRGT) because the initial Stückelberg function (T 0 (r, t)) does not appear in the final result. In fact, whenever we find time-independent quantities defined originally inside the framework of GR, when extended to the dRGT case, will not be affected by the presence of the extra-degrees of freedom in spherically symmetric configurations.
Only time-dependent quantities as they are defined originally inside the framework of GR, will suffer some changes when they are extended to the dRGT case. The definition of event horizon taken from the combination of Stückelberg function components in eq. (19) is unchanged with respect to the GR case. Then the event horizon is defined with the condition U = 0 or V = 0, when r = 2GM . This can also be verified from the metric (2), then the event-horizon position will not be affected by the presence of the extra-degrees of freedom. The complexified horizon in this case, is defined as the surface where the non-zero of the coordinates U or V is extended analytically to complex values. Any singularities for the propagator will come from the point W = 0 as in the standard case [6]. By dividing the interval of W as in [6], the propagator becomes: where K 0 (x, x ) corresponds to the propagator of the interval [0, W 0 ] for the parameter W with W 0 being infinitesimal and near the divergence point of the propagator. We can observe that there is a singularity as s c (x, x ) = −i . In massive gravity, the extra-degrees of freedom are contained inside this condition if we take into account the Stückelberg trick as has been formulated in eq. (10). Then in terms of the Stückelberg function, this result corresponds to a null geodesic, connecting x with the complexified manifold. When expressed explicitly in terms of the standard time coordinate t, this previous null condition contains some extraterms (Stückelberg fields) which will modify the analyticity properties of the propagator. All the geodesics starting from real values of x , will intersect the horizon at real sections, namely, sections where the extended versions of V and U are real. We can consider the geodesics which connect a real value of x which is outside the black-hole with the future event horizon. Then we can define again new Stückelberg functions by extending the notion of the advanced Eddington-Finkelstein coordinates in agreement with (20), covering V ≥ 0: with V representing the extended version (Stückelberg function) of the advanced-null coordinates V = v + A(r, v). Here T 0 (r, t) = St + A(r, v). At this point of the calculation, we have two options. The first one is to work by using V as the Killing time. In such a case, we will obtain exactly the same results of GR. This corresponds to the observers defining the vacuum in agreement with T 0 (r, t). The second option, is to use v = V − A(r, t) as the Killing time, defining then eq. (23) as: It is possible for some observers to define the vacuum in agreement with v, namely, the usual time-coordinate. By complexifying the geodesics assuming an imaginary affine parameter λ, then the geodesic would be a two-dimensional sheet in the complex coordinate system and then from eq. (23), we know that for complex extensions of V , there is a periodicity associated to the coordinate V and given by Imκ = 2π/V or: Note however that if we use the coordinate v as the Killing time, then any dependence on v coming from the function A(r, v), will affect the periodicity properties of the complexified variable v. For stationary solutions with A(r) being independent of v, like the one explored in this manuscript, the periodicity structure of the propagator still can be affected by a couple of effects that can appear due to the presence of extra-degrees of freedom. These effects will be explained later in detail. What is important for the moment is to note that the periodicity associated with the variable V , can be different to the one associated to the variable v. By convenience, we select the strip which contains real values of V . We suppose that for some affine parameter λ = 0, the coordinates (and Stückelberg functions) take the values V 0 , r 0 , θ 0 , φ 0 , which are assumed to be real. Then we can study which complex values of the extended coordinate (Stückelberg function) V in the strip with real values of r = 2GM , θ and φ are contained inside the two-dimensional sheet generated by the null geodesics. It is clear that the two-dimensional sheet generated by the complex values of V , when expressed in terms of v, will in general contain some arbitrary combinations of the complex values of v and A(r, v), depending on the specific form of the Stückelberg function T 0 (r, v). By extending the method proposed by Hartle and Hawking [6], the singularities of the propagator on the complexified horizon are slightly displaced from the real values of U' and V'. They correspond to the poles marked by s(x, x ) = −i , which are infinitesimal displacements with respect to the real values on the complexified horizon. The direction of the displacements, can be found in the same way as Hartle and Hawking did [6]. The analyticity properties of the propagator suggests that the periodicity of the poles will change with respect to the GR case for observers taking the time in agreement with t defined in an arbitrary direction with respect to T 0 (r, t). Observers defining the time in agreement with T 0 (r, t), will not perceive any change in the periodicity of the poles of the propagator. For these special observers, there is no difference with respect to the GR case. We can define Stückelberg functions in agreement with the extended null-coordinates U and V defined as: with the (initial) Stückelberg function T 0 (r, t) playing the role of an extended notion of time.
The results (26) are valid for the cases U > 0 and V > 0. On the other hand, we also have: for the cases U < 0 and V > 0. Analogous relations can be found for the other quadrants of the extended Penrose diagram. In the diagram (2), the region II has the Cauchy data over part of the future horizon and part of the past horizon. Then the initial data over these surfaces, determine uniquely the propagator inside this region. If we complexify the original Stückelberg function T 0 (r, t), and keeping the coordinates r, θ, φ real, then eqns. (26) can be written in the modulus-argument form: where T 0 (r, t) = γ(r, t)+iψ(r, t). Then the problem of determining the propagator is reduced to solve the wave equation (17), in the real functions |U | and |V | for fixed values of ψ(r, t). This is analogous to the standard case illustrated by Hawking and Hartle [6]. However, the main difference here is the fact that the function ψ(r, t), contains the information of the extra-degrees of freedom of the theory. The propagator is analytic on the complexified horizon on the upper-half plane of the complex variable U and on the lower half plane for the complex variable V . The Cauchy data is regular when: where ψ(r, t) = µ +Ā(r, t), withĀ(r, t) being the imaginary component of A(r, t). If we know the exact dependence of the function A(r, t), then we can separate easily the real part from the imaginary one. The extra-component of radiation, coming from the deformation of the periodicity pattern of the propagator produced by the Stückelberg function T 0 (r, t), depends on whether or not there is an imaginary component for A(r, t) after doing the analytical extension on the variable t. The observers defining the time in agreement with T 0 (r, t), will never be able to perceive any extra-component of radiation and they will define the same result of temperature as in GR. On the other hand, observers defining the time in an arbitrary direction with respect to T 0 (r, t), might be able to perceive an extra-component of radiation, depending on the explicit form of the Stückelberg function. The propagator is analytic in T 0 (r, t) in the same form as it appears in GR. The functionsŪ andV satisfy the Cauchy-Riemann (orthogonality) conditions, already demonstrated in [6]. Then given the previous arguments, the propagator K(x, x ) will be analytic in T 0 (r, t) in a strip of width 4πM . Notice however, that this means that: If we want to solve for µ this previous inequality, we have to take into account that A(r, t) can be time-dependent. This means that this function can be written as: where we have defined t = γ + iµ. The exact range of regularity for the propagator given by eq. (30), and solved in terms of the standard complex time coordinate will depend on the exact form of the function A(r, t). Only if this function is time-independent, then it is possible to solve trivially for µ in eq. (30) as follows: − 4πGM −Ā(r) < µ < −Ā(r). (32) From these previous results, we can observe that the periodicity of the poles of the propagator might change due to the presence of the extra-degrees of freedom. In other words, the periodicity associated to the variable t, is not necessarily the same periodicity associated to the function T 0 (r, t). The exact behavior of the poles of the propagator with respect to any observer, depends on the exact functional dependence of the Stückelberg function T 0 (r, t).
The modified period of the propagator, will reproduce a change on the amount of particles detected if the detectors operate by using the usual notion of time t. The Vainshtein radius marks the scale after which the extra-degrees of freedom become relevant, then it is expected that the effects explained in this manuscript are relevant after the Vainshtein radius [1]. If we use the Stückelberg function as the time-coordinate, then the Penrose diagram for the case with zero cosmological constant will not differ from the standard Schwarzschild diagram shown in Fig. (1) [1]. Then the propagator K(x, x ) with x in the region II of the diagram (2), still can be continued to T 0 (r, t) → T 0 (r, t) − i4πM . This analytical extension, simply relate the propagator connecting the points x and x , with the one which connects the points x and x in the diagram (2). Note that x is in the region III. The Penrose diagram (2), contains implicitly the extra-degrees of freedom of the theory. If the observers define the time in agreement with T 0 (r, t), then everything will look like GR. These special kind of observers will never perceive any difference between massive gravity and GR and they will describe the physics in agreement with the standard Penrose diagram. The analytical extension just mentioned, marks the relation between the rates of emission and absorption of the blackhole. If we want to understand the physical effects of the extra-degrees of freedom, then we have to express the propagator in terms of the ordinary time t. Then for detectors operating with respect to t, instead of T 0 (r, t), the analytic extension does not necessarily relates the rate of emission with the rate of absorption of the black-hole radiation. This is because the extension is related not only to the complex part of the ordinary time variable t, but also to the complex part of the function A(r, t), namelyĀ(r, t). The analyticity of the propagator is periodic as in the standard case with period 8πGM for observers describing the physics with respect to T 0 (r, t) [6]. For the case of observers describing the physics with respect to the usual time t, it is not clear at all what is the periodicity associated to the complex time. Later I will demonstrate that for the case under study in this manuscript, the periodicity structure of the propagator, depends on the values of the free-parameters of the theory.
A. The analyticity of the propagator: The case β = (3/4)α 2 The explicit Stückelberg function T 0 (r, t), can be found for the case β = (3/4)α 2 from eq. (7), the result after integration is: This result looks similar to the tortoise coordinates introduced in GR. The argument of the natural Logarithm has a different behavior in this case however. The scale at which the logarithm diverges, depends on the exact value of S. The scale factor S is a function of the free-parameters of the theory [4]. If we compare the result (33) with the definition of tortoise coordinate given by: we immediately notice the similarity in both cases. Note that if r = 0 (The singularity point), the expression (33) becomes T 0 (r, t) = St, this means that when the gravitational field becomes stronger, the extra-degrees of freedom become negligible as it should be if the Vainshtein mechanism operates. Then at scales close to the source, GR is approximately recovered. The scale factor S, as a solution of the parameters of the theory and for the case β = (3/4)α 2 can be expressed as: Note that this function becomes equivalent to: Is interesting to observe that as α → ∞ or α → −∞, then S → 1. Note that S is a well behaved function, except when α = −2 which represents a point of discontinuity. The figure  (3) illustrate the behavior of S as a function of α. If we select the parameter α to be very small, then one of the two possible branches of S (see fig. (3)) is also small and the distance where the Stückelberg function (33) diverges can be very large. The scale of divergence is given by: which is related to the scale of the black-hole. When we analyze the Hawking radiation, usually we use the advanced or retarded coordinates in GR. In massive gravity, we can define in analogous way, the Stückelberg functions as extended advanced and retarded coordinates in agreement with: where the re-scaling factor S has been used. Here there are two differences with respect to the GR case; the first one is the replacement of the time coordinate by the Stückelberg function T 0 (r, t). The second one is the inclusion of the factor S re-scaling the distance coordinate r. We then get different results depending on the combination of signs between T 0 (r, t) and r * defined in agreement with eq. (34) but taking into account the re-scaling factor S. The following are the possible cases obtained when we extend the notions of advanced and retarded coordinates in order to define the Stückelberg functions: B. Case I: C. Case II: and combinations of these previous results. The logarithmic argument in eqns. (40) and (43) is positive defined when: Otherwise the Logarithmic function provides a complex result. We can do a new re-definition of the Stückelberg function as follows: taking into accont the implicit appearance of the extra-degrees of freedom, then from eqns.
E. Case II: and combination of these cases. If we analytically extend the time-coordinate to complex values, then we have t = τ + iµ. And then the Stückelberg functions defined in agreement with V and U can be written in the form: for any of the cases illustrated before. The Cauchy data is then regular when the condition: is satisfied. If we compare this result with the notation used in eq. (32), then the periodicity of the propagator is related to the complex function: Then for the case β = (3/4)α 2 , the observers defining the physics in agreement with t, will perceive two effects, the first one is the shift of the period with respect to GR, even if S = 1. In this case, the periodicity associated to the propagator can be the same as in GR, but the shift of the positions of the poles, create the effects of time-delay for the detection of particles. In other words, if an observer inside the theory of GR has a detector at scales located after r V , he will disagree with the instant of times at which the particles appear with respect to the case of an observer located at the same distant and working inside the dRGT formulation of massive gravity. However, both observers will agree in the average temperature obtained. The second effect inside dRGT is produced by the factor S. This factor opens the possibility of having branch points if it takes a non-integer value. In fact, S 2 represents the order of the branch point with respect to the GR case. Then the order of the branch point will depend on the parameters of the theory. For the case explored in this manuscript, it will depend on α. Note that for α = −2, the branch point is of zero order and the detection of particles is suppressed. If α = ±∞, then the branch point condition disappears and the theory looks as GR except for the fact of the disagreement of the instant at which the particles appear when we compare the GR case with respect to the dRGT one. Finally, if α = 0, the branch point is of infinite order one and the particle creation process is also suppressed. Fig. (4) illustrate the behavior of the of periodicity structure of the propagator for observers defining the time in agreement with t, with respect to observers defining the time in agreement with T 0 (r, t).

VI. BLACK HOLE RADIANCE
In the GR case, when we analyze the amplitude of emission of a mode with definite positive energy E. Its time-dependence is given by f exp(−iEt). The information about emission is contained in the amplitude ε E , which is defined as [6]: In this case however, the ordinary time coordinate has been replaced by the Stückelberg function T 0 (r, t). Here, ε E is the component of energy, coming in principle from the surface C + as in the GR case. However, in massive gravity the effect of the extra-degrees of freedom will appear through the Stückelberg function T 0 (r, t). Inside the GR formulation, C + is the only surface which contributes to the radiation perceived by the detector located at large scales. In massive gravity however, since the extra-degrees of freedom appear, it is possible that some other surfaces can contribute. Here however, we can still assume that the pre-collapse surface does not contribute to the final result. This means that for a detector confined to work inside an interval (−t 1 , t 1 ), with t 1 very large, still we can assume the absence of particles in the infinite past. In massive gravity however, we would be talking about an interval defined by the Stückelberg function. Then for some observers, the pre-collapse surface given at the time T 0 (r, t) = t 1 → −∞, would be equivalent to St + A(r, t) → −∞. Then the assumption of a vacuum devoid of particles in the spacelike pre-collapse surface for any observer, is equivalent to say that the extra-degrees of freedom are also absent there. Even if still the contributions due to the extra-degrees of freedom will not appear from the pre-collapsing surface, they might appear from the timelike surface connecting the space-like surface inside the future horizon with the observation surface, where the detectors are located. When the detectors are located at scales larger than r V , the contribution coming from the extra-degrees of freedom becomes relevant. For an observer located at scales larger than r V and defining the time arbitrarily with respect to T 0 (r, t), the perceived effect is equivalent to some extra-modes not necessarily coming from the event horizon of the black-hole. In the spherically symmetric solution in massive gravity found in [4], the background solution has a symmetry under time-translations. The quantity conserved under time-translations is however a combination of the usual notion of energy inside GR, with a velocity-dependent quantity associated to the extra-degrees of freedom [2,8,9]. This conserved quantity defines the notion of energy in agreement with the observers defining the time with respect to the Stückelberg function. Then the propagator will be a function of T 0 (r, t) − T 0 (r, t) at large values of T 0 (r, t). Observers defining the time in agreement with T 0 (r, t), will not perceive any contribution coming from the time-like surface connecting the space-like surface inside the future event horizon, with the space-like surface of observation. On the other hand, the observers defining the time arbitrarily, will perceive some (possible) extra-contribution. The extra-contribution will appear as a change in the periodicity of the propagator. The extra-component of radiation, will become relevant at scales larger than r V . Now we can relate the rate of emission with the rate of absorption. For the detectors defining the time in agreement with T 0 (r, t), the result is trivial and it is the same to the one obtained in GR [6]. In this trivial situation, by distorting the contours of integration in T 0 → T 0 − 4πM i, then eq. (53), becomes: If we however expand eq. (54) by using the usual notion of time and separating real and imaginary contributions of the Stückelberg function in agreement with T 0 (r, t) = St + iµ + A(r, t) + iĀ(r, t), then there is no guarantee that the previous result relates the rate of emission with the rate of absorption for the standard Penrose diagram, namely, the one written with respect to the usual notion of time t. Then in such a case, the extended relation between the rates of emission and absorption, not only includes the radiation coming from the event horizon, but also some extra-component appearing at scales larger than r V . This extra contribution comes from the time-like surface connecting C + with the space-like observation surface. Similar conclusions were found in [1].
A. The radiation with respect to observers defining arbitrary time When the degrees of freedom inside the dynamical metric appear explicitly, and the detectors are peaked in agreement with the usual time-coordinate, then some extra-contribution to the radiation process might appear at scales larger than r V . In agreement with the previous analysis of the analyticity of the propagator, there are poles when s(x, x ) = −i . An observer defining the time in agreement with the Stückelberg function, will define the positions of the pole propagator in agreement with the expression [6]: In terms of the usual notion of time, this is equivalent to: In the case under study here, this previous condition also corresponds to the singularities of the propagator. The singularities are periodic in agreement with the analysis of the previous section. It is the periodicity of the poles what reproduces the effect of particle creation. If we imagine for a moment the function T 0 (r, t) = St + A(r, t), with A(r, t) arbitrary (initially), then the analyticity of the propagator reveals that the period associated to the Stückelberg function T 0 (r, t) is not necessarily the same period associated to the usual time-coordinate. This is expected by intuition if we observe that in agreement with the expressions (55) and (56), the periodicity of the Stückelberg function T 0 (r, t), has to be divided by two: 1). One part belonging to the complex time variable µ as has been defined previously. 2). The other fraction going to the complex part the functionĀ(r, t). This phenomena will in general change the perception of temperature for the detectors operating under the usual notion of time t because the relation between the emission and absorption of the black-hole radiation for the stationary situation is not trivial anymore. For the stationary case analyzed in this manuscript with the parameter combination β = (3/4)α 2 , the relation between the emission and absorption, can be trivially derived by using eq. (54). The extension in this case is trivial because of the time-independence of the function A(r) contained inside the original Stückelberg function. The trick now is to express eq. (54) explicitly in terms of the usual time coordinate. In general however, the result is not so trivial if the function A(r, t) is time-dependent. The reason is that the delta function coming from the time-integral and representing the conservation of energy will not appear and then there is no justification for using the simple expression (53) in this general situations [6]. However, if A(r, t) is timeindependent, then still there might appear a delta-function, and then the modification to the black-hole temperature as it is perceived by the observer working with the usual notion of time t is trivially found. Then we get: where τ represents the modified period. The modification of the period with respect to the GR case will depend on the extra-degrees of freedom as has been explained previously. It is clear from eq. (57), that the correspondence between the rates of emission and absorption related to the variable T 0 (r, t), does not necessarily implies the same correspondence for the rates of emission and absorption described with respect to t. By using these previous arguments, the temperature perceived by a detector using as a reference the time t is given by: This is in agreement with the result found in [2] where however, the modified period was expressed in terms of the complex component of the function A(r). For the case explored in this manuscript, if we use the results of Sec. (V A), the exact value of temperature is given by: This value of temperature can be higher or smaller than the GR case, depending on the relation S 2 /M , which also depends on the free-parameeters of the theory.

VII. CONCLUSIONS
In this manuscript, I explained the particle creation process of black-holes in dRGT non-linear massive gravity by using the path integral formulation. The first attempt for explaining the particle creation process in dRGT, was done in [1] by the author. In this manuscript, tha author focused in the asymptotically flat case with β = (3/4)α 2 . This case represents the situation of zero cosmological constant with a non-zero graviton mass. There is however, a non-trivial configuration of Stückelberg fields around the black-hole. They are represented by the Stückelberg function T 0 (r, t). By using the Stückelberg trick, it was possible to extend the notion of Kruskal coordinates in order to analyze the periodicity properties of the propagator. It was demonstrated that for the observers defining the time in agreement with the Stückelberg function T 0 (r, t), the periodicity properties of the propagator are the same as in the GR case, then these special observers will not perceive any difference between GR and dRGT at the moment of measuring the temperature. On the other hand, observers defining the time in agreement with t, will be able to perceive some extra-component of radiation coming from the extra-degrees of freedom. If the extra-degrees of freedom are relevant after the Vainshtein scale, then only the observers located at scales larger than r V will be able to perceive the new effects. These observers will experience one of the following situations: 1). Branch point conditions which for the cases of α → 0 or α → −2, suppress any particle creation process. For other cases, with α arbitrary, the perception of particles is affected due to the distortion on the period of the propagator. 2). Shift effect for the case S = 1 (α → ∞). In this case, although the period of the propagator is the same with respect to the GR case for any observer, the instants at which the particles appear are different if observers working under the theory of GR compare their results with the ones obtained by observers working under the theory of dRGT.