GUP Modified Hawking Radiation in Bumblebee Gravity
Sara Kanzi, \.Izzet Sakall{\i}

TL;DR
This paper investigates how Lorentz symmetry breaking and the generalized uncertainty principle influence Hawking radiation and greybody factors of a Schwarzschild-like black hole in bumblebee gravity, revealing modifications to temperature and radiation properties.
Contribution
It introduces a novel analysis of Hawking radiation under GUP and LSB effects in bumblebee gravity, providing semi-analytic greybody factors and demonstrating the impact of LSB on radiation characteristics.
Findings
GUP alters the Hawking temperature of the black hole.
LSB decreases the greybody factor, affecting radiation emission.
The effects are consistent across different particle spins.
Abstract
The effect of Lorentz symmetry breaking (LSB) on the Hawking radiation of Schwarzschild-like black hole found in the bumblebee gravity model (SBHBGM) is studied in the framework of quantum gravity. To this end, we consider Hawking radiation spin-0 (bosons) and spin- particles (fermions), which go in and out through the event horizon of the SBHBGM. We use the modified Klein-Gordon and Dirac equations, which are obtained from the generalized uncertainty principle (GUP) to show how Hawking radiation is affected by the GUP and LSB. In particular, we reveal that, independent of the spin of the emitted particle, GUP causes a change in the Hawking temperature of the SBHBGM. Furthermore, we compute the semi-analytic greybody factors (for both bosons and fermions) of the SBHBGM. Thus, we reveal that LSB is effective on the greybody factor of the SBHBGM such that its redundancy…
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
GUP Modified Hawking Radiation in Bumblebee Gravity
Sara Kanzi
İzzet Sakallı
Physics Department, Arts and Sciences Faculty, Eastern Mediterranean University, Famagusta, North Cyprus via Mersin 10, Turkey.
([)
Abstract
The effect of Lorentz symmetry breaking (LSB) on the Hawking radiation of Schwarzschild-like black hole found in the bumblebee gravity model (SBHBGM) is studied in the framework of quantum gravity. To this end, we consider Hawking radiation spin-0 (bosons) and spin- particles (fermions), which go in and out through the event horizon of the SBHBGM. We use the modified Klein-Gordon and Dirac equations, which are obtained from the generalized uncertainty principle (GUP) to show how Hawking radiation is affected by the GUP and LSB. In particular, we reveal that, independent of the spin of the emitted particle, GUP causes a change in the Hawking temperature of the SBHBGM. Furthermore, we compute the semi-analytic greybody factors (for both bosons and fermions) of the SBHBGM. Thus, we reveal that LSB is effective on the greybody factor of the SBHBGM such that its redundancy decreases the value of the greybody factor. Our findings are graphically depicted.
Hawking Radiation, Lorentz Symmetry Breaking, Bumblebee Gravity, GUP, Greybody, Quantum Gravity, Klein-Gordon Equation, Dirac Equation
††preprint: HEP/123-qed
09.03.2018]
Contents
I INTRODUCTION
In spite of their overwhelming successes in describing nature, General Relativity (GR) (i.e., detection of the gravitational waves Abbott2016v1 ; Abbott2016v2 and observation of the shadow of the M87 supermassive black hole (BH) Akiyama2019 ) and Standard Model (SM) (i.e., detection of the Higgs boson HiggsBoson ) of particle physics are incomplete theories. While Einstein’s theory of GR successfully describes gravity at a classical level, SM explains particles and the other three fundamental forces (electromagnetic, and the strong and weak nuclear forces) at a quantum level. The unification of GR and SM is a fundamental quest, and this success will necessarily lead us to a deeper understanding of nature. In the search for this unification, some quantum gravity theories (QGTs) have been proposed, but direct tests of their features are beyond the energy scale of the currently available experiments. Because, they will be observed on the Planck scale which is around (). However, it is possible that some signals of the QGT appear at sufficiently low energy scales and their effects can be observed in experiments on existing energy scales. One of these signals could be related to the LSB Nomura2013 .
The theory of LSB has been under intense research since the proposed SM Extension (SME) SME1v1 ; SME1v2 ; SME2 ; SME3v1 ; SME3v2 ; SME3v3 ; SME3v4 ; SME4v1 ; SME4v2 ; SME4v3 , which is an effective field theory that includes the SM, GR, and every possible operator that breaks the Lorentz symmetry. With the SME, further investigations of the LSB can be made in the context of high energy particle physics, nuclear physics, gravitational physics, and astrophysics. The simplest models that contain a vector field which dynamically breaks the Lorentz symmetry are called bumblebee models mybumblev1 ; mybumblev2 ; mybumblev3 ; mybumblev4 ; mybumblev5 .These models, although owning a simpler form, have interesting features such as rotations, boosts, and CPT violations. In a bumblebee gravity model (BGM), potential is included in the action , which evokes a vacuum expectation value (VEV) for the vector field. The potential is formed as a function of a scalar combination of the vector and the metric (plus the other matter fields, if there are any). The potential has a minimum at . At the , the bumblebee field incorporates a vacuum value shown by , which is the so-called vacuum vector. In fact, the vacuum vector is nothing but a background vector that gives rise to local (spontaneous) LSB myLSB . The scalar of the BGM, in general, reads as in which is a constant having dimensions of mass (). Thus, the satisfies the condition of for . Here, is spontaneously induced as a timelike vector abiding by . For instance, the aether models aether1 ; aether1v ; aether2 ; aether3 are based on a vector field, which is in the Lagrangian density of the system with a non-vanishing VEV. The vector field dynamically selects a preferred frame at each point in the considered spacetime and spontaneously breaks the Lorentz invariance. This is a mechanism reminiscent of the breaking of local gauge symmetry described by the Higgs mechanism. In general, the subclass of aether models obeys the following action aetherAct :
[TABLE]
where the parameter having dimensions of denotes the coupling between the Ricci tensor () and . is the bumblebee field strength:
[TABLE]
As mentioned above, is the potential of the bumblebee field that drives the breaking of the Lorentz symmetry of the Lagrangian by collapsing onto a non-zero minimum at or In fact, is one of the Lorentz breaking coefficients and it shows a preferred direction in which the equivalence-principle is locally broken for a certain Lorentz frame. Observations of Lorentz violation can emerge if the particles or fields interact with the bumblebee field aetherAct . It is worth noting that when a smooth quadratic potential is chosen as
[TABLE]
where is a dimensionless constant, one gets the Nambu-Goldstone excitations (massless bosons) besides the massive excitations aetherAct . Besides, the linear Lagrange-multiplier potential is given by These potentials (1) and (2) present also the breaking of the gauge invariance and other implications to the behavior of the matter sector, the photon, and the graviton. For a topical review (from experimental proposals to the test results) of the BGMs, the reader is referred to mybumblev4 and references therein. Furthermore, the studies using the bumblebee models have gained momentum for the last two decades. The vacuum solutions for the bumblebee field for purely radial, temporal-radial, and temporal-axial Lorentz symmetry breaking were obtained in myex1 . New spherically static black hole (BH) myex2 and traversable wormhole myex3 solutions in the BGM have been recently discovered. Bluhm myex4 discussed the Higgs mechanism in the BGM. The electrodynamics of the bumblebee fields was studied by myex5 in which the bumblebee field was considered as a photon field. Propagation velocity of the photon field, along with its possible effects on the accelerator physics and cosmic ray observations, was also investigated. BGMs are also used to limit the likelihood of Lorentz violation in astrophysical objects such as the Sun myex6 . For other studies demonstrating the physical effects (quasinormal modes, thermodynamics, etc.) of the bumblebee field, the reader may refer to mybumble2v1 ; mybumble2v2 ; mybumble2v3 ; mybumble2v4 ; mybumble2v5 ; mybumble2v6 ; mybumble2v7 ; mybumble2v8 ; mybumble2v9 ; mybumble2v10 and references therein.
Hawking’s ground-breaking studies myHR1 ; myHR2 can be considered as the onset of QGT myWald ; WaldGBF . Since then there have been numerous research papers on the subject of Hawking radiation (HR) in the literature (see, for instance, 5v1 ; 5v2 ; 5v3 ; 5v4 ; 5v5 ; 5v6 ; 5v7 ; 5v8 ; 5v9 ; 5v10 ; 5v11 ; 5v12 ; 5v13 ; 5v14 ; 5v15 ; 5v16 ; 5v17 ; 5v18 ; 5v19 ; 5v20 ; 5v21 ; 6v1 ; 6v2 ; 6v3 ; 6v4 ; 6v5 ; 6v6 ; 6v7 ; 6v8 ; 6v9 ; 6v10 ; 6v11 ; 6v12 ; 6v13 ). Several methods have been developed to calculate the HR of BHs 7 ; 8 ; 9 ; 10 . In this study, we mainly focus on the quantum gravity effects on the HR of SBHBGM myex2 in the tunneling paradigm. Although a number of QGTs have been proposed, however, physics literature does not as yet have a complete and consistent QGT. In the absence of a complete quantum description of the HR, we use effective models to describe the quantum gravitational behavior of the BH evaporation. In particular, string theory, loop quantum gravity, and quantum geometry predict the minimal observable length on the Planck scale 11 ; 12 , which leads to the GUP GUP13v1 ; GUP13v2 ; GUP13v3 ; GUP14v1 ; GUP14v2 ; GUP14v3 ; GUP14v4 ; GUP14v5 ; GUP14v6 ; GUP14v7 ; GUP14v8 ; GUP14v9 ; GUP14v10 ; GUP14v11 ; GUP14v12 ; GUP14v13 ; GUP14v14 ; GUP14v15 ; GUP14v16 ; GUP14v17 ; GUP14v18 ; GUP14v19 ; GUP14v20 ; GUP14v21 ; GUP14v22 ; GUP15v1 ; GUP15v2 ; GUP15v3 ; GUP15v4 ; GUP15v5 ; GUP15v6 ; GUP15v7 ; GUP15v8 ; GUP15v9 ; GUP15v10 ; GUP15v11 ; GUP15v12 :
[TABLE]
where in which denotes the Planck mass and is the dimensionless parameter, which encodes the quantum gravity effects on the particle dynamics. The upper bound for was obtained as 16 . Today, the effects of GUP on BHs have been extensively studied in the literature 17 ; 18 ; 19 . To amalgamate the GUP with the considered wave equation, the Wentzel-Kramers-Brillouin approximation 20 is generally used. Thus, one can obtain the quantum corrections to the HR of the BH 21 ; 22 .
Since the GUP and LSB effects are high energy modifications of the QGT, it is interesting to investigate their combined effects. To this end, we study the GUP-assisted HR of bosons (spin-[math]) and fermions’ (spin-) tunneling 24 ; 25 from the SBHBGM. Although the SBHBGM looks like the Schwarzschild BH, the differences in the Kretschmann scalars confirm that both BHs are physically different. The effects of spin and Lorentz-violating parameter 26 ; 27 on the quantum corrected HR are analyzed. We also study the problem of low energy greybody factors 28 ; 29 for the bosons and fermions emitted by the SBHBGM. For this purpose, we implement a method developed by Unruh 30 ; 31 . It is also worth noting that Lorentz invariant massive gravity can be obtained dynamically from spontaneous symmetry breaking in a topological Poincare gauge theory massive1 . Besides, BH radiation in massive gravity (selecting a preferred direction of time) naturally corresponds to violations of the Lorentz symmetry massive2v1 ; massive2v2 ; massive2v3 .
The outline of the paper is as follows: In Sec. 2, we briefly introduce the SBHBGM and discuss some of its basic features. Section 3 is devoted to the computation of GUP-corrected HR of the bosons’ tunneling from the SBHBGM. In Sec. 4, we compute the quantum tunneling rate for the fermions of the SBHBGM using the GUP-modified Dirac equation and derive the modified HR. In the following section, we derive the greybody factor of the SBHBGM. In Sec. 6, we summarize our results. (Throughout the paper, we use geometrized units: .)
II SBHBGM SPACETIME
The Lagrangian density of the BGM 32 ; 33 yields the following extended vacuum Einstein equations
[TABLE]
where and are the Einstein and bumblebee energy-momentum tensors, respectively. is the gravitational coupling and is given by
[TABLE]
where is the real coupling constant (having dimension ) that controls the non-minimal gravity-bumblebee interaction. From now on, the prime symbol shall denote the differentiation with respect to its argument. Meanwhile, there are other generic bumblebee models having non-zero torsion in the literature (see for instance 32 ). In Eq. (6), the potential provides a non-vanishing VEV for . As it was stated above (see also 36 ; 37 ), the VEV of the bumblebee field is determined when Taking the covariant divergence of the bumblebee Einstein equations (5) and using the contracted Bianchi identities, one gets
[TABLE]
which gives the covariant conservation law for the bumblebee total energy-momentum tensor . Thus, Eq. (5) reduces to
[TABLE]
One can immediately see that when the bumblebee field vanishes, we recover the ordinary Einstein equations. Recently, the vacuum solution in the BGM induced by the LSB has been derived by Casana et al. myex2 . The solution is obtained when the bumblebee field remains frozen in its VEV 38 ; 38n . Namely, we have
[TABLE]
Thus, the extended Einstein equations are found to be
[TABLE]
Assuming a spacelike background for as
[TABLE]
and using the condition constant, LSB parameter () is defined as 32 . A spherically symmetric static vacuum solution to Eq. (10) is obtained as follows myex2
[TABLE]
which we call it SBHBGM solution. This BH solution represents a purely radial Lorentz-violation outside a spherical body characterizing a modified BH solution. In the limit one can immediately see that the usual Schwarzschild metric is recovered. For the metric (12), the Kretschmann scalar becomes
[TABLE]
which is different than the Kretschmann scalar of the Schwarzschild BH. It means that none of the coordinate transformations link the metric (12) to the usual Schwarzschild BH. When , Eq. (12) becomes finite: the coordinate singularity can be removed by applying a proper coordinate transformation. However, in the case of , physical singularity cannot be removed. So, we see that the behaviors of the physical () and coordinate ( : event horizon) singularities do not change in the BGM.
The Hawking temperature of the metric (12) can be computed from Eq. (1), in which the surface gravity is given by myWald
[TABLE]
where is the timelike Killing vector field. Thus, the Hawking temperature of the SBHBGM (12) reads
[TABLE]
One can easily see from Eq. (15) that the non-zero LSB parameter has the effect of reducing the Hawking temperature of a Schwarzschild BH.
III GUP ASSISTED HR OF SBHBGM: BOSONS’ TUNNELING
The generic Klein-Gordon equation within the framework of GUP is given by 36
[TABLE]
where and are the GUP parameter and mass of the scalar particle, respectively. Introducing the following ansatz for the wave function
[TABLE]
where is the classically forbidden action for quantum tunneling. Substituting Eq. (17), together with the metric functions of line-element (12), into Eq. (16), we get
[TABLE]
[TABLE]
where
[TABLE]
It is easy to see that SBHBGM (12) admits two Killing vectors . The existence of these symmetries implies that we can assume a following separable solution for the action
[TABLE]
where and denote the energy and angular momentum of the radiated particle, respectively. Substituting Eq. (20) in Eq. (18), we obtain
[TABLE]
[TABLE]
We focus only on the radial trajectories in which only the sector is considered. Thus, one can set
[TABLE]
where is a constant. So, Eq. (21) becomes
[TABLE]
which can be rewritten as a bi-quadratic equation as follows
[TABLE]
where
[TABLE]
[TABLE]
[TABLE]
Eq. (24) has four roots if . We deduced from our analytical computations that only two roots () have physical meaning at the event horizon of the SBHBGM. These roots are
[TABLE]
It is worth noting that a sign represents an outgoing/ingoing wave. On the other hand, the integrand of the integral (28) has a pole at . Evaluating the integral by using the Cauchy’s integral formula, we obtain the imaginary part of the action as
[TABLE]
Thus, the tunneling rate of the scalar particles becomes
[TABLE]
Recalling the expression of the Boltzmann factor
[TABLE]
one can read the modified Hawking temperature () as follows
[TABLE]
As can be seen above, after terminating the GUP parameter i.e., , one can recover the standard Hawking temperature (15).
IV GUP-ASSISTED HR OF SBHBGM: FERMIONS’ TUNNELING
In this section, we aim to derive the modified Hawking temperature in the case of radiating fermions. To this end, we consider the Dirac equation, which is given by 37
[TABLE]
where denotes the test spinor field. The matrices for the metric (12) are given by
[TABLE]
in which ’s represent the well-known Pauli matrices 39 . One can easily ignore the terms having since is the effect of quantum gravity and it is a relatively very small quantity. For spin-up particles, the wave function can be expressed as 37
[TABLE]
where , and are functions of coordinates is the action of the emitted fermion. It is worth noting that here we only consider the spin-up case since it is physically same with the spin-down case; the only difference is the sign. Substitution of the wave function in the generalized Dirac equation (33) results in the following coupled equations
[TABLE]
and
[TABLE]
Then, one can get the following decoupled equations
[TABLE]
and
[TABLE]
By using the fact that SBHBGM spacetime has a timelike Killing vector , one can obtain the radial action by performing the separation of variables technique:
[TABLE]
where is the fermion energy. Substituting Eq. (40) into Eqs. (38) and (39), we find out the identical equations for and equations. Thus, we have
[TABLE]
or
[TABLE]
It is obvious that the expression inside the square brackets can not vanish; thus, one should have
[TABLE]
and the solution of , therefore, does not contribute to the tunneling rate. The above result helps us to simplify Eqs. (36) and (37) [with ansatz (40)] as follows
[TABLE]
and
[TABLE]
In the simple way, one can set
[TABLE]
[TABLE]
where
[TABLE]
and
[TABLE]
After making some manipulations, we see that
[TABLE]
which yields
[TABLE]
where
[TABLE]
[TABLE]
[TABLE]
[TABLE]
ignoring terms, Eq. (51) reduces to
[TABLE]
Therefore, we have
[TABLE]
in another form
[TABLE]
Recalling Eq. (30), we find the tunneling rate of fermions as follows
[TABLE]
Thus, with the help of the Boltzmann factor (31), we get the GUP-consolidated temperature of the SBHBGM via the emission of the fermions:
[TABLE]
in which represents the original Hawking temperature (15)
[TABLE]
The above result (60) shows that GUP corrected temperature deviates from the standard Hawking temperature.
V Greybody Factors of SBHBGM
In this section, we shall first derive the effective potentials of the scalar and fermion perturbations in the geometry of the SBHBGM. Then, the obtained effective potentials will be used for computing the greybody factors of the SBHBGM. The results will be depicted with some plots and discussed.
V.1 Scalar Perturbations of SBHBGM
The massless Klein- Gordon equation is given by
[TABLE]
where the D’Alembert operator is denoted by the box symbol and . For the SBHBGM (12), we have
[TABLE]
and therefore Eq. (62) reads
[TABLE]
We invoke the following ansatz for the scalar field in the above equation:
[TABLE]
so that we have
[TABLE]
When one changes the independent variable to , the angular equation is found to be
[TABLE]
where denotes the eigenvalue. The above equation is nothing but the Legendre differential equation when one sets
[TABLE]
The radial equation then becomes
[TABLE]
Introducing a new variable we get a Schrödinger-like wave equation
[TABLE]
where is the tortoise coordinate defined by
[TABLE]
The effective potential felt by the scalar field then becomes
[TABLE]
It is obvious from Fig. (1) that the effective potential vanishes both at the event horizon of the SBHBG and at spatial infinity. This behavior will help us to analytically derive the greybody factor of the scalar field emission from the SBHBG.
V.2 Fermion Perturbations of SBHBGM
In this subsection, we shall employ the Newman-Penrose formalism 40 to find the effective potential of the fermion fields propagating in the geometry of the SBHBGM. Chandrasekar-Dirac equations (CDEs) are given by Chandram
[TABLE]
[TABLE]
[TABLE]
[TABLE]
where and represent the components of the wave functions or the so-called Dirac spinors. , and are the spin coefficients, and a bar over a quantity denotes complex conjugation. The non-zero spin coefficients are found to be
[TABLE]
To have separable solutions for the CDEs (73), we introduce the following ansatzes
[TABLE]
[TABLE]
[TABLE]
[TABLE]
where denotes the azimuthal number and is the frequency of the spinor fields. Since the directional derivatives Chandram are defined by and we have
[TABLE]
[TABLE]
[TABLE]
[TABLE]
After substituting Eqs. (74-76) into the CDEs (73), one can obtain the following set of equations:
[TABLE]
[TABLE]
[TABLE]
[TABLE]
where and are the angular operators, which are known as the laddering operators:
[TABLE]
which lead to the spin-weighted spheroidal harmonics SWSH1 ; SWSH2 with the following eigenvalue myTeukolsky ; LeePRD :
[TABLE]
By considering and then we reduce the CDEs (77) to two coupled differential equations:
[TABLE]
[TABLE]
Moreover, if one sets
[TABLE]
and substitute them into Eqs. (80) and (81), after some manipulations, we get:
[TABLE]
[TABLE]
By defining and and introducing the tortoise coordinate () as , we obtain
[TABLE]
[TABLE]
One can combine the above equations by letting
[TABLE]
[TABLE]
Thus, we end up with the following pair of one dimensional Schrödinger-like wave equations:
[TABLE]
where the effective potentials for the Dirac field read
[TABLE]
V.3 Greybody Factor Computations
In general relativity, the greybody factor is one of the most important physical quantities related to the quantum nature of a BH. A high value of the greybody factor indicates a high probability that HR can reach to spatial infinity. Among the many methods (see for example WaldGBF and references therein) for obtaining the greybody factor, here we employ the method of 41 , which formulates the general semi-analytic bounds for greybody factors:
[TABLE]
where are the dimensionless greybody factors that depend on the angular momentum quantum number and frequency of the emitted particles, and
[TABLE]
in which prime denotes the derivation with respect to . We have two conditions for the certain positive function : 1) and 2) 41 . Without loss of generality, we simply set which reduces the integration of Eq. (92) to
[TABLE]
For a massless scalar field , considering the effective potential given in Eq. (72), Eq. (92) becomes
[TABLE]
Taking cognizance of the integral part of Eq. (95):
[TABLE]
the greybody factor of the SBHBGM due to scalar field radiation yields
[TABLE]
When one considers the effective potential (90) of the massless Dirac fields:
[TABLE]
the integral seen in Eq. (92) can be easily computed. Thus, we find the greybody factor expression of the SBHBGM arising from the fermion radiation:
[TABLE]
From now on, without loss of generality, we consider only . After some manipulation, one can get
[TABLE]
which recasts in
[TABLE]
or, in more compact form:
[TABLE]
We depict the greybody factors of the SBHBGM arising from the scalar (97) and fermion (102) fields in Figs. (3) and (4), respectively. As is well-known, the greybody factor of the HR must be since a BH does not perform a complete black body radiation with a 100% absorption coefficient. Our findings, as shown in Figs. (3) and (4), are in good agreement with the latter remark. Also, it can be seen from these figures that the peak values of the greybody factors decreases with increasing LSB parameter . In summary, LSB has a greybody factor-reducing effect.
VI CONCLUSION
In this paper, we studied the quantum thermodynamics QTerm of the SBHBGM. During this analysis, we had mainly two aims: 1) to obtain the modified Hawking temperature of the SBHBGM, within the framework of GUP, arising from the emission of bosons and fermions; 2) to compute the greybody factors of the scalar and fermion fields from the SBHBGM. To this end, we first derived the effective potentials of the Klein-Gordon and Dirac equations. Next, we used the obtained effective potentials in the greybody expression (92). Then, we illustrated the obtained greybody factors in Figs. (3) and (4). It was clear from those figures that as LSB effect () increases, the greybody factor decreases: low values of the greybody factor indicate a low probability that HR can reach spatial infinity. In the future, the latter observation might shed light on the LSB effects from both terrestrial experiments and astrophysical observations. Any discovery of the LSB would be an important signal beyond the SM physics Review .
In future work, we plan to extend the GUP and greybody factor analysis Visser to the various BHs in gravity’s rainbow, which is also a result of quantum gravity rainbowv1 ; rainbowv2 ; rainbowv3 ; rainbowv4 ; rainbowv5 ; rainbowv6 ; rainbowv7 ; rainbowv8 ; rainbowv9 ; rainbowv10 . The deformation of a spacetime owing to the rainbow gravity effect leads to Lorentz violations LVRGv1 ; LVRGv2 . In this way, we hope to achieve new results that will help us to understand the QGT and its effect on the LSB.
Acknowledgements
The authors are grateful to the editor and anonymous referees for their valuable comments and suggestions to improve the paper.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1(1) B. P. Abbott et al. [LIGO Scientific and Virgo Collaborations], Phys. Rev. Lett. 116 , 061102 (2016).
- 2(2) B. P. Abbott et al. [LIGO Scientific and Virgo Collaborations], Phys. Rev. Lett. 116 , 241103 (2016).
- 3(3) K. Akiyama et al. [Event Horizon Telescope Collaboration], Astrophys. J. 875 , L 1 (2019).
- 4(4) G. Aad et al. [ATLAS Collaboration], Science 338 , 1576 (2012).
- 5(5) Y. Nomura, J. Varela and S. J. Weinberg, Phys. Lett. B 733 , 126 (2014).
- 6(6) D. Colladay and V. A. Kostelecky, Phys. Rev. D 55, 6760 (1997).
- 7(7) D. Colladay and V. A. Kostelecky, Phys. Rev. D 58, 116002 (1998).
- 8(8) S.R. Coleman, and S.L. Glashow, Phys. Rev. D 59, 116008 (1999).
