Shadow of a Schwarzschild black hole surrounded by a Bach-Weyl ring

We have studied the shadows of a Schwarzschild black hole surrounded by a Bach-Weyl ring through the backward ray-tracing method. The presence of Bach-Weyl ring leads to that the photon dynamical system is non-integrable and then chaos would appear in the photon motion, which affects sharply the black hole shadow. The size and shape the black hole shadow depend on the black hole parameter, the Bach-Weyl ring mass and the Weyl radius between black hole and ring. Some self-similar fractal structures also appear in the black hole shadow, which originates from the chaotic lensing. We also study the change of the image of Bach-Weyl ring with the ring mass and the Weyl radius. Finally, we analyze the invariant manifolds of Lyapunov orbits near the fixed points and discuss further the formation of the shadow of a Schwarzschild black hole with Bach-Weyl ring.


I. INTRODUCTION
The first image of the supermassive black hole in the center of the giant elliptical galaxy M87 has been obtained by using Event Horizon Telescope [1][2][3][4][5][6], which is one of the most exciting events in astrophysics and black hole physics since it provides the first direct visual evidence that there exists exactly black hole in our Universe. The information stored in this image can help us to understand further shadow of black hole, the matter accretion process and black hole jets. A black hole shadow is a two-dimensional dark region in the observer's sky, which is caused by light rays that fall into an event horizon when propagated backwards in time. Since the shape and size of shadow carry the fingerprint of the geometry around the black hole [7,8], the shadow can be regarded as a potential tool to probe black hole parameters. For example, the shadow of a Schwarzschild black hole is a perfect black disk, but it becomes an elongated silhouette for a Kerr-like black hole due to the dragging effect arising from the black hole rotation [7,8]. Moreover, the recent investigations have shown that the cusp silhouette of shadow appears for a Kerr black hole with Proca hair [9] and for a Konoplya-Zhidenko rotating non-Kerr black hole [10] if the spacetime parameters lie in a certain range.
Interestingly, there exist self-similar fractal structures in the black hole shadow for the cases where the photon motion is non-integrable, such as, a rotating black hole with scalar hair [11][12][13][14], a Majumdar-Papapetrou binary black hole system [15,16], Bonnor black diholes with magnetic dipole moment [17], and a non-Kerr rotating compact object with quadrupole mass moment [18]. These novel structure and patterns also appear in the shadows of a two dynamic black hole system with the positive cosmological constant [19]. It is found that these self-similar fractal structures and patterns in black hole shadows are determined actually by the nonplanar bound photon orbits [9] and the invariant phase space structures [20] of the photon motion in the black hole spacetimes. The recent investigations on the shadows of black holes characterizing by other parameters  indicate that these black hole parameters yield the richer structure and patterns for shadows in various theories of gravity.
In general, a real black hole in the galaxy owns an accretion disk. This extra gravitational source around the black hole also affect the propagation of photon and further modify the shape of black hole shadow.
Considering the radiation from a thin accretion disk around a black hole, Luminet obtained a simulated photograph of a Schwarzschild black hole with thin accretion disk [46], which shows that there are primary and secondary images of the thin accretion disk outside black hole shadow. The image for a Kerr black hole with Keplerian accretion disk has been simulated in Refs. [47][48][49]. With the general relativistic ray-tracing code, F. H. Vincent et al [50] has studied the images of a thin infinite accretion disk and Ion torus around compact objects. These investigations are very important to understand the effects of accretion disk on the black hole images. It is well known that the photon motion is non-integrable in the gravity system of a black hole with extra gravitational sources including accretion disk. However, it is still an open issue how the non-integrable of photon motion arising from extra gravitational sources affect the black hole shadows.
The main purpose of this paper is to answer it through studying the shadow of a Schwarzschild black hole surrounded by Bach-Weyl ring. The Bach-Weyl ring is a concentric thin ring described by the Bach-Weyl solution [51], which can be looked as a general relativistic analog of a Newtonian ring of constant density.
For a Schwarzschild black hole with Bach-Weyl ring, the presence of the ring changes the spacetime structure and affects geodesic motion of particle in the spacetime [52][53][54]. It is shown that there exist chaotic motion for particle with certain initial conditions in this background spacetime [55]. In this paper, we will study the shadow of Schwarzschild black hole with Bach-Weyl ring and then probe the effects of the photon chaotic motion on the shadow.
The paper is organized as follows. In Sec. II, we review briefly the spacetime of a Schwarzschild black hole surrounded by Bach-Weyl ring and then analyze the null geodesics equation in this spacetime. In Sec. III, we adopt the backward ray-tracing technique and present numerically the shadows for the Schwarzschild black hole surrounded by Bach-Weyl ring. In Sec.IV, we analyze the invariant phase-space structures and further explain the formation of shadow cast by a Schwarzschild black hole with a Bach-Weyl ring. Finally, we present a summary.

II. THE SPACETIME OF A SCHWARZSCHILD BLACK HOLE WITH A BACH-WEYL RING AND NULL GEODESICS
The spacetime of a vacuum static and axially symmetric spacetime, in general, can be described by the Weyl metric where ν and λ only are the functions of ρ and z. The function ν(ρ, z) satisfies the Laplace equation and can be superposed linearly, which behaves like the gravitational potential in the Newtonian theory. However, the quantity λ does not own such a property of linear superposition. For a gravity system containing two gravitational sources with their individual functions ν 1 , λ 1 and ν 2 , λ 2 , one can write the functions ν and λ for the whole system as ν = ν 1 + ν 2 and λ = λ 1 + λ 2 + λ int , respectively. The quantity λ int is the interaction term which obeys the equations [52] λ int,ρ = 2ρ(ν 1,ρ ν 2,ρ − ν 1,z ν 2,z ) (2) λ int,z = 2ρ(ν 1,ρ ν 2,z + ν 1,z ν 2,ρ ).
Here, we focus on the spacetime generated by a Schwarzschild-type black hole with a thin ring described by the Bach-Weyl solution. In the Schwarzschild coordinates, the gravity of such a spacetime can be described by the metric [51,52,55] where λ ext = λ BW + λ int . The functions ν BW and λ BW have the forms [51,52] where M is the mass of Bach-Weyl ring, b is the Weyl radius and l 1,2 = (ρ ∓ b) 2 + z 2 . K(k) and E(k) are the 1st and the 2nd kind complete elliptic integrals with the forms For a Schwarzschild black hole with a Bach-Weyl ring, the function λ int has not an analytical form, which means that the form of λ int can be obtained only by resorting to numerical computations.
With the metric (3), one can find that the Hamiltonian of a photon propagation in the spacetime of a Schwarzschild black hole surrounded by a Bach-Weyl ring can be expressed as It is obvious that this Hamiltonian is not an explicit function of time coordinate t and angular coordinate φ.
Therefore, there exist two integration constants E and L z for the null geodesics motion Making use of these two conserved quantities, we find the equations of photon motion in this spacetime with a Bach-Weyl ring can be simplified aṡ Due to the Bach-Weyl ring, one can find the r− component and θ− component do not decouple from each other, which means that the dynamical system of photon in this spacetime is non-integrable and then chaos could occur in this case. It implies that the presence of Bach-Weyl ring will bring some novel features of the black hole shadow.

III. SHADOW CASTED BY A SCHWARZSCHILD BLACK HOLE SURROUNDED BY A BACH-WEYL RING
In this section, we adopt the backward ray-tracing technique [11-18, 20, 56] to simulate numerically the shadow of a Schwarzschild black hole surrounded by a Bach-Weyl ring. Here, it is assumed that the light rays evolve from the observer backward in time and then the information carried by each ray would be respectively assigned to a pixel in a final image in the observer's sky. Along this spirit, we solve numerically the null geodesic equation (9) for each pixel in the final image with the corresponding initial condition and obtain the image of shadow in observer's sky which is composed of the pixels connected to the light rays falling down into the horizon of black hole. Since the spacetime (3) is asymptotically flat, one can expand the observer basis {et, er, eθ, eφ} as a form in the coordinate basis [11-18, 20, 56] where the transform matrix e ν µ satisfies g µν e μ α e ν β = ηαβ, and ηαβ is the usual Minkowski metric. Considering that the spacetime of a Schwarzschild black hole with a Bach-Weyl ring (3) is static and axially symmetric, it is convenient to choice a decomposition And then, the locally measured four-momentum pμ of a photon can be obtained through the projection of its four-momentum p µ onto eμ, i.e., With these quantities, one can obtain the coordinates of photon's image in observer's sky [10-18, 20, 56] x = −r pφ pr where the spatial position of observer is set to (r obs , θ obs ).
In Fig.1, we present the dependence of the black hole shadow on the mass of Bach-Weyl ring for the fixed Weyl radius. Here we assume that Bach-Weyl ring does not emit light by itself and the observer is located on the equatorial plane of the black hole (the inclination angle θ obs = π/2). With the increase of the ring mass M, we find that the shadow becomes gradually prolate along the axis of symmetry, but turns concave in the equatorial plane, which leads to that the black hole shadow becomes a "8" type shape in the case with  Moreover, the size of the black hole shadow decreases with the Weyl radius b.
In Fig.3, we plot the shadows for the observers with different inclination angle. As θ obs = 0 • , it is obvious   [17,20].
For the spacetime of a Schwarzschild black hole surrounded by a Bach-Weyl ring (3), the fixed point x 0 = (r 0 , θ 0 , 0, 0) in phase space (r, θ, p r , p θ ) is defined by the conditioṅ which leads to Linearizing the equations (14),Ẋ = JX, one can analyze the local stability of the fixed point x 0 =(r 0 , θ 0 , 0, 0). Here the vector X ≡ (x µ ,p µ ) and J is the Jacobian. The fixed points of the dynamics for the photon motion are light rings, which are the circular photon orbits in the equatorial plane [9,20]. Considering that the function λ in the metric (3) is not analytical, we here adopt the case with M = 1.5M and b = 3M as an example to analyze numerically the formation of the shadow of a Schwarzschild black hole with a Bach-Weyl ring (3) which is shown in Fig.2

(b).
Solving equation (15) numerically, we find that there exist two fixed points in this special case. The positions of these two fixed points in phase space are overlapped at (6.869, π/2, 0, 0), but their impact parameters are Lyapunov central theorem, one can find that each purely imaginary eigenvalue yields a Lyapunov family [20], which is a one parameter family γ ǫ of periodic orbits and the orbit γ ǫ collapses into the fixed point as ǫ → 0.
In Fig. 4, we show the Lyapunov family for the fixed points (light rings) in the plane (X, θ), where X is a compacted radial coordinate defined as X = r 2 − r 2 h /(1 + r 2 − r 2 h ) as in Ref. [20]. The orbits inside the unstable invariant manifold tube can reach the observer from the region near the event horizon of black hole. Moreover, the periodic orbit touched the boundary of the black region approaches perpendicularly to the boundary V (r, θ) = 0, which is similar to that in Ref. [15]. In Fig.5(b), we present the Poincaré section in the plane (θ, p θ ) for the unstable manifold (green) with η = −2 at the observers position.
The self-similar fractal structures in the unstable manifold are responsible for the self-similar fractal behavior in black hole shadow caused by chaotic lensing. The photons starting within the green regions always move only in the unstable manifold tube. In Fig.5 (b), the intersections of the unstable manifold with the dashed line denotes the trajectories which can be detected by the observer with the corresponding inclination angle.
These intersection points also determine the positions of the photons on the image plane and the boundary of the black hole shadow. In Fig. 5 (c) and ( image in observer's sky (13) are determined by the locally measured four-momentum pμ of a photon at the observer. This is also the reason why the shadow radius of a Schwarzschild black hole at spatial infinite is 3 √ 3M rather than the radius of photon sphere 3M . From the previous discussion, we know that the black hole shadow is exactly determined by the unstable invariant manifold associated with the fixed points. In Fig. (7), we present Poincaré sections (θ, p θ ) for the unstable manifold (green) with η = −2 near the photon sphere and at spatial infinite observer, respectively. Although some similar fractal structures appear in both sections, it is obvious that their positions and sizes in the pattern depend on the position of Poincaré section. This is also shown in Fig.(8) in which we present the unstable invariant manifold (light blue line) between Poincaré sections near the photon sphere and at spatial infinite observer. Generally, the photon orbits in the unstable invariant manifold are irregular curves rather than straight lines, especially in the case the photon motion is chaotic. This leads to the difference between the basins of attraction near the photon sphere and the black hole shadow. In other words, black hole shadow depends on not only the behavior of light rays near the photon sphere, but also the propagation of photons between the photon sphere and the observer.

V. SUMMARY
In this paper we have studied the shadows of a Schwarzschild black hole surrounded by a Bach-Weyl ring.
The presence of Bach-Weyl ring leads to that the photon dynamical system is non-integrable and then chaos would appear in the photon motion. Making use of the backward ray-tracing method, we present numerically the black hole shadow. Our result indicate that the ring mass and the Weyl radius affects heavily the black hole shadow. The size of the black hole shadow increases with the ring mass M and decreases with the Weyl radius b. However, the change of the black hole shadow shape with the ring mass and the Weyl radius becomes more complicated. With the increase of the ring mass M, we find that the shadow becomes gradually prolate along the axis of symmetry, but turns concave in the equatorial plane, which yields that the black hole shadow becomes a "8" type shape in the case with the larger ring mass M. In the case with the smaller Weyl radius, the shadow is a oblate silhouette and it is convex in the equatorial plane. Moreover, we also find some bright dispersion points in the black hole shadow, which possesses self-similar fractal structures originating from the chaotic lensing. With the increase of b, the effect chaotic lensing becomes more distinct so that some bright strips appear in the black hole shadow, which yields the emergence of some eyebrow-shape shadows. But with the further increase of b, the bright dispersion points with self-similar fractal structures disappears and black hole shadow becomes gradually concave in the equatorial plane, and finally, the black hole shadow becomes a prolate silhouette with the "8" shape again in the cases with the larger b. We also study the dependence of the shadows on the inclination angle of the observers as M = 1.5M and b = 3M . When θ obs = 0 • , there are some concentric bright rings imbedded in the black disc, which is qualitatively different from that in the case of usual spherical black hole. With the increase of inclination angle θ obs , the eyebrow shape shadows appears more distinctly in this case. The presence of Bach-Weyl ring also changes the shape of Einstein ring. For some selected ring mass and Weyl radius, the Einstein ring is broken and then the Einstein cross appears.
Due to the combined action of the black hole and the Bach-Weyl ring, some parts of Einstein cross could lie inside Bach-Weyl ring and the other is located outside the ring.
Moreover, we also study the change of the image of Bach-Weyl ring with the ring mass and the Weyl radius.
In the spacetime of a Schwarzschild black hole with a Bach-Weyl ring, the image of the Bach-Weyl ring present a flying saucer shape, which distributes symmetrically in the both sides of the equatorial plane. With the increase of the ring mass M, the saucer shape of the image becomes wider for the Bach-Weyl ring. Finally, we analyze the invariant manifolds of certain Lyapunov orbits near the fixed point and discuss further the formation of the shadow of a Schwarzschild black hole with Bach-Weyl ring, which implies that the shadow of a Schwarzschild black hole with a Bach-Weyl ring is exactly determined by the unstable invariant manifold associated with the fixed points.