Influence of quintessence dark energy on the shadow of black hole

We investigate the effects of quintessence dark energy on the shadows of black hole, surrounded by various profiles of accretions. For the thin disk accretion, the images of the black hole comprises the dark region and bright region, including direct emission, lensing rings and photon rings. Although their details depend on the form of the emission, generically, the direct emission plays a major role for the observed brightness of the black hole, while the lensing ring makes a small contribution and photon ring makes a negligible contribution. The existence of cosmological horizon also plays an important role in the shadows, since the observer in the domain of outer communications is nearby the cosmological horizon. For the spherically symmetric accretion, the static and infalling matters are considered. We find that the positions of photon spheres are the same for both static and infalling accretions. However, the observed specific intensity of the image for the infalling accretion is darker than the static accretion, due to the Doppler effect of the infalling movement.


Introduction
The Event Horizon Telescope (EHT) collaboration has recently released an image of a supermassive black hole in M87 [1][2][3][4][5][6]. The image shows that there is a dark interior surrounding by a bright ring, which are called the black hole shadow and photon ring, respectively. The shadow and photon ring are resulted from the deflections of light, or gravitational lensing by the black holes, a prediction from Einstein's General Relativity (GR) [7]. Therefore, the shadow image of black hole in M87 from EHT is another strong evidence for GR.
The deflection of light by an intense star or black holes was initiated in [8,9]. Later, it was extended to a black hole surrounded by a thin accretion disk [10]. According to these papers, there exists a critical curve in the observed image of the black hole. When traced backwards from the distant observer, the light ray from this critical curve will asymptotically approach a bound photon orbit. This bound orbits for Schwarzschild black hole is r = 3M (M is the mass of the black hole), and the critical curve has the radius of b = √ 27M (b is called impact parameter). Thus, the shadow of the black hole represents the dark interior of the critical curve. However, according to [11], the radius and the profile of the shadow would depend on the location and form of the accretion surrounding the black hole. Besides the shadow, there also exists the enhanced brightness, called photon ring, around the black hole. This is because the light rays near the critical curve may orbit many times of the black hole, thus picking up extra brightness. So far, the study of shadows of black holes has become a subject of great interest . In particular, the shadows of a high-redshift supermassive black hole may serve as the standard rulers [22], which would constrain the cosmological parameters. Moreover, the EHT data may impose constraints in particle physics via the mechanism of superradiance [24,25]. More references can be referred to [15].
Recent astronomical observations show that the Universe is accelerating expansion [40][41][42], implying a state of negative pressure. Quintessence dark energy is one of the candidates to interpret the negative pressure. The state equation of quintessence dark energy is governed by the pressure p and the energy density ρ with p = wρ, where w is the state parameter. It is deemed that the range −1 < w < −1/3 causes the acceleration of the Universe. Black holes in quintessence dark energy models were extensively studied [43]. To investigate the shadows of black hole in the quintessence dark energy model will bring new insights and impose restrictions to the quintessence model. Previous work on shadows in dark energy model can be found in [44,45].
In this paper, we will investigate the shadow images of quintessence black holes with different accretions. In particular, we study the geometrically thin and optically thin disk accretion and spherically symmetric accretion. Apparently different accretions will lead to different shadow images. Specifically, for the thin disk accretion, the images of the black hole comprises the shadows, lensing rings and photon rings. The laters are resulted from the intersections of the light rays and the thin accretion for many times. Although different emission profiles of the accretion is important to the observed specific intensities, the direct image from accretion makes major contributions to the brightness of the image. Besides, the lensing ring makes small contributions and the photon ring makes negligible contributions to the observed specific intensities. For the spherically symmetric accretions, the shadows and photon rings are also spherically symmetric. The static accretion and an infalling accretion will have different impacts to the observed specific intensities, although the radii of the shadows are not changed. Specifically, the infalling accretion will have extra Doppler effect, which will lead to darker shadows than the static case. A key point in the cosmological black hole is that there exists a domain of outer communication between the event horizon and cosmological horizon [46,47]. Therefore, various quintessence state parameter w will lead to different distances between the event horizon and cosmological horizon. This will certainly have impact on the distance between the observer (near cosmological horizon) and the accretion (near the event horizon). Therefore, different quintessence state parameters will lead to different forms of the observed specific intensities. We expect that the shadow images of the quintessence black hole will impose constraints of the quintessence dark energy model in the Universe from the observation of EHT. This paper is arranged as follows: In section 2 the black hole solutions surrounded by quintessence matter and the light deflections by the quintessence black hole are introduced; In section 3 the shadow images of the black hole with thin disk accretion are given; In section 4, we show the images of the black hole with spherically symmetric accretions; Finally, we draw the conclusions and discussions in section 5.

Light deflection by a black hole in quintessence dark energy
The geometry of a black hole in a quintessence dark energy model can be expressed as where, in which, M is the mass of the black hole, a and w are the normalization factor and state parameter respectively. For the case of w = −1, the metric in Eq.(2.2) reduces to the Schwarzschild-de-Sitter black hole with which has both the event horizon r h and cosmological horizon r c . For w = − 2 3 , we can obtain the event horizon and cosmological horizon as where the restriction 8M a < 1 has been imposed. For w = −1/3, the metric in Eq.(2.2) reduces to the Schwarzschild black hole, which has only an event horizon. In this paper, we are interested in the quintessence dark energy model, in which the state parameter w is restricted in the region −1 < w < − 1 3 , and there always exists an event horizon and a cosmological horizon. The region between the two horizons is called the domain of outer communication [46,47], since any two observers in this region may communicate with each other without being hindered by a horizon.
The equation of state for the quintessence dark energy is in which p is the pressure while ρ is the energy density. In order to cause cosmic acceleration, the pressure of the quintessence dark energy p should be negative, thus the energy density ρ is positive if the normalization factor a is positive. In this paper, we will only focus on the effect of state parameter on the shadows of a quintessence black hole, and the we will set the normalization factor a = 0.05 and black hole M = 1 thereafter. In order to investigate the light deflection in the quintessence black hole, we need to find how the light ray moves in this spacetime. We will resort to the following Euler-Lagrange equation d ds in which s is the affine parameter of the light trajectory, the symbol˙is the derivative with respect to s, thusẋ µ the four-velocity of the light ray and L is the Lagrangian, taking the form as (2.7) We will also impose the initial conditions θ = π/2,θ = 0. That is, the light rays always move in the equatorial plane. In addition, since the metric coefficients do not depend explicitly on time t and azimuthal angle ψ, there are two conserved quantities correspondingly. Combining Eqs.(2.2), (2.6) and (2.7), the time, azimuthal and radial components of the four-velocity satisfy the following equations of motionṡ where + and − in Eq.(2.9) correspond to the counterclockwise and clockwise direction respectively for the motion of light rays. The parameter b = |L|/E is the impact parameter, in which E and L are the conserved quantities with respect to the time and azimuthal direction. Moreover, Eq.(2.10) can be rewritten aṡ is an effective potential. At the photon sphere, the motions of the light ray satisfyṙ = 0, andr = 0, which also means where the prime denotes the derivative with respect to the radial coordinate r. Based on Eq.(2.13), we can obtain the radius and impact parameter of the photon sphere for different w. Generically, it is difficult to obtain the analytic results of the radius and impact parameter. But for w = −2/3 the metric function can be simplified, thus the radius and impact parameter of the photon sphere can be expressed as For other w's, the numerical results of radii and impact parameters of the photon sphere, as well as the event horizons and cosmological horizons are listed in the Table 1. From this table, we see that as one increases the absolute value of w, the event horizon r h and impact parameter b ph increase as well, while the cosmological horizon r c decreases. The radius of photon sphere however is non-monotonic, it increases at first and then decreases around w = −0.7. In addition, from Table 1 we see that the locations of the event horizon and cosmological horizon become closer as |w| increases, which leads to a narrower domain of the outer communication.
From the Eq.(2.11), we see that the motion of the particle depends on the impact parameter and the effective potential. In Figure 1, the effective potential is plotted for different state parameters w = −0.5 and w = −0.7. We can see that at the event horizon, the effective potential vanishes. Then it increases and reaches a maximum at the photon sphere r ph , and later decreases and vanishes at the cosmological horizon. Let's consider a light ray moves in the radially inward direction. In Region 1, if the light ray starts its motion at r > r ph , the light ray will encounter the potential barrier and be reflected back to the outward direction. If the photon starts the motion at r < r ph , the photons will fall into the singularity. In Region 2, namely b = b ph , as the light ray approaches the photon sphere, it will revolve around the black hole infinitely many times since the angular velocity is non-zero. Actually, this orbit is circular and unstable [11]. In Region 3, the light ray will continue moving in the inward direction since it does not encounter the potential barrier. Finally it will enter the inside of the black hole and fall into the singularity.
Since we have already obtained the effective potential, it is also convenient to obtain the effective force on the photon, which is 1 The first term −3M /r 4 , called the Newtonian term, is attractive as it is negative while the second term 1/r 3 is repulsive. The contributions of the quintessence dark energy comes from the third term − 1 2 3a(w + 1)r −3w−4 , which plays the role of an attractive force. The trajectory of the light ray can be depicted on the basis of the equation of motion. Combining Eqs.(2.9) and (2.10), we have (2.17) In order to integrate it out conveniently, we set u = 1/r. Thus, Eq.(2.17) becomes The geometry of the geodesics depends on the roots of the equation Φ(u) = 0. In particular, for b > b ph , light will be deflected at the radial position u i which satisfies Φ(u i ) = 0. Therefore, it is important to find the radial position u i in order to obtain the trajectory of the light ray. In addition, the location of the observer is also important. Usually, the observer is located at infinite boundary for the asymptotically flat spacetime. However, for the model of quintessence black hole in our paper, the cosmological horizon is present. Physically, the observer should be located in the domain of outer communication which is between the event horizon and cosmological horizon, just as in de Sitter spacetime. For convenience, the observer in our paper is set to be near the cosmological horizon. By virtue of the Eq.(2.18), we can solve the trajectories of the light rays, shown in Figure 2. All the light rays approach the black hole from the right side. The green, red and blue lines correspond to b < b ph , b = b ph and b > b ph respectively. Specifically, for the case of b < b ph , the light rays (green lines) drop all the way into the black hole. For b > b ph , the light rays (blue lines) are deflected but never enter the black hole. In particular, the light rays close to the black hole can even be reflected back to the right side, i.e., the side they approach the black hole. For b = b ph , the light rays (red lines) revolve around the black hole. This conclusion is consistent with those from the analysis of effective potential in Figure 1. Indeed, region 1, 2, 3 in Figure 1 correspond to the blue, red and green lines in Figure 2, respectively.
We need to emphasize that due to the existence of the cosmological horizon, different w's will result in different geometries of the light rays. We have assumed that the observer is in the domain of the outer communications and sitting close to the cosmological horizon. Therefore, for the case of w = −0.5 (left panel in Figure 2), the cosmological horizon is much far from the event horizon (see Table 1 ). Therefore, the observed entering light rays are approximately parallel, which can be seen from the left panel of Figure 2. However, for the case of w = −0.7, the cosmological horizon is close to the event horizon (see Table  1). Therefore, the observed entering light rays behave significantly different from those of w = −0.5 (see the right panel of Figure 2), i.e., the entering light rays for w = −0.7 are not parallel.
As the light is deflected, a shadow of the black hole will be formed. The light is emitted by the accretion matters, the profiles of the accretion matters thus are important for the shadows of the black holes. Next, we will investigate the effect of the profiles of the accretions on the shadows of the black holes. We will take the thin disc accretion and spherical accretion as two examples.

Shadows, photon rings and lensing rings with thin disk accretion
In this section, we consider an optically and geometrically thin disk accretion around the black hole, viewed face-on. From [11], we know that an important feature of thin disk accretion is that there are photon rings and lensing rings surround the black hole shadow. The lensing ring consists of light rays that intersect the plane of the disk twice outside the horizon, while the photon ring 2 comprises light rays that intersect the plane of the disk three or more times. Thus, the trajectory of the photon is important to distinguish the photon ring and lensing ring.

Number of orbits of the deflected light trajectories
The trajectories of the light rays are shown on the right column of Figure 3 for w = −0.5 (top row) and w = −0.7 (bottom row), respectively. In these figures, the red lines, blue lines and green lines are direct emission, lensing rings and photon rings, respectively. This can be understood from the definition in [11] that these light rays intersect the disk plane once, twice and more than twice. It should be noted that due to the domain of outer communications, the light rays going towards the observer will behave differently for different state parameters w. For the case of w = −0.5, the observer locates near the cosmological horizon which is far away from the event horizon (see Table 1). Thus, the light rays nearby the observer is almost parallel, which is similar to the Schwarzschild black hole [11]. However, for w = −0.7 the cosmological horizon is near the event horizon (refer to Table 1). Therefore, the light rays going towards the observer are not parallel, which is shown in the bottom row of right panel of Figure 3. One can also use the total number of orbits, defined by n = ψ/2π, to distinguish the trajectories of light rays [11]. The total number of orbits are shown in the left column of Figure 3 for w = −0.5 (top row) and w = −0.7 (bottom row), respectively. We still use the red lines, blue lines and green line to represent the direct emission, lensing rings and photon rings, respectively. Therefore, from the definitions of these light rings we can explicitly see that the direct emissions correspond to n < 3/4, while lensing rings correspond to 3/4 < n < 5/4 and photon rings correspond to n > 5/4. Lensing ring : 3/4 < n < 5/4, 5.71774 < b < 5.97314 and 6.04784 < b < 7.76682 Photon ring : n > 5/4, 5.97314 < b < 6.04784 Direct emission : n < 3/4, b < 6.92848 and b > 8.36103 Lensing ring : 3/4 < n < 5/4, 6.92848 < b < 7.2156 and 7.28972 < b < 8.36103 Photon ring : n > 5/4, 7.2156 < b < 7.28972

Observed specific intensities and transfer functions
Next, we will investigate the observed specific intensity of the thin disk accretion. We assume that the thin disk emits isotropically in the rest frame of static worldlines. As in [11], we take the disk to lie in the equatorial plane of the black hole, with a static observer at the North pole. We denote the emitted specific intensity and frequency as I e (r) and ν e , the observed specific intensity and frequency as I obs (r) and ν. Since I e /ν 3 e is conserved along a light ray from Liouville's theorem, the observed specific intensity thus can be written as The total specific intensity is obtained by integrating specific intensity with different frequencies where I em (r) = I e (r)dν e is the total emitted specific intensity near the accretion. If a light ray traced backwards from the observer intersects the disk, it will pick up brightness from the disk emission. For 3/4 < n < 5/4, the light ray will bend around the black hole and hit the opposite side of the disk from the back (see the blue lines in Figure  3). Thus, it will pick up additional brightness from this second passage through the disk. For n > 5/4, the light ray will bend around the black hole more, and then hit the front side of the disk once again (see the green lines in Figure 3). This leads to additional brightness from the third passage through the disk. Therefore, the observed intensity is a sum of the intensities from each intersection, that is where r m(b) , called transfer function, is the radial position of the m th intersection with the disk plane outside the horizon. For simplicity, we have neglected absorption of the light for the thin accretion, which may decrease the observed intensity from the additional passages. The transfer function describes the relation between the radial coordinate r and the impact parameter b. The slope of the transfer function, dr/db, is the demagnification factor. For different state parameters w, the transfer functions with respect to impact parameter b are shown in Figure 4. The black line, orange line and red line represent the first (m = 1), second (m = 2) and third (m = 3) transfer functions. The first transfer function corresponds to the direct image of the disk, and its slope is about 0.85 (0.67) for w = −0.5 (w = −0.7). Therefore, the direct image profile is actually the redshifted source profile. The second transfer function corresponds to the lensing ring (including the photon ring). In this case, the observer will see a highly demagnified image of the back side of the disk. Finally, the third transfer function corresponds to the photon ring. In this case, one will see an extremely demagnified image of the front side of the disk since the slope is about infinite. It thus contribute negligibly to the total brightness of the image.

Shadows, lensing rings and photon rings
Having obtained the transfer functions, we can further obtain the specific intensity on the basis of Eq.(3.5) as the emitted specific intensity is given. In this paper, we assume the following toy-model functions which may emit by some thin matters. Firstly, we assume that I em (r) is a decay function of the power of second order, This profile is plotted in the top left panels in Figure 5 for w = −0.5 and Figure 6 for w = −0.7. This emission is sharply peaked at r = 6M , and then drops to zero abruptly as r < 6M . In this example, the region of emission is well outside the photon sphere (see Table 1). Secondly, we assume that the emission is a decay function of the power of third order as But this time the peak of the mission is at the photon sphere r ph . The profiles of this emission can be seen from the middle left panels in Figure 5 for w = −0.5 and Figure 6 for w = −0.7. Finally, we assume a moderate decay of emission as This emission is outside of the horizon r h , which is shown in the bottom left panels in Figure 5 for w = −0.5 and Figure 6 for w = −0.7.
For the case of w = −0.5, the observed specific intensities are shown in the middle and right columns in Figure 5. The middle column is about the one dimensional functions of I(r) with respect to b, while the right column is about the two dimensional density plots of I(r), viewed faced-on. In the first row the emission from the accretion is well outside the photon sphere (left panel). The corresponding observed direct image decays similarly as b > 7.5M due to the gravitational lensing (middle panel). However, the observed lensing ring emission is confined to a thin region 6.4M < b < 7.4M , since the back side image is highly demagnified for r > 6M . The photon ring emission is the spike at b ∼ 6M , which is hardly to see in the right panel of the density plots. One needs to zoom in the plot to see the photon ring. Therefore, from the observer, the lensing ring makes a small contribution to the total brightness while the photon ring makes negligible contributions.
In the second row of Figure 5, the emission starts from the photon sphere r ph = 3.216. The most important difference from the first row is that the observed lensing ring and photon ring emission are now superimposed on the direct image emission, which decays from b > 4.25M . The lensing ring has a spike in the brightness in the range of 6M < b < 6.5M , while the photon ring has an even narrower spike at b ∼ 6M which is hardly to be distinguished from the lensing ring. In this case, the lensing ring still makes a very small contribution to the total brightness and the photon ring makes a negligible contribution.
Finally, in the third row of Figure 5, the emission extends to the event horizon r h = 2.158 and the decay of the emission is much moderate compared to the first and second row. The inner edge of the observed intensity at b ∼ 3.2M is the lensed position of the event horizon. The increase of the intensity outside of the central dark area is due to the gravitational redshift. The very narrow spike at b ∼ 6M is the photon ring, while the broader bump at b ∼ 6.2M is the lensing ring. In this case, the lensing ring makes a prominent brightness to the observed intensity, but the photon ring is still entirely negligible.
From Figure 5, we see that the dark central regions are different for different profiles of the emissions. But the locations of the photon ring always stay at b ∼ 6M . In addition, for all the cases, we find that the observed brightness mainly stems from the direct emission, while the lensing ring provides only a small contribution, and the photon sphere provides a negligible contribution.
For the case of w = −0.7, we plot the profiles of the emission and observed intensities in Figure 6. They behave similarly to those of w = −0.5. The differences are the quantities of the intensities, as well as the locations of photon ring and lensing ring.

Shadows and photon spheres with spherical accretions
In this section, we are going to investigate the shadows and photon spheres with the spherical accretions for various state parameters w.

Shadows and photon spheres with a static spherical accretion
In this section, we intend to investigate the shadow and photon sphere of a quintessence black hole with a static spherical accretion. We mainly focus on the specific intensity observed by the observer (ergs −1 cm −2 str −1 Hz −1 ), which can be expressed as [12,14] in which g = ν o /ν e is the redshift factor, while ν e is the radiated photon frequency and ν o is the observed photon frequency, j(ν e ) is the emissivity per unit volume measured in the static frame of the emitter, dl p is the infinitesimal proper length, and γ stands for the trajectory of the light ray.
In the quintessence black hole (2.1), the redshift factor g = f (r) 1/2 . In addition, we assume that the radiation of light is monochromatic with fixed a frequency ν s , thus the specific emissivity takes the form as where we have adopted the 1/r 2 profile as in [14]. According to Eq.(2.1), the proper length measured in the rest frame of the emitter is in which dψ/dr is given by the inverse of Eq. (2.17). In this case, the specific intensity observed by the static observer is Next we will employ Eq.(4.4) to investigate the shadow of the quintessence black hole with the static spherical accretion. Since the intensity depends on the trajectory of the light ray, which is determined by the impact parameter b, we will investigate how the intensity varies with respect to the impact parameter b. The observed specific intensities are shown in Figure 7 with w = −0.5 (top row) and w = −0.7 (bottom row). From Figure 7, we see that as b increases, the intensity ramps up first, then reaches a peak at b ph (b ph = 5.98M for w = −0.5 and b ph = 7.23M for w = −0.7), and finally drops rapidly to smaller values. This result is consistent with Figure 1 and Figure 2. Because for b < b ph , the intensity originated from the accretion matter is absorbed mostly by the black hole, so the observed intensity is little. For b = b ph , the light ray revolves around the black hole several times, thus the observed intensity is maximal. For b > b ph , only the refracted light contributes to the observed intensity. In particular, as b becomes larger, the refracted light becomes less. The observed intensity thus vanishes for large enough b.
From Figure 7, we can also see how the state parameter of dark energy will affect the observed intensity. For a fixed a and M , we find that the larger the state parameter is, the stronger the intensity will be, i.e., the intensity of w = −0.5 is stronger than that of w = −0.7. The 2-dimensional plot of the shadows are shown in the right column of Figure  7, with the impact parameter b of the radius. We can see obviously that the shadow is circularly symmetric for the quintessence black hole, and outside the black hole, there is a bright ring, which is the photon sphere. The radius of the photon sphere for different w's have been listed in Table 1. Obviously, the result in Figure 7 is consistent with that in Table 1. That is, the larger the state parameter is, the smaller the radius of the photon sphere will be.
In addition, inside the photon sphere, we see that the specific intensity does not vanish but has a small finite value. The reason is that part of the radiation is escaped from the black hole. For r < r ph , the solid angle of the escaping radiation is 2π(1 − cos θ), while for r > r ph , the solid angle is 2π(1 + cos θ), in which θ is given by (4.5) In principle, one can calculate the escaped net luminosity by the following relation Obviously, for different state parameters, the net luminosity observed by the static observer in the quintessence black hole is different, since solid angle is dependent on the state parameter. This result is consistent with the result in Figure 7.

Shadows and photon spheres with an infalling spherical accretion
In this section, the optically thin accretion is supposed to be infalling matters. This infalling model is thought to be more realistic than the static accretion model, since most of the accretion matters are dynamical in the universe. We still employ Eq.(4.4) to investigate the shadow cast by the infalling accretion. Different from the static accretion, the redshift factor for the infalling accretion is related to the velocity of the accretion, that is in which k µ =ẋ µ is the four-velocity of the photon, u µ o = (1, 0, 0, 0) is the four-velocity of the static observer, and u µ e is the four-velocity of the accretion under consideration, given by The four-velocity of the photon has been obtained previously from Eq.(2.8) to Eq.(2.10). From these equations, we know that k t = 1/b is a constant, and k r can be inferred from the equation k α k α = 0, that is in which the +/− corresponds to the case that the photon gets close to/away from the black hole. With Eq.(4.9), the redshift factor in Eq.(4.7) can be simplified as which is different from the static accretion case. In addition, the proper distance is defined by where s is the affine parameter along the photon path γ. With regard to the specific emissivity, we also assume that it is monochromatic, so that Eq.(4.2) is still valid. Integrating Eq.(4.1) over the observed frequencies, we obtain I ∝ γ g 3 k t dr r 2 |k r | . (4.12) Now we will use Eq.(4.12) to investigate the shadow of the quintessence black hole numerically with an infalling accretion. Note that there is an absolute sign for k r in the denominator. Therefore, when the photon changes its motion direction, the sign before k r should also change. For different state parameter w's, the observed intensity with respect to b are shown in Figure 8, in which the top row is for w = −0.5 while the bottom row is for w = −0.7. From Figure 8, we find that as b increases the intensity will increase as well to a peak b = b ph , after the peak it drops to smaller values. This behavior is similar to that in the static accretion in Figure 7. We can also observe the effect of w on the intensity, and find that the intensity increases as the value of w increases, i.e., the intensity for w = −0.5 is stronger than w = −0.7. The two dimensional image of the intensity is shown on the right column in Figure 8. We see that the radii of the shadows and the locations of the photon spheres are the same as the static case. That means the motion of the accretion does not affect the radii of the shadows and the locations of the photon spheres. However, different from the static accretion, the central region of the intensity for the infalling accretion is darker, which can be accounted for by the Doppler effect. In particular, near the event horizon of the black hole, this effect is more obvious.

Conclusions and discussions
We investigated the images of a black hole in the quintessence dark energy model. In particular, we studied the quintessence black hole surrounded by a geometrically and optically thin disk accretion, as well as a spherically symmetric accretion. We found that for the observed specific intensities, there existed a dark interior (shadow), lensing ring and photon ring for the thin disk accretion. But their positions depended on the profiles of the emission from the accretion near the black hole. Despite of this distinctions, the direct image of the accretion made a major contribution to the brightness of the black hole, while the lensing ring made a minor contribution and the photon ring made a negligible contribution. The effects of the quintessence state parameter was also important to the shadow images of the black hole, since this parameter would have impact on the distances between the observer and the event horizon.
We also studied the static and infalling spherically symmetric accretions, the images of the black hole would have a dark interior and a photon sphere. The infalling accretion would have a darker interior than the static case, which was due to the Doppler effect of the infalling matter. Different quintessence state parameter would change the positions of the photon spheres of the image. We expect that this would bring insights and implications for the quintessence dark energy model from the observations of the shadow images of the black holes in future EHT experiments.