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.

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 .

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].

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

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

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

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.