On Horizonless Temperature with an Accelerating Mirror

A new solution of a unitary moving mirror is found to produce finite energy and emit thermal radiation despite the absence of an acceleration horizon. In the limit that the mirror approaches the speed of light, the model corresponds to a black hole formed from the collapse of a null shell. For speeds less than light, the black hole correspondence, if it exists, is that of a remnant.


Introduction: Some Puzzles with Moving Mirrors and Evaporating Black Holes
The discovery of Hawking radiation from a black hole came at a surprise over 40 years ago [1], since it implies that given enough time, black holes in asymptotically flat spacetimes would radiate off their mass; and if no new physics comes into play, they would eventually evaporate completely. This leads to the information paradox [2] -where did the information that falls into a black hole disappear to, if the black hole disappears completely? Attempts to recover information from evaporating black holes continue to produce new paradoxes, such as the firewall controversy [3][4][5], which threatens the conventional understanding that there should be nothing unusual -much less a diverging energy density -at the low curvature spacetime region near the event horizon of a sufficiently massive black hole. In fact, the situation is worse than previously thought: if a firewall does exist, then in principle it can be much further away from the event horizon [6], so that an unexpected space traveler would hit it and be burned to death even if he is nowhere near a black hole. While the resolution of the information paradox may require us to fully understand how to unify general relativity with quantum field theory, it is certainly possible that progress can nevertheless be made without a working theory of quantum gravity 2 . In recent years, it has become increasingly clear that we at least need to understand the subtle "physics of information" [8][9][10], and how it applies to black holes [11][12][13][14][15][16][17], in order to understand the information retrieval process from Hawking quanta.
In fact, a central piece of the puzzle regarding information loss, is the understanding of entanglement entropy, S. Let us consider the formation of a black hole from a gravitational collapse of some matter in a pure state. It is often believed that the entanglement entropy of the Hawking radiation received at null infinity, which should be zero at the beginning before any radiation arrives, should first increase, but then decrease at some point (known as the "Page time"), so that eventually the entanglement entropy vanishes. Such a "Page curve" [18][19][20] -the plot of the entanglement entropy against timeis crucial, since it gives insight into how information may be retrieved from the highly scrambled Hawking emission. There are, however, subtleties that are often overlooked. Notably, any calculation of entanglement entropy necessitates regularizing ultraviolet divergences. One notices that imposing a cutoff is a tricky procedure since modes that have sufficiently high energy at some point in the spacetime can be redshifted at some other point due to spacetime curvature. This implies that a mode that is beyond a cutoff scale can be red-shifted below the scale, so the cutoff is not well-defined [21,22]. Furthermore, the results obtained could depend on the cutoff scheme. Progress has been made recently with the introduction of the "causal-splitting regularization" scheme of Bianchi and Smerlak [23][24][25], which allows one to compute the production of entanglement entropy in a cutoff-independent manner. However, even then, there are still a few puzzles regarding the entanglement entropy of an evaporating black hole. We list two such puzzles below.
(1) Negative Energy Flux and Negative Entropy: Firstly, it has been observed, e.g. by Bianchi and Smerlak [23,24], that mass loss of a black hole, assuming unitarity, is not monotonic. In other words, at some point in time, the mass of an "evaporating" black hole actually increases. For an observer at infinity, Hawking radiation reduces the mass of a black hole by emitting particles. This means that an asymptotic observer sees a flux of positive energy coming out from the black hole. If the mass increases, albeit briefly 3 , during the course of Hawking evaporation, this would correspond to the emission of negative energy flux (hereinafter, "NEF") from the black hole. (NEF emission from evaporating black holes, at least in the case of a (1+1)-dimensional dilaton gravity, was already known in the literature for over 20 years [27].) Curiously, an observer equipped only with a particle detector would not be able to see any sign in the spectrum of the Hawking quanta to know that NEF has been emitted [28]. Even more surprising than the emission of NEF, in several models of evaporating black holes, Bianchi and Smerlak have explicitly shown that the entanglement entropy obtained from their causal-splitting scheme can also become negative at late times in the course of the evaporation. The physical interpretation of negative entropy is dependent on the particular type of entropy. Their physical relationships are subtle and not fully clear: in the literature it has been interpreted as either the radiation being less correlated than the vacuum [25], or the radiation being more correlated than the vacuum [29].
(2) Unitarity and Information Recovery: Secondly, even with the fully covariant cutoffindependent regularization of Bianchi and Smerlak, it was shown that [28] their criterion only requires that the entanglement entropy tends to a constant asymptotically (both in the far past and in the late future), instead of a more stringent requirement that the entropy should tend to zero so as to recover the pure state. In fact, it is not so surprising that the entanglement entropy can increase monotonically. Even with the causal-splitting regularization, as long as there is a cutoff, then modes can be redshifted by spacetime curvature, and consequently a mode that is below the cutoff scale at future null infinity at some late time u, when traced back to the past null infinity, could be well-above the cutoff. Since there is no longer a one-to-one map between early modes and late time modes, it is not surprising that one does not obtain S(u) → 0 at late times. Note that such a phenomenon can be shown to occur even in a simplified model involving a moving mirror in flat spacetime -as was explicitly shown in [22,28] -in which the quantum field theory is unitary. It was suggested in [28] that the corresponding black hole picture might be a black hole remnant [30,31], and that unitarity is maintained in the sense that if one takes into account both the exterior and the interior of a black hole, then the entire quantum state is pure at all times. In such a picture, information remains hidden inside the ever-shrinking black hole horizon (the interior spacetime can still have a large volume [21,[31][32][33][34][35][36], see however [37]) and the radiation is never purified. Since the end state is a remnant, the information inside is never destroyed 4 . It would be interesting to investigate if S(u) → const. = 0 at late times necessarily implies that the corresponding black hole model ends as a remnant. A main question here, however, is the following: do we understand entanglement entropy enough to base the debate of information loss paradox of black holes on it? The subtleties here suggest that we should be more careful in dealing with this issue. While Hawking evaporation has been a mainstay of the "quantum fields under external conditions" enterprise, the moving mirror model [38,39] has flourished mostly as an ancillary effect. However, we feel that it is prudent to first understand a moving mirror model in which both issues (1) and (2) can be investigated, after all, a (1+1)-dimensional flat spacetime with a mirror trajectory is much simpler than a (3+1)-dimensional black hole spacetime with non-trivial curvature 5 . A trajectory of a moving mirror in spacetime is a reflecting boundary on which field modes are constrained. A moving mirror excites the modes, thereby producing particles. The spectrum of the emission depends on the 4 Note that there are two types of black hole remnant in the literature: the "long-lived" or "metastable" remnant, and the "eternal" remnant [31] . The former has lifetime much longer than that of a black hole, that is, proportional to M n where n > 3. An eternal remnant, on the other hand, lives forever. Our model corresponds to an eternal red-shifted scenario. 5 Einstein's gravity is topological in (1+1)-dimensions, but with a suitable coupling to matter fields, gravity need not be trivial in (1+1)-dimensions. This allows one to study black hole evaporation. In our moving mirror model, of course, there is no gravity, only an accelerating mirror, so there is no complication either. exact trajectory of the mirror. With a suitably chosen trajectory, we can therefore mimic the particle production from, say, a collapsing star. (Readers who are not familiar with moving mirrors can refer to [40] for a pedagogical exposition.) Despite many decades of research, physical interpretations of external-potential type problems, like the moving mirrors, have never been made entirely clear. Moreover, serious confrontations with questions that are evaded in these non-gravitational analogs become inevitable during the study of gravitational problems. In other words, we would like to emphasize that one should not hope that a moving mirror model can fully resolve information paradox of a bona fide black hole, but understanding the subtleties of quantum field theory in a moving mirror model is a first step toward the more complicated physics of black holes.
We have come a long way in the study of moving mirrors. Notably, Moore [41], DeWitt [42], and later on, Davies and Fulling [38,39] initiated a program using field theories with external conditions, which eventually demonstrated that quantities like the expectation values of the stress-energy tensor, and the localization of particles using wave packets, can be calculated in various physical problems and used to extract significant physical consequences of the quantum fields. Indeed, there has been renewed interest [23,24,[43][44][45][46][47][48][49][50][51][52][53][54] in the moving mirror model in recent years due, in part, to claims of experimental verification of the dynamical Casimir effect [55,56]. In many of these particle production scenarios, various systems that exploit the simple mathematical set up of the (1+1)-dimensional moving mirror model have led to novel experimental designs. (See [54] for one of the latest proposals.) As previously mentioned, a central advantage to the moving mirror model is its simplicity. This is both in general, and in the context of the recent one-to-one correspondence with a black hole [57][58][59][60], which found that for one concrete example, the particle production is exactly the same in both the mirror and black hole cases in (1+1)-dimensions. The simplicity of (1+1)-dimensions allows the crux of Hawking radiation to stand out more clearly, separated from the specialized details associated with higher dimensional curved geometry and back-scattering. With the one-to-one correspondence, the moving mirror model can be treated as an even more precise analogy to Hawking's original argument, and therefore it is of interest to extend the one-to-one correspondence to more physically realistic circumstances, while holding on to this fortunate simplicity. However, even as a relatively simple theoretical model of black hole evaporation in (1+1) -dimensions, the moving mirror model, has in practice, been very hard to extend to solutions for exact trajectories where the global Bogoliubov coefficients may be evaluated. Few solutions have been found 6 and finite-nonzero-energy cases are scarce. 7 Nevertheless, mirror trajectories that produce finite amount of energy are precisely those that are physically more realistic. We are therefore interested in such trajectories.
Although there is pedagogical value associated with the reduced complexity of the direct and straightforward calculations in the moving mirror model in (1+1)-dimensions, one should note that there is a peculiarity in (1+1)-dimensions, namely that Minkowski spacetime in (1+1)-dimensions have two sets of past and future null infinities. The right past infinity will be denoted by I − R , and the right future infinity I + R . Similarly, we also have I − L and I + L . That is to say, its Penrose diagram is a diamond, whereas in higher dimensions, the Penrose diagram for Minkowski spacetime is a triangle 8 . The implication is this: if we want to understand quantum field theory of a moving mirror in (1+1)-dimensions, we should also take into account the left side of the mirror.
Yet another motivation to consider both sides of the mirror is the following: In (1+1)dimensional conformal field theory, there are left and right "temperatures", T L and T R , which are related to the amount of left and right-moving excitations in the field. In a thermal ensemble of states, the thermodynamic temperature is related to these left and right temperatures. The Cardy formula for microstate degeneracy [65], explicitly involves these temperatures, along with the left and right central charges c L and c R . So from this point of view of thermodynamics, it is natural to include both left and right-moving excitations.
The most important reason to consider both sides of the mirror, however, is related to the eventual aim to understand information loss of black holes. To do so, as we have argued, a crucial first step is to understand how entanglement works in the simpler case of a moving mirror. In particular, as far as unitarity is concerned, it is important that we include the entire spacetime, this is to ensure that information is not hiding in some part of spacetime that one might otherwise overlook.
In this work, we are mainly motivated in asking the following questions: Is there a moving mirror in (1+1)-dimensions, satisfying unitarity in the sense allowed by the Bianchi-Smerlak criterion (namely, S(u) → const. as u → ±∞), that has no acceleration horizon, produces finite amount of energy, and serves as an analog for Eddington-Finkelstein coordinate null-shell gravitational collapse in its limiting case? Furthermore, by looking at both sides of the mirror trajectory, can we understand, or at least reveal some additional features, regarding negative energy flux and negative entropy?
We shall now explain the reason we would like to consider only trajectories that have no acceleration horizon (hereinafter, we will simply refer to an acceleration horizon as "horizon", unless there is a risk of confusion.) The basic utility of various trajectories examined in the moving mirror model with respect to black hole radiation has, like black hole evaporation itself, also been well-known for over forty years. But the problem of relating Hawking's global construction calculation (see [66] for a stochastic route) to the physical mechanism responsible for the particle-energy creation effect, which is safely assumed to involve local curvature, has been mostly elucidated by various methods using the late-time Davies-Fulling black hole-moving mirror correspondence. However, with the advent of the all-time, exact, black hole-moving mirror correspondence [57], the physics of the boundary condition effect is more directly related to the curved space effect than previously supposed. We therefore hope to timely exploit this insight by extending the flat-space model to one with a boundary condition that does not contain an asymptotic horizon. It is clear from basic causality that the existence of a horizon in the future is not essential to the production of particles and energy at early times, and in this calculation we will explicitly demonstrate the particle production by localization of the spectra at early times. (In fact, the same causality argument works for the case of a black hole, i.e. the event horizon is not an essential feature for early-time particle creation. See also [67].) Beyond this, we will demonstrate that the production process at any time need not have an acceleration horizon in order to reach thermal equilibrium. The use of wave packets to obtain spectral resolution of thermal Hawking radiation in such a limiting-case horizonless model where the evaporation process stops at late times, has not been investigated before. It is important to recognize that this model is an extension of the particular moving mirror which has a one-to-one correspondence to the exactly solvable black hole case [57]. The model presented here is novel, in part, because its limit is that of the exact correspondence, which means that a certain acceleration parameter (see below) is appropriately scaled. A novel result of the removal of the horizon, for this unique trajectory, is that the finite total emission of energy is simple enough that it can be expressed analytically. Also important to this solution is the introduction of a new method to find other moving mirror solutions which may otherwise be intractable.
In this work, we are focused on the mirror model. Although the solutions and their novel features could represent various black hole collapse scenarios, it is premature to make any strong conclusion here. In addition, even though the moving mirrors that we investigate below have no acceleration horizon, this does not necessary mean that the corresponding "black holes" have no event horizon 9 (or at least a trapped surface of some sort). After all, a main purpose of moving mirror models is to reproduce particle emission of a black hole -it is the properties of the produced particles that are in correspondence -one should be very careful if one wishes to identify the mirror trajectories themselves to geometric properties of black holes such as the event horizons 10 .
The present work is organized in the following manner. A short introductory treatment of the moving mirror machinery is given in Section (2). It shows how to generalize the exactly solvable "black mirror" case to remove the horizon. The subject of interest in Section (3) is the explicit trajectory of our new mirror solution. In Section (4) we solve for the energy flux and the total energy production on both sides of the mirror. The energy flux is shown to have a spike on the left side and a plateau on the right side. We also quantify the equilibrium temperature. In Section (5) we calculate the entanglement entropy flux and confirm it remains consistent with the energy-entropy relation and unitarity. In Section (6) we investigate the stress tensor correlations for the model by solving for the correlation ratio R 1 exactly. This helps confirm thermal equilibrium. In Section 9 The possibility that black holes may not have an event horizon has been investigated by many authors, see [68] [69] and the references therein. 10 While it is true that an accelerated observer in Rindler spacetime sees Unruh temperature that is analogous to the Hawking temperature of a black hole, it does not follow that an acceleration horizon of a mirror always corresponds to a black hole horizon. Furthermore, the temperature in the mirror case, is measured by observers far away from the mirror, not a Rindler-like observer on the mirror trajectory.
(7) we calculate the particle production and investigate the spectral dynamics. Here we include consistency checks to verify the total energy produced via the stress tensor, agrees with the total energy produced from summing particle quanta. We find constant emission of particles emitted to the observer on the right for any arbitrary long period of time. We conclude with some discussions in Section (8). Throughout this work, we use the units G = c = = k B = 1.

The Machinery of Moving Mirrors
In Section (2.1), we shall first introduce some basic concepts necessary for understanding our construction of a horizonless solution that generalizes the "black mirror". We then discuss the removal of the acceleration horizon in Section (2.2). The solution satisfies three criteria: the presence of a horizonless temperature, an appropriately scaled acceleration parameter, and the termination of evaporation at late times.

Some Conventions of Moving Mirrors
As the simplest example of the dynamical Casimir effect, the moving mirror model also serves as a way to understand black hole evaporation by imposing an external boundary condition in 1+1 dimensions on the quantum field, rather than an external curved spacetime. Consider then, such a boundary that does not accelerate forever, starting and ending at time-like past infinity i − , and time-like future infinity i + , respectively; possessing asymptotically zero acceleration in both the far past and far future, and always moving slower than the speed of light. This fully asymptotically inertial mirror will contain no horizon. Thus, it will contain no pathological acceleration singularity either. The plot of such a trajectory is in Figure (1). A salient pay-off for horizon-removal is that the mirror system, in addition to being unitary (henceforth, by unitarity we always mean unitarity in the broad sense allowed by Bianchi-Smerlak criteria), produces only a finite amount of total energy, as we will demonstrate in Section (4) and Section (5).
The quantum field, Ψ, is the massless scalar of the Klein-Gordon equation Ψ = 0, whose value is zero, Ψ| z = 0, when evaluated at the position of the moving mirror, z(t). The modes, φ ω and ψ ω , are on equal footing in the sense that they can both be used to expand the field: The modes are orthonormal and complete and can be exactly solved in the (1+1)dimensional case: where the functions p(u) and f (v) are the usual notation for the ray-tracing functions, which are intimately related to the trajectory of the mirror z(t) itself, see [63] and Section (2.1.1). The famous Bogoliubov coefficients appear by expanding one set of modes in terms of the other set of modes, where with the flat space scalar product defined in null coordinates, (u, v), by The Bogoliubov coefficients α ω ω and β ω ω also give the operators a ω and a † ω in terms of the operators b ω and b † ω , while the orthonormality of the modes hold according to the usual convention, see [28] for more detail.

The Four Functions of Mirror Physics
There are four functions, v s (u), u s (v), x s (t), t s (x), which are useful for doing global calculations involving the aforementioned field modes. The first two are the ray-tracing functions (expressed in null coordinates), where u s (u) ≡ p(u) and v s (v) ≡ f (v), and the last two are the associated spacetime coordinate functions. The inverses are expressed like so: We shall collectively call all four of them, "shock wave functions" or "shock functions" for short, after the collapse of the null shell shock wavefront description to form a black hole. There are many other auxiliary functions, such as t s (v), v s (t), t s (u), u s (t). However, the original four functions of coordinates v, u, t, and x will prove efficient at calculating observables. The information about how the field modes become red-shifted due to external conditions is fully contained in these four functions. The relationships between them are demonstrated as follows.
First consider the usual null coordinates on Minkowski spacetime u ≡ t − x and v ≡ t + x, and their analogous auxiliary functions as functions of time, These contain the shock function x s (t), which is the trajectory of the mirror. The inverses of Eqs. (10), contain the shock functions, u s (v) and v s (u), Functional inverses should be obvious from the notation. Useful auxiliary inverses are: . The total energy emitted, the energy flux, and the beta Bogoliubov coefficients have expressions that are conveniently written in terms of the four shock functions.

How to Remove a Horizon
There is an easy way to remove the horizon, (recall that c = 1), from a future asymptotically null moving mirror trajectory, z(t). The idea is to modify this so that the new trajectory, z(t), has where 0 < ξ < 1 is the future asymptotically drifting speed. This can be achieved by writing the horizonless trajectory, z(t), in terms of the trajectory with a horizon (henceforth 'horizon trajectory'), z(t), This works ifż (t) = ξż(t).
Taking this approach helps answer whether the particle spectra can (1) reach equilibrium for an extended period of time, and (2) proceed to shut off. The mirror does not strictly have a null horizon, yet as we will see, it can still achieve a "thermal plateau" (i.e. the emission is virtually thermal for some arbitrary finite amount of time). This approach also ensures (3) the correct scale for the acceleration parameter κ (not to be confused with the physical acceleration, see below). A correct scale is critical for the red-shifting of the modes to correspond to the exactly solvable black hole case [57][58][59][60] in the limit ξ → 1. This automatically extends the mirror in the black hole-moving mirror correspondence by promoting it to a more physical footing where the total evaporation energy is finite and unitarity is preserved. While we have found the mirror solution that meets these strict requirements, a possible black hole counterpart calculation is beyond the scope of this work. In the model we are about to present, we do not claim that it actually corresponds to any realistic evaporating black hole spacetime. For the present work we only seek a simple mirror model in which the three conditions presented above are met, so that we may study the energy and energy flux, the entropy, the correlations, and the particle spectra, together in the absence of a horizon. It may or may not have an exactly tractable black hole correspondence. In a subsequent work this will be investigated, but as we have emphasized in the Introduction, even if it has such a black hole correspondence, the absence of horizon in the mirror model does not necessary entail the absence of any event horizon or trapped surface in the black hole geometry. It is worth pointing out that this mirror solution is new -it is the first explicit demonstration of a unitary solution with a thermal plateau 11 that has limiting red-shifting functions which correspond to the black hole-moving mirror system in [57].
The information contained in the trajectory equation of motion of the mirror is also contained in the shock functions. The exactly solvable mirror case in [57] has shock functions: The W is the product log or W Lambert function, which commonly appears in thermal equilibrium contexts, e.g. Wien's law. 12 One way to get these is as follows: Firstly, one has the simple form u s (v) as it is a simple choice for redshifting ray-tracing f (v) function in the mirror case (or from the spacetime matching solution in the null-shell case). This is given. Secondly, one takes the inverse to get v s (u). While easy, as it turns out, it was unhelpful in obtaining the other shock wave functions. The efficient approach is to notice that u s (v) has a simpler form than v s (u), so one uses u s (v) again to write down t s (v). The inverse of this can be calculated. It is, of course, v s (t). (Note that if one chooses v s (u) to write down t s (u) instead, the inverse is not quite as straight-forward to compute, in fact it is much more complicated.) So, using v s (t), one is set to write down x s (t) = v s (t) − t.
Its inverse is, fortuitously tractable, and gives the above expression for t s (x). We shall interchangeably call the horizon trajectory in the mirror analog case the "black mirror" [60] or "Omex" for short [57]. ("Om"-after the Omega constant, W (1) = Ω = 0.567 where Ωe Ω = 1, and 'ex' after the exponent in W argument.) The acceleration parameter κ in the black mirror case can be identified with the surface gravity in the black hole case, κ = (4M ) −1 , for all times.
The new moving mirror has the following more complicated shock functions: x While these expressions still depend on the primary parameter κ, the intricacy of these expressions arises from the introduction of a second parameter, ξ. Recall that v H in the black mirror case is the location of the horizon. We retain v H for completeness, but make no mistake: the mirror no longer asymptotes to infinite acceleration at a null horizon, located at v H . We shall therefore refer to v H as a "residual horizon". This mirror begins at rest in the far past, and therefore has no initial asymptotic horizon either. The absence of horizons generates the finite total energy, akin to the notion that evaporating black holes exhale only a finite energy flux [70][71][72][73]. 12 The maximum frequency of the (3+1)-dimensional Planck distribution, V i.e. the famous displacement law β ω max = 2.82144.

The Domex Trajectory
The motion of the mirror is given by the trajectory Eq. (22), where v H = 0 for simplicity, and 0 < ξ < 1 is the final speed of the mirror as t → ∞. The motion is initially asymptotically static, lim t→−∞ż (t) = 0, and most notably, the mirror does not approach a future asymptotically static resting state because its future asymptotic speed is lim making this trajectory future asymptotically coasting. The future drifting feature of this mirror means it is an exact model for a remnant [28,30,31] as described by an early anticipation of such solutions by Wilczek in [74]. The trajectory Eq. (24) is plotted in both the spacetime and Penrose diagrams in Figure  (1). The acceleration parameter, κ, is κ > 0, and to be clear, it is not the acceleration of the mirror, α(t) = κ. The rectilinear proper acceleration, α = γ 3z , is time-dependent: The negative sign on Eq.(26) gives a mirror whose motion is to the left. The acceleration has asymptotic behavior such that lim t→±∞ α(t) = 0, making this trajectory asymptotically inertial, despite the drift. As we shall now show, this solution has several analytically tractable results. The special physical aspects of this solution will be investigated in the following sections. We shall refer to this horizonless mirror as Drifting-Omex ("Domex") for short.

The Energy Flux
The energy flux of a moving mirror was first derived by Davies and Fulling [38]. Expressed in terms of the shock functions, it may be computed via Figure 1: Left: In this Penrose diagram, the color curves are "Domex", with asymptotically inertial trajectories. The black curve is "Omex", a horizon mirror (moving mirror with a horizon) [57]. The different coasting speeds correspond to ξ = 0.6, 0.7, 0.8, 0.9, for green, blue, red, purple, respectively. Right: The asymptotically inertial trajectories (Domex) with the same final coasting speeds displayed in the usual spacetime diagram. The dashed lines represent the light cone, and the dotted-dashed horizon line is at v H = 0. The trajectory example here is the same as in the conformal diagram. For comparison, the black line indicates the horizon mirror (Omex) [57], which contains a horizon coinciding with the light cone.

Right Side
The energy flux, emitted to an observer at the right side of the mirror, I + R , using the trajectory of Eq. (24) in Eq. (30), is therefore easily calculated: It contains a build-up phase, a thermal plateau, and an end-phase accompanied by negative energy flux (NEF), see Figure (2). The residual horizon location has been set to v H = 0. A period of thermal emission occurs at extremely high coasting speeds, giving a thermal plateau, which is, in the limit ξ → 1, located for some time ∆t T P , at Interestingly, this is the same as the constant flux produced by the (eternally thermal) Carlitz-Willey trajectory [61]. The Carlitz-Willey mirror radiates a thermal Planckian distribution of particles for all times, at F = F T . In our model, this value occurs because in the limit ξ → 1, this mirror has the same shock functions as the black mirror, which has thermal radiation at late times. However, now it is apparent that in this model, the evaporation eventually stops, effectively decoupling the late-time approximation from the high-frequency approximation. The left observer sees energy flux that strongly peaks the faster the coasting speed of the mirror: ξ = 0.6, 0.7, 0.8, 0.9, colored by green, blue, red, purple, respectively. Notice the initial NEF, and the above-thermal, F (t) > 1, emission. Here κ 2 = 48π. Right: Successive plots of the energy flux observed by the right observer, F (t) ≡ Tuu , Eq.(33), from smallest peak to the largest peak, with varying limiting mirror speeds, ξ = 0.6, 0.7, 0.8, 0.9. Also included are ξ = 1 − 0.1 x where x = 2, 3, 4, 6 colored by orange, pink, cyan, black, respectively. Thermal equilibrium occurs only for a very fast final coasting speed. The total NEF is qualitatively unchanged at this speed. The acceleration parameter is set to κ 2 = 48π so that the plateau levels out at F = 1, the dot-dashed line.
Allowing ξ to be nearly the speed of light, (for example, ξ = 1 − 0.1 10 , with no formal limit), the energy flux emitted to the right observer, has a simple minimal negative value, at some late time, t 0 , where which is a fairly significant proportion of the maximum magnitude amplitude of thermal emission.

Left Side
The energy flux, emitted to an observer at the left of the mirror, I + L , using the trajectory of Eq.(24) by symmetry reversing the sign on ξ, is: The energy flux contains an initial nascent NEF, a rapid reversal to positive energy flux and build-up to a non-thermal positive energy flux peak, and finally a rapid end-phase that falls to zero emission, see Figure (2). It is now clear that while one side of the mirror is approaching thermal equilibrium emission, the other side is experiencing a single nonthermal, ever-more-narrow burst, demonstrating a characteristic difference between the left and right observers. We investigate the pulse via particle spectra in Section (7).

Temperature of Domex
Domex achieves a temperature, 2πT = κ, to lowest order in where ξ ≡ 1− , via a "twice rapid acceleration" (κ(u) = |p /p | = |2η |) approximation. The rapid acceleration, η (u), is identically constant, such that κ(u) = κ, for the eternally thermal mirror (Carlitz-Willey). One finds, For large κu, so long as, κu ≪ −1 , then W (e −κu ) → 0, and to lowest order in , the rapid acceleration is constant, 2πT = |2η | = κ. From the energy flux production, we can help quantify the equilibrium condition of Domex. The simplicity of the time-space function, t s (x), allows for analytic tractability. Finding where the radiation is most near equilibrium amplitude, F ≈ F T ≡ κ 2 /(48π) is possible. Using v H = 0, and t s (x) = −κ −1 e 2κx/ξ − x/ξ, one has the flux as a function of space: Maximizing F (x, ξ, κ) with respect to x, gives the spatial location, x 0 , where the flux is maximum, F (x 0 , ξ, κ) = F max (ξ, κ). Since drift speed is high, then to lowest order in , ignoring the imaginary component of this spatial locus, the real location is The maximum flux, to lowest order in , is then Following Davies [39], Eq. 3.10, or Walker [75], Eq. 5.10, we consider the property that the energy flux of a thermal trajectory has Temperature can be expressed as, where we have taken the positive root, T > 0. To low order in , the result is The lowest order dependence on drift speed scales as ∼ (1 − ξ) 2/3 , indicating, e.g. that speeds of ξ = 1 − 0.1 9 give a millionth part deviation from equilibrium temperature. To ensure Domex is very near equilibrium for an extended period of time, we will use speeds far faster while investigating the time dependence of particle production in Section (7).

Total Energy Produced by Domex
It proves possible to calculate the finite total emitted energy analytically. In terms of the shock wave functions, the total energy to the right side of the mirror, is computed via, or after integration by parts, where the boundary term is ignored due to asymptotic inertial character, where v s ≡ v s (u), u s ≡ u s (v), x s ≡ x s (t) and t s ≡ t s (x). The primes always mean derivatives with respect to the respective function variable. For the mirror trajectory here with finite energy production we shall use the x s (t) integral over dt and confirm it with quanta summing of particles in Section (7), where the total emitted energy is E = ∞ 0 dω ω ∞ 0 dω |β ωω | 2 . Note that by "total", we mean that it is the total amount of energy that only the observer on the right side detects. The mirror emits energy on both sides to two separate observers: left and right. To find the energy emitted to the left, by symmetry, one can simply reverse the motion and compute the energy on the right side again.

Right Side
The total energy radiated to I + R is therefore: Domex does not result in the emission of infinite energy to the usual observer at I + R . Note that the solution here is monotonic for increasing coasting speed and never negative for 0 < ξ < 1. Here, the lim ξ→0 E R = 0, and the lim ξ→1 E R = +∞.

Left Side
For an observer to the left at I + L the total energy emitted is found by simply substituting, ξ → −ξ, into E R , Again, the energy is finite as long as the speed is less than the speed of light. The expression, E L (ξ) is a monotonic function of ξ.

Both Sides
For the high coasting speeds we are interested in, the energy emitted to the left is always much greater than the energy emitted to the right, E L E R . For small values of ξ one finds indicating that E L > E R . As it turns out, E L > E R for all values of the final drift speed, 0 < ξ < 1. The total energy emitted to both observers is E T = E L + E R : where η = tanh −1 ξ, is the final rapidity. See Figure (3) to see a graph of the combined total emitted energy from both sides of the mirror. Notice the divergent behavior as the coasting speed approaches the speed of light. The energy increases monotonically as a function of the coasting speed.

The Entropy of Domex
The von Neumann entanglement entropy in the unitary moving mirror case, can be found from Bianchi-Smerlak's formula 13 [24] as a function of v s (u) ≡ p(u), In terms of the mirror trajectory (see ), this is where z(t) is the trajectory motion of the moving mirror, the dot represents the time derivative and η(t) ≡ tanh −1 [ż(t)] is the time-dependent rapidity. It is simple to see that the faster the mirror moves, the greater the entropy. Unitarity in this sense strictly means that the entropy must achieve a constant value in the far past and far future. The mirror we consider has a non-zero asymptotic entropy, underscoring that eternally thermal radiation is not the end state. Since evolution from a possible initially pure state, to a final mixed state does not occur, radiation of energy flux stops and it is possible to reestablish a possible initial pure state, preserving unitarity. This mirror's trajectory is found to be consistent with the Bianchi-Smerlak [24] entropy-energy relationship where the entropy for Domex, expressed as a function of time, t, for both the right and left sides, approaches a constant value, as the energy flux approaches zero. Notice the independence of κ or the position of the horizon v H for the final entropy value. Here η is the final coasting rapidity (ξ is the final coasting speed). These results are in dramatic contrast to Omex [57] which has infinite total energy and a divergent entanglement entropy in the far future, t → +∞. A plot of the entropy (on both sides of Domex), Eq. (59), is in Figure (4). For Domex, the entanglement entropy has no turn-over of the Page curve [18][19][20] as it is monotonically increasingly approaching some constant value. Although outside the scope of this work, it is worth mentioning that the radiation in an analogous black hole context would never get purified. The preservation of unitarity is maintained only in the sense that the pure state remains pure taking into account both the exterior and interior [30] of the black hole remnant. There are unresolved challenges associated with this type of scenario stemming from the known infinite production problem, stability, and energy conservation. Detailed discussions are provided in [31]. Note however, that Domex, as a possible candidate model for a remnant, is not in conflict with Wilczek redshifting [74] and non-monotonic mass loss of Bianchi-Smerlak [23], as the former occurs even when Domex is once again inertial and in the later, the non-monotonic mass loss occurs as long as the entropy approaches a constant (not necessarily zero). We note that the entanglement entropy of the radiation on the left side of the mirror is negative. This is not unprecedented in the literature, however its interpretation is somewhat unclear. Bianchi, De Lorenzo, and Smerlak interpreted a negative entanglement entropy as the result of the radiation being less correlated than the vacuum [25]. On the other hand, Holzhey, Larsen, and Wilczek were of the opinion that negative renormalized entropy is the result of the radiation having more correlation than the vacuum [29]. We leave open the interpretation of our results, but it is interesting to note that for any fixed ξ, the positive entanglement entropy on the right side of the mirror exactly cancels the negative entanglement entropy on the left side. Therefore, the entanglement entropy of the radiation on the entire slice of any constant t is zero.

Correlation Functions
An under-appreciated lesson stressed by Ford and Roman [76], is that there is a great deal more happening in the accelerating mirror geometry than is revealed by the expectation value of the stress-energy tensor alone. There are subtle increases or reductions in correlations between the flux along rays even where the expectation value vanishes. The stress-energy tensor correlation function is of interest in our situation because it reveals information about the energy flux that demonstrates the thermal character of the radiation above and beyond that of the thermal plateau of the stress-energy tensor during equilibrium. The shock functions for the moving mirror are needed to compute the correlation functions for the stress-energy tensor. It was previously shown that the ray-tracing function p(u) is useful for delta-function pulse piece-wise mirror trajectories [76]. In this section, we extend this work to continuous trajectories and compute the correlations with an emphasis on the equilibrium period of Domex. The energy fluxes emitted by any moving mirror can be positive and negative, but they are only average values. The fluctuations around this average value are generally expected because the quantum state is not an eigenstate of the stress-energy tensor operator.
The correlation function for the stress-energy tensor is where the spacetime points are indicated by y = (u, v) and y = (u , v ). The correlation functions between two right moving rays, two left moving rays, and right and left moving rays are, respectively: Solved in terms of the ray tracing function, p(u), the results are [76] where p (u) = dp(u)/ du and p (u ) = dp(u )/ du . The above expressions deal only with correlations of distinct rays. These expressions simplify, as would be expected, in vacuum or with a static mirror present. For a static mirror we have the condition, v = p(u) = u, and In vacuum C LR (v, u ) = 0 because there can only be correlations with left and right moving fluxes with a mirror present. The correlation limits for C RR and C LL hold for either vacuum or a static mirror, hence the xor, ⊕, in the subscript. The ratios and can tell us about enhancement and suppression of correlations. For R i > 1 one interprets enhancement, for R i < 1 there is suppression.

Correlation Solutions
We focus on one ratio only in order to help confirm thermal equilibrium to the right observer. This correlation ratio, R 1 , involving C RR , associates two right moving rays. These rays come off the mirror heading to the right observer at I + R . The R 1 solutions for the ratios for the three mirrors of interest, which, illustrates thermal correlations at all times. However, for the black mirror, we have which, illustrates thermal correlations at late times. Finally, we write down the Domex's ratio, All three mirrors give the same thermal correlations when comparing a ray that occurs in the appropriate equilibrium period. This is not at very late times for Domex. When one picks a very late time ray, then Omex and Carlitz-Willey are still thermal, but Domex begins to break pattern to abide by the inevitable non-equilibrium completion of emission. This is to be expected because at very late times the mirror abandons thermal character as the radiation ceases. See Figure ( 7 The Particle Production of Domex

The Beta Bogoliubov Coefficient Integrals
The distinction between energy flux and particle flux has been well-studied by Walker-Davies [62] and Walker [75]. The key ingredient is the information of the mirror trajectory equation of motion, which is encapsulated in the shock functions. One needs the Bogoliubov coefficients because the particle emission detected by an observer is There are four ways to calculate the beta Bogoliubov coefficient using the shock functions: Here ω p ≡ ω + ω and ω n ≡ ω − ω . The integration bounds also assume all light rays hit the mirror and propagate to future null infinity on the right. For light rays that do not hit the mirror, one must stop short the integration and add up only to the last null ray as described by the relevant variable. All of the bounds are written for a mirror which starts at positive spatial infinity and proceeds left to negative spatial infinity in the far future. The bounds must be appropriately changed for mirrors which have different behaviors and/or horizons. For Domex we use Eq. (80) because of the simplicity of t s (x) function particular to Domex. Notice the negative sign analogous to Eq. (78) drops away because the mirror starts at x → +∞ at t → −∞. This integral can be used to obtain other trajectories when the other integrals are intractable. We choose this integral because it will allow investigation of particle production while avoiding integration with the product log, (see [63] for z(t) approach), and therefore insert Domex's trajectory, Eq. (24), into the integral Eq. (80) identifying, t s (x) = t(z). The solution is This exact beta solution also contains the thermal plateau, similar to the instantaneous energy flux which closely approaches the thermal line at F = F T . Both particle and energy flux approach the thermal plateau only at high final speeds making it a salient feature of the radiation. One should be confident the behavior shows up in both the particle production and energy flux because the packetized particles carry signatures of the instantaneous energy flux emission [28].
In Figure (6), we construct localized beta Bogoliubov coefficients from the global beta Bogoliubov coefficients in Eq. (81), These are the usual orthonormal complete wave packets [1] which are used to find the time-frequency localized particle count, Particles arrive at I + R in the range of frequencies j ω (j + 1) and in the range of times (2πn − π)/ u (2πn + π)/ . Details on how to construct these packets in more general situations can be found in [63].

Global Particle Distribution
The horizonless beta solution Eq. (81) of Domex has distribution, while the horizon beta Bogoliubov coefficients of Omex have particle count per modesquared, It is easy to see that as ξ → 1 that the spectra coincide. The high frequency approximation, ω ω, applied to Omex's distribution, Eq. (85), gives the usual thermal result of Carlitz-Willey and late-time Davies-Fulling trajectories, One checks that the Domex distribution is thermal, |β ωω ξ | 2 ≈ |β ωω | 2 T , by a series expansion that first approximates the distribution with a very fast end-state drifting mirror ξ ≈ 1 and then applies the high-frequency approximation ω ω. This approach explicitly decouples the late-time approximation from the high-frequency approximation.

Consistency Check
It is a fair claim that the total summation of the energies of each quanta should be equal to the integral over the energy flux: Therefore, the beta Bogoliubov particle results can be confirmed to be consistent with the stress-energy by computing the total energy, using Eq. (76) and Eq. (46), This consistency helps confirm the particles do indeed carry the energy. Explicitly, we have numerically confirmed the beta coefficients, Eq. (81), that for an observer to the right at I + R , the total energy emitted is and for an observer to the left at I + L , the total energy emitted is in agreement with the analytical results of the stress-energy tensor of Section (4). We have found the largest relative numerical error here to be less than 10 −11 using κ = 1, and various values of 0 < ξ < 1.

Discussions
In the Introduction we raised the following question: "Is there a moving mirror in (1+1)dimensions, satisfying unitarity in the sense allowed by the Bianchi-Smerlak criterion (namely, S(u) → const. as u → ±∞), that has no acceleration horizon, produces finite amount of energy, and serves as a limiting case analog for Eddington-Finkelstein coordinate null shell gravitational collapse?" We have answered this question in the affirmative, by constructing an exact mirror solution that satisfies these properties. Furthermore, we investigated both sides of the mirror trajectory, and found interesting features regarding negative entropy and negative energy flux. The hallmark trait of the solution is the fact that it is an asymptotically coasting mirror which does not have an accelerating horizon, yet approaches arbitrarily close to thermal equilibrium. Thermal radiation arises from a sufficiently fast final drifting speed. The ray-tracing function is identical to the spacetime matching condition of the black hole case in the limit that the mirror drifts to the speed of light.
The global approach to treating horizons tends to work well in fully equilibrium thermodynamics, especially so with a priori non-dynamical assumptions (i.e. constant energy flux) [61]. It is well-known that the non-equilibrium cases are not so easy to formulate using the traditional methods [66]. A practical outcome of this paper has been to show how robust the traditional methods can be when the horizon is removed from the start. Non-equilibrium dynamical conditions follow suit, however the system can still achieve equilibrium, for an arbitrary extended amount of time. With a consistent scaling (κ in Domex and Omex are the same scaled parameter relative to thermal emission), we have explicitly used the global geometric properties of the spacetime, and in the case of particle creation, only localized after solving for the global beta Bogoliubov coefficients.
The new mirror was described in terms of its energy flux, total energy, entropy flux, correlations, and particle flux. The temperature can be detected by asymptotic observers with particle detectors (the radiation demonstrates a Planckian distribution for a very fast final drift speed and the use of the high frequency approximation). The evidence for thermality is strengthened further by the long-lived steady-state stress-energy tensor and the correlations which match the eternally thermal equilibrium mirror of Carlitz-Willey [61]. However, a few remarks are in order comparing our analysis to Carlitz-Willey's seperate 1987 trajectory [80] which is not eternally thermal. A notable similarity between this apparent horizon trajectory and the one presented here is the constant rate of particle emission during a finite period of time. We confirm the locally thermal state in both unitary moving mirror trajectories. One notable difference is that their trajectory is in terms of an approximate ray-tracing function with a kink. 14 The trajectory in this paper is exact, C ∞ , for all-times, and is expressed as an explicit space-time trajectory function, Eq. (24). A consequence of this fact is that it so happens our trajectory does not come to rest at late times, while their trajectory requires the mirror eventually become stationary and consequently the entire remnant mass is radiated away to leave behind a flat region of spacetime. In other words, Carlitz-Willey consider a "meta-stable" or "long-lived" remnant that slowly evaporates away, whereas our remnant is eternal.
Unlike the eternal thermal Carltiz-Willey mirror [61] or the black mirror [57], the Domex mirror gives rise to negative energy flux, and by the result of Bianchi-Smerlak [23,24], also to the non-monotonic mass loss of the any corresponding black hole. Current efforts are being directed to explore the generalization of the tortoise coordinate from which can be evident from the generalization of the shock function of Eq. (17) to Eq. (21), in the null-shell case which matches spacetimes outside and inside the shell. For details on the null shell case, see e.g. Unruh (1976) [77], Massar (1996) [78] or Fabbri (2005) [79]. The generalization Eq. (92) and a possible coupling between the parameters ξ and M , may provide clues to understanding any corresponding black hole solution of Domex, and by necessity a different all-time collapse scenario. It is understood [1] that at very early times of gravitational collapse, the system cannot be described by the no-hair theorem. Therefore it is appropriate to consider the type of modifications that can provide various early time approaches to a thermal distribution, particularly those modifications that can afford unitarity and finite evaporation energy. The modifications that can take into account energy conservation like those of the dilaton gravity models have had significant success as a laboratory for studying black hole evaporation. The physical problem in 1+1 dilaton gravity of the evaporating black hole and its modified emission extends to complete evaporation for the Russo, Susskind, and Thorlacius (RST) model [81] and to partial evaporation leaving a remnant for the Bose, Parker, and Peleg (BPP) model [82]. The similarity of the Domex mirror to the BPP model is striking in several qualitative aspects: NEF emission as a thunderpop, a left over remnant, and finite total energy emission. It is also interesting that the mass of the remnant in the BPP model is independent of the mass M of the infalling matter, since with respect to the issue of energy conservation, there is no known physical analog for M = 1/(4κ), the initial mass of the shockwave, in the mirror model. Finally, we shall comment on the peculiar emission we find on the left side of the mirror trajectory, in particular, its possible relevance to the black hole correspondence (if one exists). We conjecture that the "left emission" corresponds to in-falling flux into the black hole. As such, its associated temperature and entropy might shed some light on the information paradox, since as we emphasized in the Introduction, unitarity is a property of the Hilbert space defined on the entire spacetime. In fact, such a "left temperature" in the context of black hole physics already exists in the literature, see e.g., [83]. It should also be emphasized that the recent result in the literature [84], concerning the study of two-dimensional model of gravitational collapse, shows that a geodesic observer on the left side measures late time thermal radiation but zero flux. This result is drastically different from ours. It might be interesting to conduct a comparative study between our model with that of [84].
Ultimately, while Domex is elementary, it embraces several surprisingly interesting traits. Since some of these traits are shared with more sophisticated systems, this solution may be a precursor for ensuing developments (the overt example being curved spacetime collapse). On the other hand, this solution has exposed several explicit general attributes which are unanticipated and must be understood in order to claim a good grasp of the dynamics of the particle creation effect in non-thermal equilibrium.
The outstanding advantage of this mirror solution is the exact expressions for quantities of interest. Since one natural speculation is the direct applicability to a curved spacetime analog, we aim to examine this pertinent and interesting follow-up topic in a later manuscript with primary consideration to energy conservation of the black hole's modified evaporation emission, metric continuity across the shock boundary, and to the Bogolubov coefficients of specific dilaton gravity models.