Coalescence of Kerr Black Holes—Binary Systems from GW150914 to GW170814

We investigate the energy of the gravitational wave from a binary black hole merger by the coalescence of two Kerr black holes with an orbital angular momentum. The coalescence is constructed to be consistent with particle absorption in the limit in which the primary black hole is sufficiently large compared with the secondary black hole. In this limit, we analytically obtain an effective gravitational spin–orbit interaction dependent on the alignments of the angular momenta. Then, binary systems with various parameters including equal masses are numerically analyzed. According to the numerical analysis, the energy of the gravitational wave still depends on the effective interactions, as expected from the analytical form. In particular, we ensure that the final black hole obtains a large portion of its spin angular momentum from the orbital angular momentum of the initial binary black hole. To estimate the angular momentum released by the gravitational wave in the actual binary black hole, we apply our results to observations at the Laser Interferometer Gravitational-Wave Observatory: GW150914, GW151226, GW170104, GW170608 and GW170814.


Introduction
The coalescence of black holes is one of the most important sources of gravitational waves. A gravitational wave occurs owing to a variation in the gravitational field, such as the motion of massive bodies. Although a gravitational wave can commonly occur owing to a small variation in the gravitational field, the magnitude of the gravitational wave in such a variation is too small to detect at our observable accuracy. Hence, an observable gravitational wave needs a sufficient magnitude so that it can be detected at observatories such as the Laser Interferometer Gravitational-Wave Observatory (LIGO). Since the magnitude of a gravitational wave becomes large along with the mass of its source, the candidate sources of detectable gravitational waves should be massive. Thus, the coalescence of black holes can be a source that releases a gravitational wave that is sufficiently large to be detected at the LIGO. Nowadays, several gravitational waves are detected at observatories. Most of their sources are the binary black hole mergers in GW150914, GW151226, GW170104, GW170608, GW170814 [1][2][3][4][5]. Thus, the coalescence of black holes can be a frequently detectable source of gravitational waves. Further, a system of binary black holes has three angular momenta-two spin angular momenta of the black holes and one orbital angular momentum-and these angular momenta play important roles in the analysis of the system. Note that the angular momenta for a binary black hole merger are described by the primary spin parameter a 1 , secondary spin parameter a 2 and orbital angular momentum L orb in the initial state, and the spin parameter a f in the final state.
Black holes have conserved quantities such as mass and angular momentum. Two types of energies are included in the mass of the black hole: a reducible energy and an irreducible mass. The roles of these two energies can be clarified in consideration of particle absorption to vary the black hole [6,7]. In particle absorption, the variations in the mass and angular momentum of a black hole can be related to a particle's momenta such as the radial and angular momenta. Here, the black hole's angular momentum can be reduced by our choice of the particle's angular momentum; thus, part of the particle's angular momentum is a reducible energy in the variation of the black hole's mass. Interestingly, the remaining part cannot be reduced with any choice of the particle's momenta; therefore, the remaining part is called the irreducible mass of a black hole [6,8]. This irreducible mass can be understood as the surface energy of the black hole [9]. For example, even if the mass of a black hole decreases during particle absorption, the irreducible mass still increases [8]. Since this behavior of the irreducible mass is very similar to the entropy in thermodynamics, its square-the area of the horizon-is related to the Bekenstein-Hawking entropy [10,11]. Further, the black hole can emit energy away from its horizon according to a quantum process. Owing to this small portion of the energy, the black hole can be assumed to be a thermal object having the Hawking temperature accounted from the emission [12,13]. For the variation in the black hole's mass, the Bekenstein-Hawking entropy and Hawking temperature act as thermodynamic variables in the first law of thermodynamics. In addition, the Bekenstein-Hawking entropy increases in an irreversible process in the second law of thermodynamics according to the behavior of the irreducible mass. Thus, the black hole itself can be considered as a thermal system obeying the laws of thermodynamics.
The black hole has a surface called the event horizon through which a particle or matter cannot escape from the gravity of the black hole. Then, the particle cannot be observed outside the horizon anymore. While absorbing a particle, the black hole also undergoes variations with respect to its physical properties. The stability of the black hole's horizon is one of its properties that are testable using particle absorption with respect to the black hole. The horizon of the black hole should be stable to cover the singularity located within it because a naked singularity, which is an observable singularity without a horizon, causes a breakdown in the causality of spacetime. This is suggested as a weak cosmic censorship conjecture that prevents a naked singularity [14,15]. Various investigations of the weak cosmic censorship conjecture have been applied to various black holes. The conjecture for a Kerr black hole was first tested by adding a particle [16]. Here, the conserved quantities carried by the particle cause a variation in the Kerr black hole with respect to its corresponding conserved quantities. By the addition of the particle's angular momentum, the angular momentum of the Kerr black hole can increase but it cannot exceed the extremal limit of the black hole; thus, its horizon still exists and covers the singularity inside it. This implies that the weak cosmic censorship conjecture is valid when adding a particle. However, the validity of the conjecture depends on the test approach. For example, in a near-extremal Kerr black hole, the conjecture was found to be invalid [17]; thus, its angular momentum can be overspun owing to the addition of a particle. Then, the horizon disappears. This can be resolved by considering the self-force effect with the addition of a particle [18][19][20]. The conjecture for a Reissner-Nordström (RN) black hole is valid when considering the back-reaction effect [21,22]. There are also various tests of the conjecture for anti-de Sitter (AdS) black holes . In particular, from a thermodynamic point of view, since the second law of thermodynamics ensures an increase in the area of the black hole's horizon, it can be a sufficient condition for the validity of the conjecture. This thermodynamic picture has been investigated for AdS black holes for particle absorption [49][50][51].
The coalescence of black holes is also understood from a thermodynamic point of view. As a thermal system, the coalescence can be an irreversible process in which the initial two black holes become a final black hole given a thermal preference. The increase in the entropy due to the irreversible process, the second law of thermodynamics, was first shown for the coalescence of two Schwarzschild black holes [52]. Further, during coalescence, the energy released by the gravitational wave is expected to be in terms of the upper limit of the energy of the gravitational wave under the second law of thermodynamics. However, the energy limit is much larger than the real energy of the gravitational wave; therefore, there are some difficulties in deriving detailed physical implications from the limit. In spite of these difficulties, when applied to the coalescence of Kerr black holes, the upper limit indicates the existence of an effective gravitational spin interaction between black holes [53], where attraction or repulsion acts on the black holes according to their alignments. For the case where one of the black holes is much smaller than the other black hole, the form of the interaction is exactly coincident with the interaction potential acting on a particle spinning around the Kerr black hole, as already obtained from the Mathisson-Papapetrou-Dixon (MPD) equations [54][55][56][57][58][59][60][61][62][63]. Hence, a gravitational spin interaction is induced from the coupling of the spin angular momenta of the two black holes. When the spin angular momentum of a black hole is sufficiently large, as in Myers-Perry (MP) black holes where there is no extremal limit for the spin angular momentum in higher-dimensional spacetime, the interaction between MP black holes plays an important role in the instability of the black hole system and the upper limit of the energy released by the gravitational wave [64]. Further, using the upper limit, the various constants of a given model of a regular black hole can be specified with physically possible ranges by matching LIGO observations [65][66][67].
In this paper, we will investigate the gravitational wave released from a binary black hole merger by the coalescence of two Kerr black holes with an orbital angular momentum. Since an astrophysical black hole can be theoretically approximated as a Kerr black hole [1,68], we consider a binary system consisting of two spin angular momenta due to Kerr black holes and one orbital angular momentum for their orbit. In particular, the orbital angular momentum is found to be important to the state of the final black hole because its angular momentum is provided by the orbital angular momentum. During coalescence, we assume conservation of the system's energy and angular momentum. In addition, the irreducible mass of Kerr black holes is assumed to be irreducible. Because the irreducible mass cannot be extracted by a physical process, including a Penrose process [7,69], it can be expected to not decrease during the coalescence of black holes, which is an irreversible process. This implies that the irreducible mass is used for the formation of the final black holes and a gravitational wave is released from the reducible mass such as the kinetic and rotational energies in the initial state. Our assumptions for coalescence are consistent with the second law of thermodynamics. In the particle limit for one of the black holes, our model is proven to be consistent with particle absorption [6,8] and the MPD equations [53]. Further, we analytically obtain that the energy released by the gravitational wave depends on the effective gravitational spin-spin and spin-orbit interactions. In particular, the spin-orbit interaction is newly obtained from an effective force and is one of the advantages of our model, which considers the orbital angular momentum. Then, we numerically compute the final black hole and the energy released by the gravitational wave to obtain the effects of various variables during coalescence for the equal-mass case. Interestingly, the obtained energy of the gravitational wave is within a very similar range as those of LIGO observations [1][2][3][4][5]; therefore, we can ensure that our results are relevant to the understanding of an actual binary black hole merger. This is another advantage of our model. By numerical computation, the relations between the alignments of the black holes and the gravitational wave are found to be consistent with the effective interactions. In addition, the orbital angular momentum also shows similar effects as the spin angular momentum. Finally, we apply our results to five recent LIGO observations to find the relevant range of angular momenta in our model. Here, we obtain the initial spin parameters, orbital angular momentum, total angular momentum and angular momentum of the gravitational wave. Their median values are set to be very consistent with the conservation of angular momentum. In particular, the orbital angular momentum constitutes the largest portion of the total angular momentum in the initial state. Hence, it can be important to the formation of a black hole with a high spin parameter during coalescence. This paper is organized as follows. In Section 2, we review Kerr black holes. Particle absorption and the MPD equations are also introduced to show the analytical forms of the effective interactions. In Section 3, we construct the basic framework of our model. Then, to prove its consistency with particle absorption and the MPD equations, we compute a variation of our model in the particle limit. In this process, we will prove that the effective interactions predicted by our model are coincident with those of particle absorption and the MPD equations. In Section 4, we numerically obtain the final states in equal-mass cases. Here, the effective interactions are obtained as analytical forms and found to work as predicted in the particle limit. In Section 5, we apply our model to recent LIGO observations of binary black hole mergers. We compute the initial spin parameters by using other parameters. Further, the orbital angular momentum and the angular momentum of the released gravitational wave are obtained by our model. In Section 6, we briefly summarize our results.

Thermodynamics and Spin Interaction of Kerr Black Hole
We will assume that the binary black hole is a system consisting of two Kerr black holes with an orbital angular momentum in their orbit. The angular momenta of the binary black hole play important roles in the effective interactions between black holes. Further, during the coalescence of the binary black hole, the energy released by the gravitational wave is significantly affected by the interactions. In this work, we will construct the coalescence of the binary black hole that satisfies the laws of thermodynamics. In this procedure, the effective interactions related to the angular momenta of the system will be shown to be consistent with what are expected from the particle absorption and MPD equations at the particle limit of one of the black holes. Hence, we will review the basics of Kerr black holes, particle absorption and the gravitational spin interaction from the MPD equations.

Basics of Kerr Black Hole
The Kerr black hole is a solution to the Einstein equations in four-dimensional spacetime. The metric of a Kerr black hole with a mass M and spin angular momentum J in Boyer-Lindquist coordinates is where the spin parameter is given as a, which is bounded when M ≥ a. For M > a, the Kerr black hole has two event horizons.
where r in and r h denote the inner and outer horizons, respectively. Since the inside of the outer horizon cannot be seen by an asymptotic observer, the properties of the Kerr black hole are defined at its outer horizon. An asymptotic observer measures the angular velocity of spacetime at the outer horizon as The Hawking temperature and Bekenstein-Hawking entropy are For a given mass, the spin parameter a is saturated at M = a, where the Kerr black hole satisfies the extremal condition. Under the extremal condition, the inner and outer horizons are coincident with each other. Then, When the spin parameter exceeds the value of the mass, M < a, there is no horizon covering the inside of the black hole. Then, the curvature singularity of spacetime is exposed to the observer located outside. This is called a naked singularity. However, according to the cosmic censorship conjecture, it is expected that there is no physical process for overspinning a Kerr black hole into a naked singularity [14,15]. This can be shown by particle absorption, among the various verifications investigated for the cosmic censorship conjecture for Kerr black holes [16].

Thermodynamics of Kerr Black Hole under Particle Absorption
In particle absorption, the acceleration of a Kerr black hole's angular velocity is considered by adding a particle. Owing to the energy and angular momentum of the particle, the corresponding conserved quantities of the Kerr black hole vary in the energy equation of the particle. In these variations, the mass of the Kerr black hole is divided into two parts: the reducible energy and the irreducible mass [6,8]. In a Kerr black hole, the reducible energy includes the rotational and kinetic energies. Further, the irreducible mass is a type of energy distributed on the surface of the horizon [9]. We will introduce the irreducible mass by using particle absorption. To obtain a relation between the conserved quantities of a particle passing through the outer horizon of the Kerr black hole, the particle's equations of motion are obtained using Hamilton-Jacobi method.
which are the Hamiltonian and Hamilton-Jacobi action of a particle, respectively, in the metric of a Kerr black hole. Using a separate variable K [70], we can write the radial and θ-direction equations of motion as In combination with Equation (7), the removal of K gives the particle's energy equations for a given location. Then, where the particle is assumed to be absorbed into the Kerr black hole when it passes through its outer horizon. At that moment, the relation between the energy and momenta of the particle is given by the energy equation, which becomes [6,8] where we choose the (+) sign among the two solutions because the particle enters the Kerr black hole in a positive time flow. Then, the solution shows the relation between the particle's energy and momenta at the outer horizon. When the particle passes through the outer horizon, the energy and angular momentum of the particle are assumed to be those of the Kerr black hole. Then, the conserved quantities of the black hole vary as much as those of the particle; thus, For particle absorption, we can write the relation between the variations in the mass and angular momentum of the Kerr black hole using Equation (9). Then, which constrains the variation in the Kerr black hole in the particle absorption process [16]. Under this constraint, the variation in the entropy of the Kerr black hole becomes which is always positive because r h > M and the equality r h = M also results in an increase in the entropy. This is consistent with the second law of thermodynamics [49]. We can obtain the first law of thermodynamics by inserting Equation (12) into Equation (11) [51]. Then, Therefore, particle absorption varies the Kerr black hole and satisfies the laws of thermodynamics. Here, we can obtain an interesting property of black holes for particle absorption. The rewritten Equation (11) becomes an inequality given as where the left-hand side is an irreducible property in the process. By integrating out the left-hand side of Equation (34), we can define a property having the same dimension as the mass [6,8].
which is called the irreducible mass [6,8,10]. The irreducible mass is assumed to be an energy distributed on the surface of the horizon [9]. Then, in terms of the irreducible mass, the mass of the Kerr black hole can be divided into the irreducible mass and a rotation energy such that [71] Hence, the mass of the Kerr black hole actually consists of irreducible and reducible masses. Then, the irreducible mass increases in a physical process such as the Penrose process [69], even if the mass of the Kerr black hole can be reduced by the extraction of the rotational energy in Equation (16).

Gravitational Spin-Spin Interaction in MPD Equations
We will consider a binary black hole system; thus, the binary system can be expected to have an effective gravitational spin interaction between two black holes with their spin angular momenta. Using the second law of thermodynamics, the increase in the entropy, the contribution of the spin-spin interaction is estimated in the energy of the gravitational wave during the coalescence of the black holes and the form of the spin-spin interaction potential derived from the MPD equations for a spinning particle is clearly coincident in the limit where one of the black holes is much smaller and slowly rotating compared with the other black hole [53]. In this section, we review these results with an introduction to the MPD equations [54][55][56] in the Kerr black hole spacetime. We consider a Kerr black hole with a mass M 1 and an angular momentum J 1 = M 1 a 1 and a spinning particle with a mass M 2 and an angular momentum J 2 . The spinning particle also has a four-velocity v µ and linear momentum p a . In the MPD equations, the momentum and four-velocity have a difference related to a proper time s due to the spinning effect. Then, where R a bcd is the Riemann curvature tensor of the Kerr metric. The spin of the particle is given in terms of the spin tensor S ab and spin vector S a . The trajectory of the spinning particle can be derived by imposing the supplementary condition [72] p a S ab = 0.
Then, we can determine the motion of the spinning particle using Equations (17) and (18). The mass and spin angular momentum of the spinning particle are defined as To obtain the form of the spin interaction potential, we assume that the axes of the spin angular momenta of both the black hole and particle are coincident with each other and the pole of the particle is aligned and moves towards that of the black hole. This way, we can prevent the precession of the spinning particle. Further, the particle slowly enters the black hole, such that the particle is assumed to be nonrelativistic, v a 1. Here, the nonrelativistic limit is possible because we can freely choose the initial condition of the particle. Then, the initial state for the velocity v a and spin vector S a of the particle is given as , v a , 0, 0 , J a 2 = 0, where the two vectors are normalized with the nonrelativistic limit of the particle. The interaction potential between the black hole and the particle is obtained from the energy of the spinning particle, which is a conserved quantity with respect to the Killing vector ξ t in terms of the time coordinate. Hence, the energy of the spinning particle is derived as where the first term is related to the kinetic energy and the second term is the spin interaction with respect to the spin tensor. Then, the second term in Equation (21) becomes which is exactly coincident with the spin interaction derived from the second law of thermodynamics [53]. The sign of the spin potential in Equation (22) depends on the alignment between J 1 and J 2 and implies effective attraction and repulsion between the black hole and the particle. The parallel alignment, which has a positive sign, has a positive potential; thus, the effective force acts in a repulsive manner. The antiparallel alignment has a negative sign; hence, the effective force acts as an attraction. This effective interaction plays an important role in the gravitational wave released in the collision of the black holes. We will derive using our approach for the particle limit and investigate the effect of the potential in the following sections. Note that we will use dimensionless coordinates and variables scaled by the solar mass M , such as where we will omit the tildes for convenience.

Basic Framework
The effects of the angular momenta in a binary black hole merger as a source of the gravitational wave detected at the LIGO will be investigated. Here, we assume an initial binary black hole in a model that consists of two Kerr black holes with an orbital angular momentum. Then, the binary black hole merger produces a final Kerr black hole with the released gravitational wave. In this model, the energy of the gravitational wave can be estimated, satisfying the laws of thermodynamics between the initial and final states. In the initial state, two Kerr black holes are located far from each other; thus, their gravitational interaction can be ignored. These Kerr black holes rotate with the orbital angular momentum L orb , which will be included in the total angular momentum. The primary and secondary black holes are (M 1 , a 1 ) and (M 2 , a 2 ) in the initial state and their axes of spin angular momenta have an angular difference ψ. These angular momenta will play an important role in explaining the final state of the black hole system. According to Equation (16), the energies of the initial state are divided into an irreducible mass and a rotation energy including the orbital angular momentum. Then, as they slowly come together with a spiral motion due to the orbital angular momentum, the two Kerr black holes merge into a Kerr black hole of (M f , a f ) in the final state. Since the total energy of the system should be conserved in the coalescence, the released gravitational wave is equivalent to the loss of mass between the initial and final states. Then, the energy and angular momentum of the gravitational wave are where, for simplicity, we assume in the initial state that the sum of J 1 and J 2 is aligned with L orb ; thus, J f and J gw are also aligned with the same rotating plane. Then, the magnitude of the sum of angular momenta in the initial state becomes where the sum of J 1 and J 2 and the final angular momentum J f are aligned with L orb according to our assumption. The ratios of the mass and angular momentum of the gravitational wave with respect to the total mass and angular momentum in the initial state are respectively defined as To estimate the value of M gw , we assume that the total irreducible mass increases during coalescence. Since the irreducible mass cannot be extracted or decreased by the Penrose process in the Kerr black hole, our assumption can be a reasonable generalization of particle absorption for the coalescence of black holes. This is our main assumption in this work. Fortunately, from a thermodynamic point of view, this will provide quite precise predictions about the energy of the gravitational wave released in the coalescence compared with its upper limits in References [52,53,64,67]. Then, the increase in the irreducible mass from the initial state to the final state is Our assumption in Equation (27) satisfies the second law of thermodynamics. According to the definition of the Bekenstein-Hawking entropy in Equation (4), the increase in the irreducible mass in Equation (27) becomes a sufficient condition; hence, Therefore, the entropy of the system increases during the process of coalescence. Then, our assumption is relevant to the second law of thermodynamics. Physically, the increase in the irreducible mass implies that a Kerr black hole is approximated as a solid body with a spin angular momentum. Hence, most of the energy of the gravitational wave is released from the reducible energy M re included in the mass M. Here, the reducible energy is simply assumed as the difference between the mass and irreducible mass of the Kerr black hole. Owing to our assumption for the initial condition, the main part of the reducible energy is the rotational energy; therefore, we define the reducible energy by the rotational energy: M rot = M − M ir . In combination with Equations (24), (26) and (27), the upper limit of the energy of the gravitational wave can be estimated as The upper limit of the energy of the gravitational wave precisely approaches the real value of the gravitational wave in LIGO observations with our model. It will be investigated in the following sections. There are two ratios related to the mass and angular momentum. We will find the constraints on M and J considering particle absorption. Further, we will obtain a type of effective interaction related to the orbital angular momentum L orb . Note that the mass of the black hole includes the rotation effect according to Equation (16). Hence, the source of the gravitational wave includes the gravitational energy of the black hole and its rotation effect. We consider herein an asymptotic observer measuring initial and final states of black holes. Then, the observer only measures masses and angular momenta of black holes, because gravitational potential is decayed at the asymptotic region. Here, we assumed no interaction between initial black holes, because black holes are assumed far from each other. For simplicity, we have removed herein the effect of the gravitational potential by setting long distance between black holes. If we consider a binary black hole system in a real case, two black holes interact with each other in a finite distance. Hence, the gravitational interaction energy can be considered to describe the system and the interaction can contribute to the energy of the gravitational wave.

Ratios Related to the Mass and Angular Momentum
Our assumption for the coalescence of Kerr black holes should be consistent with particle absorption when M 1 M 2 , where the primary black hole is fixed as a background and the secondary black hole is treated as a particle. Hence, physical constraints on M and J can be obtained from particle absorption. In particle absorption, the ratios M and J become very small because M 1 M 2 and M gw 1, where we will assume that M 1 = M and M 2 = E. The angular momentum corresponds to the orbital angular momentum; thus, L orb = L and J 1 L because the particle has only an orbital angular momentum. Then, and from Equation (24), By a combination of Equations (30) and (31), Because the primary black hole is fixed as the background and the variables are the particle's energy and momenta during particle absorption [6,8], we have to assume that (1 − M )M ≈ M and (1 − J )J ≈ J to fix the mass and angular momentum of the primary black hole. This implies that there is no radiation from the primary black hole as the background. Then, all of the variations in the black hole originate from the variables of the particle. The coalescence of black holes is not technically coincident with particle absorption; hence, we need to modify this part. Then, the particle contributes to the modified variation in the primary black hole as which implies that part of the particle's energy is absorbed into the black hole. The other part is released by the gravitational wave. This is also consistent with Equation (10) with addition of ratios. Using Equation (33), the relation in Equation (9) becomes from which the variation in the entropy becomes Owing to the second law of thermodynamics, the entropy should increase in particle absorption owing to the irreversible process. The sign of the first term in Equation (35) depends on the alignment between J and L; thus, we assume and consider the specific case that M ≈ J to remove the dependency on the alignment. Then, the entropy becomes irreducible. Hence, under our analysis, the ratio of the emitted angular momentum will be as much as that of the emitted mass between the initial and final states. Then, this ensures the increase in the irreducible mass of the system under M = J as Here, the initial condition related to the slow approach of the black holes during their coalescence can be converted to |p r | ≈ 0 in particle absorption. Therefore, we will assume that Note that our assumptions can be applied limited cases to the gravitational wave but these reduce the number of variables related to ratios to one and simplify our analysis.

Gravitational Spin-Orbit Interaction from Particle Absorption
Since most of the binary black holes observed at the LIGO are expected to have antiparallel alignment, as implied by the inspiral spin parameters χ eff having values of approximately zero instead of the black hole's spin angular momenta cancelling each other, the orbital angular momentum of the binary black hole plays an important role in the spin angular momentum of the final black hole. Here, we will investigate the effects of the orbital angular momentum on the configuration of the final black hole and the gravitational wave in the case where M 1 M 2 and L orb M 2 1 , which is consistent with particle absorption. The secondary black hole is assumed to have a spin angular momentum of zero to show the contribution of the orbital angular momentum clearly. The variation in the orbital angular momentum in the initial state affects the final black hole and gravitational wave; hence, from Equation (24), Then, we impose an irreducible mass for the equality in Equation (27) because p r ≈ 0 and its variation becomes By combining Equations (38) and (39), the variation in the final black hole's mass with respect to the orbital angular momentum is obtained as which is exactly coincident with the contribution of the particle's angular momentum for particle absorption in Equation (11). Thus, the mass of the final black hole affected by the orbital angular momentum is obtained by a partial integration of L orb as where we impose the condition M f ≈ M 1 and a f ≈ a 1 . This represents the contribution of L orb to the mass of the final black hole. Owing to Equation (38), the change in the mass of the final black hole is opposite of the energy of the released gravitational wave with respect to the orbital angular momentum. Then, the orbital angular momentum contributes to the energy of the gravitational wave by The effect of the orbital angular momentum can be considered as an effective interaction potential of the orbital angular momentum in a black hole system. The accumulated energy for the potential is released in terms of the gravitational wave. Hence, the sign of the interaction potential is opposite to that of M gw ; hence, This is the spin-orbit interaction potential, which can be obtained in our model having the orbital angular momentum. In the coalescence of a binary black hole, the energy of the gravitational wave contributes as much as the interaction potential in Equation (43). Then, the antiparallel alignment between the primary black hole's angular momentum and the orbital angular momentum releases more energy of the gravitational wave than their parallel alignment owing to the contribution of the interaction potential because the interaction potential is U orb,int < 0 in the antiparallel case and U orb,int > 0 in the parallel case. This implies that the antiparallel alignment releases potential energy owing to the attraction but the parallel alignment needs to overcome its repulsion; thus, less energy is released compared with the antiparallel case. Note that attraction or repulsion can be easily shown by the sign of the interaction potential in Equation (43). Effectively, the energy released by the gravitational wave is equal to the interaction potential.

Gravitational Spin-Spin Interaction from MPD Equations
Here, we will investigate the effects of the spin parameter of the black hole. This will show that the alignment between black holes also plays an important role in the emitted gravitational wave. The procedure is similar to that in the previous section. For the collision of two black holes, in which one of them is a slowly rotating black hole with a small mass, M 2 M 1 and a 2 M 1 in the initial state. The gravitational wave with respect to the variation in a 2 becomes The variation in Equation (27) with respect to a 2 is Then, the variation in the mass of the final black hole with respect to a 2 is obtained as which implies that the mass of the final black hole is partially affected by the spin parameter a 2 . By integrating Equation (46) with Equation (44), the energy of the gravitational wave is partially emitted from the energy depending on the spin angular momenta, which is The spin interaction can be written in terms of the spin interaction potential as This is exactly coincident with the gravitational interaction potential predicted by the MPD equations in Equation (22). Thus, our approach for the binary black hole is consistent with particle absorption and the MPD equations in the particle limit of the secondary black hole. The spin interaction potential in Equation (48) changes its sign with J 1 and J 2 . The interaction is attractive for J 1 J 2 < 0 or repulsive for J 1 J 2 > 0. Then, the antiparallel alignment releases more energy than the parallel alignment during the coalescence of the binary black hole. Note that our analysis based on the irreducible mass provides the same result as the MPD equations and entropy-based analysis provided in Reference [53].

Energy of Gravitational Wave in Coalescence of Kerr Black Holes
We will numerically investigate the gravitational wave released in the collision of two Kerr black holes when M 1 ≈ M 2 . Since our approach from Equations (24)-(27) is consistent with the particle limits when M 1 M 2 , we will apply our model to various black hole pairs to investigate the energy of the gravitational wave. Further, the effects of variables such as the spin and orbital angular momenta will be studied. Once again, we assume the coalescence of two Kerr black holes such that where the second equation related to the radial momentum means that the black holes slowly come together in the radial direction; therefore, the irreducible mass of the system is almost conserved, as shown in Equation (37). This assumption is applied in the following sections. We now investigate the energy of the gravitational wave with respect to the alignments of the spin angular momenta during coalescence. Here, the final black hole and released gravitational wave energy M gw are obtained from Equations (24)-(29) by solving the numerical method. The effects of the spin parameters in the initial state are shown in Figure 1 with respect to the second black hole's spin parameter a 2 for a given value of the first black hole's spin parameter.  The alignments with rotating axes are dependent on the sign of a 1 a 2 in Figure 1. For a plus sign, the two black holes are in a parallel alignment and for a minus sign, they are in an antiparallel alignment. In Figure 1a, the released energy M gw is the largest at the extremal values a 2 where the rotational energy of the secondary black hole is maximum because M gw is emitted from the reducible energy such as the rotational energy. Hence, the amount of energy M gw is proportional to the rotational energy in the initial state. For the alignments, the released energy is greater in an antiparallel alignment than in a parallel one owing to the contribution of the gravitational spin interaction, as expected in Equation (47). Owing to the dependence on the alignment, the minimum point of the released energy is located at the parallel alignment for a given rotational axis of the primary black hole. The ratio of the released energy with respect to the initial mass is shown in Figure 1b. Interestingly, 0 < ≤ 5%, which is similar to that of the LIGO observations of about 3-4.5% [1][2][3][4][5]. This implies that our approach based on the irreducible mass can provide results consistent with observations. Further, this supports the fact that most of the released energy originates from the reducible energy of the initial state of the binary system. Note that this is an improvement of lowering the upper limits given in previous studies [52,53,67] to a realistic level.
The detailed effects of the alignment can be obtained from the alignment angle ψ in Equation (25). The alignment angle ψ shows the angle difference between the rotating axes of the primary and secondary black holes; thus, the parallel alignment corresponds to ψ = 0 and the antiparallel alignment is ψ = π. The energy released by the gravitational wave with respect to ψ is shown in Figure 2.
The energy released by the gravitational wave is the smallest at ψ = 0 and increases as ψ increases. Then, it attains a maximum at ψ = π. Hence, the attraction or repulsion due to the spin interaction plays an important role in the coalescence process. Further, the effects of the interaction are shown to be dependent on the angle ψ. The ratio of the gravitational wave energy also depends on the rotational energy for a given initial state. Although the primary black hole in Figure 2a has a spin parameter a 1 that is two times larger than that in Figure 2b, the ratio in Figure 2b is greater than that in Figure 2a because the primary black hole in Figure 2b is an extremal black hole, which has the maximum rotational energy for a given mass. Note that the rotational energy exponentially increases as the spin parameter approaches that of the extremal black hole. Hence, we can expect that more rotational energy included in the initial state in Figure 2b is released as a gravitational wave than that in (a). Note that the start and end points of Figure 1b are coincident with those of Figure 2b.  Since the energy of the gravitational wave is assumed to be released from the reducible energy such as the rotational energy, we investigate the ratio of the energy of the gravitational wave to the rotational energies of the initial state and final black hole. The ratios of the energy of the gravitational wave and final black hole's rotational energy with respect to the initial masses are defined as where the rotational energy is defined as M rot = M − M ir . Then, gw,rot indicates the rotational energy released by the gravitational wave with respect to the initial rotational energy and f,rot is the remaining rotational energy in the final black hole. The interaction potential still plays an important role in releasing the gravitational wave; thus, the antiparallel alignment emits more energy than the parallel alignment, as shown in Figure 3.
In addition, owing to the conservation of the total energy, gw,rot + f,rot = 1; thus, gw,rot and f,rot move opposite to each other. Compared with Figure 1, the energy of the gravitational wave in Figure 3a is not maximized at the extremal value of the spin parameter a 2 . Hence, the ratio related to the rotational energy is not exactly coincident with that of the total mass . Further, in Figure 1, the maximum value of gw,rot is located at M 1 a 1 + M 2 a 2 ≈ 0, where the final black hole is almost close to a Schwarzschild black hole with a zero angular momentum because most of the rotational energy is released by the gravitational wave. As a special case, when the first black hole is extremal, if the second black hole is also extremal and antiparallel, the ratio becomes maximum, as shown in Figure 3b. Incidentally, if we set a 1 = 0 and a 2 = 0, gw,rot and f,rot will diverge because there is no rotational energy in the initial state according to our assumption. Here, most values of gw,rot are quite large; therefore, we can expect that the final black hole is slowly rotating. However, the final black hole is at a f /M f ≈ 0.7 in the LIGO observations. Thus, we need more angular momentum to fill this rotational energy gap.  This is expected because a greater spin parameter for a given mass has more rotational energy in the initial state. The ratio also behaves in the same way, as shown in Figure 4b. Only the case of a 2 = 0 increases as M 2 increases, where the secondary black hole is a Schwarzschild black hole; thus, the increase in the secondary black hole's mass induces more rotational energy from the primary black hole.
The behaviors of the nonzero spin parameter cases are due to the amount of rotational energy in the initial state, as shown in Figure 5. For a secondary black hole with a fixed spin parameter, the energy released by the gravitational wave decreases as the mass of the secondary black hole increases in Figure 5a because the rotational energy of the secondary black hole decreases when the mass of the secondary black hole increases for a fixed spin parameter, as shown in Figure 5b.
The masses of the initial states simply affect the energy released by the gravitational wave, as shown in Figure 4. For a given mass of the primary black hole M 1 , the released energy increases as a 2 increases in Figure 4a.    Here, we will investigate the effects of the orbital angular momentum L orb during the coalescence of the binary black hole. The orbital angular momentum causes spiral motions of the black holes before coalescence. This can be more appropriately written as The orbital angular momentum is included in the total angular momentum of the black hole system and it makes similar contribution to the particle limit as the spin angular momentum does, as shown in Equation (43). Owing to the orbital angular momentum, the angular momentum of the final black hole is expected to increase because most of the spin angular momenta cancel each other in the LIGO observations, where the black holes have antiparallel alignment.
Hence, the orbital angular momentum increases the rotational energy of the final black hole, as shown in Figure 6. For a given L orb , the remaining rotational energy in the final black hole almost increases but is not exactly coincident because the released energy also depends on the alignment.
In Figure 6, for a primary black hole with a 1 = 5M , the orbital angular momentum of antiparallel cases is greater than those in the parallel cases. Further, there exist end points for a given initial condition in Figure 6 due to the extremal condition of the final black hole. This is easily checked in terms of the ratio in Figure 7. The orbital angular momentum L orb for a given initial state is proportional to the final spin parameter of the final black hole in Equations (24) and (26). This saturates the final state of the extremal black hole; thus, to avoid saturation, more mass remains in the final black hole. This leads to a decrease in the ratio . However, if a sufficiently large orbital angular momentum is set in the initial state, the total angular momentum is too large to form a final black hole. Thus, the end point = 0 implies that the final state is the extremal black hole, as shown in Figure 7.    The effects of the orbital angular momentum L erb with respect to the spin parameter a 2 contributes to the energy released by the gravitational wave, as shown in Figure 8.
According to the particle limit in Equations (42) and (43), there exists a gravitational spin-orbit interaction; thus, more energy is released by the gravitational wave in the antiparallel alignment than in the parallel one. The orbital angular momentum affects the energy of the gravitational wave with respect to the variation in the primary black hole's spin parameter a 1 , as shown Figure 8.  In Figure 8, the orbital angular momentum moves the location of the minimum of the ratio to the parallel alignment. In addition, the ratio attains a maximum for an extremal black hole with an antiparallel alignment. Similarly, the rotational energy released by the gravitational wave is maximized at the extremal black hole, as shown in Figure 9. However, the ratio gw,rot becomes small for a large value of the orbital angular momentum because the increase in L orb produces a final black hole with a large spin parameter a f ; thus, it becomes difficult to release the mass of the initial black holes by the gravitational wave to form the final black hole under the extremal bound. Here, the antiparallel alignment still releases more energy than the parallel alignment.

Gravitational Wave Observations
We have investigated the contributions of various variables in our model of the coalescence of the binary black hole. Here, by combining all of the information about these variables, we will analyse the initial and final states of black holes including the orbital angular momentum in LIGO observations such as GW150914, GW151226, GW170104, GW170608 and GW170814 [1][2][3][4][5]. In particular, we will mainly focus on the magnitude of the orbital angular momentum L orb , which has not been thoroughly studied in previous studies by thermodynamics. In addition, the initial spin angular momenta, including the orbital angular momentum, will be estimated by applying our approach. By applying our approach to GW150914, we will introduce a general procedure for our investigation of the LIGO observations. According to our investigation, our analysis produces similar median values of spin and orbital angular momenta as the LIGO observations. Then, a similar analysis will be applied to other observations. Note that we will use source-frame masses related to detector-frame masses by applying the source redshift z.

GW150914
The source of GW150914 is a binary black hole merger with the effective inspiral spin parameter and a f /M f = 0.67 +0.05±0.01 −0.07±0.02 with a released gravitational wave energy M gw /M = 3.0 +0.5 −0.5 of about 4.6% of the total mass. Since the parameter ranges of a 1 and a 2 are much larger than that of a f , we will find the proper parameter ranges of a 1 and a 2 by applying our approach. Finally, in agreement with the model, the orbital angular momentum and released angular momentum will be estimated. Here, our input values including errors are set to Then, the ranges of a 1 , a 2 , L orb , J tot and J gw will be obtained with respect to the parameter range 22 . We now utilise these observation data in our model. The basic framework is the same as Equations (27)- (30) in Section 3. Since the spin parameters of the initial black holes are dependent on the effective inspiral parameter χ eff , our approach should be modified to add this parameter. In our definitions of the parameters, the effective inspiral parameter is given as where the vector direction of the orbital angular momentum is fixed and |L orb | = +1. Further, we assume in Equation (52) that the axes of the spin angular momenta are already aligned parallel or antiparallel to the axis of the orbital angular momentum. Hence, a positive spin parameter implies that the spin and orbital angular momenta rotate in the same direction, and a negative spin parameter implies that the spin and orbital angular momenta rotate in the opposite directions. In addition, we already assumed the conservation of the irreducible mass. Then, from Equation (27), In combination with Equations (52) and (53), the magnitudes and alignments of a 1 and a 2 can be obtained. Then, the magnitude of the orbital angular momentum can be calculated from the ratio of the gravitational wave ≈ M ≈ J , which is rewritten as From Equation (54), the orbital angular momentum becomes This implies that the absence of the total angular momentum between the initial and final states with the released one is supplied from the initial orbital angular momentum. This is supported from the observations of GW150914. (i) There exists a difference between the sum of the spin angular momenta of the primary and secondary black holes and the final black hole. Approximately, M f a f − (M 1 a 1 + M 2 a 2 ) (62M ) 2 · 0.67 − ((35.8M ) 2 · 0.32 + (29.1M ) 2 · 0.44) +1800M 2 , (56) where an additional angular momentum of at least 1800M 2 is needed for the final black hole to satisfy the conservation of the angular momentum. (ii) Since the median value of χ eff −0.07 is around zero, the initial spin angular momenta of the initial black holes almost cancel each other; thus, the contribution of the spin angular momenta is very limited. Then, the spin angular momentum of the final black hole has to be supplied from other angular momenta in the initial state. Except the spin angular momenta, the only remaining angular momentum is the initial orbital angular momentum of the binary black hole system.
Satisfying Equations (52)-(54), the initial angular momenta can be divided into two cases, a 1 ≥ 0 and a 1 < 0, as shown in Figures 10 and 11. Because χ eff 0, the spin parameters of the initial black holes are in antiparallel alignment in Figures 10a and 11a; therefore, the orbital angular momentum is of similar magnitude as the spin angular momentum of the final black hole in Figure 10b. The orbital angular momentum is given as a parameter with respect to the total mass of the initial state; therefore, a orb = L orb (M 1 +M 2 ) , as shown in Figures 10b and 11b. By considering the ranges of the input values, we can obtain the magnitudes of the estimated values of where the median values are at χ eff = −0.07. The estimated spin parameters of the initial black holes are slightly larger than those of Reference [73] but within their error ranges. Further, as predicted from our approach, the ranges of spin parameters are tighter than those in Reference [73] and ensure the conservation of the angular momentum. Here, we newly estimate the sum of the total angular momentum in the initial state J tot , which is a very large value compared with the initial spin angular momenta. This is based on the fact that the initial spin angular momentum is almost cancelled owing to the alignment of the initial black holes. Hence, most of the spin angular momentum of the final state including the gravitational wave is provided from the orbital angular momentum of the binary black hole. This could be one reason why the frequency of the gravitational wave is almost proportional to the orbital frequency of the binary black hole system before coalescence. Thus, the orbital angular momentum is crucial in the coalescence of the binary black hole and the released gravitational wave.    We will repeat the same procedure for other observations of gravitational waves in the following subsections.

Summary
We have investigated the coalescence of a binary black hole with a released gravitational wave by constructing a model using Kerr black holes with an orbital angular momentum. In particular, this construction is expected to provide a more detailed analysis of the spin and orbital angular momenta of the binary system. Located far from each other in the initial state, two Kerr black holes having orbital angular momenta slowly come together to form the final Kerr black hole.
In the basic framework, we apply three assumptions in our approach: (i) the conservation of energy; (ii) the conservation of angular momentum and (iii) the conservation of the irreducible mass. Since we consider the irreducible mass, the mass of the Kerr black hole is divided into irreducible and reducible masses. Since our model should be coincident with particle absorption for a Kerr black hole, we obtain the constraints M ≈ J ≈ and identify that the irreducible mass can be approximately conserved for slowly moving black holes. Owing to the conservation of the irreducible mass, the upper limit of the energy released by a gravitational wave is very close to the actual energy of the gravitational wave; hence, we assume that M gw,upper ≈ M gw . Further, analytical descriptions of the gravitational spin-orbit and spin-spin interactions are obtained for a variation of our model. Interestingly, these analytical forms of U orb,int and U spin,int exactly correspond to the results from particle absorption and the MPD equations.
Under the constraints and from an analytical analysis of particle absorption and the MPD equations, we have numerically applied our approach for the coalescence of binary black holes having equal masses. In this case, the range of energies of the released gravitational wave is about 3.0-4.5%. This is a very important advantage of our model because most of the LIGO observations support that the released energy ratio of the gravitational wave is about 4.5%. Therefore, we could provide a more realistic analysis using our approach based on a simple thermodynamic description. The energy released by the gravitational wave depends on the alignments of the spin and orbital angular momenta in the initial state. Since the dependency of the alignment is the same as those in the gravitational spin-orbit and spin-spin interactions, as we expected, the released energy in an antiparallel alignment is greater than that in a parallel alignment for a fixed one of angular momenta. In addition, owing to the formation of the final black hole, the orbital angular momentum is limited and has the maximum value.
Finally, we apply our approach to five LIGO observations of binary black hole mergers. We have estimated the parameters of the initial state using other parameters having small error ranges. Since the inspiral spin parameters are around zero in most of the observations, the large values of the spin angular momenta cancel each other. Hence, the spin angular momenta of the final black hole and gravitational wave have to be provided from the orbital angular momentum, which becomes a large value compared with the spin angular momentum. Using our approach, we have obtained a 1 , a 2 , a orb , J tot and J gw .
We have shown that the binary black hole merger can be approximated as the coalescence of two Kerr black holes having an orbital angular momentum with consideration of the irreducible mass. This implies that the energy source of the gravitational wave is the reducible energy such as the rotational energy and kinetic energy included in the mass of the Kerr black hole.