Gravitational luminosity of a hot plasma in R^2 gravity

The R^{2}-gravity contribution to energy loss of a hot plasma due to the gravitational bremsstrahlung is calculated in the linearized theory on the basis of classical Coulomb scattering of plasma constituents in small-angle scattering approximation. The explicit dependence of the gravitational luminosity on the plasma temperature is derived and its relevance to the Einstein gravity is demonstrated. The result when applied to the Sun as a hot plasma, shows very good agreement with available data.


Introduction
For a long time, the discovery of gravitational wave (GW) emissions from the compact binary system with two neutron stars PSR1913+16 [1] has been the ultimate motivation for the design, implementation, and advancement of extremely sophisticated GW detection technology. Physicists working in this field of research need this technology to conduct thorough investigations of GWs in order to advance science. The observation of GWs from a binary black hole (BH) merger (event GW150914) [2], which occurred in the 100th anniversary of Albert Einstein's prediction of GWs [3], has recently shown that this ambitious challenge has been won. The event GW150914 represented a cornerstone for science and for gravitational physics in particular. In fact, this remarkable event equipped scientists with the means to give definitive proof of the existence of GWs, the existence of BHs having mass greater than 25 solar masses and the existence of binary systems of BHs which coalesce in a time less than the age of the Universe [2]. After the event GW150914, LIGO detected a second burst of GWs from merging BHs, the event GW151226 [4]. The great hope is that such detections, also through the collaboration with other detectors [5,6], will soon become routine and part of a nascent GW astronomy. GW astronomy will be important for a better knowledge of the Universe and also to confirm or to rule out the physical consistency of the general theory of relativity (GTR) or of any other theory of gravitation [7]. This is because, in the context of extended theories of gravity (ETG), some differences between the GTR and the others theories can be pointed out starting by the linearized theory of gravity [7]. In this picture, detectors for GWs are in principle sensitive also to a hypotetical scalar component of gravitational radiation, that appears in ETG like scalar-tensor gravity and f (R) theories [7]. Let us clarify some important motivations which lead to a potential extension and generalization of the GTR.
Although Einstein's GTR [8] achieved great success (see for example the opinion of Landau who says that the GTR is, together with quantum field theory, the best scientific theory of all [9]) and withstood many experimental tests, it also displayed many shortcomings and flaws which today make theoreticians question whether it is the definitive theory of gravity, see the reviews [10,11] and references within. As distinct from other field theories, like the electromagnetic theory, the GTR is very difficult to quantize. This fact rules out the possibility of treating gravitation like other quantum theories, and precludes the unification of gravity with other interactions. At the present time, it is not possible to realize a consistent quantum theory of gravity (QTG) which leads to the unification of gravitation with the other forces. From an historical point of view, Einstein believed that, in the path to unification of theories, quantum mechanics had to be subjected to a more general deterministic theory, which he called generalized theory of gravitation, but he did not obtain the final equations of such a theory (see for example the biography of Einstein in [12]). At present, this point of view is partially retrieved by some theorists, starting from the Nobel Laureate G. 't Hooft [13].
However, one has to recall that, during the last 30 years, a strong, critical discussion about both the GTR and quantum mechanics has been undertaken by theoreticians in the scientific community. The first motivation for this historical discussion arises from the fact that one of the most important goals of modern physics is to obtain a theory which could, in principle, show the fundamental interactions as different forms of the same symmetry [10,11]. Considering this point of view, today one observes and tests the results of one or more breaks of symmetry. In this way, it is possible to say that we live in an unsymmetrical world. In the last 60 years, the dominant idea has been that a fundamental description of physical interactions arises from quantum field theory. In this tapestry, different states of a physical system are represented by vectors in a Hilbert space defined in a spacetime, while physical fields are represented by operators (i.e. linear transformations) on such a Hilbert space. The greatest problem is that such a quantum mechanical framework is not consistent with gravitation, because this particular field, i.e the metric h µν , describes both the dynamical aspects of gravity and the spacetime background. In other words, one says that the quantization of dynamical degrees of freedom of the gravitational field is meant to give a quantum-mechanical description of the spacetime. This is an unequalled problem in the context of quantum field theories, because the other theories are founded on a fixed spacetime background, which is treated like a classical continuum. Thus, at the present time, an absolute QTG, which implies a total unification of various interactions has not been obtained. In addition, the GTR assumes a classical description of the matter which is totally inappropriate at subatomic scales, which are the scales of the relic Universe [14,15].
In the unification approaches, from an initial point of view, one assumes that the observed material fields arise from superstructures like Higgs bosons or superstrings which, undergoing phase transitions, generate actual particles. From another point of view, it is assumed that geometry (for example the Ricci curvature scalar R) interacts with material quantum fields generating back-reactions which modify the gravitational action adding interaction terms (examples are high-order terms in the Ricci scalar and/or in the Ricci tensor and non minimal coupling between matter and gravity, see below). Various unification approaches have been suggested, but without palpable observational evidence in a laboratory environment on Earth. Instead, in cosmology, some observational evidences could be achieved with a perturbation approach [15]. Starting from these considerations, one can define as ETG those semi-classical theories where the Lagrangian is modified, in respect of the standard Einstein-Hilbert gravitational Lagrangian, adding high-order terms in the curvature invariants (terms like R 2 , R αβ R αβ , R αβγδ R αβγδ , R R, R k R) or terms with scalar fields nonminimally coupled to geometry (terms like φ 2 R), see [10,11] and references within. In general, one has to emphasize that terms like those are present in all the approaches to the problem of unification between gravity and other interactions. Additionally, from a cosmological point of view, such modifications of the GTR generate inflationary frameworks which are very important as they solve many problems of the standard model of the Universe [14,15,16].
In the general context of cosmological evidence, there are also other considerations which suggest an extension of the GTR. As a matter of fact, the accelerated expansion of the Universe, which is observed today, implies that cosmological dynamics is dominated by the so called Dark Energy, which gives a large negative pressure. This is the standard picture, in which this new ingredient is considered as a source on the right-hand side of the field equations. It should be some form of un-clustered, non-zero vacuum energy which, together with the clustered Dark Matter, drives the global dynamics. This is the so called "concordance model" (ΛCDM) which gives, in agreement with the CMBR, LSS and SNeIa data, a good picture of the observed Universe today, but presents several shortcomings such as the well known "coincidence" and "cosmological constant" problems [17]. An alternative approach is changing the left-hand side of the field equations, to see if the observed cosmic dynamics can be achieved by extending the GTR, see [7,10,11] and references within. In this different context, it is not required to find candidates for Dark Energy and Dark Matter, that, till now, have not been found; only the "observed" ingredients, which are curvature and baryon matter, have to be taken into account. Considering this point of view, one can think that gravity is different at various scales and there is room for alternative theories. In principle, the most popular Dark Energy and Dark Matter models can be achieved considering f (R) theories of gravity, where R is the Ricci curvature [7,10,11]. In this picture, the nascent GW astronomy could, in principle, be important. In fact, a consistent GW astronomy will be the definitive test for the GTR or, alternatively, a strong endorsement for ETG [7].
According to the GTR, a system with a time varying mass moment will loss its energy by radiating the GWs [3,9,18]. This energy loss, at the lowest order, is proportional to the 3th order time derivative of the quadrupole momentum of the mass-energy distribution [18]. In R 2 -gravity, which is the simplest extension of f (R)-gravity, because of the presence of third polarization mode arising from the R 2 curvature term, the situation is different: the extra massive mode contribution leads to an extra energy loss which is proportional to 4th order time-derivative of the quadrupole moment [19].
By comparing the theoretical considerations with the observed decay rate of binary systems PSR B1913+16 [1] and PSR J0348+0432 [20] some constraints on the strength of the R 2 -dependent term is obtained [19,21,22]. In many astrophysical situations, the hot plasma of ionized atoms emits electromagnetic and gravitational radiation through the coulomb collisions between the electrons and ions [23][24][25][26]. Thus, studying the gravitational luminosity of a plasma is of general interest and may be another test for the validity of f (R) theory of gravity. In [27] it has been derived an expression for the amount of radiated energy in a classical gravitational bremsstrahlung in R 2 -gravity, assuming the smallangle scattering approximation. In the present paper, we apply it to derive the gravitational luminosity of a hot plasma with the gravitational bremsstrahlung as a mechanism for the energy loss. In Sec. 2, we linearize the R 2 -gravity theory and, after that, we briefly discuss the quadrupole radiation in R 2 -gravity and energy loss due to gravitational bremsstrahlung in a single Coulomb collision between two charged particles. Sec. 3 is devoted to the calculation of the thermal gravitational radiation of the hydrogen plasma. In Sec. 4 we finally illustrate the correction with an application to the Sun. A summary of the main results is presented in Sec. 5.

Linearized Theory and Quadrupole Radiation in Quadratic Gravity
In the general framework of f (R) gravity, the R 2 theory, which was originally proposed by Starobinski [16], has been analysed in various interesting works, see [28 -32] for example. Specifically, the non-singular behaviour of this class of models is discussed in [28]. In [29] R 2 inflation is combined with the Dark Energy stage and in [30] an oscillating Universe, which is well tuned with some cosmological observations is discussed. Finally, in [31,32] the possibility to partially solve the Dark Matter problem in the linearized R 2 theory has been analyzed.
It is also quite important to emphasize that the R 2 is the simplest one among the class of viable models with R m terms in addition to the Einstein-Hilbert theory. In [33], it has been shown that such models may lead to the (cosmological constant or quintessence) acceleration of the universe as well as an early time era of inflation. Moreover, they seem to pass the Solar System tests, i.e. they have the acceptable newtonian limit, no instabilities and no Brans-Dicke problem (decoupling of the scalar) in the scalar-tensor version.
The field equations of the R 2 -gravity can be derived from the action [28 -33] where L M is the Lagrangian density of matter and a 2 represents the coupling constant of the R 2 term. By varying the action with respect to g µν one obtains with the associated Klein-Gordon equation obtained by taking the trace of Eq.
(2) as where E is known as curvature energy term and defined via E 2 = 1 6a2 [32]. Relation (3) implies the idea of considering the Ricci scalar as an effective scalar field [32].
Before starting the analysis, let us emphasize an important point. As one wants the R 2 -gravity theory to be viable, one needs that it passes the Solar System tests. Thus, one must assume that the constant coupling of the R 2 term in the gravitational action results much minor in respect to the linear term R. In this way, the variation from the standard GTR is very weak and the theory can pass the Solar System tests. Regarding this important issue, there are precedent works illustrating this and we need to explicitly show that the bounds are respected. The key point is that as the effective scalar field arising from curvature is very energetic, then the constant coupling of the R 2 nonlinear term → 0 [34]. In this case, the Ricci curvature, which is an extra dynamical quantity in the metric formalism, must have a range longer than the size of the Solar System. An important work is ref. [35], where it is shown that this is correct if the effective length of the scalar field l is much shorter than the value of 0.2 mm. In such a case, the presence of this effective scalar is hidden from Solar System and terrestrial experiments. Another important test concerns the deflection of light by the Sun. This effect was studied in R 2 gravity by calculating the Feynman amplitudes for photon scattering, and it was found that, to linearized order, this deflection is the same as in the standard GTR [36]. In [32] it has been shown that, in order to partially solve the Dark Matter problem the value of the curvature energy term implies a very low value of the constant coupling of the R 2 term in the gravitational action, that is a 2 ≃ 10 −34 cm 4 in natural units. In that case, the R 2 -gravity theory results viable and l ≪ 0.2 mm is guaranteed. Now, let us proceed to linearize the R 2 -gravity theory. We stress that in the following linearization process we closely follow [32] with a small difference in the definition of the effective scalar field.
Starting frome eq. (3), the identifications [37] Φ → 2a 2 R +a 1 and permit to obtain a Klein -Gordon equation for the effective scalar field Φ as To study GWs, one analyzes the linearized theory in vacuum with a little perturbation of the background, which is assumed given by a a Minkowskian background plus Φ = Φ 0 , that is, one linearizes into a background with constant curvature [32]. One also assumes Φ 0 to be a minimum for V (natural units will be used in the linearization process): where the constant m has mass dimension. Setting to first order in h µν and δΦ, one calls R µνρσ , R µν and R the linearized quantity which correspond to R µνρσ , R µν and R [32]. Thus, one writes down the linearized field equations as [32]: with R µνρσ and eqs. (8) are invariants for gauge transformations [32] h Thus, one defines [32] Let us consider the transformation for the parameter ǫ µ [32] ǫ ν = ∂ µh µν , which permits to choose a gauge analogous to the Lorenz one of electromagnetic waves [38] ∂ µh µν = 0.
Now, the field equations become [32] h µν = 0 The solutions of eqs. (14) and (15) are plan waves [32] h with Eqs. (14) and (16) represents the equation and the solution for the standard tensor GWs of the GTR [18]. Eqs. (15) and (17) are respectively the equation and the solution for the massive scalar mode instead [32]. We stress that the dispersion law for the modes of the massive scalar field h R is not linear [32]. In fact, the velocity of the tensor modesh µν is the light speed c, but the dispersion law (the second of eq. (18)) for the modes of h R is that of a massive field which is interpreted in terms of a wave-packet [32]. We recall that the group-velocity of a wave-packet of This is exactly the velocity of a massive particle with mass m and momentum − → p . From the second of eqs. (18) and eq. (19) one gets As one wants a constant speed of the wave-packet, one obtains [32].
Let us continue our analysis in the Lorenz gauge [38] with trasformations of the type ǫ ν = 0; these trasformations permit to obtain a condition of transversality for the tensor part of the field: k µ A µν = 0 [32]. On the other hand, they do not give the transversality for the total field h µν . From eq. (11) one gets [32] h µν =h µν −h 2 At this point, if being in the massless case, one could set [39] ǫ µ = 0 Eqs. (23) give the total transversality of the field. On the other hand, in the massive case this is impossible [32]. In fact, if one applies the Dalembertian operator to the second of eqs. (23) and uses the field equations (14) and (15), one gets [32] which is in contrast with the first of eqs. (23). In the same way, it is possible to show that there is no linear relation between the tensorial fieldh µν and the massive scalar field h R [32]. Thus, one cannot choose a gauge in which h µν is purely spatial (that is, one cannot set h µ0 = 0, see eq. (22)) [32]. One can set the traceless condition to the fieldh µν instead [32] ǫ µ = 0 From eqs. (25) one gets [32] ∂ µh µν = 0.
If one wants to preserve the conditions ∂ µh µν andh = 0 transformations like [32] ǫ µ = 0 can be used. Thus, by taking − → p in the z direction, one chooses a gauge in which only A 11 , A 22 , and A 12 = A 21 are different to zero [32]. Settingh = 0 one gets A 11 = −A 22 . Now, one puts these equations in eq. (22), obtaining The term A + (t − z)e (+) µν describes the two standard tensor GW polarizations which arise from the GTR [32]. The term h R (t − v G z)η µν is the massive field arising from the R 2 -gravity theory instead [32]. In other words, the Ricci scalar generates a third massive GW polarization which is not present in the standard GTR [32]. Now, the post-newtonian expansion of the theory requires to assume the space-time metric as a small perturbation expanded around the flat background metric. In this context, after a lengthy algebra one finds the energy-momentum pseudo-tensor of the gravitational field as [19] where h µν now denotes the fluctuating part of the space-time metric, k µ the 4-vector tangent to the world line of a GW andḣ αβ ≡ ∂ 0 h αβ . The rate of energy loss of a matter system coupled to the gravity is found to be [9,18] The symbol e i stands for the unite vector along the i−th axis and the quadrupole moment of mass is defined to be [9] Q ij = mx i x j − r 2 δ ij . Some efforts are devoted to determine the validity of above formula by probing the observational parameters of the binary pulsar PSR 1913+16 [19,21,22]. Setting a 1 = 4 3 in the first term of above equation [18], re-produces the well-known energy loss of the GTR [18]. The gravitational energy radiated due to the Coulomb collision between an electron with charge e and speed v and an ion with charge +Ze, in small-angle scattering regime can be obtained as [27] where B ≡ a2 30 .... Q ij (t) .... Q ij (t) and b denotes the impact parameter. In small-angle approximation, one considers the particle's trajectory as a straight line [25,26]. Let us compute B, that is the contribution of the R 2 term to the gravitational energy loss. For the time derivatives of the quadruple moment one gets ...
By the help of (35), for the time derivatives of the quadrupole moment one obtains Therefore, one gets the contribution of the R 2 term to the gravitational energy loss as The set of following integrals have been used in evaluating (48) Thus, eq. (31) becomes

Thermal Gravitational Radiation of a Hot Plasma
To obtain the gravitational luminosity of a plasma with the gravitational bremsstrahlung as a mechanism for the loss of its energy, one must multiply (31) with the electron flux vn e , ion density n i , and integrate over the impact parameter b [25,26]. Therefore, we obtain the luminosity L (energy loss per volume V ) of the plasma as This integral diverges as b → 0. Thus a cut-off, denoted by b min , is introduced to get a finite result for the luminosity. Based on either classical or quantum mechanical considerations the cut-off takes the form [25,26] respectively. The final result for the luminosity depends on which form for the cut-off is engaged. We will restrict ourself to the Hydrogen plasma. Thus, Z = 1 and n e = n i . Hence from (50), with quantum mechanical cut-off b min = mev , the energy loss takes the form L = 2π 2 e 4 n 2 e G c 5 1 24 One notes that the speed of light, c is restored in (52). By taking thermal average of the above expression, one gets the thermal luminosity of the plasma. In many astrophysical objects, the ratio of Coulomb interaction energy to thermal energy is negligible, so the hot plasma behaves like an almost ideal gas [26]. Thus, one can calculate the thermal luminosity of the plasma by averaging the electron speed in (52) over a thermal distribution of speeds. For an ensemble of particles at temperature T , obeying the Maxwell-Boltzman statistic, the thermal average is where f (v) is an arbitrary function of particle's velocity and k B denotes the Boltzman constant. In particular one obtains v 2n+1 = 2 π n+1 2π βm Thus, from (52) and (54) which is the well-known result derived earlier by Weinberg within the context of the GTR [23] 4 Gravitational Luminosity of the Sun Based on the massive scalar mode arising from the R 2 term, the coupling constant a 2 comes to be very small with respect to linear term R. Assuming the typical galactic scale for the curvature energy, E ≃ 10 45 g, we find a 2 = 10 −34 cm 4 in natural units [32]. In this way, the variation from the standard GTR is very weak. The parameters needed to obtain the above result (in CGS units) are As we can see the first term in Eq. (68) that coming from the standard GTR is in good agreement with available results [23,24,40]. Beside, due to the extra massive mode of gravitational wave arising from the R 2 term in the action, the associated coupling constant of the theory becomes very small and makes the contribution, Eq. (69), more sensible.

Conclusion remarks
The presence of the massive term in Eq. (28) is a characteristic of higher-order terms in f (R)-gravity. Thus, the R 2 -gravity theory include massive GW modes. Hence, our results beside confirming the standard GTR, stimulate the validity of f (R)-gravity. Until now there is not available data to confront our result to the experiment and fix the parameter of the R 2 -gravity contribution.