## Services on Demand

## Journal

## Article

## Indicators

- Cited by SciELO
- Access statistics

## Related links

- Cited by Google
- Similars in SciELO
- Similars in Google

## Share

## Brazilian Journal of Physics

##
*Print version* ISSN 0103-9733

### Braz. J. Phys. vol.40 no.4 São Paulo Dec. 2010

#### http://dx.doi.org/10.1590/S0103-97332010000400014

**Radiatively induced non linearity in the Walecka model **

**Rafael Cavagnoli ^{I,}^{*}; Marcus Benghi Pinto^{II,}^{†}**

^{I}Departamento de Física, Universidade Federal de Santa Catarina, 88040-900 Florianópolis, Santa Catarina, Brazil and Centro de Física Teórica - Dep. de Física - Universidade de Coimbra - P-3004 - 516 - Coimbra - Portugal

^{II}Departamento de Física, Universidade Federal de Santa Catarina, 88040-900 Florianópolis, Santa Catarina, Brazil

**ABSTRACT**

We evaluate the effective potential for the conventional linear Walecka non perturbatively up to one loop. This quantity is then renormalized with a prescription which allows finite vacuum contributions to the three as well as four 1PI Green's functions to survive. These terms, which are absent in the standard relativistic Hartree approximation, have a logarithmic energy scale dependence that can be tuned so as to mimic the effects of ø^{3} and ø^{4} type of terms present in the non linear Walecka model improving quantities such as the compressibility modulus and the effective nucleon mass, at saturation, by considering energy scales which are very close to the nucleon mass at vanishing density.

**Keywords: **Renormalization; Nuclear matter; Finite densities, Vacuum corrections; Walecka model.

**1. INTRODUCTION**

Quantum hadrodynamics (QHD) is an effective relativistic quantum field theory, based on mesons and baryons, which can be used at hadronic energy scales where the fundamental theory of strong interactions, quantum chromodynamics (QCD), presents a highly nonlinear behavior. The Walecka model [1] to be considered here represents QHD by means of a Lagrangian density formulated so as to describe nucleons interacting through the exchange of an isoscalar vector meson (ω) as well as of a scalar-isoscalar meson (ø) which is introduced to simulate intermediate range attraction due to the *s*-wave isoscalar pion pairs. The original Walecka model (QHD-I) is described by the Lagrangian density

where Ψ, ø and ω denote respectively baryon, scalar and vector meson fields (with the latter being coupled to a conserved baryonic current). The term *U*(ø,*V*), which describes mesonic self interactions was set to zero in the original model so as to minimize the many body effects while the term _{CT} represents the counterterms needed to eliminate any potential ultra violet divergences arising from vacuum computations.

It is important to recall that, roughly, counterterms are composed by two distinct parts the first being a divergent piece which exactly eliminates the divergence resulting from the evaluation of a Green function at a given order in perturbation theory. The second piece is composed by a finite part which is *arbitrary* and can be fixed by choosing an appropriated renormalization scheme [2].

The important parameters are the ratios of coupling to masses, C^{2}_{s} and C^{2}_{v}, with C_{i}^{2} = (*g _{i}M/m_{i}*)

^{2}which are tuned to fit the saturation density ρ

_{0}= 0.193 fm

^{-3}and binding energy per nucleon, BE = -15.75 MeV [1]. However, the QHD-I predictions for some other relevant static properties of nuclear matter do not agree well with the values quoted in the literature. For example, using the Mean Field Approximation (MFA) which considers only in medium contributions at the one loop level one obtains that, at saturation, the effective nucleon mass is

*M**

_{sat}~ 0.56 M, which is somewhat low, while the compression modulus, K ~ 540 MeV, is too high according to the accepted values in the literature [3-5]:

*M**

_{sat}/M ~ 0.70 to 0.80 and K ~ 200 MeV to 300 MeV. Let us point out that the effective nuclear mass at saturation density is not known accurately. According to the type of data and its analysis, the effective mass is defined in different ways, and so this mass usually is not that of the nuclear field theory, although they might be close and related to each other.

In principle, since this is a renormalizable quantum field theory, vacuum contributions (and potential ultra violet divergences) can be properly treated yielding meaningful finite results. These contributions were first considered by Chin [6], at the one loop level, in the so called Relativistic Hartree Approximation (RHA) which produced a more reasonable value for the effective mass, *M**_{sat} ~ 0.72 M. However, the compression modulus remained at a high value, K ~ 470 MeV.

One could then try to improve the situation by also considering exchange contributions since both, MFA and RHA, consider only direct terms in a nonperturbative way. When vacuum contributions are neglected this approximation is known as the Hartree-Fock (HF) approximation producing *M**_{sat} ~ 0.53 M and K ~ 585 MeV [1]. By comparing the results from MFA, RHA and HF one sees how vacuum effects can improve the values of *M**_{sat } and K. Then, the natural question is if the situation could be further improved by considering the vacuum in HF type of evaluations. The main concern now being the difficulty to deal with overall, nested, and overlapping type of divergences which surely arise due to the self consistent procedure.

The complete evaluation of vacuum contributions within direct and exchange terms was performed by Furnstahl, Perry and Serot [7] (see Ref. [8] for an early attempt in which only the propagators have been fully renormalized ). This cumbersome calculation considers the nonperturbative evaluation and renormalization of the energy density up to the two loop level showing the non convergence of the loop expansion. Later, the situation has been addressed with the alternative Optimized Perturbation Theory (OPT) which allows for an easier manipulation of divergences [9]. Two loop contributions have been evaluated and renormalized in a perturbative fashion with nonperturbative results further generated via a variational criterion. However, saturation of nuclear matter could not be achieved with the results supporting those of Ref. [7].

Meanwhile, the compressibility modulus problem has been circumvented by introducing some more parameters, in the form of new couplings [10], to the original Walecka model. One then considers *U*(ø,*V*) appearing in Eq. (1) as

where *k* = 2*bMg _{s}*

^{3}and λ = 6

*cg*

_{s}^{4}. This version of QHD is known as the nonlinear Walecka model (NLWM) and the main role of the two additional parameters,

*b*and

*c*, is to bring the compression modulus of nuclear matter and the nucleon effective mass under control. However, one may object to this course of action since the mesonic self interactions will increase the many body effects apart from increasing the parameter space. Notice also that terms like ø

^{n}(

*n*

__>__5) are not allowed since then, in 3+1 dimensions, one would need to introduce coupling parameters with negative mass dimensions spoiling the renormalizability of the original model apart from increasing the parameter space.

Heide and Rudaz [11] have then realized that it is still possible to keep *U*(ø,*V*) = 0 while improving both *K* and *M**_{sat}. The key ingredient in their approach is related to the complete evaluation (regularization and renormalization) of divergent vacuum contributions. Regularization is a formal way to isolate the divergences associated with a physical quantity for which many different prescriptions exist, e.g., sharp cut-off, Pauli-Villars, and Dimensional Regularization (DR). Within DR, which was used by Chin, one basically performs the evaluations in d - 2ε dimensions taking ε → 0 at the end so that the ultra violet divergences show up as powers of 1/ε. However, to keep the dimensionality right when doing d → d-2ε one has to introduce *arbitrary* scales with dimensions of energy (Λ, or the related^{1} ). Chin has chosen a renormalization prescription in which the final results do not depend on the arbitrary energy scale while Heide and Rudaz chose one in which such a dependence remains, as in most QCD applications. Since the latter authors also worked at the one loop level their approximation became known as the Modified Relativistic Hartree Approximation (MRHA) and their main result was to show that it is possible to substantially improve *K* and *M**_{sat} by suitably fixing the energy scale, Λ. Moreover by choosing Λ = *M* the MRHA recovers RHA. In connection with neutron stars, the MRHA has been applied to the Walecka model in Refs. [12]. The dependence of QHD results on the choice of renormalization conditions, which may be expected in phenomenological models, has been shown long ago (see Ref. [13] and references therein). At the same time, it has been discussed that approximations of the pure Walecka model are inconsistent leading to an unstable ground state [14] and ignoring the usual nonlinear contributions due to renormalization carries the danger to reenter this instability also in the quantum corrected Walecka model. In principle these facts will spoil any rigorous attempt to include vacuum corrections, in the usual Walecka model, at arbitrarily high orders. Very often however, this model is still being treated at the one loop level with the MFA and the aim of the present paper is to show that in this case the values of *K *and *M**_{sat} can be improved simply by evaluating vacuum corrections up to an energy scale very close (but not equal) to the nucleon mass at zero density.

Therefore, even if our procedure may be spoiled at higher orders by the facts mentioned above we believe that for it remains a useful, and easy to implement, procedure which improves MFA results without the need for new paremeters in the theory. This is an important feature to be considered. Secondly, as we shall see, models which lead to low effective masses at saturation are not suitable for neutron stars calculations, another important application of RMF (Relativistic Mean Field) models.

Also, one of our goals is to treat the Walecka model using a formalism which is closely related to the one used in QCD and other modern quantum field theories. Within the QHD model, the temperature and density are usually introduced using the real time formalism employed in the original work of Walecka. Instead, we use Matsubara's imaginary time formalism treating the divergent integrals with DR adapted to the modified minimal subtraction renormalization scheme [2] which constitute the framework most commonly used within QCD. To obtain the ground state energy density, e, we will first evaluate the effective potential (or Landau's free energy), , whose minimum gives the pressure, *P*. By choosing appropriate renormalization conditions we generate effective three- and four-body couplings, in , which are not present at the classical level. As we shall see the numerical values of these effective couplings run with the energy scale, , allowing for a good tuning of K and *M**_{sat} which have their values improved at energy scales of about 0.92 *M*-0.98 *M* (*M* = 939 MeV) while the usual RHA results are retrieved for the choice = *M*.

The MRHA proposed by Heide and Rudaz suggests that if one seeks to minimize many-body effects in nuclear matter *at saturation*, the choice Λ ≅ *M**_{sat} is the necessary one. Our philosophy is slightly different and perhaps simpler to implement. Since possible modifications in the behavior of *K* and *M** seem to be dictated by the presence of k_{eff}ø^{3} and λ_{eff}ø^{4} type of terms we shall use the Chin-Walecka renormalization prescription to deal with ø^{n} (n = 0,1,2) vacuum contributions by requiring that their respective contributions vanish at zero density (a requirement which was also adopted within the MRHA). However, as far as the vacuum contributions related to ø^{3} and ø^{4} are concerned we advocate that one only needs to keep the finite energy scale dependent parts in the effective three- and four-body couplings. In this way, not only k_{eff} and λ_{eff} run with but, as we shall see, one also retrieves the RHA at = M.

Considering the effective potential at zero density we will choose an appropriate renormalization prescription for this particular model. Since our main goal is to improve K and *M* _{sat}* by quantically renormalizing k = 0 → k

_{eff}( ) and λ = 0 → λ

_{eff}( ) we can keep M,m

_{s},m

_{v}as representing

*the*vacuum physical masses for simplicity. This choice means that, at k

_{F}= 0, all mass parameters (

*M,m*) represent the effective vacuum masses and shall not run with as opposed to k

_{s},m_{v}_{eff}and λ

_{eff}. In theories such as QCD the running of the couplings is dictated by the β function whose most important contributions come from the so-called leading logs, e.g. ln( /

*M*), which naturally arise in DR evaluations. The application of renormalization group (RG) equations to the effective Walecka model is beyond the scope of our work

^{2}. Nevertheless, our renormalization prescription to obtain a scale dependence so as to better control

*K*and

*M**

_{sat}is inspired by the leading logs role in the b function and the the renormalization scheme presented here proposes that one should preserve only the scale dependent leading logs which appear in the expressions for k

_{eff}ø

^{3}and λ

_{eff}ø

^{4}. As a byproduct, and contrary to the MRHA case, both quantities will display the same scale dependence. Here, this approximation will be called the Logarithmic Hartree Approximation (LHA). The numerical results show that the best LHA predictions for

*K*and

*M**

_{sat}according to the literature are obtained at energy scales which are only about 5 % smaller than that of the RHA, that is ≅ 0.95

*M*. This is a nice feature since the values of the energy scale and that of the highest mass in the spectrum are almost the same whereas in the MRHA the optimum scale, set to be close to

*M**

_{sat}is about 35 % smaller than

*M*. From the quantitative point of view, the LHA produces better results than the MRHA as will be shown.

The work is presented as follows. In the next section the one loop free energy is evaluated using Matsubara's formalism. The renormalization of the vacuum contributions is discussed in Section III and the complete renormalized energy density is presented in Section IV. Numerical results and discussions appear in Section V while our conclusions are presented in Section VI. For completeness, in the appendix, we discuss a case in which *m _{s}* does not represent the physical mass.

**2. THE FREE ENERGY TO ON LOOP**

In quantum field theories the effective potential (or Landau's free energy), , is defined as the generator of all one particle irreducible (1PI) Green's functions with zero external momentum. The standard textbook definition (for one field, ø) reads [2]

where we have absorbed non relevant factors of i and n! by defining ^{(n)}(0) = (-*i*)^{n} Γ^{(n)}(0)/*n*! with Γ^{(n)}(0) representing the 1PI n-point Green's function and ø_{c} representing the classical (space-time independent) scalar field. In practice, this quantity incorporates quantum (or radiative) corrections to the classical potential which appears in the original Lagrangian density. While the latter is always finite the former can diverge due to the evaluation of momentum integrals present in the Feynman loops. One way to obtain this free energy density is to perform a functional integration over the fermionic fields [2]. To one loop this leads to

Notice that this free energy density contains the classical potential (zero loop or tree level term) present in the Lagrangian density plus a one loop quantum (radiative) correction represented by the third term. Working in the rest frame of nuclear matter we assume that the classical fields are time-like (V_{c,}_{µ} = δ_{µ}_{,0}V_{c,}_{µ}). Then, after taking the trace one can write the free energy as

where γ = 4(2) is the spin-isospin degeneracy for nuclear (neutron) matter. To obtain finite density results one may use Matsubara's imaginary time formalism with *k*_{0} → *i*(ω_{n} - *i**µ*) where m represents the chemical potential while, for fermions, ω* _{n}* = (2

*n*+1) τ

*T*(

*n*= 0,1,...) is the Matsubara frequency with

*T*representing the temperature. Then, upon using

the free energy reads

The Matsubara's sums can be performed using

where *E*^{2}(**k**) = **k**^{2}+(*M*-g_{s}ø_{c}) ^{2} and *µ*' = *µ*- *g*_{v} V_{c,0}. Being interested in the *T *= 0 case one may take the zero temperature limit of Eq. (8) which is given by^{3}

Then, at* T* = 0 and *µ* ≠ 0, the one loop free energy for the Walecka model becomes

where

Power counting shows that Δ(ø_{c}) is a divergent quantity while the *µ* dependent term of Eq. ( 10) is convergent due to the Heaviside step function.

**3. THE RENORMALIZED VACUUM CORRECTION TERM**

In order to renormalize the vacuum correction term one must first isolate the divergences which is formally achieved by regularizing the divergent integral. Here we use DR performing the divergent integrals in 2ω = 3-2ε dimensions. Then, in order to introduce the energy scale, , commonly used within QCD one redefines the integral measure as

where γ* _{E}* = -0.5772... represents the Euler-Mascheroni constant. Note that, with this definition, irrelevant factors of γ

*and 4π are automatically canceled but the results of Refs. [6, 9, 11] can be readily reproduced by using = Λ. The integral can then be performed yielding [2]*

_{E}As one can see, by expanding the the term proportional to 1/ε, there are five potentially divergent contributions ranging from *g*^{0} to *g*^{4} while all terms of order *g ^{n}* (n

__>__5) are convergent. The divergent terms proportional to Γ

^{(n)}ø

_{c}

*(n = 0,...,4) are respectively*

^{n}and

The counterterms contained in* *_{CT} needed to render the free energy finite are [6, 9]

where the α_{n} coefficients have the general form

Now, within the renormalization scheme generally adopted within QCD one simply sets *f*_{n} = 0 and the counterterms have only the bare bones needed to eliminate the 1/ε poles while the final finite contributions depend on the arbitrary energy scale. If one adopts this scheme within the Walecka model the free energy would look like the dashed curve in Fig. 1 which shows versus ø_{c} for the values^{4} = 0.9 GeV, *M* = 1 GeV, *m _{s}* = 0.55 GeV, and

*g*= 1. As it is well known, within this scheme

_{s}*M, m*, and

_{s}*m*do not represent the measurable physical vacuum masses which are instead taken as mass parameters whose values, like the values of the couplings, run with in a way ultimately dictated by RG equation.

_{v}

If instead, like Chin, one adopts the so-called on-mass renormalization scheme the counterterms completely eliminate the total contributions represented by Eqs (14-18). Within this choice the results are scale independent while *M, m _{s}*, and

*m*represent the measurable physical masses at zero density whereas the three and four-body mesonic couplings vanish in agreement with the tree level result displayed by the original Lagrangian density. The free energy generated by this scheme is represented by the dashed line in Fig. 2. Considering the relevant k

_{v}_{F}= 0 case, let us find a hybrid alternative scheme between the and the on-mass-shell so that a residual, scale dependent, contribution survives within the three and four 1PI Green's function given by Eqs (17) and (18).

To do that, let us analyze each of the arbitrary* f*_{n} terms contained in the counterterm coefficients from the physical point of view starting with *f*_{0} which is contained in the field independent Γ^{(0)}. This contribution is renormalized by the constant counterterm a_{0} which can be referred to as the "cosmological constant" [17]. In practice, the only effect this term has is to give the zero point energy value and by its complete elimination one assures that* *(ø_{c} = 0) = 0 which, within the Walecka model, will later assure that the pressure as well as the energy density vanish at *k _{F}* = 0. Therefore, as in the on-mass shell prescription, we can impose that

*f*

_{0}be exactly equal to the finite part of the Γ

^{(0)}term. It is important to point out that even if one uses the scheme this term can be absorbed in a vacuum expectation value subtraction of the zero point energy so that the exact way in which it done is not too relevant for the present purposes.

The effect of the the linear (tadpole) term Γ^{(1)} ø_{c} is to shift the origin so that the minimum is not at the origin (* _{c}* ≠ 0) as shown by the dashed line of Fig. 1. Also any finite contribution left in the tadpole will cause direct terms to contribute to the baryon self energy which, at the present level of approximation, means that M does not represent the physical nucleon mass at k

_{F}= 0,

*M**. This can be understood by recalling that the baryon self-energy is S

_{vac}_{B}~ g

_{s}[(G)\tilde]

^{(1)}( ) so that the vacuum effective baryon mass is given by

*M**= M+S

_{vac}_{B}( ) and since

*M**= 939 MeV one sees that M, as well as g

_{vac}_{s}and m

_{s}, should depend on . However, for the purposes of controlling K and

*M**the renormalization of the baryonic vacuum mass from M to

_{vac}*M**does not generate the wanted ø

_{vac}^{3}and ø

^{4}vertices. Therefore, for simplicity, we can also set

*f*

_{1}so as to completely eliminate the tadpole vacuum contribution. This choice for

*f*

_{0}and

*f*

_{1}together with

*f*

_{2}=

*f*

_{3}=

*f*

_{4}= 0 produces the dot-dashed line of figure 1. The term G

^{(2)}represents a (momentum independent) vacuum correction to the scalar meson mass, m

_{s}. As in the previous case, getting rid of this term assures that m

_{s}be taken as the physical mass simplifying the calculations since (

*m**)

_{s, vac}^{2}=

*m*

_{s}^{2}+ S

_{s}( ) where S

_{s}( ) ~ g

_{s}

^{2}[(G)\tilde]

^{(2)}( ). Fixing

*f*

_{2}so as to completely eliminate the G

^{(2)}contribution produces the dotted line of figure 1. In summary, so far we have adopted the usual Chin-Walecka on-mass shell renormalization conditions for

*f*

_{0},

*f*

_{1}, and

*f*

_{2}so that: the vacuum energy is normalized to zero, ø

_{c}= 0 is the minimum of

*(also meaning that M = M*

^{*}

_{vac}), while m

_{s}represents the vacuum scalar meson mass. In this approach, none of the vacuum mass parameters present in the original Lagrangian density run with . Note that, physically, our choice was also inspired by the NLWM observation that the compressibility modulus is improved by the introduction of ø

^{3}and ø

^{4}terms which is consistent with our choice of neglecting any corrections to terms proportional to ø

^{2}and which are respectively related with the scalar meson and baryon masses.

Now, taking *f*_{3} = 0 and *f*_{4} = 0 would leave us with the wanted ø^{3} and ø^{4} scale dependent terms. However, inspection of Eq. (17) and Eq. ( 18) shows that these contributions would vanish at different scales, given by = *Me*^{1/3} and = *Me*^{4/3} respectively. As already emphasized the NLWM controls the compression modulus with the kø^{3} and λø^{4} terms so one can impose that, within our approach, both k_{eff} and λ_{eff} arise at the same energy scale. This can be achieved by imposing that Γ^{(3)} = Γ^{(4)} = 0 at = *M* in which case the RHA is always reproduced. Finally, in our RG-NLWM inspired prescription we also impose that any dependence should be left within the leading logs which naturally emerge within DR, as shown by Eqs (14-18), and which are the main contributing terms to the β function. In this case, the α_{3} and α_{4} counterterms also eliminate the independent constants in Eqs (17) and (18). One then obtains the continuous line of Fig. 1. The complete finite, scale dependent, vacuum contribution is then given by

The free energy obtained with this finite vacuum contribution term is shown in Fig 2 for < *M* (continuous line), > *M* (dot-dashed line) as well as for = *M* (dashed line) in which case the usual RHA is retrieved. As one can check, the first term in Eq. ( 21) is just the RHA vacuum correction [6] so that, in view of Eq. ( 2), one can write

where k_{eff} = 2*g*_{s}^{3}*Mb*_{eff}( ) and λ_{eff} = 6*g _{s}*

^{4}

*c*

_{eff}( ) with

and *b*_{eff}( ) = -3*c*_{eff}( ). In this way not only k_{eff} and λ_{eff} vanish at the same scale but an inversion of their respective signs happen at the same time. We have then achieved our goal by quantically inducing k = 0 → k_{eff}( ) and λ = 0 → λ_{eff}( ) in a way that all the scale dependence is contained in the leading logs and also achieving k_{eff}( ) = λ_{eff}( ) = 0 at = M.

For comparison purposes let us quote the MRHA result

One notices that the major difference between the MRHA and our prescription amounts to the finite contribution contained within the cubic term where the scale dependence is not restricted to the leading log being also contained in an extra linear term which does not naturally arise when expands the DR results for the loop integrals in powers of e, as shown by Eqs. ( 14-18).

**4. RENORMALIZED ENERGY DENSITY**

To obtain the thermodynamical potential, Ω, one minimizes the free energy (or effective potential) with respect to the fields. That is, Ω =(_{c}, _{0}) = - *P*, where *P* represents the pressure. Then, the LHA renormalized pressure is

where E* = (**k**^{2}+*M**)^{1/2} with M* = *M*-g_{s}* _{c}*while the Fermi momentum is given by k

^{2}

_{F}= (

*µ-g*

_{v}_{0,c})

^{2}-M*

^{2}. For the vector field one gets

where ρ_{B} = (γk_{F}^{3})/(6π^{2}) is the baryonic density whereas for the scalar field the result is

where

represents the scalar density and

To get the energy density, e, one can use the relation ε = -*P*+*µ*ρ* _{B}* obtaining

where

and

where

**5. NUMERICAL RESULTS**

Let us now investigate the numerical results furnished by LHA for the baryon mass at saturation as well as for the compressibility modulus, with the latter given by

Table shows the coupling constants and saturation properties for some values of the renormalization scale ( ) that yield *BE* = -15.75 MeV and* k _{F}* = 1.42 fm

^{-1}(280.20 MeV). These values are chosen just in order to compare with the original Walecka Model (QHD-I) [1]. The meson masses are m

_{s}= 512 MeV and

*m*= 783 MeV. This table shows that some of the best LHA values are obtained with values which are very close to M. Since at = M the RHA result is reproduced one concludes, based on our results, that a slight decrease from the RHA energy scale produces an enormous effect on the values of both,

_{v}*K*and

*M**

_{sat}.

Figures 3 (a) and (b) show the binding energy per baryon, *BE = E/A -M*, and the effective baryon mass as functions of the Fermi momentum for some values, shown in table 1. One easily sees the effect of considering the vacuum contribution and its improvements on the compressibility and the effective mass. As expected, when = M, the RHA results are recovered. From figures 4 (a) and (b) it is possible to see some properties obtained in table 1 within the LHA approach, as functions of /M. One notes from figure 4 (a) that when increases the value of the nuclear compressibility (*K*) also increases and *M**_{sat} decreases. The crossing point in figure 4 (a) represents the RHA values of K and M^{*} which occurs when we set = M. Figure 4 (b) shows the effective couplings that arise due to the LHA as functions of /M. Similarly when reaches the value M the RHA results are recovered and the effective couplings vanish.

To compare our numerical results with those provided by the MRHA let us make a remark concerning the effective nucleon mass. From a non-relativistic analysis of scattering of neutron-Pb nuclei it has been found [3] that *M* _{sat}*/M ≈ 0.74 to 0.82 which can be viewed as approximately describing the Landau effective mass [4]. The relativistic isoscalar component known as the effective mass defined in Eq. (32) can be called the Dirac effective mass and is related to the Landau effective mass. Therefore, the range expected for the Dirac effective mass at saturation density lies in the range

*M**/M ≈ 0.70 to 0.80 whereas for the nuclear compressibility at saturation the most widely accepted values are K ≈ 200 MeV to 300 MeV [5]. For this range of K and according to table II the MRHA predicts

_{sat}*M**/

_{sat}*M*≈ 0.80 to 0.85 for L/M ≈ 1.185 to 1.466. However, one should note that this MRHA energy scale range is not unique and can also be reproduced with

*L/M*≈ 0.753 to 0.778 which in turn leads to a rather low range for

*M**values,

_{sat}*M**≈ 0.65 to 0.69. Our results, shown in tables I and II, seem to produce a better agreement for this range of K giving the

_{sat}/M*unique*range

*M**/

_{sat}*M*≈ 0.76 to 0.83 for

*/M*≈ 0.920 to 0.977 with k

_{eff}> 0 and λ

_{eff}< 0.

As a last remark we would like to point out that if one chooses */M* = 0.9805, (g_{s}/m_{s})^{2} = 9.468 fm^{2} and (*g _{v}/m_{v}*)

^{2}= 4.879 fm

^{2}the LHA approach reproduces the same saturation properties as performed by the so-called GM2 parameter set according to [18]: K = 300 MeV,

*M**/

_{sat}*M*= 0.78, BE = -16.3 MeV, k

_{F}= 1.313 fm

^{-1}and ρ

_{0}= 0.153 fm

^{-3}. The resulting effective couplings are:

*b*

_{eff}= 0.005986 and

*c*

_{eff}= -0.001995.

In the appendix we show that leaving a leading log dependence also in the two point Green's function, Γ^{(2)}, only increases the numerical complexity without producing results better than the ones generated by the simplest LHA version employed so far.

**6. CONCLUSIONS**

We have considered the simplest form of the Walecka model to analyze how the values of the compressibility modulus as well as the baryon mass, at saturation, can be improved by adopting an appropriate renormalization scheme in which cubic and quartic effective couplings are radiatively generated. With this aim we have evaluated the effective potential to the one loop level using Matsubara's formalism to introduce the density dependence. The vacuum contributions have been evaluated using dimensional renormalization compatible with the renormalization scheme.

We have then chosen the renormalization conditions in such a way so that all the mass parameters appearing in the original Lagrangian density represent the physical mass at zero density and therefore do not run with the energy scale. For our purposes the most important part was to renormalize the values of the cubic and quartic terms (kø^{3} / 3! and λø^{4} / 4!) which vanish at the classical (tree) level in the original model. We have then allowed only scale dependent logarithms, which naturally arise within DR, to be present in the final finite expressions and, contrary to the MRHA prescription, we obtained that both couplings have exactly the same type of scale dependence. In other words, the parameters b and c contained in the cubic and quartic terms have been dressed by one loop vacuum contributions so that *b* = 0 →* b*_{eff} = 3/π^{2} ln( /M) and* c *= 0 → - *b*_{eff}/3.

In this approach each value of the energy scale produces only one value for K and *M* _{sat}* while two values can be obtained within the MRHA. In our case the best values for these physical quantities occur at energy scales very close to the highest mass value, M. Since the RHA is obtained for = M one concludes that a small variation around this value of the energy scale can significantly improve both K and

*M**as shown by our numerical results which predict

_{sat}*M**/M ≈ 0.76 to 0.83 and K ≈ 200 MeV to 300 MeV at /M ≈ 0.920 to 0.977. These results turn to be in excellent agreement with the most quoted estimates

_{sat}*M**/M ≈ 0.70 to 0.80 and K ≈ 200 MeV to 300 MeV. To achieve these K values the MRHA predicts either

_{sat}*M**/M ≈ 0.80 to 0.85 or

_{sat}*M**/M ≈ 0.65 to 0.69 in the two possible energy scale ranges. Recalling that at /M = 1 the (RHA) results are

_{sat}*M**/M = 0.718 and

_{sat}*K =*468.14 MeV one may further appreciate how a small tuning of the energy scale within the LHA greatly improves the situation. To compare the LHA with the MRHA we recall that the philosophy within the latter is that many-body effects in nuclear matter at saturation can be minimized by choosing the energy scale close to

*M**

_{sat}in which case the values

*M**

_{sat}/M = 0.731 and

*K =*162 MeV are reproduced. Although the former seems reasonable the latter seems too low according to the above quoted estimates. The philosophy of the LHA, proposed in the present work, is to keep only the scale dependent leading logs in the finite parts of the effective cubic and quartic couplings.

In practice, the main difference between the two approximations is reflected by the fact that the MRHA effective cubic coupling, apart from the logarithmic term, also displays a term which depends linearly on the energy scale accounting for the numerical differences cited above. It is worth pointing out that, within the LHA as well as the MRHA, a given scale sets both k_{eff} and λ_{eff} so that both K and *M* _{sat}* cannot be separately tuned as in the NLWM where k and λ can be set separately. However, even in an effective theory such as the Walecka model, an increase in the parameter space as the one generated by the NLWM can be viewed as an unwanted feature and the LHA succeeds in improving the values of

*K*of

*M**without the drawback of increasing many body effects and parameter space. The method proposed here should be easy to be implemented within many existing MFA or RHA applications where Δ

_{sat}^{LHA}can be added to the energy density (in the MFA case) or used to replace the existing Δ

^{RHA}in a RHA type of calculation.

In principle the LHA philosophy could be extended to the two loop level in a calculation similar to the one performed in Refs. [7] and [9]. Then, by tuning the energy scale appropriately one could try to reduce the size of the two loop corrections producing physically meaningful results. However, at that level of approximation one may possibly encounter other issues related to the inclusion of the vacuum.

Our method is another demonstration of how phenomenological models can be affected by the choice of renormalization conditions but instead of taking this as a drawback we have shown how a judicious choice of renormalization conditions and scale can improve the value of relevant phenomenological quantities. The LHA presented here can extend MFA and RHA applications related to neutron stars as well as to the evaluation of other physical quantities, such as the symmetry energy. Then, is very plausible that MFA (and RHA) results related to quantities such as the solution of the Tolman-Oppenheimer-Volkov equations will be readily improved as we intend to demonstrate in a forthcoming work. Note also that models and/or approximations which lead to low effective masses at saturation are not suitable for neutron stars calculations since as the density increases the effective mass vanishes so quickly that higher densities cannot be properly reached as needed [19]. In principle, the LHA has potential to correct this problem without the need to introduce extra mesonic interactions with their respective parameters. This work was partially supported by CAPES and CNPq. We are grateful to Débora Menezes, Jean-Loïc Kneur and Rudnei Ramos for comments and suggestions.

**Acknowledgments**

This work was partially supported by CAPES and CNPq. We are grateful to Débora Menezes, Jean-Loïc Kneur and Rudnei Ramos for comments and suggestions.

[1] J. D. Walecka, Ann. Phys. (N.Y.) **83**, 491 (1974); [ Links ] B. D. Serot and J.D. Walecka, Adv. Nucl. Phys. **16** (1986); [ Links ]

[2] P. Ramond, *"Field Theory: a Modern Primer"* (Westview Press, 2001); [ Links ] D. Bailin and A. Love, *"Introduction to Gauge Field Theory Revised Edition"* (Taylor and Francis, 1993). [ Links ]

[3] C. H. Johnson, D. J. Horen and C. Mahaux, Phys. Rev. C **36**, 2252 (1987); [ Links ] C. Mahaux and R. Sartor, Nucl. Phys. **A475**, 247 (1987); [ Links ] Nucl. Phys. **A493**, 157 (1989). [ Links ]

[4] N. K. Glendenning, Phys. Rev. C **37**, 2733 (1988); [ Links ] M. Jaminon and C. Mahaux, Phys. Rev. C **40**, 354 (1989); [ Links ] Zhong-Yu Ma, Jian Rong, Bao-Qiu Chen, Zhi-Yuan Zhu and Hong-Qiu Song, Phys. Lett. **B604**, 170 (2004). [ Links ]

[5] J. P. Blaizot, D. Gogny and B. Grammiticos, Nucl. Phys. A**265**, 315 (1976); [ Links ] J. P. Blaizot, Phys. Rep. **64**, 171 (1980); [ Links ] H. Krivine, J. Treiner and O. Bohigas, Nucl. Phys. **A336**, 155 (1980); [ Links ] N. K. Glendenning, Phys. Rev. Lett. **57**, 1120 (1986); [ Links ] Phys. Rev. C **37**, 2733 (1988); [ Links ] M. M. Sharma, W. T. A. Borghols, S. Brandenburg, S. Crona, A. van der Woude and M. N. Harakeh, Phys. Rev. C **38**, 2562 (1988). [ Links ]

[6] S. A. Chin, Phys. Lett. **B62**, 263 (1976); [ Links ] Ann. Phys. (N.Y.) **108**, 301 (1977). [ Links ]

[7] R. J. Furnstahl, R. J. Perry and B. D. Serot, Phys. Rev. C **40**, 321 (1989); [ Links ] D. A. Wasson and R. J. Perry, Phys. Rev. C **42**, 2040 (1990). [ Links ]

[8] A Bielajew and B. Serot, Ann. Phys. (N.Y.) **156**, 215 (1984). [ Links ]

[9] D. P. Menezes, M. B. Pinto and G. Krein, Int. J. of Mod. Phys. E **9**, 221 (2000). [ Links ]

[10] J. Boguta and A. R. Bodmer, Nucl. Phys. **A292**, 413 (1977). [ Links ]

[11] E. K. Heide and S. Rudaz, Phys. Lett. **B262**, 375 (1991). [ Links ]

[12] M. Prakash, P. J. Ellis, E. K. Heidi and S. Rudaz, Nucl.Phys **A575**, 583 (1994); [ Links ] S. S. Rocha, A. R. Taurines, C. A. Z. Vaconcellos, M. B. Pinto, and M. Dillig, Mod. Phys. Lett. A **17**, 1335 (2002). [ Links ]

[13] P.A. Henning and P. Manakos, Nucl. Phys. **A** 466, 487 (1987). [ Links ]

[14] B.L. Friman and P.A. Henning, Phys. Lett.**B** 206, 579 (1988). [ Links ]

[15] P.A. Henning, Phys. Rep. **253**, 235 (1995). [ Links ]

[16] Wang Zisheng, Ma Zhongyu and Zhuo Yizhong, Phys. Rev. C **55**, 3159 (1997). [ Links ]

[17] L. Brown, *"Quantum Field Theory"* (CUP, 1994). [ Links ]

[18] N. K. Glendenning and S.A. Moszkowski, Phys. Rev. Lett. **67**, 2414 (1991). [ Links ]

[19] A. M. S. Santos and D. P. Menezes, Phys. Rev. C **69**, 045803 (2004); [ Links ] Braz. J. Phys. **34**, 833 (2004). [ Links ]

(Received on 25 August, 2010)

* Electronic address: rafael@fsc.ufsc.br

† Electronic address: marcus@fsc.ufsc.br

1 The relation between both scales is given by a constant term, =Λ, where γ* _{E}* = -0.5772....

2 See Ref. [16] for a RG investigation of the Walecka model.

3 As discussed in Ref. [15] this procedure must be taken with care if one includes loop corrections to the meson propagators which is not the case here.

4 Note that some of the these values are close to the ones which will later be used in our numerical procedure. However, at this stage they are not intended to represent any realistic physical situation apart from letting us compare possible different shapes of .

For completeness, let us check numerically the effects of leaving a leading log dependence also in the two point Green's function with zero external momentum, Γ^{(2)}, given by Eq. ( 16). Then,

and

In this case, the effective potential gives a first (momentum independent) correction to the effective scalar mass in the vacuum, m_{s,vac}^{*}. Then, for each energy scale, apart from the BE requirement one also has to fix the parameter set so that the effective scalar meson mass is *m* _{s, vac}* = 512 MeV. This effective mass is obtained by considering one loop momentum independent self energy

which clearly indicates that *m _{s}* (as well as

*g*) must run with the energy scale. However, this more cumbersome approach has almost no effect in our best results for

_{s}*K*and

*M**as table III shows indicating the adequacy of the LHA simple prescription previously adopted.

_{sat}