Investigating shadow images and rings of the charged Horndeski black hole illuminated by various thin accretions

In this paper, we investigate the shadows and rings of the charged Horndeski black hole illuminated by accretion flow that is both geometrically and optically thin. We consider two types of accretion models: spherical and thin-disk accretion flow. We find that in both types of models, the size of the charged Horndeski black hole shadow decreases with the increase of the charge, and it decreases more slowly for the Reissner–Nordström (RN) black hole. In the spherical accretion flow model, we find that the increase of the charge of Horndeski black hole brightens the light ring around it, and it brightens more significantly in comparison with RN black hole. Due to the Doppler effect, the charged Horndeski black holes with accretion flow of radial motion have darker shadows than those with the static accretion flow, but the size of the shadow is not affected by accretion flow motion. In the thin disk-shaped accretion flow model, we find that the brightness of the light ring around the charged Horndeski black hole is dominated by the direct emission from the accretion flow, and the contribution from lensed rings is relatively small, and that from the photon rings is negligible. We also find that the ring brightness decreases as the charge of Horndeski black hole increases, and the decrease is more significant than that in the RN black hole case. Moreover, the radiation position of the accretion flow can affect the shadow size and the ring brightness of the charged Horndeski black hole.


I. INTRODUCTION
Black holes are thought of as the most mysterious one of celestial bodies predicted by General Relativity (GR).Recently, the Laser-Interferometer Gravitational Wave-Observatory (LIGO) has successfully detected the gravitational waves from the merger of two black holes [1] and the Event Horizon Telescope (EHT) collaboration directly observed the ultra-high angular resolution images of supermassive black holes in the M 87 * [2-4] and Sagittarius A * [5,6].These observations strongly confirm the existence of black holes in our universe.Particularly, the black hole images present that there is a dark central region surrounded by a bright ring, which are called shadow and photon ring of the black hole, respectively.The formation of black hole shadow is due essentially to the deflection of light in the strong gravity field [7][8][9][10], its image may carry some valuable information of the geometry around the black hole.Consequently, the investigation of the black hole shadow may provide a new window to constrain different gravity models [11][12][13][14][15].
Black holes can obtain the angular momentum via the gravitational collapse of a massive star.Thus we should think the astrophysical black holes are rotating black holes in realistic astrophysical settings.Bardeen illustrated that the shadow shape of the Kerr black hole is the deformation due to the drag effect [16].Wang et al. investigated the shadow of a Konoplya-Zhidenko rotating non-Kerr black hole with an extra deformation parameter [17].Haroon et al. studied the effects of perfect fluid dark matter and a cosmological constant on the shadow of a rotating black hole [18].Wei et al. explored the nature of Gauss-Bonnet gravity by four-dimensional rotating black hole shadow [19].The shadow of the rotating black hole has also been comprehensively investigated with many interesting results [20][21][22][23][24][25][26][27].The energy and angular momentum for a rotating black hole can be extracted by the Penrose superradiance process [28] or the Blandford and Znajek mechanism [29], which would lead to a rotating black hole degenerating into a non-rotating one.Thus we cannot exclude the possibility that static black holes exist in the universe.The shadow of a static spherically symmetric black hole also is extremely important for studying the properties of the physics of a black hole, and for recent the literature see [30][31][32][33][34][35][36].It is well known that the accretion matter exists around a black hole, so their contribution is always very important and interesting, when we in detail study shadows and optical images of such black holes in the Universe.
Considering the astrophysical black hole is surrounded by a luminous accretion flow, a few pioneering works have investigated the nature of the shadow and photon ring of the black hole.Luminet investigated the optical appearance of a spherical black hole surrounded by thin accretion disk and pointed out that the shadow and photon ring observation characteristics are related to the position and profile of accretion flow [37].Bambi studied apparent images of the Schwarzschild black hole and the static wormhole to distinguish black holes from wormholes [38].Gralla et al. primarily investigated the simple case of emission from an optically and geometrically thin disk near a Schwarzschild black hole and proposed that the bright ring near the black hole shadow consists of the direct emission, lensed ring, and photon ring [39], which are determined by the times of the intersects between the light ray and the thin accretion disk.Narayan et al. explored a simple spherical model of optically thin accretion on a Schwarzschild black hole, and shown that the size of the observed shadow is hardly influenced by accretion flow [40].Cunha et al. studied the lensing and the shadow of the Schwarzschild black hole with a thin and heavy accretion disk [41].Zeng et al. investigated influence of the Gauss-Bonnet coupling parameter or quintessence dark energy on the optical appearance of the black hole with static/infalling spherical accretion flows [42,43].Guo et al. showed the feature of the observed shadows and rings of the Hayward black hole depend on the accretion flow property and the black hole magnetic charge [44].The photon ring and observational appearances of black holes have been extensively studied in [45][46][47][48][49][50][51][52][53][54][55], which disclosed the images of black holes surrounded by various accretions beyond GR.
Although GR successfully describes the gravitational interaction at the galactic and cosmological scales, it cannot give a satisfactory description for the accelerated expansion of our Universe [56,57].Therefore, it is generally thought that GR might need to be modified to give a satisfactory explanation of this phenomena [58].The scalar-tensor theories are regarded as a simplest nontrivial modification of GR, which contain a scalar field φ nonminimal couplings with a metric tensor g µν .Horndeski has constructed the most famous four-dimensional scalar-tensor theories in 1974 [59], which is called Horndeski gravity, by inspiration of the work of Lovelock [60].Horndeski gravity has been extensively studied in P-V criticality [61] and thermodynamics [62,63] of these black holes, holographic applications [64][65][66][67] and other interesting features (please see [68] for a review).In addition, Horndeski theory might be able to evade the no-hair theorems [69][70][71][72].
In this paper, we mainly focus on the shadows and photon rings of a charged Horndeski black hole with distinctive thin accretion flows, whose metric is obtained from the Horndeski theory with the scalar field coupled to the Einstein tensor in the presence of an electric field.Its thermodynamical properties were explored in [73], the weak and strong deflection gravitational lensing by a charged Horndeski black hole were analyzed in [74].A black hole can have a positive net electric charge due to the twisting of magnetic field lines [75] and the balance between the Coulomb and gravitational forces for charged particles near the surface of the compact object [76].Hence, it is very important to study the properties of charged black holes.By the investigation of this paper, we find that the influence of the electric charge of the Horndeski black hole on event horizon, photon sphere, shadow size and brightness of light ring are more significantly in comparison with Reissner-Nordström (RN) black hole cases as the charge increases.If future observations show that black hole shadows are produced by charged black holes, this result may have important signification for verifying GR and Horndeski theory.
Our paper is organized as follows: In Sec.II, we discuss the effective potential and shadow radius of the charged Horndeski black hole, and study the trajectories of photons by utilizing the ray-tracing method.In Sec.III and IV, we consider that the charged Horndeski black hole illuminated by spherical and thin disk accretion flows, respectively, and explore the shadow contour, photon rings as well as the corresponding observed luminosity for a distant observer.The Sec. V is the conclusions and discussions.Throughout this paper we use the geometric units with G = c = 1.

II. NULL GEODESIC AND SHADOW RADIUS OF THE CHARGED HORNDESKI BLACK HOLE
The most general action of Horndeski gravity has been constructed in [59].In our paper, we only investigate a special case in Horndeski gravity with the scalar field φ couple to the Einstein tensor G µν controlled by parameter η in the presence of an electromagnetic field F µν , and the corresponding action is written as [73,74,77] The variation of the action (2.1) with respect to the metric tensor, the scalar field, and the electric field yields respectively, where T µν and E µν are defined as In this paper, we focus on an asymptotically flat solution describing the charged Horndeski black hole to (2.2)-(2.4),while the metric, scalar field and Maxwell field are obtained [73,77] where where M denotes the black hole mass, Q is the total electric charge1 , and Ψ 0 is an integration constant.The (2.9) are consistent with those of the Schwarzschild black hole as Q = 0, while they can't reduce to the RN black hole metrics due to η can not be zero2 .Utilizing the (2.9), we analytical derive the event horizon radius (2.11) The Q need to satisfy the condition to ensure the existence of the event horizon [73,74] (2.12) The charged Horndeski black hole (2.9) is distinct from the RN black hole [78-80] by A(r)B(r) ̸ = 1 as well.Therefore, it is important that we analyze the different characteristics of two charged metric solutions via studying the shadow and optical appearance of the black hole.The Euler-Lagrange equation is given as follows: where λ is the affine parameter of the light trajectory, and ẋµ denotes the four-velocity of the photon.The Lagrangian L of the photon can be written as Taking the spherical symmetry into account, we only investigate, without loss of generality, the light traveling on the equatorial plane of black hole, i.e. θ = π/2 and θ = 0.The spacetime metric in (2.7) does not depend explicitly on time t and azimuthal angle ϕ.Hence, there are two conserved quantities in this spacetime correspond to the total energy and angular momentum of the photon, which read as For the trajectory of the light, one impose the null condition ds 2 = 0.According to (2.7) and (2.16), one can easily derive the following equations of motion of photon: ) where the affine parameter λ has been replaced with λL, the impact parameter b ≡ L/E.The effective potential V ef f (r) for the motion of a photon can be defined as larger electric charge leads to a smaller peak effective potential at the larger radius in the charged Horndeski black hole, and to the opposate situation in the RN black hole.In addition, the V ef f (r) is affected more significantly in the charged Horndeski black hole in comparison with the RN black hole case with the increase of the Q in Fig. 1.
The motion of the photon should satisfy ṙ = 0 and r = 0 at the photon sphere, which result in from (2.20) where b ph ≡ b(r ph ) is the critical impact parameter, and r ph is the radius of the photon sphere.From (2.21), one easily obtain Substituting (2.9) into (2.22) and (2.23), the sphere photon and the critical impact parameter are derived respectively as follows where Using (2.11), (2.24) and (2.25), the charged Horndeski black hole outside event horizon radius r + , the radius of the photon sphere r ph and the critical impact parameter b ph as the function of the electric charge Q are plotted respectively in Fig. 2. It clearly show that the r + , r ph and b ph are gradually decrease with the Q increasingly compared with the Schwarzschild back hole (Q=0) in the Q region (2.12), and the influences of the Q of the Horndeski black hole are more significant than those in the RN black hole cases.Based on (2.18) and (2.19), the trajectory of light can be given by

.27)
For convenience later, we introduce a parameter u ≡ 1/r, and then the (2.27) is rewritten as . (2.28) Utilizing the ray-tracing code, in Fig. 3, we have plotted the trajectories of the photons around the charged Horndeski black hole for different values of Q.In Fig. 3, the black dashed line represents the photon sphere (r = r ph ), and the projection of the event horizon surface is regarded as a black disk.On the green line (b = b ph ), the light rays will make a circular motion around the black hole infinitely many times without perturbations; the area above the green line (b > b ph ), the light rays will encounter the potential barrier (with the right hand of peak in Fig. 1) and then are scattered to infinite after passing through the turn point; the area under the green line (b < b ph ), the light rays will move in the inward direction (with the left hand of peak in Fig. 1) and are absorbed by the black hole eventually.This kind of light rays cannot be observed by the distant observer, so a shadow forms in the observational sky.Compared with the Schwarzschild back hole (Q=0), the Fig. 3 shows that the photon sphere and the black disk are smaller with the increase of the charge of the Horndeski black hole, expect the deflection angle increases, and their changes are more slowly for the RN black hole case.

III. SHADOWS AND RINGS OF THE CHARGED HORNDESKI BLACK HOLE ILLUMINATED BY DIFFERENT SPHERICAL ACCRETION FLOWS
There is plenty of free-moving matters in the Milky Way that might be accreted by its central black hole.The light rays radiated by these substances are the light source that illuminates the black hole.In this section, we will continue to investigate the properties of the shadows and photon rings of the charged Horndeski (or RN) black hole with static and infalling spherical accretion flows, by regarding them as optically and geometrically thin.

A. The static spherical accretion flow model
We firstly consider a static spherically symmetric accretion flow distributed outside the event horizon of a charged Horndeski (or RN) black hole.For an observer at infinity, the observed specific intensity (measured in erg s −1 cm −2 str −1 Hz −1 ) of photon with a frequency v s o is determined by [38,81] where v s em is the intrinsic photon frequency; g s is the red-shift factor, which can be defined by g s ≡ v s o /v s em ; dl prop represents the infinitesimal proper length; j(v s em ) is the emissivity per unit volume in the rest frame of the emitter.Integrating (3.1) over all the observed frequencies, the total observed intensity is given by From the four-dimensional spherical symmetric black hole spacetime (2.7), the red-shift factor can be written as g s = A(r) 1/2 .In addition, assuming a simple case that the radiation of light is monochromatic with fixed a frequency v f , thus the specific emissivity takes the form and the proper length measured in the rest frame of the emitter can be expressed as where dϕ/dr is given from (2.27).Putting (3.3) and (3.4) into (3.2), the total photon intensity measured by the distant observer can be written as The total observed intensities I s obs as a function of the impact parameter b are shown for several representative values of electric charge in Fig. 4. One can observe that the intensity curve sharply peaks at b ph , which corresponds to the photon rings.Meanwhile, the observed intensity is stronger as a larger Q in the charged Horndeski black hole, but the corresponding b ph gets smaller.We observe that the observed intensity is more obvious in the charged Horndeski black hole than those of the RN black hole with the increase of Q in Fig. 4. The Fig. 5 demonstrates that the two-dimensional shadows cast of the charged Horndeski (or RN) black holes in the celestial coordinates.One can see clearly that the bright ring outside the shadow of black hole is the photon sphere, here the luminosity is strongest.Furthermore, the specific intensity among dark shadow does not completely vanish but has a small finite value.The reason is that the tiny fraction of radiating gas behind the black hole can escape to infinity due to spherically symmetric accretion flow.

B. The infalling spherical accretion flow model
In realistic astrophysical settings, the substances with the intrinsic initial velocity may move toward the centre of celestial body.For example, the accretion flows fall the centre of the supermassive black hole in M 87 * [82].For simplicity, we consider a charged Horndeski (or RN) black hole that is surrounded by radial infalling spherical accretion flow in this section.In the scenario, the (3.2) still is valid to calculate the total observed intensity, but the red-shift factor is related to the velocity of the accretion flow, thus it rewritten as [38] where k µ is the four-velocity of the photon; u µi o is the four velocity of an observer; u νi em = (u it em , u ir em , 0, 0) is the four-velocity of the accretion flow.Utilizing (2.17)-(2.19),one easily obtain and then their ratio in which the symbol ± corresponds to the photon is approaching (+) or away (−) from the black hole.For the simplicity of analysis, we assume that the distant observer is stationary with u µi o = (1, 0, 0, 0), and the four-velocity of the infalling spherical accretion flow u νi em are given by respectively.For the scenario, the red-shift factor in (3.6) is rewritten as and the proper distance is expressed by Finally, the total observed intensity of the infalling spherical accretion flow in (3.2) is obtained dr.
(3.12) Fig. 6 plots that the total observed intensity I i obs of the charged Horndeski (or RN) black hole surrounded by an infalling accretion flow as a function of the impact parameter b.The features of the function curves are similar to the Fig. 4, except the intensity has an extremely sharp rise before the peak.Fig. 7 shows that the two-dimensional shadows cast on the celestial coordinates in this scenario.In the same way, the observed intensity and shadow size are increase faster with the increase of the Q in the charged Horndeski black hole, while it is slower in the RN black hole.One can intuitively see that the black hole shadows with state accretion flow of the Fig. 5 are brighter than those with the infalling accretion flow of the Fig. 7 due to the Doppler effect, but the size of the black hole shadows does not change with the same electric charge.The results manifest the size of the black hole shadow depends on the spacetime geometry, while the luminosity of the black hole shadow relies on the accretion flow model.

IV. SHADOWS AND RINGS OF THE CHARGED HORNDESKI BLACK HOLE ILLUMINATED BY THIN DISK ACCRETION FLOW
The compact celestial body (i.e.white dwarf, neutron star and black hole) can accrete various matters around it that forms an accretion disk in our Universe.There are some accretion disk models near the compact object have been investigated in [83][84][85][86].The accretion substances can be regarded as an light source that illuminates the compact object.Therefore, the investigation of the optical appearance around compact object may help us to constraint the accretion disk model and explore the distribution of matters around it.In this section, considering a simple model that the substances around the black hole are regarded as an optically and geometrically thin disk-shaped accretion disk [84], and one mainly focus on the emission from the accretion disk is isotropic in the rest frame of static world-lines, the disk is located on the equatorial plane and the observer is at the north pole.

A. Direct emission, lensed ring and photon ring
The literature [39] manifested that the black hole shadow with the thin-disk accretion is surrounded by photon rings and lensed rings.Based on the definition of the total number of light orbits, i.e. n ≡ ϕ/2π, the light rings near the black hole can be divided into three types [39]: (i) n < 3/4 is the direct emission, implying the light trajectories only intersect with the thin disk once.(ii) 3/4 < n < 5/4 is lensed ring, meaning the light trajectories of intersect with the the thin disk at twice.(iii) n > 5/4 is photon ring , denoting the light trajectories intersect the equatorial plane at least three times.

TABLE I:
The range of impact parameter corresponding to direct emission, lensed ring and photon ring of the charged Horndeski (or RN) black hole, where the black hole mass as M = 1, and the electric charge taking as Q = 0, 0.1, 0.3, 0.5.For various values of the Q, the ranges of the b of the direct emission, lensed ring, and photon ring are respectively shown in Table I.In addition, in order to depict these three regions in a more intuitive way, the Q = 0, 0.3, 0.5 are chosen as examples.Fig. 8 presented the total number of orbits as a function of the impact parameter, and Fig. 9 shows that the corresponding photon trajectories around the black hole as well.The black dashed circle is the photon sphere and the black disk stands for the black hole.Combining with Table I and Fig. 9, we can find that the ranges of b of the lensed rings and photon rings in the charged Horndeski black hole spacetime are narrower compared with the Schwarzschild black hole, meaning the thickness of the lensed rings and photon rings are getting thinner with existence of the electric charge.In addition, the b regions of the lensed rings and photon rings in the charged Horndeski black hole shrink faster than those of the RN black hole as a larger Q.

B. Total observed intensity and transfer functions
The specific intensity and frequency of the emitted light from the thin accretion disk are denoted as I e and v e .Based on the Liouville's theorem, the I e /(v e ) 3 is an invariant in the direction of light propagation, thus one can give where the I o is the observed specific intensity measured by the infinity observer, and v o = A(r)v e is corresponding to the red-shifted frequency.The observed intensity for a specific frequency can be expressed as Therefore, the total observed intensity can be obtained by an integral over all the observed photon frequencies where the I em (r) ≡ I e (r)dv e denotes the total emitted specific intensity on the thin accretion disk.One merely consider that the intensity of light is emitted by the thin disk accretion flows, and the absorption and reflection of light are ignored due to optically and geometrically thin.The light can pick up an additional brightness for each intersection between the light and the accretion disk on the equatorial plane.As previous discussions, for 3/4 < n < 5/4, the light ray will bend around the black hole and falls on the back of the accretion disk (See blue lines in Fig. 9).The light will pick up the additional brightness from the second crossing between the light and the accretion flow on thin disk.For n > 5/4, the light again arrives at the front side of the thin disk (See green lines in Fig. 9), which will add more brightness from the third passage through the thin disk.Consequently, the sum of the intensity at each intersection is the total observed intensity, which should be written as where the r n (b) is the so-called transfer function, representing the radial position of the n th intersection between the light and the thin disk.The slope of the transfer function dr/db describes the demagnification factor at each b [39].
The transfer functions as a function of the impact parameter are shown with the different values of the electric charge in Fig. 10.In Fig. 10, the red lines denote the first (n = 1) transfer function, corresponding to the direct image.The slope of the red curve is approximately equal to 1, thus the direct image profile is the gravity red-shifted source profile.Then, the blue lines are associated with the second (n = 2) transfer function, representing the lensed ring image.In the case, the image of the backside of the disk will be demagnified due to the slope of the blue curve much larger than 1.Finally, the green lines are the third (n = 3) transfer function, relating to the photon ring image.Because the slope of the green curve tends to be infinite for the case, the image of the front side of the disk will be extremely demagnified.Therefore, its contribution can be hardly ignored to the total brightness of the image.In addition, compared with the transfer functions of the Schwarzschild black hole, we clearly see that the electric charge results in toward left migration of the transfer functions in Fig. 10, while the effects of the Q on the transfer function in the charged Horndeski black hole are larger than those of the RN black hole in the region (2.12).

C. Observed appearance of the charged Horndeski black hole
Utilizing the relation between the transfer functions and the total observed intensity (4.4), we move on to investigate the optical appearance of the charged Horndeski (or RN) black hole.The total emission intensity I em (r) is determined by the emission position r on the accretion disk.Thus our aim is to analyze the influence of the accretion disk radiation position on the observation appearance of the black hole.To end this purpose, we parameterize the radiations intensity profile of the accretion disk as the recently introduced Gralla-Lupsasca-Marrone (GLM) disk model [87][88][89], which has been illustrated to be in a close agreement with the observational predictions of general relativistic magnetohydrodynamics simulations of astrophysical accretion disks.The radiations intensity profile of the GLM model is written as [90] where the shape of the radiations intensity profile I em depends on these free parameters γ, α and β: γ is related to the rate of increase of the intensity profile from infinity down to the peak; α adjusts the radial position translation of the whole intensity profile; β controls the dilation of the intensity profile as a whole.Hence, we can adjust the three free parameters to select adequate radiations intensity profiles for the models under study.In this paper, we mainly investigate that three kinds of innermost radiation position of the accretion disk: the innermost stable circular orbit r isco , the photon sphere radius r ph and the outside event horizon r + .The result for these three situations are plotted in Figs.11-16, and the more detailed discussions are as follows.
Case I: γ = −2, α = r isco and β = M/4.Firstly, one assume that the innermost stable circular orbit r isco as the radiation stop position, which is the the bounder between test particles circling the black hole and test particles falling into the black hole.It can be obtained by [44] which is derived in detail in the Appendix A. Using (4.4) and (4.5), we respectively depict I em (r) as the function of r in Fig. 11 and I obs (b) as the function of b in Fig. 14 for different values of Q, and the Fig. 14 shows the two-dimensional shadows cast of the black holes in the celestial coordinates as well.From the Fig. 14, one can see that the direct emission has a big range of b, the lensed ring is limited to a small range of b, and the photon ring is an extremely narrow rang of b.Combining with the Fig. 11 and Fig. 14, the emitted positions on thin disk, the total observed intensity and the shadow size of the black hole decrease with the increase of the Q in the charged Horndeski black hole, while these results decrease more slowly in the RN black hole when the Q increases.Case II: γ = −2, α = r ph and β = M/8.Then, one assume that the region of emission is well outside the photon sphere.The Fig. 12 and Fig. 15 describe that the I em (r) as a function of r and I obs (b) as a function of b with distinctive Q, respectively.One can see that the most important difference from the first model is that the positions of direct emissions are always located inside the lensed and photon ring intensities in Fig. 15.Similar to previous analysis, the Fig. 15 shows the regions of the direct emission are much more than the regions of the lensed ring and the photon ring.From the Fig. 12 and Fig. 15, we can observe that the emitted positions on thin disk, the total observed intensity and the shadow size of the charged Horndeski black hole are smaller as a larger Q.In the same way, they are more weakly affected in the RN black hole as a larger Q value.
Case III: γ = −3, α = r + and β = M/8.Finally, one assume that the position of emission is extended all the way down to the outside horizon.For the different Q, the I em (r) as a function of r and I obs (b) as a function of b are already depicted in Fig. 13 and Fig. 16, respectively.The lensing ring and photon ring are superimposed on the direct image as well.From the Fig. 16, the lensed ring is more prominent, but the direct emission remains dominant.The photon ring continues to be entirely negligible.Similarly, the the emitted positions on thin disk, the total observed intensity and the shadow size of the charged Horndeski black hole are smaller than those of the Schwarzschild black hole.In addition, the results are more significant in the charged Horndeski black hole in comparison with the RN black hole in Fig. 13 and Fig. 16.

V. CONCLUSION AND DISCUSSION
In this paper, we have investigated the shadow images and rings of the charged Horndeski black hole illuminated by spherical and thin disk-shaped accretion flows.We found that the larger electric charge results in a smaller peak effective potential at a larger charged Horndeski black hole radius, and the corresponding event horizon radius, shadow radius and critical impact parameter are gradually decrease with the increase of electric charge.In addition, we observed that the three functions are more sensitive to Q in the charged Horndeski black hole than those of the RN black hole with the increase of Q value.Then, we studied the observation characteristics of the charged Horndeski black hole shadows and rings on two toy models of spherical accretion.For the charged Horndeski black hole illuminated by the static spherical accretion flow, we found that the total observed intensity increases at a smaller b ph with the increase of the electric charge.It has been found that the two-dimensional shadows in celestial coordinates are not a totally dark shadow with zero intensity due to part of the radiation of the accretion flow behind the black hole can escape to infinity.The shadows size and rings brightness around the charged Horndeski black hole are smaller and brighter than those of the RN black hole as a larger Q value, respectively.Similarly, for the infalling spherical accretion flow, we found that there are have similar to results as well.It has shown that the brightness of rings in infalling accretion model are much less than those of static accretion model due to the Doppler effect, and the infalling accretion flow doesn't affect the shadow size of the back hole, implying that the black hole shadow size depends on the spacetime geometry and the ring luminosity depends on the accretion flow models.
Finally, we analyzed the shadow images and rings of the charged Horndeski black hole illuminated by thin disk-shaped accretion flows with different emission profiles.It has shown that the bright region near the charged Horndeski black hole can be divided into the direct emissions, lensed rings and photon rings according to the intersection times between the light line and the thin disk, and the contribution of the direct emission on the total observed intensity is dominating, the lensed ring is small, and the photon ring is hardly neglected.We found that the shadows size and the rings brightness around the charged Horndeski black hole are decrease with the increase of the electric charge, while the results are more obvious to Q in the charged Hondeski black hole than those of the RN black hole.It has been demonstrated that the radiation positions of the accretion flows affect the shadows size and the rings brightness around the charged Horndeski black hole.In addition, one can believe that the radius of the observed shadow is related to the boundary of the direct emissions.Therefore, we further investigate the relationship between the starting point of the direct emissions and the critical impact parameter, which may help us to utilize the EHT observations to measure the critical impact parameter, and then to test GR.

. 20 )FIG. 1 :
FIG. 1: Effective potential V ef f (r) as a function of black hole radius r for different electric charges with b = M = 1.

FIG. 2 :
FIG.2:The event horizon radius r+, the radius of the photon sphere r ph and the critical impact parameter b ph as the function of the electric charge Q with M = 1.

5 FIG. 3 :
FIG. 3: Trajectories of light rays around black hole of different Q in the polar coordinates (r, ϕ) with M = 1.

FIG. 4 : 8 FIG. 5 :
FIG. 4: The total observed intensities I s obs as a function of the impact parameter b with a static spherical accretion flow for M = 1.

FIG. 6 : 4 FIG. 7 :
FIG. 6: The observed specific intensities I i obs as a function of the impact parameter b with an infalling spherical accretion flow for M = 1.

5 FIG. 9 :
FIG. 9:The selection of associated photon trajectories for the black hole in the polar coordinates (r, ϕ).The black dashed line is the photon sphere, the black disk is a black hole.These red, blue and green line represent the light trajectories intersect with the thin disk once, twice and least three times, respectively.
FIG. 14: Optical appearance of the black holes for the case I with different Q.First column: the total observed intensities I obs (b) as a function of the impact parameter.Second column: the two-dimensional shadows cast of the black holes in the celestial coordinates.Third column: the zoomed in sectors.The black hole mass is set to 1.
FIG. 15: Optical appearance of the black holes for the case II with different Q.First column: the total observed intensities I obs (b) as a function of the impact parameter.Second column: the two-dimensional shadows cast of the black holes in the celestial coordinates.Third column: the zoomed in sectors.The black hole mass is set to 1.
FIG. 16: Optical appearance of the black holes for the case III with different Q.First column: the total observed intensities I obs (b) as a function of the impact parameter.Second column: the two-dimensional shadows cast of the black holes in the celestial coordinates.Third column: the zoomed in sectors.The black hole mass is set to 1.