AbstractThis study set out to investigate charged vector particles tunneling via horizons of a pair of accelerating rotating charged NUT black holes under the influence of quantum gravity. To this end, we use the modified Proca equation incorporating generalized uncertainty principle. Applying the WKB approximation to the field equation, we obtain a modified tunneling rate and the corresponding corrected Hawking temperature for this black hole. Moreover, we analyze the graphical behavior of the corrected Hawking temperature with respect to the event horizon for the given black hole. By considering quantum gravitational effects on Hawking temperatures, we discuss the stability analysis of this black hole. For a pair of black holes, the temperature increases with the increase in rotation parameters a and , correction parameter , black hole acceleration , and arbitrary parameter k, and decreases with the increase in electric e and magnetic charges g.
A black hole (BH) is a physical object that has the ability to absorb all types of energy from its surroundings as a result of its strong gravitational pull. According to the theory of general relativity, a BH attracts all types of particles that interact with the event horizon. Hawking (1974) described that a BH acts as a black body and emits particles in the form of radiation through its horizon Considering quantum effects in the background of curved spacetime, this radiation is known as Hawking radiation  and exhibits a particular temperature, called the Hawking temperature. Quantum tunneling proved to be one of the best techniques [2-4] employed to study the Hawking radiation phenomenon. This technique depends on electron-positron pair-creation, which requires an electric field. As particles overcome the horizon, their energy changes the sign, such that a member of the pair created just inside/outside the horizon can materialize with zero total energy, while the other has the ability to tunnel to the opposite side of the horizon .
According to the tunneling phenomenon, particles are permitted to follow classically forbidden trajectories, by starting from outside the horizon to infinity [6, 7]. BH evaporation can be studied with the discharge of quantum particles in form of Hawking radiation, which enables a BH to lose its mass. When a BH loses more matter than it gains through different means it disseminates, shrinks, and eventually disappears. This phenomenon causes a change in the thermodynamical characteristic of a BH, i.e., charge, mass, and angular momentum. These particles follow the frame of outgoing/ingoing radial null geodesics. For outgoing geodesics, the particles must be imaginary, although for ingoing geodesics, they are considered to be real. This is because only a real particle that has a speed, which is not exactly or equivalent to the speed of light can exist inside the horizon.
To calculate the imaginary part of the classical action, there are two main approaches, i.e., null geodesic technique and Hamilton-Jacobi strategy. The first one proposed by Parikh and his colleagues  and second one introduced by Srinivasan and Padmanabhan . To construct a relation between classical and quantum theory, an approximation is known as WKB approximation utilized by Wentzel, Kramers, and Brillouin (WKB) . Various authors studied the tunneling of vector particles, bosonic particles, spin-2, spin-3/2, and fermionic particles to obtain the Hawking temperature for different BHs and wormholes [9-59].
Sign up for Chinese Physics C TOC Get the latest issue of Advanced Photonics delivered right to you！Sign up now
By incorporating generalized uncertainty principle (GUP) effects, it is conceivable to discuss quantum-corrected thermodynamical properties of BH . The GUP provides high-energy remedies to BH thermodynamics, which leads to the possibility of a minimal length in quantum gravity theory. According to the modified fundamental commutation relation
[61, 62], the GUP expression can be given as
denotes the dimensionless parameter, and
depicts the Plank mass. The
stand for modified position and momentum operators, respectively. They can be defined as
are standard position and momentum operators, respectively, which satisfy the usually commutation relation
. We have considered the values of
according to this condition, which also satisfies the condition of GUP relation . For corrections of the Hawking temperature, we have considered only the first order terms of
in our calculation. The idea of GUP has been utilized for various BHs. The tunneling process for Kerr, Kerr-Newman, and Reissner-Nordström (RN) BHs  provides an important contribution towards the BH physics. Jiang  calculated the tunneling of Dirac particles and analyzed the Hawking radiation spectrum for a black ring. Jian and Bing-Bing  studied the fermion particles tunneling from uncharged and charged BHs. Yale  studied the scalar, fermion, and boson particles tunneling from BHs without back-reaction effects, they also analyzed the exact Hawking temperature for these particles. Sharif and Javed  studied the Hawking radiation through the quantum tunneling process for various types of BHs and also derived the tunneling probabilities as well as their corresponding Hawking temperatures. The same authors  calculated the tunneling behavior of fermion particles from traversable wormholes. They also discussed the tunneling behavior of fermion particles for charged accelerating rotating BHs associated with the NUT parameter .
n et al.  calculated the tunneling rate of charged massive bosons for various types of BHs surrounded by the perfect fluid in Rastall theory. Anacleto et al.  studied the finite quantum-corrected entropy from non-commutative acoustic BHs. Li and Zu  analyzed GUP effects by applying the Klein-Gordon equation and also studied corrected temperature for Gibbon-Maeda-Dilation BH. Anacleto et al.  investigated the quantum-corrected entropy using a tunneling method with GUP in self-dual BH. Övgün and Jusufi  studied spin-
vector particles tunneling for charged non-commutative BHs and their GUP-corrected thermodynamical properties. Singh et al.  analyzed the particle radiation for Kerr-Newman BH using the quantum tunneling phenomenon. Sakalli et al.  investigated the scalar particles tunneling from acoustic BHs with a rotation parameter by applying the WKB approximation and Klein-Gordon equation with GUP to obtain the Hawking temperature of massive particles.
Javed et al.  analyzed the charged vector particles tunneling phenomenon for a pair of accelerating and rotating, as well as 5D gauged super-gravity BHs and investigated the corresponding Hawking temperatures. As a continuation of their work, we study the Hawking radiation phenomenon by considering quantum corrections of boson particles through the event horizon of accelerating and rotating BH with the NUT parameter. We also discuss the graphical behavior of the corrected temperature with respect to the horizon for a given BH and analyzed the effects of different parameters on temperature and visualize the stable and unstable states of the BH. The main motivation of this work is to study the charged accelerating and rotating black hole with NUT parameters under the influence of quantum gravity, and to investigate how the results reduce to the previously obtained results when we neglect quantum gravity effects and other parameters. The organization of this study is as follows: In section 2, we study the quantum-corrected tunneling rate and correct the Hawking temperature for charged BH solution with acceleration, rotation, and NUT parameter. Moreover, we analyze the graphical behavior of corrected temperature with respect to the event horizon of BH and discuss the effects of different parameters on the temperature. Section 3 is based on the conclusion and discussion of the mathematical results for these BHs.
2. Pair of charged accelerated black holes involving rotation and nut parameters
Universally, the NUT parameter is associated with the twisting behavior of the surrounding spacetime or with the gravito-magnetic monopole parameter of the central mass. The precise physical significance of the NUT parameter might not be investigated. In contrast, the higher dimensional origin of the Kerr-NUT-(anti) de-Sitter BH and its physical significance was studied [79, 80]. For a BH, the apparent dominance of the NUT parameter on the rotation parameter removes the metric free of curvature singularity, and the corresponding result is known as the NUT-like result. When the rotation parameter commands the NUT parameter, the solution is Kerr-like, and a ring curvature singularity occurs. This type of singularity is independent from the existence of the cosmological constant.
The exact understanding of the NUT parameter is conceivable when a static Schwarzschild mass is submerged in a stationary source-free electromagnetic universe . Moreover, the NUT parameter is associated with the representation of twisting behavior of the electromagnetic universe, allowing off-Schwarzschild central mass. When electric and magnetic fields vanish, the NUT parameter represents the twist of the vacuum spacetime . Consequently, for the creation of the NUT parameter, the bending of the environmental space couples with the mass of the background field.
The line-element of this BH can be defined as follows 
represent mass, magnetic charge, and electric charge, respectively;
are the NUT parameter, rotation, and acceleration of source, respectively. Moreover,
is rotation parameter (parameter of Kerr-like), and
is defined as
We observe that
. The parameter
represents the bending behavior of BHs relative to the rotation parameter (
). The parameters
are independent. As
, the line-element (2) implies the Kerr-Newman NUT results. For
, the line-element (2) yields the BHs with rotating and charged parameters. When
approach zero, we obtain a Schwarzschild BH, and if
approach zero, the C-metric is implied.
The line-element (2) can be rewritten in the following form
where the metric functions Z, B, C, D, and F are given as follows
The electromagnetic potential for this pair of BHs is defined as
The horizons are obtained from
, which implies that
is not zero, thus
is zero, such that r has real roots, i.e.,
represent the acceleration horizons, and
denotes the inner and outer horizons, such that
The angular velocity of BH at outer horizon is given as follows
To study the corrected tunneling rate for vector bosons through the pair of BHs horizon, we consider the Lagrangian equation with quantum gravity effects. Considering a spacetime with a potential due to both the electric and magnetic field, the behavior of massive (spin-1) vector is described by the wave equation. The GUP modified Lagrangian equation with vector field
for vector bosons can be expressed as [56, 83]
The modified wave equation for massive vector bosons can be defined as follows
are the matrix coefficients, mass of a particle, and anti-symmetric tensor respectively. The anti-symmetric tensor can be defined as
is the representation of the quantum gravity effect in terms of the correction parameter, if
then the above wave equation reduces to the expressions given in Ref. , and if
then the above wave equation converts into the expressions given in Ref. . Moreover, if both
are zero, then this wave equation reduces to the wave equation given in Ref. . While
is considered as the BH electromagnetic potential,
denotes the vector bosons charge, and
represents the covariant derivative. In the equation of wave motion, the
vector bosons have alike behavior, and the quantum tunneling phenomena for both types of particles is also similar. For the
field, the values of
are obtained as follows
Applying the WKB approximation 
) in the Lagrangian wave of Eq. (18), ignoring the higher terms, the following set of equations is obtained
Using the variables separation technique, the particle's action is defined as follows
where J and E denote the angular momentum and particles energy, respectively, and from the above Eqs. (26)–(29), we obtain a non-trivial matrix equation, i.e.,
matrix, whose elements are given below
. For the non-trivial solution, we use
, and by computing the radial part for resultant components, the following integral can be obtained as follows
denote the radial function of incoming and outgoing vector bosons, respectively. Although the function
can be determined as
After applying Taylor's series, the functions
near the horizon can be obtained as
Using the above relations in Eq. (48), we consider that the leading wave equation has two poles at
. Utilizing Eqs. (48) and (51), by integrating around the pole, we get
The surface gravity
is given as follows 
The corrected tunneling probability
for vector bosons can be obtained as
We calculate the corrected Hawking temperature by comparing the tunneling probability with the Boltzmann factor, i.e.,
by considering only first order quantum corrections, we can write
where the semi-classical Hawking temperature
is given as follows
The corrected tunneling probability depends on
, which are the vector potential of BHs, energy of particle, NUT parameter, kerr-like rotation parameter, charges of particles, acceleration of BHs, angular momentum of particle, and correction parameter, respectively. We can observe that the corrected Hawking temperature does not only depend upon the BH properties but also depends upon the mass and angular momentum of the radiated particles and quantum corrections
. The first order correction term is same as semi-classical original Hawking term
, while the next order correction term must be smaller than the preceding term satisfying GUP.
The corrected temperature of vector bosons given in Eq. (56) reduced to the temperature of fermion particles in Eq. (4.20) for
of Ref. . Hence, the corrected temperature depends on the correction parameter. Further, for
, the above result reduces to the Hawking temperature of accelerating and rotating BHs with electric and magnetic charges . For
, we recover the temperature for non-accelerating BHs . Moreover, for
in Eq. (56), the Hawking temperature of the Kerr-Newman BH  is recovered, and it is reduced for
to the temperature of ReissnerNordstr
m BH, and for
, the temperature reduces exactly to the Hawking temperature of the Schwarzschild BH . To calculate the residual mass of BH, the temperature can be expressed as
denotes the kinetic energy component of radiated particles related with the tangent plane at horizon. For residual mass, we approximate the kinetic energy component as
. Quantum corrections decelerate the increase in temperature during the radiation process. These corrections cause the radiation ceased at some specific temperature, leaving the remnant mass. The temperature stops increasing when this condition holds 
is the Planck mass, and
is a dimensionless parameter representing quantum gravity effects where
[89, 90], one can obtain the following constraints
Here, it is important to mention that the value of the corrected Hawking temperature is smaller than the original temperature, and BH stops radiating when the mass of the BH reaches its minimal value
2.1. Graphical analysis of
This subsection is devoted to the study of the graphical behavior of Hawking temperature
with respect to the horizon
under the influence of quantum gravity. Furthermore, we study the physical significance of these graphs under the influence of the correction parameter
, NUT parameter
, rotation parameters
, BH acceleration
, electric and magnetic charges
, arbitrary parameter
, for fixed BH mass
and arbitrary parameter
. We analyze stability and instability of accelerating and rotating BH associated with the NUT parameter.
Figure 1 indicates the behavior of
with respect to
Figure 1.(color online)
(i) An exponential increase in the temperature
is observed, and after attaining a height,
slightly decreases with the increasing horizon
for varying values of the rotation parameter
decreases as the horizon increases, and this physical behavior reflects the BH stability with positive temperature untill
. The temperature increases with the increase in
(ii) We observe the behavior of the temperature for varying values of the correction parameter
with fixed values of the other parameters. The
decreases with increasing horizon in the domain
after attaining a maximum height. The maximum temperature with non-zero horizon reflects the BH remnant. This physical behavior of
indicates the BH stability with a positive range. It is also worth notin that
increases with the increase in correction parameter
Figure 2 shows the behavior of
Figure 2.(color online)
(i) exponentially increases and slowly drops from a height for different values of electric and magnetic charges. The decrease in
with the increasing horizon exhibits the stable state of BH in the domain
decreases with the increase of BH electric and magnetic charges
(ii) The temperature eventually drops down from a height for different values of the arbitrary parameter
. There is a significant change in temperature as it decreases exponentially and attains an asymptotically flat shape, which indicates BH stability untill
increases with an increase of
Figure 3 depicts the behavior of
in the domain
Figure 3.(color online)
(i) Initially, the temperature increases with increasing horizon, and after a maximum height it exponentially decreases, which indicates the stable state of BH for different values of BH acceleration
(ii) There can be seen that
exponentially increases and eventually drops from a height and decreases as the horizon increases untill
. This physical behavior indicates the stability of BH in positive ranges of different values of
. We observe that with the increase in the value of
3. Conclusion and discussion
In this study, we have investigated the quantum gravity effects for vector bosons from charged accelerating rotating BH with the NUT parameter. By considering the GUP effects, we employed the modified Lagrangian equation incorporating quantum effects describing the motion of spin-1 particles. Subsequently, by applying the Hamilton-Jacobi technique, we have calculated the tunneling probabilities of vector bosons. Moreover, we analyzed the corrected Hawking temperatures of these BHs. We concluded that the modified tunneling probabilities are not just dependent on the BHs properties, but also on the properties of emitted vector bosons, i.e., energy, potential, surface gravity, particle charge, and total angular momentum. Moreover, it is important to note that the modified tunneling probabilities as well as the Hawking temperature depend on the quantum particles, which contribute gravitational radiation in form of massive particle (BH's energy carrier) tunneling.
When the quantum gravity effects are neglected, i.e.,
, then the corrected Hawking temperature Eq. (56) is reduced to the absolute temperature obtained by the quantum tunneling of vector bosons . If we ignore the potential effects
, the modified Hawking temperature is reduced to the temperature of vector (spin-1) particles provided in Refs. [69, 91]. Further, for
, the above results are reduced to the Hawking temperature of the accelerating and rotating BHs with electric and magnetic charges . For
, we recovered the temperature of non-accelerating BHs from the Hawking temperature of the charged accelerating and rotating BHs . Moreover, for
in Eq. (56), the Hawking temperature of the Kerr-Newman BH  is recovered, which is reduced for
to the temperature of the RN BH. For
, the temperature reduces exactly to the Hawking temperature of the Schwarzschild BH at the residual mass of the BH .
In our analysis, we have found that the quantum corrections decelerate the increase in temperature during the radiation process. This correction causes the radiation to cease at some specific temperature, leaving the remnant mass. The remnant mass is obtained at the specific condition
Here, it is important to mention that the value of the corrected Hawking temperature is smaller than the original Hawking temperature, and the BH stops radiating, when its mass reaches the minimal value
. This result remains valid if the background BH geometry is more generalized.
The results from the graphical analysis of corrected Hawking temperatures with respect to the horizon for the given BH are summarized as follows:
● For accelerating and rotating BH with NUT parameter, the
decreases with the increasing horizon, and BH reflects the stable state for varying values of the rotation parameter
and correction parameter
. The corrected temperature
also increases with the increase in
. The BH remnant can be obtained at nonzero horizon with maximum temperature for different values of
in the domain
● Corrected temperature
decreases with the increase of
. For different values of electric and magnetic charges, the BH exhibits stability in the domain
increases with the increase in the value of arbitrary parameter
increases with the increase in BH acceleration
and rotation parameter
● In our analysis, we considered the value
, then the condition of GUP must be satisfied for arbitrary values of
, the correction term is smaller than the usual term, and the positive temperature is obtained. For
, the first order correction term becomes greater than the usual term, and the condition of GUP is not satisfied. Thus, we observe the negative temperature, which is non-physical. Furthermore, for
, the semi-classical term cancels out with the first order correction term, and hence the temperature vanishes.
 . Nature, 248, 30(1974).
 . Class. Quant. Grav., 25, 095014(2008).
 . Phys. Rev. D, 73, 064003(2006).
 . Phys. Rev. D, 98, 104001(2018).
 . Phys. Rev. Lett., 85, 5042(2000).
 . Phys. Rev. D, 60, 024007(1999).
 . J. High Energy Phys., 06, 095(2008).
 . Comptes Rendus de lAcademie des Sciences, 183, 24(1926).
 . Gen. Rel. Grav., 48, 1(2016).
 . Eur. Phys. J. Plus, 130, 110(2015).
 . EPL, 110, 10008(2015).
 . EPL, 118, 60006(2017).
 . Astrophys. Space Sci., 359, 32(2015).
 . J. Exp. Theor. Phys., 121, 404(2015).
 . Eur. Phys. J. Plus, 131, 184(2016).
 . Int. J. Theor. Phys., 55, 2919(2016).
 . Eur. Phys. J. Plus, 132, 298(2017).
 . Adv. High Energy Phys., 2017, 8798657(2017).
 . Eur. Phys. J. C, 77, 613(2017).
 . Eur. Phys. J. C, 78, 840(2018).
 . Mod. Phys. Lett. A, 34, 1950057(2019).
 . Gen. Rel. Grav., 50, 62(2018).
 . Adv. High Energy Phys., 2017, 1573904(2017).
 . Int. J. Theor. Phys., 56, 1725(2017).
 . Int. J. Mod. Phys. A, 29, 1430054(2014).
 . Eur. Phys. J. C, 74, 2687(2014).
 . Mod. Phys. Lett. A, 33, 1850164(2018).
 . Phys. Lett. B, 773, 391(2017).
 . Adv. High Energy Phys., 2018, 8728564(2018).
 . Preprints, 2019030194(2019).
 . Int. J. Theor. Phys., 56, 2279-2292(2017).
 . Class. Quant. Grav., 24, 5849-5860(2007).
 . Int. J. Theor. Phys., 47, 2203-2217(2008).
 . Gen. Rel. Grav., 36, 1171-1181(2004).
 . Int. J. Theor. Phys., 48, 3517(2009).
 . Phys. Lett. B, 642, 124(2006).
 . Theor. Math. Phys., 163, 774(2010).
 . Am. J. Phys., 78, 685(2010).
 . Int. J. Mod. Phys. D, 19, 159(2010).
 . Phys. Lett. B, 666, 269(2008).
 . Phys. Lett. B, 673, 227(2009).
 . Int. J. Mod. Phys. A, 22, 1705(2007).
 . Gen. Rel. Grav., 50, 152(2018).
 . Int. J. Mod. Phys. A, 33, 1850070(2018).
 . Int. J. Mod. Phys. D, 24, 1550005(2014).
 . Chin. Phys. B, 18, 462(2009).
 . Phys. Lett. B, 671, 167(2009).
 . Phys. Lett. B, 661, 370(2008).
 . EPL, 79, 20001(2007).
 . Phys. Lett. B, 651, 58(2007).
 . Phys. Lett. B, 653, 330(2007).
 . Chin. Phys. Lett., 24, 3336(2007).
 . Phys. Lett. B, 641, 139(2006).
 . Phys. Lett. B, 751, 34(2015).
 . Int. J. Mod. Phys. Rev. A, 30, 1550083(2015).
 . Astrophys. space Sci., 48, 359(2015).
 . Astrophys. Space Sci., 361, 207(2016).
 . J. Cosmol. Astropart. Phys, 03, 036(2014).
 . JHEP, 1008, 043(2010).
 . Adv. High Energy Phys., 2013, 432412(2013).
 . Phys. Rev. D, 65, 104022(2002).
 . Phys. Rev. D, 78, 044009(2008).
 . Acta Physica Polonica, 40, 241(2009).
 . Phys. Lett. B, 697, 398(2011).
 . Can. J. Phys., 90, 903(2012).
 . Can. J. Phys., 91, 43(2013).
 . Eur. Phys. J. C, 72, 1997(2012).
 . Advances in High Energy Physics, 3131620, 11(2018).
 . Annals of Phys., 362, 436(2015).
 . Journal of Applied Mathematics and Physics, 3, 134(2015).
 . Phys. Lett. B, 749, 81(2015).
 . Eur. Phys. J. Plus, 131, 177(2016).
 . Astrophys Space Sci, 103, 361(2016).
 . Astrophys Space Sci, 361, 330(2016).
 . Eur. Phys. J. C, 77, 296(2017).
 . Class. Quantum Gravity, 23, 5323(2006).
 . Class. Quant. Gravity, 24, 1687(2007).
 . Gen. Relativ. Gravity, 38, 1729(2006).
 . Class. Quant. Gravity, 22, 3467(2005).
 . Phys. Lett. B, 763, 80(2016).
 . Eur. J. Phys. Education, 2, 1309(2011).
 . Phys. Lett. B, 699, 15(2011).
 . Phys. Lett. B, 665, 277(2008).
 . Phys. Lett. B, 665, 106(2008).
 . JCAP, 1403, 036(2014).
 . Phys. Rev. Lett., 109, 151602(2012).
 . JHEP, 04, 024(1999).