Quasinormal oscillations and late-time tail of massless scalar perturbations of a magnetized black hole in Rastall gravity

  • In this study, we investigate the quasinormal mode and late-time tail of charged massless scalar perturbations of a black hole in generalized Rastall gravity. The black hole metric in question is spherically symmetric, accompanied by a power-Maxwell field surrounded by a quintessence fluid. We show that the massless scalar field, when dressed up with the magnetic field, acquires an effective mass, which significantly affects the properties of the resultant quasinormal oscillations and late-time tails. Specifically, the quasinormal frequencies become distorted and might even be unstable for particular spacetime configurations. Additionally, the exponent of the usual power-law tail is modified according to the modification in the structure of the branch cut of the retarded Green's function. In particular, as the effective mass is generated dynamically owing to the presence of the magnetic field, we may consider a process through which the field is gradually removed from the spacetime configuration. In this context, while the quasinormal oscillations converge to the case of massless perturbations, we argue that the properties of resultant late-time tails do not fall back to their massless counterpart. The relevant characteristics are investigated using numerical and analytic approaches.
  • The black hole is one of the most intriguing predictions of Einstein's general relativity. It is a particular spacetime region typically created by the gravitational collapse of a compact object in its late stage, when its mass exceeds the upper limit of neutron stars [1]. From both theoretical and experimental perspectives, the subject has gained much interest over the past few decades. As the culmination of continuous efforts, the recent advent of empirical detection of gravitational waves emanating from black holes has elicited further interest. In particular, the merger signals of a black hole binary captured by the LIGO and Virgo collaboration [2] furnished the first direct piece of evidence of black holes and subsequently inaugurated a novel era of gravitational-wave astronomy.

    The quasinormal modes (QNMs) carry the characteristic sound of a black hole. In terms of dissipative oscillations, they bear the essential properties of the underlying spacetime metric [3]. The associated temporal profile constitutes the ringdown phase of the collapse, from which a black hole is formed. On the one hand, such signals are mathematically straightforward compared with those emanating from the merger process. On the other hand, the strength of the associated gravitational radiation is speculated to still fall within the scope of the ongoing space-borne gravitational wave detector programs, such as LISA and TianQin [4, 5]. In particular, studies have been performed to evaluate the detector signal-to-noise ratio regarding the feasibility of an eventual detection [6, 7].

    Theoretically, QNMs can be investigated using black hole perturbation theory. Among the earliest endeavors, Regge and Wheeler created the framework of the metric perturbations in spherically symmetric black holes [8]. Subsequent developments by Zerill [9] and Moncrief [10] further refined a few pertinent aspects. For rotating black holes, Teukolsky first derived the linearized scalar and metric perturbation equations for the Kerr spacetime [11, 12]. Although the physical system is dissipative and therefore non-Hermitian, the mathematical procedure to obtain QNMs can be effectively considered an eigenvalue problem. Subsequently, the emergence of discrete complex frequencies is attributed to the ingoing and outgoing boundary conditions, respectively, applied at the horizon and outer spatial bound (often located at infinity). By using the continued fraction method proposed by Leaver [13, 14], we encounter the discrete complex frequencies reminiscent of the atomic spectrum of hydrogen in quantum mechanics. The eigenvalue nature of the problem is also transparently considered in terms of the matrix method [1518], in which we explicitly solve a matrix secular equation for the quasinormal frequencies. In contrast, Leaver proposed the Green's function approach, which reformulates the general solution in terms of the spectrum decomposition [3, 19]. From this perspective, the main characteristics of the perturbation fields are attributed to the intrinsic singularities of the relevant Green's function. Specifically, the QNMs correspond to the poles of the Green's function, whereas the branch cuts essentially govern the late-time tails. In addition to the QNMs, the late-time tail is a topic of pertinent relevance. On the one hand, the time profile of QNMs decays exponentially; on the other hand, the late-time tail largely follows a power-law form. Therefore, the final stage of the temporal evolution of the perturbations is dominated primarily by the late-time tail. As mentioned earlier, mathematically, the origin of the late-time tail is attributed to the branch cut and cannot be obtained by summing the QNM poles. For instance, for massless scalar perturbations, the tail is due to a branch cut on the negative imaginary axis in the frequency domain, which stretches from the origin. Researchers have argued that the presence of the branch cut sensitively depends on the asymptotical structure of the potential [20]. Additionally, it is often interpreted in the literature as the backscattering of perturbed wave packets by the spacetime far away from the horizon [2123]. Price was the first to observe an inverse power-law tail, t(2l+3) or t(2l+2), for massless scalar perturbations in the Schwarzschild black hole metric [24]. Extensive studies for different metrics further indicated that such an inverse power-law form is a rather general characteristic for massless perturbations in spherical spacetimes that are asymptotically flat [21, 25, 26]. Moreover, the temporal profile for massive perturbations might be qualitatively different from their massless counterparts. For insignificant mass, Hod and Piran discovered that the time-domain profile at the intermediate-late time in Reissner-Nordström spacetime has an oscillatory tail with a decay rate of t(l+3/2)sin(μt) [27], where μ is the mass of the field. Moreover, Koyama and Tomimatsu analytically showed that the asymptotical late-time tail in the Schwarzschild spacetime is t5/6sin(μt) [28, 29]. For some specific cases, the asymptotical form was recenetly reported to also be t1sin(μt) [30] in the limit of vanishing scalar mass. In all the massive cases, the branch cut can be conveniently selected to be on the real axis in a finite interval determined by the scalar mass. It is also worth noting that such a decay rate does not depend on the angular momentum l. In addition to the power-law form, the exponential tails have also been observed in asymptotically de Sitter spacetimes for various types of fields [3133]. Notably, the Dirac perturbations eventually decay to a non-zero constant [23, 34]. These distinctive characteristics have significantly enriched the physical content regarding the late-time waveforms. More recently, Cardoso et al. proposed [35] that the late-time ringdown may serve to probe the horizon and discriminate between different gravitational systems, inclusively, the exotic compact objects such as gravastar and wormhole. Based on this, black hole echoes [36], and particularly, their connections with the QNMs and late-time tail [37, 38] have also attracted much interest.

    The magnetic field is another indispensable assembly piece of the universe. Experimentally, employing the rotation measures, information on the magnetic field in the disk of galaxies can be extracted for a large number of pulsars [39]. In particular, strong magnetic fields, as high as 104108 G, have been observed in the vicinity of stellar-mass or supermassive black holes [40]. Therefore, relevant studies on black hole metrics, including their perturbations, have to be placed in the context of a realistic astrophysical environment. A stationary axisymmetric black hole solution with a uniform magnetic field placed along the symmetry axis was first derived by Wald [41]. The solution for the Kerr-Newman black hole was subsequently obtained by Ernst and Wild [42]. Among many important developments, the magnetic field was observed to have an essential role in the stability of the dynamic system. The latter encompass the superradiance instability [43] and black hole QNMs [44, 45]. In particular, the Zeeman effect was recently uncovered in the QNM spectrum of charged scalar perturbations in a magnetized Schwarzschild black hole [44]. Moreover, the instability of the metric can be triggered by tuning the effective scalar mass, which is modified by the interaction between the charge and magnetic field.

    As an alternative theory of general relativity, Rastall gravity was formulated based on the possible breaking of the conservation law of the stress energy-momentum tensor [46]. The theory aims at addressing several unsettled problems in general relativity regarding the largest scale when compared with cosmological observation. In the literature, various intriguing aspects of the theory have been explored, which include black hole physics [4750] and cosmology [5155]. Among others, the theory inherently results in a particle creation mechanism, which naturally supplies an alternative implementation for dark energy [56]. It is meaningful to further investigate the impact of the magnetic field in Rastall theory on spacetime stability.

    This paper is motivated by the above considerations. We investigate the QNMs and late-time tails of charged massless scalar perturbations of a magnetized black hole solution derived from generalized Rastall gravity. First, we generalize the black hole solution in Rastall gravity into a more realistic context by including the magnetic field. We then analyze the effect of the magnetic field on the spacetime stability in the context of Rastall gravity by exploring the quasinormal oscillations. For this, the study of massless scalar perturbations is sufficient for the purpose. For stable spacetime configurations, we know that massless and massive scalar fields behave distinctively in the late-time evolution. Therefore, it would be rather interesting to investigate the implications for Rastall's theory with the presence of a magnetic field. In particular, we observe that the exponent of the power-law tail is modified. This can be attributed to the fact that the effective mass of the field is "dressed up" owing to the interaction between the magnetic field and charge.

    The remainder of the paper is organized as follows. In the Sec. II, we present the metric of a magnetized black hole investigated in this study. The QNMs are explored in Sec. III by employing the Wentzel-Kramers-Brillouin (WKB) approximation, matrix method, and finite difference method. In Sec. IV, we analyze the late-time tail of the massless scalar perturbations and the effect of the magnetic field. Both the Green's function and the finite difference method are utilized. Further discussions and concluding remarks are provided in the final section.

    According to Rastall's recipe [46], the breaking of the conservation law of the stress energy-momentum tensor is

    Tvμ;v=λR,μ ,

    (1)

    where λ is the Rastall coupling parameter, and R is the Ricci scalar. If λ=0, the stress energy-momentum tensor is conserved and Rastall gravity falls back to general relativity. Moreover, the corresponding field equation can be expressed as

    Rμν+(κλ12)gμνR=κTμν,

    (2)

    where κ is a constant related to the trace anomaly of the energy momentum tensor in Rastall gravity.

    In this paper, we consider a spherically symmetric black hole solution surrounded by a quintessence fluid and power-Maxwell field first derived in [57]. Specifically, the metric of a static, spherically symmetric four-dimensional spacetime is expressed as

    ds2=f(r)dt2+f(r)1dr2+r2dθ2+r2sin2θdϕ2.

    (3)

    The action for the electromagnetic sector is of a nonlinear form, which can be expressed as [58]

    LF=(ξF)s,

    (4)

    where F=FμνFμv, Fμν is the Faraday tensor, and s and ξ are constants. s=1 corresponds to the case of a linear electromagnetic field. Owing to the presence of the electromagnetic field and quintessence fluid, the resulting energy momentum tensor has the form

    Tμv=Eμv+Tμv,

    (5)

    where the electromagnetic part, Eμv, is given by

    Evμ=(ξ)s(F)s1(2sFσμFσv12δvμF).

    (6)

    Tμv is the energy-momentum tensor of the quintessence field. If we assume a barotropic equation of state p=wρ, it can be expressed as

    Ttt=Trr=ρ(r),Tθθ=Tφφ=12ρ(r)(3w+1).

    (7)

    By substituting the energy momentum tensor into the field equation Eq. (2), we can show that the black hole metric is [57]

    f(r)=12Mr+Q2r212s(1+2s)32s12s(32s)s+Car16κλ+w(36κλ)1+3(1+w)κλ×(13(1+w)κλ)23(1+4κλ)(κλ+w(1+κλ)),

    (8)

    where Ca is a constant of integration, and M is related to the mass of the black hole.

    For this study, we will focus on the case of the linear electromagnetic field with s=1 and the sole presence of a constant magnetic field, namely, the electric charge vanishes, Q=0. When black hole is in the magnetic field, as a good approximation, the vector potential of the electromagnetic field can be expressed as

    Aμ=12Br2sin2θ(0,0,0,1).

    (9)

    The temporal evolution of a massless charged scalar field is governed by the Klein-Gordon equation, as follows:

    gμν(μiqAμ)(νiqAν)Φ=0.

    (10)

    To simplify the above equation, we introduce some reasonable assumptions [45]. First, we obtain αAα=0 by selecting the Lorentz gauge. Additionally, we assume that q2B21 by considering the scenario in which the coupling between the scalar and electromagnetic fields is sufficiently weak. Under these conditions, Eq. (10) can be rewritten as

    1r2r(r2fΦr)(L2r2qBLz)Φ1f2Φt2=0,

    (11)

    where

    L2=[1sinθθ(sinθθ)+1sin2θ2ϕ2]

    (12)

    and

    Lz=iϕ.

    (13)

    Now, we may employ the method of separation of variables by expressing the wave function as

    Φ(t,r,θ,ϕ)=12πdωeiωtlRlm(r,ω)rYlm(θ,ϕ),

    (14)

    where Ym(θ,ϕ) are the spherical harmonics, and the radial part of the wave function Rm(r,ω) satisfies

    (d2dr2+ω2Veff(r))Rlm(r,ω)=0,

    (15)

    where the effective potential is given by

    Veff(r)=f(r)[l(l+1)r2+f(r)rmqB],

    (16)

    and r is the tortoise coordinate defined by dr=drf(r). We note that Veff=mqB as r.

    In the following, we evaluate the QNMs using the sixth order WKB approximation [5964] and matrix method [1518]. Additionally, the temporal evolutions of the perturbation are investigated using the finite difference method [20, 65]. We also study the dependence of the resulting quasinormal frequencies on the spacetime configurations. Moreover, we use the Prony method [66, 67] to extract the complex frequencies and compare them to those obtained using the WKB approximation and matrix method. The matrix method discretizes the master equation and approximates it using a matrix equation that can eventually be solved by a nonlinear equation solver. Of course, the QNMs can also be investigated using other well-known approaches, such as the continued fraction method [13]. Specifically, the parameterization given by Eq. (17) is analytic; therefore, we expect that we can derive a three-term recurrence relation and apply the method. However, for this study, we focus on the three methods discussed above. Regarding the QNM frequencies, we are interested in the role of the magnetic field. In addition to the choice of s=1 and Q=0 as discussed earlier, we consider a simplied form of the metric by assuming the following metric parameters: w=109, κ=1, Ca=9, and λ=16. The resultant metric is expressed as

    f=12Mr+310r3.

    (17)

    The resultant quasinormal frequencies are presented in Tables 12 and Figs. 14. From Tables 12, we first conclude that the results obtained using the three different approaches are in reasonable agreement. In particular, the results of the matrix method and sixth-order WKB approximation are mostly consistent. We observe that the term meff=mqB functions as an effective "mass" of the initially massless scalar field. For positive effective mass meff, as the effective mass increases, the real part of the quasinormal frequency increases, whereas the magnitude of the imaginary part decreases. This indicates that the period of the dissipative oscillations decreases, whereas the amplitude decreases slower in time. At vanishing magnetic quantum number m=0, as expected, the quasinormal frequencies are independent of the strength of the magnetic field B.

    Table 1

    Table 1.  QNMs of the massless charged scalar field in the magnetized black hole metric in Rastall gravity. The sixth order WKB approximation and matrix method have been employed for the calculations, by assuming the metric given by Eq. (17).
    lmqB=0qB=0.05qB=0.1qB=0.15
    WKBMatrix methodWKBMatrix methodWKBMatrix methodWKBMatrix method
    110.298-0.096i0.298-0.096i0.277-0.108i0.281-0.109i0.256-0.120i0.266-0.127i0.236-0.130i0.200-0.178i
    00.298-0.096i0.298-0.096i0.298-0.096i0.298-0.096i0.298-0.096i0.298-0.096i0.298-0.096i0.298-0.096i
    −10.298-0.096i0.298-0.096i0.320-0.083i0.320-0.085i0.342-0.069i0.343-0.069i0.363-0.052i0.363-0.053i
    220.492-0.095i0.492-0.095i0.462-0.105i0.464-0.104i0.433-0.115i0.438-0.115i0.404-0.124i0.438-0.115i
    10.492-0.095i0.492-0.095i0.477-0.100i0.477-0.100i0.462-0.105i0.464-0.104i0.448-0.111i0.451-0.108i
    00.492-0.095i0.492-0.095i0.492-0.095i0.492-0.095i0.492-0.095i0.492-0.095i0.492-0.095i0.492-0.095i
    −10.492-0.095i0.492-0.095i0.508-0.090i0.508-0.091i0.523-0.085i0.524-0.085i0.539-0.080i0.539-0.079i
    −20.492-0.095i0.492-0.095i0.523-0.085i0.524-0.085i0.555-0.074i0.555-0.073i0.587-0.062i0.587-0.063i
    DownLoad: CSV
    Show Table

    Table 2

    Table 2.  QNMs of the massless charged scalar field in the magnetized black hole metric in Rastall gravity obtained using the code of Prony's method developed in [66].
    lmqB=0qB=0.05qB=0.1qB=0.15
    110.293-0.098i0.272-0.113i0.247-0.108i0.229-0.206i
    00.293-0.098i0.293-0.098i0.293-0.098i0.293-0.098i
    −10.293-0.098i0.316-0.084i0.339-0.070i0.370-0.051i
    220.484-0.097i0.453-0.108i0.427-0.120i0.377-0.114i
    10.484-0.097i0.468-0.102i0.453-0.108i0.438-0.113i
    00.484-0.097i0.484-0.097i0.484-0.097i0.484-0.097i
    −10.484-0.097i0.500-0.091i0.516-0.086i0.532-0.080i
    −20.484-0.097i0.516-0.086i0.548-0.074i0.581-0.061i
    DownLoad: CSV
    Show Table

    Figure 1

    Figure 1.  (color online) Calculated temporal evolutions of charged massless scalar perturbations for different values of the effective mass. The results are obtained using the finite difference method using the paramters l=1 and M=1.

    Figure 4

    Figure 4.  (color online) Effective potentials Veff(r) for l=2 and different effective mass.

    These results are also in agreement with the temporal evolution calculated using the finite difference method shown in Fig. 1. In particular, as also shown in Fig. 2, as the value of meff increases, the imaginary part eventually becomes insignificant as the envelope of the oscillations tends to lie horizontally. It is depicted explicitly in Fig. 3, where we calculate the quasinormal frequencies as functions of the effective mass using the matrix method. As the effective increase, the real part of the quasinormal frequency increases, and the imaginary part decreases to approach zero. Subsequently, the resulting temporal evolutions are featured by the so-called quasi-resonance [68, 69]. When we further increase the effective mass until it exceeds a critical value, the quasinormal oscillation will disappear entirely. The latter is demonstrated in Fig. 3 as the quasinormal spectrum approaches the real axis as the magnetic field increases. As discussed in the next section, such a configuration becomes favorable for investigating the properties of the late-time tail in a more transparent fashion.

    Figure 2

    Figure 2.  (color online) Calculated temporal evolutions of charged massless scalar perturbations, where instability is observed for a negative effective mass. The results are obtained using the finite difference method for the parameters l=1 and M=1.

    Figure 3

    Figure 3.  (color online) QNMs of different effective masses calculated using the matrix method for the massless charged scalar field with l=1,2, m=1, and n=0.

    For negative meff<0, the WKB and matrix methods fail to produce meaningful results, and we resort to the finite difference method entirely. This is signaled by the fact that the corresponding effective potential has a region of negative values, which indicates possible instability. The form of the relevant effective potential, governed by Eq. (16), is shown in Fig. 4. Similar to Ref. [45], we observe that such a configuration indeed results in instability by exploring the temporal evolution of the initial perturbations. To further verify this instability, we use the finite difference method to determine the temporal evolution of the massless charged scalar perturbations. As demonstrated in Fig. 2, as we gradually decrease the value of the effective mass, the oscillation eventually becomes unstable. Moreover, note that there is a distinction between the quasinormal oscillation and the late-time tail in the limit when the magnetic field vanishes, which is elaborated in the next section.

    In this section, we explore the late-time tails of the charged massless scalar perturbations using the analytic method proposed in Ref. [19] and the numerical integration.

    In terms of the Green's function, the wave function of the scalar field can be given by [70]

    ψ(r,t)=[G(r,r;t)ψt(r,0)+Gt(r,r;t)ψ(r,0)]dr,

    (18)

    where the initial conditions are encoded in ψ(r,0) and ψttψ(r,0), and the Green's function G(r,r;t) is defined by

    [2t22r2+V]G(r,r;t)=δ(t)δ(rr),

    (19)

    When ˜G(r,r;ω) is obtained, the evolution of an arbitrary initial perturbation can be obtained through the integration given by Eq. (18). Indeed, both the QNM and late-time tail can be understood in terms of the singularities and branch cuts of the Fourier transform of the Green's function.

    ˜G(r,r;ω)=+0G(r,r;t)eiωtdt,

    (20)

    whose inverse transform is given by

    G(r,r;t)=12π+ic+ic˜G(r,r;ω)eiωtdω,

    (21)

    where c is a positive constant. Therefore, the analysis of the properties of the Green's function becomes an essential task.

    It is well-known that the Green's function ˜G(r,r;ω) can be constructed using two independent solutions of the corresponding homogeneous equation

    (d2dr2+ω2V)˜ψi=0.

    (22)

    These two solutions, ˜ψ1(r,ω) and ˜ψ2(r,ω), satisfy the appropriate boundary condition at both spatial boundaries [71], namely,

    ˜ψ1(r,ω){eiωrrA(ω)eiωr+B(ω)eiωrr+,

    (23)

    ˜ψ2(r,ω){C(ω)eiωr+D(ω)eiωrre+iωrr+.

    (24)

    We can easily show that the following form satisfies Eq. (19):

    ˜G(r,r;ω)=1W(ω){˜ψ1(r,ω)˜ψ2(r,ω),r>r˜ψ1(r,ω)˜ψ2(r,ω),r<r,

    (25)

    where

    W(ω)˜ψ1(r,ω)˜ψ2,r(r,ω)˜ψ2(r,ω)˜ψ1,r(r,ω)

    (26)

    is the Wronskian of ˜ψ1(r,ω) and ˜ψ2(r,ω).

    Research shows [24] that the late-time tails of quasinormal oscillations can be attributed to the branch cut of the Green's function. For the massless scalar field, the branch cut lies on the negative part of the imaginary axis [21]. For the massive scalar field, Hod and Piran [27] were the first to indicate that the branch cut lies on the real axis between the two branch points governed by the mass. Using reasonable approximation, an appropriate estimation of the contribution from the branch cut subsequently yields the main characteristic of the late-time tail. From a physical perspective, because the low-frequency waveforms are more likely to be backscattered by the curvature at the asymptotic infinity, the late-time tail is primarily dominated by the low-frequency contributions [27]. Mathematically, this is because both types of late-time tails originate from the branch cut, which remains closer to the real axis than any QNM pole.

    To proceed, we introduce

    Rlm(r,ω)=(12Mr+r(2+a))1/2˜ψ.

    (27)

    By further referring to Refs. [27, 29, 72, 73], we assume that the initial perturbation and observer are both located away from the black hole. Subsequently, the master equation can be simplified by expanding Eq. (15) in Mr. Specifically, we ignore the terms of order crα with α>2 and obtain

    [d2dr2+ω2+Bqm+4Mω2+Bqm2Mrl(l+1)r2]˜ψ=0.

    (28)

    By further introducing

    ˜ψ=zl+1ez2Φ(z),z=2Bqmω2r2ϖr,λ=M(Bqm)ϖ2Mϖ,

    (29)

    we can express two relevant solutions for homogeneous equation Eq. (22) as follows

    ˜ψ1=AMλ,(l+12)(2ϖr),˜ψ2=BWλ,(l+12)(2ϖr),

    (30)

    where A=Cϖ(l+1) and B,C are some one-valued even functions of ϖ. Mk,m and Wk,m are the Whittaker functions, which can be expressed in terms of the solutions of Kummer's equation M(a,b,z) and U(a,b,z), namely, the confluent hypergeometric functions [74],

    Mλ,(l+12)(2ϖr)=eϖr(2ϖr)(l+12)+12M(l+1λ,2l+2,2ϖr),Wλ,(l+12)(2ϖr)=eϖr(2ϖr)(l+12)+12U(l+1λ,2l+2,2ϖr).

    (31)

    Note that the specific choice of C guarantees the proper in-going wave boundary condition at the horizon. Moreover, as discussed later, ˜ψ1 satisfies Eq. (38) as it does not contain any discontinuity when crossing the branch cut [29, 72, 73]. In contrast, ˜ψ2 yields the outgoing wave at spatial infinity that has two branch points and, subsequently, a branch cut joining them. As initially proposed in [27], we may conveniently place the branching cut on the real axis of the ω-plane, as shown in Fig. 5. In this context, the late-time tail is primarily governed by the properties of ˜ψ2. Note that it is mandatory to show that the branch cut is in the lower half of the complex plane. Moreover, for the result to be physically meaningful, the specific choice of the branch cut's location must also be irrelevant. A detailed account of this point was discussed recently in Ref. [30].

    Figure 5

    Figure 5.  Relevant contour of the integration on the complex frequency plane of the Green's function. The inverse Fourier transform corresponds to the integral along the upper side of the real axis. It can be complemented by an integral along a large semicircle in the lower half of the complex plane, which is subsequently deformed to go around the quasinormal poles and brach cuts indicated in the plot.

    By using the explicit forms of Eqs. (30)–(31), we can evaluate the contribution from the difference between the two sides of the branch cut. Therefore, the inverse Fourier transform Eq. (21) yields

    GC(r,r;t)=12πBqmBqmF(ϖ)eiωtdω,

    (32)

    where

    F(ϖ)˜ψ1(r,ϖeiπ)˜ψ2(r,ϖeiπ)W(ϖeiπ)˜ψ1(r,ϖ)˜ψ2(r,ϖ)W(ϖ).

    (33)

    Regarding the integrand in Eq. (32), we first calculate the Wronskian in the denominator using the following relations [74]

    Wλ,l+12(2ϖr)=Γ(2l1)Γ(lλ)Mλ,l+12(2ϖr)+Γ(2l+1)Γ(l+1λ)Mλ,(l+12)(2ϖr).

    (34)

    and the non-vanishing Wronskian between the confluent hypergeometric functions

    W{Mλ,l+12(2ϖr),Mλ,(l+12)(2ϖr)}=(2l+1)(2ϖ).

    (35)

    We observe that

    W(ϖ)=AB(2l1)(2ϖ)Γ(2l+1)Γ(l+1λ).

    (36)

    To account for the discontinuity of the integrand across the branch cut, we use

    M±λ,l+12(eiπ2ϖr)=e(l+1)iπMλ,l+12(2ϖr),

    (37)

    to evaluate ˜ψ1,2(r,ϖeiπ) as follows:

    ˜ψ1(r,ϖeiπ)=˜ψ1(r,ϖ),

    (38)

    ˜ψ2(r,ϖeiπ)=B[e(l+1)iπΓ(2l1)Γ(l+λ)Mλ,l+12(2ϖr)+e(l)iπΓ(2l+1)Γ(l+1+λ)Mλ,(l+12)(2ϖr)].

    (39)

    In the above results, note that the subscripts λ of the Whittaker functions remain unchanged. Eq. (38) manifestly indicates that ˜ψ1 does not have any branch cut as it is also a generic physical requirement [21, 27, 29]. Additionally, the Wronskian yields

    W(ϖeiπ)=AB(2l1)(2ϖ)Γ(2l+1)Γ(l+1+λ)e(l)iπ

    (40)

    We note that in deriving the above result, we should always take note of Eqs. (38) and (39). Specifically, the discontinuity associated with the factor ϖ(l+1) embeded in A does not manifest itself, whereas the factor eliπ on the r.h.s. of Eq. (39) does.

    By consolidating all the pieces, i.e., Eqs. (38), (39), and (40), we observe that the integrand of Eq. (32) yields

    F(ϖ)=1(2l1)(2ϖ)A2˜ψ1(r,ϖ)˜F(ϖ)˜ψ1(r,ϖ)

    (41)

    where

    ˜F(ϖ)=e(2l+1)iπΓ(2l1)Γ(l+1+λ)Γ(l+λ)Γ(2l+1)Γ(2l1)Γ(l+1λ)Γ(lλ)Γ(2l+1).

    (42)

    Here, owing to Eq. (38), ˜ψ1(r,ϖ) can be readily factorized out. The remaining term, ˜ψ2(r,ϖ), can be simplified by removing irrelevant contributions that have the identical discontinuity of the Wronskian. As a result, a second ˜ψ1(r,ϖ) can be removed.

    For the intermediate-late time scale, which corresponds to the range MrtM/(qBmM2), the frequency ϖ=O(qBm/t) results in a major contribution to the integral, which implies λ1 owing to Eq. (29). In this case, the term proportional to 1r in Eq. (28), which describes the backscattering due to the curvature of asymptotic infinity, is insignificant [28, 29]. Subsequently, the inverse Fourier transform of the Green's function at the limit λ1 can be performed mostly analytically, which results in

    GC(r,r;t)=(1e(2l+1)iπ)Γ(2l1)Γ(l+1)Γ(l+32)(qBm)2l+14(rr)l+1t2l+32π(4l+2)23l52Γ(2l+1)Γ(l)×cos[qBmtπ(2l+34)].

    (43)

    In other words, we observe that at the intermediate-late time, the tail is dominated by the form

    tl32cos[qBmtπ(l2+34)].

    It oscillates while decaying through a power-law form. In addition to the magnetic field, the intermediate-late time tail depends on the multipole number l and azimuthal number m. In particular, the presence of the magnetic field results in an effective mass that significantly modifies the properties of the tail. However, if we remove the magnetic field, the tail governed by taking the limit B0 in Eq. (43) does not simply fall back to its massless counterpart. Specifically, while the synodal oscillation is suppressed, the remaining power-law exponent becomes el32, which differs from the late-time tail e2l3 observed for a massless scalar in the Schwarzschild background. Moreover, as discussed later, at the limit of the vanishing magnetic field, the power-law exponent of the late-time tail does not match that of its massless counterpart. We postpone further discussions of this intriguing characteristic to the end of this section.

    Now, we discuss the asymptotic behavior of the late-time tails. At a significant time scale qBmt1/(qBmM)2, the relevant terms in the master equation behave differently from the case of intermediate-late time tail. The backscattering from the spacetime curvature at asymptotic infinity can no longer be ignored for this scenario. Thus, because both the ˜ψ1 factors in Eq. (41) do not depend sensitively on the frequency ω, we only require to focus on the term Eq. (42). By applying the above approximation and using the asymptotical forms of Γ functions [74], we obtain

    ˜F(ϖ)Γ(2l1)Γ(2l+1)λ2l+1[e(2l+1)iπη+eiπλ+ηeiπληeiπλ+η+eiπλ],

    (44)

    where

    η±=e±iπl.

    (45)

    The resulting integral of the inverse Fourier transform is somewhat complicated, but its asymptotic properties are governed by the terms related to λ. At small frequencies, we obtain ϖ0 and λ. Their product, which appears in the argument of the Bessel functions, yields

    λϖ=(M(Bqm)ϖ2Mϖ)ϖ=M(Bqm)+O(ϖ2),

    which is finite and varies gradually in the low-frequency region ωBqm+ [72]. Both the terms λ2l+1 and e±iπλ of Eq. (44) oscillate significantly as λ, whereas the latter type is even more dramatic compared with the former. This indicates that the integral's main contribution results from these terms. In most cases, we may follow the arguments of Koyama and Tomimatsu [28] that the main contribution of the integral can be obtained by employing the method of steepest descent. The strategy of such a treatment involves isolating the part that oscillates significantly and then encountering the saddle point, where the oscillations evolve the slowest. The remaining part might be tedious in form but varies moderately. We rewrite the integrand as

    ˜F(ϖ)eiωtΓ(2l1)Γ(2l+1)λ2l+1eiϕei(2πλωt),

    (46)

    where the phase can be defined by

    eiϕ=η++ηe2iπλη++ηe2iπλ.

    (47)

    As shown in [30], the variation in the phase ϕ precisely cancels that in the term e2iπλ. In this case, the integral cannot be evaluated using the method of steepest descent. Nonetheless, we may still evaluate the most dominant term of the integral and obtain

    GC(r,r;t)t1sin(μt+φ),

    (48)

    where the phase shift φ contains some minor dependence on t. The temporal dependence of Eq. (48) implies that late-time tail evaluated regarding the contribution from the saddle point is t1. Note that the obtained exponential does not depend on either the angular momentum number or the strength of the magnetic field.

    We present the numerical results in Figs. 69. To confirm the above analytical results, we perform the numerical calculations using numerical integration. The resultant temporal evolutions of the intermediate-late time are presented in Figs. 67. In Fig. 6, we show the temporal evolutions at the intermediate-late time for different effective masses. For a given l, through Eq. (43), the period of the oscillations increases as the effective mass decreases. Moreover, Fig. 7 shows that when the effective mass is fixed, the oscillating frequency largely remains unchanged for different values of l. Furthermore, the attenuation becomes more significant as l increases. In Fig. 7, the time profiles obtained numerically are also compared against the analytic ones shown in dashed curves given by Eq. (43). The extracted exponent from the envelope of the oscillation yields tβ, where β=2.5087±0.0020 and β=3.5015±0.0020 with respect to l=1, l=2. Therefore, the asymptotic form t(l+3/2) for the intermediate-late time agrees well with the numerical results for different l values.

    Figure 6

    Figure 6.  (color online) Calculated intermediate-late time tail in the Green's function of charged massless scalar perturbations in a magnetized black hole in Rastall gravity. The results are obtained using numerical integration for different effective masses with l=1,2.

    Figure 9

    Figure 9.  (color online) Calculated quasinormal oscillations and late-time tails of charged massless scalar perturbations in a magnetized black hole in Rastall gravity. The calculations aim to explore the asymptotic behavior of the QNMs and tails as the effective mass gradually vanishes. The results are obtained using the parameters l=1 and M=1. Top: The quasinormal oscillations are shown to converge to the case of massless perturbations as the effective mass approaches zero. Bottom: The envelopes of the tails are compared with the analytic forms tl3/2 and t2l3, respectively. The results indicate that the tails' asymptotic behavior does not fall back to that of massless perturbations, as discussed in the text.

    Figure 7

    Figure 7.  (color online) Calculated intermediate-late time tail in the Green's function of charged massless scalar perturbations in a magnetized black hole in Rastall gravity. The results are obtained using numerical integration for different effective masses with l=1,2. The envelopes of dissipative oscillations shown in dashed curves are compared with the analytic ones t(l+3/2) given in the text.

    The temporal evolutions of the late-time tails are presented in Fig. 8. The calculations are performed for different effective masses and angular momenta. Again, by comparing the two plots in Fig. 8, as the effective mass increases, we observe that the oscillation frequency also increases, but the attenuation rate remains the same. Moreover, as shown in the bottom plot of Fig. 8, neither the attenuation rate nor the oscillation period depends on the angular momentum l. We also show in dashed curves the nonlinear fitting to the envelope of the oscillations. The extracted exponents from the a power-law form tβ for the two cases are β=0.9978±0.0068 and β=0.9762±0.0082, respectively, which are reasonably consistent with the analytic result t1.

    Figure 8

    Figure 8.  (color online) Calculated late-time tail in the Green's function of charged massless scalar perturbations in a magnetized black hole in Rastall gravity. The results are obtained using numerical integration for different effective masses with the parameters l=1,2 and M=1. The envelopes of dissipative oscillations shown in dashed curves are compared with the analytic ones t1 given in the text.

    Finally, we elaborate on the specific role of the magnetic field in this study. A distinct characteristic of the present spacetime configuration is that the effective mass is generated by the magnetic field; therefore, the effect on the late-time tail is dynamic. In particular, we are allowed to consider a process in which one gradually tunes the magnetic field to approach the limit B0. Based on the preceding discussions, although the quasinormal frequencies are observed to vary continuously, we conclude that the resultant tail would not converge to its massless counterpart. This phenomenon is manifestly shown in Fig. 9. The top plot of Fig. 9 shows that as the effective mass decreases, the dissipative oscillations converge to the case of massless perturbations. This is consistent with the numerical results obtained in Table 3, where the Prony method [66, 67] has again been employed to extract the complex frequencies from the time profiles. In contrast, as shown in the bottom plot of Fig. 9, the late-time tail of massless perturbations is distinct from the remaining ones. Specifically, the massless perturbation indicated by the orange curve does oscillate at a finite period, whereas the oscillation periods of massive perturbations are more significant (c.f., the green and dark red curves). Regarding the attenuation rate, all the curves of massive perturbations are essentially parallel, which differs from that of the massless one. Quantitatively, the values extracted from the numerical fitting deviate from one another, as indicated in the caption of the plot.

    Table 3

    Table 3.  QNMs for the metric Eq. (17) are calculated by gradually reducing the effective mass. The Prony method, sixth order WKB approximations, and matrix method are used for the calculations. The results are obtained for M=1,l=1.
    meffPronyWKBMatrix method
    00.29861-0.09614i0.29842-0.09636i0.29845-0.09622i
    0.000010.29861-0.09614i0.29842-0.09636i0.29845-0.09622i
    0.00010.29865-0.09612i0.29846-0.09634i0.29849-0.09620i
    0.0010.29904-0.09590i0.29885-0.09611i0.29884-0.09600i
    0.010.30288-0.09421i0.30276-0.09384i0.30246-0.09408i
    DownLoad: CSV
    Show Table

    We now argue that the case of massless perturbation corresponds to a scenario in which the magnetic field does not exist in the first place. Such an interesting difference is somewhat reminiscent of the demagnetization process in ferromagnetic materials. Mathematically, this can be understood as follows. The convergence of the quasinormal frequencies can be understood as the analytic properties of the effective potential essentially remaining unchanged as the effective mass varies continuously. However, the origin of a massive field's late-time tail is associated with the presence of two branch points on the real axis and the absence of any branch point at infinity. As the effective mass approaches zero, the two branch points become degenerate, but a new branch cut will not appear, which extends to infinity. Moreover, a single branch point will be located at the origin in the case of a massless field. The branch cut on the negative imaginary axis connects the origin with the second branch point at infinity. As a result, the relevant structure on the frequency plane due to a massless field is intrinsically different from that of a massive perturbation, although the mass is insignificant. In other words, if we gradually remove the magnetic field from a magnetized black hole, the perturbation will evolve differently from the spacetime configuration in which the magnetic field in question does not ever exist. We note that the above picture is plausible only when the time scale of the demagnetization process is significantly smaller than the inverse of the typical distance between adjacent quasinormal frequencies. This is because the resulting effective potential is now a function of time. Therefore, the notion that QNMs correspond to the zeros of the Wronskian and subsequently the poles of the frequency domain Green's function [3] is only approximately valid when the process can be considered quasistatic. In literature, the master equation with a time-dependent potential is mostly treated numerically [7577].

    In this study, we investigate the QNMs and late-time tails of charged massless scalar perturbations in a magnetized black hole in Rastall gravity. We show that the magnetic field has a significant role in both aspects of the resultant properties of the scalar perturbation. Specifically, the massless scalar acquires an effective mass through the magnetic field, which significantly affects the temporal evolution of the initial perturbations. For the quasinormal oscillations, the complex frequencies are distorted and might become either quasinormally resonant or unstable for specific parameters. For the intermediate and late-time tails, power-law forms are obtained analytically and numerically, reminiscent of those of massive scalar perturbations. Moreover, owing to the dynamic nature of the effective mass generated by the magnetic field, we argue that it has an interesting characteristic. In particular, as we gradually reduce the external magnetic field, while the quasinormal frequencies converge to its massless counterpart, the behavior of the late-time tail deviates from the latter. From a mathematical perspective, this can be understood by analyzing the structure of branch points and cuts on the complex frequency plane of the relevant Green's function.

    We wish to thank Rui-Hong Yue for enlightening discussions.

    In this Appendix, we show that when the magnetic field is small, the metric given by Eqs. (3)–(9) is a valid approximation. For simplicity, we assume that the strength of the magnetic field is insignificant and does not consider the "backreaction" to the metric. First, we readily verify that Eq. (9) satisfies the Maxwell equation for linear electromagnetic field in Schwarzschild's solution of Einstein gravity. For the Schwarzschild solution Eq. (8) of generalized Rastall gravity, we must ascertain that the deviation for Eq. (9) as linear electromagnetic field (s=1) is of higher-order, consistent with the approximation performed later. To show this, we rewrite the deviation from Eq. (17) as

    ˜f=f+δf,

    and solve the Maxwell equation for δf.

    Evidently, δf satisfies

    2δf+2rδf=65r3,

    which implies the magnitude of the correction δfO(r3), irrelevant to the effective potential in the master equation Eq. (28).

    We delegate the validity of this approximation to the Appendix.

    This approximation is also justified by the consistent results obtained from straightforward numerical calculations.

    [1] S. A. T. Stuart L. Shapiro, Black Holes, White Dwarfs, and Neutron Stars (John Wiley and Sons, Ltd, 1983)
    [2] B. P. Abbott et al. (LIGO Scientific Collaboration and Virgo Collaboration), Phys. Rev. Lett. 116, 061102 (2016) doi: 10.1103/PhysRevLett.116.061102
    [3] H.-P. Nollert, Class. Quant. Grav. 16, R159 (1999) doi: 10.1088/0264-9381/16/12/201
    [4] LISA, P. Amaro-Seoane et al., (2017), arXiv: 1702.00786
    [5] J. Luo et al., Class. Quant. Grav. 33, 035010 (2016) doi: 10.1088/0264-9381/33/3/035010
    [6] E. Berti, V. Cardoso, and C. M. Will, Phys. Rev. D 73, 064030 (2006), arXiv:gr-qc/0512160 doi: 10.1103/PhysRevD.73.064030
    [7] H. Liu, C. Zhang, Y. Gong et al., Phys. Rev. D 102, 124011 (2020), arXiv:2002.06360 doi: 10.1103/PhysRevD.102.124011
    [8] T. Regge and J. A. Wheeler, Phys. Rev. 108, 1063 (1957) doi: 10.1103/PhysRev.108.1063
    [9] F. J. Zerilli, Phys. Rev. Lett. 24, 737 (1970) doi: 10.1103/PhysRevLett.24.737
    [10] V. Moncrief, Annals Phys. 88, 323 (1974) doi: 10.1016/0003-4916(74)90173-0
    [11] S. Teukolsky, Astrophys. J. 185, 635 (1973) doi: 10.1086/152444
    [12] S. Chandrasekhar and S. Gottlober, Astron. Nachr. 306, 128 (1985)
    [13] E. W. Leaver, Proc. Roy. Soc. Lond. A 402, 285 (1985)
    [14] E. W. Leaver, Phys. Rev. D 41, 2986 (1990)
    [15] K. Lin and W.-L. Qian, (2016), arXiv: 1609.05948
    [16] K. Lin and W.-L. Qian, Class. Quant. Grav. 34, 095004 (2017), arXiv:1610.08135 doi: 10.1088/1361-6382/aa6643
    [17] K. Lin, W.-L. Qian, A. B. Pavan et al., Mod. Phys. Lett. A 32, 1750134 (2017), arXiv:1703.06439
    [18] K. Lin and W.-L. Qian, Chin. Phys. C 43, 035105 (2019), arXiv:1902.08352 doi: 10.1088/1674-1137/43/3/035105
    [19] E. W. Leaver, Phys. Rev. D 34, 384 (1986) doi: 10.1103/PhysRevD.34.384
    [20] C. Gundlach, R. H. Price, and J. Pullin, Phys. Rev. D 49, 883 (1994) doi: 10.1103/PhysRevD.49.883
    [21] E. S. C. Ching, P. T. Leung, W. M. Suen et al., Phys. Rev. D 52, 2118 (1995), arXiv:gr-qc/9507035 doi: 10.1103/PhysRevD.52.2118
    [22] S. Chen and J. Jing, Mod. Phys. Lett. A 23, 359 (2008), arXiv:gr-qc/0511098 doi: 10.1142/S0217732308023888
    [23] C.-Y. Shao et al., Mod. Phys. Lett. A 35, 2050193 (2020), arXiv:2005.00674
    [24] R. H. Price, Phys. Rev. D 5, 2419 (1972)
    [25] V. Cardoso, S. Yoshida, O. J. C. Dias et al., Phys. Rev. D 68, 061503 (2003), arXiv:hep-th/0307122 doi: 10.1103/PhysRevD.68.061503
    [26] S. Hod and T. Piran, Phys. Rev. D 58, 024017 (1998), arXiv:gr-qc/9712041 doi: 10.1103/PhysRevD.58.024017
    [27] S. Hod and T. Piran, Phys. Rev. D 58, 044018 (1998), arXiv:gr-qc/9801059 doi: 10.1103/PhysRevD.58.044018
    [28] H. Koyama and A. Tomimatsu, Phys. Rev. D 63, 064032 (2001) doi: 10.1103/PhysRevD.63.064032
    [29] H. Koyama and A. Tomimatsu, Phys. Rev. D 65, 084031 (2001), arXiv:gr-qc/0112075
    [30] W.-L. Qian, K. Lin, C.-Y. Shao et al., (2022) arXiv: 2203.04477[gr-qc]
    [31] P. R. Brady, C. M. Chambers, W. Krivan et al., Phys. Rev. D 55, 7538 (1997) doi: 10.1103/PhysRevD.55.7538
    [32] P. R. Brady, C. M. Chambers, W. G. Laarakkers et al., Phys. Rev. D 60, 064003 (1999) doi: 10.1103/PhysRevD.60.064003
    [33] C. Molina, D. Giugno, E. Abdalla et al., Phys. Rev. D 69, 104013 (2004) doi: 10.1103/PhysRevD.69.104013
    [34] N. Varghese and V. C. Kuriakose, Mod. Phys. Lett. A 29, 1450113 (2014), arXiv:1407.6292
    [35] V. Cardoso, E. Franzin, and P. Pani, Phys. Rev. Lett. 116, 171101 (2016) doi: 10.1103/PhysRevLett.116.171101
    [36] T. A. Jacobson and R. Parentani, Sci. Am. 293, 68 (2005)
    [37] H. Liu et al., Phys. Rev. D 104, 044012 (2021), arXiv:2104.11912 doi: 10.1103/PhysRevD.104.044012
    [38] Z. Mark, A. Zimmerman, S. M. Du et al., Phys. Rev. D 96, 084002 (2017) doi: 10.1103/PhysRevD.96.084002
    [39] J.-L. Han, Chin. J. Astron. Astrophys. Suppl. 6, 02211 (2006), arXiv:astro-ph/0603512
    [40] M. Y. Piotrovich, N. A. Silant'ev, Y. N. Gnedin et al., (2010), arXiv: 1002.4948
    [41] R. M. Wald, Phys. Rev. D 10, 1680 (1974) doi: 10.1103/PhysRevD.10.1680
    [42] F. J. Ernst and W. J. Wild, J. Math. Phys. 17, 182 (1976) doi: 10.1063/1.522875
    [43] R. Brito, V. Cardoso, and P. Pani, Phys. Rev. D 89, 104045 (2014), arXiv:1405.2098 doi: 10.1103/PhysRevD.89.104045
    [44] R. A. Konoplya and R. D. B. Fontana, Phys. Lett. B 659, 375 (2008), arXiv:0707.1156 doi: 10.1016/j.physletb.2007.10.065
    [45] B. Turimov, B. Toshmatov, B. Ahmedov et al., Phys. Rev. D 100, 084038 (2019), arXiv:1910.00939 doi: 10.1103/PhysRevD.100.084038
    [46] P. Rastall, Phys. Rev. D 6, 3357 (1972) doi: 10.1103/PhysRevD.6.3357
    [47] Y. Heydarzade, H. Moradpour, and F. Darabi, Can. J. Phys. 95, 1253 (2017), arXiv:1610.03881 doi: 10.1139/cjp-2017-0254
    [48] R. Kumar and S. G. Ghosh, Eur. Phys. J. C 78, 750 (2018), arXiv:1711.08256 doi: 10.1140/epjc/s10052-018-6206-1
    [49] I. P. Lobo, H. Moradpour, J. P. Morais Graça et al., Int. J. Mod. Phys. D 27, 1850069 (2018), arXiv:1710.04612
    [50] R. Kumar, B. P. Singh, M. S. Ali et al., Phys. Dark Univ. 34, 100881 (2021), arXiv:1712.09793 doi: 10.1016/j.dark.2021.100881
    [51] M. Capone, V. F. Cardone, and M. L. Ruggiero, Nuovo Cim. B 125, 1133 (2011), arXiv:0906.4139
    [52] H. Moradpour, Y. Heydarzade, F. Darabi et al., Eur. Phys. J. C 77, 259 (2017), arXiv:1704.02458 doi: 10.1140/epjc/s10052-017-4811-z
    [53] C. E. M. Batista, J. C. Fabris, O. F. Piattella et al., Eur. Phys. J. C 73, 2425 (2013), arXiv:1208.6327 doi: 10.1140/epjc/s10052-013-2425-7
    [54] H. Moradpour, Phys. Lett. B 757, 187 (2016) doi: 10.1016/j.physletb.2016.03.072
    [55] A. S. Al-Rawaf and M. O. Taha, Phys. Lett. B 366, 69 (1996) doi: 10.1016/0370-2693(95)01145-5
    [56] K. Lin and W.-L. Qian, Eur. Phys. J. C 80, 561 (2020), arXiv:2006.03229 doi: 10.1140/epjc/s10052-020-8116-2
    [57] K. Lin, Y. Liu, and W.-L. Qian, Gen. Rel. Grav. 51, 62 (2019), arXiv:1809.10075 doi: 10.1007/s10714-019-2548-8
    [58] S. H. Hendi, Eur. Phys. J. C 71, 1551 (2011), arXiv:1007.2704 doi: 10.1140/epjc/s10052-011-1551-3
    [59] B. F. Schutz and C. M. Will, Astrophys. J. 291, L33 (1985) doi: 10.1086/184453
    [60] K. D. Kokkotas and B. F. Schutz, Phys. Rev. D 37, 3378 (1988) doi: 10.1103/PhysRevD.37.3378
    [61] S. Iyer and C. M. Will, Phys. Rev. D 35, 3621 (1987)
    [62] R. A. Konoplya, Phys. Rev. D 68, 024018 (2003) doi: 10.1103/PhysRevD.68.024018
    [63] J. Matyjasek and M. Opala, Phys. Rev. D 96, 024011 (2017) doi: 10.1103/PhysRevD.96.024011
    [64] R. A. Konoplya, A. Zhidenko, and A. F. Zinhailo, Class. Quant. Grav. 36, 155002 (2019), arXiv:1904.10333 doi: 10.1088/1361-6382/ab2e25
    [65] C. Gundlach, R. H. Price, and J. Pullin, Phys. Rev. D 49, 890 (1994) doi: 10.1103/PhysRevD.49.890
    [66] W.-L. Qian, K. Lin, J.-P. Wu et al., Eur. Phys. J. C 80, 959 (2020), arXiv:2006.07122 doi: 10.1140/epjc/s10052-020-08539-x
    [67] E. Berti, V. Cardoso, J. A. Gonzalez et al., Phys. Rev. D 75, 124017 (2007), arXiv:gr-qc/0701086 doi: 10.1103/PhysRevD.75.124017
    [68] A. Ohashi and M.-a. Sakagami, Class. Quant. Grav. 21, 3973 (2004), arXiv:gr-qc/0407009 doi: 10.1088/0264-9381/21/16/010
    [69] R. A. Konoplya and A. V. Zhidenko, Phys. Lett. B 609, 377 (2005), arXiv:gr-qc/0411059 doi: 10.1016/j.physletb.2005.01.078
    [70] P. Morse and H. Feshback, Methods of Theoretical PhysicsInternational Series in Pure and Applied Physics (McGraw-Hill Book Company, 1953)
    [71] H.-P. Nollert and B. G. Schmidt, Phys. Rev. D 45, 2617 (1992) doi: 10.1103/PhysRevD.45.2617
    [72] H. Yu, Phys. Rev. D 65, 087502 (2002) doi: 10.1103/PhysRevD.65.087502
    [73] J.-l. Jing, Phys. Rev. D 70, 065004 (2004), arXiv:gr-qc/0405122 doi: 10.1103/PhysRevD.70.065004
    [74] H. Bateman and A. Erdélyi, Higher transcendental functionsCalifornia Institute of technology. Bateman Manuscript project (McGraw-Hill, New York, NY)
    [75] C.-G. Shao, B. Wang, E. Abdalla et al., Phys. Rev. D 71, 044003 (2005), arXiv:gr-qc/0410025
    [76] C. Chirenti and A. Saa, Phys. Rev. D 84, 064006 (2011), arXiv:1105.1681
    [77] K. Lin, Y. Liu, W.-L. Qian et al., Phys. Rev. D 100, 065018 (2019), arXiv:1909.04347
  • [1] S. A. T. Stuart L. Shapiro, Black Holes, White Dwarfs, and Neutron Stars (John Wiley and Sons, Ltd, 1983)
    [2] B. P. Abbott et al. (LIGO Scientific Collaboration and Virgo Collaboration), Phys. Rev. Lett. 116, 061102 (2016) doi: 10.1103/PhysRevLett.116.061102
    [3] H.-P. Nollert, Class. Quant. Grav. 16, R159 (1999) doi: 10.1088/0264-9381/16/12/201
    [4] LISA, P. Amaro-Seoane et al., (2017), arXiv: 1702.00786
    [5] J. Luo et al., Class. Quant. Grav. 33, 035010 (2016) doi: 10.1088/0264-9381/33/3/035010
    [6] E. Berti, V. Cardoso, and C. M. Will, Phys. Rev. D 73, 064030 (2006), arXiv:gr-qc/0512160 doi: 10.1103/PhysRevD.73.064030
    [7] H. Liu, C. Zhang, Y. Gong et al., Phys. Rev. D 102, 124011 (2020), arXiv:2002.06360 doi: 10.1103/PhysRevD.102.124011
    [8] T. Regge and J. A. Wheeler, Phys. Rev. 108, 1063 (1957) doi: 10.1103/PhysRev.108.1063
    [9] F. J. Zerilli, Phys. Rev. Lett. 24, 737 (1970) doi: 10.1103/PhysRevLett.24.737
    [10] V. Moncrief, Annals Phys. 88, 323 (1974) doi: 10.1016/0003-4916(74)90173-0
    [11] S. Teukolsky, Astrophys. J. 185, 635 (1973) doi: 10.1086/152444
    [12] S. Chandrasekhar and S. Gottlober, Astron. Nachr. 306, 128 (1985)
    [13] E. W. Leaver, Proc. Roy. Soc. Lond. A 402, 285 (1985)
    [14] E. W. Leaver, Phys. Rev. D 41, 2986 (1990)
    [15] K. Lin and W.-L. Qian, (2016), arXiv: 1609.05948
    [16] K. Lin and W.-L. Qian, Class. Quant. Grav. 34, 095004 (2017), arXiv:1610.08135 doi: 10.1088/1361-6382/aa6643
    [17] K. Lin, W.-L. Qian, A. B. Pavan et al., Mod. Phys. Lett. A 32, 1750134 (2017), arXiv:1703.06439
    [18] K. Lin and W.-L. Qian, Chin. Phys. C 43, 035105 (2019), arXiv:1902.08352 doi: 10.1088/1674-1137/43/3/035105
    [19] E. W. Leaver, Phys. Rev. D 34, 384 (1986) doi: 10.1103/PhysRevD.34.384
    [20] C. Gundlach, R. H. Price, and J. Pullin, Phys. Rev. D 49, 883 (1994) doi: 10.1103/PhysRevD.49.883
    [21] E. S. C. Ching, P. T. Leung, W. M. Suen et al., Phys. Rev. D 52, 2118 (1995), arXiv:gr-qc/9507035 doi: 10.1103/PhysRevD.52.2118
    [22] S. Chen and J. Jing, Mod. Phys. Lett. A 23, 359 (2008), arXiv:gr-qc/0511098 doi: 10.1142/S0217732308023888
    [23] C.-Y. Shao et al., Mod. Phys. Lett. A 35, 2050193 (2020), arXiv:2005.00674
    [24] R. H. Price, Phys. Rev. D 5, 2419 (1972)
    [25] V. Cardoso, S. Yoshida, O. J. C. Dias et al., Phys. Rev. D 68, 061503 (2003), arXiv:hep-th/0307122 doi: 10.1103/PhysRevD.68.061503
    [26] S. Hod and T. Piran, Phys. Rev. D 58, 024017 (1998), arXiv:gr-qc/9712041 doi: 10.1103/PhysRevD.58.024017
    [27] S. Hod and T. Piran, Phys. Rev. D 58, 044018 (1998), arXiv:gr-qc/9801059 doi: 10.1103/PhysRevD.58.044018
    [28] H. Koyama and A. Tomimatsu, Phys. Rev. D 63, 064032 (2001) doi: 10.1103/PhysRevD.63.064032
    [29] H. Koyama and A. Tomimatsu, Phys. Rev. D 65, 084031 (2001), arXiv:gr-qc/0112075
    [30] W.-L. Qian, K. Lin, C.-Y. Shao et al., (2022) arXiv: 2203.04477[gr-qc]
    [31] P. R. Brady, C. M. Chambers, W. Krivan et al., Phys. Rev. D 55, 7538 (1997) doi: 10.1103/PhysRevD.55.7538
    [32] P. R. Brady, C. M. Chambers, W. G. Laarakkers et al., Phys. Rev. D 60, 064003 (1999) doi: 10.1103/PhysRevD.60.064003
    [33] C. Molina, D. Giugno, E. Abdalla et al., Phys. Rev. D 69, 104013 (2004) doi: 10.1103/PhysRevD.69.104013
    [34] N. Varghese and V. C. Kuriakose, Mod. Phys. Lett. A 29, 1450113 (2014), arXiv:1407.6292
    [35] V. Cardoso, E. Franzin, and P. Pani, Phys. Rev. Lett. 116, 171101 (2016) doi: 10.1103/PhysRevLett.116.171101
    [36] T. A. Jacobson and R. Parentani, Sci. Am. 293, 68 (2005)
    [37] H. Liu et al., Phys. Rev. D 104, 044012 (2021), arXiv:2104.11912 doi: 10.1103/PhysRevD.104.044012
    [38] Z. Mark, A. Zimmerman, S. M. Du et al., Phys. Rev. D 96, 084002 (2017) doi: 10.1103/PhysRevD.96.084002
    [39] J.-L. Han, Chin. J. Astron. Astrophys. Suppl. 6, 02211 (2006), arXiv:astro-ph/0603512
    [40] M. Y. Piotrovich, N. A. Silant'ev, Y. N. Gnedin et al., (2010), arXiv: 1002.4948
    [41] R. M. Wald, Phys. Rev. D 10, 1680 (1974) doi: 10.1103/PhysRevD.10.1680
    [42] F. J. Ernst and W. J. Wild, J. Math. Phys. 17, 182 (1976) doi: 10.1063/1.522875
    [43] R. Brito, V. Cardoso, and P. Pani, Phys. Rev. D 89, 104045 (2014), arXiv:1405.2098 doi: 10.1103/PhysRevD.89.104045
    [44] R. A. Konoplya and R. D. B. Fontana, Phys. Lett. B 659, 375 (2008), arXiv:0707.1156 doi: 10.1016/j.physletb.2007.10.065
    [45] B. Turimov, B. Toshmatov, B. Ahmedov et al., Phys. Rev. D 100, 084038 (2019), arXiv:1910.00939 doi: 10.1103/PhysRevD.100.084038
    [46] P. Rastall, Phys. Rev. D 6, 3357 (1972) doi: 10.1103/PhysRevD.6.3357
    [47] Y. Heydarzade, H. Moradpour, and F. Darabi, Can. J. Phys. 95, 1253 (2017), arXiv:1610.03881 doi: 10.1139/cjp-2017-0254
    [48] R. Kumar and S. G. Ghosh, Eur. Phys. J. C 78, 750 (2018), arXiv:1711.08256 doi: 10.1140/epjc/s10052-018-6206-1
    [49] I. P. Lobo, H. Moradpour, J. P. Morais Graça et al., Int. J. Mod. Phys. D 27, 1850069 (2018), arXiv:1710.04612
    [50] R. Kumar, B. P. Singh, M. S. Ali et al., Phys. Dark Univ. 34, 100881 (2021), arXiv:1712.09793 doi: 10.1016/j.dark.2021.100881
    [51] M. Capone, V. F. Cardone, and M. L. Ruggiero, Nuovo Cim. B 125, 1133 (2011), arXiv:0906.4139
    [52] H. Moradpour, Y. Heydarzade, F. Darabi et al., Eur. Phys. J. C 77, 259 (2017), arXiv:1704.02458 doi: 10.1140/epjc/s10052-017-4811-z
    [53] C. E. M. Batista, J. C. Fabris, O. F. Piattella et al., Eur. Phys. J. C 73, 2425 (2013), arXiv:1208.6327 doi: 10.1140/epjc/s10052-013-2425-7
    [54] H. Moradpour, Phys. Lett. B 757, 187 (2016) doi: 10.1016/j.physletb.2016.03.072
    [55] A. S. Al-Rawaf and M. O. Taha, Phys. Lett. B 366, 69 (1996) doi: 10.1016/0370-2693(95)01145-5
    [56] K. Lin and W.-L. Qian, Eur. Phys. J. C 80, 561 (2020), arXiv:2006.03229 doi: 10.1140/epjc/s10052-020-8116-2
    [57] K. Lin, Y. Liu, and W.-L. Qian, Gen. Rel. Grav. 51, 62 (2019), arXiv:1809.10075 doi: 10.1007/s10714-019-2548-8
    [58] S. H. Hendi, Eur. Phys. J. C 71, 1551 (2011), arXiv:1007.2704 doi: 10.1140/epjc/s10052-011-1551-3
    [59] B. F. Schutz and C. M. Will, Astrophys. J. 291, L33 (1985) doi: 10.1086/184453
    [60] K. D. Kokkotas and B. F. Schutz, Phys. Rev. D 37, 3378 (1988) doi: 10.1103/PhysRevD.37.3378
    [61] S. Iyer and C. M. Will, Phys. Rev. D 35, 3621 (1987)
    [62] R. A. Konoplya, Phys. Rev. D 68, 024018 (2003) doi: 10.1103/PhysRevD.68.024018
    [63] J. Matyjasek and M. Opala, Phys. Rev. D 96, 024011 (2017) doi: 10.1103/PhysRevD.96.024011
    [64] R. A. Konoplya, A. Zhidenko, and A. F. Zinhailo, Class. Quant. Grav. 36, 155002 (2019), arXiv:1904.10333 doi: 10.1088/1361-6382/ab2e25
    [65] C. Gundlach, R. H. Price, and J. Pullin, Phys. Rev. D 49, 890 (1994) doi: 10.1103/PhysRevD.49.890
    [66] W.-L. Qian, K. Lin, J.-P. Wu et al., Eur. Phys. J. C 80, 959 (2020), arXiv:2006.07122 doi: 10.1140/epjc/s10052-020-08539-x
    [67] E. Berti, V. Cardoso, J. A. Gonzalez et al., Phys. Rev. D 75, 124017 (2007), arXiv:gr-qc/0701086 doi: 10.1103/PhysRevD.75.124017
    [68] A. Ohashi and M.-a. Sakagami, Class. Quant. Grav. 21, 3973 (2004), arXiv:gr-qc/0407009 doi: 10.1088/0264-9381/21/16/010
    [69] R. A. Konoplya and A. V. Zhidenko, Phys. Lett. B 609, 377 (2005), arXiv:gr-qc/0411059 doi: 10.1016/j.physletb.2005.01.078
    [70] P. Morse and H. Feshback, Methods of Theoretical PhysicsInternational Series in Pure and Applied Physics (McGraw-Hill Book Company, 1953)
    [71] H.-P. Nollert and B. G. Schmidt, Phys. Rev. D 45, 2617 (1992) doi: 10.1103/PhysRevD.45.2617
    [72] H. Yu, Phys. Rev. D 65, 087502 (2002) doi: 10.1103/PhysRevD.65.087502
    [73] J.-l. Jing, Phys. Rev. D 70, 065004 (2004), arXiv:gr-qc/0405122 doi: 10.1103/PhysRevD.70.065004
    [74] H. Bateman and A. Erdélyi, Higher transcendental functionsCalifornia Institute of technology. Bateman Manuscript project (McGraw-Hill, New York, NY)
    [75] C.-G. Shao, B. Wang, E. Abdalla et al., Phys. Rev. D 71, 044003 (2005), arXiv:gr-qc/0410025
    [76] C. Chirenti and A. Saa, Phys. Rev. D 84, 064006 (2011), arXiv:1105.1681
    [77] K. Lin, Y. Liu, W.-L. Qian et al., Phys. Rev. D 100, 065018 (2019), arXiv:1909.04347
  • 加载中

Cited by

1. Sun, Q., Zhang, Y., Xie, C.-H. et al. Shadow of Kerr black hole surrounded by a cloud of strings in Rastall gravity and constraints from M87*[J]. Physics of the Dark Universe, 2024. doi: 10.1016/j.dark.2024.101599
2. Magalhães, R.B., Masó-Ferrando, A.S., Bombacigno, F. et al. Echoes from bounded universes[J]. Physical Review D, 2024, 110(4): 044058. doi: 10.1103/PhysRevD.110.044058
3. Gogoi, D.J., Heidari, N., K̆rí̆z, J. et al. Quasinormal Modes and Greybody Factors of de Sitter Black Holes Surrounded by Quintessence in Rastall Gravity[J]. Fortschritte der Physik, 2024, 72(3): 2300245. doi: 10.1002/prop.202300245
4. Sun, Q., Li, Q., Zhang, Y. et al. Quasinormal modes, Hawking radiation and absorption of the massless scalar field for Reissner-Nordström black hole surrounded by a cloud of strings in Rastall gravity[J]. International Journal of Geometric Methods in Modern Physics, 2024. doi: 10.1142/S0219887824502426
5. Shen, S.-F., Qian, W.-L., Zhang, J. et al. Matrix method and the suppression of Runge's phenomenon[J]. SciPost Physics Core, 2024, 7(2) doi: 10.21468/SciPostPhysCore.7.2.034
6. Shen, S.-F., Qian, W.-L., Guo, H. et al. An implementation of the matrix method using the Chebyshev grid[J]. Progress of Theoretical and Experimental Physics, 2023, 2023(9): 093E01. doi: 10.1093/ptep/ptad107
7. Toshmatov, B., Stuchlík, Z., Ahmedov, B. Can electromagnetic charge inhabit in Rastall gravity?[J]. Physics of the Dark Universe, 2023. doi: 10.1016/j.dark.2023.101257

Figures(9) / Tables(3)

Get Citation
Cai-Ying Shao, Yu-Jie Tan, Cheng-Gang Shao, Kai Lin and Wei-Liang Qian. Qusinormal oscillations and late-time tail of massless scalar perturbations of a magnetized black hole in Rastall gravity[J]. Chinese Physics C. doi: 10.1088/1674-1137/ac7855
Cai-Ying Shao, Yu-Jie Tan, Cheng-Gang Shao, Kai Lin and Wei-Liang Qian. Qusinormal oscillations and late-time tail of massless scalar perturbations of a magnetized black hole in Rastall gravity[J]. Chinese Physics C.  doi: 10.1088/1674-1137/ac7855 shu
Milestone
Received: 2022-04-19
Article Metric

Article Views(1911)
PDF Downloads(57)
Cited by(7)
Policy on re-use
To reuse of subscription content published by CPC, the users need to request permission from CPC, unless the content was published under an Open Access license which automatically permits that type of reuse.
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Email This Article

Title:
Email:

Quasinormal oscillations and late-time tail of massless scalar perturbations of a magnetized black hole in Rastall gravity

  • 1. MOE Key Laboratory of Fundamental Physical Quantities Measurement, Hubei Key Laboratory of Gravitation and Quantum Physics, PGMF, and School of Physics, Huazhong University of Science and Technology, 430074 Wuhan, China
  • 2. Hubei Subsurface Multi-scale Imaging Key Laboratory, Institute of Geophysics and Geomatics, China University of Geosciences, 430074 Wuhan, China
  • 3. Escola de Engenharia de Lorena, Universidade de São Paulo, 12602-810, Lorena, SP, Brazil
  • 4. Faculdade de Engenharia de Guaratinguetá, Universidade Estadual Paulista, 12516-410, Guaratinguetá, SP, Brazil and
  • 5. Center for Gravitation and Cosmology, School of Physical Science and Technology, Yangzhou University, 225002 Yangzhou, China

Abstract: In this study, we investigate the quasinormal mode and late-time tail of charged massless scalar perturbations of a black hole in generalized Rastall gravity. The black hole metric in question is spherically symmetric, accompanied by a power-Maxwell field surrounded by a quintessence fluid. We show that the massless scalar field, when dressed up with the magnetic field, acquires an effective mass, which significantly affects the properties of the resultant quasinormal oscillations and late-time tails. Specifically, the quasinormal frequencies become distorted and might even be unstable for particular spacetime configurations. Additionally, the exponent of the usual power-law tail is modified according to the modification in the structure of the branch cut of the retarded Green's function. In particular, as the effective mass is generated dynamically owing to the presence of the magnetic field, we may consider a process through which the field is gradually removed from the spacetime configuration. In this context, while the quasinormal oscillations converge to the case of massless perturbations, we argue that the properties of resultant late-time tails do not fall back to their massless counterpart. The relevant characteristics are investigated using numerical and analytic approaches.

    HTML

    I.   INTRODUCTION
    • The black hole is one of the most intriguing predictions of Einstein's general relativity. It is a particular spacetime region typically created by the gravitational collapse of a compact object in its late stage, when its mass exceeds the upper limit of neutron stars [1]. From both theoretical and experimental perspectives, the subject has gained much interest over the past few decades. As the culmination of continuous efforts, the recent advent of empirical detection of gravitational waves emanating from black holes has elicited further interest. In particular, the merger signals of a black hole binary captured by the LIGO and Virgo collaboration [2] furnished the first direct piece of evidence of black holes and subsequently inaugurated a novel era of gravitational-wave astronomy.

      The quasinormal modes (QNMs) carry the characteristic sound of a black hole. In terms of dissipative oscillations, they bear the essential properties of the underlying spacetime metric [3]. The associated temporal profile constitutes the ringdown phase of the collapse, from which a black hole is formed. On the one hand, such signals are mathematically straightforward compared with those emanating from the merger process. On the other hand, the strength of the associated gravitational radiation is speculated to still fall within the scope of the ongoing space-borne gravitational wave detector programs, such as LISA and TianQin [4, 5]. In particular, studies have been performed to evaluate the detector signal-to-noise ratio regarding the feasibility of an eventual detection [6, 7].

      Theoretically, QNMs can be investigated using black hole perturbation theory. Among the earliest endeavors, Regge and Wheeler created the framework of the metric perturbations in spherically symmetric black holes [8]. Subsequent developments by Zerill [9] and Moncrief [10] further refined a few pertinent aspects. For rotating black holes, Teukolsky first derived the linearized scalar and metric perturbation equations for the Kerr spacetime [11, 12]. Although the physical system is dissipative and therefore non-Hermitian, the mathematical procedure to obtain QNMs can be effectively considered an eigenvalue problem. Subsequently, the emergence of discrete complex frequencies is attributed to the ingoing and outgoing boundary conditions, respectively, applied at the horizon and outer spatial bound (often located at infinity). By using the continued fraction method proposed by Leaver [13, 14], we encounter the discrete complex frequencies reminiscent of the atomic spectrum of hydrogen in quantum mechanics. The eigenvalue nature of the problem is also transparently considered in terms of the matrix method [1518], in which we explicitly solve a matrix secular equation for the quasinormal frequencies. In contrast, Leaver proposed the Green's function approach, which reformulates the general solution in terms of the spectrum decomposition [3, 19]. From this perspective, the main characteristics of the perturbation fields are attributed to the intrinsic singularities of the relevant Green's function. Specifically, the QNMs correspond to the poles of the Green's function, whereas the branch cuts essentially govern the late-time tails. In addition to the QNMs, the late-time tail is a topic of pertinent relevance. On the one hand, the time profile of QNMs decays exponentially; on the other hand, the late-time tail largely follows a power-law form. Therefore, the final stage of the temporal evolution of the perturbations is dominated primarily by the late-time tail. As mentioned earlier, mathematically, the origin of the late-time tail is attributed to the branch cut and cannot be obtained by summing the QNM poles. For instance, for massless scalar perturbations, the tail is due to a branch cut on the negative imaginary axis in the frequency domain, which stretches from the origin. Researchers have argued that the presence of the branch cut sensitively depends on the asymptotical structure of the potential [20]. Additionally, it is often interpreted in the literature as the backscattering of perturbed wave packets by the spacetime far away from the horizon [2123]. Price was the first to observe an inverse power-law tail, t(2l+3) or t(2l+2), for massless scalar perturbations in the Schwarzschild black hole metric [24]. Extensive studies for different metrics further indicated that such an inverse power-law form is a rather general characteristic for massless perturbations in spherical spacetimes that are asymptotically flat [21, 25, 26]. Moreover, the temporal profile for massive perturbations might be qualitatively different from their massless counterparts. For insignificant mass, Hod and Piran discovered that the time-domain profile at the intermediate-late time in Reissner-Nordström spacetime has an oscillatory tail with a decay rate of t(l+3/2)sin(μt) [27], where μ is the mass of the field. Moreover, Koyama and Tomimatsu analytically showed that the asymptotical late-time tail in the Schwarzschild spacetime is t5/6sin(μt) [28, 29]. For some specific cases, the asymptotical form was recenetly reported to also be t1sin(μt) [30] in the limit of vanishing scalar mass. In all the massive cases, the branch cut can be conveniently selected to be on the real axis in a finite interval determined by the scalar mass. It is also worth noting that such a decay rate does not depend on the angular momentum l. In addition to the power-law form, the exponential tails have also been observed in asymptotically de Sitter spacetimes for various types of fields [3133]. Notably, the Dirac perturbations eventually decay to a non-zero constant [23, 34]. These distinctive characteristics have significantly enriched the physical content regarding the late-time waveforms. More recently, Cardoso et al. proposed [35] that the late-time ringdown may serve to probe the horizon and discriminate between different gravitational systems, inclusively, the exotic compact objects such as gravastar and wormhole. Based on this, black hole echoes [36], and particularly, their connections with the QNMs and late-time tail [37, 38] have also attracted much interest.

      The magnetic field is another indispensable assembly piece of the universe. Experimentally, employing the rotation measures, information on the magnetic field in the disk of galaxies can be extracted for a large number of pulsars [39]. In particular, strong magnetic fields, as high as 104108 G, have been observed in the vicinity of stellar-mass or supermassive black holes [40]. Therefore, relevant studies on black hole metrics, including their perturbations, have to be placed in the context of a realistic astrophysical environment. A stationary axisymmetric black hole solution with a uniform magnetic field placed along the symmetry axis was first derived by Wald [41]. The solution for the Kerr-Newman black hole was subsequently obtained by Ernst and Wild [42]. Among many important developments, the magnetic field was observed to have an essential role in the stability of the dynamic system. The latter encompass the superradiance instability [43] and black hole QNMs [44, 45]. In particular, the Zeeman effect was recently uncovered in the QNM spectrum of charged scalar perturbations in a magnetized Schwarzschild black hole [44]. Moreover, the instability of the metric can be triggered by tuning the effective scalar mass, which is modified by the interaction between the charge and magnetic field.

      As an alternative theory of general relativity, Rastall gravity was formulated based on the possible breaking of the conservation law of the stress energy-momentum tensor [46]. The theory aims at addressing several unsettled problems in general relativity regarding the largest scale when compared with cosmological observation. In the literature, various intriguing aspects of the theory have been explored, which include black hole physics [4750] and cosmology [5155]. Among others, the theory inherently results in a particle creation mechanism, which naturally supplies an alternative implementation for dark energy [56]. It is meaningful to further investigate the impact of the magnetic field in Rastall theory on spacetime stability.

      This paper is motivated by the above considerations. We investigate the QNMs and late-time tails of charged massless scalar perturbations of a magnetized black hole solution derived from generalized Rastall gravity. First, we generalize the black hole solution in Rastall gravity into a more realistic context by including the magnetic field. We then analyze the effect of the magnetic field on the spacetime stability in the context of Rastall gravity by exploring the quasinormal oscillations. For this, the study of massless scalar perturbations is sufficient for the purpose. For stable spacetime configurations, we know that massless and massive scalar fields behave distinctively in the late-time evolution. Therefore, it would be rather interesting to investigate the implications for Rastall's theory with the presence of a magnetic field. In particular, we observe that the exponent of the power-law tail is modified. This can be attributed to the fact that the effective mass of the field is "dressed up" owing to the interaction between the magnetic field and charge.

      The remainder of the paper is organized as follows. In the Sec. II, we present the metric of a magnetized black hole investigated in this study. The QNMs are explored in Sec. III by employing the Wentzel-Kramers-Brillouin (WKB) approximation, matrix method, and finite difference method. In Sec. IV, we analyze the late-time tail of the massless scalar perturbations and the effect of the magnetic field. Both the Green's function and the finite difference method are utilized. Further discussions and concluding remarks are provided in the final section.

    II.   POWER-MAXWELL CHARGED BLACK HOLE SOLUTIONS IN RASTALL GRAVITY
    • According to Rastall's recipe [46], the breaking of the conservation law of the stress energy-momentum tensor is

      Tvμ;v=λR,μ ,

      (1)

      where λ is the Rastall coupling parameter, and R is the Ricci scalar. If λ=0, the stress energy-momentum tensor is conserved and Rastall gravity falls back to general relativity. Moreover, the corresponding field equation can be expressed as

      Rμν+(κλ12)gμνR=κTμν,

      (2)

      where κ is a constant related to the trace anomaly of the energy momentum tensor in Rastall gravity.

      In this paper, we consider a spherically symmetric black hole solution surrounded by a quintessence fluid and power-Maxwell field first derived in [57]. Specifically, the metric of a static, spherically symmetric four-dimensional spacetime is expressed as

      ds2=f(r)dt2+f(r)1dr2+r2dθ2+r2sin2θdϕ2.

      (3)

      The action for the electromagnetic sector is of a nonlinear form, which can be expressed as [58]

      LF=(ξF)s,

      (4)

      where F=FμνFμv, Fμν is the Faraday tensor, and s and ξ are constants. s=1 corresponds to the case of a linear electromagnetic field. Owing to the presence of the electromagnetic field and quintessence fluid, the resulting energy momentum tensor has the form

      Tμv=Eμv+Tμv,

      (5)

      where the electromagnetic part, Eμv, is given by

      Evμ=(ξ)s(F)s1(2sFσμFσv12δvμF).

      (6)

      Tμv is the energy-momentum tensor of the quintessence field. If we assume a barotropic equation of state p=wρ, it can be expressed as

      Ttt=Trr=ρ(r),Tθθ=Tφφ=12ρ(r)(3w+1).

      (7)

      By substituting the energy momentum tensor into the field equation Eq. (2), we can show that the black hole metric is [57]

      f(r)=12Mr+Q2r212s(1+2s)32s12s(32s)s+Car16κλ+w(36κλ)1+3(1+w)κλ×(13(1+w)κλ)23(1+4κλ)(κλ+w(1+κλ)),

      (8)

      where Ca is a constant of integration, and M is related to the mass of the black hole.

      For this study, we will focus on the case of the linear electromagnetic field with s=1 and the sole presence of a constant magnetic field, namely, the electric charge vanishes, Q=0. When black hole is in the magnetic field, as a good approximation, the vector potential of the electromagnetic field can be expressed as

      Aμ=12Br2sin2θ(0,0,0,1).

      (9)
    III.   QUASINORMAL MODES OF CHARGED MASSLESS SCALAR PERTURBATIONS
    • The temporal evolution of a massless charged scalar field is governed by the Klein-Gordon equation, as follows:

      gμν(μiqAμ)(νiqAν)Φ=0.

      (10)

      To simplify the above equation, we introduce some reasonable assumptions [45]. First, we obtain αAα=0 by selecting the Lorentz gauge. Additionally, we assume that q2B21 by considering the scenario in which the coupling between the scalar and electromagnetic fields is sufficiently weak. Under these conditions, Eq. (10) can be rewritten as

      1r2r(r2fΦr)(L2r2qBLz)Φ1f2Φt2=0,

      (11)

      where

      L2=[1sinθθ(sinθθ)+1sin2θ2ϕ2]

      (12)

      and

      Lz=iϕ.

      (13)

      Now, we may employ the method of separation of variables by expressing the wave function as

      Φ(t,r,θ,ϕ)=12πdωeiωtlRlm(r,ω)rYlm(θ,ϕ),

      (14)

      where Ym(θ,ϕ) are the spherical harmonics, and the radial part of the wave function Rm(r,ω) satisfies

      (d2dr2+ω2Veff(r))Rlm(r,ω)=0,

      (15)

      where the effective potential is given by

      Veff(r)=f(r)[l(l+1)r2+f(r)rmqB],

      (16)

      and r is the tortoise coordinate defined by dr=drf(r). We note that Veff=mqB as r.

      In the following, we evaluate the QNMs using the sixth order WKB approximation [5964] and matrix method [1518]. Additionally, the temporal evolutions of the perturbation are investigated using the finite difference method [20, 65]. We also study the dependence of the resulting quasinormal frequencies on the spacetime configurations. Moreover, we use the Prony method [66, 67] to extract the complex frequencies and compare them to those obtained using the WKB approximation and matrix method. The matrix method discretizes the master equation and approximates it using a matrix equation that can eventually be solved by a nonlinear equation solver. Of course, the QNMs can also be investigated using other well-known approaches, such as the continued fraction method [13]. Specifically, the parameterization given by Eq. (17) is analytic; therefore, we expect that we can derive a three-term recurrence relation and apply the method. However, for this study, we focus on the three methods discussed above. Regarding the QNM frequencies, we are interested in the role of the magnetic field. In addition to the choice of s=1 and Q=0 as discussed earlier, we consider a simplied form of the metric by assuming the following metric parameters: w=109, κ=1, Ca=9, and λ=16. The resultant metric is expressed as

      f=12Mr+310r3.

      (17)

      The resultant quasinormal frequencies are presented in Tables 12 and Figs. 14. From Tables 12, we first conclude that the results obtained using the three different approaches are in reasonable agreement. In particular, the results of the matrix method and sixth-order WKB approximation are mostly consistent. We observe that the term meff=mqB functions as an effective "mass" of the initially massless scalar field. For positive effective mass meff, as the effective mass increases, the real part of the quasinormal frequency increases, whereas the magnitude of the imaginary part decreases. This indicates that the period of the dissipative oscillations decreases, whereas the amplitude decreases slower in time. At vanishing magnetic quantum number m=0, as expected, the quasinormal frequencies are independent of the strength of the magnetic field B.

      lmqB=0qB=0.05qB=0.1qB=0.15
      WKBMatrix methodWKBMatrix methodWKBMatrix methodWKBMatrix method
      110.298-0.096i0.298-0.096i0.277-0.108i0.281-0.109i0.256-0.120i0.266-0.127i0.236-0.130i0.200-0.178i
      00.298-0.096i0.298-0.096i0.298-0.096i0.298-0.096i0.298-0.096i0.298-0.096i0.298-0.096i0.298-0.096i
      −10.298-0.096i0.298-0.096i0.320-0.083i0.320-0.085i0.342-0.069i0.343-0.069i0.363-0.052i0.363-0.053i
      220.492-0.095i0.492-0.095i0.462-0.105i0.464-0.104i0.433-0.115i0.438-0.115i0.404-0.124i0.438-0.115i
      10.492-0.095i0.492-0.095i0.477-0.100i0.477-0.100i0.462-0.105i0.464-0.104i0.448-0.111i0.451-0.108i
      00.492-0.095i0.492-0.095i0.492-0.095i0.492-0.095i0.492-0.095i0.492-0.095i0.492-0.095i0.492-0.095i
      −10.492-0.095i0.492-0.095i0.508-0.090i0.508-0.091i0.523-0.085i0.524-0.085i0.539-0.080i0.539-0.079i
      −20.492-0.095i0.492-0.095i0.523-0.085i0.524-0.085i0.555-0.074i0.555-0.073i0.587-0.062i0.587-0.063i

      Table 1.  QNMs of the massless charged scalar field in the magnetized black hole metric in Rastall gravity. The sixth order WKB approximation and matrix method have been employed for the calculations, by assuming the metric given by Eq. (17).

      lmqB=0qB=0.05qB=0.1qB=0.15
      110.293-0.098i0.272-0.113i0.247-0.108i0.229-0.206i
      00.293-0.098i0.293-0.098i0.293-0.098i0.293-0.098i
      −10.293-0.098i0.316-0.084i0.339-0.070i0.370-0.051i
      220.484-0.097i0.453-0.108i0.427-0.120i0.377-0.114i
      10.484-0.097i0.468-0.102i0.453-0.108i0.438-0.113i
      00.484-0.097i0.484-0.097i0.484-0.097i0.484-0.097i
      −10.484-0.097i0.500-0.091i0.516-0.086i0.532-0.080i
      −20.484-0.097i0.516-0.086i0.548-0.074i0.581-0.061i

      Table 2.  QNMs of the massless charged scalar field in the magnetized black hole metric in Rastall gravity obtained using the code of Prony's method developed in [66].

      Figure 1.  (color online) Calculated temporal evolutions of charged massless scalar perturbations for different values of the effective mass. The results are obtained using the finite difference method using the paramters l=1 and M=1.

      Figure 4.  (color online) Effective potentials Veff(r) for l=2 and different effective mass.

      These results are also in agreement with the temporal evolution calculated using the finite difference method shown in Fig. 1. In particular, as also shown in Fig. 2, as the value of meff increases, the imaginary part eventually becomes insignificant as the envelope of the oscillations tends to lie horizontally. It is depicted explicitly in Fig. 3, where we calculate the quasinormal frequencies as functions of the effective mass using the matrix method. As the effective increase, the real part of the quasinormal frequency increases, and the imaginary part decreases to approach zero. Subsequently, the resulting temporal evolutions are featured by the so-called quasi-resonance [68, 69]. When we further increase the effective mass until it exceeds a critical value, the quasinormal oscillation will disappear entirely. The latter is demonstrated in Fig. 3 as the quasinormal spectrum approaches the real axis as the magnetic field increases. As discussed in the next section, such a configuration becomes favorable for investigating the properties of the late-time tail in a more transparent fashion.

      Figure 2.  (color online) Calculated temporal evolutions of charged massless scalar perturbations, where instability is observed for a negative effective mass. The results are obtained using the finite difference method for the parameters l=1 and M=1.

      Figure 3.  (color online) QNMs of different effective masses calculated using the matrix method for the massless charged scalar field with l=1,2, m=1, and n=0.

      For negative meff<0, the WKB and matrix methods fail to produce meaningful results, and we resort to the finite difference method entirely. This is signaled by the fact that the corresponding effective potential has a region of negative values, which indicates possible instability. The form of the relevant effective potential, governed by Eq. (16), is shown in Fig. 4. Similar to Ref. [45], we observe that such a configuration indeed results in instability by exploring the temporal evolution of the initial perturbations. To further verify this instability, we use the finite difference method to determine the temporal evolution of the massless charged scalar perturbations. As demonstrated in Fig. 2, as we gradually decrease the value of the effective mass, the oscillation eventually becomes unstable. Moreover, note that there is a distinction between the quasinormal oscillation and the late-time tail in the limit when the magnetic field vanishes, which is elaborated in the next section.

    IV.   LATE-TIME TAILS OF CHARGED MASSLESS SCALAR PERTURBATIONS
    • In this section, we explore the late-time tails of the charged massless scalar perturbations using the analytic method proposed in Ref. [19] and the numerical integration.

      In terms of the Green's function, the wave function of the scalar field can be given by [70]

      ψ(r,t)=[G(r,r;t)ψt(r,0)+Gt(r,r;t)ψ(r,0)]dr,

      (18)

      where the initial conditions are encoded in ψ(r,0) and ψttψ(r,0), and the Green's function G(r,r;t) is defined by

      [2t22r2+V]G(r,r;t)=δ(t)δ(rr),

      (19)

      When ˜G(r,r;ω) is obtained, the evolution of an arbitrary initial perturbation can be obtained through the integration given by Eq. (18). Indeed, both the QNM and late-time tail can be understood in terms of the singularities and branch cuts of the Fourier transform of the Green's function.

      ˜G(r,r;ω)=+0G(r,r;t)eiωtdt,

      (20)

      whose inverse transform is given by

      G(r,r;t)=12π+ic+ic˜G(r,r;ω)eiωtdω,

      (21)

      where c is a positive constant. Therefore, the analysis of the properties of the Green's function becomes an essential task.

      It is well-known that the Green's function ˜G(r,r;ω) can be constructed using two independent solutions of the corresponding homogeneous equation

      (d2dr2+ω2V)˜ψi=0.

      (22)

      These two solutions, ˜ψ1(r,ω) and ˜ψ2(r,ω), satisfy the appropriate boundary condition at both spatial boundaries [71], namely,

      ˜ψ1(r,ω){eiωrrA(ω)eiωr+B(ω)eiωrr+,

      (23)

      ˜ψ2(r,ω){C(ω)eiωr+D(ω)eiωrre+iωrr+.

      (24)

      We can easily show that the following form satisfies Eq. (19):

      ˜G(r,r;ω)=1W(ω){˜ψ1(r,ω)˜ψ2(r,ω),r>r˜ψ1(r,ω)˜ψ2(r,ω),r<r,

      (25)

      where

      W(ω)˜ψ1(r,ω)˜ψ2,r(r,ω)˜ψ2(r,ω)˜ψ1,r(r,ω)

      (26)

      is the Wronskian of ˜ψ1(r,ω) and ˜ψ2(r,ω).

      Research shows [24] that the late-time tails of quasinormal oscillations can be attributed to the branch cut of the Green's function. For the massless scalar field, the branch cut lies on the negative part of the imaginary axis [21]. For the massive scalar field, Hod and Piran [27] were the first to indicate that the branch cut lies on the real axis between the two branch points governed by the mass. Using reasonable approximation, an appropriate estimation of the contribution from the branch cut subsequently yields the main characteristic of the late-time tail. From a physical perspective, because the low-frequency waveforms are more likely to be backscattered by the curvature at the asymptotic infinity, the late-time tail is primarily dominated by the low-frequency contributions [27]. Mathematically, this is because both types of late-time tails originate from the branch cut, which remains closer to the real axis than any QNM pole.

      To proceed, we introduce

      Rlm(r,ω)=(12Mr+r(2+a))1/2˜ψ.

      (27)

      By further referring to Refs. [27, 29, 72, 73], we assume that the initial perturbation and observer are both located away from the black hole. Subsequently, the master equation can be simplified by expanding Eq. (15) in Mr. Specifically, we ignore the terms of order crα with α>2 and obtain

      [d2dr2+ω2+Bqm+4Mω2+Bqm2Mrl(l+1)r2]˜ψ=0.

      (28)

      By further introducing

      ˜ψ=zl+1ez2Φ(z),z=2Bqmω2r2ϖr,λ=M(Bqm)ϖ2Mϖ,

      (29)

      we can express two relevant solutions for homogeneous equation Eq. (22) as follows

      ˜ψ1=AMλ,(l+12)(2ϖr),˜ψ2=BWλ,(l+12)(2ϖr),

      (30)

      where A=Cϖ(l+1) and B,C are some one-valued even functions of ϖ. Mk,m and Wk,m are the Whittaker functions, which can be expressed in terms of the solutions of Kummer's equation M(a,b,z) and U(a,b,z), namely, the confluent hypergeometric functions [74],

      Mλ,(l+12)(2ϖr)=eϖr(2ϖr)(l+12)+12M(l+1λ,2l+2,2ϖr),Wλ,(l+12)(2ϖr)=eϖr(2ϖr)(l+12)+12U(l+1λ,2l+2,2ϖr).

      (31)

      Note that the specific choice of C guarantees the proper in-going wave boundary condition at the horizon. Moreover, as discussed later, ˜ψ1 satisfies Eq. (38) as it does not contain any discontinuity when crossing the branch cut [29, 72, 73]. In contrast, ˜ψ2 yields the outgoing wave at spatial infinity that has two branch points and, subsequently, a branch cut joining them. As initially proposed in [27], we may conveniently place the branching cut on the real axis of the ω-plane, as shown in Fig. 5. In this context, the late-time tail is primarily governed by the properties of ˜ψ2. Note that it is mandatory to show that the branch cut is in the lower half of the complex plane. Moreover, for the result to be physically meaningful, the specific choice of the branch cut's location must also be irrelevant. A detailed account of this point was discussed recently in Ref. [30].

      Figure 5.  Relevant contour of the integration on the complex frequency plane of the Green's function. The inverse Fourier transform corresponds to the integral along the upper side of the real axis. It can be complemented by an integral along a large semicircle in the lower half of the complex plane, which is subsequently deformed to go around the quasinormal poles and brach cuts indicated in the plot.

      By using the explicit forms of Eqs. (30)–(31), we can evaluate the contribution from the difference between the two sides of the branch cut. Therefore, the inverse Fourier transform Eq. (21) yields

      GC(r,r;t)=12πBqmBqmF(ϖ)eiωtdω,

      (32)

      where

      F(ϖ)˜ψ1(r,ϖeiπ)˜ψ2(r,ϖeiπ)W(ϖeiπ)˜ψ1(r,ϖ)˜ψ2(r,ϖ)W(ϖ).

      (33)

      Regarding the integrand in Eq. (32), we first calculate the Wronskian in the denominator using the following relations [74]

      Wλ,l+12(2ϖr)=Γ(2l1)Γ(lλ)Mλ,l+12(2ϖr)+Γ(2l+1)Γ(l+1λ)Mλ,(l+12)(2ϖr).

      (34)

      and the non-vanishing Wronskian between the confluent hypergeometric functions

      W{Mλ,l+12(2ϖr),Mλ,(l+12)(2ϖr)}=(2l+1)(2ϖ).

      (35)

      We observe that

      W(ϖ)=AB(2l1)(2ϖ)Γ(2l+1)Γ(l+1λ).

      (36)

      To account for the discontinuity of the integrand across the branch cut, we use

      M±λ,l+12(eiπ2ϖr)=e(l+1)iπMλ,l+12(2ϖr),

      (37)

      to evaluate ˜ψ1,2(r,ϖeiπ) as follows:

      ˜ψ1(r,ϖeiπ)=˜ψ1(r,ϖ),

      (38)

      ˜ψ2(r,ϖeiπ)=B[e(l+1)iπΓ(2l1)Γ(l+λ)Mλ,l+12(2ϖr)+e(l)iπΓ(2l+1)Γ(l+1+λ)Mλ,(l+12)(2ϖr)].

      (39)

      In the above results, note that the subscripts λ of the Whittaker functions remain unchanged. Eq. (38) manifestly indicates that ˜ψ1 does not have any branch cut as it is also a generic physical requirement [21, 27, 29]. Additionally, the Wronskian yields

      W(ϖeiπ)=AB(2l1)(2ϖ)Γ(2l+1)Γ(l+1+λ)e(l)iπ

      (40)

      We note that in deriving the above result, we should always take note of Eqs. (38) and (39). Specifically, the discontinuity associated with the factor ϖ(l+1) embeded in A does not manifest itself, whereas the factor eliπ on the r.h.s. of Eq. (39) does.

      By consolidating all the pieces, i.e., Eqs. (38), (39), and (40), we observe that the integrand of Eq. (32) yields

      F(ϖ)=1(2l1)(2ϖ)A2˜ψ1(r,ϖ)˜F(ϖ)˜ψ1(r,ϖ)

      (41)

      where

      ˜F(ϖ)=e(2l+1)iπΓ(2l1)Γ(l+1+λ)Γ(l+λ)Γ(2l+1)Γ(2l1)Γ(l+1λ)Γ(lλ)Γ(2l+1).

      (42)

      Here, owing to Eq. (38), ˜ψ1(r,ϖ) can be readily factorized out. The remaining term, ˜ψ2(r,ϖ), can be simplified by removing irrelevant contributions that have the identical discontinuity of the Wronskian. As a result, a second ˜ψ1(r,ϖ) can be removed.

      For the intermediate-late time scale, which corresponds to the range MrtM/(qBmM2), the frequency ϖ=O(qBm/t) results in a major contribution to the integral, which implies λ1 owing to Eq. (29). In this case, the term proportional to 1r in Eq. (28), which describes the backscattering due to the curvature of asymptotic infinity, is insignificant [28, 29]. Subsequently, the inverse Fourier transform of the Green's function at the limit λ1 can be performed mostly analytically, which results in

      GC(r,r;t)=(1e(2l+1)iπ)Γ(2l1)Γ(l+1)Γ(l+32)(qBm)2l+14(rr)l+1t2l+32π(4l+2)23l52Γ(2l+1)Γ(l)×cos[qBmtπ(2l+34)].

      (43)

      In other words, we observe that at the intermediate-late time, the tail is dominated by the form

      tl32cos[qBmtπ(l2+34)].

      It oscillates while decaying through a power-law form. In addition to the magnetic field, the intermediate-late time tail depends on the multipole number l and azimuthal number m. In particular, the presence of the magnetic field results in an effective mass that significantly modifies the properties of the tail. However, if we remove the magnetic field, the tail governed by taking the limit B0 in Eq. (43) does not simply fall back to its massless counterpart. Specifically, while the synodal oscillation is suppressed, the remaining power-law exponent becomes el32, which differs from the late-time tail e2l3 observed for a massless scalar in the Schwarzschild background. Moreover, as discussed later, at the limit of the vanishing magnetic field, the power-law exponent of the late-time tail does not match that of its massless counterpart. We postpone further discussions of this intriguing characteristic to the end of this section.

      Now, we discuss the asymptotic behavior of the late-time tails. At a significant time scale qBmt1/(qBmM)2, the relevant terms in the master equation behave differently from the case of intermediate-late time tail. The backscattering from the spacetime curvature at asymptotic infinity can no longer be ignored for this scenario. Thus, because both the ˜ψ1 factors in Eq. (41) do not depend sensitively on the frequency ω, we only require to focus on the term Eq. (42). By applying the above approximation and using the asymptotical forms of Γ functions [74], we obtain

      ˜F(ϖ)Γ(2l1)Γ(2l+1)λ2l+1[e(2l+1)iπη+eiπλ+ηeiπληeiπλ+η+eiπλ],

      (44)

      where

      η±=e±iπl.

      (45)

      The resulting integral of the inverse Fourier transform is somewhat complicated, but its asymptotic properties are governed by the terms related to λ. At small frequencies, we obtain ϖ0 and λ. Their product, which appears in the argument of the Bessel functions, yields

      λϖ=(M(Bqm)ϖ2Mϖ)ϖ=M(Bqm)+O(ϖ2),

      which is finite and varies gradually in the low-frequency region ωBqm+ [72]. Both the terms λ2l+1 and e±iπλ of Eq. (44) oscillate significantly as λ, whereas the latter type is even more dramatic compared with the former. This indicates that the integral's main contribution results from these terms. In most cases, we may follow the arguments of Koyama and Tomimatsu [28] that the main contribution of the integral can be obtained by employing the method of steepest descent. The strategy of such a treatment involves isolating the part that oscillates significantly and then encountering the saddle point, where the oscillations evolve the slowest. The remaining part might be tedious in form but varies moderately. We rewrite the integrand as

      ˜F(ϖ)eiωtΓ(2l1)Γ(2l+1)λ2l+1eiϕei(2πλωt),

      (46)

      where the phase can be defined by

      eiϕ=η++ηe2iπλη++ηe2iπλ.

      (47)

      As shown in [30], the variation in the phase ϕ precisely cancels that in the term e2iπλ. In this case, the integral cannot be evaluated using the method of steepest descent. Nonetheless, we may still evaluate the most dominant term of the integral and obtain

      GC(r,r;t)t1sin(μt+φ),

      (48)

      where the phase shift φ contains some minor dependence on t. The temporal dependence of Eq. (48) implies that late-time tail evaluated regarding the contribution from the saddle point is t1. Note that the obtained exponential does not depend on either the angular momentum number or the strength of the magnetic field.

      We present the numerical results in Figs. 69. To confirm the above analytical results, we perform the numerical calculations using numerical integration. The resultant temporal evolutions of the intermediate-late time are presented in Figs. 67. In Fig. 6, we show the temporal evolutions at the intermediate-late time for different effective masses. For a given l, through Eq. (43), the period of the oscillations increases as the effective mass decreases. Moreover, Fig. 7 shows that when the effective mass is fixed, the oscillating frequency largely remains unchanged for different values of l. Furthermore, the attenuation becomes more significant as l increases. In Fig. 7, the time profiles obtained numerically are also compared against the analytic ones shown in dashed curves given by Eq. (43). The extracted exponent from the envelope of the oscillation yields tβ, where β=2.5087±0.0020 and β=3.5015±0.0020 with respect to l=1, l=2. Therefore, the asymptotic form t(l+3/2) for the intermediate-late time agrees well with the numerical results for different l values.

      Figure 6.  (color online) Calculated intermediate-late time tail in the Green's function of charged massless scalar perturbations in a magnetized black hole in Rastall gravity. The results are obtained using numerical integration for different effective masses with l=1,2.

      Figure 9.  (color online) Calculated quasinormal oscillations and late-time tails of charged massless scalar perturbations in a magnetized black hole in Rastall gravity. The calculations aim to explore the asymptotic behavior of the QNMs and tails as the effective mass gradually vanishes. The results are obtained using the parameters l=1 and M=1. Top: The quasinormal oscillations are shown to converge to the case of massless perturbations as the effective mass approaches zero. Bottom: The envelopes of the tails are compared with the analytic forms tl3/2 and t2l3, respectively. The results indicate that the tails' asymptotic behavior does not fall back to that of massless perturbations, as discussed in the text.

      Figure 7.  (color online) Calculated intermediate-late time tail in the Green's function of charged massless scalar perturbations in a magnetized black hole in Rastall gravity. The results are obtained using numerical integration for different effective masses with l=1,2. The envelopes of dissipative oscillations shown in dashed curves are compared with the analytic ones t(l+3/2) given in the text.

      The temporal evolutions of the late-time tails are presented in Fig. 8. The calculations are performed for different effective masses and angular momenta. Again, by comparing the two plots in Fig. 8, as the effective mass increases, we observe that the oscillation frequency also increases, but the attenuation rate remains the same. Moreover, as shown in the bottom plot of Fig. 8, neither the attenuation rate nor the oscillation period depends on the angular momentum l. We also show in dashed curves the nonlinear fitting to the envelope of the oscillations. The extracted exponents from the a power-law form tβ for the two cases are β=0.9978±0.0068 and β=0.9762±0.0082, respectively, which are reasonably consistent with the analytic result t1.

      Figure 8.  (color online) Calculated late-time tail in the Green's function of charged massless scalar perturbations in a magnetized black hole in Rastall gravity. The results are obtained using numerical integration for different effective masses with the parameters l=1,2 and M=1. The envelopes of dissipative oscillations shown in dashed curves are compared with the analytic ones t1 given in the text.

      Finally, we elaborate on the specific role of the magnetic field in this study. A distinct characteristic of the present spacetime configuration is that the effective mass is generated by the magnetic field; therefore, the effect on the late-time tail is dynamic. In particular, we are allowed to consider a process in which one gradually tunes the magnetic field to approach the limit B0. Based on the preceding discussions, although the quasinormal frequencies are observed to vary continuously, we conclude that the resultant tail would not converge to its massless counterpart. This phenomenon is manifestly shown in Fig. 9. The top plot of Fig. 9 shows that as the effective mass decreases, the dissipative oscillations converge to the case of massless perturbations. This is consistent with the numerical results obtained in Table 3, where the Prony method [66, 67] has again been employed to extract the complex frequencies from the time profiles. In contrast, as shown in the bottom plot of Fig. 9, the late-time tail of massless perturbations is distinct from the remaining ones. Specifically, the massless perturbation indicated by the orange curve does oscillate at a finite period, whereas the oscillation periods of massive perturbations are more significant (c.f., the green and dark red curves). Regarding the attenuation rate, all the curves of massive perturbations are essentially parallel, which differs from that of the massless one. Quantitatively, the values extracted from the numerical fitting deviate from one another, as indicated in the caption of the plot.

      meffPronyWKBMatrix method
      00.29861-0.09614i0.29842-0.09636i0.29845-0.09622i
      0.000010.29861-0.09614i0.29842-0.09636i0.29845-0.09622i
      0.00010.29865-0.09612i0.29846-0.09634i0.29849-0.09620i
      0.0010.29904-0.09590i0.29885-0.09611i0.29884-0.09600i
      0.010.30288-0.09421i0.30276-0.09384i0.30246-0.09408i

      Table 3.  QNMs for the metric Eq. (17) are calculated by gradually reducing the effective mass. The Prony method, sixth order WKB approximations, and matrix method are used for the calculations. The results are obtained for M=1,l=1.

      We now argue that the case of massless perturbation corresponds to a scenario in which the magnetic field does not exist in the first place. Such an interesting difference is somewhat reminiscent of the demagnetization process in ferromagnetic materials. Mathematically, this can be understood as follows. The convergence of the quasinormal frequencies can be understood as the analytic properties of the effective potential essentially remaining unchanged as the effective mass varies continuously. However, the origin of a massive field's late-time tail is associated with the presence of two branch points on the real axis and the absence of any branch point at infinity. As the effective mass approaches zero, the two branch points become degenerate, but a new branch cut will not appear, which extends to infinity. Moreover, a single branch point will be located at the origin in the case of a massless field. The branch cut on the negative imaginary axis connects the origin with the second branch point at infinity. As a result, the relevant structure on the frequency plane due to a massless field is intrinsically different from that of a massive perturbation, although the mass is insignificant. In other words, if we gradually remove the magnetic field from a magnetized black hole, the perturbation will evolve differently from the spacetime configuration in which the magnetic field in question does not ever exist. We note that the above picture is plausible only when the time scale of the demagnetization process is significantly smaller than the inverse of the typical distance between adjacent quasinormal frequencies. This is because the resulting effective potential is now a function of time. Therefore, the notion that QNMs correspond to the zeros of the Wronskian and subsequently the poles of the frequency domain Green's function [3] is only approximately valid when the process can be considered quasistatic. In literature, the master equation with a time-dependent potential is mostly treated numerically [7577].

    V.   CONCLUDING REMARKS
    • In this study, we investigate the QNMs and late-time tails of charged massless scalar perturbations in a magnetized black hole in Rastall gravity. We show that the magnetic field has a significant role in both aspects of the resultant properties of the scalar perturbation. Specifically, the massless scalar acquires an effective mass through the magnetic field, which significantly affects the temporal evolution of the initial perturbations. For the quasinormal oscillations, the complex frequencies are distorted and might become either quasinormally resonant or unstable for specific parameters. For the intermediate and late-time tails, power-law forms are obtained analytically and numerically, reminiscent of those of massive scalar perturbations. Moreover, owing to the dynamic nature of the effective mass generated by the magnetic field, we argue that it has an interesting characteristic. In particular, as we gradually reduce the external magnetic field, while the quasinormal frequencies converge to its massless counterpart, the behavior of the late-time tail deviates from the latter. From a mathematical perspective, this can be understood by analyzing the structure of branch points and cuts on the complex frequency plane of the relevant Green's function.

    ACKNOWLEDGMENTS
    • We wish to thank Rui-Hong Yue for enlightening discussions.

    APPENDIX A: APPENDIX
    • In this Appendix, we show that when the magnetic field is small, the metric given by Eqs. (3)–(9) is a valid approximation. For simplicity, we assume that the strength of the magnetic field is insignificant and does not consider the "backreaction" to the metric. First, we readily verify that Eq. (9) satisfies the Maxwell equation for linear electromagnetic field in Schwarzschild's solution of Einstein gravity. For the Schwarzschild solution Eq. (8) of generalized Rastall gravity, we must ascertain that the deviation for Eq. (9) as linear electromagnetic field (s=1) is of higher-order, consistent with the approximation performed later. To show this, we rewrite the deviation from Eq. (17) as

      ˜f=f+δf,

      and solve the Maxwell equation for δf.

      Evidently, δf satisfies

      2δf+2rδf=65r3,

      which implies the magnitude of the correction δfO(r3), irrelevant to the effective potential in the master equation Eq. (28).

Reference (77)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return