The Standard Model Coupled to Quantum Gravitodynamics

We show that the renormalizable SO(4) X U (1) X SU (2) X SU (3) Yang Mills coupled to matter and the Higgs field fits all the experimentally observed differential cross sections known in nature. This extended Standard Model reproduces the experimental gravitational differential cross sections without resorting to the graviton field and instead by exchanging SO(4) gauge fields. By construction, each SO(4) generator in quantum gravitodynamics does not commute with the Dirac gamma matrices. This produces additional interactions absent to non-Abelian gauge fields in the Standard Model. The contributions from these new terms yield differential cross sections consistent with the Newtonian and post Newtonian interactions derived from General Relativity. Dimensional analysis of the Lagrangian shows that all its terms have total dimensionality four or less and therefore that all physical quantities in the theory renormalize by finite amounts. These properties make QGD the only renormalizable 4-dimensional theory describing gravitational interactions.


Introduction
The structure and success of General Relativity have lead us to the believe that the spin-2 graviton field mediates quantum gravitational interactions; but to date this believe has eluded experimental detection. Furthermore, the quantum limit of General Relativity fails to renormalize properly forcing us to further believe that General Relativity comes from a classical projection of a 10 dimensional string. On the other hand, we know that the geometrical object that describes the geodesics in General Relativity is the connection and not the metric which enters that construction only after arbitrarily imposing the metric compatibility constraint. Then, we can in practice construct a quantum theory of gravity where the fundamental field is the connection as in [1] or a similar object as we do here. Nothing but experimental observation constrains these possibilities. Therefore, we should not force the graviton upon ourselves and instead focus on matching the proposed model's differential cross sections to experimental evidence; more so, when we must yet detect the particle mediating gravitational interactions.
Here we present quantum gravitodynamics, QGD, the only renormalizable 4-dimensional theory that reproduces the interactions in General Relativity to O( v 2 c 2 ) and couples to the Standard Model through a covariant derivative like the other three interactions in nature. As we show below, all diagrams in QGD match those of it counterpart in General Relativity to O( v 2 c 2 ) and the potential that determines a test particle's classical trajectories, are identical for both theories to that same order. As opposed to General Relativity which does not renormalize, QGD has a metric void of any dynamics, explicitly lacks a graviton field or any spin-2 state, and instead uses SO(4) Yang Mills fields coupled to spinor singlets to mediate the gravitational interactions. Furthermore, QGD reproduces all diagrams in General Relativity to O( v 2 c 2 ) and therefore all experimental evidence to that same order. However, as opposed to General Relativity, the terms in the QGD Lagrangian involve products of fields and their derivatives with total dimensionality of at most four and therefore renormalization follows [2]. Thus, QGD is renormalizable, and, like General Relativity, constructed in 4 dimensions. Because QGD is defined in terms of Yang Mills fields, it also couples naturally to the Standard Model to produce a Unified Model with SO(4)×U(1)×SU(2)×SU(3) gauge symmetry that reproduces all the experimentally observed differential cross sections in nature, including those produced by gravitational interactions to O( v 2 c 2 ). QGD also has the ability to fit those diagrams contributions from gravity that may become available in the future such as Higgs-Gluon gravitational scattering. In addition, QGD has gravitational SU(2) Yang Mills BPS states that match the entropy of black holes and type-II strings. Pure Yang's gravity [1] has similar classical equations to those of Einstein's Relativity and includes Schwarzschild-like solutions, whose singularity at the horizon, as opposed to black hole solutions, cannot be removed through a coordinate transformation [3]. QGD and Yang's gravity are both constructed in terms of the same connection defined over a fiber bundle, However, Yang's gravity is not renormalizable because it couples its Yang-Mills field directly to the energy-momentum tensor. Instead, in QGD, the Yang Mills fields couple to matter through a SO(4) × U(1) covariant derivative, in the same manner as the other forces in nature and which preserves renormalization. Therefore, the couplings in QGD ensure that all terms in the Lagrangian have dimensionality four or less which guarantees renormalization [2].
String theory is the only other theory, besides QGD, to incorporate gravitational interactions while preserving renormalization. It has been a successful gravitational theory, because like QGD, scattering processes associated with the β-function constrain the metric to be a solution of Einstein's equation. Therefore, like QGD it reproduces all General Relativity tests to O( v 2 c 2 ). Like strings, QGD has the necessary BPS states to resolve the black hole entropy problem. However, unlike QGD, string theory cannot be defined without ambiguity in 4 dimensions and neither can it couple to the Standard Model which instead must be reformulated along side gravity.
The QGD connection defined on the SO(4) fiber bundle can be projected to the connection defined on the SO(1,3) tangent bundle. This projection allows the generators of the SO(4) gauge fields to be defined in terms of Dirac gamma matrices and therefore these generators do not commute with the spacetime Dirac γ µ matrices. This non-commutativity contrasts QGD from the other gauge theories in the Standard Model because it introduces terms proportional to [γ µ , T ab ] which couple to fermionic matter fields with arbitrary couplings. The map from the SO(4) gravitational gauge field, ω ab µ , to the spin-connection defined in [4] motivate calling ω ab µ the connecton. In order to test QGD, we must reproduce the experimental evidence supporting General Relativity. As discussed at the end of section 4, this also suffices to show that QGD will reproduce all effects of General Relativity to O( v 2 c 2 ). The procedure used to fit QED to experimental evidence, see for example [5], transfers well to QGD and permits the fitting of all the QGD parameters to experimental observation. The T-matrix calculation in [5] associated with the tree-level e − e − → e − e − scattering process produces a differential cross section; comparison of that differential cross section with the one obtained from the Schroedinger equation with a Coulomb potential allows the fitting of the electron coupling found in QED to the experimentally measured charge used in the Coulomb potential. This procedure extends to the quantized but non-renormalizable General Relativity theory. In [6] the tree-level scattering process for two Bosons, b b ′ → b b ′ , produces a differential scattering matrix which fits the Schroedinger equation with a Newton potential while the 1-loop scattering correction fits the differential cross section produced by the Schroedinger equation with the leading classical post Newtonian correction to the potential. The strategy implemented in [5] and [6] was also successfully implemented to the gravitational sector in string theory [7]- [10]. Thus, matching the model's differential cross sections to experimental/expected differential cross sections allows us to test if a model correctly describes nature.
Here we follow the same strategy as that used in [5]- [10] to show that QGD coupled to the Standard Model fits all the known experimental differential cross sections, including those gravitational in nature. Section 2 details the symmetries and action derived for a single matter field and the U(1) electromagnetic field, both coupled to gravity. In section 3 we calculate the propagators and vertices of the model presented in section 2. In section 4 we calculate T-matrix elements to produce differential cross sections which are then fitted to the differential cross section of the Schroedinger equations with the Coulomb, the Newtonian and post Newtonian potentials as well as the differential cross section for the deflection of light by a point particle in the small angle approximation found in General Relativity. Section 5 presents the full QGD coupling to the Standard Model while section 6 fits the differential cross sections to experimentally observed differential cross sections. Section 7 moves outside the realm of experimental observation and shows that QGD can include the expected gravitational interactions of the SU(2) and SU(3) gauge Bosons as well as the Higgs doublet. In section 8 we consider the isomorphic SU(2) × SU(2) Yang Mills representation of QGD, which when N = 2 supersymmetrized has BPS states whose entropy corresponds to that of a black hole, and that can bridge to type-II strings [13]. The conclusions appear in section 9.

Symmetries and Actions
The fundamental geometrical object of General Relativity, the connection on the tangent bundle, does not require a metric to describe a particle's geodesics. Consistent with General Relativity and in the spirit of [1], QGD chooses a connection over the metric as the fundamental field. QGD does not use a yet to be detected graviton field or spin-2 particle and it uses instead a constant metric void of any dynamics. Instead, connections defined on the SO(4) fiber bundle mediate all gravitational interactions; these SO(4) gauge fields alone mediates all the gravitational interactions between particles. These fields suffice to obtain all the differential cross sections expected from General Relativity for a single matter-photon system. More importantly, this construction does not have any experimental impediments because we have yet to detect the graviton or any other spin-2 particle.
We present the Lagrangian We use the following conventions. The diagonal and flat space-time metric, η µν , does not carry any dynamics and has signature (+, −, −, −). ǫ µνρσ represents the Levi-Civita tensor in 4 dimensions. The vielbeins e a ρ connect space-time coordinates (Greek indices) with the SO(4) fiber bundle coordinates (Latin indices). In particular e a ρ e bρ = η ab and e aρ e σ a = η ρσ . The metric η ab , defined on the fiber bundle, has signature (+, +, +, +).
We define the field strengths appearing in (3) A µ represents the U(1) electromagnetic gauge field and the SO(4) gauge field ω ab µ , the connecton, mediates of gravitational interactions. The couplings e, g, g ′ and α denote coupling constants and β = √ 1 + 6α 2 . Withω µab defined as a copy of a textbook SO(N=4) gauge field necessarily antisymmetric in indices a and b, then ω µab =ω µab + e c µ e ν aω νcb , so that the indispensable property ω µab e a ν e b ρ ǫ µνρσ = 0 (8) holds.
The covariant derivative operator D µ defined in (4) acts on a 4 dimensional Dirac spinor field ψ and depends on the SO(4) generators T ab defined in (1). The U(1) local gauge transformation Γ = exp(iλ), with λ a scalar function, acts on the following fields which also transforms the covariant derivative in the standard way but leaves (2) unchanged. The SO(4) local gauge transformations withΓ = exp(iΛ ab T ab ) and Λ ab a 2-Tensor, act on the following fields which also transforms the covariant derivative in the standard way but also leaves (2) However, the SO(1,3) gauge transformations, can not act on the pure Yang-Mills Lagrangian because the Cartan-Killing metric for that gauge group, given by the trace over two group elements, is not positive definite which yields states not bounded from below [2]. This property therefore distinguishes theories defined over the tangent bundle with local SO(1,3) symmetry from those defined over the fiber bundle with local SO(4) symmetry.
The former cannot be used to define Yang-Mills theories even when the fiber bundles can be projected to the tangent bundle. For the avoidance of doubt, gravitational theories like spin-gravity can be defined over the tangent bundle only because, as a result of the different construction from Yang-Mills, their spectrum does not become unbounded from below when the Cartan-Killing metric is not positive definite. The Lorenz gauge fixings and property (7) simplify the quadratic component in (3) to The structure constant defined in [4] C abcdef = −η af η ce η bd + η ad η ce η bf − η ae η cb η df + η ac η eb η df simplifies the relevant terms of the cubic terms in (3) to The terms in (23) take the form ωω∂ω and ∂AAω. In the low energy limit the photon momentum, r, and the connecton momentum, k, satisfy or equivalently This condition accounts for the omission of all other terms. For example, for matter-photon scattering the vertex generated by Aω∂ω would only participate at 1-loop level, would necessarily appear twice and thus takes values proportional to the square of the connecton momentum. Instead, the vertices produced by the second term in (23) would appear first at tree-level and then at 1-loop level where it would contribute terms proportional to the photon momentum and its square respectively. Dimensional analysis of General Relativity shows that it does not renormalize properly. However, it follows from (9) that the connecton has the same dimensionality as the U(1) gauge field. Therefore, the terms in the Lagrangian (21) and (23) involve products of fields and their derivatives with total dimensionality four or less. Renormalization then follows as shown in [2].
Note that (28) does not depend on ∂ µ . This means that gauge invariance constrains (27) but not (28). Thus the rescaling of (27) by a constant redefines the coupling. However, the redefinition guarantees gauge invariance only after a subsequent redefinition of the spinor ψ. On the other hand, a recalling by a constant, even a complex one, of (28) does not affect gauge invariance. This contrasts with SU(2) of SU(3) couplings where the Dirac matrices commute with their generators and therefore (28) necessarily vanishes. The relations (27) and (28) simplify the Lagrangian (4) to Here the real coupling constants e, g, α and g ′ adjust the theory to observation. In particular e represents the electric charge and gauge invariance constrains g and coincides with the coupling constant found in the purely bosonic sector (23). Instead, gauge invariance does not constrain g ′ which can take any value, real, purely imaginary or complex.
Again, the terms in the Lagrangian (29) involve products of fields and their derivatives with total dimensionality four or less. Therefore, the terms in the Lagrangians (21), (23) and (29), and thus Lagrangian (2), involve products of fields and their derivatives with total dimensionality four or less. Renormalization then follows as shown in [2].
The Lagrangian L f has an complex coupling constant g ′ while L gauge has a real coupling constant g. Then, self energy corrections can have imaginary contributions and therefore it would be tempting to assume that the probability density determined by the wavefunction decreases over time and unitarity is lost. This is not the case.
We choose a pure real g and a pure imaginary g ′ and will show when fitting theory to experiment in the subsequent sections that where G represents the gravitational constant. The self-energy contributions to O(G 2 ) for both Bosons A µ and ω ab µ come from purely bosonic diagrams and those with fermionic loops. We construct the purely bosonic diagrams with vertex (c) and vertex (d) defined in section 3. Therefore, the diagrams depend only on the real valued coupling constant g Σ (G) Thus, those diagrams always contribute real amounts to the self energies Σ . As we show below, real contributions to the self energy do not affect unitarity.
The diagrams involving fermionic loops always involve an even number of vertices (a) or an even number of vertices (b) defined in section 3. Therefore, corrections to O(G 2 ) to the self energies of both A µ and ω ab µ always contribute even powers of g ′ Σ (G) for some real constants a 1 and a 2 . Even though the coupling g ′ is imaginary, all corrections to the self energies for both Bosons A µ and ω ab µ contribute real amounts to O(G 2 ), Even more, because Fermions always appear in loops regardless of the order of perturbation, the self energy corrections to the Bosons A µ and ω ab µ always remain real. On the other hand, the contributions to the self energies of the Fermion contribute real amounts to O(G) but complex amounts to O(G 2 ). The contribution to the Fermion self energy to O(G) comes from diagrams constructed from a pair of vertices (a) or a pair of vertices (b). Therefore for some real constants b 1 and b 2 . The first imaginary contribution comes from the diagram involving one vertex (c) and three vertices (a). Therefore Therefore, the leading imaginary correction to the Fermion self energy contributes at most of O(G 2 ). The time evolution of the wavefunction for the Fermion where |ψ(t in ) > describes the wavefunction at the emitter and |ψ(t out ) > represents the wavefunction at the detector, implies that the probability density at time t out , < ψ(t out )|ψ(t out ) >, decays as When < ψ(t out )|ψ(t out ) ><< ψ(t in )|ψ(t in ) > particles are lost as they travel between emitter and detector. This experimentally unobserved reduction in flux prevents large imaginary values of Σ p . However, when Im(Σ f ) satisfies then, the self energy Σ f can take on imaginary values without loss of unitarity.
Since the leading imaginary correction to the self energy of the Fermion contributes at most O(G 2 ), then, whenever t out − t in < 10 34 years. Bearing in mind that the universe has existed less than 10 11 years and that estimates put the the proton's half life at 10 32 years, QGD can take imaginary couplings without affecting unitarity. This contrasts with other theories where the coupling constants are much larger than g ′ leading to a loss of unitarity if not of gauge invariance as well since g ′ = g. This conclusion is further supported by the absence of diagrams contributing to O(G 2 ) to the self energies constructed with vertex (a) and vertices involving the electromagnetic field.

Propagators and Vertices
The Lagrangian densities (21) and (29) yield the following propagators For the avoidance of doubt, in the limit considered here or equivalently m >> r >> k. The Lagrangian densities (23) and (29) generate the following vertices: In the following section we relate the tree and 1-loop level diagrams to physical processes in order to determine the four parameters in QGD.

Tree and 1-loop Diagrams and Their Relationship to Observed Quantities
In this section we calculate 4 scattering amplitudes: the tree-level photon exchange between two particles with equal mass, the tree-level connecton exchange and the 1-loop connecton exchange between two particles with equal mass and the tree-level connecton exchange between a photon and a particle with mass. These scattering amplitudes yield a total of 4 differential cross sections; we equate the first 3 to the differential cross section for a Schroedinger equation with 3 different potentials: the Coulomb potential, the Newton potential and the post-Newtonian potential while the last cross section matches the cross section for gravitational lensing of a photon by a massive point particle. The Lagrangian provides 4 coupling constants: e, g, g ′ and α. Therefore by the end of this section, the 4 coupling constants appearing in the differential cross sections calculated from QGD match those of the Schroedinger equation as well as that obtained from gravitational lensing in General Relativity.
The contribution from the photon exchange between two mass particles in diagram (a) produces the T-matrix suppresses the second term in (42). A further limiting to the spinless Schroedinger limit requires summation over the spin of all in-states and out-states so that < J Aµ >=< J ′ Aµ >= 0. Thus, in the non-relativistic quantum limit, the third term in (42) vanish exactly reducing (42) to This T-matrix yields the differential cross section while the Schroedinger equation with a Coulomb potential, with e exp , the measured electric charge, and R, the distance between the two particles, produces a differential cross section Comparing the differential cross section (45) and (47) yields Diagram (b) contributes to the differential scattering cross section of the connecton exchange between two mass particles. The T-matrix element produces the quantum gravitodynamic differential cross section The Schroedinger equation with a Newtonian potential, with G representing the gravitational constant, m exp representing the measured mass and R representing the distance between the two particles, produces a differential cross section Comparing the differential cross section (50) and (52) yields and The post-Newtonian correction comes from diagrams (c) and (d). All other 1-loop diagrams contribute to either analytic terms, or quantum corrections of order O(ln(−q 2 )) and do not contribute to the low energy limit considered here [6]. In addition, diagrams involving photons and connectons simultaneously do not contribute because the transformation ω ′ ab µ = ω ab µ − α β ǫ abν µ A ν leaves the contribution of the above diagram unchanged while all terms but the first two in (23) remain present in the action. These two diagrams contribute to the quantum gravitodynamic post Newtonian T-matrix This result as those below require the following steps. First, we focus only on the "electric" form factor which means that Equivalently, we disregard terms proportional to [γ µ , γ ν ] that contribute to the "magnetic" form factor. Second, the on-shell external momenta imply the relations Third, we suppress terms of O(q 4 ) using the low energy limit property q << p. Fourth, we use the approximationψγ µ ψ = 2p µ andψ ′ γ µ ψ ′ = 2p ′ µ . Finally, we use the appendix in [6] with the expressions for the several Feynman integrals.
The T-matrix (55) produces the differential scattering cross section The Schroedinger equation with a post Newtonian correction to the gravitational potential (see equation (39) of [6] ), where the constant a depends on the post-Newtonian expansion, produces, after using (53) and (54), the differential scattering cross section Comparing the differential cross section (59) and (61) yields Therefore, by fixing e, g and g ′ through equations (48), (53) and (62), QGD reproduces the low energy limit of General Relativity for the matter sector. Therefore, the SO(4) Yang Mills theory coupled to matter and an inert metric is a valid renormalizable relativistic quantum field theory to describe gravity's experimental observations.
Diagram (e) describes the scattering of a Fermion f with a photon γ that produces the T-matrix which in turn, in the small angle approximation where q = rsin(θ/2), yields the differential cross section where the photon energy, E γ , satisfies Derivation of (64) requires the photon on-shell conditions, along with the steps used for the matter-matter scattering above. Substitution of (53) and (62) into (64) yields In the low energy limit, equation (68) describes the Rutherford scattering of a mass m projectile off of a mass E γ target.
In the small angle deflection approximation and after setting the speed of light to 1, the gravitational lensing of a photon by a point particle of mass m in General Relativity relates the impact parameter b L with the deflection angle The relation between the impact parameter, b L , and the differential cross section, simplifies the latter to or equivalently Equation (72) describes the Rutherford scattering of mass E γ projectile off of a mass m target. The exchange m ↔ E γ modifies this differential cross section to which then describes the scattering of a mass m particle of a photon target with energy E γ . Equating (68) with (73) requires that 16 27 which determines α as a function of the known constant a. With (74) satisfied, Fourier transform in the Born approximation of (68) yields the low energy potential which describes the motion of a photon as it had mass equal to its energy in the Newtonian approximation.
Diagram (e) also describes the gravitational redshift. When the impact parameter vanishes, a photon starting at −∞ approaches the mass target and experiences the potential (75) which accelerates the photon and blueshifts its frequency. Thus the energy of the photon at point r 1 reads and when ω 1 >> ω 2 − ω 1 as required.
Diagram (e) has the same value as its counterpart diagram in General Relativity obtained following the methods of [6] which produces the exact same potential. However, renormalization aside, it does so in the frame of reference where the photon scatters off the Fermion, and again confirms that we can model photons as particles with mass E γ and (75). Like the gravitational redshift, the Shapiro delay also stems from diagram (e) and it is just the length of the path traveled by the photon when bending about a mass particle and after reflection from some other body retraces its path back to the source. This feature transpires the quantum limit were diagrams in General Relativity match the diagrams for QGD to O( v 2 c 2 ). In the classical limit, the trajectories of General Relativity on a Schwarzschild background also follow from consideration of a particle moving about a potential. Therefore, since we have shown that both theories have the same potential to O( v 2 c 2 ), any classical or quantum test that General Relativity passes will also be passed by QGD to that same order.
The QGD Lagrangian (2) has 4 free parameters e, g, g ′ and α. Equations (48), (53), (62) and (74), repeated here: 16 27 fix these parameters to experimentally observed quantities from QED, and General Relativity: the electric charge, the Newtonian mass and gravitational constant, the precession of a mass particle about another mass particle and the gravitational lensing of light by mass particles. Therefore, QGD suffices to describe the gravitational effects observed in nature without resorting to a spin-2 particle and using instead a quantum connection defined over an SO(4) fiber bundle. The dimensionality of the connecton is the same as that of the U(1) gauge field. Therefore, the terms in the Lagrangians (21), (23) and (29) involve products of fields and their derivatives with total dimensionality four or less. Then all scattering amplitudes considered in this section renormalize by finite amounts [2].
The matching of the above 4 differential cross sections with those of General Relativity coupled to matter and light are more than simple tests. By matching the differential cross sections we have, at least for two particle interactions, shown that the equations of motion for matter and light coupled to General Relativity exactly match those for QGD to O( v 2 c 2 ). This is manifested by the fact that the Schroedinger equation for both theories are the same. Therefore, any test that follows from General Relativity coupled to matter and light will be reproduced by QGD to order O( v 2 C 2 ). While it may very well be the case that General Relativity and QGD diverge beyond O( v 2 c 2 ), the reader is reminded that QGD needs to replicate nature and not General Relativity. We know that because QGD and General Relativity are indistinguishable to O( v 2 c 2 ), it replicates nature to that order. But General Relativity, safe for [12], has not been experimentally proven beyond that order, and therefore it may very well be that if QGD and General Relativity diverge beyond O( v 2 c 2 ), it will be up to nature to judge which most closely resembles it, always taking into account that General Relativity is at best and effective theory because it is not renormalizable like QGD.
It should be further noted that while the differential cross sections pertain two particle interactions and omit single particle and three-, four-particle interactions, these do not impose further constraints. The single particle diagrams involve loops which General Relativity cannot handle due to lack of renormalization. Furthermore, Three and four particle interactions are constrained by two particle scattering; see [6] for the diagrams used in General Relativity. Therefore, by considering only the equivalency of two particle scattering processes, we fix the equations of motions for all other scattering processes in both theories.

The Standard Model
Lack of evidence supporting the existence of a spin-2 particle along with the existence of a mapping between an SO(4) fiber bundle to an SO(1,3) tangent bundle motivate the formulation of a gravitational theory based on the SO(4) connecton instead of one based on the graviton. Without a graviton field and with a constant diagonal metric void of any dynamics, QGD again incorporates all the necessary gravitational interactions to the Standard Model using instead N × M SO(4) gauge fields ω suffices to obtain all the experimentally verified differential cross sections, and those expected from General Relativity.
The experimental evidence for General Relativity exists only for the matter and U(1) gauge sector. At present, no experimental evidence exists supporting that the SU(2), SU(3) and the Higgs sectors of the Standard Model couple to gravity; however, here we show that QGD can incorporate such couplings to satisfy the intuitive expectation of General Relativity.

The Lagrangian
µνab Ω (i,m)µνab (83) We define eachω µab as a copy of a textbook SO(4) gauge field necessarily having antisymmetry in indices a and b, then, so that ω (i,m) µab e a ν e b ρ ǫ µνρσ = 0. Since where φ µ = e a µ φ a and φ a , a = 1, ..., 4 describes the Higgs in the fundamental of SO(4), the constant couplings α im must take the form and 1 β 2 When representing a quark doublet, the spinors ψ (j) , acting on a quark doublet ψ (j) transforms covariantly under SO(4) × U(1) × SU(2) × SU(3). The U(1) transformation generated by Γ = e iλ acts on the following fields and leaves (82) SU (3) transformations generated by Γ = e iΛ a T a act on the following fields where C a bc represents the SU(3) structure constant, to leave (82) invariant while each D (i,m) µ ψ (j) transforms covariantly. Similarly, SU(2) transformations generated by Γ = e iΛ a τ a act on the following fields where C a bc represents the SU(2) structure constant, and also leave (82) invariant while each D (i,m) µ ψ (j) transforms covariantly. SO(4) transformations generated by Γ = e iΛ ab T ab act on the following fields and also leave (82) invariant while each D (109) The dimensionality of the connecton is the same as that of the U(1) gauge field. Therefore, the terms in the Lagrangian (109) involve products of fields and their derivatives with total dimensionality four.
The pure Standard Model SU(2) and SU(3) cubic and quartic bosonic interactions present in (83) do not mix with the SO(4) gauge fields and do not merit further treatment here. For the same reasons presented for (23), the relevant cubic terms in Lagrangian (83) include As in (23), the terms in (110) take the form ωω∂ω and ∂AAω. In the low energy limit the gauge Boson momentum, r, and the connecton momentum, k, satisfy r >> k. The terms in the Lagrangian (110) also involve products of fields and their derivatives with total dimensionality four. Therefore any diagrams constructed from the terms in (110) will renormalize by finite amounts [2].
After adding the contributions from (83) dependent on A 13 to (85) we obtain the Standard Model Higgs Lagrangian along with QGD interactions where φ d describes the Higgs in the SU(2) doublet representation instead of the SO(4) vector representation and 1 β 2 (112) The dual representation of the Higgs further strengthens the proposal of the SO(4) connecton as the carrier of quantum gravitational interactions. The terms in the Lagrangian (111) again involve products of fields and their derivatives with total dimensionality four. Thus, any corrections involving only these terms will renormalize by finite amounts [2].
When ψ (j) represents a quark doublet, the Lagrangian (84) simplifies to where e (u) and e (d) are the electric charges of each respective quark and Similarly, when ψ (j) represents a lepton doublet, the Lagrangian (84) simplifies to where e (ν) and e (e) are the electric charges of each respective lepton. In addition to the transformation of the terms in the pure Standard Model, the electroweak rotation which produces the U(1) gauge field A µ only transforms the last term in (113) and (116) where θ w describes the weak mixing angle and takes the same values as in the pure Standard Model. Thus, for a quark doublet with L f,sm representing the fermionic sector of the pure Standard Model after the rotation between B µ and W 3 µ . A similar expression follows for the lepton doublet.
The dimensionality of the connecton is the same as that of the U(1) gauge field. Therefore, the terms in the Lagrangian (118) involve products of fields and their derivatives with total dimensionality of at most four. Renormalization then follows [2]. More importantly, the whole Lagrangian, the sum of (109), (110), (111) and (118), involve products of fields and their derivatives also with total dimensionality at most four and therefore any corrections involving their terms renormalize by finite amounts [2].
The electroweak rotation transforms L (2) gauge into and it transforms L At this point we can apply the Higgs mechanism by selecting the usual VEV for the Higgs field which will produce the same effects that take place in the Standard Model. The rotation between the fields does not affect the dimensionality of the fields. Therefore, the rotated terms in the Lagrangian involve products of fields and their derivatives with total dimensionality at most four. Renormalization then follows [2].
with j = 1, ..., J running over the number of singlet spinors, g i = g i0 and α i = α i0 . The Lagrangian (119) reduces to and the Lagrangian (120) reduces to The dimensionality of the connecton is the same as that of the U(1) gauge field. Therefore, the terms in the Lagrangian (121), (122) and (123) involve products of fields and their derivatives with total dimensionality four. Renormalization then follows as shown in [2].
In sections 2 to 4 we considered a single particle, like the electron, and successfully reproduced the necessary post-Newtonian corrections to fit theory to experiment. Expanding the particle spectrum increases the number of scattering cross sections well beyond 4; This leads to a lack of solution when a single connecton describes the gravitational interactions. Instead we introduced N connectons to increase the number of couplings from g ′ ∈ I to g ji ∈ C and from g ∈ R to g i ∈ R. With this increase in the number of connectons the number of parameters in the model also increases sufficiently to fit all these differential cross sections.
The Lagrangian densities (121)-(123) yield the following propagators For the avoidance of doubt, in the low energy limit m >> ∂A >> ∂ω such that m >> r >> k. These same Lagrangian densities generate the following vertices: After assuming as before that the axial current vanishes, the scattering process in diagram (f) describing the photon exchange between particle j and k produces a T-matrix element, which in turn yields a differential cross section contribution, with µ jk defined as the reduce mass of particles j and k. On the other hand, the Schrodinger equation for a Coulomb potential with e j,exp and e k,exp describing the particles measured electric charge and R describing the distance between the two particles, yields the differential cross section dσ jk Comparing the differential cross section (126) and (128) constrains the parameters e j and e k to e j = e j,exp , e k = e k,exp .
Diagram (g), exchanging a connecton between particle j and k produces the T-matrix, which in turn produces the differential cross section, On the other side, the differential cross section obtained from the Schrodinger equation for a Newton potential, where m j,exp and m k,exp describe the particles measured mass and R the distance between them, produces the differential cross section, Comparing the differential cross sections (131) and (133) imposes the constraint Diagrams (h) and (i) produce the 1-loop correction to the connecton exchange between particles j and k to the T-matrix, after using the same steps used in section 4, and in turn it produces the differential cross section, (136) The Schroedinger equation with the Einstein correction to the gravitational potential between two particles, yields the differential scattering cross section Comparing the differential cross section (138) and (136) yields N i g i g ji g ki (g ji (m j +17m k )+g ki (17m j +m k )) = 16 3 π 3 6 aG 2 m 2 j,exp m 2 k,exp (m j,exp +m k,exp ). (139) The T-matrix contribution from diagram (j) scattering a mass particle off of a photon simplifies to After using the small angle approximation, the differential cross section obtained from this T-matrix equates to In the low energy limit, equation (141) describes the Rutherford scattering of a mass m j projectile off of a mass E γ target. For a small angle deflection in General Relativity the gravitational lensing of a photon by a point particle of mass m j , the differential cross section calculation follows in the same manner as in section 4.3 and resulting in The exchange m j ↔ E γ transforms the differential cross section to describe the scattering of a particle of mass m j from a photon target with energy E γ and as a result takes the form Equating (141) with (143) requires that Renormalization follows by construction because the structure added to the Standard Model is just N standard Yang-Mills gauge fields with SO(4) symmetry. Dimensional analysis of the Standard Model coupled to QGD shows that the connectons ω (i,m)ab µ have the same dimension as the U(1), SU(2) and SU(3) gauge fields. Therefore, all integrands in the Standard Model coupled to QGD, including those involving the connectons ω (i,m) have the dimension four or less which guarantees renormalization of all corrections to physical quantities [2].
If a solution also exists for the system of equations (134), (139) and (144), then the renormalizable Lagrangian (121), (122) and (123) correctly describes experimental observation without resorting to the graviton and instead using SO(4) gauge fields.
If a solution exists to the system of equations then, a solution also exists for the system of equations (134) We require the complex system of equations (145) with g j ∈ C N , g j = (g j1 , ..., g jN ) instead of (134) because solutions to (145)  To show existence of a solution to the system (145), (146) and (147) first note that (146) and (147) is a complex linear equation system in g i ∈ R and then assume that α i = 1. Our model contains 12 Fermions with j = 1, ..., J = 12, so that the linear system (146) has (12 × 13)/2 real equations and (12 × 13)/2 imaginary equations. The linear system (147) has 12 real and 12 imaginary equations. We can also represent the complex linear system (146) and (147) instead as a real system of 180 linear equations where g = (g 1 , ..., g N ) and A represents a square matrix with dimension N=180. Then, a solution to (149) exists whenever the determinant(A) = 0 or equivalently when the vector rows in A are not parallel among themselves. A simple inspection of the coefficients of A produced by (146) and (147) shows that when the complex vectors g j are not parallel among themselves, the different vector rows in A are not parallel. A solution to the complex system of equations (145) can be obtained for the set g j ∈ C 180 , j = 1, .., 12 following these steps: 1) assign a real random number between 0 and 1 to each Re(g ji ) and Im(g ji ), j = 1, ..., 12; i = 1, ..., N, 2) use the Fletcher-Reeves-Polak-Ribiere algorithm to minimize the equation Im(g ji )Re(g ki ) + Im(g ki )Re(g ji )) 2 .
This method produces a minimum of zero for this equation which amounts to solving (145). Note that the first squared term in (150) is the real part of (145) while the second squared term is the imaginary part of (145). It is not surprising that a minimum of zero exists for (150); there are 12 × (180 + 180) variables and only 88 equations. Furthermore, the vectors produced by this method are not parallel and explicit calculation of A shows that its determinant does not vanish 3 . Therefore, a solution can always be found to the system of 12 Fermions and one U(1) gauge Boson provided there are 180 connectons. This does not preclude the existence of solutions whenever N < 180, specially after abandoning the assumption that α i = 1.
The reader should be concerned by the number of parameters necessary to fit QGD coupled to the Standard Model when compared with General Relativity. However, two points should be kept in mind. First, the intention of this section was to show the existence of a solution rather than to find a solution with the smallest possible number of couplings. Given the non-linear nature of the system of equations, it is very likely that non-linear solutions with a significantly smaller number of couplings may very well exist; for example solutions obtained through the gradient method or simulated annealing. Second, QGD should only be compared with theories that successfully incorporate gravity and more importantly renormalize. Thus, when compared with heterotic strings, which have 496 gauge fields before even considering compactification, the number of gauge fields required by QGD is similar. Or, alternatively, like string theory, the gauge fields may represent some compactification of M-theory.
Nevertheless, we can cast QGD in a more symmetric manner where each of the 25 particles, whether a Fermion or a gauge Boson, couples to 25 different connectons; each particle carries charge with respect to 25 non-Abelian gauge fields; each particle connects to a different particle through a unique connecton, including one which self couples. Then in this construction a picture emerges where each of the 325 connectons mediates the interaction between two particular particles only (for example between a Higgs and a ν τ ). Then, the number, 25 × 13 + 25 × 13 = 650, of differential cross sections to O( v 2 c 2 ) is smaller than the number of couplings which, when omitting the α's, counts to 25 × 25 + 325 = 950: the first term counts the covariant derivative couplings and the second term counts the Yang Mills couplings.

Bending the Remaining Bosons
Experimental data evidencing the gravitational interactions does not yet exist for the W ± , the Z and the Higgs particles or the Gluons. However, we expect for them to become available in the future. Nevertheless, QGD can incorporate gravitational interactions for all these particles through the second term in (110). We present the propagators for the remainder of the bosonic sector r m, µ m ′ , ν where the first propagator describes the Gluons, the second describes the electroweak vector Bosons W and Z, and the third describes the Higgs in the SO(4) representation. Lagrangian (110) produces the following vertices: In general we expect that Gluon couplings satisfy g im = g i,SU (3) , m = 5, ..., 12; but, due to the mass difference between the W and Z particles, we expect g iW = g iZ and g jiW = g jiZ . These additional vertices produce the following tree-level and O(G) diagrams: These diagrams produce the T-matrix elements, which in turn produce the differential cross sections, and and whose derivation follows exactly along the lines in section 6 above. On the other hand the W, Z and Higgs Bosons are all massive and therefore we must compare their differential cross sections to those obtained from the Schroedinger equation with a standard Newtonian potential. These constraints yield the following conditions Finally, diagrams (o)-(q) above give the desired correction to O(G 2 ) for the W, Z, and Higgs particles interacting with Fermions: In the low energy limit, these diagrams produce the T-matrix elements, with B = W, Z, H. These T-matrices in turn produce the following differential cross sections, Comparison with the differential cross sections derived for a Schrodinger with a post Newtonian potential (138) imposes Inclusion of these differential cross sections further increases the number of differential cross section by 12+2×(12×3) = 84 for a total of 252 differential cross sections with a consequent increase in the total number of equations to be fitted. However, a solution can be shown to exist following the method in the previous section. We could go on and construct the O(G) amplitude for the scattering described by the exchange of a connecton between different Bosons as well as the O(G 2 ) amplitude for the corresponding 1-loop corrections. This requires modifying Ω (i,m) µab in (82) to read The requirement that the quadratic Lagrangian (109) remain invariant requires that the constant real symmetric matrices K The vertices (h)-(k) for Bosons emitting a connecton will then have a quadratic dependence on the matrices K (i) . Then, diagrams for Fermion-Fermion scattering remain unchanged; diagrams for Fermion-Boson scattering will produce tree-level T-matrices with quadratic dependence on the matrices K (i) ; their 1-loop corrections will produce T-matrices with quadratic and quartic dependence on the matrices K (i) ; diagrams for Boson-Boson scattering will produce tree-level T-matrices with quartic dependence on the matrices K (i) ; their 1-loop corrections will produce T-matrices with sextic dependence on the matrices K (i) . The coefficients of these matrices will be determined in the same manner as in the previous section: equating the differential cross sections of QGD to those obtained from the Schroedinger equation and gravitational lensing. Given the lack of experimental evidence, the task of writing each equation which reproduce General Relativity expectations is left as a curiosity for the reader. However we note that it may very well turn out that some processes not described by General Relativity, like those analogous to s-channel Bhabha scattering produced by the term ω (i,m) K mb K ma W b G a , are confirmed through experiment or that some processes described by General Relativity, like those produced by the term ω (i,m) K ma K ma G a G a , are confirmed to be absent through experiment. We stress that the matrices K (i) allow all these possible experimental scenarios.

Black Hole Entropy and Strings
QGD also explains the black hole entropy problem in the same way that strings do. To construct the map between QGD and strings we need first break down SO(4) into its isomorphic representation SU(2) × SU(2). We N=2 supersymmetrize at least one of QGD's SU(2) gauge groups and consider for example a K3 spacetime manifold. Each SU(2) Yang Mills theory then has BPS states which have the same moduli space as that of a supersymmetric sigma model which in turn can be used to calculate the degeneracy of BPS states for large energy states [13]. The degeneracy thus obtained determines the entropy which coincides with that of certain black hole solutions in string theory [13]. Thus QGD has states whose entropy is that of certain black hole solutions, and therefore, those states represent the quantum description of the black hole solutions in the exact same way they do in String theory.
The solution for SU(2) Yang-Mills found in [3] differs from the Schwarzschild black hole solution because the singularity of its horizon cannot be removed through a coordinate transformation. However, both singularities can be located at 2GM. Therefore in SU(2) Yang Mills, the horizon area and the black hole mass are related in the same way as for the Schwarzschild black hole: A = 16πM 2 . This result implies that the SU(2) Yang Mills black hole has the same entropy as the Schwarzschild black hole which in turn can be expressed in terms of BPS solutions. Therefore one has a dual representation of the same physical objects in the N = 2 supersymmetric limit: the macroscopic system is described by black hole solution found in [3], while the microscopic system is described by BPS states.

Conclusions
QGD, as presented here, reproduces the expected Newtonian and post Newtonian interaction evidenced by experimental measurements for both matter fields and photon fields. In addition to matching the differential cross sections from QGD to those of the Schroedinger equation with Newtonian and post Newtonian corrections as well as the differential cross section for the deflection of light by a point mass particle in General Relativity, QGD couples to the Standard Model in a straight forward and natural manner. Thus, QGD does not require us to re-invent the Standard Model to have a unified theory. This contrasts other attempts at describing quantum gravity which require a recasting of the already successful Standard Model.
The fundamental theory, a straight forward SO(4) Yang-Mills, couples to matter in the same manner as SU(3) and SU(2) theories describe the other fundamental forces and thus bears a high degree of symmetry with the nongravitational forces found in nature. Dimensional analysis of QGD and QGD coupled to the Standard Model show that all terms in the Lagrangians have dimensionality four or less and therefore renormalization follows [2]. QGD accommodates the full particle spectrum currently known and yields in the low energy limit the expected relativistic forces provided a sufficiently large number of SO(4) symmetries exists. The absence in QGD or QGD coupled to the Standard Model of the graviton should not concern the reader; after all, we have yet to detect such a particle.
The theory here proposed not only reproduces the observed experimental results of General Relativity, exhibits renormalization and couples in a straight forward manner to the Standard Model, it also incorporates yet to be observed, but highly expected, gravitational interactions between matter and the SU(3), SU (2) and Higgs Bosons as well as those between these Bosons. These features make QGD coupled to the Standard Model with SO(4) × U(1) × SU(2) × SU(3) local symmetry the most attractive theory to describes all interactions in nature.

Acknowledgments
I am most grateful to J. Lee for providing Ricci.m which was indispensable for all calculations. Also indispensable was the encouragement and funding by W. Israel during the early attempts to construct this model as well as funding by E. Braun in its latter stages.