Light from Schwarzschild black holes in de Sitter expanding universe

A new method is applied for deriving simultaneously the redshift and shadow of a Schwarzschild black hole moving freely in the de Sitter expanding universe as recorded by a remote co-moving observer. This method is mainly algebraic, focusing on the transformation of the conserved quantities under the de Sitter isometry relating the black hole co-moving frame to observer’s one. Hereby one extracts the general expressions of the redshifts and shadows of the black holes having peculiar velocities but their expressions are too extended to be written down here. Therefore, only some particular cases and intuitive expansions are presented while the complete results are given in an algebraic code (Cotăescu in Maple code BH01, https://physics.uvt.ro/~cota/CCFT/codes, 2020).


Introduction
The light emitted by cosmic objects is one of the principal sources of empirical data in astrophysics. An important accessible observable is the redshift which encapsulates information about the cosmic expansion and possible peculiar velocity of the observed object. For separating these two contributions one combined so far the Lemaître rule [2,3] of Hubble's law [4], governing the cosmological effect, [5][6][7] with the usual theory of the Doppler effect of special relativity [8] even though there is evidence that our universe is expanding.
Recently, we proposed an improvement of this approach replacing special relativity with our de Sitter relativity [9,10] where local charts of the same type, playing the role of inertial frames, are related among themselves through de Sitter isometries. In the case of the longitudinal Doppler effect, when a point-like source is moving along the axis observersource, we obtained a redshift formula having a new term a e-mail: i.cotaescu@e-uvt.ro (corresponding author) combining the cosmological and kinetic contributions in a non-trivial manner [11].
The next step is to extend this method to the black holes which, in general, are no point-like sources. The light emitted by a black hole comes from an apparent source situated on a sphere, surrounding the black hole, which is observed as the black hole shadow. When this is not negligible, as in the case of the object M87 [12,13], and the black hole may have a peculiar velocity, the Doppler effect is no longer longitudinal such that the transverse contributions due to the black hole shadow must be evaluated. This can be done only by studying simultaneously the Doppler effect and black hole shadow in the same theoretical framework.
As in this paper we would like to study the relation between the redshifts and shadows of the Schwarzschild black holes having peculiar velocities in de Sitter expanding universe we must abandon the geometric method adopting the algebraic approach of Ref. [11] which is suitable for deriving the redshift. This is based on the de Sitter relativity where we may relate the moving black hole proper frames to those of remote co-moving observers freely falling in the de Sitter expanding universe.
We start supposing that our expanding universe is satisfactory described by the expanding portion of a (1 + 3)dimensional de Sitter manifold. As the actual observations show with reasonable accuracy that this universe is spatially flat, we consider only local charts with Painlevé coordinates [44] since these have flat space sections. These local charts are de Sitter co-moving frames [45] where the coordinates are the cosmic time and Cartesian or spherical space coordinates. These frames may carry observers related among themselves through the de Sitter isometries which transform simultaneously the coordinates and the conserved quantities [9,10].
The Schwarzschild-de Sitter black holes are usually considered in proper frames with static coordinates and Kottler metric [46]. However, here we prefer the corresponding comoving frames with Painlevé coordinates whose metrics have the same asymptotic behavior as the metric of the observer co-moving frame. Then for a remote observer the black hole co-moving frame appears as an empty de Sitter one which can be related to observer's co-moving frame through a de Sitter isometry, in accordance with the relative motion of the black hole with respect to observer.
We assume that at the initial moment when the black hole emits the photon this is translated and has a relative velocity with respect to the remote observer known as peculiar velocity. For avoiding extremely complicated calculations we assume that this velocity is longitudinal, in the black holeobserver direction. Thus we start with precise initial conditions determining the suitable isometry relating the black hole and observer proper frames. This will give the conserved quantities measured in the observer's frame we need for extracting physical results without resorting to geodesics or other geometric objects. For this reason we say that our method is algebraic, observing that there are some advantages among them the principal one is of a coherent framework offered by the de Sitter relativity which prevent us of using supplemental hypotheses.
In this approach we derive the related redshift and shadow of a Schwarzschild black hole freely moving in the de Sitter expanding universe. Our principal new results are a closed formula of the black hole shadow depending on its peculiar velocity and the corrections to our new redshift formula derived in Ref. [11] due to the dimension of the black hole shadow. These corrections are too complicated to be written down here but can be derived with the help of an algebraic code on computer [1].
We start in the second section with a brief review of the metrics of the black hole and observer co-moving frames with Cartesian or spherical coordinates, revisiting the equation giving the geodesic shapes in the black hole co-moving frames. In the next section we present the solutions of this equation representing the null geodesics around the black hole. These are the circular geodesics on the photon sphere and the Darwin [47,48] spiral geodesics which determine the black hole shadow. The fourth section is devoted to the de Sitter isometries relating the conserved quantities measured in different co-moving frames. In the next section we obtain our new results assuming that a remote observer sees that the light is emitted by an apparent source on a null de Sitter geodesic whose conserved quantities can be determined. Furthermore, by using an isometry formed by a translation followed by a Lorentzian isometry we obtain the conserved quantities in the observer's co-moving frame from which we extract the observed redshift and shadow of the moving black hole. More specific, the redshift results from the observed energy while the angular radius of the black hole shadow is derived by using the components of the photon momentum in the observer's origin, where the photon angular momentum must vanish. These results are elementary but with a large number of terms that cannot be written here in the general case of a moving black hole. Consequently, we restrict ourselves to presenting here only their series expansions with respect to a common small parameter, the particular case when the relative velocity vanishes and the flat limit. As mentioned, the complete results which cannot be written here are given in an algebraic code on computer [1]. In order to convince oneself that the flat limit is correct we derive in Appendix B the results that can be obtained by applying our algebraic method to a Schwarzschild black hole in Minkowski flat space-time. Finally we present some concluding remarks.
As our approach may be applied even in quantum theory we introduce a special notation denoting by ω H = 3 c the de Sitter Hubble constant (frequency) since H is reserved for the Hamiltonian operator [49]. Moreover, the Hubble time t H = 1 ω H and the Hubble length l H = c ω H will have the same form in the natural Planck units with c =h = G = 1 we use here.
Here we focus on a Schwarzschild black hole of mass M embedded in the de Sitter expanding universe for which the metric (2) of its static frame has the Kottler [46] (or Schwarzschild-de Sitter) form with where, as mentioned before, ω H is the de Sitter Hubble constant in our notation. The corresponding frames with Painlevé coordinates have the asymptotic behavior of the de Sitter co- For this reason we say that the black hole frames with Painlevé coordinates, denoted by {t, x} BH and {t, r, θ, φ} BH , are the comoving frames of the Schwarzschild black hole in de Sitter expanding universe.
The frames of the remote observers, {t, x} and {t, r, θ, φ}, located in the asymptotic zone, are genuine de Sitter comoving frames where the astronomical observations are performed and recorded. The observers stay at rest in the origins of their own frames evolving along the unique time-like Killing vector field of the de Sitter geometry which is not time-like everywhere but has this property just in the null cone where the observations are allowed [49].
Here we use simultaneously Cartesian and spherical coordinates since the Cartesian coordinates are suitable for studying the conserved quantities and the transformation rules under isometries while the spherical coordinates help one to integrate the geodesic equations. For example, the spherical symmetry is obvious in Cartesian coordinates where the metrics (1) and (3) are invariant under the global rotations x i → R i j x j such that we can use the vector notation. On the other hand, only in spherical coordinates one can integrate the geodesics equations in the black hole co-moving frame we revisit briefly in the next.
In the frame {t, r, θ, φ} BH with the line element (4) the conserved quantities along geodesics are the energy E and angular momentum L. These give rise to the prime integrals of a geodesic of a particle of mass m moving in the equatorial plane of the black hole (with fixed θ = π 2 ) as where 'dot' denotes the derivatives with respect to the affine parameter λ which satisfies ds = m dλ. The third prime integral comes from the line element in the equatorial plane, which reads as it results from Eq. (4). Hereby one may derive the function After a little calculation, combining the above prime integrals, one obtains the well-known equation giving the geodesic shapes but which is not enough for finding the time behavior of the functions φ(t) and for which one must apply special methods [50]. Note that Eq. (11) derived in the co-moving frame is the same as that of the static frame since this equation is static giving only the shape of trajectory in the same space coordinates. In fact, the time evolution on geodesics is quite different in the static and co-moving frames.

Light around black holes
The problem of the gravitational lensing which has a long history [51] was studied in general relativity first by Einstein and Eddington [52] but was solved by Darwin [47,48] which derived the null geodesics around the photon sphere of a Schwarzschild black hole in the flat Minkowski space-time. Applying the same commonly used method [53][54][55][56][57][58] we may inspect briefly the null geodesics in the co-moving frame {t, r, θ, φ} B H of the Schwarzschild-de Sitter system.
The shapes of the photon geodesics are given by the functions r (φ) which satisfy Eq. (11) with m = 0 that now reads This equation has two types of solutions, namely circular geodesics on the photon sphere and associated spiral geodesics [47]. The circular geodesics satisfy simultaneously the conditions Fig. 1 The functions r (±) (φ) of the spiral photon geodesics closest to the photon sphere of radius 3M giving the radius of the photon sphere r ph = 3M and the mandatory condition derived in Ref. [15]. Thus the photons with circular geodesics are trapped on the photon sphere without escaping outside. Furthermore, by substituting the condition (14) in Eq. (12) we obtain the equation which is independent on the Hubble de Sitter constant ω H . Apart from the circular geodesics, this equation allows the solutions known as the spiral geodesics [47]. These are determined up to a rotation, φ → φ − φ 0 , fixing the origin of this angular coordinate. For example, if we translate the arguments of the functions r (±) as φ → φ ± = φ ± ln 6M then we recover the elegant Darwin form [47] 1 of the spiral geodesics. Thus we may conclude that the presence of de Sitter gravity is encapsulated only in Eq. (14), while the photon sphere and the shapes of the spiral geodesics remain the same as in Minkowski flat space-time (when ω H = 0). The spiral geodesics are symmetric, In Fig. 1 we see that the functions , since between the vertical asymptotes their values are negative having no physical meaning. It is interesting that this opaque window is independent on the black hole mass, On the physical domain these trajectories remain outside the photon sphere, r (±) (φ) > 3M, but approaching to this for large |φ| since This gives us the image of the spiral geodesics rolled out around the photon sphere (as in Fig. 2) escaping outside only when φ is approaching to the values (18) where the functions r (±) can take larger values near singularities. The geodesics r (±) (φ) are the closest trajectories to the photon sphere of the first photons that can be observed at the limit of the black hole shadow. Therefore, for studying this shadow and the associated redshift we have to consider only these photons.

de Sitter isometries
The de Sitter co-moving frames play the role of inertial frames being related among themselves through de Sitter isometries as in our de Sitter relativity [9,10]. Moreover, the black hole frames have the asymptotic de Sitter symmetry which governs the relative motion of the black hole with respect to remote observers such that we may use these isometries for relating the observer co-moving frames to the black hole one.
The de Sitter isometries can be studied easily since this manifold is a hyperboloid of radius 1/ω H embedded in the five-dimensional flat space-time (M 5 , η 5 ) of coordinates z A (labelled by the indices A, B, . . . = 0, 1, 2, 3, 4) and metric η 5 = diag(1, −1, −1, −1, −1). The local charts can be introduced giving the set of functions z A (x) which solve the hyperboloid equation, giving the line element The functions that introduce our Painlevé coordinates are The de Sitter isometry group is just the stable group SO(1, 4) of the embedding manifold (M 5 , η 5 ) that leave invariant its metric and implicitly Eq. (21). Therefore, given a system of coordinates defined by the functions z = z(x), each transformation g ∈ SO(1, 4) gives rise to the isometry x → x = φ g (x ) derived from the system The local charts related through these isometries play the same role as the inertial frames of special relativity. The classical conserved quantities under de Sitter isometries are given by the Killing vectors k (AB) of the de Sitter manifold [49] that are related to those of (M 5 , η 5 ) as allowing us to derive the covariant components of the Killing vectors in an arbitrary chart {x} of the de Sitter space-time as where z A = η AB z B . The conserved quantities along the timelike geodesic of a particle of mass m have the general form K (AB) (x, P) = ω H k (AB) μẋ μ . The conserved quantities with physical meaning are the energy E, momentum P, angular momentum L and a specific vector Q that we call the adjoint momentum [49]. A geodesic in the co-moving frame {t, x} [59], depends only on the momentum P (P = |P|) and the initial condition x(t 0 ) = x 0 fixed at the time t 0 . The conserved quantities in an arbitrary point (t, x(t)) of this geodesic read [9,59] satisfying the obvious identity corresponding to the first Casimir invariant of the SO(1, 4) algebra [49]. In the flat limit, when ω H → 0, we have Q → P such that this identity becomes just the usual mass-shell condition, E 2 − P 2 = m 2 , of special relativity. The conserved quantities E, P and the new ones, form a skew-symmetric tensor on M 5 , whose components transform under the isometries x → x = φ g (x ) defined by Eq. (24) as where g = η 5 g η 5 [9]. Summarizing, we can say that the de Sitter isometries are generated globally by the SO (1, 4) transformations which determine the transformations of the coordinates and conserved quantities. We have thus a specific relativity on the de Sitter space-time allowing us to study different relativistic effects in the presence of the de Sitter gravity. In what follows we use the Lorentzian isometries defined in Ref. [9] and the translations presented in the Appendix A.

Observing light from black holes
Le us consider now a mobile black hole in its proper comoving frame {t , x } B H with the origin in O B H and a fixed remote observer in his own co-moving frame {t, x} having the origin in O. We consider that the space Cartesian axes of these frames remain parallel with the basis of unit vectors (e 1 , e 2 , e 3 ) such that the geodesic of the emitted photon is in the plane (e 1 , e 2 ). In this geometry we assume that the photon is emitted at the initial moment t = t = 0 when the origin O B H is translated with d and has the relative velocity V = e 1 V with respect to O. Note that the velocity V = P M is conserved depending on the conserved momentum P of the black hole geodesic observed by O.

Related conserved quantities
A remote observer sees the photon of momentum k = n k k, energy E ph = |k| and angular momentum (14), as emitted from an apparent source S of position vector n S r S on the sphere of radius which is just the apparent radius of the black hole shadow. Therefore, r S is the radius of the sphere hosting different photon sources S that can be observed nearest to the black hole shadow (as in Fig. 3). When ω H → 0 this becomes just the shadow radius 3 √ 3M derived in special relativity [47]. The apparent trajectory of the emitted photon is a de Sitter null geodesic of momentum k that, according to Eq. (27), reads complying with the initial condition x ph (0) = n S r S . This geodesic depends on the orthogonal unit vectors, which can be represented as n k = −e 1 cos α − e 2 sin α, (37) n S = −e 1 sin α + e 2 cos α, (38) where the angle α, giving the apparent direction of the photon, will depend on observer's position. We have thus the opportunity of defining the de Sitter conserved quantities on this geodesic as in an apparent de Sitter empty comoving frame {t , x } associated to {t , x } B H . Here one vector is missing, namely the adjoint momentum Q that can be derived simply at the time t = 0 according to Eq. (30). We complete thus the set of conserved quantities, which satisfy the condition (31). Now we can deduce how these conserved quantities are measured by the fixed observer O since the observer frame {t, x} and the apparent black hole one, {t , x }, are related through an isometry, x = φ g (x ), of the de Sitter relativity [9]. According to our hypotheses, this is generated by the SO(1,4) transformation, formed by a translation (A.1) of parameter a = e 1 d = (d, 0, 0), having the form [9,60] g(a) = followed by the Lorentz boost of the particular Lorentzian isometry we need here [9]. Applying then the transformation (34) with g given by Eq. (43) we obtain the conserved quantities observed by O. This calculation is elementary but complicated, involving many terms that can be manipulated only by using suitable algebraic codes on computer. For presenting the final result it is convenient to introduce the notation which allows us to write down the conserved quantities in the observer's frame as while L 1 = L 2 = P 3 = Q 3 = 0. As expected, these quantities satisfy the invariant identity (31).
In other respects, we observe that all the vector components we derived above can take any real values in contrast with the energy which must remain positive definite. This condition is fulfilled only if This is in fact the mandatory condition for observing the photon in O at finite time. When the relative velocity V exceeds this limit then the photon cannot arrive in O at finite time because of the background expansion. Thus V lim defines a new velocity horizon restricting the velocities such that for very far sources with α = 0 and δ = ω H d ∼ 1 this limit vanishes.

Shadow and redshift
The angle α, which depends on the relative position between black hole and observer, can be found simply imposing the condition L 3 = 0 when the photon is passing through the point O. Solving this equation for L 3 given by Eq. (52) we find Substituting then this angle in Eqs. (47)-(53) we obtain all the conserved quantities measured by O and the maximal velocity V lim . Hereby we can extract the quantities of general interest, namely sin α, measured in the black hole proper frame, the angular radius sin α obs = |P 2 | P , P = P 1 2 + P 2 2 (55) of the shadow, measured in the observer's frame, and the redshift z defined as Unfortunately, the exact expressions of these quantities have a huge number of terms that cannot be written here but can be manipulated on computer [1] for extracting significant particular cases or intuitive approximations. The simplest particular case is when the black hole does not have a initial relative velocity with respect to O. Then by setting V = 0 we obtain the simple formulas showing how the observations of the shadow and redshift are related each other.
In the general case of V = 0, we observe that the expansions around ξ = 0 are useful since this is the only parameter which remains very small in astronomical observations as long as the other ones have larger ranges, 0 < δ = ω H d < 1 and 0 < V < V lim . Computing these series we write down here only the lowest terms which are still comprehensible and can be interpreted. Similarly, for the velocity limit we have The principal novelty here is that sin α obs depends on the relative velocity V as in Eq. (60) just from the first order of the expansion which may be observationally accessible. In contrast, the expansion (61) has a first term independent on ξ recovering just the redshift formula we derived recently for a point-like source moving along the e 1 axis [11]. Therefore, the influence of the black hole dimensions could be observed only when the second term of the order O(ξ 2 ) could be measured with a satisfactory accuracy.

Flat limit
The limit of the de Sitter relativity when the de Sitter-Hubble constant ω H vanishes is just the usual version of special relativity. Then all the co-moving frames of the de Sitter relativity become inertial frames in Minkowski space-time without affecting the black hole geometry in its proper frame. Now a remote observer sees a photon emitted in S as having an apparent rectilinear trajectory with momentum k and energy E ph = k. Then all the measured quantities can be obtained from the de Sitter ones in the limit ω H → 0. In this limit our principal parameters becomes while the other quantities have the limits where As expected, hereby we recover the usual aberration formulas In addition, we obtain the limit velocitŷ which does not make sense since this can exceed the speed of light. Thus the limitation of velocities disappears together with the de Sitter event horizon. For interpreting these results a good choice is Δ 1 as the parameter ξ remains very small. Then Eq. (66) is just the redshift due to the Doppler effect in special relativity. Moreover, we may convince ourselves that all the limits derived above are just the results that may be obtained by applying our method in special relativity, as presented briefly in Appendix B.

Concluding remarks
We studied how a co-moving observer measures simultaneously the shadow and redshift of a Scwarzschild black hole freely falling in the de Sitter expanding universe. For this purpose we used a new algebraic method offered by our de Sitter relativity that provides us with suitable isometries transforming the conserved quantities of the emitted photon into those recorded by a remote co-moving observer. In this manner we obtained the closed formula (60) of the shadow depending on the peculiar velocity and the corrections to the new redshift formula (61) that can be calculated on computer by using the code [1].
Another advantage of our method is that this is somewhat independent on the coordinates which are involved only in imposing the initial conditions. For example, the choice of the co-moving frames with Painlevé coordinates simplifies the calculations since then we use the translation (2) which does not affect the time. In contrast, in static coordinates, defined by Eq. (4), the same translation gives the transformation (6) which affects the time such that it is more difficult to synchronise the clocks by setting common initial conditions when V = 0. However, this is not a real impediment as long as we know how the coordinates transform among themselves.
This relative independence on coordinates can be tested in the case of V = 0 applying our method to a static black hole. Our preliminary calculations indicate that the shadow of the static black hole is given by Eq. (57) just as in the case of the co-moving frames with Painlevé coordinates. This stability comes from the fact that in both cases we start with the same conserved quantities (39)-(42) transformed by the same translation (44) whose parameter d is the physical distance between black hole and observer.
Under such circumstances we may compare how the algebraic and geometric methods work in determining the black hole shadow at least in the case of V = 0. The shadow formula derived in Ref. [15] by using the geometric method can be written in our notation as sin α obs = r S d f (d ) where d is now the radial coordinate of the fixed observer. Thus we see that in our approach the factor f (d ) is missing. This means that between these two methods there are some minor differences that may come from the approximation of remote observers on which the algebraic method is based and from the fact that in the static frame the radial coordinate d does not coincide with the physical distance d. However, now it is premature to say more about this relationship before analyzing many examples in de Sitter relativity.
In other respects, we must specify that the study of the conserved quantities is not enough for understanding the entire information carried out by the light emitted by moving black holes. There are important observable quantities resulted from the coordinate transformations under isometries as, for example, the photon propagation time or the real distance between observer and black hole at the time when the photon is measured. In Ref. [11] we derived such quantities in the longitudinal case of a point-like source moving along the observer-black hole direction. Therefore, when we apply the algebraic method the coordinate transformations under isometries or other geometric tools may complete our investigation.
We hope that the algebraic method proposed here will improve the general geometric approach for getting over the difficulties in analyzing the light emitted by various cosmic objects moving in the de Sitter expanding universe.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: There are no experimental or numerical data. The code [1] is exclusively algebraic.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .
Solving again Eq. (24) for these functions and the transformation (A.1) we find the transformation rules under translations, which we present here for the first time.

Appendix B: Moving black holes in special relativity
Let us consider now a Schwarzschild black hole embedded in the flat space-time where the remote observers stay in inertial frames. We apply our algebraic method assuming that at the initial time t  (45). Performing this isometry we obtain the four-vectors x = x and p = p whose components give the energy and angular momentum observed in O as The condition L 3 = 0 of the photon passing through O gives