The energy spectrum of graviton emitted by the black hole binary is calculated in the first part of the chapter. Then, the total quantum loss of energy is calculated in the Schwinger theory of gravity. In the next part, we determine the electromagnetic shift of energy levels of H-atom electrons by calculating an electron coupling to the black hole thermal bath. The energy shift of electrons in H-atom is determined in the framework of nonrelativistic quantum mechanics. In the last section, we determine the velocity of sound in the black hole atmosphere, which is here considered as the black hole photon sea. Derivation is based on the thermodynamic theory of the black hole photon gas.
- Schwinger source theory
- spectrum of H-atom
- Coulomb potential
- black hole spectrum
- energy shift
1. The graviton spectrum of the black hole binary
In 1916, Schwarzschild published the solution of the Einstein field equations  that were later understood to describe a black hole [2, 3], and in 1963, Kerr generalized the solution to rotating black holes . The year 1970 was the starting point of the theoretical work leading to the understanding of black hole quasinormal modes [5, 6, 7], and in the 1990s, higher-order post-Newtonian calculations  were performed and later the extensive analytical studies of relativistic two-body dynamics were realized [9, 10]. These advances, together with numerical relativity breaks through in the past decade [11, 12, 13]. Numerous black hole candidates have now been identified through electromagnetic observations [14, 15, 16]. The black hole binary and their rotation and mergers are open problem of the astrophysics, and it is the integral part of the binary black hole physics.
The binary pulsar system PSR B1913+16 (also known as PSR J1915+1606) discovered by Hulse and Taylor  and subsequent observations of its energy loss by Taylor and Weisberg  demonstrated the existence of gravitational waves .
By the early 2000s, a set of initial detectors was completed, including TAMA 300 in Japan, GEO600 in Germany, the Laser Interferometer Gravitational-Wave Observatory (LIGO) in the United States, and Virgo in Italy. In 2015, Advanced LIGO became the first of a significantly more sensitive network of advanced detectors (a second-generation interferometric gravitational wave detector) to begin observations .
Taylor and Hulse, working at the Arecibo Radiotelescope, discovered the radio pulsar PSR B1913+16 in a binary, in 1974, and this is now considered as the best general relativistic laboratory .
Pulsar PSR B1913+16 is the massive body of the binary system where each of the rotating pairs is 1.4 times the mass of the Sun. These neutron stars rotate around each other in an orbit not much larger than the Sun’s diameter, with a period of 7.8 h. Every 59 ms, the pulsar emits a short signal that is so clear that the arrival time of a 5 min string of a set of such signals can be resolved within 15 s.
A pulsar model based on strongly magnetized, rapidly spinning neutron stars was soon established as consistent with most of the known facts ; its electrodynamical properties were studied theoretically  and shown to be plausibly capable of generating broadband radio noise detectable over interstellar distances. The binary pulsar PSR B1913+16 is now recognized as the harbinger of a new class of unusually short-period pulsars, with numerous important applications.
Because the velocities and gravitational energies in a high-mass binary pulsar system can be significantly relativistic, strong-field and radiative effects come into play. The binary pulsar PSR B1913+16 provides significant tests of gravitation beyond the weak-field, slow-motion limit [24, 25].
We do not repeat here the derivation of the Einstein quadrupole formula in the Schwinger gravity theory . We show that just in the framework of the Schwinger gravity theory, it is easy to determine the spectral formula for emitted gravitons and the quantum energy loss formula of the binary system. The energy loss formula is general, including black hole binary, and it involves arbitrarily strong gravity.
Since the measurement of the motion of the black hole binaries goes on, we hope that sooner or later the confirmation of our formula will be established.
1.1 The Schwinger approach for the problem
Source methods by Schwinger are adequate for the solution of the calculation of the spectral formula of gravitons and energy loss of binary. Source theory [27, 28] was initially constructed to describe the particle physics situations occurring in high-energy physics experiments. However, it was found that the original formulation simplifies the calculations in the electrodynamics and gravity, where the interactions are mediated by photon and graviton, respectively. The source theory of gravity forms the analogue of quantum electrodynamics because, while in QED the interaction is mediated by the photon, the gravitational interaction is mediated by the graviton . The basic formula in the source theory is the vacuum-to-vacuum amplitude :
where the minus and plus symbols refer to any time before and after the region of space–time with action of sources. The exponential form is postulated to express the physically independent experimental arrangements, with result that the associated probability amplitudes multiply and the corresponding expressions add [27, 28].
where the dimensionality of has the dimension of the Planck constant and is the momentum and energy tensor that, for a particle trajectory , is defined by the equation 
where is the relativistic four-momentum of a particle with a rest mass and
and the relativistic energy is defined by the known relation
where v is the three-velocity of the moving particle.
Symbol in Eq. (2) is defined as , and is the graviton propagator whose explicit form will be determined later.
1.2 The power spectral formula in general
It may be easy to show that the probability of the persistence of vacuum is given by the following formula :
where the so-called power spectral function has been introduced . For the extraction of the spectral function from Im , it is necessary to know the explicit form of the graviton propagator . This propagator involves the graviton property of spreading with velocity . It means that its mathematical form is identical with the photon propagator form. With regard to Schwinger et al. , the -representation of in Eq. (2) is as follows:
1.3 The power spectral formula for the binary system
In the case of the binary system with masses and , we suppose that they move in a uniform circular motion around their centre of gravity in the plane, with corresponding kinematical coordinates:
For the tensor of energy and momentum of the binary, we have
where we have omitted the tensor , which is associated with the massless, gravitational field distributed all over space and proportional to the gravitational constant .
1.4 The quantum energy loss of the binary
Using the following relations
we get for , with being restricted to positive:
Using the definition of the Bessel function
from which the derivatives and their integrals follow, we get for and the following formulae:
Using , where is supposed to be small in comparison with radii and , we obtain
So, instead of Eq. (19), we get
From  we have Kapteyn’s formula:
After differentiating the last relation with respect to , we have
From  we learn other Kapteyn’s formulae:
Instead of using Kapteyn’s formulae for the interference term, we will perform a direct evaluation of the energy loss of the interference term by the -integration in (27) . So, after some elementary modification in the -integral, we get
Using the definition of the -function and its derivative, we have, instead of Eq. (34), with
According to the Schwinger article , we express the delta function as follows:
and it means that
Now, we can write Eq. (36) in the following form after some elementary operations:
where can be written as follows:
After the application of the per partes method, we get from Eq. (40) the following mathematical object:
We get after some elementary operations
where are functions which must be determined.
So we get instead of Eq. (41) the following final form:
Let us remark that we can use simple approximation in Eq. (41) as follows: . Then, after using the well-known formula
So, instead of Eq. (46), we have
2. Energy shift of H-atom electrons due to the black hole thermal bath
We here determine the electromagnetic shift of energy levels of H-atom electrons by calculating an electron coupling to the black hole thermal bath. The energy shift of electrons in H-atom is determined in the framework of nonrelativistic quantum mechanics.
The Gibbons-Hawking effect is the statement that a temperature can be associated to each solution of the Einstein field equations that contain a causal horizon.
Schwarzschild space-time involves an event horizon associated with temperature of a black hole of mass . We consider here the influence of the heat bath of the Gibbons-Hawking photons on the energy shift of H-atom electrons.
The analogical problems are solved in the scientific respected journals. There is a general conviction of an analogy between the black hole and the hydrogen atom. Corda  used the model where Hawking radiation is a tunneling process. In his article the emission is expressed in terms of the black hole quantum levels. So, the Hawking radiation and black hole quasinormal modes by Corda  are analogical to hydrogen atom by Bohr.
In this model  the corresponding wave function is written in terms of a unitary evolution matrix. So, the final state is a pure quantum state with no information loss. Black hole is defined as the quantum systems, with discrete quantum spectra, with Hooft’s assumption that Schrödinger equations are universal for all universe dynamics.
Thermal photons by Gibbons and Hawking are blackbody photons, with the Planck photon distribution law [40, 41, 42], derived from the statistics of the oscillators inside of the blackbody. Later Einstein  derived the Planck formula from the Bohr model of atom where photons and electrons have the discrete energies related with the Bohr formula , being the initial and final energies of electrons.
Now, we determine the modification of the Coulomb potential due to blackbody photons. At the start, the energy shift in the H-atom is the potential , generated by nucleus of the H-atom. The potential at point is [44, 45]
The average of the last equation in space enables the elimination of the so-called the effective potential:
where is the average value of the square coordinate shift caused by the thermal photons. The potential shift follows from Eq. (51):
The shift of the energy levels is given by the standard quantum formula :
In case of the Coulomb potential, which is the case of the H-atom, we have
Then for the H-atom we can write
where we used the following equation for the Coulomb potential
The motion of electron in the electric field is evidently described by elementary equation:
which can be transformed by the Fourier transformation into the following equation
where the index concerns the Fourier component of the above functions.
Using Bethe idea  of the influence of vacuum fluctuations on the energy shift of electron, the following elementary relations were applied by Welton , Akhiezer et al.  and Berestetzkii et al. :
and in case of the thermal bath of the blackbody, the last equation is of the following form :
because the Planck law in (60) was written as
where the term
is the average energy of photons in the blackbody and
is the number of electromagnetic modes in the interval .
where involves the number of frequencies in the interval .
So, after some integration, we get
where is the primitive function of the omega-integral with
which is not elementary, and it is not in the tables of integrals.
Frequencies and can be determined from the field of thermal photons. It was performed for the Lamb shift [44, 47] caused by the interaction of the Coulombic atom with the field fluctuations. The Bethe-Welton method is valid here too and so we take Bethe-Welton frequencies. It means an electron does not respond to the fluctuating field if the frequency is much less than the atom binding energy given by the Rydberg constant  . So, the lower frequency limit is
where is so-called the fine structure constant.
The second frequency follows from the cutoff, determined by the neglection of the relativistic effect in our theory. So, we write
If we express the thermal function in the form of the geometric series
and the first thermal contribution is
then, with Eq. (55)
where according to Sokolov et al. 
Let us only remark that the numerical form of Eq. (72) has deep experimental astrophysical meaning.
Haroche  and his group performed experiments with the Rydberg atoms in a cavity. We used here Gibbons-Hawking black hole for the determination of the energy shift of H-atom electrons in the black hole gas.
3. Velocity of sound in the black hole photon gas
We have seen that the black hole can be modeled by the blackbody, and it means that there is the velocity of sound in the Gibbons-Hawking black hole thermal bath. So, let us derive the sound velocity from the thermodynamics of photon gas and energy mass relation.
In order to be pedagogically clear, we start with the derivation of the speed of sound in the real elastic rod.
Let be the cross-section of the element of a rod on the axis . Let be the deflection of at point at time . The shift of the at point is evidently
Now, we suppose that the force tension acting on the of the rod is given by Hooke’s law:
where is Young’s modulus of elasticity. We easily derive that
The mass of is , where is the mass density of the rod and the dynamical equilibrium is expressed by Newton’s law of force:
is the velocity of sound in the rod.
The complete solution of Eq. (79) includes the initial and boundary conditions. We suppose that Eq. (80) is of the universal validity also for gas in the cylinder tube. If is the relative prolongation of a rod, then an analogue for the tube of gas is , , where is the volume of a gas and is gas pressure. Then, the modulus of elasticity as the analogue of Eq. (76) is
The sound in ideal gas is the adiabatic thermodynamic process with no heat exchange. This is the model of the sound spreading in the gas of blackbody photons. Such process is described by the thermodynamic equation:
where is the Poisson constant defined as , with being the specific heat under constant pressure and under constant volume.
After differentiation of Eq. (82), we get the following equation:
with being the gas mass density.
The equilibrium radiation density has the Stefan-Boltzmann form:
Then, with regard to the thermodynamic definition of the specific heat,
Similarly, with regard to the general thermodynamic theory,
because for photon gas, and in such a way, for photon gas. According to the theory of relativity, there is a relation for mass and energy, namely, . At the same time, the pressure and the internal energy of the blackbody gas are related as . So, in our case
which was derived by Partovi  using the QED theory of the photon gas. We correctly derived .
So, we have performed the derivation of the velocity of sound in the relic photon sea. It is not excluded that the relic sound can be detected by the special microphones of Bell Laboratories. If we use van der Waals equation of state or the Kamerlingh Onnes virial equation, the obtained results will be modified with regard to the basic results.
Our derivation of the light velocity in the blackbody photon gas was based on the classical thermodynamic model with the adiabatic process (), controlling the spreading of sound in the gas. Partovi  derived additional radiation corrections to the Planck distribution formula and the additional correction to the speed of sound in the relic photon sea. His formula is of the form
where is the fine structure constant and G Kelvin. We see that our formula is the first approximation in the Partovi expression.
There is the Boltzmann statistical theory of transport of sound energy in a gas . After the application of this theory to the photon gas or relic photon gas, we can obtain results involving the cross-section of the photon-photon interaction :
where cm is the classical radius of electron and is the fine structure constant with numerical value .
4. Discussion and summary
We have derived the spectral density of gravitons and the total quantum loss of energy of the black hole binary. The energy loss is caused by the emission of gravitons during the motion of the two black hole binaries around each other under their gravitational interaction. The energy loss formulae of the production of gravitons are derived here by the Schwinger method. Because the general relativity and theory of gravity do not necessarily contain the last valid words to be written about the nature of gravity and it is not, of course, a quantum theory , they cannot give the answer on the production of gravitons and the quantum energy loss, respectively. So, this article is the original text that discusses the quantum energy loss caused by the production of gravitons by the black hole binary system. It is evident that the production of gravitons by the binary system forms a specific physical situation, where a general relativity can be seriously confronted with the source theory of gravity.
This article is an extended version of an older article by the present author , in which only the spectral formulae were derived. Here we have derived the quantum energy loss formulae, with no specific assumption concerning the strength of the gravitational field. We hope that future astrophysical observations will confirm the quantum version of the energy loss of the binary black hole.
In the next part of the chapter, the electromagnetic shift of energy levels of H-atom electrons was determined by calculating an electron coupling to the Gibbons-Hawking electromagnetic field thermal bath of the black hole. The energy shift of electrons in H-atom is determined in the framework of nonrelativistic quantum mechanics.
In the last section, we have determined the velocity of sound in the blackbody gas of photons inside of the black hole. Derivation was based on the thermodynamic theory of the photon gas and the Einstein relation between energy and mass. The spectral form for the n-dimensional blackbody was not here considered. The text is based mainly on the author articles published in the international journals of physics [33, 54, 55].
There is the fundamental problem concerning the maximal mass of the black hole. The theory of the space–time with maximal acceleration constant was derived by authors [56, 57]. In this theory the maximal acceleration constant is the analogue of the maximal velocity in special theory of relativity. Maximal acceleration determines the maximal black hole mass where the mass of the black hole is restricted by maximal acceleration of a body falling in the gravity field of the black hole.
Another question is what is the relation of our formulae to the results obtained by LIGO (Laser Interferometer Gravitational-Wave Observatory)? LIGO is the largest and most sensitive interferometer facility ever built. It has been periodically upgraded to increase its sensitivity. The most recent upgrade, Advanced LIGO (2015), detected for the first time the gravitational wave, with sensitivity far above the background noise. The event with number GW150914 was identified with the result of a merger of two black holes at a distance of about 400 Mpc from Earth . Two additional significant detections, GW151226 and GW170104, were reported later. We can say that at this time it is not clear if the LIGO results involve information on the spectrum of gravitons calculated in this chapter.