This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Black-Body Radiation in an Accelerated Frame

Seramika Ariwahjoedi1,2 [email protected]    Apriadi Salim Adam2 [email protected]    Hadyan Luthfan Prihadi3,4 [email protected]    Freddy Permana Zen3,4 [email protected] 1Asia Pacific Center for Theoretical Physics,
Pohang University of Science and Technology, Pohang 37673, Gyeongsangbuk-do, South Korea,
2Research Center for Quantum Physics,
National Research and Innovation Agency (BRIN), South Tangerang 15314, Banten, Indonesia,
3Theoretical High Energy Physics Group, Department of Physics,
FMIPA, Institut Teknologi Bandung, Jl. Ganesha 10 Bandung 40132, West Java, Indonesia,
4Indonesia Center for Theoretical and Mathematical Physics (ICTMP), Indonesia. 
(December 30, 2024)
Abstract

We derive Planck’s radiation law in a uniformly accelerated frame expressed in Rindler coordinates. The black-body spectrum is time-dependent by its temperature and Planckian at each instantaneous time, but it is scaled by an emissivity factor that depends on the Rindler spatial coordinate and the acceleration magnitude. The observer in an accelerated frame will perceive the black-body as black, hyperblack, or grey, depending on its position with respect to the source (moving away or towards), the acceleration magnitude, and the case of whether it is accelerated or decelerated. For an observer accelerating away from the source, there exists a threshold on the acceleration magnitude beyond which it stops receiving radiation from the black-body. Since the frequency and the number of modes in Planck’s law evolve over time, the spectrum is continuously red or blue-shifted towards lower (or higher) frequencies as time progresses, and the radiation modes (photons) could be created or annihilated, depending on the observer’s position and its acceleration or deceleration relative to the source of radiation.

preprint: APS/123-QED

I Introduction

The status of temperature in relativistic thermodynamics remained a subject of ongoing debate, marked by the lack of a universal agreement in the physics community regarding the transformation of temperature for moving bodies (Farias, ; Derakshani, ). In their works, both Planck and Einstein independently derived that the temperature of a moving body would transform as T=TγT^{\prime}=\frac{T}{\gamma}, with γ\gamma representing the Lorentz factor associated with the velocity of the moving body (or moving observer, if the body is at rest) (Planck, ; Einstein, ; Einstein2, ). However, under different assumptions, Ott and Arzelie`\mathrm{\grave{e}}s derived the transformation of temperature to be T=γTT^{\prime}=\gamma T (Ott, ; Arzelies, ). Furthermore, Landsberg proposed the invariance of temperature under transformations between inertial frames, namely, T=TT=T^{\prime} (Landsberg, ; Lansberg2, ). Each of these different formulation for temperature transformations in inertial frames is supported by compelling arguments, hence, although the most widely-accepted theory is the one by Einstein-Planck (Einstein, ; Planck, ), the attempt to reconcile these different approaches remains an ongoing research. A recent advancement in the field of special relativistic thermodynamics was achieved by Nakamura in (Nakamura, ; Nakamura2, ) who revised the Israel-van Kampfen covariant, inverse, 4-temperature (vanKampen, ; Israel, ). This revision allows the derivation of the three distinct temperature transformations within the Israel-van Kampen formulation, depending on the definition of the 3-dimensional volume and the chosen decomposition of the 4-momentum (Nakamura, ; Nakamura2, ).

Moreover, Tolman pioneered the study of general relativistic thermodynamics in (Tolman, ; Tolman2, ), by proposing a law that in a gravitational field, the temperature of a system in thermal equilibrium is inversely proportional to the square root of the gravitational potential at that location. This law can be expressed as T=TgttT^{\prime}=T\sqrt{g_{tt}}, where gttg_{tt} represents the purely-time component of the metric tensor gμνg_{\mu\nu}. In general relativity, the metric tensor describes the gravitational field as the curvature of spacetime, and the gttg_{tt} component corresponds to the gravitational potential in the Newtonian limit (Einstein-1, ; Carrol, ).

On the other hand, the concept of temperature can be defined from several theoretical framework. In thermodynamics, temperature is defined as the parameter shared by two bodies in thermal equilibrium, a state characterized by the absence of net heat transfer between the bodies (Salinger, ); this definition is widely accepted among 19th19^{\mathrm{th}}-century physicists including Clausius, Carnot, Kelvin, and many others. From the statistical mechanics approach, it is defined as the average of kinetic energy of a many-body system, and this is described by the relation EκBT\left\langle E\right\rangle\sim\kappa_{B}T, with κB\kappa_{B} is the Boltzmann constant (Boltzmann, ). Statistical mechanics provides a procedure for bridging macroscopic observables of a system, such as temperature and entropy, to its microscopic properties, i.e., the position, velocity, and kinetic energy of the individual degrees of freedom (Liboff, ). Combining these two perspectives, temperature can be defined as an observable that emerge from the statistical properties of individual degrees of freedom (i.e., the kinetic energy) of a many-body system that is in a thermal equilibrium with its surrounding. As far as we are concern, this is the current status of the temperature conceptualization widely-accepted by physics community at the moment.

In our study, we adopt a more pragmatic approach to the concept of temperature -specifically, through the use of Planck’s law on black-body radiation, originally proposed in 1900 (Planck2, ). This law initiated the birth of quantum mechanics and remains valid through experiments to the present day, and it founds wide applications in pyrometric measurements. As explained in the preceding paragraphs, the thermodynamic definition of temperature -a parameter indicating two systems are in thermodynamic equilibrium- is particularly applicable for temperature measurement in a rest frame. The measurement apparatus must establish direct contact with the system and afford sufficient time for thermalization until a (thermal) equilibrium is reached. In this case, the temperature is well-defined. However, for these two systems to be in thermal equilibrium, in prior, they need to be in a mechanical equilibrium, while maintaining a direct contact to undergo thermalization. This requirement is challenging to fulfill if one of the systems is not at rest relative to the other (vanKampen, ; vanKampen3, ). In principle, at least by the definition originating from thermodynamical approach, the notion of temperature can only emerge as a consequence of local measurement (Mares, ). Furthermore, a quantum-field theory calculation, using the Unruh-DeWitt detector as a thermometer moving through a thermal bath, gives a result that the particle distribution is non-Planckian, which made it difficult to define temperature in this context (Costa, ). Based on this work, then Landsberg and Matsas proposed an argument on the impossibility of a universal relativistic temperature transformation, due to the fact that there is no continuous function that could map the non-Planckian distribution from (Costa, ) to a Planckian one (Matsas, ; Matsas2, ). These reasons supports the argument that the concept of temperature of a moving body is not well-defined (Mares, ).

However, using the Planckian spectrum, one can predict the temperature of a distant object without necessitating direct contact, thereby skipping the need to reach thermal equilibrium with the object under observation. This methodology, treating temperature as a derived parameter, finds widespread application – from the infrared thermometers to the calculation of the temperatures of celestial bodies such as stars and the cosmic microwave background (Bbracewell, ). The process involves the collection of the complete radiation spectrum emitted by the object and, -under the assumption that the object behaves as a black-body-, the comparation of the spectrum data with the theoretical black-body spectrum. Temperature prediction can be carried out using either Wien’s displacement law (Wien, ; Wien2, ) or the Stefan-Boltzmann law (Stefan, ). In the case of a black-body source, these two distinct temperature calculation methods coincide. However, it is important to note that assuming all objects behave as perfect black-bodies is overly restrictive. For a more realistic approach, the collected data should be compared to a calibrated grey-body spectrum derived from the black-body spectrum, but with emissivity or absorptivity values set below 1. This method, to the best of our knowledge, represents the sole viable approach for ’measuring’ the temperature of a moving object. Furthermore, experiments based on this measurement technique are relatively straightforward to conduct. These facts emphasize the importance of Planck’s law in the determination of temperatures of moving bodies.

We derived Planck’s law and calculated the black-body spectrum in a uniformly accelerated frame expressed in Rindler coordinates. The reason behind our quest on this subject can be explained as follows: In the framework of general relativity, an inertial frame is an idealization, hence one needs to move to a more general reference of frame. As a first realistic step towards achieving this objective, we consider a uniformly accelerated frame, expressed in Rindler coordinates. The calculation of the black-body spectrum in the inertial frame had been done in (Peebles, ; Heer, ; Henry, ). These works demonstrated that the Planckian spectrum is invariant under Lorentz transformation, and the effective temperature of the moving object depends on its velocity and the direction of observation (Peebles, ; Heer, ; Henry, ). Another different approach shows only the zero-point temperature term is invariant under Lorentz transformation (Bouyer, ). Attempts to calculate the black-body spectrum in an accelerated frame had been done in (Lee, ), but to the best of our knowledge, the calculation of the spectrum in Rindler coordinates is still lacking in the existing literature. In the present work, we follow the (effective) ’directional-temperature’ approach as in (Henry, ), but instead of using the temperature TT in Planck’s law, we substitute it with variables that maximize the energy in the spectrum (in our case, we use the ’maximal’ frequency). With this substitution, we avoid the problem associated with temperature transformation, as we are certain about the Lorentz transformation for the ’maximal’ variables.

The main result presented in this article is the black-body spectrum in an accelerated frame. The spectrum has an explicit time dependency on its temperature, Planckian at each instantaneous time, but it is scaled by a factor dependent on the spatial coordinate and the acceleration magnitude. The spatial, coordinate-dependent scale factor is proportional to eαξ\sim e^{\alpha\xi}, while the scale factor related to the acceleration is proportional to eαξ(1±α)\sim e^{\alpha\xi}\left(1\pm\alpha\right), where the ±\pm sign depending on the observer is either deccelerated or accelerated. This scale factor could be physically interpreted as the emissivity factor of the source, hence, in an accelerated frame, a black-body could be perceived as grey or ’hyper-black’, depending on the magnitude of the acceleration. Furthermore, the variables in Planck’s law, specifically the number of modes and frequency, evolve over time. However, since we are considering all the possible value of frequencies from zero to infinity, the time dependence of the frequencies does not explicitly shown in the law. The Planckian spectrum is continuously red or blue-shifted towards lower (or higher) frequencies as time progresses. In the accelerated frame, we demonstrate that the radiation modes (photons) could be positive or negative, depending on the acceleration or deceleration of the observer, and zero for vanishing acceleration. In the end, assuming the validity of Wien’s displacement law in an accelerated frame, the time-dependent, (directional)-temperature of a body in an accelerated frame is given by T¯[τ]=Teαξ(coshς¯[τ]+cosθ¯[τ]sinhς¯[τ])1\bar{T}\left[\tau\right]=Te^{-\alpha\xi}\left(\cosh\bar{\varsigma}\left[\tau\right]+\cos\bar{\theta}\left[\tau\right]\sinh\bar{\varsigma}\left[\tau\right]\right)^{-1}. It is worth noting that the spatial, coordinate-dependent scale factor eαξ\sim e^{-\alpha\xi} resembles the temperature scale factor in Tolman’s theory of general relativistic thermodynamics (Tolman, ), and for α=0,\alpha=0, the transformation of the temperature returns to the ones obtained in (Peebles, ; Heer, ; Henry, ).

The paper is structured as follows: Section II contains the derivation of black-body radiation in an accelerated frame. First, we introduce the properties of the Rindler space, followed by the transformation from Minkowski to Rindler coordinates for each variable contained in Planck’s law. With these transformed variables, we derive Planck’s distribution law in the accelerated frame. Section III contains discussions and the physical interpretation of the results, including the aberration of light, relativistic beaming, and red/blueshifts due to the relativistic Doppler effect, creation/annihilation of modes in the accelerated/decelerated frame, the emissivity factor of the black-body source, and the transformation of temperature, assuming Wien’s law holds in the accelerated frame. Finally, in Section IV, we provide a conclusive summary of our work and outline some insights into potential further research. We have included appendices to ensure the self-contained nature of this paper. The appendices contain the derivation of the black-body spectrum in a frame moving with a constant velocity. This supplementary content aids readers in understanding each step of the spectrum derivation in the accelerated frame. One crucial aspect is that we have aimed to keep this paper as classical as possible. The only quantum assumptions we adopt are the explicit use of the de Broglie postulate p3=k\,{}^{3}\mathrm{p}=\hbar\mathrm{k} for photon momentum (with k\mathrm{k} is the photon wave-vector and \hbar is the Planck constant) and the assumption of energy discreteness E=nωE=n\hbar\omega (with ω\omega is the photon angular frequency and nn are integers) as the requirement to derive Planck’s law.

II Blackbody Radiation in Accelerated Frame

II.1 Rindler Coordinates

Let (,𝐱)\left(\mathcal{M},\mathbf{x}\right) be a Minkowski space equipped with a metric η\eta. The (global) coordinate 𝐱\mathbf{x} =(x,y,z,t)\left(x,y,z,t\right), with t:=c𝔱t:=c\mathfrak{t} is the time coordinate and (x,y,z)\left(x,y,z\right) is the spatial part, is used to parametrize \mathcal{M}. In this coordinate, the metric η\eta is a diagonal metric with signature (1,1,1,1)\left(-1,1,1,1\right). Naturally, one could attach an inertial reference of frame to the coordinate 𝐱\mathbf{x}, and let 𝒪\mathcal{O} be an observer at rest with respect to this inertial frame. Let us call this frame/observer as 𝒪.\mathcal{O}.

Suppose we have another observer 𝒜\mathcal{A} that moves with a constant acceleration of magnitude aa in the direction of the xx-axis, with respect to 𝒪\mathcal{O}. The trajectory of 𝒜\mathcal{A} in \mathcal{M} is 𝐱[λ]=xμ[λ]\mathbf{x}\left[\lambda\right]=x^{\mu}\left[\lambda\right], with λ\lambda is a parameter (a proper time with respect to 𝒜\mathcal{A}), satisfying:

t[λ]\displaystyle t\left[\lambda\right] =1asinhaλ,\displaystyle=\frac{1}{a}\sinh a\lambda, (1)
x[λ]\displaystyle x\left[\lambda\right] =1acoshaλ,\displaystyle=\frac{1}{a}\cosh a\lambda, (2)

and y,zy,z being constant. The 4-acceleration is aμ=d2xμdλ2,a^{\mu}=\frac{d^{2}x^{\mu}}{d\lambda^{2}}, where the non-zero components are:

at\displaystyle a^{t} =asinhaλ,\displaystyle=a\sinh a\lambda,
ax\displaystyle a^{x} =acoshaλ,\displaystyle=a\cosh a\lambda,

such that the magnitude of aμa^{\mu} is aμaμ=a\sqrt{a^{\mu}a_{\mu}}=a. Note here that aa is positive definite: a0a\geq 0. From (1) and (2), the trajectory of 𝒜\mathcal{A} satisfies x2t2=a2x^{2}-t^{2}=a^{-2}, describing a collection of hyperboloids with asymptotes at the lines x=tx=-t and x=tx=t. The asymptotes divide \mathcal{M} into 4 regions, and let us focus on one of the region where 0<x<0<x<\infty and x<t<x-x<t<x, usually labeled as the (right) Rindler wedge, see FIG. 1.

Refer to caption
Figure 1: The left and right Rindler wedge. The hyperboloids describe trajectories of constant ξ\xi, while the straight lines of varying slopes describe trajectories of constant τ\tau. The point (x,t)=(0,0)\left(x,t\right)=\left(0,0\right) in the Cartesian coordinates coincides with the points (τ,ξ)=(0,)\left(\tau,\xi\right)=\left(0,-\infty\right) in the Rindler coordinate. The (Killing) horizons are located at x=±tx=\pm t.

One could choose a new coordinate patch 𝐱¯\bar{\mathbf{x}} =(ξ,y¯,z¯,τ)\left(\xi,\bar{y},\bar{z},\tau\right), to parametrize the Rindler wedge as follows:

τ=12α(lnα(x+t)lnα(xt)2ς),ξ=12α(lnα(x+t)+lnα(xt)),\begin{array}[]{cc}\tau=&\frac{1}{2\alpha}\left(\ln\alpha\left(x+t\right)-\ln\alpha\left(x-t\right)-2\varsigma\right),\\ \xi=&\frac{1}{2\alpha}\left(\ln\alpha\left(x+t\right)+\ln\alpha\left(x-t\right)\right),\end{array} (3)

and their inverse:

t=1αeαξsinh(ατ+ς),x=1αeαξcosh(ατ+ς),\begin{array}[]{cc}t&=\frac{1}{\alpha}e^{\alpha\xi}\sinh\left(\alpha\tau+\varsigma\right),\\ x&=\frac{1}{\alpha}e^{\alpha\xi}\cosh\left(\alpha\tau+\varsigma\right),\end{array} (4)

with y¯=y\bar{y}=y , z¯=z\bar{z}=z, and α,\alpha, ς\varsigma are constant parameters. The (right) Rindler wedge is covered by the coordinate patch with <ξ,τ<-\infty<\xi,\tau<\infty. For ς=0\varsigma=0, the coordinate is usually known as the Rindler coordinates in the Lass/radar representation (Carrol, ). In this coordinate, (1)-(2) becomes:

τ[λ]\displaystyle\tau\left[\lambda\right] =1α(aλς),\displaystyle=\frac{1}{\alpha}\left(a\lambda-\varsigma\right),
ξ[λ]\displaystyle\xi\left[\lambda\right] =1αlnαa.\displaystyle=\frac{1}{\alpha}\ln\frac{\alpha}{a}.

Notice that along the trajectory 𝐱[λ]\mathbf{x}\left[\lambda\right], ξ\xi is constant and τ\tau is proportional to the proper time λ\lambda, although it is shifted by a constant amount ςα-\frac{\varsigma}{\alpha}. An observer moving along constant ξ\xi will experience an acceleration of magnitude:

a=αeαξ.a=\alpha e^{-\alpha\xi}. (5)

Since aa is positive definite, then so does α\alpha: α0\alpha\geq 0. Naturally, the coordinate 𝐱¯\bar{\mathbf{x}} is adopted by the observer/frame 𝒜\mathcal{A}.

One can obtain the infinitesimal transformation of (3), which is:

dτ\displaystyle d\tau =xdttdxα(x2t2)=eαξ(coshς¯dtsinhς¯dx),\displaystyle=\frac{xdt-tdx}{\alpha\left(x^{2}-t^{2}\right)}=e^{-\alpha\xi}\left(\cosh\bar{\varsigma}dt-\sinh\bar{\varsigma}dx\right), (6)
dξ\displaystyle d\xi =xdxtdtα(x2t2)=eαξ(coshς¯dxsinhς¯dt),\displaystyle=\frac{xdx-tdt}{\alpha\left(x^{2}-t^{2}\right)}=e^{-\alpha\xi}\left(\cosh\bar{\varsigma}dx-\sinh\bar{\varsigma}dt\right),

where we write ς¯=ς¯[τ]=ατ+ς\bar{\varsigma}=\bar{\varsigma}\left[\tau\right]=\alpha\tau+\varsigma, and their inverses:

dt\displaystyle dt =eαξ(coshς¯dτ+sinhς¯dξ)=α(xdτ+tdξ),\displaystyle=e^{\alpha\xi}\left(\cosh\bar{\varsigma}d\tau+\sinh\bar{\varsigma}d\xi\right)=\alpha\left(xd\tau+td\xi\right),
dx\displaystyle dx =eαξ(coshς¯dξ+sinhς¯dτ)=α(xdξ+tdτ),\displaystyle=e^{\alpha\xi}\left(\cosh\bar{\varsigma}d\xi+\sinh\bar{\varsigma}d\tau\right)=\alpha\left(xd\xi+td\tau\right), (7)

with dy¯=dyd\bar{y}=dy and dz¯=dzd\bar{z}=dz. ς¯\bar{\varsigma} is the time-dependent rapidity that is linear to the “proper time” τ\tau. For α=0\alpha=0, the acceleration aa vanishes, and (6)-(7) will return to the standard Lorentz transformation (81)-(82).

The line element in Minkowksi space could be written in Rindler coordinate using (7) as follows:

ds2=gμνdx¯μdx¯ν=e2αξ(dτ2+dξ2)+dy¯2+dz¯2;ds^{2}=g_{\mu\nu}d\bar{x}^{\mu}d\bar{x}^{\nu}=e^{2\alpha\xi}\left(-d\tau^{2}+d\xi^{2}\right)+d\bar{y}^{2}+d\bar{z}^{2}; (8)

The coefficient of the line element gives the Rindler metric gμνg_{\mu\nu} in Lass representation. One could also consider the left Rindler wedge (the region IV in FIG. 1), where it could be covered by a coordinate patch similar to the RRW, but with the flip in the sign of (4) as follows:

t=1αeαξsinh(ατ+ς),x=1αeαξcosh(ατ+ς).\begin{array}[]{cc}t&=-\frac{1}{\alpha}e^{\alpha\xi}\sinh\left(\alpha\tau+\varsigma\right),\\ x&=-\frac{1}{\alpha}e^{\alpha\xi}\cosh\left(\alpha\tau+\varsigma\right).\end{array} (9)

In this region, the future-directed time-like Killing vector is τ-\partial\tau, instead of τ\partial\tau in region I.

II.2 The Aberration of Light in Accelerated Frame

Suppose with respect to an inertial frame 𝒪\mathcal{O} we have an object moving with a constant velocity 𝐯\mathbf{v} as follows:

𝐯=(vx,vy,vz)=(dxdt,dydt,dzdt).\mathbf{v}=\left(v_{x},v_{y},v_{z}\right)=\left(\frac{dx}{dt},\frac{dy}{dt},\frac{dz}{dt}\right).

According to the accelerated frame 𝒜\mathcal{A}, the velocity of the object is:

𝐯¯[τ,ξ]=(vξ,v¯y,v¯z)=(dξdτ,dy¯dτ,dz¯dτ).\bar{\mathbf{v}}\left[\tau,\xi\right]=\left(v_{\xi},\bar{v}_{y},\bar{v}_{z}\right)=\left(\frac{d\xi}{d\tau},\frac{d\bar{y}}{d\tau},\frac{d\bar{z}}{d\tau}\right).

Notice that the components of 𝐯¯[τ,ξ]\bar{\mathbf{v}}\left[\tau,\xi\right] are functions of the Rindler coordinate (τ,ξ)\left(\tau,\xi\right). Using the transformation (6)-(7), one can obtain the relation between the instantaneous 𝐯¯[τ,ξ]\bar{\mathbf{v}}\left[\tau,\xi\right] and 𝐯\mathbf{v} as follows:

vξ=vxtanhς¯1vxtanhς¯=(xvxt)(xtvx),v_{\xi}=\frac{v_{x}-\tanh\bar{\varsigma}}{1-v_{x}\tanh\bar{\varsigma}}=\frac{\left(xv_{x}-t\right)}{\left(x-tv_{x}\right)}, (10)
v¯y,z=eαξvy,zcoshς¯(1vxtanhς¯)=α(x2t2)vy,z(xtvx).\bar{v}_{y,z}=\frac{e^{\alpha\xi}v_{y,z}}{\cosh\bar{\varsigma}\left(1-v_{x}\tanh\bar{\varsigma}\right)}=\frac{\alpha\left(x^{2}-t^{2}\right)v_{y,z}}{\left(x-tv_{x}\right)}. (11)

Relation (10)-(11) could be derived from the standard velocity addition formula (or, similarly, from the relativistic acceleration formula). Their inverses could also be obtained as follows:

vx=vξ+tanhς¯1+vξtanhς¯,v_{x}=\frac{v_{\xi}+\tanh\bar{\varsigma}}{1+v_{\xi}\tanh\bar{\varsigma}},
vy,z=v¯y,zeαξcoshς¯(1+vξtanhς¯).v_{y,z}=\frac{\bar{v}_{y,z}}{e^{\alpha\xi}\cosh\bar{\varsigma}\left(1+v_{\xi}\tanh\bar{\varsigma}\right)}.

In the following paragraphs, we will derive the light aberration formula for an accelerated observer 𝒜\mathcal{A}. For a brief review on the light aberration formula for an inertial, moving observer, one could refer to Appendix B, or consult (Johnson, ). Let us consider a black-body source (photons inside a cavity) at rest with respect to an inertial frame 𝒪.\mathcal{O}. Electromagnetic radiations (photons) are emitted by the black-body with propagation velocity cc in the direction of n^=(cosθ,sinθ,0)\hat{n}=\left(\cos\theta,\sin\theta,0\right), with θ\theta is the (polar) observation (or inclination) angle between the axis xx and n^\hat{n} on plane dxdydx\wedge dy. The velocity of the photon with respect to frame 𝒪\mathcal{O} is:

𝐯=(vx,vy,vz)=cn^=(ccosθ,csinθ,0).\mathbf{v}=\left(v_{x},v_{y},v_{z}\right)=c\hat{n}=\left(c\cos\theta,c\sin\theta,0\right). (12)

The 4-velocity of the photon satisfies the null-vector condition, see Appendix X. To derive the 3-velocity in 𝒜\mathcal{A}, let us first write the 4-velocity in Rindler coordinate:

𝐮=𝐮¯=τλ(1,𝐯¯),\mathbf{u}=\bar{\mathbf{u}}=\frac{\partial\tau}{\partial\lambda}\left(1,\bar{\mathbf{v}}\right),

where λ\lambda is a real parameter (usually the proper time with respect to the moving object) and 𝐯¯\bar{\mathbf{v}} is the 3-velocity in 𝒜\mathcal{A}:

𝐯¯=(vξ,v¯y,v¯z)=(ξτ,y¯τ,z¯τ).\bar{\mathbf{v}}=\left(v_{\xi},\bar{v}_{y},\bar{v}_{z}\right)=\left(\frac{\partial\xi}{\partial\tau},\frac{\partial\bar{y}}{\partial\tau},\frac{\partial\bar{z}}{\partial\tau}\right).

Using the null-vector condition for light, we have:

u¯μu¯μ=gμνu¯μu¯ν=(τλ)2(e2αξ+e2αξvξ2+v¯y2+v¯z2)=0.\bar{u}^{\mu}\bar{u}_{\mu}=g_{\mu\nu}\bar{u}^{\mu}\bar{u}^{\nu}=\left(\frac{\partial\tau}{\partial\lambda}\right)^{2}\left(-e^{2\alpha\xi}+e^{2\alpha\xi}v_{\xi}^{2}+\bar{v}_{y}^{2}+\bar{v}_{z}^{2}\right)=0.

In 𝒪\mathcal{O}, the light only propagates in the direction n^\hat{n} on plane dxdydx\wedge dy, so even in 𝒜\mathcal{A}, the z¯\bar{z} component of the velocity of light is zero. The null condition then becomes:

vξ2+(eαξv¯y)2=1.v_{\xi}^{2}+\left(e^{-\alpha\xi}\bar{v}_{y}\right)^{2}=1.

Writing vξ2=cos2θ¯v_{\xi}^{2}=\cos^{2}\bar{\theta} and (eαξv¯y)2=sin2θ¯,\left(e^{-\alpha\xi}\bar{v}_{y}\right)^{2}=\sin^{2}\bar{\theta}, the 3-velocity of the photon in 𝒜\mathcal{A} is:

𝐯¯=(vξ,v¯y,v¯z)=cn¯^=(ccosθ¯,ceαξsinθ¯,0),\bar{\mathbf{v}}=\left(v_{\xi},\bar{v}_{y},\bar{v}_{z}\right)=c\hat{\bar{n}}=\left(c\cos\bar{\theta},ce^{\alpha\xi}\sin\bar{\theta},0\right), (13)

with n¯^\hat{\bar{n}} and θ¯\bar{\theta} are the propagation direction and the inclination angle of the photon at 𝒜\mathcal{A}, respectively. Notice that since the observer 𝒜\mathcal{A} is accelerated in the direction +x+x, it will receive only photons that are moving towards the observer, namely, the ones that have the +x+x component in the velocity. The relation between 𝐯¯\bar{\mathbf{v}} and 𝐯\mathbf{v} can be obtained using the velocity addition formula in the accelerated frame (10)-(11). Inserting (12)-(13) to (10)-(11) will result in the light aberration for an accelerated frame 𝒜\mathcal{A}:

cosθ¯\displaystyle\cos\bar{\theta} =cosθtanhς¯1cosθtanhς¯,\displaystyle=\frac{\cos\theta-\tanh\bar{\varsigma}}{1-\cos\theta\tanh\bar{\varsigma}}, (14)
sinθ¯\displaystyle\sin\bar{\theta} =sinθcoshς¯(1cosθtanhς¯),\displaystyle=\frac{\sin\theta}{\cosh\bar{\varsigma}\left(1-\cos\theta\tanh\bar{\varsigma}\right)}, (15)
tanθ¯2\displaystyle\tan\frac{\bar{\theta}}{2} =1coshς¯(1tanhς¯)tanθ2,\displaystyle=\frac{1}{\cosh\bar{\varsigma}\left(1-\tanh\bar{\varsigma}\right)}\tan\frac{\theta}{2}, (16)

where the last equation is obtained from the trigonometry identity tanθ2=sinθ1+cosθ\tan\frac{\theta}{2}=\frac{\sin\theta}{1+\cos\theta}. One could also obtained their inverses:

cosθ\displaystyle\cos\theta =cosθ¯+tanhς¯1+cosθ¯tanhς¯,\displaystyle=\frac{\cos\bar{\theta}+\tanh\bar{\varsigma}}{1+\cos\bar{\theta}\tanh\bar{\varsigma}}, (17)
sinθ\displaystyle\sin\theta =sinθ¯coshς¯(1+cosθ¯tanhς¯),\displaystyle=\frac{\sin\bar{\theta}}{\cosh\bar{\varsigma}\left(1+\cos\bar{\theta}\tanh\bar{\varsigma}\right)}, (18)
tanθ2\displaystyle\tan\frac{\theta}{2} =1coshς¯(1+tanhς¯)tanθ¯2.\displaystyle=\frac{1}{\cosh\bar{\varsigma}\left(1+\tanh\bar{\varsigma}\right)}\tan\frac{\bar{\theta}}{2}. (19)

The aberration formula in 𝒜\mathcal{A} has exactly the same form with the ones in a constant moving frame 𝒪\mathcal{O}^{\prime} (see relation (89)-(91)), however, it should be kept in mind that θ¯\bar{\theta} is time-dependent, namely θ¯=θ¯[τ]\bar{\theta}=\bar{\theta}\left[\tau\right], while θ\theta^{\prime} is constant.

II.3 The Doppler Effect in Accelerated Frame

In this paper we will only consider the longitudinal relativistic Doppler effect (LongitudDoppler, ), where the observer/source velocity has component parallel to the wave propagation. To derive this, let us first consider the 4-momentum constraint (107) which is valid for any reference frame. The photon is massless, so using the de Broglie postulate E=ωE=\hbar\omega and 𝐩=𝐤\mathbf{p}=\hbar\mathbf{k}, (107) becomes the dispersion relation:

ω2=|𝐤|2,\omega^{2}=\left|\mathbf{k}\right|^{2}, (20)

for the wave vector of the photon satisfies 𝐤=2πn^λ.\mathbf{k}=\frac{2\pi\hat{n}}{\lambda}. In 𝒪,\mathcal{O}, the dispersion relation becomes:

ω2=kx2+ky2.\omega^{2}=k_{x}^{2}+k_{y}^{2}.

Since the coordinate axis xx is perpendicular to yy, one could write kx=ωcosθk_{x}=\omega\cos\theta and ky=ωsinθk_{y}=\omega\sin\theta, with θ\theta is the inclination angle at 𝒪.\mathcal{O}. Therefore, the 4-momentum of a photon with frequency ω\omega according to the rest frame 𝒪\mathcal{O} (with respect to the black-body source) is 𝐩=(pt,px,py,pz)\mathbf{p}=\left(p_{t},p_{x},p_{y},p_{z}\right), satisfying:

𝐩=ωc(1,n^)=ωc(1,cosθ,sinθ,0),\mathbf{p}=\frac{\hbar\omega}{c}\left(1,\hat{n}\right)=\frac{\hbar\omega}{c}\left(1,\cos\theta,\sin\theta,0\right), (21)

Now, let us obtain the photon’s 4-momentum according to the accelerated frame 𝒜\mathcal{A}, that is 𝐩¯=(pτ,pξ,p¯y,p¯z).\bar{\mathbf{p}}=\left(p_{\tau},p_{\xi},\bar{p}_{y},\bar{p}_{z}\right). The dipersion relation (20) is also satisfied in 𝒜\mathcal{A}, however, written in different coordinate, it becomes:

e2αξω¯2=e2αξkξ2+ky¯2,e^{2\alpha\xi}\bar{\omega}^{2}=e^{2\alpha\xi}k_{\xi}^{2}+k_{\bar{y}}^{2}, (22)

(notice that the dispersion relation comes from pμpμ=0,p_{\mu}p^{\mu}=0, hence the metric components needs to be taken accounted). Since the coordinate axis ξ\xi and y¯\bar{y} are also perpendicular to one another, one could defined that kξ=ω¯cosθ¯k_{\xi}=\bar{\omega}\cos\bar{\theta} and ky¯=ω¯eαξsinθ¯k_{\bar{y}}=\bar{\omega}e^{\alpha\xi}\sin\bar{\theta}. Notice that θ¯\bar{\theta} is the inclination angle according to 𝒜\mathcal{A}. With this, then the 4-momentum at the accelerated frame is:

𝐩¯=ω¯(1,n¯^[ξ,τ])=ω¯(1,cosθ¯,eαξsinθ¯,0).\bar{\mathbf{p}}=\hbar\bar{\omega}\left(1,\hat{\bar{n}}\left[\xi,\tau\right]\right)=\hbar\bar{\omega}\left(1,\cos\bar{\theta},e^{\alpha\xi}\sin\bar{\theta},0\right). (23)

ω¯\bar{\omega} is the frequency of the photon, observed by 𝒜\mathcal{A}. Here, we assume that de Broglie postulate is still valid in 𝒜,\mathcal{A}, namely the relation 𝐩¯[τ,ξ]=𝐤¯[τ,ξ]\bar{\mathbf{p}}\left[\tau,\xi\right]=\hbar\bar{\mathbf{k}}\left[\tau,\xi\right] is satisfied for every position in space and each instantaneous time. One could refer to Appendix B for a more detailed explanation on the photon’s momentum.

It needs to be kept in mind that there exists another solution to (20), namely: ω=k,\omega=-k, giving another possible value for kx=ωcosθk_{x}=-\omega\cos\theta and ky=ωsinθk_{y}=-\omega\sin\theta in 𝒪\mathcal{O} and kξ=ω¯cosθ¯k_{\xi}=-\bar{\omega}\cos\bar{\theta} and ky¯=ω¯eαξsinθ¯k_{\bar{y}}=-\bar{\omega}e^{\alpha\xi}\sin\bar{\theta} in 𝒜\mathcal{A}. However, since we are considering only the right Rindler wedge (RRW), the reasonable case for an observer in RRW to receive the signal is by the waveform that moves to the right (in +x+x direction). For an observer in the left Rindler wedge (LRW, which we will consider in Section IV), we should consider the case where ω=k.\omega=-k.

The (generalized) momentum is naturally a covector, however, in this paper, we use its contravariant counterpart to be consistent with the 4-velocity (12)-(13). The entire results do not depends on the choice of vector/covector for the derivation. Since the 4-momentum is a 4-vector, the transformation between 𝐩\mathbf{p} and 𝐩¯\bar{\mathbf{p}} follows (6)-(7), where we write the cofficient of the transformation in terms of hyperbolic functions:

pτ=eαξ(ptcoshς¯pxsinhς¯),pξ=eαξ(pxcoshς¯ptsinhς¯),\begin{array}[]{cc}p_{\tau}&=e^{-\alpha\xi}\left(p_{t}\cosh\bar{\varsigma}-p_{x}\sinh\bar{\varsigma}\right),\\ p_{\xi}&=e^{-\alpha\xi}\left(p_{x}\cosh\bar{\varsigma}-p_{t}\sinh\bar{\varsigma}\right),\end{array} (24)

with their inverses:

pt=eαξ(pτcoshς¯+pξsinhς¯),px=eαξ(pξcoshς¯+pτsinhς¯),\begin{array}[]{cc}p_{t}&=e^{\alpha\xi}\left(p_{\tau}\cosh\bar{\varsigma}+p_{\xi}\sinh\bar{\varsigma}\right),\\ p_{x}&=e^{\alpha\xi}\left(p_{\xi}\cosh\bar{\varsigma}+p_{\tau}\sinh\bar{\varsigma}\right),\end{array} (25)

p¯y=py\bar{p}_{y}=p_{y} and p¯z=pz.\bar{p}_{z}=p_{z}. Moreover, by inserting (21)-(23) to the transformation (24), we can obtain the 3 equivalent forms of the relativistic Doppler effect in the accelerated frame as follows:

ω¯ω\displaystyle\frac{\bar{\omega}}{\omega} =eαξ(coshς¯cosθsinhς¯),\displaystyle=e^{-\alpha\xi}\left(\cosh\bar{\varsigma}-\cos\theta\sinh\bar{\varsigma}\right), (26)
=eαξ(cosθcoshς¯sinhς¯)cosθ¯,\displaystyle=e^{-\alpha\xi}\frac{\left(\cos\theta\cosh\bar{\varsigma}-\sinh\bar{\varsigma}\right)}{\cos\bar{\theta}}, (27)
=sinθsinθ¯eαξ,\displaystyle=\frac{\sin\theta}{\sin\bar{\theta}}\,e^{-\alpha\xi}, (28)

together with their inverses:

ωω¯\displaystyle\frac{\omega}{\bar{\omega}} =eαξ(coshς¯+cosθ¯sinhς¯),\displaystyle=e^{\alpha\xi}\left(\cosh\bar{\varsigma}+\cos\bar{\theta}\sinh\bar{\varsigma}\right), (29)
=eαξ(cosθ¯coshς¯+sinhς¯)cosθ.\displaystyle=e^{\alpha\xi}\frac{\left(\cos\bar{\theta}\cosh\bar{\varsigma}+\sinh\bar{\varsigma}\right)}{\cos\theta}.

Notice that ω¯=ω¯[τ,ξ]\bar{\omega}=\bar{\omega}\left[\tau,\xi\right] is a function of time and position in the Rindler coordinate. From (26) and (29), one could retrieve the aberration formula (14)-(18).

Let us consider one of the light aberation formula (17). Differentiating (17) with respect to any parameter λ\lambda give:

sinθdθdλ=sinθ¯(dθ¯dλ+αsinθ¯dτdλ)(coshς¯+cosθ¯sinhς¯)2,-\sin\theta\frac{d\theta}{d\lambda}=\frac{\sin\bar{\theta}\left(-\frac{d\bar{\theta}}{d\lambda}+\alpha\sin\bar{\theta}\frac{d\tau}{d\lambda}\right)}{\left(\cosh\bar{\varsigma}+\cos\bar{\theta}\sinh\bar{\varsigma}\right)^{2}},

while using (18), (28), and the chain rule dτdλ=dτdθ¯dθ¯dλ\frac{d\tau}{d\lambda}=\frac{d\tau}{d\bar{\theta}}\frac{d\bar{\theta}}{d\lambda}, gives:

dθdθ¯=ω¯ωeαξ(1αθ¯˙sinθ¯),\frac{d\theta}{d\bar{\theta}}=\frac{\bar{\omega}}{\omega}e^{\alpha\xi}\left(1-\frac{\alpha}{\dot{\bar{\theta}}}\sin\bar{\theta}\right), (30)

with θ¯˙=dθ¯dτ\dot{\bar{\theta}}=\frac{d\bar{\theta}}{d\tau}. θ¯˙\dot{\bar{\theta}} is the rate of change of the inclination angle θ¯\bar{\theta} as observed by the accerelated observer 𝒜\mathcal{A}. In 𝒜\mathcal{A}, the inclination angle θ¯\bar{\theta} is time-dependent, satisfying equation (14); differentiating (14) with respect to τ\tau, we obtain:

θ¯˙=dθ¯dτ=sinθ¯,\dot{\bar{\theta}}=\frac{d\bar{\theta}}{d\tau}=\sin\bar{\theta}, (31)

and hence (30) can be simplified into:

dθdθ¯=ω¯ωeαξ(1α).\frac{d\theta}{d\bar{\theta}}=\frac{\bar{\omega}}{\omega}e^{\alpha\xi}\left(1-\alpha\right). (32)

Using the transformation of the inclination (polar) angle (30), one could obtain the transformation for the solid angle dΩ=sinθdθdϕd\Omega=\sin\theta d\theta d\phi as follows:

dΩdΩ¯=(ω¯ω)2e2αξ(1α),\frac{d\Omega}{d\bar{\Omega}}=\left(\frac{\bar{\omega}}{\omega}\right)^{2}e^{2\alpha\xi}\left(1-\alpha\right), (33)

where we use the Doppler effect (28) and the fact that the transformation of the azimuth angle satisfies dϕ=dϕ¯d\phi=d\bar{\phi}.

As explained in the preceeding subsection, observer 𝒜\mathcal{A} is accelerated along the +x+x direction; causing it to receive only photons approaching from that direction, namely, those with the +x+x component in their velocity. Moreover, 𝒜\mathcal{A} is moving away from the black-body source at rest with respect to 𝒪.\mathcal{O}. As a consequence, the frequency of the photons received by the accelerated frame 𝒜\mathcal{A} is redshifted by equation (26). However, in constrast with the inertial case presented in Appendix B where the redshift remains constant, the redshift in the accelerated frame 𝒜\mathcal{A} increases in time, with the wavelength of the photons shift progressively towards the infrared range. To obtain the blueshift case, one could invert the situation by flipping the sign of τ\tau (or tt), effectively moving backward in time.

II.4 Density of State in Accelerated Frame

Another important parameter we must determined in the accelerated frame 𝒜\mathcal{A} is the modes distribution of the photons in the cavity. To determine how this quantity transform from the inertial to accelerated frame, first we need to determine transformation of the phase-space volume. In the inertial frame 𝒪\mathcal{O} and the accelerated frame 𝒜\mathcal{A}, the infinitesimal 3-volume elements are defined as, respectively:

d3𝐱\displaystyle d^{3}\mathbf{x} =dxdydz,\displaystyle=dx\wedge dy\wedge dz, (34)
d3𝐱¯\displaystyle d^{3}\mathbf{\bar{x}} =dξdy¯dz¯.\displaystyle=d\xi\wedge d\bar{y}\wedge d\bar{z}.

To obtain the density of states of our black-body system, we need to calculate how much modes are inside a volume element. Let the finite volume element in 𝒪\mathcal{O} be ΔxΔyΔz\Delta x\Delta y\Delta z. The spatial length Δx\Delta x is defined by 2 simultaneous events 𝐩=(ti,xi)\mathbf{p}=\left(t_{i},x_{i}\right) and 𝐪=(tf,xf)\mathbf{q}=\left(t_{f},x_{f}\right) along the xx-axis where ti=tf=0t_{i}=t_{f}=0, xi=0x_{i}=0 and xf=x_{f}=\ell. According to 𝒪,\mathcal{O}, the length between these 2 event is Δx=xfxi=.\Delta x=x_{f}-x_{i}=\ell. For an observer at 𝒜,\mathcal{A}, whose accelerated in the xx-direction with respect to 𝒪,\mathcal{O}, the 2 events 𝐩=(τi,ξi)\mathbf{p}=\left(\tau_{i},\xi_{i}\right) and 𝐪=(τf,ξf)\mathbf{q}=\left(\tau_{f},\xi_{f}\right) are not simultaneous, but are separated by a time interval Δτ=τfτi\Delta\tau=\tau_{f}-\tau_{i}. Let us labeled the spatial length of these 2 events in 𝒜\mathcal{A} as Δξ=ξfξ,i\Delta\xi=\xi_{f}-\xi{}_{i},where τf,τi,ξ,iξf,\tau_{f},\tau_{i},\xi{}_{i},\xi_{f}, could be obtained from (3). Taking the infinitesimal limit Δξdξ\Delta\xi\rightarrow d\xi, the transformation of the infinitesimal spatial length satisfies (6). Notice that the measurement of length in 𝒪\mathcal{O} is obtained from 2 simultaneous event, hence in 𝒪\mathcal{O}, dt=0;dt=0; this is not the case in 𝒜.\mathcal{A}. Therefore, the (infinitesimal) spatial length of 𝐩,𝐪\mathbf{p},\mathbf{q} in 𝒜\mathcal{A} is:

dξ=eαξcoshς¯dx=xα(x2t2)dx.d\xi=e^{-\alpha\xi}\cosh\bar{\varsigma}\,dx=\frac{x}{\alpha\left(x^{2}-t^{2}\right)}\,dx. (35)

Furthermore, the infinitesimal volume element dxdydzdx\wedge dy\wedge dz is perceived by an observer in 𝒜\mathcal{A} as the volume swept by the plane dy¯dz¯=dydzd\bar{y}\wedge d\bar{z}=dy\wedge dz from ξi\xi{}_{i} to ξf,\xi_{f}, along an infinitesimal time interval dτ.d\tau. This is the physical interpretation of d3𝐱¯=dξdy¯dz¯,d^{3}\mathbf{\bar{x}}=d\xi\wedge d\bar{y}\wedge d\bar{z}, see also discussion on Appendix B. Inserting (35) to (34) gives the volume transformation:

d3𝐱¯=eαξcoshς¯d3𝐱=xα(x2t2)d3𝐱.d^{3}\mathbf{\bar{x}}=e^{-\alpha\xi}\cosh\bar{\varsigma}\,d^{3}\mathbf{x}=\frac{x}{\alpha\left(x^{2}-t^{2}\right)}\,d^{3}\mathbf{x}. (36)

For the next step, we need to determine the infinitesimal volume element in the momentum space, this could be obtained from (21)-(23); in the inertial frame 𝒪\mathcal{O} and the accelerated frame 𝒜\mathcal{A}, they are, respectively:

d3𝐩\displaystyle d^{3}\mathbf{p} =dpxdpydpz,\displaystyle=dp_{x}\wedge dp_{y}\wedge dp_{z}, (37)
d3𝐩¯\displaystyle d^{3}\mathbf{\bar{p}} =dpξdp¯ydp¯z.\displaystyle=dp_{\xi}\wedge d\bar{p}_{y}\wedge d\bar{p}_{z}.

The 4-momentum is subjected to the energy-momentum constraint, see Appendix B:

|𝐩4|=pμpμ=pt2+px2+py2+pz2=m2c2.\left|{}^{4}\mathbf{p}\right|=p_{\mu}p^{\mu}=-p_{t}^{2}+p_{x}^{2}+p_{y}^{2}+p_{z}^{2}=m^{2}c^{2}. (38)

Taking the differential of (38) gives:

dpt=1pt(pxdpx+pydpy+pzdpz).dp_{t}=\frac{1}{p_{t}}\left(p_{x}dp_{x}+p_{y}dp_{y}+p_{z}dp_{z}\right). (39)

Furthermore, using (4), one could show that:

x2t2=1α2e2αξ,x^{2}-t^{2}=\frac{1}{\alpha^{2}}e^{2\alpha\xi}, (40)

so that the transformation of the 4-momentum (24) can be written as follows:

pτ\displaystyle p_{\tau} =1α(x2t2)(xpttpx),\displaystyle=\frac{1}{\alpha\left(x^{2}-t^{2}\right)}\left(xp_{t}-tp_{x}\right), (41)
pξ\displaystyle p_{\xi} =1α(x2t2)(xpxtpt),\displaystyle=\frac{1}{\alpha\left(x^{2}-t^{2}\right)}\left(xp_{x}-tp_{t}\right),

with p¯y=py\bar{p}_{y}=p_{y} and p¯z=pz.\bar{p}_{z}=p_{z}. Differentiating (41) with any parameter will give the infinitesimal version of 4-momentum transformation as follows:

dpτ\displaystyle dp_{\tau} =1α(x2t2)[(2t[xpttpxx2t2]px)dt+(pt2x[xpttpxx2t2])dx+xdpttdpx],\displaystyle=\frac{1}{\alpha\left(x^{2}-t^{2}\right)}\left[\left(2t\left[\frac{xp_{t}-tp_{x}}{x^{2}-t^{2}}\right]-p_{x}\right)dt+\left(p_{t}-2x\left[\frac{xp_{t}-tp_{x}}{x^{2}-t^{2}}\right]\right)dx+xdp_{t}-tdp_{x}\right], (42)
dpξ\displaystyle dp_{\xi} =1α(x2t2)[(2t[xpxtptx2t2]pt)dt+(px2x[xpxtptx2t2])dx+xdpxtdpt].\displaystyle=\frac{1}{\alpha\left(x^{2}-t^{2}\right)}\left[\left(2t\left[\frac{xp_{x}-tp_{t}}{x^{2}-t^{2}}\right]-p_{t}\right)dt+\left(p_{x}-2x\left[\frac{xp_{x}-tp_{t}}{x^{2}-t^{2}}\right]\right)dx+xdp_{x}-tdp_{t}\right].

Note that in equation (42), there exist infinitesimal coordinates components, namely dtdt and dxdx. This is due to the fact that the 4-momentum transformation from inertial to accelerated frame (41) varies with coordinates (x,t)\left(x,t\right), in constrast to the transformation between inertial frames (98) in Appendix B, which is independent from the temporal and spatial coordinates.

Similar to the inertial case in Appendix B, using the hypersurface constraint dt=0dt=0 and inserting the infinitesimal momentum constraint (39), (42) becomes:

dpτ\displaystyle dp_{\tau} =1α(x2t2)[(pt2x[xpttpxx2t2])dx+(xpxptt)dpx+x1pt(pydpy+pzdpz)],\displaystyle=\frac{1}{\alpha\left(x^{2}-t^{2}\right)}\left[\left(p_{t}-2x\left[\frac{xp_{t}-tp_{x}}{x^{2}-t^{2}}\right]\right)dx+\left(x\frac{p_{x}}{p_{t}}-t\right)dp_{x}+x\frac{1}{p_{t}}\left(p_{y}dp_{y}+p_{z}dp_{z}\right)\right], (43)
dpξ\displaystyle dp_{\xi} =1α(x2t2)[(px2x[xpxtptx2t2])dx+(xtpxpt)dpxt1pt(pydpy+pzdpz)],\displaystyle=\frac{1}{\alpha\left(x^{2}-t^{2}\right)}\left[\left(p_{x}-2x\left[\frac{xp_{x}-tp_{t}}{x^{2}-t^{2}}\right]\right)dx+\left(x-t\frac{p_{x}}{p_{t}}\right)dp_{x}-t\frac{1}{p_{t}}\left(p_{y}dp_{y}+p_{z}dp_{z}\right)\right],

dp¯y=dpy,d\bar{p}_{y}=dp_{y}, and dp¯z=dpzd\bar{p}_{z}=dp_{z}. Now, one can construct altogether the infinitesimal phase-space volume element in frame 𝒜\mathcal{A} as follows:

d3𝐱¯d3𝐩¯=dξdy¯dz¯dpξdp¯ydp¯z.d^{3}\bar{\mathbf{x}}\wedge d^{3}\bar{\mathbf{p}}=d\xi\wedge d\bar{y}\wedge d\bar{z}\wedge dp_{\xi}\wedge d\bar{p}_{y}\wedge d\bar{p}_{z}.

Inserting (36) and (43) into the equation above, we obtain the transformation of the phase-space volume element from the inertial frame 𝒪\mathcal{O} to the accelerated frame 𝒜\mathcal{A}:

d3𝐱¯d3𝐩¯=x(xtpxpt)α2(x2t2)2d3𝐱d3𝐩.d^{3}\bar{\mathbf{x}}\wedge d^{3}\bar{\mathbf{p}}=\frac{x\left(x-t\frac{p_{x}}{p_{t}}\right)}{\alpha^{2}\left(x^{2}-t^{2}\right)^{2}}\,d^{3}\mathbf{x}\wedge d^{3}\mathbf{p}. (44)

Note that by the definition of the wedge product ,\wedge, the terms containing equivalent components of differential forms, i.e. dxdxdx\wedge dx, vanishes. Using (40) and then (21), one could simplify (44) into:

d3𝐱¯d3𝐩¯d3𝐱d3𝐩=(1pxpttanhς¯)e2αξcosh2ς¯.\frac{d^{3}\bar{\mathbf{x}}d^{3}\bar{\mathbf{p}}}{d^{3}\mathbf{x}d^{3}\mathbf{p}}=\left(1-\frac{p_{x}}{p_{t}}\tanh\bar{\varsigma}\right)e^{-2\alpha\xi}\cosh^{2}\bar{\varsigma}. (45)

One the other hand, we have, from (24):

pτpt=(1pxpttanhς¯)eαξcoshς¯,\frac{p_{\tau}}{p_{t}}=\left(1-\frac{p_{x}}{p_{t}}\tanh\bar{\varsigma}\right)e^{-\alpha\xi}\cosh\bar{\varsigma},

and therefore, using (21)-(23):

d3𝐱¯d3𝐩¯d3𝐱d3𝐩=pτpteαξcoshς¯=ω¯ωeαξcoshς¯,\frac{d^{3}\bar{\mathbf{x}}d^{3}\bar{\mathbf{p}}}{d^{3}\mathbf{x}d^{3}\mathbf{p}}=\frac{p_{\tau}}{p_{t}}e^{-\alpha\xi}\cosh\bar{\varsigma}=\frac{\bar{\omega}}{\omega}e^{-\alpha\xi}\cosh\bar{\varsigma}, (46)

where we write the transformation (45) in terms of the Doppler factor ω¯ω\frac{\bar{\omega}}{\omega}. (46) is the relation between the phase-space volume element in 𝒪\mathcal{O} and 𝒜.\mathcal{A}.

Let us define the relativistic distribution function, or density of state f(𝐱,𝐩)f\left(\mathbf{x},\mathbf{p}\right) as the number of world-lines that cross the phase-space, i.e, the states dNdN, per phase-space volume element (Liboff, ):

f(𝐱,𝐩)=dNd3𝐱d3𝐩.f\left(\mathbf{x},\mathbf{p}\right)=\frac{dN}{d^{3}\mathbf{x}d^{3}\mathbf{p}}. (47)

We assume that the density of state f(𝐱,𝐩)f\left(\mathbf{x},\mathbf{p}\right) is invariant under coordinate transformation, namely:

dNd3𝐱d3𝐩=dN¯d3𝐱¯d3𝐩¯,\frac{dN}{d^{3}\mathbf{x}d^{3}\mathbf{p}}=\frac{d\bar{N}}{d^{3}\mathbf{\bar{x}}d^{3}\mathbf{\bar{p}}},

with dN¯d\bar{N} is the number of states per phase-space volume d3𝐱¯d3𝐩¯d^{3}\mathbf{\bar{x}}\,d^{3}\mathbf{\bar{p}} of the accelerated frame 𝒜\mathcal{A}. This is a reasonable assumption, and for a detailed explanation, one could refer to Appendix B. Some previous works in the existing literature have assume the invariance of the number of states dNdN instead of f(𝐱,𝐩)f\left(\mathbf{x},\mathbf{p}\right), however, this assumption is not suitable for our case because our treatment in this works relies only on a coordinate transformations at some part of the Minkowski space, i.e., the Rindler wedge. While these transformations leave the world-lines invariant, they alter the unit volume due to the coordinate transformation, leading to a different count of world-lines crossing the unit volume. Similar arguments are implicitly employed in (Peebles, ; Heer, ; Henry, ) as well.

By the invariance of the density of state under (infinitesimal) coordinate transformation, one could obtain the relation between the number of states/modes in inertial frame 𝒪\mathcal{O} and the accelerated frame 𝒜\mathcal{A} as follows:

dN¯dN=d3𝐱¯d3𝐩¯d3𝐱d3𝐩=ω¯ωeαξcoshς¯.\frac{d\bar{N}}{dN}=\frac{d^{3}\bar{\mathbf{x}}d^{3}\bar{\mathbf{p}}}{d^{3}\mathbf{x}d^{3}\mathbf{p}}=\frac{\bar{\omega}}{\omega}e^{-\alpha\xi}\cosh\bar{\varsigma}. (48)

II.5 Black-body Radiation in Accelerated Frame

We already had all the Rindler-transformed variables to derive the Planck’s law in a uniformly-accelerated frame. Similar to the inertial frame case that we derived in Appendix B, in the rest frame 𝒪\mathcal{O}, we use the Planck’s law version of equation (77) as follows:

dN=ω24π3c31ewω/ωmax1dωdΩdV,dN=\frac{\omega^{2}}{4\pi^{3}c^{3}}\frac{1}{e^{\nicefrac{{w\omega}}{{\omega_{\mathrm{max}}}}}-1}d\omega d\Omega dV, (49)

where the term containing the temperature of the black-body, i.e., kBTk_{B}T, is replaced by ωmax\omega_{\mathrm{max}}, the angular frequency that maximize the radiation energy EE of the black-body. This is possible by the Wien’s displacement law, see a detailed explanation on the Appendix A. The reason for this substitution is to avoid the problem of the temperature in moving bodies. As explained in the Introduction, currently, there is no consensus regarding the transformation of temperature in moving bodies (Farias, ). However, we convincingly agree on the transformation of the black-body’s frequency, hence, if ω\omega satisfies the relativistic Doppler effect, then so does the ωmax\omega_{\mathrm{max}}. The dimensionless constant ww in (49) is the Wien coefficient, obtained from solving the non-linear equation that arise from the maximization of energy with respect to ω\omega. One could refer to Appendix B for a detailed derivation.

In this section, our objective is to ascertain the validity of the distribution (49) in the accelerated frame 𝒜\mathcal{A}. Inserting the relativistic Doppler effect (29), the transformation of solid angle (33), the volume contraction (36), and the transformation of the number of modes distribution (48) to (49), we obtain:

dN¯=(1α)e2αξ4π3c3ω¯2ewω¯/ωmax¯1dω¯dΩ¯dV¯.d\bar{N}=\left(1-\alpha\right)\frac{e^{2\alpha\xi}}{4\pi^{3}c^{3}}\frac{\bar{\omega}^{2}}{e^{\nicefrac{{w\bar{\omega}}}{{\overline{\omega_{\mathrm{max}}}}}}-1}d\bar{\omega}d\bar{\Omega}d\bar{V}. (50)

One can observe that in the right hand side of (50) there is a scale factor:

ε[α,ξ]=e2αξ(1α),\varepsilon\left[\alpha,\xi\right]=e^{2\alpha\xi}\left(1-\alpha\right), (51)

that prevents the distribution (50) to be the original Planckian (49). The scale factor depends on the acceleration α\alpha and the spatial coordinate ξ\xi. For the case with α=0\alpha=0 where the acceleration vanishes, the distribution becomes Planckian as in (49). The fundamental distinction between the distribution (50) and the distribution (52) in the moving inertial frame lies in the fact that all the variables in the (50) are coordinate-dependent, namely, they vary with τ\tau and ξ\xi, which represent the proper time and position in 𝒜\mathcal{A}. This is due to the fact that the transformations of these variables, namely (29), (33), (36), and (48) are coordinate-dependent. However, since we are considering all the possible value of frequencies from 0ω<0\leq\omega<\infty, the time dependence of ω¯\bar{\omega} does not explicitly shown in the law. Consequently, the ’black-body’ spectrum (50) in the accelerated frame 𝒜\mathcal{A} experiences a redshift that increase over time and varying with position, where the wavelength of the photon increases toward the infrared direction. This behavior distinguishes (50) from the black-body spectrum in the inertial frames, where the spectrum shifts similarly occur, but the magnitudes of these shifts remain constant, as indicated by equation (100).

Similar with the inertial case in Appendix B, ωmax¯\overline{\omega_{\mathrm{max}}} is the “maximal” frequency of the distribution (49) of a black-body in the rest frame, viewed by the accelerated observer 𝒜\mathcal{A}. Let us check if ωmax¯\overline{\omega_{\mathrm{max}}} also maximizes the energy in distribution (50). By multiplying the left and right hand side of (50) with ω¯,\hbar\bar{\omega}, one obtains:

dE¯=ε4π3c3ω¯3ewω¯/ωmax¯1ρ¯[ω¯]Edω¯dΩ¯dV¯.d\bar{E}=\underset{\bar{\rho}{}_{E}\left[\bar{\omega}\right]}{\underbrace{\frac{\varepsilon}{4\pi^{3}c^{3}}\frac{\hbar\bar{\omega}^{3}}{e^{\nicefrac{{w\bar{\omega}}}{{\overline{\omega_{\mathrm{max}}}}}}-1}}}d\bar{\omega}d\bar{\Omega}d\bar{V}. (52)

E¯\bar{E} is the energy of modes with frequency ω¯\bar{\omega} observed in the accelerated frame 𝒜\mathcal{A}. The value of ω¯\bar{\omega} that maximize E¯\bar{E} can be obtained by solving the following equation:

dρ¯Edω¯=ε4π3c3ω¯2exp(wω¯ωmax¯)1(3wω¯ωmax¯exp(wω¯ωmax¯)exp(wω¯ωmax¯)1)=0.\frac{d\bar{\rho}{}_{E}}{d\bar{\omega}}=\frac{\varepsilon}{4\pi^{3}c^{3}}\frac{\hbar\bar{\omega}^{2}}{\exp\left(w\frac{\bar{\omega}}{\overline{\omega_{\mathrm{max}}}}\right)-1}\left(3-w\frac{\bar{\omega}}{\overline{\omega_{\mathrm{max}}}}\frac{\exp\left(w\frac{\bar{\omega}}{\overline{\omega_{\mathrm{max}}}}\right)}{\exp\left(w\frac{\bar{\omega}}{\overline{\omega_{\mathrm{max}}}}\right)-1}\right)=0. (53)

Note that the scale factor ε\varepsilon does not affect the equation we need to solve. Let us define w¯=wω¯ωmax¯\bar{w}=w\frac{\bar{\omega}}{\overline{\omega_{\mathrm{max}}}}, the equation from (53) that we need to solve becomes:

w¯=3(1ew¯),\bar{w}=3\left(1-e^{-\bar{w}}\right), (54)

which is the same as in the inertial cases, namely (74) and (115). The solution is then similar:

w¯=3+W0(3e3),\bar{w}=3+W_{0}\left(-3e^{-3}\right),

with W0W_{0} is the Lambert WW-function (Euler, ). Hence, the Wien’s coefficient is invariant, even under the transformation from the inertial to accelerated frame: w¯=w\bar{w}=w. Suppose ω¯=ω¯max\bar{\omega}=\bar{\omega}_{\mathrm{max}} is the solution to (53), therefore:

w¯w=ω¯maxωmax¯=1,\frac{\bar{w}}{w}=\frac{\bar{\omega}_{\mathrm{max}}}{\overline{\omega_{\mathrm{max}}}}=1,

or simply ω¯max=ωmax¯\bar{\omega}_{\mathrm{max}}=\overline{\omega_{\mathrm{max}}}: the transformation (26) sends “maximal” frequency in an inertial frame to a “maximal” frequency in an accelerated frame at an instantaneous time. (50) becomes:

dN¯=ε4π3c3ω¯2ewω¯/ω¯max1dω¯dΩ¯dV¯.d\bar{N}=\frac{\varepsilon}{4\pi^{3}c^{3}}\frac{\bar{\omega}^{2}}{e^{\nicefrac{{w\bar{\omega}}}{{\bar{\omega}_{\mathrm{max}}}}}-1}d\bar{\omega}d\bar{\Omega}d\bar{V}. (55)

As previously mentioned, the presence of the scale factor ε\varepsilon introduces deviations from the Planckian distribution of the black-body radiation in the accelerated frame 𝒜\mathcal{A}. For an observer 𝒜\mathcal{A} moving along ξ=0\xi=0, it will observe a Planckian distribution but with a redshift increasing over time. For another accelerated observer moving along a constant ξ>0\xi>0 , it will experience a similar increasing redshift, but now the distribution will be scaled by a factor of e2αξ(1α)e^{2\alpha\xi}\left(1-\alpha\right). Such an observer will interpret this as a change in the emissivity of the black-body source; this facts will be analyzed in detail in the next section.

It needs to be kept in mind that the Planck formula (55) applies in the right Rindler wedge (RRW). With a similar derivation, one could show that the same formula is also valid for the left Rindler wedge, but with the Doppler formula (26)-(27) changed into:

ω¯ω\displaystyle\frac{\bar{\omega}}{\omega} =eαξ(coshς¯+cosθsinhς¯),\displaystyle=e^{-\alpha\xi}\left(\cosh\bar{\varsigma}+\cos\theta\sinh\bar{\varsigma}\right), (56)
=eαξ(cosθcoshς¯+sinhς¯)cosθ¯.\displaystyle=e^{-\alpha\xi}\frac{\left(\cos\theta\cosh\bar{\varsigma}+\sinh\bar{\varsigma}\right)}{\cos\bar{\theta}}. (57)

This is due to the fact that there exists another solution to the dispersion relation (20), namely ω=k.\omega=-k.

III Physics in Rindler Coordinate

III.1 Acceleration and Deceleration

The "Planck’s law" (55) describes the frequency spectrum measured by an observer 𝒜\mathcal{A} whose accelerating away from the radiation source along the +x+x direction, as described in Section II. This observer moves with an initial velocity β=tanhς\beta=\tanh\varsigma and undergoes an acceleration with magnitude aa. Let us add three possible scenarios besides the ones we had considered in the previous section. First, note that the origin of the rest frame 𝒪\mathcal{O} does not coincide to the origin of frame 𝒜,\mathcal{\mathcal{A}}, namely, the point (t,x)=(0,0)\left(t,x\right)=\left(0,0\right) at 𝒪\mathcal{O} is equivalent to point of (τ,ξ)=(ςα,)\left(\tau,\xi\right)=\left(-\frac{\varsigma}{\alpha},-\infty\right) at 𝒜,\mathcal{\mathcal{A}}, while the point of (τ,ξ)=(0,0)\left(\tau,\xi\right)=\left(0,0\right) at 𝒜,\mathcal{\mathcal{A}}, is the point of (t,x)=1α(sinhς,coshς)\left(t,x\right)=\frac{1}{\alpha}\left(\sinh\varsigma,\cosh\varsigma\right) at 𝒪,\mathcal{O}, see FIG. 2 and FIG. 3.

Refer to caption
Figure 2: The accelerated away case as observed by the rest frame 𝒪\mathcal{O}. Let us set α=1\alpha=1 and ς=0\varsigma=0, so that the origin of 𝒜\mathcal{A} is at (t,x)=(0,1)\left(t,x\right)=\left(0,1\right). The source is located at x=1x=1 (or ξ=0\xi=0), at rest with respect to 𝒪\mathcal{O}, and its worldline with respect to 𝒪\mathcal{O} is described by the orange-dashed line. 𝒜\mathcal{A} is accelerated with respect to 𝒪\mathcal{O}, hence its wordline is described by the blue hyperboloid. Note that the distance between the source and 𝒜\mathcal{A} increases, and so does the time interval to receive the signal. The source signal at the horizon (and above) will never reach 𝒜\mathcal{A} since it needs infinite time to reach the observer.
Refer to caption
Figure 3: The accelerated away case as observed by the rest frame 𝒜\mathcal{A}. The origin of 𝒪\mathcal{O} is at (τ,ξ)=(0,)\left(\tau,\xi\right)=\left(0,-\infty\right) (not shown). Now the worldline of 𝒜\mathcal{A} is constant in the Rindler coordinate, described by the blue line ξ=0\xi=0. The source, located at ξ=0\xi=0 when τ=0\tau=0, is moving away from 𝒜,\mathcal{A}, described by the orange-dashed hyperboloid. This is exactly the same situation as in FIG. 2 but described in the Rindler coordinate.

To simplify the analysis, let us set the black-body source to be located at any point in the range of 0x1αcoshς0\leq x\leq\frac{1}{\alpha}\cosh\varsigma or 1αcoshς<x\frac{1}{\alpha}\cosh\varsigma<x\leq\infty, depending on the scenarios. With respect to 𝒪\mathcal{O}, let frame 𝒜\mathcal{A} have an initial point at (t,x)=1α(sinhς,coshς)\left(t,x\right)=\frac{1}{\alpha}\left(\sinh\varsigma,\cosh\varsigma\right), that means the clock in 𝒜\mathcal{A} is slightly delayed by the amount of 1αsinhς\frac{1}{\alpha}\sinh\varsigma. Then 𝒜\mathcal{A} moves with the time-dependent rapidity ς¯[τ]\bar{\varsigma}\left[\tau\right], which could satisfy one of the four following case:

ς¯[τ]=±ς±ατ,\bar{\varsigma}\left[\tau\right]=\pm\varsigma\pm\alpha\tau, (58)

noting that the initial rapidity ς0\varsigma\geq 0. The sign ±\pm in front of the initial rapidity ς\varsigma depends on the rate of change of the position of observer 𝒜\mathcal{A} with respect to the source; it is positive if their distance increases in time, and negative if it decreases. Meanwhile, the sign ±\pm in front of ατ\alpha\tau depends on the rate of change of ς\varsigma. If it increases in time (accelerating), then it has the same sign with ς\varsigma, and vice versa. Let us list all the possible cases. The first case, where ς¯[τ]=ς+ατ\bar{\varsigma}\left[\tau\right]=\varsigma+\alpha\tau, corresponds to the case described in Section II, used to derive the Planck’s law (55), see FIG. 4.

Refer to caption
Figure 4: Case 1: The “accelerating away” case could be realized in the right Rindler wedge by locating the source (orange dashed-line) at 0x1αcoshς0\leq x\leq\frac{1}{\alpha}\cosh\varsigma, where 1αcoshς\frac{1}{\alpha}\cosh\varsigma is the origin of 𝒜.\mathcal{A}. In this case, the source is located at x=0.6x=0.6 and the origin of 𝒜\mathcal{A} is at x=1,x=1, given α=1\alpha=1 and ς=0\varsigma=0. The trajectory of observer 𝒜\mathcal{A} is described by the blue curve.

The second case, where ς¯[τ]=ςατ\bar{\varsigma}\left[\tau\right]=\varsigma-\alpha\tau, describes an observer decelerating away from the source; it moves in the +x+x direction, but the speed is decreasing in time; we will only consider this case until the observer reach zero velocity as a result of deceleration. See FIG. 5.

Refer to caption
Figure 5: Case 2: The second case could be similarly realized in the right Rindler wedge by locating the source at 1αcoshς<x\frac{1}{\alpha}\cosh\varsigma<x\leq\infty, in this case x=2.6x=2.6. The observer’s trajectory coincides with the source at some time in the past τi<0\tau_{i}<0. At τi,\tau_{i}, the (initial) velocity of 𝒜\mathcal{A} is ς\varsigma; which decreases as time progresses by the amount of ατ-\alpha\tau, and eventually, reaching zero at τf=0\tau_{f}=0. At this time, the observer is at a finite distance from the source. The collection of events in the time interval of τiττf\tau_{i}\leq\tau\leq\tau_{f} is considered as the “decelerating away” case.

The third case with ς¯[τ]=ς+ατ\bar{\varsigma}\left[\tau\right]=-\varsigma+\alpha\tau describes an observer that is decelerating towards the source, with a decreasing speed. Similar to the second case, we will only consider the case until the velocity is zero, see FIG. 6.

Refer to caption
Figure 6: Case 3: To describe the “decelerated towards” case in the right Rindler wedge, the source is located at 0x1αcoshς0\leq x\leq\frac{1}{\alpha}\cosh\varsigma, in this case x=0.6x=0.6. The observer approaches the source from x=x=\infty with an initial velocity ς-\varsigma, starts to decelerate, and then stops at τ=0\tau=0.

The last case with ς¯[τ]=ςατ\bar{\varsigma}\left[\tau\right]=-\varsigma-\alpha\tau describes an observer accelerating towards the source, see FIG. 7.

Refer to caption
Figure 7: Case 4: To describe the “accelerated towards” in the right Rindler wedge, the source is located at 1αcoshς<x\frac{1}{\alpha}\cosh\varsigma<x\leq\infty, namely x=2.6x=2.6. The observer approaches the source from the origin of 𝒜\mathcal{A} with an increasing velocity and then coincides with the source at a time τ>0\tau>0 in the future.

The concept of acceleration and deceleration is formally defined as the rate of change of the speed, specifically, an acceleration if d|ς¯|dτ>0\frac{d\left|\bar{\varsigma}\right|}{d\tau}>0 and deceleration if d|ς¯|dτ<0\frac{d\left|\bar{\varsigma}\right|}{d\tau}<0. It is also important to highlight that the notion of acceleration/deceleration do not necessarily correlate to redshift/blueshift, as elaborated in Subsection III C.

III.2 The Relativistic Beaming

The alteration of the inclination angle in the moving frame leads to a peculiar phenomenon known as relativistic beaming where, in moving a frame, the rays of light shows a tendency to either converge or diverge in the direction of motion, depending on whether the observer is approaching or receding away from the source (Cohen, ). For a frame 𝒪\mathcal{O}^{\prime} moving away with constant velocity of β\beta , the inclination angle θ\theta will transform by (90). For the case where the observer is moving towards the source, one could flip the sign in front β\beta to become positive. The photons that reach the observer will have the inclination angle within the range of π2θπ2-\frac{\pi}{2}\leq\theta\leq\frac{\pi}{2}. It can be demonstrated that for the “moving towards” case (+β+\beta), the condition 1γ(1+βcosθ)1\leq\gamma\left(1+\beta\cos\theta\right)\leq\infty holds, implying |θ||θ|\left|\theta^{\prime}\right|\leq\left|\theta\right|. It implies that for an observer approaching the source, the light rays tend to converge to the direction of motion, see FIG. 8.

Refer to caption
Figure 8: The alteration of the inclination angle on a frame approaching the source with velocity ς\varsigma. The slopes of the collection of curves θ[θ]\theta^{\prime}\left[\theta\right] never exceed 1, except for the case where ς=0\varsigma=0, described by the dashed-line. As the rapidity ς\varsigma increases, the slope of the curve decreases, hence |θ||θ|.\left|\theta^{\prime}\right|\leq\left|\theta\right|.

For the “moving away” case (β-\beta), there are 2 conditions that affect the light rays bending. First, if 1γ(1βcosθ)1\leq\gamma\left(1-\beta\cos\theta\right)\leq\infty holds, then |θ||θ|\left|\theta^{\prime}\right|\leq\left|\theta\right|, implying that in this interval, the light rays converge. Second, if the condition 0γ(1βcosθ)10\leq\gamma\left(1-\beta\cos\theta\right)\leq 1 holds, the light rays diverge in the direction of motion, namely |θ||θ|\left|\theta^{\prime}\right|\geq\left|\theta\right|. This case occurs when cosθββ2\cos\theta\leq\frac{\beta}{\beta-2}, see FIG. 9.

Refer to caption
Figure 9: The alteration of the inclination angle in a frame moving away from the source with velocity ς\varsigma. The dashed-line is the case with the velocity ς=0\varsigma=0, hence θ=θ\theta^{\prime}=\theta. As the velocity increases, so does the difference between θ\theta^{\prime} and θ\theta. For a specific θ\theta and ς\varsigma, a threshold exists at θ=±π2\theta^{\prime}=\pm\frac{\pi}{2}. For |θ||±π2|\left|\theta^{\prime}\right|\leq\left|\pm\frac{\pi}{2}\right|, the angle diverge, represented by the slope 1\geq 1. However, there exist combination values of θ\theta and ς\varsigma that gives |θ||±π2|\left|\theta^{\prime}\right|\geq\left|\pm\frac{\pi}{2}\right|, which implies that the angle spread “backward” in the opposite direction of the observer’s velocity ς\varsigma. Thus the observer will perceive the angle as converging.

It is worth noting that in a moving inertial frame, θ\theta^{\prime} is time-independent, resulting in a constant angle deviation.

For the case of an accelerated frame 𝒜\mathcal{A}, the inclination angle θ\theta will be perceived as time-dependent, satisfying (15). Therefore, a similar beaming phenomenon as in the moving inertial cases occurs. In the scenario where the accelerated observer approaching the source, θ¯\bar{\theta} decreases in time, resulting in a convergence at τ\tau\rightarrow\infty, see FIG. 10.

Refer to caption
Figure 10: The evolution of the inclination angle perceived by an approaching accelerated observer 𝒜\mathcal{A}, namely θ¯\bar{\theta}, for different value of θ\theta. θ¯0\bar{\theta}\rightarrow 0 as τ\tau\rightarrow\infty.

In the scenario where the accelerated observer receding from the source, θ\theta will diverge until it reaches the maximal values, namely π/2\pi/2 for 0θπ/20\leq\theta\leq\pi/2 and π/2-\pi/2 for π/2θ0-\pi/2\leq\theta\leq 0. From this point, it will stop diverging and start to converge as time increases, as illustrated in FIG. 11.

Refer to caption
Figure 11: The evolution of the observation angle θ¯\bar{\theta} with respect to proper time τ\tau, as observed by an accelerated frame 𝒜\mathcal{A}, moving away from the source. The light rays initially dispersed at some time interval but subsequently stopped diverging and started to converge as time increased. Different initial angle θ\theta will results in different peaking-time, but the light rays will eventually converge at τ\tau\rightarrow\infty.

It is important to emphasize that our analysis of the relativistic beaming phenomenon is restricted to specific cases, namely Cases 1 (accelerated away) and 4 (accelerated toward). The decelerated cases (namely, Case 2 and 3) can be derived by inserting the condition (58) into (15). In the context of relativistic beaming, Case 3 (decelerating toward) will experience a similar phenomenon to Case 1 (accelerating away), while Case 2 (decelerating away) is similar to Case 4 (accelerating towards). Observing only the beaming of lightrays, these two pairs of cases are indistinguishable.

To be more precise, one could obtain the rate of change of the inclination angle by inserting (31) to (15)- as follows:

dθ¯dτ=sinθcoshς¯cosθsinhς¯.\frac{d\bar{\theta}}{d\tau}=\frac{\sin\theta}{\cosh\bar{\varsigma}-\cos\theta\sinh\bar{\varsigma}}. (59)

Note that the inclination angle satisfies π2θπ2-\frac{\pi}{2}\leq\theta\leq\frac{\pi}{2}; in these ranges, coshς¯\cosh\bar{\varsigma}, cosθsinhς¯\cos\theta\sinh\bar{\varsigma} are positive definite, and the quantity (coshς¯cosθsinhς¯)1\left(\cosh\bar{\varsigma}-\cos\theta\sinh\bar{\varsigma}\right)^{-1} could be positive or negative, depending on whether coshς¯>cosθsinhς¯,\cosh\bar{\varsigma}>\cos\theta\sinh\bar{\varsigma}, or vice versa. For the case where the inclination angle at the rest frame satisfies 0θπ20\leq\theta\leq\frac{\pi}{2}, sinθ\sin\theta will be positive definite, so if the quantity (coshς¯cosθsinhς¯)1\left(\cosh\bar{\varsigma}-\cos\theta\sinh\bar{\varsigma}\right)^{-1} is positive, then dθ¯dτ0\frac{d\bar{\theta}}{d\tau}\geq 0; the angle at 𝒜\mathcal{A} will diverge. Otherwise, dθ¯dτ0\frac{d\bar{\theta}}{d\tau}\leq 0, and the angle will converge. Meanwhile, for the case of π2θ0-\frac{\pi}{2}\leq\theta\leq 0, sinθ\sin\theta will be negative definite, so the angle at 𝒜\mathcal{A} will converge if (coshς¯cosθsinhς¯)1\left(\cosh\bar{\varsigma}-\cos\theta\sinh\bar{\varsigma}\right)^{-1} is positive -resulting in dθ¯dτ0\frac{d\bar{\theta}}{d\tau}\leq 0-, otherwise it will diverge.

III.3 The Redshift and Blueshift

Considering the quantity 𝐩,𝐱=pμxμ\left\langle\mathbf{p},\mathbf{x}\right\rangle=p_{\mu}x^{\mu} with xμx^{\mu} is a 4-vector and pμp_{\mu} is a 4-momentum for a single photon, one could show that 𝐩,𝐱\left\langle\mathbf{p},\mathbf{x}\right\rangle is invariant under Lorentz transformation. A simple way to prove this is to calculate pμxμp_{\mu}x^{\mu} at frame 𝒪\mathcal{O} using (96), and pμxμp^{\prime}_{\mu}x^{\prime\mu} at frame 𝒪\mathcal{O}^{\prime} using (97). With the help of (78), (89)-(90), and (101), one could show that pμxμ=pμxμp_{\mu}x^{\mu}=p^{\prime}_{\mu}x^{\prime\mu}.

Furthermore, with the invariance of 𝐩,𝐱\left\langle\mathbf{p},\mathbf{x}\right\rangle in inertial frames, one can conclude that a plane wave moving in n^\hat{n}-direction as follows:

ψ[x,t]=ei𝐩,𝐱=eiω(t+x),\psi\left[x,t\right]=e^{\frac{i}{\hbar}\left\langle\mathbf{p},\mathbf{x}\right\rangle}=e^{i\omega\left(-t+x\right)}, (60)

is invariant under Lorentz transformation. Hence, a plane wave in 𝒪\mathcal{O} is also a plane wave in 𝒪\mathcal{O}^{\prime} but with the change in frequency satisfying the Doppler effect (101). The Doppler effect causes the frequency in the inertial moving frame 𝒪\mathcal{O}^{\prime} to undergo a redshift if it is moving away from, -and blueshift if it is approaching, the source. In the inertial frames, the redshift and blueshift are constant in time.

However, this is not the case in an accelerated frame. Let us take only the temporal part of (60) at 𝒪\mathcal{O} for simplicity, namely eiΘ,e^{-i\Theta}, with Θ=ωt\Theta=\omega t is the phase-angle. The accelerated observer 𝒜\mathcal{A} will measure the photon of constant frequency ω\omega at 𝒪\mathcal{O} as having frequency ω¯[τ]\bar{\omega}\left[\tau\right] satisfying the relativistic Doppler effect (29). Without loosing of generality, let us set θ=0\theta=0 so that ω¯[τ]\bar{\omega}\left[\tau\right] satisfies:

ω¯[τ]=ωeα(ξ±τ)ς,\bar{\omega}\left[\tau\right]=\omega e^{-\alpha\left(\xi\pm\tau\right)-\varsigma}, (61)

considering only the case of ς¯=ς±ατ\bar{\varsigma}=\varsigma\pm\alpha\tau from (58). The phase-angle of the wave of the single photon with respect to 𝒜\mathcal{A} is:

Θ¯[τ]=ω¯𝑑τ=ωαeα(ξ±τ)ς=ω¯[τ]α.\bar{\Theta}\left[\tau\right]=\int\bar{\omega}\,d\tau=\mp\frac{\omega}{\alpha}e^{-\alpha\left(\xi\pm\tau\right)-\varsigma}=\mp\frac{\bar{\omega}\left[\tau\right]}{\alpha}. (62)

Θ¯\bar{\Theta} is time-dependent, since ω¯\bar{\omega} is a function of time τ\tau. This is different with the phase-angle in the inertial case, where Θ\Theta could be written simply as ωt\omega t, since ω\omega in the inertial case is independent of time. Therefore, a single mode of (the temporal part of) the radiation wave with respect to 𝒜\mathcal{A} is:

ψ¯[τ,ξ]=eiΘ¯=e±iωαeα(ξ±τ)ς.\bar{\psi}\left[\tau,\xi\right]=e^{-i\bar{\Theta}}=e^{\pm i\frac{\omega}{\alpha}e^{-\alpha\left(\xi\pm\tau\right)-\varsigma}}. (63)

(63) is not necessarily a plane-wave, in contrast to its inertial counterpart (60), since the frequency measured by 𝒜\mathcal{A} experiences a shift increasing in time. The direction of the frequency’s shift depends on the rate of change of the frequency with respect to the proper time, namely:

dω¯dτ=αωeα(ξ±τ)ς;\frac{d\bar{\omega}}{d\tau}=\mp\alpha\omega e^{-\alpha\left(\xi\pm\tau\right)-\varsigma}; (64)

dω¯dτ<0\frac{d\bar{\omega}}{d\tau}<0 gives a redshift while dω¯dτ>0\frac{d\bar{\omega}}{d\tau}>0 gives a blueshift.

Now, having the definition of redshift and blueshift, let us consider the 4 cases of accelerated observers in Section III A. Note that by the relativistic Doppler effect, all these 4 cases will experienced time-dependent shifts on their frequencies: the first and the third cases will have their frequency shifting in time towards the infrared direction (redshift), while the second and the last cases will have the shift toward the ultraviolet direction (blueshift). These pairs of cases could be distinguished by their instantaneous frequency with respect to the frequency of the source at rest: the first case has ω¯[τ]ω\bar{\omega}\left[\tau\right]\leq\omega, while the third has ω¯[τ]ω\bar{\omega}\left[\tau\right]\geq\omega, although they experience increasing time-dependent redshift on their frequencies. A similar condition occurs for the second case, with ω¯[τ]ω\bar{\omega}\left[\tau\right]\leq\omega, and the last case, with ω¯[τ]ω\bar{\omega}\left[\tau\right]\geq\omega; where both experience increasing blueshift on their frequencies. Comparing the definition on frequency’s shift with the definition of acceleration and deceleration in Section IIIA, it is clear that that they are not correlated, i.e., acceleration does not always correspond to redshift, while deceleration does not always correspond to blueshift.

The "maximal" frequency ωmax\omega_{\mathrm{max}} similarly transformed as (61), therefore the rate of change of ω¯max\bar{\omega}_{\mathrm{max}} satisfies (64). The peak of the spectrum will shift towards higher or lower frequencies as time progresses, depending on which one from the 4 cases is satisfied. See FIG. 12.

Refer to caption
Figure 12: The 3D diagram of the ”Planck” spectrum evolving in time, as perceived by an accelerated observer 𝒜\mathcal{A}, approaching the source (Case 4). The (x,y,z)\left(x,y,z\right)-axes are, respectively, the time τ\tau, frequency ω¯\bar{\omega}, and energy density ρ¯E\bar{\rho}{}_{E}. Notice that the dynamics in (55) is only observed in the parameter of "maximal" frequency ω¯max.\bar{\omega}_{\mathrm{max}}.

One may wonder the relationship between the frequency shift discussed in this work and the gravitational red/blueshift phenomena (Eddington, ). Gravitational frequency shift is an alteration in the photon’s frequency as it travels a gravitational field. The physical clock at different location along the field has different rate of time, therefore, photon appears to be red/blue-shifted relative to the frequency of the clock (Okun, ). On the other hand, the frequency shift considered in this work is due to its velocities; it has the same origin with the relativistic Doppler effect, with the main difference lies in the time-dependent nature of this shift. Now, by the equivalence principle, that asserts the indistinguishability of gravity with acceleration, one might expect an equivalency between the frequency shift due to acceleration and gravitational red/blueshift. However, given a stationary source of gravity, the gravitational redshift remains constant over time, while the frequency shift due to acceleration is time-dependent. It is worth noting that the equivalence principle holds only locally in space and in time, therefore if one considers a gravitational field in a small region in space, it is approximately indistinguishable with a uniform acceleration. This is also valid in reverse: if one consider an acceleration in a (nearly) instantaneous time, it is indistinguishable from a gravitational field. In this limit, the frequency shift due to acceleration is practically indistinguishable to the gravitational red/blueshift.

III.4 Annihilation and Creation of Modes

The number of modes/quanta of inertial moving frame, in general, will differ with the number in a rest frame. This is due to the fact that different inertial coordinates will measure different size of phase-space volume element, hence, they will calculate different number of worldlines that crosses the phase-space volume element. However, there will be no creation or annihilation of quanta/modes, if their number NN is conserved in time. In contrast, the number of states dN¯d\bar{N} in the accelerated frame evolves in time. One could obtain the rate of the quanta/modes production as follows. Let us consider only the accelerated and decelerated away case where ς¯=ς±ατ\bar{\varsigma}=\varsigma\pm\alpha\tau. Inserting (26) to (48), and setting θ=0\theta=0 for simplicity gives:

dN¯dN=e2αξ2(1+e2(ς±ατ)):=σ[α,ξ,τ].\frac{d\bar{N}}{dN}=\frac{e^{-2\alpha\xi}}{2}\left(1+e^{-2\left(\varsigma\pm\alpha\tau\right)}\right):=\sigma\left[\alpha,\xi,\tau\right].

From here we can obtain the rate of quanta/modes production in 𝒜\mathcal{A}, with respect to τ\tau:

dσdτ=αe2αξe2(ς±ατ)=α(ω¯ω)2.\frac{d\sigma}{d\tau}=\mp\alpha e^{-2\alpha\xi}e^{-2\left(\varsigma\pm\alpha\tau\right)}=\mp\alpha\left(\frac{\bar{\omega}}{\omega}\right)^{2}. (65)

For the case where 𝒜\mathcal{A} is accelerated away from the black-body source at rest, the number of modes/quanta observed by 𝒜\mathcal{A} is decreasing exponentially in time, as described by the minus sign in (65). This corresponds to the increasing redshift experienced by 𝒜\mathcal{A}. For the case where 𝒜\mathcal{A} is decelerated away, the number of modes increases exponentially: Note that for dσdτ<0\frac{d\sigma}{d\tau}<0, the modes will be annihilated, while for dσdτ>0\frac{d\sigma}{d\tau}>0, the modes will be created in time.

III.5 The Emissivity Factor

Let us return to the Planck’s law at an accelerated frame (50). The scale factor ε[α,ξ]\varepsilon\left[\alpha,\xi\right] could be interpreted physically as the emissivity resulted from the acceleration of the observer. The emissivity is the factor that describes the imperfectness of a black-body, i.e., a physical property that describes how efficiently an object emits thermal radiation. The classical range of the emissivity is 0ε10\leq\varepsilon\leq 1, where 1 is the emissivity of a black-body and 0 is the emissivity of a perfect thermal mirror, where no radiation is emitted. In this interval, are the region of the grey-body. Classically, ε>1\varepsilon>1 and ε<0\varepsilon<0 are not defined, except for special cases for particles smaller than the dominant radiation wavelength (Golyk, ).

There are two factors that affect the emissivity due to acceleration as in (51), namely, the spatial coordinate ξ\xi, and the magnitude of the acceleration α\alpha. Let us consider the case where ς¯=ς±ατ\bar{\varsigma}=\varsigma\pm\alpha\tau, namely the case where the observer are respectively, accelerated and decelerated away from the source. One could derive that for these cases, the scale factor is:

ε[α,ξ]=e2αξ(1α),\varepsilon\left[\alpha,\xi\right]=e^{2\alpha\xi}\left(1\mp\alpha\right), (66)

where the minus sign is for the accelerated away case, and the plus for the decelerated. For the accelerated away case, the region with negative emissivity is obtained when α>1\alpha>1; at this condition, no radiation will be detected by the accelerated observer. Meanwhile, positive emissivity factor is obtained when 0α10\leq\alpha\leq 1. For this case, if ξ>0,\xi>0, it is possible to have the emissivity ε>1\varepsilon>1, hence the black-body will be perceived as “hyperblack” in the accelerated frame; otherwise, for ξ<0\xi<0, it will be perceived as grey. See FIG. 13.

Refer to caption
Figure 13: The plot of the source’s emissivity ε\varepsilon with respect to the acceleration magnitude α\alpha, for different value of spatial position ξ\xi. There is a strict threshold α1\alpha\geq 1 where there is no radiation detected by the accelerated observer due to ε0\varepsilon\leq 0, independent from the value of ξ\xi. For 0α10\leq\alpha\leq 1, the source could be perceived as grey, black, or hyperblack, depending on ξ\xi. The dashed-line is the condition of the black-body with ε=1\varepsilon=1; above and below the line are, respectively, the ’hyperblack’ and grey regions. For ξ=0\xi=0, the relation between emissivity and acceleration becomes linear: ε=1α\varepsilon=1-\alpha.

For the decelerated away case, the emissivity factor ε0\varepsilon\geq 0, so the accelerated frame will observe a range of hyperblack, black, and greybody, depending on the value of α\alpha and ξ\xi. See FIG. 14.

Refer to caption
Figure 14: The plot of the source’s emissivity ε\varepsilon with respect to the deceleration magnitude α\alpha, for different value of spatial position ξ\xi. The emissivity for this case is positive definite: ε0\varepsilon\geq 0, hence the observer will detect radiation for any value of α\alpha. The source could be perceived as grey, black, or hyperblack, depending on ξ\xi. The dashed-line is the condition of the black-body with ε=1\varepsilon=1, above and below the line are, respectively, the ’hyperblack’ and grey region. For ξ=0\xi=0, the relation between emissivity and acceleration becomes linear: ε=1+α\varepsilon=1+\alpha.

While in the decelerated away case one can perceive the source with different ’blackness’ depending on the value of α\alpha and ξ\xi, an interesting phenomenon occurs in the case of an accelerated away scenario. For this case, if 0α10\leq\alpha\leq 1 is satisfied, the emissivity is 0ε10\leq\varepsilon\leq 1; the black-body is perceived as a greybody in 𝒜\mathcal{A}, where the signal becomes dimmer as α1\alpha\rightarrow 1. At α=1\alpha=1, ε=0\varepsilon=0: in the frame of 𝒜\mathcal{A}, the black-body stops emitting radiation as if it is a perfect thermal mirror. From this point, the emissivity can only be negative definite as α\alpha increases. Remarkably, for a constant ξ\xi, the emissivity is independent of time and the position of the observer from the source, which means that no matter how far or close the object is from the observer, theoretically, they will measure the same emissivity. There could be several possible explanations for this phenomenon, which might include the relativistic beaming and the existence of the (Killing) horizon at x=±tx=\pm t, however, further research is needed to understand this phenomenon completely.

III.6 Wien’s Displacement Law and Temperature

At the moment, we are not able to derive the transformation of Wien’s displacement law in a moving frame without knowing how the temperature TT in relation (76) transform under Lorentz transformation. However, if we assume that Wien’s law is valid also in a moving frame, we can define the directional temperature using the Wien’s law (76), simply as a quantity proportional to the frequency that maximies the energy of the source. This, at least, will be useful for the calculation of a temperature for a black-body in moving frames. As we had mentioned in the Introduction, any measurement of temperature in moving bodies are measured via it’s radiation, since we are not able to do a direct measurement of the temperature in moving bodies. If a body is moving, it needs to be at rest with respect to an observer so that it could be thermalize/ in a thermal equilibrium with the measurement apparatus. Here, satisfied or not, the temperature for moving bodies, is treated as a derived variables. Moreover, if we make an assumption that (76) is valid for inertial and moreover, uniformly accelerated frame, then the directional temperature TT must transform in the same way as the frequency, namely:

T¯=Teαξcoshς¯[τ]+cosθ¯[τ]sinhς¯[τ].\bar{T}=\frac{Te^{-\alpha\xi}}{\cosh\bar{\varsigma}\left[\tau\right]+\cos\bar{\theta}\left[\tau\right]\sinh\bar{\varsigma}\left[\tau\right]}. (67)

(67) is a generalization of the ’directional’-(effective) temperature in (Henry, ) which is valid in any inertial frame; it is an extension of the result in (Henry, ) to a uniformly accelerated frame. One could observe that in (67), there is a coordinate-dependent scale factor proportional to eαξ,\sim e^{-\alpha\xi}, similar to the temperature scale factor in Tolman’s theory of general relativistic thermodynamics (Tolman, ).

IV Discussion and Conclusion

IV.1 Discussion

The Planck’s formula in the accelerated frame (50) are based on several assumptions in the derivation, hence, if one or more of these assumptions are not valid, this will also affect the validity of (50). Let us list the assumptions we use for the derivation:

  1. 1.

    E=nωE=n\hbar\omega. This is implicitly used because we started from the Planck’s law in rest frame (the relation is use to derive Planck’s law).

  2. 2.

    The validity of the de Broglie postulate p=kp=\hbar k in the accelerated frame, namely: p[τ,ξ]=k[τ,ξ]p\left[\tau,\xi\right]=\hbar k\left[\tau,\xi\right], for every Rindler time τ\tau and Rindler position ξ\xi.

  3. 3.

    The invariance of the relativistic distribution function (number of world-lines that crosses a hypersurface) under coordinate transformation.

Starting from these assumptions, together with some standard definitions in the theory of relativity, all the equations we derived, most importantly, the relativistic Doppler effect (29), the transformation of solid angle (33), the volume contraction (36), and the transformation of the number of modes distribution (48), are consequences of the assumptions and definitions we used. Let us check if all these assumptions are reasonable. The first assumption is used to derive Planck’s law in the original form, and since it is not use explicitly to derive the results, this should not concern our work. The second assumption is used explicitly to derive the dispersion relation (20) from the momentum constraint pμpμ=0.p_{\mu}p^{\mu}=0. (20) is crucial for deriving the momentum in 𝒜,\mathcal{A}, namely (23). However, one can also obtained the dispersion relation (20) without the de Broglie postulate, if we consider the electromagnetic wave equation 𝐄=0\boxempty\mathbf{E}=0, with \boxempty is the d’Alembertian operator. The form of the wave equation and its solution are invariant under Rindler coordinate transformation, and so does the dispersion relation. Therefore, without the de Broglie postulate, (20) is still valid in 𝒜\mathcal{A} as long as the Maxwell equation is satisfied by the black-body radiation.

Finally, opinions on the last assumption are divided among the researchers in the field. The relativistic distribution function f(𝐱,𝐩)f\left(\mathbf{x},\mathbf{p}\right) is defined as (47). Since our work follows closely the work of (Peebles, ; Heer, ; Henry, ), we use the same assumption that the relativistic distribution function f(𝐱,𝐩)f\left(\mathbf{x},\mathbf{p}\right) is coordinate-invariant, instead of the invariance of the number of states dNdN, which is also employed in many existing works in the literature. The reason for this difference is because there are 2 distinct measurements of volumes in moving frames: the ones that measure all the points in the volume simultaneously at an instant time, and the one that measure all the points in the volume at a time interval Δτ=τfτi\Delta\tau=\tau_{f}-\tau_{i}. Our work use the second type of measurement, and hence the number of states counted by the ’moving volume’ at time interval Δτ\Delta\tau is not the same as dNdN (the number of states in rest frame); hence we use assumption 3.

Some readers may wonder if a simple coordinate transformation from inertial to Rindler coordinate is sufficient to obtain the physics in an accelerated frame. As far as we understand, in the theory of relativity (special or general), the use of different coordinates on spacetime, in general, will result in different observational perspectives. This is not only valid for coordinates related by Lorentz transformation, but also for more general coordinate systems, including the Rindler coordinate. Therefore, to obtain observational perspective from an accelerated frame, it is sufficient to do a transformation from inertial/Cartesian coordinate to Rindler coordinate. Furthermore, one might also wonder why the Unruh effect is not included in our work. This is out of the scope of the subject in our work, since the Unruh effect is derived from the quantum field theoretic derivation, and we have tried to write our paper as classical as possible. The reason to write the paper as classical as possible is because the subject is full of controversy, so in our opinion, it will be an advantage if we could, at least, clearly understand the subject in the (semi)-classical level. Although it might be possible to include the Unruh effect in the discussion, in the classical level, it can not be obtained just by coordinate transformation from Minkowski/Cartesian coordinate to Rindler coordinate. There is another assumption needed to obtain the Unruh effect in the classical setting, and this is related to the definition of the classical vacuum. This will be an insteresting subject to pursue.

IV.2 Conclusion

Let us conclude our work in this article as follows. We have derived Planck’s law and calculated the black-body spectrum in a uniformly accelerated frame using Rindler coordinates. The spectrum is time-dependent, Planckian at each instantaneous time, but it is scaled by an emissivity factor ε=eαξ(1±α)\varepsilon=e^{\alpha\xi}\left(1\pm\alpha\right) that depends on the Rindler spatial coordinate ξ\xi and the acceleration magnitude α\alpha. The spatial, coordinate-dependent scale factor is proportional to eαξ,\sim e^{\alpha\xi}, while the scale factor related to the acceleration is proportional to eαξ(1±α)\sim e^{\alpha\xi}\left(1\pm\alpha\right), depending on the observer is either accelerated or decelerated. An observer decelerating away from the source will perceive the black-body as hyperblack, black, or grey, while for the accelerating-away observer, there is a limit in the acceleration magnitude in receiving the radiation, namely, if α1\alpha\geq 1, the accelerated observer will stop receiving radiation from the black-body. Outside of this limit, the black-body is perceived either as grey or hyperblack, depending on the spatial Rindler coordinate ξ\xi.

The variables in the Planck’s law, specifically the number of modes and frequency, evolve over time. The Planckian spectrum is continuously red-shifted towards lower frequencies as time progresses for the case where the observer is accelerating-away or decelerating towards, and blue-shifted towards lower (or higher) frequencies for the case where the observer is accelerating-towards or decelerating-away. In the accelerated frame, the production of radiation modes (photons) can be positive or negative, depending on the acceleration or deceleration of the observer, and zero for vanishing acceleration. Besides this, there exists a peculiar phenomenon perceived by the accelerating observer, the relativistic beaming, where the rays of light tends to converge, or diverge in the direction of motion, depending on whether the observer is moving towards, or away from the source. In the end, assuming Wien’s displacement law also holds in the accelerated frame, the time-dependent, (directional)-temperature of a body within an accelerated frame is given by T¯[τ]=eαξ(coshς¯[τ]+cosθ¯[τ]sinhς¯[τ])1T.\bar{T}\left[\tau\right]=e^{-\alpha\xi}\left(\cosh\bar{\varsigma}\left[\tau\right]+\cos\bar{\theta}\left[\tau\right]\sinh\bar{\varsigma}\left[\tau\right]\right)^{-1}T.

Our last comment is on the temperatures of the systems. For the moment, we do not have a universal definition of temperature. However, the (effective) temperature (67) is derived from their frequency spectra, this is also the case for the Hawking-Unruh temperature. Since for a moving body we do not have other choice than treating (effective) temperature as a derived quantity, we argue that one needs to take seriously the directional approach of temperature. This is supported by the argument in (Aldrovandi, ) that the inertial version of effective temperature (67) is not merely a mathematical parameter, but a real transformation law. We are optimistic that this will be a useful approach to understanding the nature of temperature.

Acknowledgements.
S. A. was supported by an appointment to the Young Scientist Training Program at the Asia Pacific Center for Theoretical Physics (APCTP) through the Science and Technology Promotion Fund and Lottery Fund of the Korean Government. This was also supported by the Korean Local Governments - Gyeongsangbuk-do Province and Pohang City. H. L. P. would like to thank Ganesha Talent Assistanship (GTA) Institut Teknologi Bandung for financial support.

Appendix A Black-Body Radiation in Rest Frame

A.1 Planck’s Radiation Law

To make this article self-contained, in this subsection of the Appendix A, we derive the Planck’s law in the form of relation (49) from the original equation (Planck2, ):

B[f,T]=2hf3c21ehf/(kBT)1.B\left[f,T\right]=\frac{2hf^{3}}{c^{2}}\frac{1}{e^{\nicefrac{{hf}}{{\left(k_{B}T\right)}}}-1}. (68)

BB is the intensity or spectral radiance, defined as:

B[f]=dEcosθdfdAdΩdt,B\left[f\right]=\frac{dE}{\cos\theta dfdAd\Omega dt}, (69)

with ff is the frequency of the electromagnetic radiation and TT is the temperature of the black-body source emitting the radiation. h,h, cc, and kBk_{B} are respectively, the Planck constant, the speed of light, and the Boltzmann constant. The intensity B[f]B\left[f\right] is defined as the infinitesimal energy dEdE that passes through an infinitesimal area dAdΩdA\,d\Omega of the surface of a sphere, per time interval dtdt, in the frequency range [f,f+df]\left[f,f+df\right]. θ\theta is the inclination angle: the angle between the velocity of the photon (that is parallel to the normal to dAdA) with the direction of the observation, but with frame 𝒜\mathcal{A} replaced with frame 𝒪\mathcal{O}^{\prime}. Inserting (69) to (68), and multiplying the right hand side of the equation with cc=1\frac{c}{c}=1 gives:

dE=2hf3c31ehf/(kBT)1dfdΩdAccosθdt.dE=\frac{2hf^{3}}{c^{3}}\frac{1}{e^{\nicefrac{{hf}}{{\left(k_{B}T\right)}}}-1}dfd\Omega dA\,c\cos\theta dt. (70)

The quantity dAccosθdtdA\,c\cos\theta dt is the infinitesimal volume dVdV swept out by the radiation. Using the angular frequency ω=2πf\omega=2\pi f instead of ff, the Planck’s law can be written as:

dE=ω34π3c31eω/(kBT)1ρE[ω]dωdΩdV.dE=\underset{\rho_{E}\left[\omega\right]}{\underbrace{\frac{\hbar\omega^{3}}{4\pi^{3}c^{3}}\frac{1}{e^{\nicefrac{{\hbar\omega}}{{\left(k_{B}T\right)}}}-1}}}d\omega d\Omega dV. (71)

Moreover, given the energy of each frequency mode ω\omega as dE=ωdNdE=\hbar\omega dN, with NN is the number (or distribution) of the radiation modes (quanta of radiation/photon, however, we avoid such terminologies because we want this paper to be as classical as possible) having an angular frequency ω\omega, we could obtain the Planck’s distribution:

dN=ω24π3c31eω/(kBT)1dωdΩdV,dN=\frac{\omega^{2}}{4\pi^{3}c^{3}}\frac{1}{e^{\nicefrac{{\hbar\omega}}{{\left(k_{B}T\right)}}}-1}d\omega d\Omega dV, (72)

with the famous Planckian spectrum. (72) is the form of Planck’s law used in (Henry, ).

A.1.1 Wien’s Displacement Law

To derive Wien’s displacement law from the Planck’s law (71), one needs to find the angular frequency ωmax\omega_{\mathrm{max}} that maximize the radiation energy EE. This can be obtained by solving the equation dρE[ω]dω=0\frac{d\rho_{E}\left[\omega\right]}{d\omega}=0 as follows:

4π3c3ω2eω/(kBT)1(3ωkBTeω/(kBT)eω/(kBT)1)=0.\frac{\hbar}{4\pi^{3}c^{3}}\frac{\omega^{2}}{e^{\nicefrac{{\hbar\omega}}{{\left(k_{B}T\right)}}}-1}\left(3-\frac{\hbar\omega}{k_{B}T}\frac{e^{\nicefrac{{\hbar\omega}}{{\left(k_{B}T\right)}}}}{e^{\nicefrac{{\hbar\omega}}{{\left(k_{B}T\right)}}}-1}\right)=0. (73)

Let us define w=ωkBTw=\frac{\hbar\omega}{k_{B}T}, the equation from (73) that we need to solve becomes:

w=3(1ew),w=3\left(1-e^{-w}\right), (74)

which can be solved by the Lambert function W0W_{0} as follows:

w=3+W0(3e3).w=3+W_{0}\left(-3e^{-3}\right). (75)

With this, one could obtain the Wien’s displacement law in terms of angular frequency:

kBT=1(3+W0(3e3))ωmax,k_{B}T=\frac{1}{\left(3+W_{0}\left(-3e^{-3}\right)\right)}\hbar\omega_{\mathrm{max}}, (76)

with ωmax\omega_{\mathrm{max}} is the angular frequency that maximize the radiation energy EE, the ’maximal’ frequency. The statement that the temperature of a black-body source is proportional to the frequency that maximizes the energy of the radiation comes directly from (76).

One can insert the Wien’s law (76) to (72) and obtain:

dN=ω24π3c31ewω/ωmax1dωdΩdV.dN=\frac{\omega^{2}}{4\pi^{3}c^{3}}\frac{1}{e^{\nicefrac{{w\omega}}{{\omega_{\mathrm{max}}}}}-1}d\omega d\Omega dV. (77)

We called the dimensionless coefficient ww satisfying (74) as the Wien constant. Notice that for a different representation of the Planck’s law, i.e., if the law written as a function of the wavelength λ\lambda instead of ω\omega, the Wien ’constant’ will differ.

The Planck’s law in the form of equation (72) is equivalent with the two equations: the “Planck’s law” in terms of maximal frequency (77), together with the Wien’s law (76). To avoid the problem of the temperature in moving bodies (Farias, ; Derakshani, ), we will use the “Planck’s law” (77) in our analysis for the black-body spectrum in a uniformly accelerated frame.

Appendix B Black-Body Radiation in Inertial Frames

B.1 Lorentz Transformation

In this subsection, we review the basic properties of Lorentz transformation and define the notations used in our article. Let 𝒪\mathcal{O} be a frame at rest, with spatial coordinates (x,y,z)\left(x,y,z\right) and time coordinate 𝔱.\mathfrak{t}. In the covariant approach, the time and space are regarded in an equal manner, so let us define the ’time’ t:=c𝔱t:=c\mathfrak{t} such that tt has the same dimension as length; cc is the speed of light. Let 𝐱\mathbf{x} =(x,y,z,t)\left(x,y,z,t\right) be an inertial coordinate that parametrizes the Minkowski space.

Suppose another inertial frame 𝒪\mathcal{O}^{\prime} is moving with a constant velocity uu in the xx-direction. Let 𝐱=(x,y,z,t)\mathbf{x}^{\prime}=\left(x^{\prime},y^{\prime},z^{\prime},t^{\prime}\right) be another inertial coordinate of the Minkowski space related to frame 𝒪\mathcal{O}^{\prime}. Frames 𝒪\mathcal{O} and 𝒪\mathcal{O}^{\prime} are related by the coordinate (Lorentz) transformation as follows:

t=γ(tβx),x=γ(xβt),t=γ(t+βx),x=γ(x+βt),\begin{array}[]{cc}t^{\prime}&=\gamma\left(t-\beta x\right),\\ x^{\prime}&=\gamma\left(x-\beta t\right),\end{array}\quad\begin{array}[]{cc}t&=\gamma\left(t^{\prime}+\beta x^{\prime}\right),\\ x&=\gamma\left(x^{\prime}+\beta t^{\prime}\right),\end{array} (78)

where y,zy,z are not affected by the transformation, namely y=yy^{\prime}=y and z=z.z^{\prime}=z. Here, we use β=uc\beta=\frac{u}{c} and γ=11β2.\gamma=\frac{1}{\sqrt{1-\beta^{2}}}. Let us define the rapidity as ς=cosh1γ\varsigma=\cosh^{-1}\gamma, then we have the following relations:

γ=coshς,βγ=sinhς,β=tanhς.\gamma=\cosh\varsigma,\;\beta\gamma=\sinh\varsigma,\;\beta=\tanh\varsigma. (79)

With relations (79), the Lorentz transformation (78) can be written in terms of hyperbolic functions:

t=tcoshςxsinhς,x=xcoshςtsinhς,t=tcoshς+xsinhς,x=xcoshς+tsinhς.\begin{array}[]{cc}t^{\prime}&=t\cosh\varsigma-x\sinh\varsigma,\\ x^{\prime}&=x\cosh\varsigma-t\sinh\varsigma,\end{array}\;\begin{array}[]{cc}t&=t^{\prime}\cosh\varsigma+x^{\prime}\sinh\varsigma,\\ x&=x^{\prime}\cosh\varsigma+t^{\prime}\sinh\varsigma.\end{array} (80)

The infinitesimal (4-vector) transformation related to (80) can be obtained as follows:

dt=γ(dtβdx)=dtcoshςdxsinhςdx=γ(dxβdt)=dxcoshςdtsinhς,\begin{array}[]{ccc}dt^{\prime}&=\gamma\left(dt-\beta dx\right)&=dt\cosh\varsigma-dx\sinh\varsigma\\ dx^{\prime}&=\gamma\left(dx-\beta dt\right)&=dx\cosh\varsigma-dt\sinh\varsigma\end{array}, (81)

together with their inverses:

dt=γ(dt+βdx)=dtcoshς+dxsinhςdx=γ(dx+βdt)=dxcoshς+dtsinhς,\begin{array}[]{ccc}dt&=\gamma\left(dt^{\prime}+\beta dx^{\prime}\right)&=dt^{\prime}\cosh\varsigma+dx^{\prime}\sinh\varsigma\\ dx&=\gamma\left(dx^{\prime}+\beta dt^{\prime}\right)&=dx^{\prime}\cosh\varsigma+dt^{\prime}\sinh\varsigma\end{array}, (82)

with dy=dydy^{\prime}=dy and dz=dzdz^{\prime}=dz.

B.2 The Relativistic Aberration of Light

Suppose with respect to the rest frame 𝒪\mathcal{O} we have an object moving with velocity 𝐯\mathbf{v} as follows:

𝐯=(vx,vy,vz)=(dxdt,dydt,dzdt).\mathbf{v}=\left(v_{x},v_{y},v_{z}\right)=\left(\frac{dx}{dt},\frac{dy}{dt},\frac{dz}{dt}\right).

In frame 𝒪\mathcal{O}^{\prime}, the velocity of the object is:

𝐯=(vx,vy,vz)=(dxdt,dydt,dzdt).\mathbf{v}^{\prime}=\left(v^{\prime}_{x},v^{\prime}_{y},v^{\prime}_{z}\right)=\left(\frac{dx^{\prime}}{dt^{\prime}},\frac{dy^{\prime}}{dt^{\prime}},\frac{dz^{\prime}}{dt^{\prime}}\right).

Using the transformation (81)-(82), one could obtain the relation between 𝐯\mathbf{v} and 𝐯\mathbf{v}^{\prime} as follows

vx=vx+β1+βvx=vx+tanhς1+vxtanhς,v_{x}=\frac{v^{\prime}_{x}+\beta}{1+\beta v^{\prime}_{x}}=\frac{v^{\prime}_{x}+\tanh\varsigma}{1+v^{\prime}_{x}\tanh\varsigma}, (83)
vy,z=vy,zγ(1+βvx)=vy,zcoshς(1+vxtanhς),v_{y,z}=\frac{v^{\prime}_{y,z}}{\gamma\left(1+\beta v^{\prime}_{x}\right)}=\frac{v^{\prime}_{y,z}}{\cosh\varsigma\left(1+v^{\prime}_{x}\tanh\varsigma\right)}, (84)

and their inverses:

vx=vxβ1βvx=vxtanhς1vxtanhς,v^{\prime}_{x}=\frac{v_{x}-\beta}{1-\beta v{}_{x}}=\frac{v_{x}-\tanh\varsigma}{1-v_{x}\tanh\varsigma}, (85)
vy,z=vy,zγ(1βv)x=vy,zcoshς(1vxtanhς).v^{\prime}_{y,z}=\frac{v_{y,z}}{\gamma\left(1-\beta v{}_{x}\right)}=\frac{v_{y,z}}{\cosh\varsigma\left(1-v_{x}\tanh\varsigma\right)}. (86)

Let a black-body source be at rest with respect to an inertial frame 𝒪.\mathcal{O}. The black-body are emitting electromagnetic radiation (photons) with propagation velocity cc in the direction of n^=(cosθ,sinθ,0)\hat{n}=\left(\cos\theta,\sin\theta,0\right), with θ\theta is the (polar) inclination angle between the axis +xx and n^\hat{n} on plane dxdydx\wedge dy, but with frame 𝒜\mathcal{A} replaced with frame 𝒪\mathcal{O}^{\prime}. The electric and magnetic part of the radiation lie, respectively, in plane dxdydx\wedge dy and n^dz\hat{n}\wedge dz, so the propagation velocity will not affect their amplitude. The velocity of the photon with respect to frame 𝒪\mathcal{O} is:

𝐯=(vx,vy,vz)=cn^=(ccosθ,csinθ,0).\mathbf{v}=\left(v_{x},v_{y},v_{z}\right)=c\hat{n}=\left(c\cos\theta,c\sin\theta,0\right). (87)

Meanwhile, at frame 𝒪\mathcal{O}^{\prime} , the velocity of the photon is:

𝐯=(vx,vy,vz)=cn^=(ccosθ,csinθ,0),\mathbf{v}^{\prime}=\left(v^{\prime}_{x},v^{\prime}_{y},v^{\prime}_{z}\right)=c\hat{n}^{\prime}=\left(c\cos\theta^{\prime},c\sin\theta^{\prime},0\right), (88)

with n^\hat{n}^{\prime} and θ\theta^{\prime} are the propagation direction and the inclination angle of the photon according to 𝒪\mathcal{O}^{\prime}, respectively. This could be derived from the null (light-like) vector condition for light, where the norm of its 4-velocity 𝐮=(γ,γ𝐯)\mathbf{u}=\left(\gamma,\gamma\mathbf{v}\right) is always zero in any coordinate system. This gives:

uαuα=γ2(𝐯21)=0,u^{\alpha}u_{\alpha}=\gamma^{2}\left(\left\|\mathbf{v}\right\|^{2}-1\right)=0,

where the 3-velocities satisfy (87) and (88).

Using the velocity addition formulas in (83)-(85), one could obtain the transformation between velocities of the photon seen by 𝒪\mathcal{O} and 𝒪\mathcal{O}^{\prime}, which results in the polar angle formulas as follows:

cosθ\displaystyle\cos\theta^{\prime} =cosθβ1βcosθ,\displaystyle=\frac{\cos\theta-\beta}{1-\beta\cos\theta}, (89)
sinθ\displaystyle\sin\theta^{\prime} =sinθγ(1βcosθ).\displaystyle=\frac{\sin\theta}{\gamma\left(1-\beta\cos\theta\right)}. (90)

Inserting (89) and (90) to a trigonometric identity, tanθ2=sinθ1+cosθ\tan\frac{\theta^{\prime}}{2}=\frac{\sin\theta^{\prime}}{1+\cos\theta^{\prime}}, gives:

tanθ2=1γ(1β)tanθ2=1+β1βtanθ2.\tan\frac{\theta^{\prime}}{2}=\frac{1}{\gamma\left(1-\beta\right)}\tan\frac{\theta}{2}=\frac{\sqrt{1+\beta}}{\sqrt{1-\beta}}\tan\frac{\theta}{2}. (91)

Equation (89)-(91) are the aberration of light formulas (Johnson, ). Their inverses are:

cosθ\displaystyle\cos\theta =cosθ+β1+βcosθ,\displaystyle=\frac{\cos\theta^{\prime}+\beta}{1+\beta\cos\theta}, (92)
sinθ\displaystyle\sin\theta =sinθγ(1+βcosθ),\displaystyle=\frac{\sin\theta^{\prime}}{\gamma\left(1+\beta\cos\theta^{\prime}\right)}, (93)
tanθ2\displaystyle\tan\frac{\theta}{2} =1β1+βtanθ2.\displaystyle=\frac{\sqrt{1-\beta}}{\sqrt{1+\beta}}\tan\frac{\theta^{\prime}}{2}. (94)

B.3 The Relativistic Doppler Effect

The next step is to derive the relativistic Doppler effect. First, we need to define the 4-momentum 𝐩\mathbf{p} of a moving body with respect to 𝒪\mathcal{O} as follows:

𝐩=(pt,px,py,pz𝐩3),\mathbf{p}=\left(p_{t},\underset{{}^{3}\mathbf{p}}{\underbrace{p_{x},p_{y},p_{z}}}\right), (95)

with pt=E/cp_{t}=\nicefrac{{E}}{{c}} is the energy of the moving body and 𝐩3{}^{3}\mathbf{p} is the relativistic 3-momentum. For the black-body radiation case, the moving bodies are photons, which, by de Broglie postulate, has energy E=ωE=\hbar\omega and momentum 𝐩3=𝐤{}^{3}\mathbf{p}=\hbar\mathbf{k}. Here, =h/2π\hbar=\nicefrac{{h}}{{2\pi}} is the Planck constant, ω=2πf\omega=2\pi f is the (angular) frequency, and 𝐤=(kx,ky,kz)\mathbf{k}=\left(k_{x},k_{y},k_{z}\right) is the wave vector satisfying 𝐤=2πn^λ.\mathbf{k}=\frac{2\pi\hat{n}}{\lambda}. Inserting these information to (95) together with the components of n^\hat{n}, one could obtain the 4-momentum of a photon with respect to an inertial observer 𝒪\mathcal{O}:

𝐩=ωc(1,n^)=ωc(1,cosθ,sinθ,0).\mathbf{p}=\frac{\hbar\omega}{c}\left(1,\hat{n}\right)=\frac{\hbar\omega}{c}\left(1,\cos\theta,\sin\theta,0\right). (96)

According to the moving frame 𝒪\mathcal{O}^{\prime}, the 4-momentum of the photon is:

𝐩=(pt,px,py,pz)=ωc(1,cosθ,sinθ,0),\mathbf{p}^{\prime}=\left(p_{t}^{\prime},p_{x}^{\prime},p_{y}^{\prime},p_{z}^{\prime}\right)=\frac{\hbar\omega^{\prime}}{c}\left(1,\cos\theta^{\prime},\sin\theta^{\prime},0\right), (97)

with ω\omega^{\prime} is the frequency of the photon, with respect to 𝒪\mathcal{O}^{\prime}. One could also obtain (96) and (97) by inserting the de Broglie postulate to the zero mass condition for the photon:

pμpμ=E2+𝐩32=0;p^{\mu}p_{\mu}=-E^{2}+\left\|{}^{3}\mathbf{p}\right\|^{2}=0;

with this, we only consider the positive solution from the dispersion relation ω2=c2k2.\omega^{2}=c^{2}k^{2}.

The 4-momentum is an element of the Minkowski space and therefore transform under Lorentz transformation, so the relation of 𝐩\mathbf{p}^{\prime} and 𝐩\mathbf{p} is:

pt=γ(ptβpx),px=γ(pxβpt),pt=γ(pt+βpx),px=γ(px+βpt),\begin{array}[]{cc}p_{t}^{\prime}&=\gamma\left(p_{t}-\beta p_{x}\right),\\ p_{x}^{\prime}&=\gamma\left(p_{x}-\beta p_{t}\right),\end{array}\quad\begin{array}[]{cc}p_{t}&=\gamma\left(p_{t}^{\prime}+\beta p_{x}^{\prime}\right),\\ p_{x}&=\gamma\left(p_{x}^{\prime}+\beta p_{t}^{\prime}\right),\end{array} (98)

with py=pyp_{y}^{\prime}=p_{y} and pz=pzp_{z}^{\prime}=p_{z}. Naturally, the (generalized) momentum is a co-vector instead of a (contravariant) vector, however, here we use the vector transformation of momentum. This will not affect the result as long as the calculations are consistent. Using the Lorentz transformation (98) on the momenta (96) and (97), one could obtain the relations between the frequency of photon ω\omega as seen by 𝒪\mathcal{O} and ω\omega^{\prime} as seen by 𝒪\mathcal{O}^{\prime}, written in 3 equivalent forms as follows (Johnson, ):

ωω\displaystyle\frac{\omega}{\omega^{\prime}} =sinθsinθ,\displaystyle=\frac{\sin\theta^{\prime}}{\sin\theta}, (99)
ωω\displaystyle\frac{\omega}{\omega^{\prime}} =γ(1+βcosθ),\displaystyle=\gamma\left(1+\beta\cos\theta^{\prime}\right), (100)
ωω\displaystyle\frac{\omega}{\omega^{\prime}} =γ(cosθ+β)cosθ.\displaystyle=\frac{\gamma\left(\cos\theta^{\prime}+\beta\right)}{\cos\theta}.

From these equations, one could retrieve the aberration formulas (89)-(91) and their inverses (92)-(94). Inserting (89) to (100), one could rewrite ωω\frac{\omega}{\omega^{\prime}} in terms of only variable θ\theta, the inclination angle seen by observer 𝒪\mathcal{O}:

ωω\displaystyle\frac{\omega}{\omega^{\prime}} =1γ(1βcosθ).\displaystyle=\frac{1}{\gamma\left(1-\beta\cos\theta\right)}. (101)

Equation (100) and (101) are the formula describing the relativistic Doppler effect for the frequency of the photon in a moving observer 𝒪\mathcal{O}^{\prime} with respect to observer at rest 𝒪\mathcal{O} (Johnson, ). For our case, the photon frequency observed by the moving observer 𝒪\mathcal{O}^{\prime} is redshifted, since observer 𝒪\mathcal{\mathcal{O}} is moving away from the black-body source. To obtain the blueshifted Doppler effect, one could flip the sign on the velocity of 𝒪\mathcal{O}^{\prime} with respect to the source to become u-u.

Since the Planck’s law (77) contains the solid angle term dΩd\Omega, we need to know how it transforms under Lorentz transformation. First, from (100) and (101) we obtain:

γ(1+βcosθ)=1γ(1βcosθ).\gamma\left(1+\beta\cos\theta^{\prime}\right)=\frac{1}{\gamma\left(1-\beta\cos\theta\right)}.

Second, differentiating equation (90) or (93) with respect to a parameter will give:

dθdθ=ωω.\frac{d\theta^{\prime}}{d\theta}=\frac{\omega}{\omega^{\prime}}. (102)

With (102), the solid angle dΩ=sinθdθdϕd\Omega=\sin\theta d\theta d\phi, will transform as follows (Johnson, ):

dΩdΩ=(ωω)2,\frac{d\Omega^{\prime}}{d\Omega}=\left(\frac{\omega}{\omega^{\prime}}\right)^{2}, (103)

using (99), (102), and the fact that the azimuth angle are not affected by the transformation, dϕ=dϕd\phi=d\phi^{\prime}.

B.4 The Phase-Space Volume Transformation

The next variable in the Planck’s law that transformed under the Lorentz transformation is NN, the number or distribution of modes with frequency ω\omega. To know how this variable changes under different inertial frames, one needs to consider the phase-space volume transformation. The phase-space volume element is not a Lorentz scalar, see a detailed explanation in (Debbasch, ). First, let us obtain the transformation for the 3-volume element of the spatial part of the phase-space. In 𝒪\mathcal{O} and 𝒪\mathcal{O}^{\prime}, the infinitesimal 3-volume element are defined as, respectively:

d3𝐱\displaystyle d^{3}\mathbf{x} =dxdydz,\displaystyle=dx\wedge dy\wedge dz, (104)
d3𝐱\displaystyle d^{3}\mathbf{x}^{\prime} =dxdydz.\displaystyle=dx^{\prime}\wedge dy^{\prime}\wedge dz^{\prime}.

To derive the transformation relation between these 2 elements, let us consider a measurement of spatial length in 𝒪.\mathcal{O}. To measure spatial length, 2 events must be simultaneous in time with respect to the observer. Let us consider 2 simultaneous events 𝐩=(ti,xi)\mathbf{p}=\left(t_{i},x_{i}\right) and 𝐪=(tf,xf)\mathbf{q}=\left(t_{f},x_{f}\right) along the xx-axis at time ti=tf=0t_{i}=t_{f}=0, where xi=0x_{i}=0 and xf=x_{f}=\ell. According to 𝒪,\mathcal{O}, the length between these 2 event is:

Δx\displaystyle\Delta x =xfxi=.\displaystyle=x_{f}-x_{i}=\ell.

Now, let these 2 events be perceived in frame 𝒪\mathcal{O}^{\prime}, moving with respect to 𝒪.\mathcal{O}. Using Lorentz transformation (78):

ti=γ(tiβxi)=0,xi=γ(xiβti)=0,tf=γ(tfβxf)=βγ,xf=γ(xfβtf)=γ;\begin{array}[]{cc}t^{\prime}_{i}&=\gamma\left(t_{i}-\beta x_{i}\right)=0,\\ x^{\prime}_{i}&=\gamma\left(x_{i}-\beta t_{i}\right)=0,\end{array}\;\begin{array}[]{cc}t^{\prime}_{f}&=\gamma\left(t_{f}-\beta x_{f}\right)=-\beta\gamma\ell,\\ x^{\prime}_{f}&=\gamma\left(x_{f}-\beta t_{f}\right)=\gamma\ell;\end{array}

it is clear that in 𝒪\mathcal{O}^{\prime} these 2 events are not simultaneous. The spatial length between 𝐩\mathbf{p} and 𝐪\mathbf{q} in 𝒪\mathcal{O}^{\prime} is:

Δx\displaystyle\Delta x^{\prime} =xfxi=γ.\displaystyle=x^{\prime}_{f}-x^{\prime}_{i}=\gamma\ell.

Now for our black-body case, we need to calculate how much modes are inside a volume element. Let the finite volume element in 𝒪\mathcal{O} be ΔxΔyΔz\Delta x\Delta y\Delta z. This volume will be perceived simultaneously by 𝒪\mathcal{O} at an instant time tt, while in 𝒪,\mathcal{O}^{\prime}, it will be perceived as the volume swept by the plane ΔyΔz=ΔyΔz\Delta y^{\prime}\Delta z^{\prime}=\Delta y\Delta z from xix^{\prime}_{i} to xfx^{\prime}_{f} along a time interval Δt=tfti;\Delta t^{\prime}=t^{\prime}_{f}-t^{\prime}_{i}; this is the physical interpretation of ΔxΔyΔz.\Delta x^{\prime}\Delta y^{\prime}\Delta z^{\prime}. Taking the infinitesimal value Δ𝐱d𝐱\Delta\mathbf{x}\rightarrow d\mathbf{x}, we have:

d3𝐱=γd3𝐱.d^{3}\mathbf{x}^{\prime}=\gamma d^{3}\mathbf{x}. (105)

Notice that the standard volume contraction formula is d3𝐱=γ1d3𝐱,d^{3}\mathbf{x}^{\prime}=\gamma^{-1}d^{3}\mathbf{x}, which has different interpretation with (105). The standard volume contraction is the comparison between volumes that are both measured instantaneously with respect to each times on each frames, while in (105), the measurement of d3𝐱d^{3}\mathbf{x}^{\prime} is done along a time interval Δt\Delta t^{\prime}. For an alternative derivation for this formula, one could consult (Peebles, ; Heer, ; Henry, ).

Second, let us obtain the transformation for the 3-volume element in the momentum space, which are defined as follows:

d3𝐩\displaystyle d^{3}\mathbf{p} =dpxdpydpz,\displaystyle=dp_{x}\wedge dp_{y}\wedge dp_{z}, (106)
d3𝐩\displaystyle d^{3}\mathbf{p}^{\prime} =dpxdpydpz,\displaystyle=dp_{x}^{\prime}\wedge dp_{y}^{\prime}\wedge dp_{z}^{\prime},

respectively for 𝒪\mathcal{O} and 𝒪\mathcal{O}^{\prime}. The 4-momentum are constrained such that it’s norm |𝐩4|=pμpμ\left|{}^{4}\mathbf{p}\right|=p_{\mu}p^{\mu} is constant:

pμpμ=pt2+px2+py2+pz2=m2c2.p_{\mu}p^{\mu}=-p_{t}^{2}+p_{x}^{2}+p_{y}^{2}+p_{z}^{2}=-m^{2}c^{2}. (107)

Differentiating this constraint with respect to any parameter gives:

dpt=1pt(pxdpx+pydpy+pzdpz).dp_{t}=\frac{1}{p_{t}}\left(p_{x}dp_{x}+p_{y}dp_{y}+p_{z}dp_{z}\right). (108)

Let us consider the infinitesimal version of the transformation of 4-momentum; such transformations have similar forms with (98). Inserting the infinitesimal version of (98) to (106), and then using the constraint (108), gives:

d3𝐩\displaystyle d^{3}\mathbf{p}^{\prime} =γ(1βpxpt)dpxdpydpz.\displaystyle=\gamma\left(1-\beta\frac{p_{x}}{p_{t}}\right)dp_{x}\wedge dp_{y}\wedge dp_{z}. (109)

However, coefficient in (109) is simply ptpt\frac{p_{t}^{\prime}}{p_{t}} by (98), so one has:

d3𝐩\displaystyle d^{3}\mathbf{p}^{\prime} =ptptd3𝐩=ωωd3𝐩,\displaystyle=\frac{p_{t}^{\prime}}{p_{t}}d^{3}\mathbf{p}=\frac{\omega^{\prime}}{\omega}d^{3}\mathbf{p}, (110)

by the photon momentum (97) and the Doppler effect (102). Now, using (105) and (110), we can obtain the transformation of the phase-space volume element between 2 inertial observers as follows:

d3𝐱d3𝐩\displaystyle d^{3}\mathbf{x}^{\prime}\wedge d^{3}\mathbf{p}^{\prime} =γωωd3𝐱d3𝐩,\displaystyle=\gamma\frac{\omega^{\prime}}{\omega}d^{3}\mathbf{x}\wedge d^{3}\mathbf{p},
=γ2(1βcosθ)d3𝐱d3𝐩.\displaystyle=\gamma^{2}\left(1-\beta\cos\theta\right)d^{3}\mathbf{x}\wedge d^{3}\mathbf{p}.

B.5 The Relativistic Distribution Function

The trajectory of the photon in spacetime is described by a curve in the Minkowski space. At an instantaneous (constant) time tt, one could define a hypersurface σt\sigma_{t} and extend the hypersurface to the phase-space by attaching the momentum space TpσtT_{p}^{*}\sigma_{t} (the cotangent bundle of σt.\sigma_{t}.) to σt.\sigma_{t}. The worldine that cross the phase-space will mark a point on the phase-space σt×Tpσt\sigma_{t}\times T_{p}^{*}\sigma_{t}, describing the state of the photon at time tt. Moreover, one can construct the phase-space volume element in the phase-space σt×Tpσt\sigma_{t}\times T_{p}^{*}\sigma_{t}, namely d3𝐱d3𝐩d^{3}\mathbf{x}\wedge d^{3}\mathbf{p}. The relativistic distribution function f(𝐱,𝐩)f\left(\mathbf{x},\mathbf{p}\right) is defined as the number of wordlines that cross the phase-space, i.e, the states, per volume element (Liboff, ):

f(𝐱,𝐩)=dNd3𝐱d3𝐩.f\left(\mathbf{x},\mathbf{p}\right)=\frac{dN}{d^{3}\mathbf{x}d^{3}\mathbf{p}}.

Now, if we have another coordinate patch for the Minkowski space, namely (t,x,y,z)\left(t^{\prime},x^{\prime},y^{\prime},z^{\prime}\right), it will perceive different time-constant hypersurface, different phase-space, and different phase-space volume element, namely d3𝐱d3𝐩d^{3}\mathbf{x}^{\prime}\wedge d^{3}\mathbf{p}^{\prime}. Since the phase-space volume element changes, the number of wordlines that cross the volume element will also change accordingly. However, the number of worldlines (and the worldlines themselves) in \mathcal{M} do not change by coordinate transformation. Hence, it is reasonable to assume that the worldline density f(𝐱,𝐩)f\left(\mathbf{x},\mathbf{p}\right) is invariant under coordinate (Lorentz) transformation, namely (Liboff, ):

dNd3𝐱d3𝐩=dNd3𝐱d3𝐩.\frac{dN}{d^{3}\mathbf{x}d^{3}\mathbf{p}}=\frac{dN^{\prime}}{d^{3}\mathbf{x}^{\prime}d^{3}\mathbf{p}^{\prime}}.

With this reasonable assumption, one could have the number of states/modes transformation between two inertial frames 𝒪\mathcal{O} and 𝒪\mathcal{O}^{\prime}, which is:

dNdN=d3𝐱d3𝐩d3𝐱d3𝐩=γ2(1βcosθ).\frac{dN^{\prime}}{dN}=\frac{d^{3}\mathbf{x}^{\prime}d^{3}\mathbf{p}^{\prime}}{d^{3}\mathbf{x}d^{3}\mathbf{p}}=\gamma^{2}\left(1-\beta\cos\theta\right). (111)

Finally, with (101), (103), (105), and (111), we have all the ingredients to show that the Planck’s law are invariant under Lorentz transformation.

B.6 Black-body Radiation in Moving Frame

The objective in this section is to know if the Planckian distribution is invariant under Lorentz transformation. To avoid the problem of temperature in moving bodies, we use the Planck’s Law in the form of equation (77), that is, the one with the term containing ωmax\omega{}_{\mathrm{max}}^{\prime}, instead of the original form (72) with the term containing κBT\kappa_{B}T. The derivation in this section is based on the derivation in (Henry, ), with some slight modifications. Inserting the relativistic Doppler effect (101), the transformation of solid angle (103), the volume contraction (105), and the transformation of the number of modes distribution (111) to (77), we obtain:

dNγ2(1βcosθ)=(ω)24π3c31γ2(1βcosθ)21exp(wωγ(1βcosθ)γ(1βcosθ)ωmax)1dωγ(1βcosθ)(ωω)2dΩ1γd3𝐱,\frac{dN^{\prime}}{\gamma^{2}\left(1-\beta\cos\theta\right)}=\frac{\left(\omega^{\prime}\right)^{2}}{4\pi^{3}c^{3}}\frac{1}{\gamma^{2}\left(1-\beta\cos\theta\right)^{2}}\frac{1}{\exp\left(w\frac{\omega^{\prime}}{\gamma\left(1-\beta\cos\theta\right)}\frac{\gamma\left(1-\beta\cos\theta\right)}{\omega{}_{\mathrm{max}}^{\prime}}\right)-1}\frac{d\omega^{\prime}}{\gamma\left(1-\beta\cos\theta\right)}\left(\frac{\omega^{\prime}}{\omega}\right)^{2}d\Omega^{\prime}\frac{1}{\gamma}d^{3}\mathbf{x}^{\prime},

that can be simplified as:

dN=(ω)24π3c31exp(wωωmax)1dωdΩdV,dN^{\prime}=\frac{\left(\omega^{\prime}\right)^{2}}{4\pi^{3}c^{3}}\frac{1}{\exp\left(w\frac{\omega^{\prime}}{\omega^{\prime}_{\mathrm{max}}}\right)-1}d\omega^{\prime}d\Omega^{\prime}dV^{\prime}, (112)

which is exactly the form of Planck’s law (77). All the variables inside the equation transform according to the Lorentz transformation, but the relation between these variables is invariant, hence, the observer in a moving frame will still observe the Planckian spectrum. However, there is a subtlety in the relation (112). ωmax\omega{}_{\mathrm{max}}^{\prime}, the ’maximal’ frequency, also transformed as the angular frequency via the relativistic Doppler effect (101). However , ωmax\omega{}_{\mathrm{max}}^{\prime} is the ’maximal’ frequency in the distribution (77) of a black-body in the rest frame, observed by a moving observer. It is not necessarily the ’maximal’ frequency obtained from the distribution (112) of a black-body observed by the moving observer. We will see that they are equivalent as follows. Let us multiply the LHS and RHS of (112) with ω\hbar\omega^{\prime} to obtain:

dE=(ω)34π3c31exp(wωωmax)1ρE[ω]dωdΩdV.dE^{\prime}=\underset{\rho^{\prime}_{E}\left[\omega^{\prime}\right]}{\underbrace{\frac{\hbar\left(\omega^{\prime}\right)^{3}}{4\pi^{3}c^{3}}\frac{1}{\exp\left(w\frac{\omega^{\prime}}{\omega{}_{\mathrm{max}}^{\prime}}\right)-1}}}d\omega^{\prime}d\Omega^{\prime}dV^{\prime}. (113)

EE^{\prime} is the energy of modes with frequency ω\omega^{\prime} observed in the moving frame 𝒪\mathcal{O}^{\prime}. Let us find the value of ω\omega^{\prime} that maximize EE^{\prime} from the equation dρE[ω]dω=0\frac{d\rho^{\prime}_{E}\left[\omega^{\prime}\right]}{d\omega^{\prime}}=0:

(ω)24π3c31exp(wωωmax)1(3wωωmaxexp(wωωmax)exp(wωωmax)1)=0.\frac{\hbar\left(\omega^{\prime}\right)^{2}}{4\pi^{3}c^{3}}\frac{1}{\exp\left(w\frac{\omega^{\prime}}{\omega{}_{\mathrm{max}}^{\prime}}\right)-1}\left(3-w\frac{\omega^{\prime}}{\omega{}_{\mathrm{max}}^{\prime}}\frac{\exp\left(w\frac{\omega^{\prime}}{\omega{}_{\mathrm{max}}^{\prime}}\right)}{\exp\left(w\frac{\omega^{\prime}}{\omega{}_{\mathrm{max}}^{\prime}}\right)-1}\right)=0. (114)

Defining w=wωωmaxw^{\prime}=w\frac{\omega^{\prime}}{\omega{}_{\mathrm{max}}^{\prime}} as in the previous subsection, the equation from (114) that we need to solve is then:

w=3(1ew),w^{\prime}=3\left(1-e^{-w^{\prime}}\right), (115)

which is exactly similar to (74). Hence it posses a same solution:

w=3+W0(3e3),w^{\prime}=3+W_{0}\left(-3e^{-3}\right),

namely, the Wien coefficient is invariant under Lorentz transformation: w=ww^{\prime}=w. Then, if ω=ωmax\omega^{\prime}=\omega^{\prime}_{\mathrm{max}} is the solution to (114), we have:

ww=ωmaxωmax=1,\frac{w^{\prime}}{w}=\frac{\omega^{\prime}_{\mathrm{max}}}{\omega{}_{\mathrm{max}}^{\prime}}=1,

or namely ωmax=ωmax\omega^{\prime}_{\mathrm{max}}=\omega{}_{\mathrm{max}}^{\prime}: the Lorentz transformation sends ’maximal’ frequency in one inertial frame, to a ’maximal’ frequency in another inertial frame. (113) becomes:

dN=(ω)24π3c31exp(wωωmax)1dωdΩdV.dN^{\prime}=\frac{\left(\omega^{\prime}\right)^{2}}{4\pi^{3}c^{3}}\frac{1}{\exp\left(w\frac{\omega^{\prime}}{\omega{}_{\mathrm{max}}^{\prime}}\right)-1}d\omega^{\prime}d\Omega^{\prime}dV^{\prime}. (116)

With this, we can state that the Planck’s radiation law in the form of (77) is invariant under Lorentz transformation. This result is similar to (Henry, ).

References

  • (1) C. Farias, V. A. Pinto, P. S. Moya. What is the temperature of a moving body?. Sci. Rep. 7: 17657. 2017.
  • (2) K. Derakhshani. Black Body Radiation in Moving Frames. 2019. https://hal.science/hal-02299780.
  • (3) M. Planck. Zur Dynamik bewegter Systeme. Sitzungsber. Preuss. Akad. Wiss. pp. 867–904. 1907.
  • (4) A. Einstein. Zur Elektrodynamik bewegter Ko¨\ddot{o}rper. Ann. Phys. 17: 891. 1905.
  • (5) A. Einstein. Uber das Relativitatsprinzip und die aus demselben gezogene Folgerungen. Jahrbuchder Radioaktivita¨\ddot{\mathrm{a}}t 4: 411. 1907.
  • (6) H. Ott. Lorentz-Transformation der Wa¨\ddot{a}rme und der Temperatur. Zeit. Phys. 175: 70. 1963.
  • (7) H. Arzelie`\mathrm{\grave{e}}s. Transformation relativiste de la temperature et de quelque autre grandeurs thermodynamique. Nuov. Cim. 35 (3): 792-804. 1965.
  • (8) P. T. Landsberg. Does a Moving Body Appear Cool?. Nature 212: 571. 1966.
  • (9) P.T. Landsberg. Special relativistic thermodynamics. Proc. Phys. Soc. 89: 1007. 1966.
  • (10) T. K. Nakamura. Thermodynamics of extended bodies in special relativity. Phys. Lett. A 352: 175. 2006.
  • (11) T. K. Nakamura. Three Views of a Secret in Relativistic Thermodynamics. Prog. Theor. Phys. 128: 463-475. 2012.
  • (12) N. G. van Kampen. Relativistic Thermodynamics of Moving Systems. Phys. Rev. 173 (295). 1968.
  • (13) W. Israel. Transient relativistic thermodynamics and kinetic theory. Annal. Phys. 100 (310). 1976.
  • (14) R. C. Tolman. On the Extension of Thermodynamics to General Relativity. Proc. Nat. Acad. Sci. 14: 268-272. 1928.
  • (15) R. C. Tolman. Effect of Inhomogeneity on Cosmological Models. Proc. Nat. Acad. Sci. 20: 169-176. 1928.
  • (16) A. Einstein. Relativity: The Special and General Theory. Henry Holt and Co. New York. 1920.
  • (17) S. M. Carroll. Spacetime and Geometry: An Introduction to General Relativity. Cambridge University Press. 2019.
  • (18) F. W. Sears, G. L. Salinger. Thermodynamics, Kinetic Theory, and Statistical Thermodynamics. 3rd Edition. Addison Wesley. 1975.
  • (19) L. Boltzmann. On the Relationship between the Second Fundamental Theorem of the Mechanical Theory of Heat and Probability Calculations Regarding the Conditions for Thermal Equilibrium. Wiss. Abhandlungen II 42: 64-223. 1909.
  • (20) R. L. Liboff. Kinetic Theory Classical, Quantum, and Relativistic Descriptions. Springer. 2003.
  • (21) W. G. Unruh. Notes on black-hole evaporation. Phys. Rev. D 14 (4): 870–892. 1976.
  • (22) S. A. Fulling. Nonuniqueness of Canonical Field Quantization in Riemannian Space-Time. Phys. Rev. D 7 (10): 2850–2862. 1973.
  • (23) P. C. W. Davies. Scalar production in Schwarzschild and Rindler metrics. J. Phys. A 8 (4): 609–616. 1975.
  • (24) P. M. Alsing, P. W. Milonni. Simplified derivation of the Hawking-Unruh temperature for an accelerated observer in vacuum. Am. J. Phys. 72: 1524-1529. 2004.
  • (25) S. W. Hawking. Black hole explosions?. Nature 248 (5443): 30–31. 1974.
  • (26) S. W. Hawking. Particle creation by black holes. Commun. Math. Phys. 43 (3): 199–220. 1975.
  • (27) S. Deffner, S. Campbell. Quantum Thermodynamics: An introduction to the thermodynamics of quantum information. Morgan & Claypool Publishers. 2019.
  • (28) F. Binder, L. A. Correa, C. Gogolin, J. Anders, G. Adesso. Thermodynamics in the Quantum Regime. Fundamental Theories of Physics. Springer. 2018.
  • (29) C. Rovelli. Statistical mechanics of gravity and the thermodynamical origin of time. Class. Quant. Grav. 10: 1549–1566. 1993.
  • (30) A. Connes, C. Rovelli. Von Neumann algebra automorphisms and time-thermodynamics relation in generally covariant quantum theories. Class. Quant. Grav. 11 (12): 2899–2917. 1994.
  • (31) M. Planck. On the Theory of the Energy Distribution Law of the Normal Spectrum. Verhandl. Dtsch. Phys. Ges. 2 (237). 1900.
  • (32) N. G. van Kampen. Relativistic thermodynamics. J. Phys. Soc. Japan 26: 316–321. 1969.
  • (33) J. J. Maresˇ\mathrm{\check{s}}, P. Hubik, J. Sˇ\mathrm{\check{S}}esta´\mathrm{\acute{a}}k, V. Sˇ\mathrm{\check{S}}picˇ\mathrm{\check{c}}ka, J. Krisˇ\mathrm{\check{s}},tofik, J. Sta´\mathrm{\acute{a}}vek. Relativistic transformation of temperature and Mosengeil–Ott’s antinomy. Physica E 42 (3): 484-487. 2010.
  • (34) S. S. Costa, G. E. A. Matsas. Temperature and relativity. Phys. Lett. A 209 (155). 1995.
  • (35) P. T. Landsberg, G. E. A. Matsas. Laying the ghost of the relativistic temperature transformation. Phys. Lett. A 223 (6): 401-403. 1996.
  • (36) P. T. Landsberg, G. E. A. Matsas. The impossibility of a universal relativistic temperature transformation. Physica A 340 (1–3): 92-94. 2004.
  • (37) R. Bracewell, E. Conklin. An Observer moving in the 3o3^{o} K Radiation Field. Nature 219: 1343–1344. 1968.
  • (38) W. Wien. U¨\ddot{U}ber die Differentialgleichungen der Elektrodynamik fu¨\ddot{u}r bewegte Ko¨\ddot{o}rper. I. Ann. der Phys. 318 (4): 641–662. 1904.
  • (39) W. Wien. U¨\ddot{U}ber die Differentialgleichungen der Elektrodynamik fu¨\ddot{u}r bewegte Ko¨\ddot{o}rper. II. Ann. der Phys. 318 (4): 663–668. 1904.
  • (40) J. Stefan. U¨\ddot{U}ber die Beziehung zwischen der Wa¨\ddot{a}rmestrahlung und der Temperatur. Sitzungsber. Preuss. Akad. Wiss. 79: 391–428. 1879.
  • (41) L. Boltzmann. Ableitung des Stefan’schen Gesetzes, betreffend die Abha¨\ddot{a}ngigkeit der Wa¨\ddot{a}rmestrahlung von der Temperatur aus der electromagnetischen Lichttheorie. Ann. der Phys. Chem. 258 (6): 291–294. 1884.
  • (42) P. J. E. Peebles, D. T. Wilkinson. Comment on the Anisotropy of the Primeval Fireball. Phys. Rev. 174 (2168). 1968.
  • (43) C. V. Heer, R. H. Kohl. Theory for the Measurement of the Earth’s Velocity through the 303^{0}K Cosmic Radiation. Phys. Rev. 174 (1611). 1968.
  • (44) G. R. Henry, R. B. Feduniak, J. E. Silver, M. A. Peterson. Distribution of Blackbody Cavity Radiation in a Moving Frame of Reference. Phys. Rev. 176 (1451). 1968.
  • (45) T. H. Boyer. Derivation of the Blackbody Radiation Spectrum without Quantum Assumptions. Phys. Rev. 182 (1374). 1969
  • (46) J. S. Lee, G. B. Cleaver. The relativistic blackbody spectrum in inertial and non-inertial reference frames. New Astronomy 52: 20-28. 2017.
  • (47) M. H. Johnson, E. Teller. Intensity changes in the Doppler effect. Proc. Natl. Acad. Sci. 79 (4) :1340. 1982.
  • (48) R. P. Feynman, R. B. Leighton, M. Sands. Relativistic Effects in Radiation. The Feynman Lectures on Physics: Volume 1. Addison-Wesley. pp. 34–7. 1977.
  • (49) L. Euler. De serie Lambertina Plurimisque eius insignibus proprietatibus. Acta Acad. Scient. Petropol. 2: 29-51. (1783).
  • (50) W. Rindler. Relativity: Special, General, and Cosmological. 2nd Ed. Oxford University Press. 2006.
  • (51) K. S. Thorne, R. D. Blandford. Modern Classical Physics: Optics, Fluids, Plasmas, Elasticity, Relativity, and Statistical Physics. Princeton University Press, 2017.
  • (52) F. Debbasch, J. P. Rivet, W. A. van Leeuwen. Invariance of the relativistic one-particle distribution function. Physica A 301 (1–4): 181-195. 2001.
  • (53) M. H. Cohen, M. L. Lister, D. C. Homan, M. Kadler, K. I. Kellermann, Y. Y. Kovalev, R. C. Vermeulen. Relativistic Beaming and the Intrinsic Properties of Extragalactic Radio Jets. The Astro. Jour. 658 (1). 2007.
  • (54) A. S. Eddington. Einstein Shift and Doppler Shift. Nature. 117 (2933): 86. (1926).
  • (55) L.B. Okun (ITEP, Moscow), K.G. Selivanov (ITEP, Moscow), V.L. Telegdi (EP Division, CERN). On the Interpretation of the Redshift in a Static Gravitational Field. Am.J.Phys.68:115,2000.
  • (56) V. A. Golyk, M. Kru¨\mathrm{\ddot{u}}ger, M. Kardar. Heat radiation from long cylindrical objects. Phys. Rev. E 85 (4): 046603. 2012.
  • (57) M. Planck. U¨\ddot{U}ber die Begru¨\ddot{u}ndung des Gesetzes der schwarzen Strahlung. Ann. Phys. 342: 642. 1912.
  • (58) S. Weinberg. The Quantum Theory of Fields, Volume 1: Foundations. Cambridge University Press. 2005 .
  • (59) R. Jaffe. Casimir effect and the quantum vacuum. Phys. Rev. D. 72 (2): 021301. 2005.
  • (60) R. Aldrovandi, J. Gariel. On the riddle of the moving thermometers. Phys. Lett. A 170 (1): 5-10. 1992.