Renormalization in a Lorentz-violating model and higher-order operators

The renormalization in a Lorentz-breaking scalar-spinor higher-derivative model involving $\phi^4$ self-interaction and the Yukawa-like coupling is studied. We explicitly de- monstrate that the convergence is improved in comparison with the usual scalar-spinor model, so, the theory is super-renormalizable, with no divergences beyond four loops. We compute the one-loop corrections to the propagators for the scalar and fermionic fields and show that in the presence of higher-order Lorentz invariance violation, the poles that dominate the physical theory, are driven away from the standard on-shell pole mass due to radiatively induced lower dimensional operators. The new operators change the standard gamma-matrix structure of the two-point functions, introduce large Lorentz-breaking corrections and lead to modifications in the renormalization conditions of the theory. We found the physical pole mass in each sector of our model.


I. INTRODUCTION
It is well known that the Lorentz-breaking field theory models can be introduced in several ways. We can list some of the most popular approaches. First, one can introduce small Lorentz-breaking modifications of the known theories through additive terms, thus implementing the Lorentz-breaking extensions of the standard model [1]. In principle, the most known extensions of the QED follow this way. A very extensive list of the possible Lorentz-breaking additive terms in different field theory models including QED is given by [2]. Second, one can start with the modified dispersion relations [3], and, in principle, try to find a theory yielding such relations. Third, the Lorentz-breaking theories can be treated as a low-energy limit of some fundamental theories, for example string theory [4] and loop quantum gravity [5]. Finally, the Lorentz symmetry can be broken spontaneously, see f.e. [6]. The main motivation behind all these approaches, however, is the same, and resides in the expectation that any experimental evidence of departure from Lorentz symmetry may provide the first germs towards the construction of a theory amalgamating both General Relativity and the Standard Model of particle physics.
At the same time, it is natural to consider one more aspect of studying the Lorentz-breaking extensions of the field theory models. It consists in introducing the essentially Lorentz-breaking terms, that is, those ones proportional to some constant vectors or tensors, involving higher derivatives. As a result, the corresponding theory will yield an essentially different quantum dynamics. The first known example of such a theory is the Myers-Pospelov extension of the electrodynamics [7] where the three-derivative term essentially involves the Lorentz symmetry breaking. Another important example of such a theory is a four-dimensional Chern-Simons modified gravity with the Chern-Simons coefficient chosen in special form Θ(x) = k µ x µ [8], which, in the weak field limit, also involves third order in derivatives of the dynamical field (that is, the metric fluctuation). Moreover, the importance of the Myers-Pospelov-like term, and analogous terms for scalar and spinor fields which can be easily introduced, is also motivated by the fact that a special choice of the Lorentz-breaking vector will allow to eliminate the presence of higher time derivatives thus avoiding the arising of the ghosts which are typically present in theories with higher time derivatives (see f.e. [9]). Also, this term was shown to arise as a quantum correction in different Lorentz-breaking extensions of QED [10] and has been studied for causality and stability [11]. In the case of including higher time derivatives, it has been shown recently that the unitarity of the S-matrix can be preserved at the one-loop order in a Myers-Pospelov QED [12]. The proof has been accomplished using the Lee-Wick prescription for quantum field theories with negative metric [13]. For other studies on unitarity at tree level for minimal and nonminimal Lorentz violations, see [14,15] respectively. It is important to notice that the Myers-Pospelov-like modifications of QED are actually experimentally studied as well within different contexts [16].
We emphasize that, up to now, the quantum impact of the Myers-Pospelov-like class of terms being introduced already at the classical level, where the higher-derivative additive term should carry a small parameter which can enforce large quantum corrections [17], almost was not studied except of the QED [18] and superfield case [19]. The presence of such effect raises the question how to define correctly the physical parameters in the renormalized theory. On the other hand, for studies in the context of semiclassical quantization it is natural to consider the presence of higher derivative terms in order to implement a consistent renormalization program [20]. With these considerations, the natural question is -what are the possible consequences of including the Lorentz-breaking higher-derivative terms into the classical action?
It is well known that loop corrections in Lorentz-invariance violating quantum field theory may lead to new kinetic operators absent in the original Lagrangian. Recently, the consequences of these radiatively induced operators have been studied in relation with the finiteness of the S-matrix and the identification of the asymptotic state space [21]. These new terms introduce modifications in the propagation of free particles and change drastically the physical content of the space of in and out states. In particular, the Kallen-Lehmann representation [22] and the LSZ reduction formalism [23] are modified in the presence of Lorentz symmetry violation [24]. An important finding is that spectral densities which in the standard case are functions of momentum-dependent observer scalars such as p 2 , in the Lorentz violating scenario may depend on other scalars such as couplings of Lorentz-violating tensor coefficients with momenta [24]. This has led to modifications in the renormalization procedure, in the definition of the asymptotic Hilbert space and in general in the treatment for external-leg physics [21]; for other studies of the renormalization in Lorentz-breaking theories, see also [25]. A natural extension for these studies is to consider the nonminimal framework of Lorentz invariance violation, that is, when the Lorentz-breaking is performed with higher-order operators [26]. It is well known that the inclusion of higher-order operators in quantum field theory will generate, via radiative corrections, all the lower dimensional operators allowed by the symmetries of the Lagrangian. For the case of breaking the Lorentz symmetry, let us say in QED and with a preferred four-vector n µ , the induced operators may involve contractions of n µ with matrices other than just γ µ , together with scalars such as (n · p). The new terms force to modify the renormalization conditions in order to extract the correct pole mass from the two-point functions. In particular, the renormalization condition for the renormalized fermion self-energy Σ R ( / p = m P ) = 0, with m P being the physical pole mass, has to be generalized, which ultimately will depend on the form of Lorentz breaking. In this work we continue these studies in order to carry out the renormalization in a theory with higher-order operators and in addition we study the possible effects of large Lorentz-violating corrections. Within our study, we consider the renormalization of the higher-derivative Lorentz-breaking generalizations of λφ 4 and Yukawa model.
The structure of the paper looks like follows. In Sec. II, we consider the classical actions of our models, write down the dispersion relations, find the poles and describe their analytical behavior in complex p 0 -plane. In Sec. III, we compute the quantum corrections corresponding to the self interaction λφ 4 . In Sec. IV we discuss the coupling of scalar and spinor fields and provide a study of the degree of divergencies in our model. In Sec. V, we compute the two-point functions in purely scalar and scalar-spinor sectors thus exhausting possible divergences and showing explicitly the radiatively induced operators with new gamma-matrix structure and large Lorentz-violating terms. In Sec. VI, we perform the mass renormalization in both sectors and find the physical masses in the theory. In the last section, we discuss our results, and in the Appendices A, B, we provide some details of the calculations.

II. THE EFFECTIVE MODELS AND POLE STRUCTURE
We are interested in the higher-order Lagrangian density describing two sectors of Lorentz-breaking theory: The first sector involves a scalar sector with a fourth derivative together with a self-interaction potential term and the second one the fermionic Myers-Pospelov model [7], with dimension-five operators and the Yukawa coupling vertex The constants g 1 = κ M 2 P l and g 2 = η M P l parametrize the higher-order Lorentz invariance violation with M P the Planck mass representing itself as a natural mass scale, κ, η are dimensionless parameters, whose presence describes the intensity of the higher-derivative terms, and n µ is a dimensionless four-vector defining a preferred reference frame. The propagators in momentum space read We begin an analysis of the dispersion relations in both sectors. A further motivation for its study, and consequently, the finding of the poles and their analytical behavior in complex p 0 -plane, consists first in the fact that in our models namely using of the residues of the propagators is a most convenient approach for calculating the quantum corrections. Second, in the presence of higher-order time-derivative terms a direct implementation of the i prescription may lead to a wrong four-momentum representation for the propagator which may spoil any attempt to preserve unitarity or causality. Let us start with the scalar dispersion relation, namely which for a purely time-like four-vector n µ = (1, 0, 0, 0), the solutions are given by and where E( p) = p 2 + M 2 . The dispersion relation can also be written as (p 2 0 − p 2 1 )(p 2 0 + P 2 2 ) = 0, hence one has the solutions p 0 = ±p 1 and p 0 = ±iP 2 so that Their exact location in complex p 0 -plane and also the contour of integration C (s) F are shown in Fig. 1. The solutions can be classified according to their perturbative behavior when taking the Lorentz violation to zero. We identify two standard solutions ±p 1 which are perturbative solutions to the usual ones ±E and two complex ones (and moreover, actually tachyonic) ±iP 2 which diverge as g 1 → 0. The extra solutions that appear ±P 2 are associated to negative-metric states in Hilbert space and have been called Lee-Wick solutions [13].
Alternatively, we can write the scalar propagator as for the fermion propagator which is defined to round the negative pole from below and the positive poles from above. At higher energies than 1/(4g2), the two solutions ω1 and W1 become complex and move in opposite directions along the imaginary line starting at 1/(2g2), we deform the contour continuously avoiding any crossing or singularity with the poles.
which agrees with the usual propagator in the limit g 1 → 0.
In the fermion sector we have the dispersion relation Again for the time-like n µ we have the equation whose standard, that is, non-singular at g 2 → 0, solutions are and the Lee-Wick ones where E( p) = p 2 + m 2 .
In the region of energies satisfying the condition 4g 2 E < 1 the four solutions are real and obey the inequality where ω 2 is a negative number. However, beyond the critical energy 1/(4g 2 ) both ω 1 and W 1 become complex and move in the opposite imaginary line at 1/(2g 2 ) as shown in Fig. 2, while the other two solutions ω 2 , W 2 remain real.
To define the contour C (f ) F we use an heuristic argument, specially to go beyond the critical energy at which complex solutions appear. We implement a correct low energy limit by considering the prescription given in [27] which has been well tested to give a suitable correspondence with the normal theory when g 2 → 0 and also to preserve the unitarity of the S matrix. In this effective region the integration contour C (f ) F is defined to round the negative pole from below and the three positive ones from above. Now we increase the energy to values at which the two solutions ω 1 and W 1 become complex, and define the new contour as the one obtained by continuously deforming the curve by avoiding any crossing and singularity with the poles, as shown in Fig. 2.
With this consideration in mind, the fermion propagator reads which differs from the direct i prescription in the quadratic terms, but allows in particular to define a consistent Wick rotation which we use later.

III. THE INTERACTION λφ 4
In this section we explore the potentially divergent one-loop radiative correction in the scalar propagator which is generated by the well-known tadpole graph given by Fig. 3. To proceed with it, we need to evaluate the basic integral Note that a naive power counting gives a logarithmic divergence for the integral which, however, as shown below using dimensional regularization is found to have a finite result in four dimensions; in a similar fashion of what happens with the Riemann zeta function ζ(s) for negative values of s. We go to d dimensions and choose the Lorentz-breaking four-vector to be timelike n µ = (1, 0, 0, 0) which yields where p 1 , P 2 are given in (7). We perform the integration in the complex p 0 -plane by closing the contour C (s) F upward and enclosing the two poles −p 1 and iP 2 as depicted in Fig. 1, yielding where Note that iF 1 − F 2 has the correct limit at g 1 → 0, recovering the usual result − 1 2E . Now, it is convenient to change variables z = 1 + 8g 1 E 2 (p) yielding which allows to write the integral (16) as with where z 0 = 1 + 8g 1 M 2 and we have used the definition of solid angle (A2).
Considering both contributions through the relation Σ 2 = Σ (1) + Σ (2) , and after some algebra with d = 4 − and expanding in , we find at lowest order Here, 2 F 1 R (0,0,1,0) 1 4 , 3 4 , 2, 1 is a hypergeometric function, whose exact value is −0, 71. Note that there is a fine tuning in this case, that is, the expression is singular at g 1 → 0. However the correction to the two-point function Σ 2 is UV finite.

IV. COUPLING OF SCALAR AND SPINOR FIELDS
Let us consider the theory involving both the quartic interaction vertex V 1 = − λ 4! φ 4 and the Yukawa coupling vertex V 2 = gψψφ. We note, that, in principle, the second time derivatives in a free action of a spinor field are present also in specific Lorentz invariant theories, for example, the known ELKO model [28]. However, our theory essentially differs from that model. To classify the possible divergences, we should calculate the superficial degree of divergence ω of this theory. The naive result for it is where E ψ is a number of spinor legs. However, this manner yields incorrect results because of the strong anisotropy between time and space components of the momenta (for example, in this case one can naively suggest that the two-point function of the spinor field can yield only the renormalization of the mass of the spinor field). So, let us proceed in the manner similar to that one used for Horava-Lifshitz-like theories (cf. [29]). Since n µ is purely time-like, we can write (n · p) 4 = p 4 0 , so, we have from (4) Following the methodology developed for the Horava-Lifshitz theories (see f.e. [29]), we suggest that the denominators of the propagators are the homogeneous functions with respect to higher orders in corresponding momenta, and the canonical dimension of the spatial momentum p is 1. Taking into account only the leading degrees, we easily conclude that the canonical dimension of the momentum p 0 is 1/2 (we note that this case does not occur in usual Horava-Lifshitz-like theories where the canonical dimension of time momenta are always more than one, cf. [29]). Therefore, the spinor propagator has the canonical dimension (and the contribution to the superficial degree of divergence) equal to (−1), and the scalar one -equal to (−2) just as in the usual case. Nevertheless, the dimension of the integral measure, that is, d 4 k = d 3 kdk 0 in this case is different from the usual one, being equal to 7/2 rather than 4. Hence the superficial degree in our theory is where L is a number of loops, and P φ and P ψ are the numbers of scalar and spinor propagators respectively. Then, let V 1 will be the number of φ 4 vertices, and V 2 -of Yukawa-like vertices. One has the identities for numbers of scalar and spinor fields in an arbitrary Feynman diagram: where E φ , E ψ are the numbers of external scalar and spinor legs respectively. We use the topological identity As a result, we eliminate numbers of loops and propagators from ω and rest with A straightforward verification shows that the superficially divergent diagrams (that is, those ones with ω > 0) can be of the following types: . This is the one-loop renormalization of the mass and kinetic terms for the spinor.
(iv): E ψ = 0, E φ = 2, V 2 = 4, ω = 1. This is the two-loop renormalization of the mass and kinetic terms for the scalar. ( . This is the one-loop renormalization of the mass term for the scalar. Actually, we already showed in the previous section that, due to the specific structure of poles of the propagator, this contribution is finite. Actually, in the cases (iii) and (iv) the divergence will be not linear but logarithmic, by the reasons of symmetry of integrals over momenta. The diagrams with odd numbers of E φ will vanish due to an analogue of the Furry theorem. So, our theory is super-renormalizable. Moreover, we note that since the kinetic term for the scalar involves two derivatives acting to the external fields, its superficial degree of freedom should be decreased at least by 1, if these derivatives are the time ones, and by two for space derivatives; actually, in one-loop case in a purely scalar sector the kinetic term simply does not arise. Also, in the cases (i) and (ii) one will have the only divergent contribution to the mass of the spinor. So, taking into account the previous section as well, we conclude that at the one-loop order one could have only the renormalization of the masses of the spinor and the scalar arisen from the Yukawa-like coupling.
We note that namely this degree of divergence correctly explains why the self-energy of the fermion diverges, as we will see further (indeed, the naive calculation yields a finite result for it). To study the renormalization, we can restrict ourselves by the lower order, that is, one loop.
So, we rest with only three potentially divergent graphs -with V 1 = 1, that is, the purely scalar tadpole we studied above, and with V 2 = 2 and E ψ = 0, E φ = 2 or E ψ = 2, E φ = 0 we study below.

V. THE YUKAWA-LIKE THEORY
In the next subsections we compute the radiative corrections to the scalar and fermion two-point function in the Yukawa-like theory which arises by considering the self-interaction term V 1 → 0 and g 1 → 0 in (1). The Lagrangian is and additionally, we impose the simplification of considering m = M and the preferred four-vector to be purely timelike n = (1, 0, 0, 0).
A. Scalar self-energy Π(p) As a first example of quantum corrections in our Yukawa-like model, we study the contribution with two external scalar legs depicted at It is represented by the integral where we define Calculating the trace gives Let us write the corresponding contribution to the effective action as Π(p) = −2g 2 φ(−p) Π(p)φ(p) and study the typical low-energy behavior of this contribution by expanding it into Taylor series The zeroth-order contribution follows directly from (30) It is convenient to rewrite as where and where the integrals (34) and (35) have been solved in the Appendix A. For the next term, it is clear that ( ∂ Π ∂p µ )| p=0 can be proportional to n µ only since there is no other vectors, the corresponding contribution to the effective action will yield d 4 xφ(n·∂)φ, that is, a surface term. So, we can disregard it. Further, one would need then to find the second derivative, that is, ( ∂ 2 Π ∂pµ∂pν )| p=0 which may contain naturally terms of higher-order in g 2 . To find it, consider Integrating by parts and neglecting the surface terms, we obtain the symmetric expression where we have used the identity ∂f (k+p) ∂pα p=0 and after some algebra we arrive at with By using the relations one obtains Considering the tensors available in our model, which are the flat metric η µν and the preferred four-vector n µ we can write and Now, consider the relation and multiplying by n µ n ν we to arrive at and by contracting with the metric η µν Solving the algebraic equation we have A similar analysis gives We find the second-order contribution Reorganizing this expression, we can write the correction to the scalar propagator up to second-order in p as where q 0 = K + 2m 2 P , (52) q 1 = −4 P + 2m 2 n 2 M , q n = 8 2g 2 n 2 S − m 2 L .

Finally one has
We provide details of the computation of q 0 , q 1 and q n in the Appendix A. The two-point function is finite and involves a fine-tuning term proportional to g −2 2 . The Lee-Wick modes improve the convergence of the theory such to make the two-point function of the scalar field essentially UV finite and involves the aether term [30].

B. Fermion self-energy Σ(p)
Now we focus on the contribution of the fermion self-energy graph depicted in Fig. 4, and recall we are considering M = m. The fermion self-energy graph is represented by the integral To find it, let us consider a Taylor expansion of the first denominator term up to second-order in p and rewrite the diagram as With the notation Σ(p) = g 2 I (0) + g 2 I (1) + g 2 I (2) , we introduce the the zeroth-order contribution the linear-order contribution and the second-order contribution C. The gamma-matrix structure of I (0) , I (1) I (2)

Zeroth-order I (0)
Let us start with Eq. (58) and rewrite it as where we have defined and From tensor analysis considerations one should have Replacing the expression (64) in Eq. (61) produces the zeroth-order contribution with We carry out the calculations of f 0 and f n 1 following the lines given in Appendix B 2. The first coefficient f 0 is naturally finite and the second one f n 1 is divergent and contains a large Lorentz-breaking correction term of the order of g −1 2 .

Linear-order I (1)
The linear-order integral (59) can be rewritten by introducing where By considering = n µ n ν B + n 2 η µν C , and after some manipulations one finds Introducing the new notation we can write where The terms f 1 , f n 2 and f n 3 are convergent and are explicitly calculated in Appendix B 3.

Second-order I (2)
Following the same methodology the integral I (2) can be written as with , We have that I (2) µν and f 2 are of the order g 2 2 so they can be neglected. For the remaining we define them with the following tensor structure = n µ n ν n αD + n 2 η µν n αĒ + n 2 η να n µF .
We rewrite as with where the explicit form of the relevant constants is given in the Appendix B 4.

VI. MASS RENORMALIZATION
A. The on-shell subtraction scheme Having computed the one-loop correction to the scalar and fermion propagators, we now proceed to determine the pole mass in our Yukawa-like model. The radiative corrections essentially involve extra terms such as the scalar (n · p), the matrix / n and their combinations. The theory is finite in the scalar sector and in the fermion sector the only divergency appears in the term f n 1 . Since the mass corrections are finite in both sectors, in some sense we are working in the scheme in which the finite parts of the renormalized mass correspond to the pole mass, and hence one can say that we are working in the on-shell subtraction scheme. Most of the ideas and method of calculation in the derivation of the pole mass, which are given below, are along the lines developed in the work of Ref. [21].
Let us remind the leading quantum corrections in the scalar sector where we have restored the constant −2g 2 . In the fermion sector we have found The coefficients q and f have been computed in the Appendix A and B.

B. The scalar pole mass
Let us start with the renormalized scalar Lagrangian and without loss of generality we absorb a term proportional to g −2 2 in the definition of renormalized mass m 2 R as appears in the correction of q 0 in (53).
Consider the renormalized two-point function for which Let us write where A φ , B φ and C φ are constants that can be deduced from the expressions (53) , (54), (55) being In order to find the pole defined by the condition atP 2 φ = 0 on the renormalized two-point function we defineP where M ph andȳ are unknown constants we should find. From (87) replacing the value of p 2 given in (92) and using the condition (91), we arrive at the equation Due to the independence of each term, and after some algebra, we find (95)

C. The fermion pole mass
In the fermion sector we consider the renormalized Lagrangian and a corresponding counterterm in order to make the quantum corrections finite and we associate the counterterm with the interaction term The renormalized two-point function is with We can write the finite contribution to the renormalized two-point function up to second-order in p as where the explicit Lorentz violation coefficients at linear order in g 2 , and coded in Σ 2 , are given by We work in the minimal subtraction scheme and hence we fix the divergent term to be (101) Within this scheme and according to the previous calculations we make the identification To find the pole we consider the ansatzP We include the large Lorentz violating terms in the pole extraction process which has been included explicitly inx.
Considering the condition forP to be a pole, that is to say, provides us with the six terms m ph , m n ,x (0) ,x (1) ,x (2) ,x (3) . We proceed as follows: we replaceP in (97), and arrive at At lowest-order we have the first condition which has the solution Considering that all coefficients are independent, we have at linear-order for the factor accompanying (n · p): and consequently we have The next conditions determinex which follows from the equations with factor accompanying 1 g2 / n, (n · p)/ n , g 2 (n · p) 2 / n, g 2 p 2 / n respectively. Therefore, one arrives at This finishes the pole extraction in both sectors.

VII. SUMMARY
We considered the Myers-Pospelov-like higher-derivative extensions of the Yukawa model which incorporates possible new physics from the Planck scale through dimension five operators coupled to a preferred four vector n µ which breaks the Lorentz symmetry. We have selected a particular configuration of Lorentz symmetry violation in which the preferred four-vector n µ is purely timelike. This choice produces higher-order time-derivative terms and leads to extra solutions and new poles in the model. We have found and identified the poles associated to standard solutions and those ones corresponding to negative-metric states or Lee-Wick solutions. Some of these poles move in the real axis, as in the usual Yukawa model, however, above a specific energy called the critical energy some solutions become complex introducing a extra difficulty to define a consistent prescription for the contour of integration in the complex plane [31]. To solve this problem we have analyzed the motion of the poles in the complex plane. In the scalar sector we have defined the usual Feynman integration contour which rounds the negative pole from below and the positive from above; the new poles that appear are located in the imaginary axis and hence they do not present problems in this sense. However, in the fermion sector we have defined the prescription with two considerations in mind: first, the ability to recover the standard location of the poles relative to the real axis when the Lorentz invariance violation is turned off. With this consideration we have imposed the minimal requirement on the contour C (f ) F to round the negative pole ω 2 from below and ω 1 from above. Second, to account for the other poles we have used the prescription putted forward in [27], which has been shown to lead to a unitary S-matrix in the regime of real solutions. Using this Lee-Wick prescription, the three positive solutions W 1 , W 2 and ω 1 are defined to lie below the contour and the negative one ω 2 above. At higher energies beyond the critical energy, we have called for the following heuristic construction. We start at lower energies in which all the poles are real and where the basic requirements for consistency of the theory are satisfied. Next, we begin to increase the energy at which complex solutions appear and define the new contour as the one obtained by continuously deforming the first one in such a way to avoid any crossing and singularity with the complex poles.
A central part of this work has been the study of mass renormalization in the presence of both Lorentz-invariance violation and higher-order operators. We have computed the one-loop corrections in the model and shown how they lead to modifications in the renormalization procedure by pushing the pole mass to a sector involving other gamma matrices besides of / p. The significance of our results, in particular, is based on the fact that this is one of the first works on the renormalization of higher-derivative Lorentz-breaking theories within the context of modifications of asymptotic Hilbert space, whereas even without introducing the higher derivatives the renormalization issues in Lorentz-breaking theory were considered only in a few papers, in particular [25,32,33]. For other approaches in renormalization of higher-order Lorentz breaking theories, see [34]. We found that the theory is super-renormalizable (in principle, increasing the value of N in the added kinetic operator (n · ∂) N , we can get a completely finite theory, the same result can be achieved through a supersymmetric extension of the theory, see [19]). The kinetic terms of the both fields are explicitly finite. As a by-product, we conclude that the aether-like terms originally introduced in [30,35] are generated within our studies both in scalar and spinor sectors, and they are essentially finite. It is interesting to note that in the spinor sector, the aether terms dominate in the low-energy limit in comparison with the second-derivative term we introduced into the classical action.
At the same time, the large quantum corrections or, in other words, fine-tuning arise in our theory. We have found in the scalar sector a large correction which in principle should be added to the usual fine tuning proportional to the square of the fermion mass. In consequence the effect of UV sensitivity of the scalar mass is maintained in our model. In the fermion sector we have that the large corrections affect a term which do not involve any physical observable protecting the theory in this sector to a genuine fine tuning. In a certain sense, this effect of large quantum corrections taking place in our theory resembles the effect of UV/IR mixing taking place in noncommutative field theories responsible for arising infrared singularities in a small noncommutativity limit [36]. Let us mention that in principle these large corrections are natural to expect since the higher-derivative extensions of the classical action can be treated as a kind of the higher-derivative regularization, hence, removing the regularization we return to singular results. We note that this effect is rather generic since the large quantum corrections are present even in the supersymmetric extensions of the Myers-Pospelov-like theories [19]. The complete elimination of large quantum corrections could consist in employing the extended supersymmetry which is well known to achieve complete finiteness of the field theory models as occurs for example for N = 4 super-Yang-Mills theory. Therefore, the natural problem could consist in study of some Myers-Pospelov-like extensions of N = 4 super-Yang-Mills theories. Another relevant problem consists in studying of the two-loop approximation thus exhausting the possible divergences. We plan to carry out these studies in one of our next papers.

Acknowledgments
This work was partially supported by Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq). The work by A. Yu. P. has been supported by the CNPq project No. 303783/2015-0. CMR wants to thank the hospitality of the Universidade Federal da Paraiba (UFPB) and acknowledges support by FONDECYT grant 1140781. A. Yu. P. thanks the group of Física de Altas Energías of UBB for the hospitality. We also want to thank Markos Maniatis for valuable comments and helpful discussions.
Appendix A: The calculation of q0, q1, qn We start to compute q 0 in (52), so we focus on K and P as given in Eqs. (34) and (35) and promote the integrals to d dimensions and consider where we have performed an integration in the angles producing the solid angle in d − 1 dimensions We obtain We rewrite the contour integrals making explicit their poles, and to compute them we use the method of residues with the contour C (f ) F closed from above enclosing the pole ω 2 , see Fig. 2. We find The next integral in the variable E is direct, and gives at lowest order in = 4 − d for both pieces It is simple to show that by combining the two contributions through q 0 = K + 2m 2 P , produces The second-order contributions q 1 and q n are given in terms of P , M , S and L as shown in Eqs. (52). We begin with M , and since it is obtainable in four dimension, we just set d = 4, giving Solving the first integral gives and calculating then the E-integral, we arrive at With the result of P given in (A5), one has We continue to compute q n . From (49) we have In the same way, we arrive at To compute L given in the first equation (48), we note that with We find and therefore L = 0. Considering the leading contributions for q n we obtain Appendix B: The calculation of Σ(p)

The general strategy
In this subsection, we describe the main steps to derive the basic integral that will appear in the computation of the fermion self-energy correction (57). Let us start to consider the integral below where each order is labelled with α = 1, 2, 3 and where F (k 0 , k) is an arbitrary function of k 0 and k. We promote the integrals to d dimensions and recall the phase space measure given in (A1) in order to write Working the denominator, we can rewrite the last contour integral as Now, using the Feynman parametrization and using β = 1 we write the integral as Next, we carry out the change of variables k 0 = k 0 − (1−x) 2g2 and perform a subsequent Wick rotation k 0 → ik 0E , as a result, and dropping the tilde, we arrive at Replacing the above expression in (B2) and again changing variables by the rule t = g 2 E produces the final expression with which is the basic integral we need to solve in the next subsections.
2. The zeroth-order terms f0 and f n 1 We start to compute f 0 in d dimensions which follows from (63) We use the Eqs. (B7) and (B8) with the identification α = 1 and F (k 0 , k) = 1, which yield To proceed, we consider some useful integrals (B14) By using the integral (B12) for r = 2, we are able to solve the time integral in (B11) arriving at Introducing the new variable ε = t(1 − x), we write We use the approximation (B20) The function AppellF 1 (a, b 1 , b 2 , c, x 1 , x 2 ) that appear above, is the hypergeometric function of two variables defined by AppellF 1 (a, b 1 , b 2 where (q) n = q(q + 1) . . . (q + n − 1), or equivalently (q) n = Γ(q + n) Γ(q) .
We continue with f n 1 , given in (66), and again promote to d dimensions and using our Eqs. (B7) and (B8), with α = 1 and F (k 0 , k) = k 0 − g 2 k 2 0 , we write with l (1) 2 Using the integrals (B12) and (B13) we arrive at Just as in the previous calculation, we define ε = t(1 − x) and consider the approximation (B17) with n = 3 for the first and second term and n = 1 for the third and fourth term in (B26). Together with this, we replace d = 4 in all the finite terms and consider the identity 1− x − t(1−x) With ε = t(1 − x), k 2 = 1 g 2 2 (t 2 − g 2 2 m 2 ) and using the approximation (B17) gives