arXiv:1711.03837v2 [physics.gen-ph] 27 Feb 2019 o lc oewstetdt eaqatt hc spootoa to proportional is which quantity a as be expressed to usually treated is previou form[5] was the to hole Bekenstein-Hawking’s analogous thermodynamic black was hole a study black and for The of evaporate calculated. strong aspects can were different as holes quantities and treated black Black people that be like shown attracts can objects has mergers compact Hawking hole theoretically the black of regarding detections[1] theories wave many exist There Introduction 1 pr Thermodynamic Thermodynamics; physics; hole Black : Keywords adpDutta Sandip Assum Hole Time Black With Schwarzschild Increasing a is of Mass Interior the for Entropy it ub uda,Sae:Ws egl India Bengal, West : State Burdwan, Gola Purba Burdwan, of : University Dist. The Mathematics, of Department h nes eprtr a ae.()Fnly eueteflxrela obt flux we the use free we Finally, the (c) From (b) taken. was energy. temperature the inverse as the identified is E where iedpnetbcgon,tesaa edwscoe as chosen In was Kerr[ (a) field for scalar description: done the clear was very background, calculation without dependent hole This paper[10] time black Schwarzschild the volume). the in taken CR of ho which volume as black inside (known extremal the Rovelli do defining and of renewed: path recently new has microsta a which entropy a question, the entropy this the to of by corresponds source question that the a is microstates Yet, what classical 90’s[7]. for holes: the known black in well of executed A thermodynamics later was thermodynamics[6]. charge of Noether laws of ordinary the to similar laws where where nrae ihtm.Aohrntrlnto fetoyi h o N Von the is entropy of notion natural Another time. with increases 2 1 [email protected] [email protected] hntevlct ftecniee ytmi oprbet h spe the to comparable incorporated. is thermo is system in charge considered change scalar the of The of velocity considered. the temperatur are when inverse thermo quantities and of temperature thermodynamic law between second to relation the a Bekenste get satisfies the we entropy to law th evaluated proportional calculate the is entropy that we the t see that Here assuming find modes, We massless horizons. time. for event with hole black hole Schwarzschild a black of of interior existence the with ρ A stedniymti faqatmsse.I lsia ee,w know we level, classical In system. quantum a of matrix density the is steae ftebakhl’ vn oio.Tescn a ftherm of law second The horizon. event hole’s black the of area the is lc oetemdnmc steae fsuyta ek orecon to seeks that study of area the is thermodynamics hole Black 1 n iart Biswas Ritabrata and exp 2 {− S [ iET ρ S = ] BH Abstract } − exp = 1 tr 4¯ { ( lnρ ρ hG iI A ( e?S,w on ogsadn su nliterature in issue long-standing a found we So, tes? bgAaei ope,,Ct Burdwan-713104, : City , Complex, Academic pbag ,θ φ θ, λ, , ) , h lsia ee?Mr eeal,wa r the are what generally, More level? classical the t a rie suasee nteae fclassical of area the in unanswered as arrived has xadn h li-odneuto ne the under equation Klein-Gordon the expanding inHwigexpression)[10]. tion(Hawking []hsbe icse.Mn thermodynamic Many discussed. been has s[4] vdneo hi xsec.Frtefis time, first the For existence. their of evidence tde fcascltemdnmc.Entropy thermodynamics. classical of studies s e aezr rnnzr nrp[] Recently, entropy[8]? non-zero or zero have les uanetoygvnby given entropy eumann i h nrp,tedrvtv ihrsetto respect with derivative the entropy, the ain ) } cse;Etoy yai lc Hole Black Dynamic Entropy; ocesses; , ev hra aito[,3.Ti concept This 3]. radiation[2, thermal leave oe.Ee h rsn a gravitational day present the Even holes. yais sn h thermodynamic the Using dynamics. .Teseilrltvsi corrections relativistic special The e. nHwigepeso.As,w can we Also, expression. in-Hawking uafrteetoyo lc oei terms in hole black of entropy the for mula a enitoue[4 yChristodoulou by introduced[14] been has t vricesn rao vn horizon. event of area increasing ever its ] o nrdcn h olwn inputs following the introducing Now 9]. ems ftebakhl increasing hole black the of mass he do ih.Teeeto presence of effect The light. of ed yai rprisaeanalyzed are properties dynamic ietelw fthermodynamics of laws the cile nrp orsodn othe to corresponding entropy e htbakhl ehnc olw the follows mechanics hole black that dnmc ttsta h entropy the that states odynamics n the ing Here we tried to get a clean knowledge about the above inputs. Express the interior volume as integral form[10] and defined an effective metric. Then we evalute the Hamiltonian for a particle moving in this background. Finally, we use the Gibb’s free energy for a massless particle and find the entropy as time dependent which is a very small impact on entropy changes. The structure of the paper is as follows: firstly, we write a brief review of the results, given in [10]. In section 3, we calculate the energy for a massless particle within the interior of the black hole. In section 4, we calculate the energy and in section 5, we consider a system with speed comparable to that of light and study the dependence v of the corresponding thermodynamics as a function of the ratio c . Section 6 will study the effect of scalar charge present and finally we conclude in section 7.

2 Review of Previous Works

First we write the well known Schwarzschild metric in Eddington-Finkelstein coordinates

ds2 = fdv2 +2dvdr + r2dΩ2 , (1) − 2 n 2M(1+a1t+a2t +...+ant ) where the function f = f(r, t)=1 r (where, 0

2 considering the cross term as zero by taking the transformation properly. Assuming that the condition f ∂v + − ∂T 2 ∂v ∂r = 1 is enforced, if the spherically symmetric hypersurface is considered as the direct product of a 2-sphere ∂T ∂T  and an arbitrary− curve parametrized by λ in v r plane then we are able to find the hypersurface Σ : T = constant [10, 14] where we have − ds2 = dT 2 + fv˙ 2 + 2v ˙r˙ dλ2 + r2dΩ2 . Σ − − The interior volume within the horizon can be written by the surface Σ γ S2 on which metric can be written as ≡ × ds2 = fv˙ 2 + 2v ˙r˙ dλ2 + r2dΩ2 , − where, r = r(λ) and v = v(λ), λ being an arbitrary parameter[14]. The volume can be written as V =4π dλ r4 fv˙ 2 + 2v ˙r˙ . (2) Σ {− } Z Now, the metric takes the form of an integrand of Lagrangian,p

dS2 = r4 fv˙ 2 + 2v ˙r˙ dλ2 (3) eff − 3choosing these variables negative such that the second law of thermodybamics is violated is nothing but to chose an unphysical system.

2 For this Lagrangian, the coordinates are (r, v) and the momenta are (Pr, Pv) respectively. Now, drdvdPr dPv represent the phase space volume. In this paper we evaluate the entropy inside this phase space volume. In the process of evaluation of entropy in a quantum statistical way, firstly we have to calculate the HamiltonianR of a particle, restricted inside the volume.

3 Calculation of the Hamiltonian of a Particle

Considering m as the mass of a particle which moves in a space time with a back ground metric given by

ds2 = g dxadxb = dt2 + r4 f(r, t)dv2 +2dvdr , (4) ansatz ab − − the action (which supposed to have the reparametrization symmet ry) seems to be 

2 2 1 a b 2 S = m dSansatz = m gabdx dx . (5) Z1 Z1  a dxa a a The velocities of the particle is given by u = dτ , where, τ is an arbitrary parameter and x = x (τ). The path of the particle is then given by 1 2 2 dxa dxb 2 S = dτ = m g dτ . (6) L ab dτ dτ Z1 Z1   1 a b 2 Comparison of both sides clear states the Lagrangian to have the form = m g dx dx . By using the Euler L ab dτ dτ Lagrange equation, we can easily find the equation of the motion of the system as  d2xa dxb dxc dua +Γa =0 +Γa ubuc =0 , dτ 2 bc dτ dτ ⇒ dτ bc a where Γbc is the Christoffel symbol. Also the above geodesic equation is true for any space time. To evaluate the Hamiltonian which describes the whole dynamical system, firstly, we have to calculate the momenta of the system given by ∂ m2 dxb P = L = g . (7) a ∂x˙ a ab dτ L Therefore, the canonical Hamiltonian is dxa H = P =0 . (8) c a dτ − L In reparametrization invariant theory, canonical Hamiltonian was zero which is a typical signature of this theory; for Minkowski space time the same thing can be observed. For the present article, we have to analyse the Hamiltonian of the system given by (4). For chronological gauge [15], in flat space time the analysis is done and for proper time gauge in [16]. For calculation of entropy, we present the constraint which are independent. Since the momenta are not independent,

2 ab 2 P = g PaPb = m . (9) Then we get primary constraint Φ= P 2 m2 0 (10) − ≈ Using Dirac’s algorithm [17], the primary constraint and the Hamiltonian are proportional to each other, therefore

H = ξ(τ)Φ = ξ(τ) P 2 m2 , (11) T − where ξ depends on τ and is known as proportionality constant.  Hence dxa x˙ a = = ua = xa, H =2ξP a (12) dτ { T } and ∂H ∂gbc P˙ = P , H = = ξ P P (13) a { a T } −∂xa − ∂xa b c 2 also using (7) and (12) we get, ua =2ξP a =2ξ m ua. So we have L ξ = L . (14) 2m2

3 Therefore, the total Hamiltonian is H = L P 2 m2 . (15) T 2m2 − So far, we have considered one constraint only (given by (10)). Th is characterises the system as of first class and hence it has gauge freedom [18]. Imposing some conditions on the arbitrary parameter τ which we are going to interpret as proper time, the gauge freedom can be removed. A proper time gauge will be imposed throughout the subsequent analysis. Therefore we are now with P 0 ψ = τ x0 0 (16) 2 m − ≈ and for primary constraint (11) we have, ψ = P 2 m2 0 , (17) 1 − ≈ this makes the system as second class. Therefore, ∂ψ P 0 τ ψ˙ = 2 + ψ , H =0 x0, H + P 0, H =0 . (18) 2 ∂t { 2 T } ⇒ m − T m T Also   0 0 2 ab 0 0 x , HT = x , ξP =2ξg x , Pa Pb =2ξP (19) and    ∂g0b ∂gbc P 0,H = ξ 2P 0 P g0a P P . (20) { T } ∂xa b − ∂xa b c   For the metric (4), we get 0 P ,HT =0 . Therefore equation (18) and (19) gives  1 ξ = . (21) 2m From equations (14) and (21) we can obtain = m. Also from equation (11) and (12) we get,L 1 P a x˙ = xa, P 2 = . (22) a 2m{ } m and 1 1 ∂gab ∂gdb P˙a = P a, P 2 = 2 gcd gca P P . (23) 2m{ } 2m ∂xc − ∂xc d b   1 P˙ a = Γa P bP c . (24) −m bc here we useing the identity ∂ gab = gadg ∂ g c − be c de for two dynamical variables, caculating the Dirac braket [9, 19]

f ,f ∗ { 1 2} is 1 2 2 f ,f ∗ = f ,f + f , P ψ ,f f , ψ P ,f . (25) { 1 2} { 1 2} 2P 0 { 1 }{ 2 2}−{ 1 2}{ 2} for a dynamical variable, the equation of motion can be obtained by the following relation [17]

f˙ = f ,H ∗ . (26) 1 { 1 } In fact H = P 0 serves our motive.

ab a a 0 a0 1 a τ a0 0 b c a a 0 ∂g 0c 1 a b c 0a 0 b c x˙ = x , P ∗ = g + P + g Γ P P and P˙ = P , P ∗ = g P Γ P P + g Γ P P { } P 0 m bc { } ∂xc b − P 0 bc bc   (27)

4 By the gauge fixing constraint a =0 has been eliminated, then the equation of motion for the space component of

a(a = µ) is P µ 1 x˙ µ = and P˙ µ = Γµ P aP b (28) P 0 P 0 ab From the above two equation, we can eliminate P 0 and get the desired geodesic equation. Now,

ab 0 2 4 4 2 2 g P P = P +2r− P P fr− P = m (29) a b − r v − r The energy of a particle, i.e., the Hamiltonian is given by

1 fP 2 2P P 2 ǫ = P 0 = r + r v m2 (30) − r4 r4 −   We are calculating the energy for a massless particle (like photon). So we can consider m 0, then the equation (30) can be written as → 1 fP 2 2P P 2 ǫ = r + r v (31) − r4 r4   4 Entropy Calculation for Classical Gibb’s free energy

Now, we will construct the expression for entropy w.r.t. the energy for the massless particle. This entropy is defined 1 inside the black hole. Also, no chemical potential terms are present, hence Gibb’s free energy is G0 = β ln Z = 1 − β Σǫ ln (1 exp βǫ ), where Z is the grand canonical partition function and β is the inverse temperature. Gibb’s free energy− is not{− a relativistic} invariant. The increment in mass can be considered as a slow process. So, we can consider that the effect of Gibb’s free energy is negligible, i.e., the differences between Gibb’s free energy before and after the change of hypersurface will be almost the same. In horizon, for maximum volume of the interior, we get the fixed value of the radial coordinate, i.e.,r ˙ = 0, where r can be obtained from the equation (28) which implies 4 (Pv +fPr ) r˙ = r− P 0 =0 Pv + fPr = 0. Again r and ingoing null coordinates are depending on λ(a parameter)[14], i.e., r = r(λ) and v =⇒v(λ), v = F (r), a function of r. Therefore the Gibb’s free energy

1 G = dP dP drdv ln (1 exp βǫ ) δ (P + fP ) δ(v F (r)) . 0 β r v − {− } × v r − Z Using Dirac-Delta functions

1 1 fP 2 2 G = dP dr ln 1 exp β r , 0 β r − − − r4 Z " (   )# where P varies from 0 to . r ∞ So 1 2 2 1 ∞ fP G = dr ln 1 exp β r 0 β − − − r4 Z Z0 " (   )#

βPr√ f − incorporating r2 = x we have π2 r2 G0 = dr (32) −6β2 √ f Z − Let us consider as time goes the radius of the black hole increases, then the radial coordinate varies from 2M(1 + 2 n a1t + a2t + ... + ant ) to zero. So, we have

π2 0 r2 G0 = dr . 2 2 n −6β 2 n 2M(1+a1t+a2t +...+ant ) Z2M(1+a1t+a2t +...+ant ) 1 r − q 5 r Incorporating y = 2 n , we have 2M(1+a1t+a2t +...+ant )

5π3 G = M 3(1 + a t + a t2 + ... + a tn)3 0 12β2 1 2 n

Hawking showed that quantum effects allow black holes to emit exact black body radiation[2]. The electromag- netic radiation is produced as if emitted by a black body with a temperature inversely proportional to the mass of the black hole. Also, our system is like black body which encloses massless particles. Hence the particles are also at inverse temperature β =8πM. Therefore, the Gibb’s free energy 5 G = β(1 + a t + a t2 + ... + a tn)3 (33) 0 6144 1 2 n Hence, the entropy is ∂G 5β2 S = β2 0 = (1 + a t + a t2 + ... + a tn)3 (34) 0 ∂β 6144 1 2 n So, we can say that the entropy is monotonically increasing as time increases but the increment of entropy depends on time which is very small. In this connection we wish to reestablish the result given by [20], where they have proved that Planck entropy is a Lorentz invariant. The empirical temperature in relativistic thermodynamics is a Lorentz invariant[21]. Let us assume in the proper frame, T0, δQ0 and δL0 are the absolute temperature, the amount of heat entering the system and the thermodynamic work done over the system respectively. For reversible processes

δQ0 = T0dS0. (35)

At the beginning if we had assumed that the efficiency of the Carnot cycle is an invariant, then the entropy is an invariant, i.e., δQ δQ0 dS = = = dS0. (36) T T0

By the second law of thermodynamics the Pfaffian form δQ0 depending only on the empirical temperature. This leads to the Clausius equality δQ 0 =0 (37) T I 0 for every cyclic reversible process. Let us assume that the second law of thermodynamics is in force in relativistic thermodynamics in any inertial frame. Hence, the absolute temperature of the form

T = T (T0, v) , lim T (T0, v)= T0 , (38) v 0 → where the thermodynamic system moving with a constant velocity ~v with respect to a frame K. Also, the Clausius equality in K reads δQ =0. (39) T I Again, the relativistic transformation of heat found by H. Ott[22] and independently by H. Arzelies[23] and C. Mφller[24, 25] δQ = γ(v)δQ0 (40) (see also [26]). Substituting (40) and (38) into (39), we obtain,

δQ 0 =0. (41) γ(v) 1T (T , v) I { }− 0 Comparing (41) and (37), we get, 1 γ(v) − T (T , v)= bT , (42) { } 0 0 where b is an arbitrary constant. Taking v 0 in the equation (42) and applying (38), we get, b = 1. So finally →

T = γ(v)T0. (43)

6 This is the second law of relativistic thermodynamics for reversible processes and the relativistic transformation of absolute temperature consistent with (40). This transformation rule and the transformation law are different[27, 20, 28, 29, 30] (P lanck) 1 T = γ(v) − T (44) { } 0 following from the transformation rule[31, 32, 33, 34] (L) T = T0 (45) consistent with (L) δQ = δQ0. (46) Now, from the equation (35), (40) and (43), we get,

δQ = T dS0 (47) which implies that according to Planck[20] entropy is a Lorentz invariant

S = S0 (48) Some important things we have to mention here : firstly, we have considered that the massless modes inside the horizon are at a temperature whose value is same as that of the event horizon. We should remember that r varies from 0 to 2M but our consideration was slightly different than this. We have assumed the mass of the black hole is increasing as time increases because matters get accreted towards a black hole with time. But the average density of the universe is very small. So the increment of the mass of the black hole is very slow. For this we assumed that the radial coordinate’s increment also is very minute. This leaded us to the consideration that the mass of the black 2 n hole is M(1 + a1t + a2t + ... + ant ) after a time t while M = M(t = t0 = 0). We have also assummed that r varies 2 n from 2M(1 + a1t + a2t + ... + ant ) to zero. Also, Hawking showed that black holes emit thermal Hawking radiation[5, 35] corresponding to a certain tem- perature(Hawking temperature)[2, 36]. Using the thermodynamic relationship between energy, temperature and entropy, Hawking was able to confirm Bekenstein’s conjecture and fix the constant of proportionality at 1/4[37]: A 2 SBH = 4 is the entropy on the horizon of the black hole, where A = 16πM is the area of the event horizon. Hence, the relation between S0 and SBH is 5π S = (1 + a t + a t2 + ... + a tn)3S 0 384 1 2 n BH

Thus we can see that the entropy on the horizon is greater than the entropy inside the black hole i.e. SBH >S0. Also the entropy inside the black hole is proportional to the time(t), i.e., if we increase the time then the entropy(S0) inside the black hole increase. The first law of thermodynamics provides the basic definition of internal energy, associated with all thermody- namics systems, and states the rule of conservation of energy. The second law is concerned with the direction of natural processes. It asserts that a natural process runs only in one sense, and is not reversible. Also, we can see from the above dS 0 > 0 dt . Hence, the total entropy can never decrease over time and the process is irreversible. According to the first law of thermodynamics

dM = T0dS0 (49) 768 1 β = 2 n 2 , ⇒ 15πT0 (1 + a1t + a2t + ... + ant ) where β is the reciprocal of the thermodynamic temperature of the system[38, 39]. Thermodynamic beta is essentially the connection between the information of a physical system through its entropy and the thermodynamics associated with its energy. Case:I 2 n If ai = 0 for i> 1 then (1 + a1t + a2t + ... + ant ) (1 + a1t) and 5π 3 768 1 ≈ 2 S0 = 384 (1 + a1t) SBH and β = 15πT0 (1+a1t) If we increase the entropy on the horizon of the black hole (SBH ) then we show that the entropy inside the black hole increases here the time assumed as a constant. On the otherhand, if we variate the time then the rate of increasing of entropy inside the black hole become higher than the above case. The relation between β and T0 is hyperbolic through out the figure, when time is not constant the eccentricity increase with time.

7 Fig.1a

Relation between S0, SBH and t

Fig.1b

Relation between β, T0 and t

Case:II 2 n 2 If ai = 0 for i> 2 then (1 + a1t + a2t + ... + ant ) (1 + a1t + a2t ) and 5π 2 3 768 1 ≈ 2 2 S0 = 384 (1 + a1t + a2t ) SBH and β = 15πT0 (1+a1t+a2t ) When time constant the entropy inside the black hole increases with the increament with the horizon entropy. But if the time is not constant the rate of increasing entropy inside the black hole is distincity high. Here the rate of increasing of the eccentricity is higher than the Case:I.

8 Fig.2a

Relation between S0, SBH and t

Fig.2b

Relation between β, T0 and t

Case:III 2 n 2 3 If ai = 0 for i> 1 then (1 + a1t + a2t + ... + ant ) (1 + a1t + a2t + a3t ) and 5π 2 3 3 768≈ 1 2 3 2 S0 = 384 (1 + a1t + a2t + a3t ) SBH and β = 15πT0 (1+a1t+a2t +a3t ) . If the time is not constant the rate of increasing entropy inside the black hole is distincity high from the Case:I and Case:II.

9 Fig.3a

Relation between S0, SBH and t

Fig.3b

Relation between β, T0 and t

Here the rate of increasing of the eccentricity is higher than the Case:I and Case:II.

5 Phenomenological Relativistic Thermodynamics

In proper inertial frame (i.e., say K0), the relativistic counter part of famous Maxwell’s distribution for ideal gas was proposed by F. Juttner[30]. Throughout a long period, Juttner distribution had been accepted by many literatures[29, 40, 41, 42, 43, 44]. When a vessel, kept at labaratory, contains gas moving with constant velocity ~v with respect to the labaratory frame K, the similar assumptions and deductions like Juttner has been constructed. Some literature based on computer simulations[45, 46, 47, 48, 49], supports partially the correctness of the Juttner formula and its generalization. Some other distributions also have been proposed. Though Juttner distribution is found to be an efficient one[50]. The relativistic Gibb’s distribution for an ideal gas contained in the vessel moving j ~v with a constant velocity ~v and 4-velocity u = γ(v) c ,γ(v) turns to be 1  dω = exp βcu j d3N pd3N q , (50) (2π)3N N!Z {− jP }

10 where N is the number of particles, β = 1 with κ being the , j is the total momentum of κT0 P the gas, d3N pd3N q is the phase space volume element and Z denotes the partition function

N V j 3N Z = 3N exp βcuj d p . (51) (2π) N! 3N {− P } ZR Again, the measure d3p transforms as follows[41]

3 ~v ~p0 3 d p = γ(v) 1+ · d p0 . (52) 2 2 2 c ~p0 + m c !

Hence, p N N V j 3 Z = 3N exp βcujp d p (2π) N! 3N {− } ZR  N N V 2 2 2 ~v ~p0 3 = exp βc ~p0 + m c γ(v) 1+ · d p0 3N 3 2 2 2 (2π) N! R N − " c ~p + m c !# ! Z n q o 0 N p N V0 2 2 2 3 = 3N exp βc ~p0 + m c d p0 = Z0, (53) (2π) N! R3N − Z n q o  where the subindex ‘0’ corresponds to the proper frame K0. From the above equation, we can say that the partition function Z is a Lorentz invariant[41]. Finally, the relativistic Gibbs distribution (50) turns to be 1 ω = exp βcu j , n =1, 2, ... (54) n Z {− jPn} where the subindex ‘n’ denotes a quantum state. Also, from the first law of relativistic thermodynamics (for ideal liquid)[51, 21] v2 dE = T dS pdV + γ(v) 2 V dp, (55) − { } c2 From (55), one can easily find v2 d E γ(v) 2 V = T dS γ(v) 2pdV (56) −{ } c2 p −{ }   Hence, ∂E v2 ∂p T = γ(v) 2 V , (57) ∂S −{ } c2 ∂S  V  V v2 ∂p ∂E p γ(v) 2 V = (58) −{ } c2 ∂V − ∂V  S  S and ∂T ∂p = γ(v) 2 (59) ∂V −{ } ∂S  S  V Therefore, from the equation (55) and (56), we obtain,

∂E T = γ(v) 2 (60) { } ∂S  V,P

Also, pressure is a Lorentz invariant[21], i.e., p = p0 and as the volume V0 = γ(v)V we can write

= E + pV = γ(v) (E + p V )= γ(v) , (61) H 0 0 0 H0 where is the enthalpy. UsingH the definition of enthalpy and equation (56), we get,

∂H T = (62) ∂S  p

11 Therefore, the free energy F v2 F = E TS γ(v) 2 pV = γ(v)F , (63) − −{ } c2 0 where F0 = E0 T0S0 is the free energy of the system in its proper frame K0. Using the equation− (63) into the equation (56), we get,

dF = SdT γ(v) 2pdV (64) − −{ } Therefore, ∂F ∂F ∂F v2 ∂F S = , γ(v) 2p = , E = F T V (65) − ∂T { } − ∂V − ∂T − c2 ∂V  V  T  V  T and the relativistic Maxwell identities ∂S ∂p = γ(v) 2 (66) ∂V { } ∂T  T  T We define G as the Gibbs function in a standard way

G = E TS + pV = γ(v)G (67) − 0 with G0 = E0 T0S0 + p0V0 Again, some− invariant thermodynamical quantities are represented by the Massieu function[52],

G G = 0 (68) T T0 Using the equations (33) and (67), we get,

5γ(v) G = β(1 + a t + a t2 + ... + a tn)3 (69) 6144 1 2 n Therefore, the entropy is

∂G 5γ(v)β2 S = β2 = (1 + a t + a t2 + ... + a tn)3 = γ(v)S (70) ∂β 6144 1 2 n 0

According to the first law of thermodynamics and using the equation (70), we get,

1 dM = T dS = γ− T0γdS0 = T0dS0, (71)

1 since T = γ(v)− T0 . The above equation is same as the equation (49). (see also [52]). A 2 SBH = 4 is the entropy[37] on the horizon of the black hole, where A = 16πM is the area of the event horizon. Hence, the relation between S and SBH is 5πγ(v) S = (1 + a t + a t2 + ... + a tn)3S , (72) 384 1 2 n BH where γ(v) = 1 and v is the velocity between both the rest and moving reference frames. Let us assume √1 v2/c2 − v = α, where 0 α< 1. So, from the equation (72) c ≤ 5π√1 α2 S = − (1 + a t + a t2 + ... + a tn)3S (73) 384 1 2 n BH Case:I 2 n If ai = 0 for i> 1 then (1 + a1t + a2t + ... + ant ) (1 + a1t) and 5π√1 α2 3 ≈ S = 384− (1 + a1t) SBH

12 Fig.4a

Relation between S, SBH and t for α = 0.5

Fig.4b

Relation between S, SBH and t for α = 0.9

We have drawn the newly calculated entropy with respect to SBH and time for α = 0.5 and 0.9 for a1 = 0.05 in fig 4a and 4b respectively. We see as we increase α, the value of entropy decreases. This means as the system is moving with a speed near to the speed of light, the variation of entropy decreases lower than the classical terminal results. Case:II 2 n 2 If ai = 0 for i> 2 then (1 + a1t + a2t + ... + ant ) (1 + a1t + a2t ) and 5π√1 α2 2 3 ≈ S = 384− (1 + a1t + a2t ) SBH

13 Fig.5a

Relation between S, SBH and t for α = 0.5

Fig.5b

Relation between S, SBH and t for α = 0.9

We have plotted the newly calculated entropy for α = 0.5 and 0.9 for a1 = 0.05, a2 = 0.04 and we can see in this case with time the increment of S is higher than the classical case. Case:III 2 n 2 3 If ai = 0 for i> 2 then (1 + a1t + a2t + ... + ant ) (1 + a1t + a2t + a3t ) and 5π√1 α2 2 3 3 ≈ S = 384− (1 + a1t + a2t + a3t ) SBH We have plotted S vs SBH and t for α = 0.5 and α = 0.9 in the figures 6a and 6b respectively. When a1, a2 and a3 all are non zero we observe that the rate of increment of S with time is higher than the previous cases.

14 Fig.6a

Relation between S, SBH and t for α = 0.5

Fig.6b

Relation between S, SBH and t for α = 0.9

Fig.4a, fig.5a, fig.6a and fig.4b, fig.5b, fig.6b represents the relations between S, SBH and t for α = 0.5 and α = 0.9 respectively. Also, the figures which we will get in Case:I, Case:II and CaseIII for α = 0, are identical as fig.1a, fig.2a and fig.3a respectively in section 4. The medium moves with a constant velocity ~v with respect to the inertial frame K. If the velocity ~v is negligible to the speed of light that case represented by taking α = 0. Also, taking α =0.5 for those medium which moves with comparable velocity ~v to the speed of light. Again, if the medium moves with a velocity near to the speed of light then consider the case for α =0.9. From the figures we can say that if we increase the velocity of the medium then the entropy decrease very fast, i.e., if we decrease the value of α then the entropy increase.

6 Black Holes with Scalar Charge

Let us consider a new type of black hole solution, namely, massive gravity black hole[53]. The ansatz for the static spherically symmetric black hole solutions can be written in the following form: dr2 ds2 = g(r, t)dt2 + + r2 dθ2 + sin2θdφ2 , (74) − g(r, t)

2 n  2M(1+a1t+a2t +...+ant ) Q where g(r, t)=1 r rλ and λ is a parameter of the model which depends on the scalar charge Q. The presence− of the scalar charge− represents a modification of the Einstein’s gravitational theory. In this

15 paper we shall focus on the case Q> 0 and λ = 1. From the equation (32), we get, the Gibb’s free energy as

3 π aβ + Q 4π 1 aβ 1 Q 1 Q G = 60 sin− + Q 60 sin− 45 sin 2 sin− 0  2  576β " r 4π ! −  aβ  −   aβ  4π + Q  4π + Q  q  q    1 Q 1 Q 1 aβ 9 sin 4 sin− sin 6 sin− + 45 sin 2 sin− + Q −   aβ  −   aβ  ( r 4π !)  4π + Q   4π + Q  q  q   1 aβ  1  aβ +9 sin 4 sin− + Q + sin 6 sin− + Q , (75) ( r 4π !) ( r 4π !) # 2 n where a = (1+ a1t + a2t + ... + ant ). Hence, the entropy is

1 Q 3 45aQ cos 2 sin− aβ 1 aβ 15aQ 4 +Q 15a   π  S0 = π + Q 3/2 + 3/2 + 576 4π " aβ Q2 aβ p Q2 2π 1 aβ Q aβ + Q   2π 4π + Q 1 aβ 4π 4π + Q 1 aβ 4π 4π − 4π +Q − 4π +Q − −   q   q q q 1 Q 1 Q 1 aβ 9aQ cos 4 sin− 3aQ cos 6 sin− 45a cos 2 sin− + Q aβ +Q aβ +Q 4π   4π    4π     + 3/2 + 3/2 + q aβ p Q2 aβ p Q2 4π 1 aβ Q aβ + Q 2π 4π + Q 1 aβ 4π 4π + Q 1 aβ 4π 4π − 4π +Q − 4π +Q − −   q   q q q 1 aβ 1 aβ 9a cos 4 sin− 4π + Q 3a cos 6 sin− 4π + Q +  q  +  q  2π 1 aβ Q aβ + Q 4π 1 aβ Q aβ + Q # − 4π − 4π − 4π − 4π q q q q 2 1 aβ 1 aβ 1 Q 1 Q + a + Q 60 sin− + Q 60 sin− 45 sin 2 sin− 768 4π " r 4π ! −  aβ  −   aβ    4π + Q  4π + Q  q  q    1 Q 1 aβ 1 Q 9 sin 4 sin− + 45 sin 2 sin− + Q sin 6 sin− −   aβ  ( r 4π !) −   aβ   4π + Q   4π + Q  q  q   1 aβ 1 aβ   +9 sin 4 sin− + Q + sin 6 sin− + Q ( r 4π !) ( r 4π !) # 3 π aβ + Q 4π 1 aβ 1 Q 1 Q 60 sin− + Q 60 sin− 45 sin 2 sin−   − 288β " r 4π ! −  aβ  −   aβ  4π + Q  4π + Q  q  q    1 Q 1 Q 1 aβ 9 sin 4 sin− sin 6 sin− + 45 sin 2 sin− + Q −   aβ  −   aβ  ( r 4π !)  4π + Q   4π + Q  q  q   1 aβ 1 aβ +9 sin 4 sin− + Q + sin 6 sin− + Q (76) ( r 4π !) ( r 4π !) #

The relation between S0, SBH and t looks like

1 Q 45aQ cos 2 sin− 3 1 15aQ √a√SBH +Q S0 = π a SBH + Q 2 +  2   576 Q 3/2 Q 3/2 "2π 1 a√SBH + Q 4π 1 a√SBH + Q  p  − a√SBH +Q − a√SBH +Q q  q 

16 1 Q 1 Q 9aQ cos 4 sin− 3aQ cos 6 sin− 15a √a√SBH +Q √a√SBH +Q + +    +    Q2 3/2 Q2 3/2 2π 1 a√SBH Q a√SBH + Q 2π 1 a√SBH + Q 4π 1 a√SBH + Q − − − a√SBH +Q − a√SBH +Q p p q q 1  1  45a cos 2 sin− a√SBH + Q 9a cos 4 sin− a√SBH + Q + + 4π 1 na√S pQ a√S +oQ 2π 1 na√S pQ a√S +oQ − BH − BH − BH − BH p 1 p p p 3a cos 6 sin− a√SBH + Q 2 1 1 + + a a S + Q 60 sin− a S + Q n  o BH BH 4π 1 a√SBH pQ a√SBH + Q# 768 " − −  p  q p  p p 1 Q 1 Q 1 Q 60 sin− 45 sin 2 sin− 9 sin 4 sin− − a√SBH + Q! − ( a√SBH + Q!) − ( a√SBH + Q!) p p p 1 Q 1 1 sin 6 sin− + 45 sin 2 sin− a SBH + Q + 9 sin 4 sin− a SBH + Q − ( a√SBH + Q!)  q p   q p  p 3 1 1 1 + sin 6 sin− a SBH + Q π a SBH + Q 60 sin− a SBH + Q # − 288 "  q p   p  q p 

1 Q 1 Q 1 Q 60 sin− 45 sin 2 sin− 9 sin 4 sin− − a√SBH + Q! − ( a√SBH + Q!) − ( a√SBH + Q!) p p p 1 Q 1 1 sin 6 sin− + 45 sin 2 sin− a SBH + Q + 9 sin 4 sin− a SBH + Q − ( a√SBH + Q!)  q p   q p  p 1 + sin 6 sin− a SBH + Q (77) #  q p  Case:I If a = 0 for i> 1 then (1 + a t + a t2 + ... + a tn) (1 + a t) i 1 2 n ≈ 1 Fig.7a

Relation between S0, SBH and t

Case:II If a = 0 for i> 2 then (1 + a t + a t2 + ... + a tn) (1 + a t + a t2) i 1 2 n ≈ 1 2

17 Fig.7b

Relation between S0, SBH and t

Case:III If a = 0 for i> 3 then (1 + a t + a t2 + ... + a tn) (1 + a t + a t2 + a t3) i 1 2 n ≈ 1 2 3 Fig.7c

Relation between S0, SBH and t

Relation between S , S and t for Case-I, Case-II and Case-III are shown in the figures 7a 7c. 0 BH − 7 Conclusion

A black hole is usually formed from the collapse of a quantity of matter or radiation, both of which carry entropy. However, the black hole’s interior and contents are veiled to an exterior observer. Thus a thermodynamic description of the collapse from that observer’s view point cannot be based on the entropy of that matter or radiation because these are unobservable. A stationary black hole is parametrized by just a few numbers (Ruffini and Wheeler 1971): its mass, electric charge and angular momentum (and magnetic monopole charge, except its actual existence in nature has not been demonstrated yet). For any specific choice of these parameters one can imagine many scenarios for the black hole’s formation. Thus there are many possible internal states corresponding to that black hole. In thermodynamics one meets a similar situation: many internal microstates of a system are all compatible with the one observed (macro) state. Thermodynamic entropy quantifies the said multiplicity. Thus by analogy one needs to associate entropy with a black hole .

18 In this paper, at first we calculated the entropy for the Schwarzchild black hole contained by the CR volume for massless modes. The approach are considered as statistical one. We first introduced the integrand of the expression of interior volume as an effective metric. Then easily identified the energy of the modes. To handle the situation, one can use the method of constraint analysis, since the canonical Hamiltonian vanishes. The exact calculation of the Hamiltonian is itself new. By using the Gibb’s free energy we can calculate the entropy. The effect is very interesting. As the mass of the black hole increasing as time increases, the average density of the universe is very small. So, the increment of the mass of the black hole is very slow. So we consider the mass of the black hole is 2 n 2 n M(1 + a1t + a2t + ... + ant ) after a time t. Then we assumed r varies from 2M(1 + a1t + a2t + ... + ant ) to zero. Finally we found that the entropy is monotonically increasing as time increases but the increment of entropy depends on time which is very small. Also we can see that the entropy on the horizon is greater than the entropy inside the black hole i.e. SBH > S. The entropy inside the black hole is proportional to the time(t) i.e. if we increase the time then the entropy(S) inside the black hole increase. Hope this paper will give a imagination knowledge about entropy.

Acknowledgement: This research is Supported by the project grant of Goverment of West Bengal, Department of Higher Education, Science and Technology and Biotecnology (File no:- ST/P/S&T/16G 19/2017). SD thanks West Bengal State Govt. for awarding Non-NET Fellowship. RB thanks Inter University Center− for Astronomy and Astrophysics(IUCAA), Pune, India for Visiting Associateship.

References

[1] Abbott, B. P. et al. :- “Observation of Gravitational Waves from a Binary Black Hole Merger”, Phys. Rev. Lett. 116, 061102 (2016); arXiv: 1602.03837v1[gr-qc]. [2] Hawking, S. W. :- “Particle creation by black holes”, Commun. Math. Phys. 43, 199-220 (1975); Preprint: Not available. [3] Hawking, S. W. :- “Breakdown of Predictability in gravitational collapse”, Phys. Rev. D 14, 2460 (1976); Preprint: Not available. [4] Bardeen, J. M., Carter, B. and Hawking, S. W. :- “The Four Laws of Black Holes Mechanics ”, Commun. Math. Phys. 31, 161-170 (1973); Preprint: Not available. [5] Bekenstein, J. D. :- “Black Holes and entropy”, Phys. Rev. D 7, 2333 (1973); Preprint: Not available. [6] Wald, R. M. :- “The thermodynamics of black holes”, Living Rev. Relativity 4, 6 (2001); arXiv: 9912119[gr-qc]. [7] Wald, R. M. :- “Black hole entropy is the Noether charge”, Phy. Rev. Relativity D 48, 3427-3431 (1993); arXiv: 9307038v1[gr-qc]. [8] Carroll, S., Jahnson, M. C. and Randall, L. :- “Extremal limits and black hole entropy”, JHEP 0911, 109 (2009); arXiv:0901.0931[hep-th]. [9] Bengtsson, I., Jakobson, E. :- “Black holes: Their large interiors”, Mod. Phys. Lett. A 30, 21 1550103 (2015); arXiv : 1502.01907[gr-qc]. [10] Zhang, B. :- “Entropy in the interior of a black hole and thermodynamics”, Phys Rev. D 92, 081501(R) (2015); arXiv: 1510.02182[gr-qc]. [11] Lange, F., Lenari, Z., Rosch, A. :- “Time-dependent generalized Gibbs ensembles in open quantum systems ”, Phys. Rev. B 97, 165138 (2018); arXiv:1801.07646 [cond-mat.str-el]. [12] Bak, D., Gutperle, M., Karch, M. :- “Time dependent black holes and thermal equilibration”, JHEP 0712 034, (2007); arXiv:0708.3691 [hep-th]. [13] Majhi, B. R., Samanta, S. :- “Entropy corresponding to the interior of a Schwarzschild black hole” Physics Letters B770, 314(2017). [14] Christodoulou, M., Rovelli, C. :- “How big is a black hole?”, Phys Rev. D 91, 064046 (2015); arXiv : 1411.2854[gr-qc].

19 [15] Hanson, A., Reggc, T., Teitelboim, C. :- “Constrained Hamiltonian System”, Academia Nazionale Dei lincei, (1976); Preprint: Not available. [16] Fulop, G., Gitman, D. M., Tyutin, I. V. :- “Reparametrization Invariance as Gauge Symmetry”, Int. J. Theor. Phys. 38, 1941-1968 (1999); arXiv :9805040[hep-th]. [17] Dirac, P. A. M. :- “Lectures on Quantum Mechanics”, Yeshiva University Press, New York (1964); Preprint: Not available. [18] Kawai, H., Yokokura, Y.:- “Interior of Black Holes and Information Recovery”, Phys. Rev. D 93, 044011 (2016); arXiv: 1509.08472[hep-th]. [19] Majhi, B. R., Samanta, S. :- “Entropy corresponding to the interior of a Schwarzschild black hole”, Phys. Lett. B 770, 314-318 (2017); arXiv: 1703.00142[gr-qc]. [20] Planck, M. :- “Zur Dynamik bewegter Systeme”, Ann. Physik 26, 1 (1908); Preprint: Not available. [21] Przanowski, M., Tosiek, J. :- “Notes on thermodynamics in special relativity”, Phys. Scr. 84, 055008 2011; arXiv: 1010.5701[cond-mat.stat-mech]. [22] Ott, H. :- “Lorentz-Transformation der Wrme und der Temperatur”, Z. Physik 175, 70 (1963); Preprint: Not available. [23] Arzelies, H. :- “Transformation relativiste de la temperature et de quelques autres grandeurs thermody- namiques”, Nuovo Cim. 35, 792 (1965); Preprint: Not available. [24] Mller, C. :- “Matematisk-fysiske meddelelser”, Det. Kong. Danske Videnskab. Selskab. Mat.-Fys. Medd. 36,1 (1967); Preprint: Not available. [25] Mller, C. :- “ Thermodynamics in the Special and the General Theory of Relativity”, Academic Press, 202221 1968; Preprint: Not available. [26] Mller, C. :- “The Theory of Relativity”, Oxford University Press, Oxford, 2nd edition (1972); Preprint: Not available. [27] Einstein, A. :- “ On the Relativity Principle and the conclusion drawn from it”, Jahrbuch der Radioaktivitat und Elektronik 4, 411 (1907); Preprint: Not available. [28] Mosengeil, K. V. :- Ann. Physik 22, 876 (1907); Preprint: Not available. [29] Pauli, W. :- Relativitatstheorie, Teubner, Berlin, (1921); Preprint: Not available. [30] Laue, M. V. :- Die Relativitatstheorie, Friedrich Vieweg und Sohn, Braunschweig, Germany, (1921); Preprint: Not available. [31] Landsberg, P. T. :- Nature 212, 571 (1966); Preprint: Not available. [32] Landsberg, P. T. :- Nature 214, 903 (1967); Preprint: Not available. [33] Landsberg, P. T., John, K. A. :- Ann. Phys. 56, 299 (1970); Preprint: Not available. [34] Landsberg, P. T. :- Europ. J. Phys. 2, 203 (1981); Preprint: Not available. [35] Hawking, S. W. :- “Black hole explosion?”, Nature 248, 30-31 (1974); Preprint: Not available. [36] Hawking, S. W. :- “The Large Scale Structure of Space-Time”, New York: Cambridge University Press. ISBN 0-521-09906-4, Ellis, G.F.R. (1973); Preprint: Not available. [37] HAwking, S., W. :- “The Nature of Space and Time”, Princeton University Press, 160 (1994); ariXiv: 9409195[hep-th]. [38] Castle, J., Emmenish, W., Henkes, R., Miller, R., and Rayne, J. :- “Science by Degrees: Temperature from Zero to Zero”, Westinghouse Search Book Series (Walker and Company, New York), (1965); Preprint: Not available. [39] Garrod, C. :- “ and Thermodynamics”, Oxford U. Press, (1995); Preprint: Not available. [40] Pathria, R. K. :- Proc. Nat. Inst. Sci. India 21 331 (1955); Preprint: Not available. [41] Pathria, R. K. :- Proc. Nat. Inst. Sci. India 23 168 (1957); Preprint: Not available. [42] Synge, J. L. :- The Relativistic Gas, North-Holland, Amsterdam, (1957); Preprint: Not available.

20 [43] Haar, D. T., Wegeland, H. :- “Thermodynamics and statistical mechanics in the special theory of relativity”, Phys. Rep. 1 31 (1971); Preprint: Not available. [44] Groot, S. R. D., Leeuwen, W. W. V., Weert, C. G. V. :- “Relativistic Kinetic Theory: Principles and Appli- cations”, North-Holland, Amsterdam, (1980); Preprint: Not available. [45] Cubero, D., Pascual, J. C., Dunkel, J., Talkner, P., Hanggi, P. :- “Thermal Equilibrium and Statistical Ther- mometers in Special Relativity”, Phys. Rev. Lett. 99 170601 (2007); Preprint: Not available. [46] Dunkel, J., Hnggi, P., Hilbert, S. :- “Nonlocal observables and lightcone-averaging in relativistic thermodynam- ics”, Nature Physics 5, 741 2009; arXiv: 0902.4651v2[cond-mat.stat-mech]. [47] Rasinariu, C. :- “Numeric experiments in relativistic thermodynamics: a moving system appears cooler”; arXiv: 0804.3836. [48] Montakhab, A., Ghodrat, M., Barati, M. :- “Statistical thermodynamics of a two-dimensional relativistic gas”, Phys. Rev. E 79 031124 (2009); arXiv: 0809.1517v2[cond-mat.stat-mech]. [49] Ghodrat, M., Montakhab, A. :- “Time parametrization and stationary distributions in a relativistic gas”, Phys. Rev. E 82 011110 (2010); arXiv: 0908.3753[cond-mat.stat-mech]. [50] Peano, F., Marti, M., Silva, L. O., Coppa, G. :- “Statistical kinetic treatment of relativistic binary collisions”, Phys. Rev. E 79 025701 (2009); arXiv: 0902.1762v1[physics.plasm-ph]. [51] Landau, L. D., Lifshitz, E. M. :- Statistical Physics, Pergamon Press, Oxford, (1980); Preprint: Not available. [52] Carrera, B. L., Granados, V., Parga, G. A. D. :- “An alternative deduction of relativistic transformations in thermodynamics”, Rev. Mex. Fs. 54 1519 (2008); Preprint: Not available. [53] Bebronne, M. V. and Tinyakov, P. G. :- “Black hole solutions in massive gravity”, JHEP 0904, 100 (2009); arXiv: 0902.3899[gr-qc].

21