<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">JMP</journal-id><journal-title-group><journal-title>Journal of Modern Physics</journal-title></journal-title-group><issn pub-type="epub">2153-1196</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/jmp.2012.329166</article-id><article-id pub-id-type="publisher-id">JMP-23118</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Time Evolution of Horizons
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>rundhati</surname><given-names>Dasgupta</given-names></name><xref ref-type="aff" rid="aff1"><sub>1</sub></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff1"><label>1</label><addr-line>Department of Physics and Astronomy, University of Lethbridge, Lethbridge, Canada</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>arundhati.dasgupta@uleth.ca</email></corresp></author-notes><pub-date pub-type="epub"><day>28</day><month>09</month><year>2012</year></pub-date><volume>03</volume><issue>09</issue><fpage>1289</fpage><lpage>1297</lpage><history><date date-type="received"><day>June</day>	<month>24,</month>	<year>2012</year></date><date date-type="rev-recd"><day>July</day>	<month>25,</month>	<year>2012</year>	</date><date date-type="accepted"><day>August</day>	<month>2,</month>	<year>2012</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  Finding the origin of Hawking radiation has been a puzzle to researchers. Using a loop quantum gravity description of a black hole slice, a density matrix is defined using coherent states for space-times with apparent horizons. Evolving the density matrix using a semi-classical Hamiltonian in the frame of an observer outside the horizon gives the origin of Hawking radiation.
 
</p></abstract><kwd-group><kwd>Black Hole Physics; Quantum Gravity; Hawking Radiation; Loop Quantum Gravity; Coherent States</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>A new theory is expected to take over at Planck distances as “quantum effects” of gravity start dominating. One of the promising approaches to the theory of quantum gravity is the theory of Loop Quantum Gravity (LQG), which is by formulation non-perturbative and background independent [1-3]. LQG has a well defined kinematical Hilbert space, and though the Hamiltonian constraint remains unsolved, the theory allows for a semiclassical sector of the theory. This includes “coherent states” [4,5] which are peaked at classical phase space elements. Using these as a starting point, I defined in a series of papers [6-8] coherent states for the Schwarzschild spacetime, and derived an origin of entropy using quantum mechanical definition of entropy from density matrices. The exact entropy is a function of the graph used to obtain the LQG phase space variables [<xref ref-type="bibr" rid="scirp.23118-ref9">9</xref>]. The zeroeth order term is proportional to the area of the horizon signifying a universality of the Bekenstein-Hawking entropy. The proportionality constant and the correction terms bring out the details of the graph [<xref ref-type="bibr" rid="scirp.23118-ref8">8</xref>].</p><p>In this paper we take this new way of finding the origin of entropy a step further by evolving the spatial slice in time [<xref ref-type="bibr" rid="scirp.23118-ref10">10</xref>], and observing the evolution of the density matrix in the process. This state as of now does not satisfy the Hamiltonian constraint, but one is allowed to take an arbitrary initial state, or a wavepacket with appropriate properties, representing a macroscopic configuration. The evolution discussed in this paper is semiclassical, i.e. no attempt is made to use the full Hamiltonian.</p><p>The quasilocal energy (QLE) of an outside observer, defined in [<xref ref-type="bibr" rid="scirp.23118-ref11">11</xref>] is used as the Hamiltonian to evolve the system. As the time clicks in the observers clock, the Hamiltonian evolves the coherent state such that the area of the horizon remains the same as predicted by classical physics. However, classically forbidden regions become accessible quantum mechanically, and vertices of the graph hidden behind the horizon in one slice emerge outside the horizon in the next slice. This gives a net change in area, and the mass deficit is emitted from the black hole. This evolution is not unitary, and the quasi-local energy which is used to evolve the slice is not mapped to a Hermitian operator. When matter is coupled to the gravitational system, a net flux emerges causing a decay of the horizon.</p><p>In Section 2 we introduce the formalism by describing the coherent state, the black hole time slice, the apparent horizon equation, and the density matrix. Section 3 describes the time evolution of the system and gives a derivation of the change in entropy. In Section 4 we give a description of a matter current emergent from behind the horizon. Finally in the concluding section we include a discussion about the implications of the nonHermitian evolution.</p></sec><sec id="s2"><title>2. The Coherent State in LQG</title><p>For gravity, finding the canonical variables which describe the physical phase space is an odd task as there is no unique time. Nevertheless a fiducial time coordinate can be chosen, which breaks the manifest diffeomorphism invariance, restored in the Hilbert space of states by imposing constraints.</p><p>The constant time slices are described by the intrinsic metric <img src="23-7500852\7422fd29-974e-4f6e-ab29-d24d2b9a7ce4.jpg" /> and the extrinsic curvature <img src="23-7500852\7b7c7862-44cf-4f67-bf77-5357a6f8e5d7.jpg" /> (a,b = 1,2,3). The theory can be formulated in terms of the square root of the metric, the triads <img src="23-7500852\a78842fb-6f2a-479f-843c-81f7bbdc0293.jpg" /> defined thus:</p><disp-formula id="scirp.23118-formula70891"><label>(1)</label><graphic position="anchor" xlink:href="23-7500852\57d0e0c8-f416-472d-9b87-1bc9b73d66ea.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="23-7500852\a7d43f42-dd5f-4a97-847b-fd3dbe7f89fb.jpg" /> represents the internal index for the rotation group SO(3) of the tangent space and<img src="23-7500852\85c03032-271a-46d1-ba46-d0c45c90d134.jpg" />. The internal group is taken to be SU(2), as this is locally isomorphic to SO(3). The theory is then defined in terms of the “spin connection” <img src="23-7500852\f1fc650d-a37c-46af-9571-50681edb86e3.jpg" />and the triads. However, a redefinition of the variables in terms of tangent space densitised triads <img src="23-7500852\a3a9ced2-00d5-411e-81d4-3cde45e3a86c.jpg" /> and a corresponding gauge connection <img src="23-7500852\635711fb-4836-44ff-8930-34860069d5d9.jpg" /> where I represents the SU(2) index simplifies the quantisation considerably.</p><disp-formula id="scirp.23118-formula70892"><label>(2)</label><graphic position="anchor" xlink:href="23-7500852\ac068906-a6c5-4c60-84bc-f1ab55fdea38.jpg"  xlink:type="simple"/></disp-formula><p>(<img src="23-7500852\bee00e3e-ac4b-47f7-97ae-de48fb45c0d3.jpg" />are the usual triads, <img src="23-7500852\706904a0-a718-4ab5-993e-72e06af25ece.jpg" />is the extrinsic curvature, <img src="23-7500852\60520239-1884-494c-a346-3933a7aa1087.jpg" />the associated spin connection, <img src="23-7500852\f656ce80-3cb5-4a06-8dbc-52b3bcd8f36e.jpg" />the one parameter ambiguity which remains named as the Immirzi parameter). The quantisation of the Poisson algebra of these variables is done by smearing the connection along one dimensional edges <img src="23-7500852\bc56a4f3-afc6-46b2-aab6-28f729fcd582.jpg" /> of length <img src="23-7500852\bf648ae9-6fb3-41e1-9e32-d0fd07c1a1c6.jpg" /> of a graph <img src="23-7500852\8cb04f83-aac5-4c5a-bcd8-0e04c677bf60.jpg" /> to get holonomies<img src="23-7500852\b11502c5-5288-4a86-8f65-00cdba178233.jpg" />. The triads are smeared in a set of 2-surface decomposition of the three dimensional spatial slice to get the corresponding momentum<img src="23-7500852\e1f6344d-1548-4e16-b34d-ac5c54269117.jpg" />. The algebra is then represented in a kinematic “Hilbert space”, in which the physical constraints have been “formally” realised [<xref ref-type="bibr" rid="scirp.23118-ref12">12</xref>]. Once the phase space variables have been identified, one can write a coherent state for these [<xref ref-type="bibr" rid="scirp.23118-ref4">4</xref>] i.e. minimum uncertainty states peaked at classical values of<img src="23-7500852\13371314-d2c8-452d-9691-b6f285b21df2.jpg" />. In analogy with the harmonic oscillator coherent states, where the coherent state is a function of the complexified phase space element<img src="23-7500852\b3d838f1-8ec0-4cba-9cbd-b8a1f5d73fec.jpg" />, the SU(2) coherent states are peaked at the complexified phase space element<img src="23-7500852\5227f0d0-4067-4cca-ab69-3c48434a4b7e.jpg" />. These <img src="23-7500852\2ea385a8-1d3f-4006-adf2-5f971d940a47.jpg" /> are thus elements in the complexification of SU(2) as <img src="23-7500852\2791f29f-e41f-442c-b8a0-61574ff77361.jpg" /> (<img src="23-7500852\2c07b608-672f-4404-bc60-1304dedb28a2.jpg" />being the generator matrices of SU(2)) is a Hermitian matrix and <img src="23-7500852\f499067d-f51c-486d-a1c0-1d8e9e3d6a35.jpg" /> is the unitary SU(2) matrix. Whether these are physical coherent states, or have appropriate behavior under the action of the constraints has to be examined carefully [<xref ref-type="bibr" rid="scirp.23118-ref13">13</xref>]. The coherent state in the momentum representation for one edge is defined to be</p><disp-formula id="scirp.23118-formula70893"><label>(3)</label><graphic position="anchor" xlink:href="23-7500852\89c117db-05c5-40e2-994d-333789c9a87d.jpg"  xlink:type="simple"/></disp-formula><p>In the above <img src="23-7500852\559c0c74-4e22-42a8-bb98-c050807b3f72.jpg" /> is a complexified classical phase space element<img src="23-7500852\13ff192b-6432-4390-a4b4-9269ffe84af0.jpg" />, (the <img src="23-7500852\fb3dbd6f-afa9-44d6-9721-55ca9e85cb67.jpg" /> and the <img src="23-7500852\69f0de0b-bd91-49bf-b132-2dad55e51cb9.jpg" /> represent classical momenta and holonomy obtained by embedding the edge in the classical metric). The <img src="23-7500852\03c974c7-a5ac-446a-8348-d782790218c6.jpg" /></p><p>are the usual basis spin network states given by<img src="23-7500852\6f63d2aa-4b83-45b3-be3e-72a781d897aa.jpg" />which is the jth representation of the SU(2) element<img src="23-7500852\72b3ef5b-d84d-43e6-b1b1-0da6d043baa6.jpg" />. Similarly, <img src="23-7500852\f68dda39-26d7-4b86-a02a-d8a1f1252b23.jpg" />dimensional representations of the <img src="23-7500852\5a49bf05-d581-46ac-a84c-b623380e24c3.jpg" /> matrix <img src="23-7500852\6c30f529-201f-4040-b91b-7fce5e35e3ac.jpg" /> are denoted as<img src="23-7500852\6b432488-141a-4158-9683-076029505c4d.jpg" />. The j is the quantum number of the SU(2) Casimir operator in that representation, and <img src="23-7500852\b2b1b698-dc2c-48bd-9c82-40337ed4bb81.jpg" /> represent azimuthal quantum numbers which run from<img src="23-7500852\beedf4e5-b829-4f18-ac30-dc60d1285291.jpg" />. The coherent state is precisely peaked with maximum probability at the <img src="23-7500852\8aa7b22a-0d2f-4215-8b85-9617fb4b0c1f.jpg" /> for the variable <img src="23-7500852\7ea5fc46-d4fc-4406-9923-ed486ad6b5a0.jpg" /> as well as the classical momentum <img src="23-7500852\b37d7fc6-021d-4b47-95e6-ef11ce38e3d3.jpg" /> for the variable<img src="23-7500852\ee96c676-59a2-4a48-9a27-55591e2fa596.jpg" />. The fluctuations about the classical value are controlled by the parameter t (the semiclassicality parameter). This parameter is given by <img src="23-7500852\2234ac28-d088-48fe-ad6e-afa7de5ebb4a.jpg" /> where <img src="23-7500852\5f188178-508f-4d51-974c-5d5447503172.jpg" /> is Planck’s constant and a a dimensional constant which characterises the system. The coherent state for an entire slice can be obtained by taking the tensor product of the coherent state for each edge which form a graph<img src="23-7500852\7435adf3-6499-4d51-b374-de8f70c0d126.jpg" />,</p><disp-formula id="scirp.23118-formula70894"><label>(4)</label><graphic position="anchor" xlink:href="23-7500852\5c4b81f6-470f-4a75-a52c-5381982ac21b.jpg"  xlink:type="simple"/></disp-formula><p>In [<xref ref-type="bibr" rid="scirp.23118-ref7">7</xref>] the <img src="23-7500852\d2d0b8d1-7e44-45c0-8870-c94399f3f3fc.jpg" /> was evaluated for the Schwarzschild black hole by embedding a graph on a spatial slice with zero intrinsic curvature. The particular graph which was used had the edges along the coordinate lines of a sphere. This simplistic graph, was very useful in obtaining the description of the space-time in terms of discretised holonomy and momenta. A particularly interesting consequence of this was that the phase space variables were finite and well defined even at the singularity.</p><p>Given that the area of a surface in gravity is measured as the integral of the square root of the metric over the surface, the area operator can be written simply as <img src="23-7500852\0d1f364f-52fe-4ba8-8359-82cd290444ed.jpg" />. The expectation value of the area operator in the coherent state emerges as [<xref ref-type="bibr" rid="scirp.23118-ref9">9</xref>]</p><disp-formula id="scirp.23118-formula70895"><label>(5)</label><graphic position="anchor" xlink:href="23-7500852\08ba2452-3c90-438e-b0d1-e8d32807f955.jpg"  xlink:type="simple"/></disp-formula><p>Thus we are considering a semiclassical state, which is a state such that expectation values of operators are closest to their classical values. The information of the classical phase space variables are encoded in the complexified SU(2) elements labeled as<img src="23-7500852\58659058-6073-4f7a-bbf9-650f0c0c2566.jpg" />. The fluctuations over the classical values are controlled by the semiclassical parameter t.</p><p>The density matrix which describes the entire black hole slice is obtained as</p><disp-formula id="scirp.23118-formula70896"><label>(6)</label><graphic position="anchor" xlink:href="23-7500852\960e427a-20f7-474a-b6fa-e7a65b0fc3c3.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="23-7500852\205affa6-0404-4f8b-ba12-e7e1ba0152f0.jpg" /> is the coherent state wavefunction for the entire slice, a tensor product of coherent state for each edge.</p><sec id="s2_1"><title>2.1. Apparent Horizons</title><p>We concentrate on the coherent state near the apparent horizon contained in the spatial slice. We find that motivated from the apparent horizon equation the graph across the horizon can be taken to be populated by radial edges, linking vertices outside and inside the horizon. One then traces over the coherent state within the horizon. Initially we take a particular time slicing of the black hole, which has the spatial slices with zero intrinsic curvature [<xref ref-type="bibr" rid="scirp.23118-ref7">7</xref>]. One such metric which has the time slices as flat is the Lemaitre metric</p><disp-formula id="scirp.23118-formula70897"><label>(7)</label><graphic position="anchor" xlink:href="23-7500852\0f906bf8-1a04-4e79-b6dc-fb6b7f112bbe.jpg"  xlink:type="simple"/></disp-formula><p>The<img src="23-7500852\ac0ef48e-4d8d-44f7-84c6-5aaa7cf3b7d7.jpg" />, (in units of c = 1) and in the <img src="23-7500852\fd371ca5-b6da-4ac9-8c37-5b0f38bdf2d5.jpg" /> constant slices one can define the induced metric in terms of a “r” coordinate defined as</p><p><img src="23-7500852\fe70981a-23d2-40f5-909c-1951da31e203.jpg" /></p><p>(<img src="23-7500852\9ef583ed-74ad-411d-b532-60bf722f698f.jpg" />) on the slice. One gets the metric of the three slice to be</p><disp-formula id="scirp.23118-formula70898"><label>(8)</label><graphic position="anchor" xlink:href="23-7500852\b9da7589-7fa9-4da8-a3d3-e77073223be0.jpg"  xlink:type="simple"/></disp-formula><p>The entire curvature of the space-time metric is contained in the extrinsic curvature or <img src="23-7500852\2aa690ee-ff98-42d3-ba7c-16b8d74ed401.jpg" /> tensor of the <img src="23-7500852\8736b40e-f8ed-4202-9965-97ad39e6ece4.jpg" /> constant slices. Now if there exists an apparent horizon somewhere in the above spatial slice, then that is located as a solution to the equation</p><disp-formula id="scirp.23118-formula70899"><label>(9)</label><graphic position="anchor" xlink:href="23-7500852\1f3d4283-b304-4f6f-977f-65262476c604.jpg"  xlink:type="simple"/></disp-formula><p>where<img src="23-7500852\2103d6df-f9bf-4f0c-9130-ec71861c45a0.jpg" />, (<img src="23-7500852\4c4fbf4f-4cf2-4662-a171-a8ed770546c3.jpg" />denote the spatial indices) is the normal to the horizon, <img src="23-7500852\6a031da5-2a49-4d92-b7fd-b8e3476ac723.jpg" />the extrinsic curvature in the induced coordinates of the slice, and <img src="23-7500852\9f8b1dfe-9347-4e01-a9eb-ba34c5efac1d.jpg" /> the trace of the extrinsic curvature. If the horizon is chosen to be the 2-sphere, then in the coordinates of (8), <img src="23-7500852\fee155ac-47be-4748-9b8f-06cea1a1b70c.jpg" />, the apparent horizon equation as a function of the metric reduces to:</p><disp-formula id="scirp.23118-formula70900"><label>(10)</label><graphic position="anchor" xlink:href="23-7500852\ada3900e-e661-4f63-8742-f0beaa40190b.jpg"  xlink:type="simple"/></disp-formula><p>Note that the first term of the equation disappears trivially as <img src="23-7500852\961a4b9c-57a6-4f92-a723-b37c2f0ff17e.jpg" /> for any point in the spatial slice. Even at the operator level the <img src="23-7500852\cb1a16f7-9586-4079-9831-d013d912d198.jpg" /> can be set to the identity operator in the first approximation, as <img src="23-7500852\cb64c1ae-a7fd-425c-a8cd-cc1282f1c1c8.jpg" /> (<img src="23-7500852\d1996c1a-d940-4a28-ad4a-bc89f5915a53.jpg" />being the volume operator) upto normalisations, and in the spherically symmetric metric <img src="23-7500852\f9adb9cd-c22d-4f90-bf12-0cb4c1e5dce7.jpg" /> (upto discretisation constants). Thus the operators in the numerator and denominator cancel and the normalisation conspire, leaving<img src="23-7500852\163539bd-27ac-426f-9f4e-f53cddbb5a10.jpg" />. To understand the rest of the equation in terms of the holonomy and momentum variables of LQG, which are classically measured in the same metric as (8), we use the following regularisation</p><disp-formula id="scirp.23118-formula70901"><label>(11)</label><graphic position="anchor" xlink:href="23-7500852\d8ff6164-1bac-4997-895c-8dd78cb8d02c.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.23118-formula70902"><label>(12)</label><graphic position="anchor" xlink:href="23-7500852\0ed3b182-85dd-4602-9d32-75e64abf5cc9.jpg"  xlink:type="simple"/></disp-formula><p>(N is a constant, a function of the edge lengths and the area bits of the discretisation) and V is the volume operator.</p><disp-formula id="scirp.23118-formula70903"><label>(13)</label><graphic position="anchor" xlink:href="23-7500852\3c0926ff-f736-44d1-8ce7-330bc04fb52f.jpg"  xlink:type="simple"/></disp-formula><p>Here <img src="23-7500852\04cd8a0a-4584-430e-9d77-ac1b1e0c159f.jpg" /> has been used as a parameter to identify the <img src="23-7500852\a4ee68d3-ea58-4133-9459-4e84241d86ef.jpg" /> operator, and this is mainly a trick. In the continuum limit</p><disp-formula id="scirp.23118-formula70904"><label>(14)</label><graphic position="anchor" xlink:href="23-7500852\44dcdab2-80e6-47da-aafc-f8980dc6e070.jpg"  xlink:type="simple"/></disp-formula><p>As the gauge connection is a function of the Immirzi parameter due to (2), the expectation value of this operator in a coherent state will be a function of the Immirzi parameter. By taking the derivative wrt to the Immirzi parameter we are giving the same status to the parameter as is given to “dimension” in a dimensional regularisation of Feynman diagrams. We let the parameter vary by an infinitesimal amount from its value in the particular quantisation sector, take the derivative, and put its original value in the final answer for the <img src="23-7500852\0160a88e-a369-440f-bf0d-7e335480bd88.jpg" /> operator. The Formula (13) is facilitated by the fact that the dependence of <img src="23-7500852\7555c6f6-acc8-449a-bac2-fd0c459d01b6.jpg" /> on the <img src="23-7500852\b0f47d02-6d3c-4f6e-9955-8ab06263297d.jpg" /> is linear. One way to check whether this gives the proper answer is to take a solved quantum mechanical system and use a similar method there. The most useful example is the Harmonic Oscillator Hamiltonian, which can be written as</p><disp-formula id="scirp.23118-formula70905"><label>(15)</label><graphic position="anchor" xlink:href="23-7500852\af2bbadd-a122-44e6-949c-5d46827ff6d7.jpg"  xlink:type="simple"/></disp-formula><p>The ground state is a coherent state, so we take that as an example. We define the operator</p><disp-formula id="scirp.23118-formula70906"><label>(16)</label><graphic position="anchor" xlink:href="23-7500852\6e7140d4-6e3f-4656-b25e-504e10bf8cad.jpg"  xlink:type="simple"/></disp-formula><p>Thus</p><disp-formula id="scirp.23118-formula70907"><label>(17)</label><graphic position="anchor" xlink:href="23-7500852\cc7e5e9f-9ac1-434e-9b69-b058d5eff60a.jpg"  xlink:type="simple"/></disp-formula><p>The regularisation (13) is thus an allowed approximation.</p><p>The terms involving the Christoffel connections like <img src="23-7500852\9692ce77-59ef-48c1-875b-7395326d1b30.jpg" /> include derivatives in the regularised version, the derivatives appear as difference of triads across two vertices. Thus</p><disp-formula id="scirp.23118-formula70908"><label>(18)</label><graphic position="anchor" xlink:href="23-7500852\29963217-2db2-423f-a030-32476ee76069.jpg"  xlink:type="simple"/></disp-formula><p>As a result of this if we impose restrictions on the Christoffel connections and one of the vertices <img src="23-7500852\73e899c1-0ece-4fe6-8a89-1de7f1c87cd3.jpg" /> is within the horizon, whereas <img src="23-7500852\61f45bb9-7425-48c7-98df-d1d95a29b928.jpg" /> is outside the horizon, there will be correlations across the horizon.</p><p>If one evaluates the expectation value of the apparent horizon equation using the regularised variables in the coherent states, then one would obtain</p><disp-formula id="scirp.23118-formula70909"><label>(19)</label><graphic position="anchor" xlink:href="23-7500852\04e2e5dd-2700-4653-9c4f-462e6ac93d49.jpg"  xlink:type="simple"/></disp-formula><p>(<img src="23-7500852\8b78466f-7f74-4d68-b458-223b840066dc.jpg" />is a constant).</p></sec><sec id="s2_2"><title>2.2. Density Matrix</title><p>The density matrix is obtained as</p><disp-formula id="scirp.23118-formula70910"><label>(20)</label><graphic position="anchor" xlink:href="23-7500852\566bdbde-20f9-413c-a2cc-554573dda2fd.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="23-7500852\ec7004ab-05ed-4568-9d6d-3ebdbbff3ee5.jpg" /> is the coherent state wavefunction for the entire slice, a tensor product of coherent state for each edge.</p><p>But given this, we concentrate in a “local” region to see the behavior of the horizon</p><disp-formula id="scirp.23118-formula70911"><label>(21)</label><graphic position="anchor" xlink:href="23-7500852\4c6c93ba-0d4f-4241-813a-850017c084df.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="23-7500852\75ce5a5b-196b-48dc-a85f-1c2309e62256.jpg" /> covers a band of vertices surrounding the horizon one set on a sphere at radius <img src="23-7500852\2cf7eb3c-95c1-40e8-89fa-46bf4e927f51.jpg" /> and one set on a sphere at radius <img src="23-7500852\fd60f515-38d1-406c-aef3-67be74339584.jpg" /> within the horizon, as described in [<xref ref-type="bibr" rid="scirp.23118-ref9">9</xref>], and in the figure enclosed. This local density matrix and the correlations due to the apparent horizon equation (19) was used to derive entropy [<xref ref-type="bibr" rid="scirp.23118-ref6">6</xref>]. This entropy counts the number of ways to induce the horizon area using the spin networks, though the constraints have not been appropriately imposed as was obtained using a Chern-Simons theory in [<xref ref-type="bibr" rid="scirp.23118-ref14">14</xref>]. However, the entropy calculation using the coherent states provides a tracing mechanism, and a method to obtain correlations across the horizon which are gravitational in origin. We will henceforth deal with<img src="23-7500852\5e88b053-72b5-4ae5-9272-0f7598b53d92.jpg" />, but we will drop the local label for brevity.</p><p><img src="23-7500852\0d1b28d3-e5ef-4e8c-ab69-49e0cdeb88b5.jpg" /></p></sec></sec><sec id="s3"><title>3. Time Evolution</title><p>In physical systems, the Hamiltonian generates time evolution, but in General Theory of Relativity, the Hamiltonian is a constraint and generates diffeomorphisms in the time direction. So the question is, what is physical time, and if that exists, what would be the operator evolving the system in that direction? In case of space-times with time like Killing vectors, notion of time can be identified with the Killing direction, and a notion of “quasilocal energy” (QLE) defined using the same. The QLE then generates translations in the Killing time. In case of the Schwarzschild space-time, the QLE has been defined in [<xref ref-type="bibr" rid="scirp.23118-ref11">11</xref>]. We build the Hamiltonian which evolves the horizon from one time slice to the next by appropriately regularising the QLE. Note the “Killing time” and QLE are classical concepts, and thus regularising QLE gives us a “semiclassical” Hamiltonian.</p><sec id="s3_1"><title>3.1. Change in Entropy</title><p>Before we get into the analysis of what QLE evolution means, we take a simple system made up of two subsystems, and examine the consequences of a Hamiltonian evolution. Let the density matrix be defined for a system whose states are given in the tensor product Hilbert space <img src="23-7500852\5c4ed5af-0078-4227-ad29-b6227d355c37.jpg" /> and given by</p><disp-formula id="scirp.23118-formula70912"><label>(22)</label><graphic position="anchor" xlink:href="23-7500852\a50359f4-0e62-4802-bd9d-87c109ece04d.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="23-7500852\14cb19d5-a62d-48c2-a316-a3fdefaea8bb.jpg" /> is the basis in <img src="23-7500852\09da8d8c-83c7-4524-8d6a-3202eb0ebd17.jpg" /> and <img src="23-7500852\530af0a8-39e8-40bd-a093-ca529fe5d74b.jpg" /> is the basis in <img src="23-7500852\83e21066-4701-41b0-9bc6-e2708805a3be.jpg" /> and <img src="23-7500852\d564ff0e-8fe8-4636-bf77-abebb5fe6c59.jpg" /> are the non-factorisable coefficients of the wavefunction in this basis. Let us label the wavefunction at time <img src="23-7500852\b7e225f2-1bd7-401b-87d8-90f2049cacb7.jpg" /> to be given by the coefficients<img src="23-7500852\5397b6c0-da22-44b0-bc4d-2e34e1295345.jpg" />. The density matrix is</p><disp-formula id="scirp.23118-formula70913"><label>(23)</label><graphic position="anchor" xlink:href="23-7500852\390737d4-16d3-4f58-830c-68cfd75d3691.jpg"  xlink:type="simple"/></disp-formula><p>The reduced density matrix if one traces over <img src="23-7500852\21cba904-e2d5-45e8-8fc1-3e7dd65dfc63.jpg" /> is:</p><disp-formula id="scirp.23118-formula70914"><label>(24)</label><graphic position="anchor" xlink:href="23-7500852\55593a53-07cb-43c6-a769-d4f3c5e67157.jpg"  xlink:type="simple"/></disp-formula><p>We now evolve the system using a Hamiltonian which has the matrix elements<img src="23-7500852\4192d169-2228-43e8-b98f-36c0d1d1f567.jpg" />, we assume that the Hamiltonian does not factorise, that is there exists interaction terms between the two Hilbert spaces. The evolution equation is:</p><disp-formula id="scirp.23118-formula70915"><label>(25)</label><graphic position="anchor" xlink:href="23-7500852\284acaa2-43b6-481c-b517-4ee0d00fb428.jpg"  xlink:type="simple"/></disp-formula><p>which in this particular basis gives the density matrix elements at a infinitesimally nearby slice to be</p><disp-formula id="scirp.23118-formula70916"><label>(26)</label><graphic position="anchor" xlink:href="23-7500852\649eedcd-b9dd-4120-86fe-3147df639db1.jpg"  xlink:type="simple"/></disp-formula><p>Thus we evolve the “unreduced” density matrix and then trace over the <img src="23-7500852\bcf52228-5569-4be6-8ee5-536f34b849d7.jpg" /> in the evolved slice. The reduced density matrix in the evolved slice is:</p><disp-formula id="scirp.23118-formula70917"><label>(27)</label><graphic position="anchor" xlink:href="23-7500852\fe6b35df-2322-4e7a-b668-9a7565c27fd2.jpg"  xlink:type="simple"/></disp-formula><p>This gives:</p><disp-formula id="scirp.23118-formula70918"><label>(28)</label><graphic position="anchor" xlink:href="23-7500852\b6085a0c-1daf-41b9-9406-59eda5094677.jpg"  xlink:type="simple"/></disp-formula><p>where A represents the commutator. Clearly the entropy in the evolved slice evaluated as</p><p><img src="23-7500852\b180c960-5415-4c95-b4ee-a4ccef448134.jpg" /></p><p>can be found as</p><disp-formula id="scirp.23118-formula70919"><label>(29)</label><graphic position="anchor" xlink:href="23-7500852\30b7cae1-45d4-4498-b0ad-dac6a4b450c1.jpg"  xlink:type="simple"/></disp-formula><p>Given the definition of<img src="23-7500852\443c8bee-4245-4a10-a9dc-d1164348f06e.jpg" />, one gets</p><disp-formula id="scirp.23118-formula70920"><label>(30)</label><graphic position="anchor" xlink:href="23-7500852\13dadce4-25af-489c-9d65-14a37bf77bff.jpg"  xlink:type="simple"/></disp-formula><p>In case both the Hamiltonian and the density operator are Hermitian, one obtains</p><disp-formula id="scirp.23118-formula70921"><label>(31)</label><graphic position="anchor" xlink:href="23-7500852\8099d9f3-bf63-498a-b83d-9a9bfdd4c04d.jpg"  xlink:type="simple"/></disp-formula><p>This is clearly calculable, and gives the change in entropy<img src="23-7500852\de505a11-2018-4b34-b197-e59e3b927da4.jpg" />. The <img src="23-7500852\bac4805f-14dd-458c-ba61-9de4d63ba448.jpg" /> term yields corrections, and we ignore it in the first approximation.</p></sec><sec id="s3_2"><title>3.2. The Hamiltonian</title><p>To trace the origin of Horizon fluctuations, we must take an observer who is stationed outside the horizon, or in other words is not a freely falling observer. The quasilocal energy is defined using a “surface” integral of the extrinsic curvature with which the surface is embedded in three space. In our case, we take the bounding surface to be the horizon and the quasilocal energy is given by the surface term [11,15].</p><disp-formula id="scirp.23118-formula70922"><label>(32)</label><graphic position="anchor" xlink:href="23-7500852\6cb2040e-15aa-43c8-9ea6-28760bd93184.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="23-7500852\6149569c-42e5-492c-bc04-97b63e3dffca.jpg" /> is the extrinsic curvature with which the 2- surface, which in this case is the horizon <img src="23-7500852\78188ba9-4fcb-4d49-8125-0d52ac5568b7.jpg" /> is embedded in the spatial 3-slices, and <img src="23-7500852\d99c4e83-6a21-4076-8fa4-d3bbb46b5fe5.jpg" /> is the determinant of the two metric <img src="23-7500852\fc3853f9-9f26-45f3-9f78-10d8a5f8fb5b.jpg" /> defined on the 2-surface. This “quasilocal energy” is measured with reference to a background metric. Thus<img src="23-7500852\2955b621-e244-4ada-93ff-3dccf645ce52.jpg" />. We concentrate on the physics observed in an observer stationed at a r = constant sphere.</p><p>The metric in static <img src="23-7500852\211258e5-3d52-4b55-b73a-295c484fae9e.jpg" /> observer’s frame is</p><disp-formula id="scirp.23118-formula70923"><label>(33)</label><graphic position="anchor" xlink:href="23-7500852\fa076fa5-15ff-431a-a9e4-57cd7521d371.jpg"  xlink:type="simple"/></disp-formula><p>The <img src="23-7500852\e775e3b7-82c1-4560-baa0-1bbebc8b342d.jpg" /> where <img src="23-7500852\ab235ab0-cda3-4373-876a-ffaacd59eda6.jpg" /> is the Schwarzschild radius. If we take <img src="23-7500852\5ff0e70e-7c0b-40a6-a678-0a6c75dc5235.jpg" /> to be the space-like vector, normal to the 2-surface, then the extrinsic curvature is given by:</p><disp-formula id="scirp.23118-formula70924"><label>(34)</label><graphic position="anchor" xlink:href="23-7500852\8751427b-33a5-442d-aabb-a31be917680e.jpg"  xlink:type="simple"/></disp-formula><p>and the trace is obviously</p><disp-formula id="scirp.23118-formula70925"><label>(35)</label><graphic position="anchor" xlink:href="23-7500852\0262c13e-0462-4a71-8cb4-231f47b2aa1c.jpg"  xlink:type="simple"/></disp-formula><p>In the special slicing of the of the stationary observer the normal to the horizon 2-surface is given by <img src="23-7500852\cab5f91a-5222-4bc6-b66a-f4c8ebcbd96a.jpg" />. However, we built the coherent state on the Lemaitre slice. The Lemaitre and the Schwarzschild observer’s coordinates are related by the following coordinate transformations,</p><p><img src="23-7500852\60def415-b0ba-4ade-a109-e978e5cfa0c9.jpg" /></p><disp-formula id="scirp.23118-formula70926"><label>(36)</label><graphic position="anchor" xlink:href="23-7500852\1489df42-901d-4089-85b2-dd839724f061.jpg"  xlink:type="simple"/></disp-formula><p>The r = const cylinder of the Schwarzschild coordinate corresponds to <img src="23-7500852\cd98872c-63bc-4d3b-8b80-b972cbbd425d.jpg" /> of the Lemaitre coordinates, and for these<img src="23-7500852\dd4c88ee-1475-454e-a523-f20be96b616b.jpg" />. Thus unit translation in the t coordinate coincides with unit translation in the <img src="23-7500852\090b1df2-f9fa-428a-967c-c039c00ae09a.jpg" /> coordinate. Further, the intersection of the r = constant cylinder with a t = constant surface coincides with the intersection of r = constant and the <img src="23-7500852\1b415cce-aae0-400d-a6d5-86fab4fce9a9.jpg" /> = constant surface. Thus in the initial slice, the QLE Hamiltonian can be written as</p><disp-formula id="scirp.23118-formula70927"><label>(37)</label><graphic position="anchor" xlink:href="23-7500852\b8534288-56e8-4495-b83a-9d9680d5f86d.jpg"  xlink:type="simple"/></disp-formula><p>The reference frames’ quasilocal energy is a number, it just defines the zero point Hamiltonian. Thus, we replace the classical expressions by operators evaluated at the <img src="23-7500852\847466f4-86d1-4f44-9d4a-f84498a9aaa6.jpg" /> = constant slice. In the first approximation we simply take the <img src="23-7500852\e4b56863-74bd-4add-9745-2f5a3502b0fd.jpg" /> as classical</p><p><img src="23-7500852\441a940a-2500-43db-ab48-6ce7e4721cac.jpg" />as this arises due to the coordinate transformation and the norm of the vector <img src="23-7500852\f1435924-b4da-4a19-9338-521226dd59a3.jpg" /> in the previous frame. In the rewriting of (37) in regularised LQG variables the Hamiltonian appears rather complicated.</p><p>One can rewrite these in a much simpler form, using the apparent horizon equation. Since the Hamiltonian is an integral over the horizon, the variables will satisfy the apparent horizon Equation (10) upto quantum fluctuations. Thus the Hamiltonian operator is then re-written as</p><disp-formula id="scirp.23118-formula70928"><label>(38)</label><graphic position="anchor" xlink:href="23-7500852\adf4168a-907a-4398-803e-791df7de7baa.jpg"  xlink:type="simple"/></disp-formula><p>where we have used the classical apparent horizon equation (10) (with<img src="23-7500852\4285f334-0e12-466e-abf4-2eb20c78181c.jpg" />).</p><disp-formula id="scirp.23118-formula70929"><label>(39)</label><graphic position="anchor" xlink:href="23-7500852\0b51a4bd-6e59-47a3-97a1-f0d7533d547d.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="23-7500852\f0b4de08-9c07-4d58-94c8-738c123318b3.jpg" /> consists of some dimensionless constants <img src="23-7500852\4761a4b2-9b91-41b5-bdc6-dbf00533c8be.jpg" /> is the 2-dimensional area bit over which <img src="23-7500852\68351ea2-3f9f-4913-acf6-2da552737b1e.jpg" /> is smeared, <img src="23-7500852\ca6143a3-7780-4c9a-8251-fa024801a7f9.jpg" />is a dimensionfull constant which appears to get the <img src="23-7500852\a236051a-7ac6-4218-8847-a01ad278a02d.jpg" /> dimension less. <img src="23-7500852\dbaf7ef5-0f68-460f-bfb9-ba0d510fb5fe.jpg" />is the length for the angular edge <img src="23-7500852\161ed4d0-e1da-46ee-b8f3-1977425950d6.jpg" /> over which the gauge connection is integrated to obtain the holonomy. The sum over <img src="23-7500852\e3c8a72d-303c-4ead-8d2a-731866e022f3.jpg" /> is the set of vertices immediately outside the horizon. The (39) can then be lifted to an operator.</p><p>This regularised expression for QLE is for the horizon 2-surface only and would not apply for any other spherical surface in the Schwarzschild space-time.</p></sec><sec id="s3_3"><title>3.3. U(1) Case</title><p>Let us take the U(1) case to make the calculations easier and observe the action of the QLE Hamiltonian on the evolution of the coherent state. The spin network states are replaced by<img src="23-7500852\3a8af1d4-86f4-438a-8471-54a239377e68.jpg" />, <img src="23-7500852\4b2fe273-559f-483b-8c7b-c2e1367ed248.jpg" />, n is an integer and the coherent states are:</p><disp-formula id="scirp.23118-formula70930"><label>(40)</label><graphic position="anchor" xlink:href="23-7500852\dcd4899e-833b-4ac3-9e81-8b275a5569cc.jpg"  xlink:type="simple"/></disp-formula><p><img src="23-7500852\0e3dbb45-76f6-4a09-be1e-5cf4bcfdab28.jpg" />is the complexified phase space element in the “n-th” representation.</p><p>The QLE operator also takes the simplified form</p><disp-formula id="scirp.23118-formula70931"><label>(41)</label><graphic position="anchor" xlink:href="23-7500852\7cc910ec-eda0-405f-b19d-6d65a5832d16.jpg"  xlink:type="simple"/></disp-formula><p>The prefactors have been clubbed into<img src="23-7500852\b66e0010-a29d-41e5-94c2-ca0791fd2bcd.jpg" />.</p><p>In the calculation of the matrix elements, we drop the label of the edges <img src="23-7500852\2b31067b-a204-4cf4-b86c-182d546938a4.jpg" /> for the Hamiltonian.</p><disp-formula id="scirp.23118-formula70932"><label>(42)</label><graphic position="anchor" xlink:href="23-7500852\6df19604-9487-4da6-92bf-333872bf2569.jpg"  xlink:type="simple"/></disp-formula><p>This calculation can be done by putting an assumption that the<img src="23-7500852\2cb1ed1f-9183-4542-8818-defaa199e88a.jpg" />. In this <img src="23-7500852\4582920c-72d6-4487-b0b2-68ed6402781d.jpg" /> are completely independent of<img src="23-7500852\c97a901f-71e3-4c81-a990-d0a951f59e5f.jpg" />. It is an allowed assumption, and identifies the <img src="23-7500852\654f9ac9-640a-4000-b10a-466c1727fc7a.jpg" /> dependence of the operator matrix elements, which are otherwise “hidden”. The calculation however introduces an arbitrariness in the formula, which can be fixed by requiring that the expectation value of the Hamiltonian agrees with the classical QLE [<xref ref-type="bibr" rid="scirp.23118-ref10">10</xref>]. However, in this paper we use the “annihilation” operators defined in [<xref ref-type="bibr" rid="scirp.23118-ref16">16</xref>].</p><p>This is done by observing that the U(1) coherent states are eigenstates of an annihilation operator defined thus:</p><disp-formula id="scirp.23118-formula70933"><label>(43)</label><graphic position="anchor" xlink:href="23-7500852\55cad82b-f683-4331-9c0b-8975a0235b16.jpg"  xlink:type="simple"/></disp-formula><p>The holonomy operator can thus be written as</p><disp-formula id="scirp.23118-formula70934"><label>(44)</label><graphic position="anchor" xlink:href="23-7500852\b3064f8d-372c-4d04-8b9c-8968a0b2d765.jpg"  xlink:type="simple"/></disp-formula><p>And the derivative wrt Immirzi parameter of the holonomy which appears in the definition of the Hamiltonian replaced by</p><disp-formula id="scirp.23118-formula70935"><label>(45)</label><graphic position="anchor" xlink:href="23-7500852\326763c1-ee5b-469f-92bd-d0d26ee2c9cc.jpg"  xlink:type="simple"/></disp-formula><p>The dependence of the operator <img src="23-7500852\0d177c93-8551-4946-a9ba-153face87a74.jpg" /> on the Immirzi parameter is known (2), and thus we could evaluate the derivative</p><p><img src="23-7500852\f54ac54b-f52e-419b-bff1-8669b6a7465e.jpg" /></p><p>The term</p><disp-formula id="scirp.23118-formula70936"><label>(46)</label><graphic position="anchor" xlink:href="23-7500852\38253e3a-e005-4157-ad8f-e7f074772dc2.jpg"  xlink:type="simple"/></disp-formula><p>is then computable. Let us take the first term of (41) and find (46). As<img src="23-7500852\f788515f-85da-4bf6-b00d-9f4ec4dff2ff.jpg" />, (46) gives simply (we drop the “e” label for brevity)</p><p><img src="23-7500852\d515ba31-67aa-40ac-a21a-b3337c23b4e0.jpg" /></p><p>We then concentrate on the 2nd term of the above</p><disp-formula id="scirp.23118-formula70937"><label>(47)</label><graphic position="anchor" xlink:href="23-7500852\d689e7f1-4063-4ec9-8a89-e987bdd50a29.jpg"  xlink:type="simple"/></disp-formula><p>where we have used the fact that coherent states resolve unity. It can be shown that the expectation value of the operators in the <img src="23-7500852\fcb75815-baf0-43df-aa71-b204780ce913.jpg" /> collapses the integral to <img src="23-7500852\24a6d937-d39d-4663-8098-c1c134b8c6c6.jpg" /> point [<xref ref-type="bibr" rid="scirp.23118-ref16">16</xref>]. Thus one obtains from the above</p><disp-formula id="scirp.23118-formula70938"><label>(48)</label><graphic position="anchor" xlink:href="23-7500852\9e13bd97-bc0c-4c86-92bc-fc2db4e18a71.jpg"  xlink:type="simple"/></disp-formula><p>which is real, and thus</p><disp-formula id="scirp.23118-formula70939"><label>(49)</label><graphic position="anchor" xlink:href="23-7500852\110086b9-be0d-415a-a681-7f02e32f5648.jpg"  xlink:type="simple"/></disp-formula><p>this is actually the classical QLE as it should be from<img src="23-7500852\2d6d257f-a2b2-44ae-ad2b-5b9247d46ea2.jpg" />.</p><p>This is obvious, as the way the Hamiltonian is defined, this is simply a function of the Hilbert space outside the horizon, and the matrix elements of this will not yield anything new. We approximated the horizon sphere by summing over <img src="23-7500852\70baaae3-b13f-45ff-9c68-00f256fca095.jpg" /> vertices immediately outside the horizon. We could do the same by summing over <img src="23-7500852\68d78bb4-0940-42e7-bc1b-c1c35f991f78.jpg" /> vertices immediately within the horizon. For the Lemaitre slice, the metric is smooth at the horizon, and one can take the “quantum operators” evaluated at the vertex<img src="23-7500852\50f2b721-4e00-4627-893e-51cf0faf9f84.jpg" />. In this case however, as the region is within the classical horizon, the norm of the Killing vector is negative, and <img src="23-7500852\47bc592d-31d2-45e4-bed0-40bf013914ab.jpg" /> has components which are imaginary. The<img src="23-7500852\f116ea40-4dec-45c4-a2b1-d979d7a29e0f.jpg" />.</p><p>Thus<img src="23-7500852\5513569e-8ee0-41e1-a7f1-74459f070720.jpg" />. In the evaluation of the QLE, the energy would emerge correct in the <img src="23-7500852\20dd2de8-d38e-4e7b-a7b3-0aa32ca3c4e5.jpg" /> limit as <img src="23-7500852\26998939-471d-4739-8f4a-1afc4e68f162.jpg" /> The regularised Hamiltonian is not Hermitian, and the evolution equation is</p><disp-formula id="scirp.23118-formula70940"><label>(50)</label><graphic position="anchor" xlink:href="23-7500852\029d7149-5fc2-4681-9fb5-b622b37d9670.jpg"  xlink:type="simple"/></disp-formula><p>And thus the operator which appears in the change of entropy equation is</p><disp-formula id="scirp.23118-formula70941"><label>(51)</label><graphic position="anchor" xlink:href="23-7500852\aafc73f8-8560-4701-ac8f-e95cf75c15c5.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.23118-formula70942"><label>(52)</label><graphic position="anchor" xlink:href="23-7500852\d858bc9f-52dd-4dce-8af5-83346afa2326.jpg"  xlink:type="simple"/></disp-formula><p>The “rate of change” of entropy is thus</p><disp-formula id="scirp.23118-formula70943"><label>(53)</label><graphic position="anchor" xlink:href="23-7500852\ab9ac171-24bb-4c86-b119-8b8ad6f10b5d.jpg"  xlink:type="simple"/></disp-formula><p>we extracted the <img src="23-7500852\b7b16adf-a66a-49a9-99a7-eb8a1d09179a.jpg" /> from <img src="23-7500852\cb239d52-a068-408b-bd0c-98012e5cbcc7.jpg" /> to get <img src="23-7500852\906ed22f-2ae4-44bb-ae7e-6b300c01a09a.jpg" /> and rewrote the rest of the constants as<img src="23-7500852\f3e3c011-73b4-418c-921e-a30ad5b33f45.jpg" />.</p><p>Thus there is a net change in entropy, but, to see if this is Hawking radiation, we have to couple matter to the system.</p></sec><sec id="s3_4"><title>3.4. SU(2) Case</title><p>The SU(2) case is easily reduced to the U(1) case in the actual calculation due to the gauge fixing. This is achieved by making the following observations: To retain the metric as in the same form as the classical metric, we impose the conditions at the operator level</p><disp-formula id="scirp.23118-formula70944"><label>(54)</label><graphic position="anchor" xlink:href="23-7500852\b3736f5c-09b7-4159-a512-88d016f1b1af.jpg"  xlink:type="simple"/></disp-formula><p>such that the corresponding metric has only the diagonal terms as non-zero. With these additional “constraints” on the operators, we can put the <img src="23-7500852\28d7312d-1742-4828-afb7-46bad307ef5d.jpg" /> such that each has only one component surviving, let’s say<img src="23-7500852\959f1367-dd5a-4251-b994-b98102805642.jpg" />. This also makes the holonomy restricted to the U(1) case, as the gauge connection <img src="23-7500852\00b6324b-4cf2-4b0a-b553-15b9c0b31ede.jpg" /> gets restricted to the <img src="23-7500852\fdfcb76c-ca1d-4956-a079-a705fd07c71a.jpg" /> and other directions can be put to zero. Thus we can take the holonomy to be diagonal. If the holonomy matrix is off-diagonal the U(1) projection still works out to be the same</p><disp-formula id="scirp.23118-formula70945"><label>(55)</label><graphic position="anchor" xlink:href="23-7500852\206dde22-fce5-4e3b-9397-fb5de07dafbe.jpg"  xlink:type="simple"/></disp-formula><p>The operator is then obtained as</p><p><img src="23-7500852\88fbadbb-918c-4a48-94a6-df80c5dac141.jpg" /></p><p><img src="23-7500852\ce0c5037-1c23-4f8c-8523-bc52e174d2ba.jpg" /></p><disp-formula id="scirp.23118-formula70946"><label>(56)</label><graphic position="anchor" xlink:href="23-7500852\1bd8b200-0bbc-4f89-9138-3480edd32eca.jpg"  xlink:type="simple"/></disp-formula><p>This is same as the U(1) Hamiltonian (upto normalisations). The spin network states also project on to U(1) subgroup, thus giving us the same techniques to use in the calculation of the U(1) states as for this one. To observe this, the non-zero elements for the holonomy matrix</p><disp-formula id="scirp.23118-formula70947"><label>(57)</label><graphic position="anchor" xlink:href="23-7500852\17c64933-9e69-4bc6-9f0e-7e37640c4349.jpg"  xlink:type="simple"/></disp-formula><p>in the j-th representation is given by:</p><disp-formula id="scirp.23118-formula70948"><label>(58)</label><graphic position="anchor" xlink:href="23-7500852\e4b49c7b-c90b-4fae-91e3-da2e781622e7.jpg"  xlink:type="simple"/></disp-formula><p>Clearly in the particular case we are considering, the<img src="23-7500852\8fff13c1-f8b5-4a0e-96fe-e55f229171a0.jpg" />, and m, n = –j and j. Thus the two non-zero elements are</p><disp-formula id="scirp.23118-formula70949"><label>(59)</label><graphic position="anchor" xlink:href="23-7500852\372a07cb-53d6-459f-b881-02454f46f805.jpg"  xlink:type="simple"/></disp-formula><p>The sum over j in the <img src="23-7500852\4dfda8c0-3e9d-447d-a294-1a95ac25f335.jpg" /> with the coherent state defined in (3) thus reduces to the U(1) case in the computation of the change in entropy. Thus the rate of change in entropy of a classically spherically symmetric black hole is given by</p><disp-formula id="scirp.23118-formula70950"><label>(60)</label><graphic position="anchor" xlink:href="23-7500852\268de829-5e3d-4e63-9065-60abf2adf243.jpg"  xlink:type="simple"/></disp-formula><p>where the classical holonomies<img src="23-7500852\823dbd0c-3fe8-4ebe-ba42-0f177e526f8e.jpg" />. If we plug in the actual values, we get this to be</p><disp-formula id="scirp.23118-formula70951"><label>(61)</label><graphic position="anchor" xlink:href="23-7500852\c010df67-3140-46c5-a1d7-132b14cfa562.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="23-7500852\ee2266a9-8494-45e5-8253-011b8bcacd8c.jpg" /> the area element at vertex <img src="23-7500852\2bfca7e0-45b0-4115-a50b-0ebbbc8fa788.jpg" /> on the sphere. This change in entropy is totally gravitational in origin, and seems to signify the emergence of “geometry” from within the horizon.</p><p>In fact, if we some over the area, we get the <img src="23-7500852\0f4914c8-3fe4-4a10-82dd-d047fef53604.jpg" />(if we set<img src="23-7500852\c32fe04e-8593-4a15-b7f2-a0865a77c61b.jpg" />), which would be the change in entropy when the radius of the horizon changes by<img src="23-7500852\76e2f36d-3904-4171-baa6-98d61db6251c.jpg" />!</p></sec></sec><sec id="s4"><title>4. Outgoing Flux of Radiation</title><p>In the previous section we found that as the system evolves in time, the horizon fluctuates and the area decreases. But is this Hawking radiation? Adding matter to a “coherent state” description of semiclassical gravity has been discussed [<xref ref-type="bibr" rid="scirp.23118-ref17">17</xref>]. Thus, given a massless scalar field Lagrangian coupled to gravity, whose Hamiltonian is given by</p><disp-formula id="scirp.23118-formula70952"><label>(62)</label><graphic position="anchor" xlink:href="23-7500852\b285ffd1-ec37-49bf-959f-d78f263f2e13.jpg"  xlink:type="simple"/></disp-formula><p>the “gravity” in the Hamiltonian can be regularised in terms of the <img src="23-7500852\1960d99f-1515-4634-9680-f6ae6d939531.jpg" /> operators in the coherent state formalism. The integral over the three volume gets converted to a sum over the vertices dotting the region. Thus</p><disp-formula id="scirp.23118-formula70953"><label>(63)</label><graphic position="anchor" xlink:href="23-7500852\ca75b701-8f6c-4a49-9a85-309a73a79b0f.jpg"  xlink:type="simple"/></disp-formula><p>This Hamiltonian is an operator, and one evaluates an expectation value of the Hamiltonian in the reduced density matrix of the initial slice, to find the classical behavior of the scalar field as observed by an observer outside the horizon. Thus</p><disp-formula id="scirp.23118-formula70954"><label>(64)</label><graphic position="anchor" xlink:href="23-7500852\2d223230-8e5e-4d2a-a1a0-6e6d55165aac.jpg"  xlink:type="simple"/></disp-formula><p>This Hamiltonian and the density matrix are then both evolved according to the time-like observers frame. One gets</p><disp-formula id="scirp.23118-formula70955"><label>(65)</label><graphic position="anchor" xlink:href="23-7500852\5583b456-8aad-4cce-a42f-09dc1d0aa201.jpg"  xlink:type="simple"/></disp-formula><p>This gives</p><disp-formula id="scirp.23118-formula70956"><label>(66)</label><graphic position="anchor" xlink:href="23-7500852\ea72b403-5147-43cc-8888-bbc8ef387746.jpg"  xlink:type="simple"/></disp-formula><p>It is very clear thus that the order <img src="23-7500852\e05f708a-aea7-46b8-8369-460eaa7e498d.jpg" /> terms are zero for this. However, allowing for the non-unitary evolution using the non-Hermitian Hamiltonian, the <img src="23-7500852\bccf6eb0-1fbe-4571-a512-b04d037303d4.jpg" /> terms survive. In fact the terms are</p><disp-formula id="scirp.23118-formula70957"><label>(67)</label><graphic position="anchor" xlink:href="23-7500852\279108ca-aea0-4052-8060-beeb4e04dba5.jpg"  xlink:type="simple"/></disp-formula><p>The first term is remarkable, it shows that the term giving rise to entropy change teams up with the expectation value of the scalar Hamiltonian. The second term yields corrections, so we ignore it in the first approximation. The exact details of the computation have to be obtained using the coherent state of the matter and gravity coupled system [17,18]. If one simple takes the matter + gravity system in a tensor product form, and one has matter quanta of energy <img src="23-7500852\c094fcc2-0a21-4841-bf0b-ecf12f0d5ac2.jpg" /> sitting at one vertex, then the first term would give new matter in the evolved slice as<img src="23-7500852\94cc79ae-6a04-45a7-b95e-20bc616886ca.jpg" />. The “rate” of particle creation thus has the form <img src="23-7500852\dc345634-debe-4c92-b3a4-60f1f3bff96b.jpg" /> where <img src="23-7500852\85ba3e95-57de-41b5-8313-32d59133b685.jpg" /> is the Hawking temperature for the signs <img src="23-7500852\466b3e82-c6b2-4555-a6e7-08783cec9679.jpg" /> and negative (positive)<img src="23-7500852\85dcadfb-ab1c-45c8-8637-b1962321422e.jpg" />.</p><p>Thus from the above it seems 1) One has found emission of matter quanta from a black hole but from a “semiclassical” description rooted in a theory of quantum gravity, beyond quantum fields in curved space-time.</p><p>2) The results indicate a non-unitary evolution which allows space to emerge from within the horizon.</p><p>3) The emission is perceived by a static or an accelerating observer as anticipated, and the non-unitary flow might be due to the semiclassical approximations. A quantum evolution using the quantum Hamiltonian might still be unitary.</p><p>The above derivation seems to be a “quantum gravity” description of the tunneling mechanism for describing Hawking radiation [<xref ref-type="bibr" rid="scirp.23118-ref19">19</xref>]. However, the results are preliminary and further investigation has to be done.</p></sec><sec id="s5"><title>5. Conclusion</title><p>In this paper we showed a method to obtain the origin of Hawking radiation using a coherent state description of a black hole space-time. We took a quantum wavefunction defined on an initial slice, peaked with maximum probability at classical phase space-variables. We then evolved the slice using a Hamiltonian, which is the “quasilocal energy” at the horizon. This QLE evolved the system in the time and the entropy was shown to change, indicating a change in black hole mass and hence an emergence of interesting non-unitary dynamics. One of course has to investigate further to see what is the endpoint of this time evolution. The time flow indicates one might have to formulate quantum theory of gravity rooted in irreversible physics. The presence of additional degrees of freedom in the form of “graphs” also indicates that the classical phase space might not be described by deterministic physics, but by distributions, a manifestation of microscopic irreversible physics in complex systems.</p></sec><sec id="s6"><title>6. Acknowledgements</title><p>This research is supported by NSERC; research funds of University of Lethbridge. I would like to thank B. Dittrich for useful discussion; J. Supina for proofreading the manuscript.</p></sec><sec id="s7"><title>REFERENCES</title></sec></body><back><ref-list><title>References</title><ref id="scirp.23118-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">A. Ashtekar, New Perspectives in Canonical Gravity,” Bibliopolis, Napoli, 1988. </mixed-citation></ref><ref id="scirp.23118-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">C. Rovelli, “Quantum Gravity,” Cambridge University Press, Cambridge, 2004.  
doi:10.1017/CBO9780511755804</mixed-citation></ref><ref id="scirp.23118-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">T. Thiemann, “Introduction to Modern Canonical Quantum General Relativity,” Cambridge University Press, Cambridge, 2007. doi:10.1017/CBO9780511755682</mixed-citation></ref><ref id="scirp.23118-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">B. Hall, “The Segal-Bargmann ‘Coherent State’ transform for Compact Lie Groups,” Journal of Functional Analysis, Vol. 122, No. 1, 1994, pp. 103-151.  
doi:10.1006/jfan.1994.1064</mixed-citation></ref><ref id="scirp.23118-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">T. Thiemann, “Gauge Field Theory Coherent States (GCS): 1. General Properties,” Classical and Quantum Gravity, Vol. 18, No. 1, 2001, pp. 2025-2064.  
doi:10.1088/0264-9381/18/11/304</mixed-citation></ref><ref id="scirp.23118-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">A. Dasgupta, “Semiclassical Quantisation of Space-Times with Apparent Horizons,” Classical and Quantum Gravity, Vol. 23, No. 3, 2006, pp. 635-671.  
doi:10.1088/0264-9381/23/3/007</mixed-citation></ref><ref id="scirp.23118-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">A. Dasgupta, “Coherent States for Black Holes,” Journal of Cosmology and Astroparticle Physics, Vol. 3, No. 8, 2004, pp. 1-36. </mixed-citation></ref><ref id="scirp.23118-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">A. Dasgupta, “Semiclassical Horizons,” Canadian Journal of Physics, Vol. 86, No. 4, 2008, pp. 659-662. </mixed-citation></ref><ref id="scirp.23118-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">A. Dasgupta and H. Thomas, “Correcting Gravitational Entropy of Spherically Symmetric Horizons,” 2006.  
arXiv:gr-qc/0602006 </mixed-citation></ref><ref id="scirp.23118-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">A. Dasgupta, “Entropic Origin of Hawking Radiation,” Proceedings of 12th Marcel Grossman Meeting, Paris, July 2009. </mixed-citation></ref><ref id="scirp.23118-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">J. Brown and J. York, “Quasilocal Energy and Conserved Charges Derived from the Gravitational Action,” Physical Review D, Vol. 47, No. 4, 1993, pp. 1407-1419.  
doi:10.1103/PhysRevD.47.1407</mixed-citation></ref><ref id="scirp.23118-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">B. Dittrich and T. Thiemann, “Testing the Master Constraint Program for Loop Quantum Gravity: I. General Framework Class,” Classical and Quantum Gravity, Vol. 23, No. 4, 2006, pp. 1025-1066.  
doi:10.1088/0264-9381/23/4/001</mixed-citation></ref><ref id="scirp.23118-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">B. Bahr and T. Thiemann, “Gauge-Invariant Coherent States for Loop Quantum Gravity. II. Non-Abelian Gauge Groups,” Classical and Quantum Gravity, Vol. 26, No. 4, 2009, p. 045012. </mixed-citation></ref><ref id="scirp.23118-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">A. Ashtekar, J. Baez, K. Krasnov and A. Corichi, “Quantum Geometry and Black Hole Entropy,” Physical Re- view Letters, Vol. 80, No. 5, 1998, pp. 904-907.  
doi:10.1103/PhysRevLett.80.904</mixed-citation></ref><ref id="scirp.23118-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">S. W. Hawking and G. Horowitz, “The Gravitational Hamiltonian, Action, Entropy and Surface Terms,” Classical and Quantum Gravity, Vol. 13, No. 6, 1996, pp. 1487- 1498. doi:10.1088/0264-9381/13/6/017</mixed-citation></ref><ref id="scirp.23118-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">T. Thiemann and O. Winkler, “Gauge Field Theory Coherent States (GCS) III: Ehrenfest Theorems,” Classical and Quantum Gravity, Vol. 18, No. 21, 2001, pp. 4629-46841.  
doi:10.1088/0264-9381/18/21/315</mixed-citation></ref><ref id="scirp.23118-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">H. Sahlmann and T. Thiemann, “Towards the QFT on Curved Space-Time Limit of QGR. 1. A General Scheme,” Classical and Quantum Gravity, Vol. 23, No. 3, 2006, pp. 867-908. doi:10.1088/0264-9381/23/3/019</mixed-citation></ref><ref id="scirp.23118-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">H. Sahlmann and T. Thiemann, “Towards the QFT on Curved Space-Time Limit of QGR. 2. A Concrete Implementation,” Classical and Quantum Gravity, Vol. 23, No. 3, 2006, pp. 909-954. doi:10.1088/0264-9381/23/3/020</mixed-citation></ref><ref id="scirp.23118-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">M. K. Parikh and F. Wilczek, “Hawking Radiation as Tunneling,” Physical Review Letters, Vol. 85, No. 24, 2000, pp. 5042-5045. doi:10.1103/PhysRevLett.85.5042</mixed-citation></ref></ref-list></back></article>