Exorcising Ghosts in Induced Gravity

Unitarity of scale-invariant coupled theory of higher-derivative gravity and matter is investigated. A scalar field coupled with dirac fermion is taken as matter sector.Following the idea of induced gravity Einstein-Hilbert term is generated via dynamical symmetry breaking of scale-invariance. The renormalisation group flows are computed and one-loop RG improved effective potential of scalar is calculated. Scalar field develops a new minimum via Coleman-Weinberg procedure inducing Newton's constant and masses in the matter sector. The spin-2 problematic ghost and the spin-0 mode of the metric fluctuation gets a mass in the broken phase of theory. The energy-dependence of VeV in the RG improved scenario implies a running for the induced parameters. This sets up platform to ask whether it is possible to evade the spin-2 ghost by keeping its mass always above the running energy scale? In broken phase this question is satisfactorily answered for a large domain of coupling parameter space where the ghost is evaded. The spin-0 mode can be made physically realisable or not depending upon the choice of initial parameters. Induced Newton's constant is seen to vanishes in ultraviolet. By properly choosing parameters it is possible to make the matter fields physically unrealisable.


I. INTRODUCTION
Finding a well-defined and mathematically consistent theory of quantum gravity is one of the most important problems of theoretical physics.Moreover, finding experimental evidence validating or falsifying one is equally hard.Presently there are several models of quantum gravity which are aimed at studying the quantum nature of space-time and investigating physics at ultra-high energies.Recently a minimalistic model in the framework of four dimension quantum field theory (QFT) in lorentizian space-time was investigated, which was shown to be renormalizable to all loops [1,2], and was recently shown to be unitary [5,6] (see also references therein).This then offers a sufficiently good and simple model of quantum field theory of gravity whose arena can be used to investigate physics at ultra-high energies.
Here in this paper, motivated by the results of [2][3][4][5][6][7][8] we study the scale-invariant higherderivative gravitational system coupled with matter fields.These constitute interesting systems.The scale-invariant purely gravitational sector consist of only dimensionless couplings.This makes the theory perturbatively renormalizable to all loops in four space-time dimensions by power-counting [9,10] (for classical picture of these theories see [1,11]).Coupling this with scale-invariant matter sector doesn't change the picture.The resulting theory is still perturbatively UV renormalizable in four dimensions due to lack of any dimensionful parameter.Classically the matter sector however has local conformal invariance, which is broken under quantum corrections due to conformal anomalies (local and non-local) [12,13].The interesting thing to note here is that in this quantum theory the counter terms generated still possess scale invariant structure (due to lack of any dimensionful parameter in theory) [14].This therefore preserves the renormalizability of theory [2], even though trace anomalies are present.
Scale-invariant gravitational systems coupled with matter have been investigated in past.Some of the first studies were conducted in [4,9,10,15,16], where the renormalisation group running of various couplings was computed and fixed point structure was analysed.Further investigation for more complicated systems were done in [17][18][19][20][21][22] (see also the book [23]).Matter coupling with conformal quantum gravity along with Gauss-Bonnet term were investigated in [24,25].Recently it has gained some momentum and these models have been reinvestigated [26,[30][31][32][33].The purpose in these papers were to see if it is possible to generate a scale dynamically starting from scale-invariant system.In [26] the authors called their model 'Agravity', where Planck scale is dynamically generated from the vacuum-expectationvalue(VeV) of potential in Einstein frame (not in Jordan frame).By this they achieve a negligible cosmological constant, generates Planck's mass, and addresses the naturalness [26,27] and inflation [28], but unitarity issues was not explored 1 .In [30][31][32][33] the authors studies the issue of dynamical generation of scale via dimensional transmutation in the presence of background curvature.This also induces Einstein-Hilbert gravity and generates Newton's constant, but unitarity problem was not addressed.An interesting idea has been suggested in [34,35] by assuming an analogy with QCD, where the authors addresses the problems of ghost and tachyon using the wisdom acquired from non-perturbative sector of QCD, as is argued that the gravitational theory enters a non-perturbative regime below the Planck scale.
The ideas of induced gravity goes long back.It was first proposed in [36,37], where the quantum matter fluctuations at high energy generate gravitational dynamics at low energy inducing cosmological and Newton's gravitational constant.Another proposal suggested in [38][39][40] induces Einstein gravity spontaneously via symmetry breaking along the lines of Higgs mechanism.Later in [41][42][43][44][45] the idea of generation of Einstein gravity via dynamical symmetry breaking was considered, following the methodology of Coleman-Weinberg [46].In [45], metric fluctuations were also incorporated in the generation of induced Newton's constant.Around the same time induced gravity from weyl-theory was also studied [47][48][49].Phase-transitions leading to generation of Einstein-Hilbert gravity due to loop-effects from conformal factor coupled with scalar field were studied in [50].In [51,52] the renormalization group improved effective-potential of the dilaton leads to running of VeV inducing masses (along with Einstein-Hilbert gravity).Furthermore the authors make a proposal along lines of [3,4] to tackle ghost and tachyons.Cosmological consequences of these models were explored in [53,54].
In this paper we therefore study scale-invariant gravitational and matter coupled systems in (4 − ǫ) dimensional regularisation scheme.The beta-functions are computed and compared with the past literature.The one-loop RG improved effective potential for the scalar field is computed by incorporating the quantum fluctuations of both matter and gravity [55].The scale invariance is broken dynamically when the scalar-field φ acquires a VeV via Coleman-Weinberg mechanism [46].This in turn induces gravitational Newton's constant, cosmological constant and masses in the matter sector.We work in lorentzian signature and take the sign of C 2 µνρσ (Weyl tensor square) to be negative, keeping the sign of R 2 term (where R is Ricci-scalar) to be always positive (this is done to avoid tachyonic instabilities).These choice of signs further allows necessary convergence in feynman +iǫ-prescription by suppressing fields modes with large action in the lorentzian path-integral.The sign of the Rφ 2 term is taken negative, so to generate the right sign of Newton's constant and to avoid tachyonic instabilities in the broken phase.The negative sign of C 2 µνρσ term is taken in order to satisfy the unitary criterion as stated in [5,6].In this case we no longer have asymptotic freedom as has been observed in euclidean case in [9,10,16,26,[30][31][32][33].The VeV generated in RG improved effective potential has running, and therefore induces running in Newton's constant and masses of matter fields.Due to the generation of Einstein-Hilbert term in the action, the propagator of metric fluctuations in the broken phase gets modified, and the modes acquires mass.In this broken phase we investigate the problem of spin-2 ghost by probing the running of its mass along the lines of [5,6] 2 .
The outline of paper is following.In section II the divergent part of the effective action is computed in (4 − ǫ) dimensional regularisation scheme.The beta-function are obtained from it.In section III one-loop renormalisation group improved effective potential for the scalar field is computed by incorporating quantum corrections from gravitational and matter degrees of freedom.In section IV the breaking of scale-invariance is studied via Coleman-Weinberg mechanism, which in turn induces gravitational Newton's constant and masses in the broken phase.RG equation for VeV is derived, which induces a running in the generated Newton's constant and masses.In section V a prescription to avoid spin-2 massive ghost is given, where a procedure to pick the right set of initial conditions is stated.In section VI numerical analysis is done to give evidence showing that there exist a large domain of coupling parameter space where spin-2 massive ghost can be made physically unrealisable.Finally in section VII conclusions are presented.

II. RG RUNNING
In this section we compute the renormalisation group running of the various coupling parameters present in our scale-invariant theory.We start first with the formalism and compute the various diagrams that contain UV divergences.These are then used to write the beta-function of the various coupling parameters.We start by considering the pathintegral of the coupled system ( = c = 1) where S GR and S matter are given by where the coupling parameters f 2 , ω, λ, ξ and y t are all dimensionless, and the geometric quantities (curvature and covariant-derivative) depends on metric γ µν .In the fermionic part of action, the Dirac gamma matrices are defined via tetrads and inverse tetrads (γ µ = e a µ γ a , γ µ = e a µ γ a ), and D µ is the spin-connection covariant derivative.
Here greek indices denote the space-time index, while the latin indices denote the internal lorentz index, and The internal indices are raised and lowered using internal metric η ab .For torsionless manifold the spin-connection can be expressed in terms of the christoffel connection Γ µ α ν (which can be re-expressed again in terms of tetrads) as, The dimensionless nature of coupling f 2 and ω/f 2 results in a fully dimensionless scaleinvariant coupled action.
We study the diffeomorphism invariant action of the coupled system using background field method [70,71] in (4 − ǫ) dimensional regularisation scheme.It is advantageous, as by construction it preserves background gauge invariance.The field is decomposed into background and fluctuation.Keeping the background fixed the path-integral is then reduced to an integral over the fluctuations.The gravitational metric field is decomposed into background and fluctuation, while the tetrads (and its inverse) are expressed in powers of this fluctuation field.The matter fields are similarly decomposed.The gauge invariance of the full metric field is then transformed into the invariance over the fluctuation field.To prevent over-counting of gauge-orbits in the path-integral measure, a constraint is applied on this fluctuation field, which results in appearance of auxiliary fields called ghosts.The effective action generated after integrating over the fluctuation and auxiliary fields still enjoys invariance over the background fields.
The quantum metric is written as where ḡµν is some arbitrary (but fixed) background and h µν is the metric fluctuation.The full action can be expanded in powers of h µν .The path-integral measure over γ µν is then replaced with measure over h µν .Integrating over the fluctuation field implies that in some sense they will appear only as internal legs and never as external legs.The background gauge invariant effective action formalism allows to choose a particular background for the ease of computation.In particular writing ḡµν = η µν + H µν (while still keeping H µν generic) allows one to use the machinery of the flat space-time QFT, thereby giving a particle interpretation to the internal (h µν ) and external (H µν ) legs, in the sense that the former behaves as virtual particle, while the later is the corresponding external particle.In this manner one can compute the effective action for the external leg H µν .Alternatively one can expand the full action around flat space-time directly, calling the fluctuation to be h ′ µν though.This is a highly non-linear gauge theory with infinite number of interactions terms (however their couplings are related to each other by diffeomorphism invariance).Then following the usual strategy of background field method and writing h ′ µν = H µν + h µν , it is quickly seen that H µν is the external leg corresponding to the virtual particle given by h µν .Integrating over quantum fluctuations h µν then gives the effective action in terms of H µν field.
One can then set-up Feynman perturbation theory by expanding the original action in powers of h µν and H µν .Similarly writing the scalar and fermion fields as one can expand the action in powers of fluctuations χ, η and η.The piece which is quadratic in only fluctuations (h µν , χ, η and η) gives the propagator while all the other terms gives the interactions vertices.In one-loop approximation the terms which are exclusively quadratic in fluctuations are retained, all other terms which involve higher powers of fluctuations (h µν , χ, η and η) contribute in higher-loops and will be ignored here.In-fact for computing the running of all matter couplings (except ξ) it is sufficient to consider the situation with H µν = 0.However, in the case for computing running of ξ, terms up-to (and at-least) linear in H µν should be retained.Similarly if one is interested in studying behaviour or R 2 and R µν R µν then one should at-least retain terms up-to quadratic in H µν .

A. Gauge Fixing and Ghosts
The path-integration over the gravitational field is ill defined.This is a general feature of gauge theories where the gauge invariance (diffeomorphism invariance for gravity) relates two field configuration by gauge transformation.Such field configuration will contribute equally to the path-integral.However this will lead to over-counting.To prevent such over-counting, gauge-invariance needs to be broken by constraining the gauge field.This procedure of systematically applying the constraint leads to ghost, which are elegantly taken care of by the Faddeev-Popov prescription [72].
However in this style of breaking the invariance one may wonder whether the gauge (or diffeomorphism) invariance emerges in the effective action.To make sure that the effective action obtained after integrating out the fluctuation field is gauge invariant, background field method is followed.It is a method (and procedure) which guarantees that the effective action constructed using it will be background gauge invariant.Below we describe the procedure for gauge-fixing in the background field method.
The diffeomorphism invariance of the full action in eq (2) implies that for arbitrary vector field ǫ ρ , the action should be invariant under the following transformation of the metric field variable, where L ǫ γ µν is the Lie derivative of the quantum metric γ µν along the vector field ǫ ρ .Decomposing the quantum metric γ µν into background (ḡ µν ) and fluctuation (h µν ) allows one to figure out the transformation of the fluctuation field while keeping the background fixed.This will imply the following transformation of h µν .
where ∇ is the covariant derivative whose connection is constructed using the background metric.This is the full transformation of the metric fluctuation field.Ignoring terms which are linear in h µν allows one to investigate only the one-loop effects.These ignored terms are however mandatory when dealing with higher-loop effects.The invariance of the action is broken by choosing an appropriate gauge-fixing condition implemented via Faddeev-Popov procedure.
The gauge fixing action chosen for fixing the invariance under the transformation of the metric fluctuation field is given by, where α and ρ are gauge parameters, while Y µν is either a constant or a differential operator depending upon the gravitational theory under consideration.For the theory considered here in eq. ( 2), we consider higher-derivative type gauge fixing with Y µν = (−ḡ µν ¯ + β ∇µ ∇ν ), where ¯ = ∇µ ∇µ .Choosing ρ = −1 and β = 0 correspond to Landau gauge condition.Taking α → 0 imposes the gauge condition sharply.
The ghost action is obtained following the Faddeev-Popov procedure [72].In general if the gauge-fixing condition on the gravitational field h µν is written as F µ = 0 (which here is F µ = ∇ρ h ρµ − 1+ρ d ∇µ h), we introduce it in the path-integral by multiplying the later with unity in the following form, where F ǫ µ is the gauge transformed F µ .As Y µν contains derivative operator, therefore its determinant is non-trivial.The original path-integral (without gauge-fixing) being invariant under transformation eq. ( 8) of the field h µν implies that a change of integration variable from h µν to h ǫ µν doesn't give rise to any jacobian in the path-integral measure.However replacing the measure over F ǫ µ with measure over ǫ ρ introduces a non-trivial jacobian in the path-integral.This is obtained as follows, In the background field formalism this jacobian consist of background covariant derivative, background and fluctuation fields, and is independent of the transformation parameter ǫ ρ .This implies that it can be taken out of the functional integral over ǫ ρ .Changing the integration variable from h ǫ µν to h µν , and ignoring the infinite constant generated by integrating over ǫ ρ , gives us the gauge fixed path integral including the determinant.
The functional determinant appearing in eq. ( 11) can be exponentiated by making use of appropriate auxiliary fields.Writing the functional determinant (det Y ) 1/2 as (det Y ) × (det Y ) −1/2 , allows to combine the former with the Faddeev-Popov determinant in eq. ( 11), which is then exponentiated by making use of anti-commuting auxiliary fields, while the later (det Y ) −1/2 is exponentiated by making use of commuting auxiliary fields.The former auxiliary fields are known as Faddeev-Popov ghosts, while those in later case are known as Nielsen-Kallosh ghosts [73,74].The path integral of the full ghost sector is given by, where Cµ and C ν are Faddeev-Popov ghost fields arising from the gauge fixing in the gravitational sector, and θ µ is the commuting ghost arising due to fact that Y µν contains derivatives.
In the case when F µ is given as in eq. ( 9), the Faddeev-Popov ghost action is given by, where, Here the last two lines contains terms linear in h µν .These are not relevant in doing oneloop computations, but at higher-loops they are important.In the following we will ignore ghost contributions completely as they are not relevant in the computation of the running of matter couplings, while the running of gravitational couplings are taken from past literature [5,6,16,26,78].

B. Gravitational Field Propagator
The propagator for the gravitational field is obtained by expanding the gravitational action around the flat space-time up-to second order in the fluctuation field h µν .By decomposing the fluctuation field in terms of various components and writing them using the ) as given in appendix B, we note that this second variation can be expressed in a neat form in momentum space in the following manner, Moreover the gauge-fixing action can be similarly expressed by using the projection operators, By writing the gauge-fixing action in terms of the projection operators allows us to see clearly which modes of the field are affected by the gauge-fixing.For example the spin-2 mode is not affected at all by the gauge-fixing condition.Interestingly it should be noted that there is another gauge-invariant mode of the field which arises due to the action of spin-s projection operator on the h µν field (see appendix B).However under harmonic type gaugefixing condition this mode doesn't remain completely unaffected.Only for some particular gauge choices this mode is not affected by the gauge-fixing condition.Landau gauge being one such choice ρ = −1, β = 0, α = 0.In this gauge choice only the purely longitudinal modes are gauge fixed.In this gauge the propagator for the metric fluctuation field is the following, where Y i are the propagators for the various spin-components: Here ∆ µναβ G is the inverse propagator for the h µν field including the gauge fixing and is symmetric in µν and αβ.As the propagator is 1/q 4 , it doesn't allow to be decomposed further via partial fractions.Here the first term in eq. ( 17) arises due the presence of C 2 µνρσ part of action, while the later comes from the R 2 part.In this form it is not clear how the unitarity will be satisfied.

C. Formalism
We employ the background field formalism and decompose the metric and matter fields as in eq.(5 and 6) respectively, where we choose the background space-time to be flat.In order to do the one-loop computation the action is expanded up-to second powers of the all fluctuation field (h µν , χ, η and η).This will result in various vertices and propagators that are required for the one-loop analysis.The second variation of the matter action is given by the following, The various vertices and matter propagators are written in detail in appendix C. Having obtained the second variation giving propagator and the vertices, we set forth by considering the path-integral over the fluctuation fields.In this case the zeroth order term will be independent of the fluctuation fields and can be taken out of the path-integral.The linear term can be removed by doing field redefinition.In general, terms proportional to equation of motion can be removed by doing field redefinition.Such a redefinition will give rise to a trivial jacobian from the functional measure.The quadratic piece can now be investigated easily by putting together all the field fluctuations in the form of a multiplet Φ T = (h µν , χ, η T , η).Using this the path-integral can be written in a more compact form as, where J = {t µν , t, ρ, ρT } is the source multiplet which couples with the fluctuation field multiplet Φ = h µν , χ, η T , η .The super matrix M is given by From the generating functional Z, one can define the one-particle-irreducible (1PI) generat- and Φ T is the expectation value of Φ T field.The 1PI generating functional is also the effective action containing the quantum corrections.In the one-loop approximation (which we are considering here), one can perform the functional integral over the super-field Φ thereby giving an expression for the one-loop effective action to be, where the first two terms correspond to tree-level diagrams while the last term contains one-loop quantum corrections.The appearance of generalised trace 'STr' means that In the following we will be computing the divergent pieces present in the STr ln M.There are various ways to compute the one-loop quantum corrections.The most common methodology to do is via Feynman diagrams after computing vertices and propagator.Here we will follow a slightly different strategy of computation via evaluation of functional determinant.We start by writing M = ∆ + V, where the former ∆ contains the various propagator while the later V contains various vertices.They are given by, ) Pulling out ∆ from the expression for M allows to expand the residual expression (I+∆ −1 •V) (where I is a generalised identity in super-field space) under the logarithm in a perturbative manner as follows, Here (26) It should be mentioned here that so far we took background metric to be flat with H µν = 0.This is enough to compute the counter-terms involving quantum gravity corrections to all matter couplings including their anomalous dimensions.If we had included terms with H µν = 0, then it is also possible to compute the counter-term proportional to Rϕ 2 .But here for simplicity we keep H µν = 0, and the counter-term proportional to Rϕ 2 will be computed using methodology of heat-kernels (HK) later.Heat-kernel method is quick, as the HK coefficients have already been computed in past [75][76][77][78].Besides, it also gives an alternative check on the computation done using feynman diagrams.For flat background the term STr ln ∆ is irrelevant, but it is not so if the background is non-flat for which case this gives purely curvature dependent divergent contributions.Such contributions have been computed elsewhere in literature [16,78] and here we will just take their results.
In the following we will be computing the various graphs that are giving quantum gravity corrections to the running of the matter couplings and the fields anomalous dimensions.

D. Graphs
Here we will be computing the various graphs that contain divergent contributions.These are basically the terms in the series expansion given in eq. ( 25), which will be evaluated one by one.The first term of the series contain tadpole graphs (those having single vertex), the second term of series has bubble graphs (those having two vertices), the third term of the series are triangle graphs containing three vertices while the fourth term of series are square graphs with four vertices.The series has infinite number of graphs, but the divergent ones are only present in the first four terms of the series expansion, and below we will be computing them.

Tadpole
These set of graphs arises from the first term of the series in eq. ( 25) which is STr(V) = STr(∆ −1 • V).Here the super-trace takes care of trace over bosonic and fermionic part, and includes the trace not only over field space but also over lorentz indices.This will imply the following, Here the first term contains graphs having an internal graviton line, while the next three terms contains the usual diagrams which are present without gravity.The former gives quantum gravity contribution.The set of graphs present in the tadpole order are shown in figure 1.
Each of these diagrams can be evaluated using the vertices given in appendix C 2. Here we will write their contribution.However, the last three terms in eq. ( 27) vanish in scaleinvariant theory.The gravitational ones are complicated and lengthy as the vertices are cumbersome.Below we write this For the other one there is more algebra as it involves Dirac-matrices.Here we will write the expression after performing the lorentz and Dirac matrix algebra.This is given by, The momentum integrals can be evaluated in the (4 − ǫ) dimensional regularisation scheme and the divergent piece can be singled out easily.The divergent piece of all the above tadpole contribution is, where M 2 /Z = 16πf 2 is introduced for convenience.

Bubble
These set of graphs arise from second term in eq. ( 25) which is Here again the super trace is evaluated as before.This will imply, Here each of the term will contain several diagrams, but only few contain the divergences that are relevant for our purpose.These diagrams contain two vertices.They can be classified in three categories: (a) those with two internal graviton lines, (b) those with one internal graviton and one internal matter line and (c) those with two internal matter lines.The set of diagrams are shown in figure 2.
Each of these diagrams can be evaluated using the vertices given in appendix C 2. The super-trace given in eq. ( 31) contains lot of diagrams, but not all contain UV divergence.Here we will mention only the ones having the UV divergences.These come from, From the various terms written in eq. ( 32) the first two lines contain diagrams which are purely matter ones and correspond to the diagrams shown in first row of figure 2, while the last two lines contain diagrams having quantum contributions and correspond to the diagrams shown in last two rows of the figure 2. The trace is over the Lorentz and spacetime indices.After performing the algebra over Dirac matrices and doing the contraction of the tensors we get the simplified expression involving momentum integrals.The divergent contributions of the purely matter diagrams is, The diagrams containing the internal graviton legs are bit complicated, as it involves lengthy Dirac matrix algebra and tensor manipulations.For doing these we have used various tricks to extract the divergent piece and also used MATHEMATICA packages (xAct [79], xTras [80] and FEYNCALC [81]).Below for simplicity we will mention only the divergent piece of these diagrams to evade unnecessary complications, while some of the details will be mentioned in the appendix.The diagrams with internal graviton line has the following contributions Here the first row contains contributions to the scalar sector, while the second row contain contributions to the fermion sector.The former correspond to diagrams of the second row in figure 2, while the later correspond to diagrams in the third row of figure 2 respectively.Putting together the full contribution of the bubble kind of diagrams, we get contribution to the one-loop effective action of the diagrams having two vertices.This is given by,

Triangular Graphs
These diagram are generated from the third order terms in the series of eq. ( 25), 1/3STrV 3 where STrV These diagrams have three vertices and consist of graphs which are either purely matter oriented or ones which include a mixture of matter and gravity.These graphs give correction to vertex: either to ϕ 4 or to yukawa vertex ϕ θθ.On expansion we will see that there are many diagrams but we will consider those which carry divergent pieces and give correction to running couplings.Here we will only mention the trace pieces which will be carrying the divergences, however in principle there will be many more diagrams.The relevant terms in the trace which will be of relevance can be guessed by looking at the set of third order diagrams in figure 3.These are given by, Here the first line correspond to set of diagrams of purely matter type, while the rest of terms contain quantum gravity corrections.The second line contribute to running of ϕ 4 coupling while the rest of the terms correspond to the running of yukawa coupling.Interestingly while doing the computation involving fermions it is noticed that not all of terms are non-zero.
The divergent contributions of these diagrams and corresponding their contribution to the effective action is given by, FIG. 4: Various diagrams containing divergences with four vertices.Here there are only two diagrams.The first one is purely matter while the second one contain quantum gravity correction.Both give correction to the ϕ 4 coupling.

Square Graphs
These set of contribution arise at fourth order of the perturbative expansion given in the series eq. ( 25) and comes from −1/4STrV 4 .Here again the super-trace is given by, These diagrams consist of four vertices and all of them contribute to the running of ϕ 4 coupling.There are only two diagrams at this order.One is purely matter type, contains a fermion loop with four external scalar legs, while the other one contains quantum gravity correction.The trace can be expanded as before and consists of large number of graphs but the ones containing the divergences are only two.These are the following, Here the former term is purely matter and contains a fermion loop, while the second term contain virtual gravitons.Considering as before just the divergent part and their corresponding contribution to the effective action, we find that, Here we compute the divergence proportional to Rϕ 2 .There are two ways to compute it.The first is via computation of feynman graphs and second via heat kernel.Conceptually both are same and give same results, however the later is quicker.In each case we break the metric γ µν appearing in path-integral is written as in eq. ( 5).In the former the background is further written as η µν + H µν (see also paragraph following eq.( 5)).Here H µν will act as external graviton for the corresponding internal leg h µν .The action is then expanded under this decomposition.This way we get additional vertices.The vertices in the previous section are the ones for which H µν = 0.If H µν = 0 then we get contributions which contains dependence on external graviton leg, and if there are derivatives acting on H µν , then those will give terms proportional to background curvature.This was employed in [26].
Alternatively, one can take the background metric ḡµν to be maximally symmetric and compute the contributions proportional to background curvature using Heat-Kernel.This will project directly the contribution proportional to Rϕ 2 .We use the heat-kernel methodology to compute the one-loop divergence proportional to Rϕ 2 .The matter fields are decomposed as in eq. ( 6), but this time we take background matter fields to be constant.The fluctuation metric h µν is further decomposed into various components as where This decomposition of h µν introduces jacobians in the path-integral, which can be cancelled by redefining the fields as Under this decomposition the Hessian for the fluctuation fields is obtained on a maximally symmetric background.This will be same as in eq. ( 20), except now M will be different.
The multiplet Φ also gets modified Φ T = h T T µν , ξµ , σ, h, χ, η T , η .As the background matter fields are constant, therefore the fermion sector is completely decoupled with mixing of fermion fluctuation field with metric and scalar fluctuations.This is not so when background matter is not constant.On the maximally symmetric background metric with constant background matter the matrix M is given by (in the Landau gauge ρ = −1) where We will not be considering the contribution for the ghost here, as they will not contribute at one-loop to the term proportional to Rϕ 2 .The one-loop effective action is given by, These functional traces can be tackled using heat-kernel [75][76][77][78].One can compute the divergent part of the effective action from this.Since the background matter fields are not constant therefore this will however not be able to give the anomalous dimensions of the matter fields.However, the anomalous dimension has already been computed earlier using feynman diagram therefore it will not be considered again here.Here we will just look the divergent contribution proportional to Rϕ 2 , which is given by, Here the first two terms arise due to matter loop while the rest of terms contain quantum gravity corrections.This is in agreement with [26].

E. Effective action and Beta Functions
Once we have computed all the relevant graphs at various order of the perturbation theory and their divergent contributions, it is easy to put them together to write the divergent part of effective action and collect all the pieces together.The divergent part of the full effective action is given by, Once the divergent part of the effective action is written it is easy to compute the betafunction of the various couplings by incorporating the effect of the wave-function renormal-isation.These are given by, where t = ln(µ/µ 0 ) (µ 0 is a reference scale) and d/dt = µd/dµ.Here Z ϕ and Z ψ is the wave-functional renormalisation of scalar and fermion field respectively, while η ϕ and η ψ is the corresponding anomalous dimension.The beta-functions obtained here agree fully with [26,27], while there is partial agreement with [23,[82][83][84][85][86].For completeness we also include the running of gravitational couplings which has been taken from past literature [16,78].These are given by, d dt We will be doing the RG analysis of the couplings and exploring the issue of unitarity later.

III. EFFECTIVE POTENTIAL
Here we compute the effective potential for the background scalar field ϕ which gets contributions not only from matter fields but also from gravitational sector.
To compute the effective potential for scalar, the background scalar field is taken to be constant.This is sufficient to compute the effective potential.The quantum gravitational fluctuations are considered around a flat background.The fermion fields are likewise decomposed into a constant background (which for simplicity is taken to be zero θ = θ = 0) plus fluctuations.This simplifies the computation very much.As the ghost action doesn't depend on the background scalar field ϕ, therefore there is no contribution by the ghost to the effective potential, and hence will be ignored in the following.Once the full second variation of the action is performed, we have the hessian needed to compute the one-loop effective potential.This can be obtained directly from eq. ( 45) by putting background R = 0 and replacing background covariant derivative with partial derivative.Being on flat background allows the freedom to work directly in momentum space.
Moreover, in flat space-time the transverse-traceless decomposition of h µν given in eq. ( 42) can be rewritten in an alternative form.In this new decomposition the field components σ and h of h µν are replaced by s and w.These new fields s and w are related to old ones in the following manner The advantage of doing this transformation is to bring out the scalar mode which remains invariant under diffeomorphism transformation stated in eq. ( 8).The field s is therefore gauge invariant, while the field w is longitudinal.So the decomposition of h µν has two gaugeinvariant components h T T µν and s, with two longitudinal components ξµ and w.Furthermore, on flat space-time one can also use the set of orthogonal projectors to project h µν on various components h T T µν , ξµ , s and w (see appendix B).In terms of new field variables, the hessian mentioned in eq. ( 45) can be rewritten (for R = 0) in a more transparent manner to see clearly the gauge-dependent and gauge-independent part.The hessian for h T T µν and ξµ remans same, while the mixing matrix of σ, h and χ gets rotated due to field transformation stated in eq.(57).The new mixing between the field variables s, w and χ has a simplified structure.Moreover, this transformation of field variable is unaccompanied by any non-trivial jacobian in the path-integral.The one-loop effective potential is therefore obtained from a simplified hessian, where the entries of the scalar mixing matrix are, From the entries of mixing matrix we clearly notice that S ss , S sw , S sφ doesn't depend on gauge parameter.The only gauge dependence is in S ww .For a generic case with an arbitrary field variable, the one-loop hessian can be written in the form (− − m 2 ) (where is m 2 contain background field contributions and couplings).
In this case the effective potential is given by the general formula where d s is the factor coming due to the degrees of freedom of the field.The term 3/2 in the bracket can be absorbed by rescaling µ 2 as μ2 = µ 2 e 3/2 .By exploiting this generic formula one can write the contribution to the effective potential from the various field modes of the metric fluctuation field, the scalar and the fermion field.In the case for spin-2, the differential operator responsible for the contribution can be factored and has the form where A 1 and A 2 are given by, respectively.Here both A 1 and A 2 are dimensionless.It should be noted that for positive λ, A 1 is negative while A 2 are positive.If the sign of ξ is reversed, the roles of A 1 and A 2 gets interchanged.A negative A 1 is tachyonic in nature.This is a source of instability in the effective potential and will give imaginary contribution to the effective potential.Plugging A 1 ϕ 2 and A 2 ϕ 2 for the m 2 in the expression for the effective potential in eq. ( 60) and summing the two, we get the contribution of the spin-2 mode to the effective potential.This imaginary piece though is an infrared effect.It is an indication that background chosen for doing the computation is not stable, and is a generic feature of gravity coupled with scalar field in flat space-time at zero temperature [19,23].This is same as the instability seen in the gas of graviton at finite temperature, an indication that flat space-time is unstable.This issue has been throughly investigated in past in [87,88].This kind of tachyonic mode will create issues with unitarity, but this one is different from the unitarity issue caused by the ghost of higher-derivative gravity, in the sense that the former is an IR problem and has no affect on the UV physics, while the later does affect the UV physics also.Since we are interested in sorting out the problem of unitarity caused by ghosts of higher-derivative, therefore we study only this by focusing on the real part of the effective potential, as the imaginary piece is relevant in IR and deals with tachyonic instability only.This is an important realisation as it decouples the two problems: (a) problem of unitarity caused by higher-derivative ghosts, and (b) problem of unitarity caused by tachyons.This paper deals with the former problem.
The contribution of the spin-0 mode is a bit complicated as it involves the scalar mixing matrix.We need to compute the determinant of this mixing matrix and then compute the effective potential of the operator so obtained from this matrix determinant.The operator obtained after matrix determinant is following, This operator is a cubic polynomial in − and will therefore have three roots.The nature of roots can be analysed using the discriminant ∆ of the equation formed by putting the cubic polynomial in (62) to zero.We write If ∆ > 0 then all roots are real, if ∆ = 0 then there is a multiple root, and ∆ < 0 then roots are complex.The operator in eq. ( 62) can be factorised as where B 1 , B 2 and B 3 are dimensionless.As the product of roots B 1 B 2 B 3 is positive and B 1 B 2 + B 2 B 3 + B 3 B 1 is negative, therefore this will imply that when ∆ > 0 then two roots will be negative.If ∆ = 0, then there is a multiple root with negative sign.In the case when ∆ < 0 there is a pair of complex conjugate root with negative real part and a positive real root.The cases with ∆ ≥ 0 has roots carrying negative sign, while for ∆ < 0 the complex conjugate pair has a negative real part.In all these cases the roots can be written as In the case when ∆ > 0, θ is positive and real, in ∆ = 0 case θ = 0, while in ∆ < 0 case θ is imaginary.The factor of −1 in the parametrisation of roots can also be exponentiated as e iπ .This factor is the source of tachyonic instability and will give rise to an imaginary contribution in the effective potential.This is similar to the instability caused in spin-2 case and is an indication that flat space-time is not a true vacuum [87,88].The contribution to the effective potential from the scalar sector is now easily computed using the generalised expression given in eq. ( 60).This is done by replacing m 2 in eq. ( 60) with B 1 ϕ 2 , B 2 ϕ 2 and B 3 ϕ 2 , and summing all together.Using the parametrisation for the roots written in eq. ( 64) and employing the properties of exponential functions, one can write the effective potential in simple terms 3 .The contribution of the fermions needs to be done in a different manner.It arises from −iTr ln(iγ µ ∂ µ −y t ϕ).This can be written in an alternative form by squaring the operator and by making use of the gamma-matrix properties.This then become −i/2Tr ln(− − y 2 t ϕ 2 ).Then using the generalised expression in eq. ( 60) one gets the contribution for the fermions.The full effective potential involves the tree-level contributions plus the quantum corrections.The one-loop RG improved full effective action is then given by, where A i 's, B 1 , r and θ are dimensionless and RG-time t dependent.When θ → iθ, cosh(2θ) → cos(2θ) and sinh(2θ) → i sin(2θ), thereby preventing the switching between real and imaginary part.
In the following we will study the real part of the effective potential.We ignore the imaginary part, as the imaginary part arises from the tachyonic modes of theory and is relevant in IR.We are interested in investigating unitarity issues caused by higher-derivatives ghosts.It should be noticed that the effective potential still posses the Z 2 symmetry, as ϕ 2 = 0 is an extrema.But due to radiative corrections the real part of quantum corrected effective potential might develop a vacuum expectation value (VeV) away from zero.

IV. SYMMETRY BREAKING
Due to RG corrections a VeV is generated in the effective potential, which then becomes the new vacuum.The original ϕ 2 = 0 vacuum becomes unstable under RG corrections and the field migrates to the new vacuum which occur at ϕ 2 = κ 2 .It is given by, At the tree level our original action of the theory is scale-invariant and there is no massparameter to begin with.However the mass parameter enters the system via RG running thereby breaking scale-invariance.The effective potential so generated not only breaks scale invariance but also breaks the Z 2 symmetry.The generation of VeV consequently gives mass to scalar and fermion fields.It also generates an effective Newton's constant, beside generating newer interactions.The generated mass and Newton's coupling can be expressed in terms of VeV κ 2 and all the other couplings as The generation of mass and Newton's constant makes the propagators for various fields massive.In particular the graviton propagator after the symmetry breaking is following, where now G is the induced Newton's constant and is defined using eq.( 67).The masses M 2 and M 2 /ω are given by The interesting thing about the generation of Newton's constant is that now as the propagator becomes massive, so there is a spin-2 massive ghost that appears in the system, which in the original theory was massless.In the original theory we cannot do the partial-fraction trick in the h µν propagator, which is possible in the broken phase due the induced Newton's constant G.Not only the spin-2 ghost becomes massive but also the scalar mode acquires mass through symmetry breaking.We call this massive scalar mode 'Riccion'.It should be pointed out that if we had taken ξ to be negative, then there will be tachyons in the broken phase.So the presence of higher-derivatives terms and requiring no tachyons to be generated in broken phase fixes the sign of ξ.This also generates right sign for induced Newton's constant.The sign of various couplings in the broken phase is then in accordance with the sign of parameters taken in [5,6].
At this point we compare the propagator of metric fluctuation written in eq. ( 17) with the one written in eq. ( 68).The former is before symmetry breaking while the later is after symmetry breaking.The former has no mass, while later contains masses.Although the appearances of the two are different one should however be careful while counting the propagating degrees of freedom in the two.In the later case (broken phase) it is easy to count: two massless graviton modes, five massive-tensor ghost modes and one massive scalar mode, thereby making eight propagating degrees of freedom.In the former case (unbroken phase) one should count carefully.The pure C 2 µνρσ -gravity has six massless propagating degrees of freedom [89,90].For pure R 2 gravity, the theory two massless propagating degrees of freedom as the linearised field equation (∂ µ ∂ ν − η µν ) h = 0 shows to have fourth-order time derivatives.Thereby totalling the degrees of freedom in unbroken phase to be eight.This implies that the propagating degrees of freedom in both phases is same, except in broken phase some of the modes acquire mass due to symmetry-breaking.
The generation of mass for the spin-2 ghost and scalar-mode gives us a hope to investigate unitarity by using the criterion stated in [3,[5][6][7][8].In the RG improved effective potential the VeV has t-dependence.This arises because at each energy scale the effective potential has a VeV.This translates into t-dependence of VeV.The RG running of VeV depends on the running of the other couplings in a complicated manner.This running of VeV then translates into running of generated Newton's constant.The running of the VeV κ 2 can be computed using the expression of the real part of effective potential given in eq. ( 65).When ϕ 2 = κ 2 , then we are at the minima.The minima condition written in eq. ( 66) then translates into following, d dϕ 2 Re(V eff ) where, As κ 2 = 0 and Z φ = 0 therefore these overall factors goes away and the residual condition simplifies to the expression in the square bracket written in eq. ( 70).As κ 2 (t)/μ 2 is dimensionless, we call it K(t).One can then directly solve for K(t) using eq.( 70) in terms of all couplings of theory.This also gives the flow of K which is generated due to the flow of various couplings present in the theory.We however take a t-derivative of the expression in the square-bracket of eq. ( 70) to compute the beta-function of the K(t).This is needed in checking and locating extrema of K(t).Such extrema are crucial points as will be seen later.
This is a linear first order differential equation for the ln K(t).Plugging the running of various couplings from eq. (52,53,54,55,56) in RHS of eq. ( 72) we get the beta function of K(t).This will be a very complicated function of various couplings.Using the running of K(t) we can compute the running of the effective Newton's constant by exploiting the expression for induced G mentioned in eq. ( 67).This is given by, In order to investigate the issues of unitarity caused by higher-derivative we consider the following combination M 2 /µ 2 .We first note that this is Taking t-derivative of this yields, Similarly the expression for induced M 2 /ωµ 2 is, and the RG flow of this combination is given by, The generation of VeV also induce masses for the scalar and the fermion fields, which is mentioned in eq. ( 67).Due to the running of VeV, these masses inherits a running.Then to investigate whether these fields are physically realisable or not, we consider the flow of combinations m 2 s /µ 2 and m 2 f /µ 2 .The running of these combinations is given by, V. UNITARITY PRESCRIPTION In this section we dictate the algorithm to choose the set of initial conditions for which the theory will have a unitary flow.We start by analysing the RG equations given in eq.(52,53,54,55,56).The first thing we do is to extract the running ω using the eq.(55 and 56).This is given by, From this running we notice that as the RHS is always positive therefore ω is a monotonically increasing function of t.In [5][6][7][8] it was shown that in order to avoid tachyonic instability, we should demand that ω ≥ 0.Here in the present scale-invariant theory we should demand the same.This is done in order to prevent the occurrence of tachyons in the broken phase.
For every value of ξ, ω will have two fixed points.
ω 1 is repulsive while ω 2 is attractive.For ξ small both these fixed points lie in the unphysical tachyonic regime.For large ξ the fixed points are complex conjugate with negative real part.
Since ω is monotonically increasing with t, therefore one can alternatively study the RG flows of various parameters in terms of ω.Prevention of tachyonic instability restricts ω to lie between zero and infinity.This then serves as a good candidate in terms of which the the RG flows can be analysed.In [5][6][7][8] the RG flows were studied in terms of ω.
The crucial problem in overcoming the issue of unitarity is to choose the right set of initial conditions so that throughout the RG evolution the flow remains unitary in the sense that the ghost mass remains always above the energy scale, and the effective potential doesn't develop any further instability (other than the ones already present).To prevent the occurrence of this instability requires that the coupling λ remains positive throughout the RG flow (as negative λ will result in tachyonic instability for scalar field φ).This particularly depends on the choice of initial condition for yukawa coupling.If the yukawa coupling is above a certain threshold then λ becomes negative too soon during the RG evolution, making the effective potential unstable.In standard model of particle physics this is an important instability problem where the electroweak vacuum becomes metastable [91] (see references therein).In present case of scale invariant gravity, we have freedom to explore the set of initial conditions which will give unitary evolution.So we just consider those domains where this instability can be avoided.
In [5][6][7][8] it was observed that the RG evolution of M 2 /µ 2 is such that its flow has a unique minima.This was a crucial feature which allowed us to seek those RG trajectories for which this minima is above unity.These RG trajectories are the ones for which the flow is unitary (massive tensor mode is not physically realisable).In the present case of gravity being induced from scale-invariant theory we seek a similar behaviour of induced M 2 /µ 2 , where now M 2 is given by eq. ( 69), and the flow of M 2 /µ 2 is given in eq. ( 75).The flow of M 2 /µ 2 is much complicated in the present case and it is difficult to give a rigorous analytic proof that there exist a minima in the RG evolution of M 2 /µ 2 .From various numerical investigations we realised that a minima does exist in the evolution of induced M 2 /µ 2 .We choose this minima to be our reference point and choose the initial conditions at this point for all other couplings.The appearance of a minima in the flow induced M 2 /µ 2 implies that at this minima the beta-function of induced M 2 /µ 2 will vanish, Plugging the RG-flows of various coupling in this, will result in a condition satisfied by all the couplings of theory at this minima.This will act as a constraint in choosing some of the initial parameters.We first choose the value of M 2 /µ 2 at this point, we call it ρ * .We require ρ * > 1.
where f 2 * , ξ * and K * are initial values of f 2 , ξ and K respectively.The imposition of this constraint makes sure that the mass of the spin-2 ghost mode remains above the running energy scale.This will imply that one of the three unknowns f 2 * , ξ * and K * can be expressed in terms of other two.We choose to write K * in terms of f 2 * and ξ * .At this point we also choose f 2 * ≪ 1.Now the left unknowns are λ * , (y t ) * , ω * and ξ * .In order to choose the matter couplings λ * we use our knowledge of non-gravitational system.In such system the running λ always hits the Landau pole if the initial value of yukawa coupling is below certain threshold, beyond which λ becomes negative leading to instability.We accordingly choose λ * 0.1.
At this point we analyse the beta-function for the coupling ξ.In this theory we have the freedom to choose ξ to be very large ( 10).This is primarily because in the perturbation theory the coupling strength of vertex containing n-gravitons and two scalars is ∼ ξ( f 2 ) n and ∼ ξ( f 2 /ω) n .Since f 2 ≪ 1, so this give us freedom to choose ξ to be very large while still being in the realm of perturbation theory4 .For ξ large the beta-functions of various coupling acquires a simplified form.Although the beta functions become simplified but still they are complicated enough that it requires the analysis to be done numerically.We tend to explore numerically this regime of parameters.
We choose to work in regime where −∆/ϕ 6 = ǫ ≪ 1, where ∆ is the discriminant mentioned in eq.(63).In this regime there will one positive root for − and a complex conjugate pair with negative real part.Under the RG flow ∆/ϕ 6 will also run.We choose the initial parameters in such a way so that at the initial point ǫ ≪ 1.By reversing this argument we say that we start with ǫ ≪ 1 and solve for the initial parameters under this constraint.This fixes the initial value problem completely.With the chosen f 2 * ≪ 1 ( 10 −6 ), ρ * > 1, λ * 0.1 and ξ * ≫ 1 ( 102 ), we use the constraint dictated by ǫ ≪ 1 to solve for ω * .From the four different solution for ω * , we choose the one which real and positive (to avoid tachyons) 5 .Knowledge of ξ * gives the initial value of K * by using the relation given in eq. ( 82).We then plug these into the minima constraint given in eq. ( 81).This constraint contains the yukawa coupling in quadratic form, and therefore on solving gives two equal and opposite values for y t * .One can choose either of the sign of yukawa coupling for the initial condition.The flow of all the other couplings doesn't depend on this sign.Once the initial parameters are known we can solve the RG flows and compute the flow of induced M 2 /µ 2 to see if it remains above unity throughout the RG evolution.

VI. NUMERICAL ANALYSIS
We tried several possible values of various parameters in order to see how the flows are for various initial conditions, and did the analysis case by case systematically.

A. Fixed λ *
We first considered case with fixed λ * , while we took different values for f 2 * , and for each f 2 * we explored a range of ξ * .Throughout this we took ρ * = 1.5 (there is nothing special about this number, as long as long as ρ * > 1).We considered three different values for f 2 * = 10 −6 , 10 −7 and 10 −8 .We have freedom over the choice of −∆ * /ϕ 6 = ǫ.It is seen that with f 2 * fixed, when ǫ is made smaller, then ω * increases.However the y t * obtained first decreases to a minima before rising again and becoming stable.We choose ǫ near this minima, so that we have more number of e-folds in the RG flows.It turns out that for each value of f 2 * the position of occurrence of this minima will be different.For smaller f 2 * , the minima occurs at a smaller value of ǫ.Thus for f 2 * = 10 −6 , 10 −7 and 10 −8 , the minima for ǫ occurs around 10 −12 , 10 −15 , and 10 −16 respectively.We consider these cases in succession.
The number of e-folds from the Planck's time to current galactic scale is ∼ 130.This stands then as another guiding principle to choose set of initial parameters.It is noticed that when f 2 * is made more smaller then the allowed upper value of ξ * (which is chosen so that we have 100 e-folds) increases.This can be understood by considering the strength of vertices.For vertex containing one graviton leg and two scalar leg, the interaction strength ∼ ξ f 2 .Demanding perturbation theory to remain valid implies ξ * f 2 * 1, which explains the behaviour.We keep ξ * large so that there is sufficient communication between the matter and gravity sector.In table I we tabulate our findings for f 2 * = 10 −6 , 10 −7 and 10 −8 .We then plot the flows of various parameters for the choice of initial parameters written in table I.Each flow is interesting to analyse.The flow of the coupling ξ for various choices of the initial conditions is shown in figure 5.The plot shown in figure 5    only, for other values we observe similar qualitative behaviour which will not be shown here.
From the running of ξ shown in figure 5 we notice that the parameter ξ runs to smaller values in the UV regime.This might be an indication of the possible existence of an UV stable fixed point, but it is hard to give a robust answer in this paper.This however can be justified by looking at the beta-function of ξ given in eq (53) whose r.h.s.can be seen to vanish for a certain choice of coupling parameters.The flow of matter couplings λ and y t for various initial conditions is shown in left and right of figure 6 respectively.For smaller values of ξ * the flow of these couplings remain almost same, while deviations are seen for large ξ * .This is again plotted for f 2 * = 10 −8 , while for other values of f 2 * qualitatively similar behaviour is seen.In the UV the flow of λ is seen to bend and run toward smaller values, which is caused by the yukawa coupling.
The flow of the VeV induces a flow in the Newton's constant.The flow of the induced Newton's constant for various initial conditions is shown in figure 7. The induced Newton's constant goes to zero in the UV and in IR.In UV it is seen to go to zero at a finite energy scale.This is similar to the flow of Newton's constant observed in [5][6][7][8], where the original action was not scale-invariant and contained Einstein-Hilbert piece in the higher-derivative action.This is somewhat interesting to note.Again this is just a numerical observation and not a rigorous analytic argument.By varying the value of f 2 * we notice that the qualitative features of the graph remains same.
Figure 8 shows the flow of parameters M 2 /µ 2 and M 2 /ωµ 2 in left and right respectively.The flow of M 2 /µ 2 is such that it is always above unity (M 2 /µ 2 > 1).This means that the propagator of metric fluctuation after symmetry breaking doesn't witness the ghost pole, as the problematic ghost mode is never realised and never goes on-shell.We observe this to happen for a large domain of coupling parameter space.A similar running of the parameter M 2 /µ 2 was also observed in [5][6][7][8], and was used to establish unitarity criterion for the higher-derivative gravity.The flow of the parameter M 2 /ωµ 2 is different from the one observed in [5][6][7][8], where a monotonic behaviour was seen.In the present case we see a convex structure with a single minima in the flow.If we demand that M 2 /µ 2 | * > 1, then it doesn't imply that M 2 /ωµ 2 | * > 1. However the reverse is always true i.e.
By choosing ρ * to be large enough one can make the scalar mode also physically unrealisable.
We then consider the induced masses in the matter sector and consider the flow of combinations m 2 s /µ 2 and m 2 f /µ 2 , where m s and m f is given in eq. ( 67).The induced RG running of them is shown in figure 9.The running of these is interesting in the sense that both of them has a minima.The value at the minima depends crucially on the initial parameters chosen to make the flow unitary.If we choose ρ * large enough then it is possible that flow of m 2 s /µ 2 and m 2 f /µ 2 will be such that the scalar and fermion will never be realised during the whole RG flow, and they never go on-shell.In that sense they affect the theory indirectly and only gravitationally but they never go on-shell.
In the previous subsection the case with fixed λ * was investigated.It is worth checking the robustness of the qualitative features when other parameters are varied.One particular important parameter is the λ * .It is important to see how the situation changes when λ * is increased.For this we fix the value of f 2 * = 10 −8 , ξ * = 10 5 ρ * = 1.5.We took a range of value of λ * = 0.2, 0.25, 0.3, 0.35 and 0.45.Although the qualitative features are same but there are minor differences.In each case we always witness that the running M 2 /µ 2 has a minima, and the flow is always above unity.This further establishes that there always exist a minima in the flow of induced M 2 /µ 2 , and it also implies that by choosing right set of initial condition it is possible to make the massive tensor ghost innocuous.
As the system contain a mixture of several coupling which are all evolving in different manner, therefore the dynamics of system is rich and interesting.This becomes more apparent when we plot the running of various parameters.The flow of λ, y t and ξ is shown in figure 10.These flows are very much similar to the ones shown for fixed λ * in the previous sub-section.It is seen that as λ * is increased the flow of ξ decrease more sharply in the UV, and in IR the flow goes to higher values, even though starting point is same.The running of yukawa coupling is simple, in the sense that when λ * increases, so does y t * and accordingly the whole RG trajectory for yukawa coupling.The flow of λ is interesting.For higher λ * , ξ flows to higher values in the IR.This makes the λ to run faster toward zero in the IR.In the UV the RG trajectories for λ has self-similarity.
The flow of G, M 2 /µ 2 and M 2 /ωµ 2 is shown in left, centre and right respectively in figure 11.The qualitative behaviour is the same in the sense that the induced Newton's constant goes to zero in the UV at a finite energy scale.It goes to zero in the IR.The RG flows for M 2 /µ 2 and M 2 /ωµ 2 have same qualitative features, and tensor ghost is physically unrealisable even when we increase λ * .Choosing appropriate ρ * will make the Riccion also physically unrealisable.The RG flow of Riccion mass is different from the one seen in [5][6][7][8].The flow of induced masses in the matter sector has similar qualitative features and is shown in figure 12.

C. Fixing Planck's scale
The renormalisation group invariance insures that the flow of couplings doesn't depend on the reference point µ 0 .This gave us freedom to choose the reference point without any conditions.As a result for sake of convenience we choose it to be the point where the flow of M 2 /µ 2 has a minima, where the initial conditions for the flow are imposed.However an interesting thing to ask is to how to relate it with phenomenology?In the sense how does the running of various parameter look like when compared to Planck's scale M P l , whose value is around 1.22 × 10 19 GeV?This is interesting point to ponder on.For this we study the running of induced G, which runs strongly in the UV and goes to zero.From observations of astrophysics and cosmology we know G N ewton remains constant for a large energy range.However it is usually expected that it will undergo strong running near the Planck's scale.For this reason we choose M P l in the regime where induced G witnesses a strong running i.e. near the point where induced G goes to zero.Once this is fixed one can plot the flow of various coupling parameters and induced masses.These are presented in figure 13 and 14.

VII. CONCLUSION
Here in this paper the idea of gravity being induced from scale-invariant theory is considered.The fundamental theory is a coupled system of scale-invariant matter and higherderivative gravity.The lorentzian path-integral of this fundamental theory incorporates quantum fluctuations from both matter and gravity.The matter sector is taken to be simple (a scalar and a dirac fermion).
The effective action of the theory is computed in the 4 − ǫ dimensional regularisation procedure.The divergent part of which gives the RG flow of the various coupling parameters of the theory.These have been computed in past also.We did it again in order to verify the past results.We agree fully with the past literature [26,27].We then compute the one-loop RG improved effective potential of the scalar field on the flat space-time.This gets contribution from the both the gravitational and matter degrees of freedom.This effective potential however contains an instability which comes up as an imaginary piece in the effective potential.The straightforward interpretation of the appearance of an imaginary piece in the effective potential is an indication that the background (flat space-time and constant scalar background) is not stable, and will decay.
The reason for the occurrence of this instability is probed.It is found that this is purely gravitational in nature, in the sense that it arose from occurrence of tachyonic modes in the spin-2 gravitational sector and the gravitational scalar sector.These kind of tachyonic instabilities have been investigated in past [87,88], and is a characteristic feature of gravitational theories coupled with scalar in flat space-time.At finite temperature this instability (also known as Jeans instability) results in collapse of gas of gravitons.While occurrence of this instability is a disturbing feature of the theory and is unavoidable, it is however an IR problem and has no effect on the UV physics.In this paper we considered a different feature of theory.We investigate the issue of ghost appearing due to the presence of higher-derivative terms in the theory, which affect also the UV physics.This is done by investigating only the real part of the effective potential and ignoring the instability caused by the tachyonic modes.
The real part of effective potential develops a VeV.This breaks the scale symmetry and induces mass scale, which in turn generates Newton's constant and masses for matter fields.The propagator of the metric fluctuation field now has mass, and in the broken phase it is easy to see problematic massive tensor ghost of the theory which remains shrouded in the symmetric phase.The scalar mode (Riccion) of the metric fluctuation also picks a mass in the broken phase.The VeV has a running, which in turn induces a running in the various parameters that are generated in the broken phase.It therefore seems sensible to ask question on behaviour of massive tensor ghost under this running in the broken phase, which is the most important aim of the paper.The induced running in the various parameters generated in the broken phase allows to investigate the running of M 2 /µ 2 (where M is the induced mass of tensor ghost).The crucial task of the paper was to see whether there exist a domain of coupling parameters space where it is possible to make M 2 /µ 2 > 1 throughout the whole RG trajectory.Satisfactory arrangement of M 2 /µ 2 > 1 will imply that the massive tensor ghost is never physically realised and never goes on-shell.This issue is however studied numerically, as the complexity of the beta-functions and the complicated running of the induced parameters in broken phase hinders to make analytic progress.
The last part of paper is devoted to numerically studying this issue.A prescription to choose the set of initial conditions so that M 2 /µ 2 > 1 for whole RG trajectory is stated.This involves solving certain constraints.It is realised that the flow of M 2 /µ 2 has a unique minima at a certain point along the RG trajectory.The existence of such a point was analytically proved in the context of higher-derivative gravity including Einstein-Hilbert term [5,6], as the RG equation were simpler.In the present paper however it is not possible to achieve this analytically and numerical support was taken to get evidence for the existence of such a minima.We do see that for a large domain of parameter space such a unique minima does exist, and by requiring that M 2 /µ 2 > 1 at this minima implies that M 2 /µ 2 > 1 for whole RG trajectory.This although is not a robust analytic proof but stands as a strong evidence.We considered different set of initial conditions by varying various parameters in a systematic fashion.In each case it was see that the minima in the flow of M 2 /µ 2 always exits and unitarity criterion can be met.
In this domain of coupling parameter space where M 2 /µ 2 > 1 for the whole RG trajectory, the behaviour of other parameters are studied.The first important thing that is noticed is the existence of a finite UV cutoff in the theory, which was also noticed in [5][6][7][8].Even though we do the analysis in dimensional regularisation, still a cutoff emerges dynamically from the RG equations.Beyond this point the flow cannot be continued and knowledge of higher-loop contributions are needed.In [5][6][7][8] we showed analytically that at this point the coupling ω diverges.In the present context we noticed this numerically.The behaviour of matter coupling λ has an interesting flow.It is seen to increase monotonically, but in the UV this stops and starts to decrease.This is due to yukawa coupling.The flow of yukawa coupling y t increases monotonically, and stops when the cutoff is reached.In the case of Rϕ 2 coupling, the coupling starts to run only near the UV, where it is seen to go toward smaller values, hinting that there might perhaps exists a stable fixed point.
The flow of the induced Newton's constant G is interesting.It approaches zero both in UV and IR.In UV it vanishes at finite energy scale.This was something which was also observed in [5][6][7][8], where Einstein-Hilbert was present in the bare action of the theory and was not induced.In that respect this is surprising that in the present picture of Einstein-Hilbert term being induced from scale-invariant theory, the flow of the induced gravitational coupling is qualitatively similar to the case where EH term is present in the theory to begin with.Such vanishing of induced G is although unexpected but a welcome feature.This is opposite to the widely known feature in Einstein-Hilbert gravity (without higher-derivatives) where Newton's constant becomes very large in UV.However those results cannot be trusted as they appear in non-renormalizable theories.In the presence of higher-derivatives the situation changes in UV.Such vanishing of Newton's constant means that in UV, gravity decouples from matter, although gravitational self-interactions continue to exists.Such a behaviour will have consequences in the early universe, and also justifies the use of flat background in the UV.This softening can also be used in addressing the Higgs naturalness problem [27].
The flow of combination M 2 /ωµ 2 is however a bit different than what has been witnessed in [5][6][7][8], in the sense that the function M 2 /ωµ 2 is no longer a monotonically decreasing function of RG time t.On the contrary it has a flow similar to M 2 /µ 2 , having a single minima.But there is a region of RG time t over which M 2 /ωµ 2 < 1.This is because of the choice of initial condition M 2 * /µ 2 * = ρ * we made.This will imply that there is a range of t where this scalar mode will be realised and will go on-shell, and outside this region it will remain unrealised.This can play a role in early universe to drive inflation.On the other hand the parameter ρ * can be chosen appropriately large in order to make this scalar mode ghost (never goes on-shell), while still having unitarity.
The running of VeV also induces a running in the generated masses for the scalar and fermion.To analyse whether they are physically realised or not, we studied the behaviour of m 2 s /µ 2 and m 2 f /µ 2 respectively.It is seen that if we choose ρ * appropriately, then it is possible to make them not physically realisable.They never go on-shell but do effect the theory gravitationally.
The RG flow equations for the dimensionless couplings are gauge independent at oneloop however at higher loop gauge dependence is expected to enter.The effective potential is gauge-dependent which is because the hessian carries gauge dependence.This gauge dependence then enters the VeV and any quantity which is related to VeV (induced Newton's constant and induced masses).In the paper we studied the problem in Landau gauge which is a physical gauge allowing propagation of only transverse modes and suppressing longitudinal ones.However such gauge dependence is expected not to change qualitative features.This was explicitly seen in the case of pure higher-derivative gravity without matter [5,6].
In the action appearing in the integrand for the lorentzian path-integral, the sign of coefficient of C 2 µνρσ is taken to be negative while the sign of coefficient of R 2 is taken positive.This is done to avoid tachyons and make the poles (in broken phase) lie on real axis, the inspiration of which comes from past study done in [5][6][7][8].Such a choice further offers necessary convergence in the lorentzian path-integral in the feynman +iǫ-prescription.This sign choice then implies that the coupling f 2 is no longer asymptotically free different from what is seen in path-integral defining a positive-definite euclidean theory [9,10,16,[30][31][32][33] (which is an entirely different theory), but instead has a Landau pole.This Landau singularity however appears way beyond the point where the RG flow of all couplings stops.Also the occurrence of Landau singularity is very possibly a one-loop effect, as at higher loops the running of f 2 gets correction thereby hinting the occurrence of fixed point [5,6].Moreover in this theory the dimensionless perturbative parameters f 2 , f 2 /ω, ξ f 2 , ξ f 2 /ω, λ and y t remain small throughout the RG flow, thereby justifying usage of perturbative approximation and we don't enter non-perturbative regimes in our analysis.The analysis done in the paper is on a flat background.This is because any generic background has locally flat regions allowed by (strong) equivalence principle.Also when one is probing ultra-short distances, one can do the analysis on flat space-time, as the perturbative UV divergences are independent of the background.Moreover, the chances of tensor-ghost becoming physically realisable is more in UV (and nowhere else), therefore its avoidance is important in UV, which can be investigated on a flat space-time.However extrapolating flat space-time analysis in deep infrared can lead to erroneous conclusions.In performing this study cosmological constant was put to zero, and was argued that it can be maintained to be zero in a supersymmetric framework.However supersymmetry is broken in IR and this will generate cosmological constant back again.Moreover current observations also favours the existence of cosmological constant in order to drive the accelerated expansion of the universe.Therefore a proper treatment of IR physics in a field theoretic language is needed.A possible direction would be to formulate the theory on a deSitter background [92][93][94][95] (see also [96] and references therein).This will give more accurate description of the theory in the infrared.
The existence of tachyonic instability is a further indication that the chosen background of flat space-time is unstable.While this is an IR effect and an unavoidable outcome of gravitational theories coupled with scalar on flat space-time, it signals the breakdown of flat space-time as the background.It is expected that performing the study on a curved background might address these issues.For this, one would require a more accurate description of the formulation of field theory on curved background, and obtain the results in low energy limit.Alternatively one may have to incorporate non-localities appropriately to deal with IR physics [97,98].This is a future direction and will be considered later.
It is interesting to wonder whether the RG equations gets modified when the decoupling of massive spin-2 ghost mode occurs, in the same manner as in flat space-time non gravitational QFTs where decoupling theorem exists [99].In flat space-time QFTs a systematic order-by-order computation leads to decoupling of heavier modes in process occurring at energies less than the mass of heavy particle.This theorem has been suitably extended to the case of matter theories on curved background [83][84][85], where the beta-functions gets a correction in mass-dependent scheme.For spin-2 fields the situation is a bit more involved, as incorporating mass in a deffeomorphism invariant manner is a tricky task.A possible way to achieve is by including higher-derivative terms in the action, which immediately brings in ghosts.If the ghost mass however is always above the energy scale, then ghosts get avoided due to decoupling.But this occurs in the quantum theory where RG running of the ghost mass is always above the energy scale.This implies an effective decoupling in the sense this spin-2 ghost mode never goes on-shell and off-shell it doesnt contribute to imaginary part of amplitudes [5,6].But currently it is unclear how such a decoupling will modify the RG flow equations of various parameters.This in really worthy of investigations and will be considered in future works.

Propagators for Matter fields
Here we write the inverse propagators for the various matter fields.These are obtained by doing the second variation of the action of the theory with respect to various fields.The mixed terms in such kind of variation will be treated as interaction terms.From the second variation of the action given in eq. ( 19) one can pick the terms corresponding to the scalar and fermion propagator.The operator whose inverse correspond to scalar propagator is, (C1) In the case of fermions the relevant inverse operator is given by, (∆ F ) ab = iγ ρ ab ∂ ρ . (C2)

Vertices
Here we specify the various vertices that are relevant for our one-loop computations.These can be categorised in 3 parts: (a) vertex with two internal graviton lines, (b) vertex with one internal graviton line and one internal matter line, and (c) vertex with two internal matter lines.

a. Gravity-gravity
In these vertices there are two internal graviton lines.In the following the term V µνρσ comes from scalar field action, while the term U µνρσ comes from fermion field action.The vertices are depicted in figure 15.Here we write the vertex that contain one internal graviton line and one internal fermion line.These vertex comes only from the fermion field action.These vertices are depicted in figure 17.Here we write the vertices which has two internal matter lines.These will be either both scalar lines, one scalar and one fermion line or both fermion lines.These vertices are depicted in figure 18.  Various vertices containing two internal matter legs (scalar-scalar, scalar-fermion, fermion-fermion).Here the dashed line is scalar while solid line with arrow is fermion.The lines ending with circle containing cross are external legs.potential from the scalar sector of theory.In order to compute the contribution of the scalar to the effective potential we need analyse a cubic equation in − operator written in eq.(62).Here in this section we will consider a generic cubic equation of the form where ∆ was also defined in eq. ( 63).The two roots of the quadratic equation will be given by, This can be solved easily by taking cube-root.Here there will be three roots for v. Corresponding to each cube-root we have a root for the eq.(D1), which is obtained by plugging v back in to u and z.The three roots will be given in terms of l 1 and l 2 (D7) In the case when the discriminant ∆ 2 1 − 4∆ 3 0 > 0, we have one real root and a complex conjugate pair.In this case we can write the exponentials in terms of sine and cosine functions.One can then write the real and imaginary part of these roots.This complex conjugate pair can be written in polar form also with a) − Tr(b) .

FIG. 1 :
FIG.1: Various diagrams containing divergences at the tadpole level.Here the dashed line represent scalar field, solid line with arrow represent fermion field while double-line represent h µνfield.The first two graphs are purely matter ones, while the other four graphs contain quantum gravity corrections.

FIG. 2 :
FIG. 2: Various diagrams containing divergences with two vertices.Here the dashed line represent scalar field, solid line with arrow represent fermion field while double-line represent h µν -field.Here the graphs on the first line are purely matter ones.The second and third line contain graphs having quantum gravity corrections.

FIG. 3 :
FIG.3: Various diagrams containing divergences with three vertices.Here the first graph is purely matter oriented and gives correction to yukawa coupling.The next two diagrams are giving correction to ϕ 4 coupling.They are only present in the quantum gravity context.The last two diagrams are giving quantum gravity correction to the yukawa coupling.

FIG. 17 : 2 η
FIG.17: Various vertices containing one internal graviton line and one internal fermion leg.Here the dashed line is scalar, solid line with arrow is fermion, while double line depicts graviton.The lines ending with circle containing cross are external legs.
) (V ψψ ) ab = −y t ϕδ ab , (V T ψψ ) ab = y t ϕδ ab , (C12) (V ψφ ) a = −y t θ a , (V T ψφ ) a = y t θ T a (C13) (V φψ ) b = −y t θb , (V T φψ ) b = y t θTb .(C14)Appendix D: Cubic EquationHere we will consider the roots of generic cubic equation with real coefficients.Such an equation emerges in section.III while computing the contribution to the scalar effective