Influence of inhomogeneities on holographic mutual information and butterfly effect

We study the effect of inhomogeneity, which is induced by the graviton mass in massive gravity, on the mutual information and the chaotic behavior of a 2+1-dimensional field theory from the gauge/gravity duality. When the system is near-homogeneous, the mutual information increases as the graviton mass grows. However, when the system is far from homogeneity, the mutual information decreases as the graviton mass increases. By adding the perturbations of energy into the system, we investigate the dynamical mutual information in the shock wave geometry. We find that the greater perturbations disrupt the mutual information more rapidly, which resembles the butterfly effect in chaos theory. Besides, the greater inhomogeneity reduces the dynamical mutual information more quickly just as in the static case.


Introduction
Over the past two decades, there have been a large overlap studies among high energy theory, condensed matter physics, quantum information and geometry based on the conjecture of gauge/gravity correspondence [1][2][3]. The typical example is the holographic entanglement entropy (HEE) proposed by Ryu and Takayanagi that relates the minimized codimension-2 bulk surface to the quantum information in the boundary field theory [4]. The Ryu-Takayanagi conjecture greatly simplifies the computations of the entanglement entropy in the field theory, and numerous papers investigates the relation between the bulk geometry and the quantum information in the boundary field theory. For recent reviews, please refer to Refs. [5,6]. Among these studies, one interesting thing is to adopt the HEE to probe the thermalization process or phase transitions in the strongly coupled systems [7][8][9][10][11][12][13][14][15][16][17].
Mutual information measures the correlations between two subsystems in quantum information theory [18]. The holographic setup of mutual information can be obtained by calculating the length of the wormhole which connects the two sides of an eternal AdS black hole [19][20][21]. In [20], the authors defined the thermofield double states (TFD) to describe the states of the entangled two sides of the black hole, i.e., where β is the inverse of the temperature while |i L and |i R are the identical quantum states on the two-sided AdS black holes. Supposing there are two identical space-like subregions where S(A), S(B) are the entanglement entropy of A and B respectively, while S(A ∪ B) is the entanglement entropy of the union of A and B. Holographically, S(A) and S(B) can be calculated by the areas of the minimal surfaces in the bulk geometry associated to A and B independently from the Ryu-Takayanagi conjecture, while S(A ∪ B) can be computed from the minimal surface which crosses the horizon and stretches through the wormhole connecting the two subregions A and B [22]. The disruption of mutual information was related to the butterfly effect by Shenker and Stanford [23]. Specifically, as a small perturbation, such as a light-like energy perturbation, is added to one side of the eternal black hole, the mutual information between the two sides may be disrupted after an amount of time t * , which means there is no dependence or entanglement between the two sides of the black hole. The time t * , usually called the scrambling time in black hole systems, is proportional to the logarithm of the entropy of the black hole, e.g., t * ∼ β log(S)/2π, in which S is the Bekenstein-Hawking entropy. The disruption of the mutual information is a piece of evidence of the system's high sensitivity to the initial conditions, which reminds us of the terminology in chaos theory-butterfly effect. Recent literatures relevant to the butterfly effect can be found in [24][25][26][27][28][29].
In this paper, we intend to study the static and dynamical mutual information in the background of massive gravity [30,31]. As we know from [32][33][34] that the graviton mass in the bulk breaks the diffeomorphism invariance, which makes the stress-energy tensor non-conserved in the dual field theory. The non-conservation of the stress-energy tensor causes the momentum dissipation in the boundary field theory. Therefore, from this aspect the graviton mass plays the role of inhomogeneity in the boundary field theory, i.e., greater graviton mass is dual to greater inhomogeneity in the boundary 1 . There have been a number of papers investigating the inhomogeneous effects caused by massive gravity so far, see for example [35][36][37][38][39][40][41][42][43].
Firstly, we study the holographic mutual information of two identical strips in a static background of a 3+1-dimensional massive gravity. For the strips with larger length, the mutual information between them decreases monotonically as the graviton mass increases; However, the mutual information between two shorter strips first increases and then decreases with respect to the graviton mass. A plausible reason is that the graviton mass (or equivalently spatial inhomogeneity) would have greater effects on the strips with larger length. When the system is near-homogeneous (or equivalently with small graviton mass), the mutual information for the strips with shorter length will increase with respect to the temperature of the black hole (referring to the right panel of Fig.4 and discussions in the main context), which is similar to the relation between mutual information and temperature in the pure homogeneous case [24]. However, as the graviton mass grows big enough, i.e., far from homogeneity, the mutual information will instead decrease with respect to the graviton mass, which behaves distinctly from the homogenous case. Therefore, we argue that the greater graviton mass (or stronger inhomogeneous effects) would play a dominant role in decreasing the mutual information compared with the temperature. On the other hand, if the strips are longer, the mutual information decreases monotonically according to the graviton mass regardless of the system being near or far from homogeneity. It intuitionaly suggests that the spatial inhomogeneity would have greater impacts to the strips with larger length by destroying the mutual information monotonically. A detailed discussion of the relation between mutual information and the graviton mass in the static background case is given in Section 3.
In all the parameter regimes we considered in this paper, we find that the critical width of the strip, which renders the mutual information vanishing, always decreases according to the graviton mass. Moreover, as the width of the boundary strip increases, it is found that the critical charge which disrupts the mutual information increases as well. As we know, in the massive gravity the temperature of the black hole will decrease as the charge increases. Therefore, we can infer that the critical width of the strip increases when the temperature of the black decreases, which is consistent with the results obtained in [24].
In order to study the dynamical behavior of mutual information in the boundary field theory, one of the approaches is to add extra energy perturbations into the bulk, which will lead to a shift on the horizon of the black hole. This shift will affect the mutual information on the two sides of the black hole. The dynamics of the bulk after adding the perturbations can be conveniently investigated in the shock wave geometry with Kruskal coordinates [44]. We find that as the shift grows (more added perturbed energy at the initial time), the mutual information will be reduced more significantly, which reminds us of the phenomenon in chaos theory -butterfly effect. By turning on the graviton mass, we also find that the greater the mass is, the smaller the mutual information will be, which indicates that the spatial inhomogeneity will reduce the mutual information just like in the static case mentioned above. This paper is organized as follows: In Section 2 we briefly review the dynamical holographic mutual information in the shock wave geometry. We investigate the static mutual information in the massive gravity background in Section 3. In Section 4 we study the dynamical mutual information by adding energy perturbations into the bulk of the massive gravity. Finally we draw our conclusions and discussions in Section 5. Through out this paper we use natural units (G = c = = 1) for simplicity.

Reviews of shock wave geometry and holographic mutual information
The butterfly effect of a black hole is usually studied in the shock wave geometry, therefore, we are going to briefly introduce the key ingredients of the shock wave geometry at first. For simplicity, we will adopt a planar symmetric black hole in 3 + 1 dimensions with the line element, in which (x, y) are the transverse directions in the bulk while r represents the radial direction. The Hawking temperature of this black hole is T = κ/2π, in which κ = f (r)| r h /2 is the surface gravity with r h the horizon radius. From the AdS/CFT dictionary, the temperature T can be regarded as the temperature of the dual field theory.

Shock wave geometry
The shock wave geometry is conveniently discussed in the Kruskal coordinates [44]. The metric in Eq.(2.1) can be rewritten as where U = t − r and V = t + r are the Eddington-Finkelstein coordinates, which can be defined by the tortoise coordinate r = dr/f (r). Therefore, as r approaches the event horizon and boundary, r tends to −∞ and 0 respectively. The Penrose diagram of the geometry (2.2) is shown on the left panel of Fig.1, where a dot represents a two dimensional space in (x, y) directions. We suppose µ > 0, ν < 0 in the left exterior region while µ < 0, ν > 0 in the right region as in [24]. Thus from Eq.(2.4) we know that the event horizon and boundary are located at µν = 0 and µν = −1, respectively. The light is going along µ = constant and ν = constant. The shock wave geometry can be obtained by adding a small perturbation of energy E into one side, for instance left side of the black hole [44]. Supposing that at the boundary time t w we add a light-like perturbation of the energy, which goes along a constant µ trajectory, into the left boundary. We label the Kruskal coordinates on the left side and right side as µ L , ν L and µ R , ν R , respectively. The constant µ trajectory of the perturbation implies (2.5) In order to find the relation between ν L and ν R , we employ the relation Generally speaking, κ L = κ R = κ for the energy E of the perturbation is much smaller than the mass of the black hole M . On the other hand, we are interested in the case t w → ∞, which implies r → r h . In this case, we can approximate r ≈ 1 2κ (log(r − r h ) + c) with c a constant of integration. Hence, e 2κr = C(r − r h ), where C = e c . Therefore, we have the identification in which the relation C L = C R = C has been used. It should be stressed that even (r hL − r hR ) → 0, the formula e κtw (r hL − r hR ) is still finite in Eq. (2.7). The difference h between ν L and ν R is the shift close to the horizon. The Penrose diagram of the shock wave geometry is shown on the right panel of Fig.1.
To get the standard shock wave geometry as showed in [23], one often employs the replacement ν → ν + h(θ)Θ(µ), where Θ(µ) is a step function. The shock wave geometry is strictly a solution to the Einstein equation is the boost energy arising from the null perturbation at the initial time.
The scrambling time t is defined as the value of t w when h ∼ O(1) for in this case the mutual information vanishes [23]. On the basis of Eq. (2.7) and the first law of the black hole thermodynamics, the scrambling time can be written as t ∼ β 2π log[c(r h )S], in which c(r h ) is a function of the black hole horizon. In various gravity models, this form is universal and the only difference is embodied in the function c(r h ) [45].

Holographic mutual information
As depicted on the left panel of Fig.1, an eternal black hole has two asymptotic AdS regions, which can be holographically described by the so-called TFD states of the two identical, non-interacting conformal field theories [20]. Our objective is to compute the mutual information of a subregion A on the left asymptotic boundary and its partner B on the right asymptotic boundary. For simplicity, we will let A = B so that the left and right boundaries are identical.
We are interested in the 3+1-dimensional planar black holes, thus the AdS boundary has a 2-dimensional space parameterized by coordinates (x, y). We will consider the subregion A or B as a strip, which has the width x ∈ (0, x 0 ) and extends along the y direction with length Y 2 . Therefore, the entanglement entropy S A of the subregion A is where Area A is the area of the minimal surface in the bulk, i.e., where r = dr/dx. If regarding the integrand in Eq. (2.8) as a 'Lagrangian' L, one can define a conserved quantity associated to translation in x-direction, that is where r min is the turning point of the surface with r = 0. According to the symmetry of the surface, the turning point locates at x = x 0 /2. With the help of the conserved equation (2.9), x 0 can be written as and the minimal area in Eq.(2.8) can be rewritten as Since B is identical with A, Area B thus takes the same form as Area A . As stressed in the introduction, we will employ the mutual information, defined by I( to study the correlation between regions A and B. Therefore, our next step is to find S A∪B or the area Area A∪B , which is the minimal surface connecting regions A (left) and B (right) by passing through the horizon of the black hole. From [24], the total area Area A∪B including both sides of the horizons can be expressed as Combining all the Eqs.(1.2), (2.11) and (2.12) together, we have We are interested in how the width of the strip x 0 affects the mutual information. Substituting Eq. (2.10) into Eq. (2.13), we obtain 14) It is of great interest to find the critical value of the width x 0c where the mutual information vanishes, i.e., I(x 0c ) = 0, which leads to

Holographic mutual information in the static background case
Massive gravity is a deformation of the Einstein gravity with graviton mass [30,31]. Diffeomorphism invariance is broken in massive gravity because of the graviton mass. Therefore, the stress energy tensor is not conserved any more in the dual field theory. The nonconservation of the stress energy tensor corresponds to a momentum dissipation on the dual boundary field theory [34]. Therefore, from this sense the graviton mass plays the role of inhomogeneity on the boundary field theory. It is of great interest to investigate the effect of inhomogeneity on the holographic mutual information, in particular, we are going to explore the effects of the graviton mass on the mutual information with and without energy perturbations in the following context. The action of the massive gravity is as follows [32,35] where m is the graviton mass parameter, F µν = ∂ µ A ν − ∂ ν A µ is the field strength, f µν is the reference metric, c i are constants, and U i are symmetric polynomials of the eigenvalues of the 4 × 4 matrix K µ ν ≡ √ g µα f αν : The square root in K stands for ( In 3+1dimensions, a black hole solution with line element in Eq.(2.1) can have the gravitational potential as in which M and Q are the mass and charge of the black hole respectively; c 0 , c 1 and c 2 are constant parameters associated to the graviton mass. In this paper, we will fix the parameters c 0 = c 1 = 1, c 2 = −1/2 in order to render the background thermodynamically stable [35,39]. The temperature of the black hole is readily obtained as 3.1 Effect of the width of the strip on the mutual information Firstly, we are going to study the relation between the mutual information and the width of the strip in the background of massive gravity. One can readily read off the relation between the width of the strip and the position of the turning point r min from Eq.(2.10). The relation is shown on the left panel of Fig.2. From Eq.(2.10) one finds that as r min → r h , the integral for the width x 0 diverges, which can be seen as well from the left plot of Fig.2. Intuitively the left panel of Fig.2 is also correct, since as we know greater r min is closer to the infinite boundary, therefore, it is obvious that the greater r min corresponds to the smaller width of the strip. From Eq.(2.13) or (2.14), one can see that as r min → r h the mutual information diverges, since in this case the widths of the strips on the two boundaries are nearly divergent (cf. the left panel of Fig.2). From [24] we know that the mutual information for divergent strips will be divergent too. This phenomenon can also be found on the right panel of Fig.2. In addition, from this subgraph we observe that the mutual information vanishes where r min ≈ 1.25, which indicates that there is a critical value for the position of the turning point (or equivalently a critical width of the strip x 0c from the left panel of Fig.2) rendering the mutual information vanishing. Combining the two panels in Fig.2, we plot the relation between the mutual information and the width of the strip x 0 explicitly in Fig.3, from which we can clearly see that the critical width of the strip is roughly x 0c ≈ 1.13. In addition, we find that the mutual information always grows as the width of the strip increases, which is easy to understand for in this case the subsystems on the two asymptotic boundaries are larger and their entanglements becomes greater as well.

Effects of the graviton mass and black hole charge on the mutual information
The effects of graviton mass m and charge Q of the black hole on the mutual information I(x 0 ) are shown in Fig.4. From the left panel of Fig.4 we see that for each curve, the mutual information decreases as the charge increases. Besides, there is a critical charge Q c that makes the mutual information vanishing, which means that there is no entanglement between the paired subregions we considered. For a small fixed charge, we find that the mutual information is smaller for greater r min . As we know from the preceding subsection, bigger r min actually corresponds to smaller width of the strip on the boundary. Therefore, the left panel of Fig.4 also indicates that smaller subregions have smaller mutual information between them, which is consistent with the preceding subsection. From Ref. [24] we learn that when the temperature of the black hole grows, the critical width of the strip decreases. It also means that for a fixed width of the strip, as temperature grows the mutual information of the two strips will grow as well. In the case of the massive gravity, the temperature will decrease as the charge Q increases when fixing other parameters of the background, please refer to Eq.(3.4). Therefore, from the left panel of Fig.4 we see that as charge grows (temperature decreases) the mutual information decreases monotonically, which matches the conclusions mentioned before. Incidentally, we also checked other cases of the graviton mass, similar results were obtained.
However, it is interesting to see that the mutual information does not have monotonic decreasing behavior to the graviton mass on the right panel of Fig.4. In particular, when r min is bigger (x 0 is smaller) the mutual information first increases to a maximum value and then decreases with respect to the graviton mass. Let's call the critical graviton mass as m c which corresponds to the maximum value of mutual information. From Eq.(3.4) we see that as m grows the temperature of the black hole grows as well when fixing other parameters. Therefore, from the conclusions in the preceding subsection it seems that as m grows the mutual information should increase as well, since the mutual information increases with respect to the temperature. However, as we see from the right panel of Fig.4, the mutual information decreases as m grows in the regime m > m c , which contradicts the above conclusions.
The possible way coming to rescue is that the conclusion in the previous subsection is mainly valid in the (near-)homogeneous case. However, when taking into account of the inhomogeneous effects induced by the graviton mass, the conclusions should be different. Therefore, we can infer from the right panel of Fig.4 that when m < m c the inhomogeneous effects are tiny. Hence, in this case the mutual information still grows as the temperature grows (i.e., m increases); However, when m > m c the inhomogeneous effects are significant, which will disrupt the mutual information and finally render it vanishing as m is big enough. Therefore, we can infer that greater inhomogeneity will spoil the mutual information.
Finally, let's come to the top curve on the right panel which shows that the mutual information decreases monotonically with respect to m. In fact the top curve corresponds to a smaller r min (bigger width x 0 of the strip). Therefore, it makes us speculate that the inhomogeneity will have greater effects on longer strips than shorter ones to reduce the mutual information. Therefore, for a long strip the mutual information only decreases according to m since now the inhomogeneous effects dominate. The critical width x 0c of the strip is the width that renders the mutual information vanishing, i.e., I(x 0c ) = 0. On the left panel of Fig.5, by fixing the charge Q = 2, we see that the critical width x 0c decreases as the graviton mass m grows for various black hole horizons r h . This means that the greater inhomogeneity will reduce the critical width x 0c . For a fixed value of m, the critical width x 0c decreases with respect to the increasing horizon r h . We know that the greater horizon corresponds to higher temperature of the black hole, therefore, the left plot of Fig.5 also indicates that the higher temperature also reduces the critical width x 0c .

The critical width
On the right panel of Fig.5, we plot the relation between the critical width x 0c and the charge Q for various horizons r h . For a fixed horizon, we find that when Q increases the critical width x 0c increases as well. In particular, when r h is small the critical width increases more rapidly; However, when r h is big, x 0c increases more mildly. As we know from Eq.(3.4) that as Q increases the temperature T M G decreases. Therefore, the right plot of Fig.5 indicates that as temperature of the black hole decreases the critical width x 0c will increase, which is consistent with the previous analysis. Moreover, for a fixed charge Q, the critical width x 0c decreases as r h increases, which states that the higher temperature reduces the critical width x 0c . This statement matches the conclusions above.

Holographic mutual information in the dynamical background case
As we know from Section 2, when a small perturbation is added from the left boundary to the bulk, there will be a shift h(x) near the horizon in the ν direction for a long enough time t w . A shock wave geometry thus forms and the wormhole connecting the left and right regions may be destroyed in some circumstances. Therefore, the mutual information between the two subregions will be disrupted. In the previous section we studied the holographic mutual information in the static background of massive gravity; In this section, we are going to investigate the dynamical behavior of the holographic mutual information in the shock wave background of massive gravity. As in Section 3, we suppose that the strip A sits in the left asymptotic boundary and its identical partner B in the right boundary. After adding the light-like perturbations, the areas of minimal surfaces Area A and Area B are unaffected by the shock wave since they do not cross the horizon, while the Area A∪B is affected by the shock wave since it passes across the horizon and stretches through the wormhole. The sketchy plot of the surface Area A∪B is shown in Fig.6. Our main goal thus is to calculate the Area A∪B and to study how it changes with respect to the shift h(x) for a fixed boundary separation.
Because of the symmetry of the minimal surface, we should only calculate the areas for the regions 1, 2 and 3 in Fig.6 [24,26]. At a surface with constant x, the induced metric can be written as in whichṙ = dr/dt. The area of minimal surface for the regions 1, 2 and 3 in Fig.6 is then given by If regarding the integrand in Eq. (4.2) as a 'Lagrangian' L, we can define a conserved 'Hamiltonian' H as in which f 0 = f (r 0 ) and r 0 is the radial position behind the horizon satisfyingṙ = 0. From In contrast, the coordinate ν 2 can be determined by choosing a reference surface r =r for which r = 0 in the black hole interior. Thus, we reach The third segment stretches from (µ 2 , ν 2 ) to (µ 3 = 0, ν 3 = h/2). Therefore, with the relation we can express h as where 14) It is obvious that h depends on the location of r 0 for a fixed r h . Combining Eq. (4.8) with Eq. (4.13), one can get the relation between I(h, x 0 ) and h. Next, we will study the relation between I(h, x 0 ) and h numerically in the background of the shock wave geometry in the massive gravity.

Black hole charge and butterfly effect
We are going to study how the charge affects the shock wave geometry and the dynamical behavior of mutual information in the following context. As we already knew from the previous subsection that as r 0 → r h (we set r h = 1 in numerics), the shift h vanishes, which can also be found in the plot (a) of Fig.7. The vanishing of h means there is no added perturbation into the bulk, therefore, the shock wave geometry goes back to the unperturbed or the static case which we studied in the previous section. From the formula of Ξ 3 in Eq.(4.16), we can see that if the denominator in the integrand vanishes, i.e., 1 + H −2 f r 2 = 0, the shift h diverges. Since H is a conserved quantity, the above equation can be readily transformed to dr r=r crit = f (r crit )r crit + 2f (r crit ) = 0,  vanishing. In fact, one can see from plot (a) of Fig.7 that for a fixed r 0 , greater values of the charge correspond to smaller values of shift h. As we know that the shift h is proportional to the energy of the added perturbation, thus we can deduce from panel (b) of Fig.7 that the more the added energy perturbation is, the smaller the mutual information will be. The relationship between mutual information and the shift h can be more clearly seen from the panel (c) of Fig.7. As the shift grows to a critical value h crit , the mutual information will vanish, i.e., when the added perturbation is big enough, it will finally disrupt the entanglement between the two strips. The disruption of the mutual information depending on the added perturbation in the initial time reminds us of the phenomenon in chaos theory, i.e., butterfly effect. For a fixed value of shift, the mutual information decreases as the charge grows, which is consistent with the statements in the previous section. In particular, when the shift is zero (or equivalently r 0 → 1), i.e., there is no added perturbation into the bulk, it will go back to the static case discussed in the previous section.

Graviton mass and butterfly effect
Plot (a) of Fig.8 shows the relation between the shift h and r 0 for different graviton masses. Just like the case in the preceding subsection, as r 0 → r h = 1, h vanishes; Besides, there is also a critical value of r 0 that makes the shift h diverge. The critical values of r 0 match the ones from Eq.(4.17). For a fixed r 0 , larger m corresponds to larger h according to the panel (a) of Fig.8; Meanwhile we find from the panel (b) that larger m corresponds to smaller mutual information, which indicates that larger inhomogeneity will disrupt the mutual information. This is consistent with the statements in the static case in the previous section.
The relationship between the mutual information and the shift h is shown in the plot (c) of Fig.8, which can be obtained from the rest two plots (a) and (b). Just like in the plot (c) of Fig.7, the mutual information decreases according to the shift h, which means the added perturbation will destroy the mutual information between the two sides of the black hole. And for a fixed value of the shift, the mutual information decreases with respect to the graviton mass as we already discussed above. There are also critical shifts h c that render the mutual information vanishing, and it is found that larger m corresponds to smaller h c .
We have also checked other values of m for these relations and found that they have similar behaviors as in Fig.8.

Conclusions and Discussions
In this paper we studied the holographic mutual information in the background of massive gravity. In the static case, we found that for shorter strips the near-homogeneous mutual information would increase as the temperature of the black hole grows, which is consistent with the conclusions in the pure homogeneous case studied previously. However, for a larger graviton mass, which corresponds to greater inhomogeneity in the boundary field theory, we found that the mutual information would decrease with respect to the graviton mass. For longer strips the mutual information would decrease monotonically with respect to the graviton mass, which implies that the spatial inhomogeneity has larger influence to the mutual information of strips with larger length. With the above results, we argued that when the system is far from the homogeneity, the spatial inhomogeneity plays a dominant role in affecting the mutual information than the temperature of the black hole. We also investigated the effect of the charge on the mutual information and found that the mutual information decreases as the charge increases, which is independent of the width of the strip.
By adding the light-like perturbations into the bulk, we studied the dynamical mutual information in the shock wave geometry of the massive gravity. The added perturbations produce a shift on the horizon in the Kruskal coordinates. From the existing studies of the shock wave geometry, we know that the shift is proportional to the added energy. We found that the more the added energy was, the smaller the mutual information would be, which suggests that the added perturbations would reduce the mutual information between the two sides of the black hole. We also investigated the effect of the charge and graviton mass on the critical values of the shift, where makes the dynamical mutual information vanish. We found that the larger the values of the charge and graviton mass are, the smaller the critical values of the shift are, which indicates that both the charge and inhomogeneity would reduce the mutual information.
It would be interesting to study the holographic mutual information in the background of real spatial inhomogeneity, such as adding lattice structures in the bulk [46] or for simplicity with the spatially dependent sources on the boundary [47]. The advantage of the latter is that the bulk spacetime is homogeneous and isotropic; Moreover, turning on the sources of the massless scalar fields on the boundary would provide more physical meanings of the momentum relaxation in the boundary field theory. It was shown in [47] that the parameters related to the scalar sources are similar to the graviton mass in some sectors of the massive gravity. Therefore, it would be very interesting to check whether the parameters of the scalar sources really play the same role as the graviton mass to the holographic mutual information in the massive gravity. We will leave this as our future study.