1. Introduction

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 [1] 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 [5].

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 [5] and second one introduced by Srinivasan and Padmanabhan [6]. To construct a relation between classical and quantum theory, an approximation is known as WKB approximation utilized by Wentzel, Kramers, and Brillouin (WKB) [8]. 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].

By incorporating generalized uncertainty principle (GUP) effects, it is conceivable to discuss quantum-corrected thermodynamical properties of BH [60]. 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

where
,
denotes the dimensionless parameter, and
, the
depicts the Plank mass. The
and
stand for modified position and momentum operators, respectively. They can be defined as
where
and
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 [63]. 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 [64] provides an important contribution towards the BH physics. Jiang [65] calculated the tunneling of Dirac particles and analyzed the Hawking radiation spectrum for a black ring. Jian and Bing-Bing [66] studied the fermion particles tunneling from uncharged and charged BHs. Yale [67] 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 [68] 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 [69] 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 [70].

vg
n et al. [71] calculated the tunneling rate of charged massive bosons for various types of BHs surrounded by the perfect fluid in Rastall theory. Anacleto et al. [72] studied the finite quantum-corrected entropy from non-commutative acoustic BHs. Li and Zu [73] analyzed GUP effects by applying the Klein-Gordon equation and also studied corrected temperature for Gibbon-Maeda-Dilation BH. Anacleto et al. [74] investigated the quantum-corrected entropy using a tunneling method with GUP in self-dual BH. Övgün and Jusufi [75] studied spin-
vector particles tunneling for charged non-commutative BHs and their GUP-corrected thermodynamical properties. Singh et al. [76] analyzed the particle radiation for Kerr-Newman BH using the quantum tunneling phenomenon. Sakalli et al. [77] 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. [78] 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 [81]. 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 [70]. 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 [82]

where

where
,
, and
represent mass, magnetic charge, and electric charge, respectively;
,
, and
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
depends on
and
. The parameter
represents the bending behavior of BHs relative to the rotation parameter (
). The parameters
,
,
,
,
, and
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
and
approach zero, we obtain a Schwarzschild BH, and if
and
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.,

Here,
and
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

where
,
, and
are the matrix coefficients, mass of a particle, and anti-symmetric tensor respectively. The anti-symmetric tensor can be defined as

and

where
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. [78], and if
then the above wave equation converts into the expressions given in Ref. [71]. Moreover, if both
and
are zero, then this wave equation reduces to the wave equation given in Ref. [76]. 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
and
vector bosons have alike behavior, and the quantum tunneling phenomena for both types of particles is also similar. For the
field, the values of
and
are obtained as follows

Applying the WKB approximation [84]

where

and (for
) 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.,

where
is a
matrix, whose elements are given below

where
,
, and
. For the non-trivial solution, we use
, and by computing the radial part for resultant components, the following integral can be obtained as follows

where
and
denote the radial function of incoming and outgoing vector bosons, respectively. Although the function
and
can be determined as

and

After applying Taylor's series, the functions
and
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

where
.

The surface gravity
is given as follows [70]

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
, *E*,
,
,
,
,
, and
, 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. [70]. Hence, the corrected temperature depends on the correction parameter. Further, for
and
, the above result reduces to the Hawking temperature of accelerating and rotating BHs with electric and magnetic charges [85]. For
, we recover the temperature for non-accelerating BHs [86]. Moreover, for
,
,
, and
in Eq. (56), the Hawking temperature of the Kerr-Newman BH [87] 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 [88]. To calculate the residual mass of BH, the temperature can be expressed as

where
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 [63]

For
,
, and
where
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
versus

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
and
, BH acceleration
, electric and magnetic charges
and
, 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
for fixed
,
and
[89, 90].

Figure 1.(color online)
vs.
for
,
&
.

(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
. The
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
w.r.t
for fixed
,
&
.

Figure 2.(color online)
vs.
for
,
, &

(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
and
, respectively.

(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
. The
increases with an increase of
.

Figure 3 depicts the behavior of
w.r.t
for fixed
,
&
for varying
and
in the domain
.

Figure 3.(color online)
vs.
for
,
&

(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
,
increases.

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 [78]. 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
and
, the above results are reduced to the Hawking temperature of the accelerating and rotating BHs with electric and magnetic charges [85]. For
, we recovered the temperature of non-accelerating BHs from the Hawking temperature of the charged accelerating and rotating BHs [86]. Moreover, for
,
,
and
in Eq. (56), the Hawking temperature of the Kerr-Newman BH [87] 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 [88].

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
and
. 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
and
. For different values of electric and magnetic charges, the BH exhibits stability in the domain
.

● The
increases with the increase in the value of arbitrary parameter
.

● The
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.