Pauli–Zeldovich cancellation of the vacuum energy divergences, auxiliary fields and supersymmetry

We have considered the Pauli–Zeldovich mechanism for the cancellation of the ultraviolet divergences in vacuum energy. This mechanism arises because bosons and fermions give contributions of the opposite signs. In contrast with the preceding papers devoted to this topic wherein mainly free fields were studied, here we have taken their interactions into account to the lowest order of perturbation theory. We have constructed some simple toy models having particles with spin 0 and spin 1 / 2, where masses of the particles are equal while the interactions can be quite non-trivial.


Introduction
Many years ago Pauli [1] suggested that the vacuum (zeropoint) energies of all existing fermions and bosons compensate each other. This possibility is based on the fact that vacuum energy of fermions has a negative sign whereas that of bosons has a positive one. As is well known, such a cancellation indeed takes place in supersymmetric models (see e.g. [2]). Subsequently in a series of papers Zeldovich [3,4] related vacuum energy to the cosmological constant. However rather than eliminating divergences through the bosonfermion cancellation, he suggested the Pauli-Villars regularization of all divergences by introducing a number of massive regulator fields. Covariant regularization of all contributions then leads to finite values for both the energy density ε and (negative) pressure p corresponding to a cosmological constant, i.e. connected by the equation of state p = −ε. In our preceding paper [5] we examined the conditions for the cancellation of the ultraviolet divergences of the vacuum energy to the leading order inh, i.e. by considering free theories and neglecting interactions. Such conditions are reduced to some sum rules involving the masses of particles present in the model. We formulated these conditions not only for the Minkowski spacetime, but also for the de Sitter one. In the latter case, the radius of the de Sitter universe also enters into the mass sum rules. In paper [6] we applied such considerations to observed particles of the Standard Model (SM) and also studied the finite part of vacuum energy. This last contribution should be very small, so as to obtain a result compatible with the observed value of the cosmological constant (almost zero with respect to SM particle masses). We showed [6] that it was impossible to construct a minimal extension of the SM by finding a set of boson fields which, besides canceling ultraviolet divergences, could compensate residual huge contribution of known fermion and boson fields of the Standard Model to the finite part of the vacuum energy density.
On the other hand, we found that addition of at least one massive fermion field was sufficient for the existence of a suitable set of boson fields which would permit such cancellations and obtained their allowed mass intervals. On examining one of the simplest SM extensions satisfying the constraints, we found that the mass range of the lightest massive boson was compatible with the Higgs mass bounds which were known at the time of the publication of the paper [6]. As is well known, later the Higgs boson was discovered at the LHC [7,8]. For some time it appeared that there might exist an the observed diphoton excess at 750 GeV [9]. This excess, had it been confirmed, could be interpreted as an indication for the existence of a new heavy elementary or composite particle with a mass of the order of 750 GeV. Later this phenomenon disappeared, nonetheless inspiring in the meanwhile quite a few theoretical works. In particular, we also studied in our preprint how the presence of such a new particle could be included into our scheme of the cancellation of ultraviolet divergences of vacuum energy [10].
In our preceding papers [5,6,10] we studied only free theories without interactions. It is also interesting to take interactions into account, at least to the lowest order of perturbation theory. This is not easy, and in the present paper we shall concentrate on the construction of relatively simple toy models where the "Pauli-Zeldovich cancellation" of ultraviolet divergences still takes place.
We wish to emphasize that the approach employed in the present paper represents a whole direction in quantum field theory which goes well beyond effective low energy field theory and, although based on some hypothesis, has not been proven to be wrong. We also wish to mention the paper by Ossola and Sirlin [11], where contributions of fundamental particles to the vacuum energy density were discussed with a special attention to relations between different regularization schemes and to the appearance of power divergences in different contexts. Other related approaches are presented in Refs. [12][13][14].
In the recent paper [15], it was noticed that under certain circumstances (in particular, but not limited to finite QFTs), the Pauli cancellation mechanism would survive the introduction of particle interactions. It was pointed out there that for the mass sum rules to be valid at different mass scales, it is necessary to impose some relations on mass runnings with energy. Thus, the corresponding relations between anomalous mass dimensions were formulated [15]. However, concrete examples were not constructed.
In the present paper we discuss some relatively simple examples of models where the Pauli-Zeldovich cancellation takes place to the first order of perturbation theory. Being inspired by the famous supersymmetric Wess-Zumino model [16], we consider models with spinor, scalar and pseudoscalar fields only. We hope to treat vector (gauge) fields in future works. The models which we discuss are not supersymmetric, but they have one important feature which makes them akin to supersymmetric models: the number of the fermion and boson degrees of freedom in them is the same. That implies an unexpected feature: the necessity to take so called auxiliary fields into account. Such fields are necessary in the supersymmetric models, because they allow one to formulate supersymmetry transformations in a coherent way.
But their role is even more ubiquitous. To conserve supersymmetry, it is necessary to have the balance between fermion and boson degrees of freedom not only on shell, but also off shell. However, the number of degrees of freedom of a spinor field doubles when it is off shell. For example, a Majorana spinor has two complex components, i.e. four degrees of freedom off shell. When we require the satisfac-tion of the first-order Dirac equation, the number of degrees of freedom becomes equal to two. Thus, for example, in the Wess-Zumino [16] model one has two fermion degrees of freedom of the Majorana spinor and two boson degrees of freedom associated with the scalar and pseudoscalar fields. Off shell the number of fermion degrees of freedom becomes equal to four, while the role of two additional boson fields is played by two auxiliary fields which become in a sense independent off shell. If we consider non-supersymmetric models with the Pauli-Zeldovich mechanism of cancellation of ultraviolet divergences for vacuum energy in the presence of interactions, then the number of the boson and fermion degrees of freedom should be equal not only on shell, but also off shell. This means that we should introduce auxiliary fields. Further, when we consider a model with interactions, we should not only take into account running of masses of the fields, but also consider cancellations of contributions coming from the potential terms in the Lagrangians. It is there that the introduction of the auxiliary fields becomes very convenient. Fortunately, we shall see that, at least in the considered class of spinor-scalar models, the introduction of auxiliary fields is equivalent to a simple rule for the calculation of some contribution to the scalar fields self-interaction.
Here we can add that, in principle, one can perform all calculations and show that in the formalism where auxiliary fields are excluded, vacuum energy in the supersymmetric models is equal to zero. However, in this case there are no separate cancellations of the potential energy and of the kinetic energy between bosons and fermions. Thus, verification of the analogous cancellation in non-supersymmetric models becomes more complicated. Hence, it is better to implement rather simple rules, equivalent to the explicit introduction of auxiliary fields, which will be used in the present paper.
Here we present a model consisting of a Majorana fermion and two scalar fields with the same mass and with different kinds of interactions, and we show that for such a model, one can find a family of coupling constants such that the Pauli-Zeldovich mechanism for the cancellation still works. Then we find an analogous family of models with a Majorana fermion, a scalar field and a pseudoscalar field. Obviously, the Wess-Zumino model belongs to this family. We also discuss briefly models where particles with different masses are present.
The structure of the paper is as follows: in the second section we briefly discuss the mass sum rules for theories without interactions; in the third section we formulate rules for the conservation of the mass sum rules when interactions are switched on. In Sect. 4 we discuss the vacuum expectation values of the potential terms and the role of auxiliary fields. In Sect. 5 we present a model with one Majorana field and two scalar fields. In the sixth section we consider a model with one Majorana field, one scalar field and one pseudoscalar field. Section 7 is devoted to the discussion of models with non-degenerate masses, the last section contains some concluding remarks.

Vacuum energy and the balance between the fermion and boson fields
One knows that vacuum energy of the harmonic oscillator is equal toh ω 2 . If one has a massive field with mass m, then ω = √ k 2 c 2 + m 2 c 4 , where k is the wave number. In the following we shall seth = 1 and c = 1. The energy density of vacuum energy of a scalar field treated as free oscillators with all possible momenta is given by the divergent integral [3]: We can regularize this integral by introducing a cutoff Λ. In this case On expanding this expression with respect to the small parameter m Λ , one obtains The contribution of one fermion degree of freedom coincides with that of Eq. (1) with the opposite sign. It now follows from Eq. (3) that to cancel the quartic ultraviolet divergences proportional to Λ 4 , one has to have equal numbers of boson and fermion degrees of freedom: The conditions for the cancellation of quadratic and logarithmic divergences are and respectively. Here the subscripts S, V and F denote scalar, massive vector and massive spinor Majorana fields respec-tively (for Dirac fields it is sufficient to put a 4 instead of 2 on the right-hand sides of Eqs. (5) and (6)). For the case in which the conditions (4), (5) and (6) are satisfied, the remaining finite part of the vacuum energy density is equal to Let us now calculate the vacuum pressure. This pressure is given by the formula [3] On introducing the cutoff Λ we have On expanding this expression with respect to the small parameter m Λ , we obtain On then comparing the expressions (3) and (10), we see that the quartic divergence satisfies the equation of state for radiation p = 1 3 ε, the quadratic divergence satisfies the equation of state p = − 1 3 ε, which sometimes is identified with the so called string gas (see e.g. [17,18]), while the logarithmic divergence behaves as a cosmological constant with p = −ε. If all these divergences cancel, then the finite part of the pressure is which also behaves as a cosmological constant. We can emphasize that the real (renormalized) vacuum energy does not have radiation-like (or, p = −ρ/3) equation of state. Instead, following Pauli and Zeldovich we investigate conditions under which it (and the whole average value of the vacuum energy-momentum tensor) can be finite without any renormalization in the absence of an exact supersymmetry (which is not realized in nature). If it is finite indeed, then any method of calculation -covariant or non-covariant -will show it.

Running masses and anomalous mass dimensions
In what follows we shall consider models having only particles with spin zero and spin 1/2. If we include the interactions, the masses begin their running and the conservation of the relations (5) and (6) implies some new restrictions on the masses and on the coupling constants. Namely, the conservation of the relation (5) gives where γ m is the mass anomalous dimension defined as where as usual μ is the renormalization mass parameter. The conservation of the relation (6) gives These relations coincide with those presented in paper [15]. We shall here derive the expressions for these anomalous mass dimensions. Generally the technique of such calculations was developed many years ago [19][20][21][22][23][24][25]. However, for convenience and completeness we shall perform all the derivations from the start. On considering our toy models with degenerate masses, we shall not really use them explicitly. It will be enough to study shifts of masses induced by radiative corrections for different fields present in the models under consideration. However, when one considers models where particles with different masses are present, the formulas given in this section become necessary.
Our treatment of the anomalous mass dimensions in the presence of quadratic divergences is based on the approach presented in paper [26], which in turn uses the version of renormalization group formalism connected with dimensional regularization [27].
Let us consider the model, including a Dirac spinor with a mass M and a scalar field with a mass m.
The full propagator of the fermion field is given by whereÂ ≡ γ μ A μ and where Σ is the self-energy operator of the fermion field. This operator in the one-loop approximation is given by the formula Here d is the dimensionality of the spacetime such that Let us first calculate the term Σ 2 : On making a Wick rotation, we obtain the integral on the Euclidean momenta: We are interested only in the divergent part of this integral. Thus, we can neglect the masses in the denominator. On using the formula the Gaussian integration, the formula connecting the Euler B and Γ functions and the fact that we arrive to the expression To find Σ 1 , we shall take the 1 4 Tr(pΣ). Then Using the identity we transform the expression (25) as This expression only contains a logarithmic divergence. On making a Wick rotation, integrating in the Euclidean momentum space and keeping only the poles in ε, we obtain Thus, On substituting the formula (29) into Eq. (16) we see that the fermion propagator in the one-loop approximation is (30) In the same approximation, this propagator can be rewritten as Thus, we see that the shift of the mass M is To compensate this shift, we should introduce a counterterm into the Lagrangian, or in other terms, we should introduce a bare mass M B which is connected with the renormalized mass M through the relation where Further, to have a canonically normalized fermion field, or in other words, to compensate a non-trivial divergent factor in the numerator of the formula (31), we should introduce a bare fermion field where On now, following the scheme elaborated in paper [27], we introduce and calculate the anomalous mass dimension for the fermion mass M. Let us remember that when we use the dimensional regularization, the renormalized quantities depend on the renormalization mass parameter μ. At the same time the bare quantities depend on the regularization parameter ε, but do not depend on the renormalization mass parameter μ. Thus, we can write down a general equation Generally, the renormalization constant Z M has the following structure: On introducing we can rewrite Eq. (37) as follows: For the case wherein the residues a n depend only on the Yukawa coupling constant g, Eq. (40) becomes We now introduce the β -function for the Yukawa constant g: Then Eq. (41) reads: The above equation should be correct in any order in ε. To the zeroth order it gives: From Eq. (34) we immediately obtain The calculation of the analogous quantity for the scalar field is more complicated because the mass renormalization in this case includes quadratic divergences. To treat them, we shall follow the approach developed in paper [26]. The full propagator of the scalar field is where Π is the self-energy operator. The contribution of the scalar field self-interaction to the one-loop order in the operator Π is Let us note that we here include the factor μ ε to provide the correct dimensionality of Π . We did not include such a factor on calculating the self-energy of the fermion, because there only logarithmic divergences were present and this factor disappeared in the limit d → 4. Here, in the presence of quadratic divergences the factor μ ε becomes crucial. A direct calculation gives One can see that this expression has the pole at d = 4 and also the pole at d = 2, corresponding to quadratic divergence [26]. Indeed, it is well known that in the theory with the Lagrangian (15) the index of divergence of a diagram G, ω (G) is where E F is a number of the external fermion lines and E B is a number of the external boson lines. Thus, the diagrams with E F = 0, E B = 2 are quadratically divergent. Let us now consider d-dimensional spacetime. In this case, the formula (49) is replaced by where L is the number of loops. Let us again consider a diagram with E F = 0, E B = 2. This diagram, which is quadratically divergent at d = 4 becomes logarithmically divergent (ω(G) = 0) at d = 4 − 2 L . That means that the quadratic divergence is represented as a pole of the quantity and in the case of the one-loop approximation Thus, expanding the expression (48) around d = 4, we have while expansion of the same expression at d → 2 gives Thus, the infinite shift of the mass squared in the full propagator of the scalar field due to its self-interaction is We can analogously calculate the contribution of the fermion loop to the self-energy of the scalar field.
The scalar field propagator corrected by the fermion loop is (1) .
Thus, the mass squared of the scalar field is shifted as On combining the last equation with Eq. (55), we obtain the full mass shift: To compensate this shift, we introduce a bare scalar field mass following the procedure elaborated in the paper [26]: where the renormalization constants in the one-loop approximation are and (1) .
Further, to have a canonical normalization of the scalar field, we introduce a bare field as follows: where On now introducing the anomalous mass dimension and requiring the independence of the bare mass (63) on the renormalization mass parameter μ and using the explicit expressions (64) and (65), we obtain in the one-loop approximation Let us also include some pseudoscalar fields into our model. The interaction between a scalar and the fermion is described by the following term in the Lagrangian: Here On repeating the preceding calculations and using the formulae (71), we see that the contribution of the interaction (70) to the anomalous mass dimension of the fermion field is The contribution of the fermion loop to the anomalous mass dimension of the pseudoscalar field χ is

Contribution of potential terms into the vacuum energy and the auxiliary fields
When we switch on the interactions and require the cancellation of ultraviolet divergences in the expression for vacuum energy, we should consider not only the mass sum rules, but also the potential terms. The contributions of the potential terms to the vacuum energy density have the following structure Here, L int is the interaction Lagrangian and the exponent should be expanded up to necessary order in the perturbation theory while V represents potential terms. For the term to obtain the result in the two-loop approximation which we study in the present paper, it is enough to take only the zeroth order of the expansion of the exponent of the action in the formula (74). The corresponding contribution is equal to where the integral I is defined as The contribution of the Yukawa interaction term is given by the structure where the second factor (−ig)ψψφ comes from the firstorder term in the expansion of the T -exponent. This contribution (for the case of a Majorana spinor) is equal to where M is the fermion mass and m is the scalar mass. A simple calculation shows that for the case of the Wess-Zumino model, when m = M and there are well-known relations between the coupling constants [16], the quartic divergences present in the contributions (76) and (79) do not cancel each other (we shall present detailed calculations in the next section). Namely, the contribution of the spinors is twice that of the scalars. The reason for this mismatch was already discussed in the Introduction. The point is that the number of fermion degrees of freedom is doubled off shell. To compensate this effect, we should introduce the auxiliary scalar fields as is done in supersymmetric models. A simple example shows that this exactly gives the doubling of the leading contribution to vacuum energy. Indeed, let us consider a model with the Lagrangian On shell this theory is equivalent to the theory where the auxiliary field is excluded F by means to the equation of motion and which has a Lagrangian It follows from Eq. (76) that vacuum energy in the theory with the Lagrangian (82) is equal to Let us calculate an analogous (two-loop) contribution to vacuum energy in the model with the Lagrangian (80). It is equal to (84) Here we have used the fact that the propagator of the auxiliary field in the massless theory is given [28] by We see that in this case the result is doubled because of the effective doubling of the number of degrees of freedom. One can check also that the contributions to the self-energy operator of the scalar field to the order of h 2 coincide in the models with the Lagrangians (80) and (82): Thus, the requirement of the explicit account of auxiliary fields arises only in the diagrams possessing quartic ultraviolet divergences and including only boson propagators, because their contribution is proportional to the number of degrees of freedom present off shell in the model under consideration. This fact gives us a practical recipe: when one calculates vacuum energy contribution of the scalar field diagrams, having the shape of "eight", one should multiply it by the factor 2.
Concluding this section, we wish to make one more comment. In the action (80) the term F 2 2 is also present. One can consider this term as a part of the kinetic energy. The contribution of this term into vacuum energy is given by the formula Thus, on summing (87) and (84) we reproduce the result (83). This means that, in the end, the results for vacuum energy in the model (80) with an auxiliary field and in the model (82), where the auxiliary field is eliminated, coincide. However, the expressions for the contributions of the potential energy and of the kinetic energy do not coincide separately.
As we have already mentioned in the Introduction, a similar effect can be observed in the supersymmetric Wess-Zumino model. If we consider the formalism in the absence of auxiliary fields, vacuum energy is still equal to zero, but the potential and kinetic energy are not equal to zero separately.

A model with one Majorana and two scalar fields
Let us consider a model with a Majorana field ψ and two scalar fields A and B. All the fields have the same mass m and the interaction is given by The two tadpole diagrams for fields A and B should be cancelled to avoid the necessity of introducing linear in fields terms into the Lagrangian. The tadpole for the field A arises due to the contraction of this field with the vertices A 3 , AB 2 andψψ A. All these contributions are proportional to the integral (77). The corresponding combinatorial factors are 3mh 1 , for A 3 , mh 4 for AB 2 and −4mg 1 for the vertexψψ A.
The last contribution arises due to the trace of the fermion propagator which is proportional to the mass m. Thus, the cancellation of the tadpole diagram for A requires Similarly the vanishing of the tadpole for the field B requires Now the self-energy operator for the propagator of the field A obtains the contributions from the vertex A 4 , from the vertex A 2 B 2 and from the pair of vertexesψψ A, A 3 , A 2 B and AB 2 . The contributions of two quartic vertexes are both proportional to the integral I . The corresponding coefficients are 12λ 1 and 2λ 3 . The contribution of the fermion loop is where the factor 2 arises due to the Majorana nature of the fermion. Then where Quadratic divergences present in the integral I should be canceled because such divergences do not arise in the selfenergy correction to the fermion propagator. Thus, we have and, analogously, Now the contribution of the pairs of the triple scalar vertices is Thus, the propagator of the scalar field A in the one-loop approximation has the form Normalizing as usual the wave function, i.e. making the coefficient at p 2 equal to 1, we obtain Thus, this effective shift of the mass squared given by defines the running of the mass for the scalar field A. The analogous shift for the second scalar field is The self-energy contribution to the fermion propagator is where the factor 4 arises due to the Majorana nature of the fermion. The fermion propagator is now On normalizing the term atp, we obtain The shift of the mass squared is The running of the masses and, hence, the shifts (99), (100) and (104) should be equal and we obtain two equations: and 18h 2 2 + 4h 2 4 + 2h 2 3 − 12g 2 2 = 12(g 2 1 + g 2 2 ).
Let us now consider the contribution of the potential term (88) to vacuum energy. The contribution of the quartic terms is The contribution coming from the two scalar-fermion vertices is given by the integral where The contribution to vacuum energy of the triple scalar interactions is We can now observe that the sum provided Eqs. (94), (95), (105) and (106) are satisfied. The coefficient 2 in front of the term E 1 is introduced to take into account the fact that the number of boson and fermion degrees of freedom should be equal also off shell. It is equivalent to the introduction of two auxiliary fields in supersymmetric models, as was explained in the preceding section. On now substituting the expressions for g 1 and g 2 from Eqs. (89) and (90) into Eqs. (105) and (106), we obtain the following pair of the consistency conditions on the constants h 1 , h 2 , h 3 and h 4 : These equations are homogeneous in h 1 , h 2 , h 3 and h 4 . Thus, we can fix h 1 = 1 and we can then change the value of h 2 . Then we shall have a system of two quadratic equations for h 3 and h 4 . This system, which is equivalent to one quartic equation for one variable, is solvable analytically but the solutions are very cumbersome. Thus, we shall just present some numerical solutions. In any case we have four solutions, two of them are complex and two of them are real. We shall only take real solutions into account. Let us note that the negative values of the coupling constants are not essential because they are in front of the odd (third) power of fields A and B. The lower bound of the scalar field potential exists and is determined by the quartic terms with positive constants λ 1 , λ 2 and λ 3 .
One can meanwhile introduce the quartic interactions using auxiliary fields in a manner similar to that used in the Wess-Zumino model. It is enough to introduce the following terms into the Lagrangian instead of terms with quartic interactions: 6 Model with a Majorana field, a scalar field and a pseudoscalar field Let us consider another toy model where the field B is a pseudoscalar. In this case h 2 = h 3 = 0 and the interaction between the pseudoscalar and the fermion is described by the Lagrangian In this case we have only one condition for the tadpole cancellation for the scalar field A which coincides with that given by Eq. (89). The conditions for the cancellation of quadratic divergences in the propagators of the scalar and pseudoscalar fields are also the same (94) and (95). However, the shifts of the mass squared for the fields A, B and ψ are different. They are proportional to respectively. Correspondingly, on requiring that the running of these three masses are the same, we obtain the following couple of equations: and From these two equations we obtain immediately If g 1 = h 1 then from Eq. (89) we obtain which, in turn, implies If, instead, then which implies a negative value for g 2 2 or for g 2 1 as follows from the couple of equations (116) and (117). Thus, the choice (119) should be discarded.
We have seen that for the case with one Majorana field, one scalar and one pseudoscalar we have less freedom in the choice of the coupling constants than in the case of two scalar fields and one Majorana field, but this choice is still broader than that in the Wess-Zumino model.
In principle, one can also consider a model with one Majorana field and two pseudoscalar fields. In this case the triple scalar interactions do not exist and the constants h 1 , h 2 , h 3 and h 4 all are equal to zero. The requirement of the absence of quadratic divergences in the shifts of the mass squared of two pseudoscalar propagators is again given by Eqs. (94) and (95). However, the shifts of the mass squared for this propagators are equal to zero, while the shift of the mass squared of the fermion propagator is given by Thus, the Yukawa interactions should vanish as well. Now, the conditions (94) and (95) can be satisfied if but the corresponding quartic potential of these two pseudoscalar fields is unbounded from below and is hardly interesting.

Models with non-degenerate masses
It is interesting to find toy models with masses which are not degenerate. In this case it is necessary to consider at least four boson and four fermion degrees of freedom [5]. The simplest models of this kind are those which include a certain number of "triplets" of the types described in two preceding sections, i.e. with degenerate masses inside any triplet and with coupling constants (again, describing interactions within a triplet) which satisfy the relations obtained in the Sects. 5 and 6. Naturally, in this case, if there are no interactions between the fields belonging to different triplets, then the Pauli-Zeldovich mechanism does work. If we introduce interactions between different triplets with different masses, then the coupling constants should satisfy some constraints. We shall illustrate this by considering a model, where there are two triplets. Both triplets contain a Majorana fermion and two scalar fields. The mass of all the particles in the first triplet is equal to m 1 , in the second tripletm 2 and the coupling constants describing interactions within the triplets are chosen in such a way that vacuum energy is equal to zero. Let us then introduce the following interaction Hamiltonian between fields belonging to different triplets: H = λ AC A 2 C 2 + λ AD A 2 D 2 + λ BC A 2 C 2 + λ B D A 2 D 2 +g Aχχ χ A + g Bχχ χ B + g Cψψ ψC + g Dψψ ψ D where the scalar fields A and B and the Majorana spinor ψ belong to the first triplet, while the scalar fields C and D and the Majorana spinor χ belong to the second triplet.
It is now possible to find relations which constrain the choice of the coupling constants given in the interaction Hamiltonian (120). However, this task is rather cumbersome and we shall postpone it and the construction of other models for future work [29].

Concluding remarks
In this paper we have studied the Pauli-Zeldovich mechanism for the cancellation of ultraviolet divergences in vacuum energy which is associated with the fact that bosons and fermions produce contributions to it having opposite signs. In contrast with the preceding papers devoted to this topic where only free fields were considered, here we have taken interactions up to the lowest order of perturbation theory into account. We have constructed a number simple toy models having particles with spin 0 and spin 1/2, wherein masses of the particles are equal while interactions can be quite non-trivial. To make calculations simpler and more transparent, it was found useful to introduce some auxiliary fields. It appears that the presence of these fields is equivalent to the modification of some contributions of the physical fields to the vacuum expectation of the potential energy. We hope to construct more complicated models including particles with different masses and in the presence of vector fields in future work [29].