
The phase diagram of Einstein-Weyl gravity
S. Silveravalle1,2and A. Zuccotti3
1Universit`a degli Studi di Trento,Via Sommarive, 14, IT-38123, Trento, Italy
2INFN - TIFPA,Via Sommarive, 14, IT-38123, Trento, Italy
3Ghent University, Technologiepark-Zwijnaarde 126, Be-9052 Gent, Belgium
Thanks to their interpretation as first order correction of General Relativity at high energies,
quadratic theories of gravity gained much attention in recent times. Particular attention has been
drawn to the Einstein-Weyl theory, where the addition of the squared Weyl tensor to the action
opens the possibility of having non-Schwarzschild black holes in the classical spectrum of the theory.
Static and spherically symmetric solutions of this theory have been studied and classified in terms of
their small scales behaviour; however, a classification of these solutions in terms of the asymptotic
gravitational field is still lacking. In this paper we address this point and present a phase diagram of
the theory, where the different types of solutions are shown in terms of their mass and the strength
of a Yukawa-like correction to the gravitational field. In particular we will show that, in the case of
compact stars, different equations of state imply different Yukawa corrections to the gravitational
potential, with possible phenomenological implications.
I. INTRODUCTION
General Relativity is one of the most successful theory
of last century, however there is still no general consensus
on how gravity should be described at the quantum level.
It is known that the Einstein-Hilbert action of General
Relativity is not renormalizable within standard pertur-
bative methods, and modifications of such action are ex-
pected at high energies. The study of modified gravity
theories is largely used to address the high energy limit
of gravity while preserving General Relativity at lower
energies.
The first corrections expected are quadratic terms in the
curvatures [1], which in 4 dimensions modify the action
as
SQG =Zd4x√−g[γR +βR2−αCµνρσCµνρσ +Lm].
(1)
Indeed such action has a long history: in [2] it is proved
that the action (1) in the vacuum is perturbatively renor-
malizable, in [3] it appears in the low energy limit of
string theory, and recently it has emerged in the frame-
work of the renormalization group flow [4–6], Asymptot-
ically Safe program [7, 8] and fakeons theory [9, 10]. The
physical content of (1) can be resumed in the standard
massless graviton, plus a massive scalar mediator and a
massive spin two ghost mediator, which at the quantum
level implies the loss of unitarity. While various authors
proposed solutions to the ghost problem, in this work we
want to focus on the classical content of the quadratic ac-
tion. In particular we are interested in studying the case
of static spherically symmetric solutions without a cos-
mological constant in order to describe the gravitational
field of isolated objects. Despite the classical solutions
of (1) have been largely investigated in recent works [11–
15], given the non linear nature of the field equations,
an exact form of the general solution is lacking, and nu-
merical methods have to be used in order to understand
the link between the asymptotic field and the physical
nature of the solution. In the following we consider the
Einstein-Weyl action
SEW =Zd4x√−g[γ R −α Cµνρσ Cµνρσ +Lm],(2)
i.e. the quadratic theory restricted to the β= 0 case.
Both the theories defined in (1) and (2) present a large va-
riety of solution families in addition to the Schwarzschild
one, but while the effect of the R2term has been largely
studied in astrophysical and cosmological contexts [16–
19], the C2term received relatively less attention un-
til recent times. The Einstein-Weyl theory indeed gives
interesting insight on the physical content of quadratic
gravity. It is proved that the solution space of Einstein-
Weyl gravity coincides with the one of (1) with the con-
straint R= 0, and a no-hair theorem presented in [20],
and corrected in [14], states that under certain condi-
tions the Ricci scalar must vanish. As main consequence
this theorem implies that all the asymptotically flat black
hole solutions of quadratic gravity are present only in
the Einstein-Weyl restriction. Moreover, the results in
[14] show that, together with black holes, the main new
families of solutions appearing in quadratic gravity, also
appear in the Einstein-Weyl theory. Finally, it has also
been shown that the C2has a much stronger effect than
the R2one in the properties of compact stars [21].
While attempts to link the different families of solutions
to the asymptotic field have been made in the full the-
ory [22, 23], the presence of two massive mediators in
the full quadratic case brings numerical issues when in-
tegrating the asymptotic field, namely some non-linear
terms can be larger than linear ones, which affects the
physical results obtained. Such numerical instabilities
have not been encountered in the Einstein-Weyl theory,
so we found sensible to restrict our results to this case
where the numerical procedure is much more reliable. In
this paper we use the analytical approximation found in
previous works, together with numerical shooting tech-
niques, in order to show the complete solution space of
the Einstein-Weyl gravity in form of phase diagram of
arXiv:2210.13877v1 [gr-qc] 25 Oct 2022