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.


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 geomea e-mail: csb3752@hunnu.edu.cn b e-mail: jcwang@hunnu.edu.cn c e-mail: jljing@hunnu.edu.cn try 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 Kerrlike 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 non-planar 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, 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 Sect. 2, 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 Sect. 3, we adopt the backward ray-tracing technique and present numerically the shadows for the Schwarzschild black hole surrounded by Bach-Weyl ring. In Sect. 4, we analyze the invariant phasespace structures and further explain the formation of shadow cast by a Schwarzschild black hole with a Bach-Weyl ring. Finally, we present a summary.

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 ) λ int,z = 2ρ(ν 1,ρ ν 2,z + ν 1,z ν 2,ρ ).
Here, we focus on the spacetime generated by a Schwarzschildtype 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 and 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.
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][12][13][14][15][16][17][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 fourmomentum p μ onto eμ, i.e., With these quantities, one can obtain the coordinates of photon's image in observer's sky [10][11][12][13][14][15][16][17][18]20,56] 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 the larger ring mass M. Moreover, the size of the black In Fig. 1, we also marked the images of Bach-Weyl ring in magenta. It is shown that 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. The presence of Bach-Weyl ring also affects the shape of Einstein ring. Especially, in the cases where the ring mass is set to M = M or 1.5M and the Weyl radius is set to b = 10M or 20M, the Einstein ring is broken into four parts, and then the Einstein ring turns into the so-called Einstein cross. From Fig. 1, due to the combined action of the black hole and the Bach-Weyl ring, for some selected b, some parts of Einstein cross lie inside Bach-Weyl ring and the other is located outside the ring. In Fig. 2, we show the change of the black hole shadow on the Weyl radius b of Bach-Weyl ring for fixed ring mass M. When the Weyl radius b is small, the shadow is a oblate silhouette and it is convex in the equatorial plane, but there exist 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, 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. 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. Finally, the black hole shadow becomes a prolate silhouette with the "8" shape again in the cases with the larger b. 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 that the shadow is center symmetric. Moreover, we find that there are concentric bright rings imbedded in the black disc, which is qualitatively different from that in the case of usual spherical black hole. These distinct features in the shadow can be attributed to the effect of Bach-Weyl ring on the space- time structures. For the observer with the other inclination angle θ obs , the center symmetry of shadow is broken and the shadows becomes axially symmetric. With the increase of inclination angle θ obs , we find the eyebrow shape shadows appears more distinctly. And finally, as θ obs = 90 • , the pattern of shadow owns the symmetry along the equatorial plane.

Invariant phase-space structures and the shadow cast by a Schwarzschild black hole surrounded by a Bach-Weyl ring
Let us now in position to analyze the invariant phasespace structures of the photon system in the spacetime of a Schwarzschild black hole with a Bach-Weyl ring and then further explain the formation of the black hole shadow. The invariant phase space structures contain fixed points, periodic orbits and invariant manifolds. As one of the most important features for dynamical systems, the invariant phase space structures are applied extensively in the design of space trajectory for various of spacecrafts including a low energy transfer from the Earth to the Moon and a "Petit Grand Tour" of the moons of Jupiter [57][58][59][60][61]. Recently, it is shown that the invariant phase space structures play an important role in the formation of black hole shadows [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), 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. 2b. 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 η 1 = −14.29 and η 2 = 14.29, respectively. This distribution of two fixed points is determined by a fact that the Schwarzschild black hole with a Bach-Weyl ring (3) is a non-rotating spacetime. The further analysis indicates that the eigenvalues of the Jacobian in the equation (16) are λ 0 = ±0.42779, ν = ±0.42779 i. From Lyapunov central theorem, one can find that each purely imaginary eigen-value 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 [20]. The thick dot represent the light ring, and the solid lines denote the corresponding family of periodic Lyapunov orbits arising from these light rings. These periodic orbits can be parameterized by the impact parameter η = L z /E on the interval [−14.29, 14.29]. The periodic orbits with the same absolute value of η is overlapped in the plane (X, θ) because that the spacetime (3) is static and axially symmetric. All of these periodic Lyapunov orbits in Fig. 4 are nearly spherical orbits with radius r = 6.869, which play an important role in determining the boundary of shadow of a Schwarzschild black hole with a Bach-Weyl ring (3). The sign of real eigenvalue λ 0 determinate the stability of the invariant manifold associated with the fixed point. For the positive real eigenvalue, the invariant manifold is unstable and points in the manifold exponentially approach the fixed point in backward time. The invariant manifolds for each Lyapunov orbit are two dimensional surfaces forming tubes in the three dimen-sional reduced phase space (r ; θ ; p θ ). In Fig. 5a, we present a projection of the unstable invariant manifolds ( the green lines) related to the periodic orbit ( the red line) for η = −2 in the plane (X, θ) for the Schwarzschild black hole with a Bach-Weyl ring with M = 1.5M and b = 3M. 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. 5b, 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. 5b, 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. 5c, d, we  Fig. 5, we find that the bright stripes and bright region in Fig. 5c, d originate from the white regions in the Poincaré section, and the fractal-like structure shown in Fig. 5b accounts for the fractal shadow structure in Fig. 3. These correspondence shows that the shadow of the Schwarzschild black hole with a Bach-Weyl ring is exactly determined by the unstable invariant manifold associated with the fixed points, which is consistent with those in the cases of Kerr black hole with scalar hair [20] and of the Bonnor black dihole spacetime [17]. It is well known that photon sphere plays an important role in the formation of black hole shadow. In the Schwarzschild spacetime, the light rays which do not form the photon sphere are either escape to infinity or get captured by the black hole. In this sense, the photon sphere can be treated as a basin boundary that separates the basins of escape and capture of the light rays. Thus, the motion of light rays near the photon sphere are crucial to understanding the properties of black hole shadows.
Recently, Shoom et al. [62,63] studied the properties of photon sphere of a Schwarzschild black hole distorted by an external gravitational field with a quadrupole moment. It is shown that there exists a fractal basin boundary due to the chaotic behavior of null geodesics around the distorted black hole. In the left panel in Fig. 6, we present the basins of attraction near the photon sphere for a Schwarzschild black hole surrounded by a Bach-Weyl ring with M = 1.5M and b = 3M. Black points correspond to the trajectory got captured by the black hole and magenta points denote the trajectory falling into the Bach-Weyl ring. While, blue and green points correspond to trajectories escaped from the black hole in the upward and downward directions, respectively. Obviously, the basin boundary near the photon sphere possesses some fractal structures, which is similar to that obtained in the distorted Schwarzschild black hole [62].
From Fig. 6, we find that basins of attraction near the photon sphere possess many similar structures to black hole shadow at observer for a Schwarzschild black hole surrounded by a Bach-Weyl ring with M = 1.5M and b = 3M. It shows that photon sphere is very important in the formation of black hole shadow. However, we also note that there exist some distinct differences between the basins of attraction near photon sphere and the black hole shadow. For example, some bright dispersion points appear in the middle of the black hole shadow, but vanish in the basins of attraction near photon sphere. Moreover, the magenta points corresponded to the trajectory falling into the Bach-Weyl ring have the distinct distribution regions in the both cases. The main reason is that the black hole shadow also depend on the propagation of light ray in the spacetime and the position of observer because the coordinates of photon's 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.

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 of 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.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: The datasets generated during and/or analysed during the current study available from the corresponding author on reasonable request.] 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 .