Print version ISSN 0103-9733
Braz. J. Phys. vol.40 no.4 São Paulo Dec. 2010
IDepartamento 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
IIDepartamento de Física, Universidade Federal de Santa Catarina, 88040-900 Florianópolis, Santa Catarina, Brazil
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.
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  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 .
The important parameters are the ratios of coupling to masses, C2s and C2v, with Ci2 = (giM/mi)2 which are tuned to fit the saturation density ρ0 = 0.193 fm-3 and binding energy per nucleon, BE = -15.75 MeV . 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 , 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 . 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  (see Ref.  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 . 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. .
Meanwhile, the compressibility modulus problem has been circumvented by introducing some more parameters, in the form of new couplings , to the original Walecka model. One then considers U(ø,V) appearing in Eq. (1) as
where k = 2bMgs3 and λ = 6cgs4. 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  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 related1 ). 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. . 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.  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  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  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 keffø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 keff 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 → keff( ) and λ = 0 → λeff( ) we can keep M,ms,mv as representing the vacuum physical masses for simplicity. This choice means that, at kF = 0, all mass parameters (M,ms,mv) represent the effective vacuum masses and shall not run with as opposed to keff 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 work2. 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 keffø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 ms 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 
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 . 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 (Vc,µ = δµ,0Vc,µ). 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 k0 → i(ωn - iµ) where m represents the chemical potential while, for fermions, ωn = (2n+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 E2(k) = k2+(M-gsøc) 2 and µ' = µ- gv Vc,0. Being interested in the T = 0 case one may take the zero temperature limit of Eq. (8) which is given by3
Then, at T = 0 and µ ≠ 0, the one loop free energy for the Walecka model becomes
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 γE 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 
As one can see, by expanding the the term proportional to 1/ε, there are five potentially divergent contributions ranging from g0 to g4 while all terms of order gn (n > 5) are convergent. The divergent terms proportional to Γ(n)øcn (n = 0,...,4) are respectively
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 fn = 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 values4 = 0.9 GeV, M = 1 GeV, ms = 0.55 GeV, and gs = 1. As it is well known, within this scheme M, ms, and mv 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.
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, ms, and mv 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 kF = 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 fn terms contained in the counterterm coefficients from the physical point of view starting with f0 which is contained in the field independent Γ(0). This contribution is renormalized by the constant counterterm a0 which can be referred to as the "cosmological constant" . 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 kF = 0. Therefore, as in the on-mass shell prescription, we can impose that f0 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 kF = 0, M*vac. This can be understood by recalling that the baryon self-energy is SB ~ gs[(G)\tilde](1)( ) so that the vacuum effective baryon mass is given by M*vac = M+SB( ) and since M*vac = 939 MeV one sees that M, as well as gs and ms, should depend on . However, for the purposes of controlling K and M*vac the renormalization of the baryonic vacuum mass from M to M*vac does not generate the wanted ø3 and ø4 vertices. Therefore, for simplicity, we can also set f1 so as to completely eliminate the tadpole vacuum contribution. This choice for f0 and f1 together with f2 = f3 = f4 = 0 produces the dot-dashed line of figure 1. The term G(2) represents a (momentum independent) vacuum correction to the scalar meson mass, ms. As in the previous case, getting rid of this term assures that ms be taken as the physical mass simplifying the calculations since (m*s, vac)2 = ms2+ Ss( ) where Ss( ) ~ gs2 [(G)\tilde](2)( ). Fixing f2 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 f0,f1, and f2 so that: the vacuum energy is normalized to zero, øc = 0 is the minimum of (also meaning that M = M*vac), while ms 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 f3 = 0 and f4 = 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 = Me1/3 and = Me4/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 keff 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  so that, in view of Eq. ( 2), one can write
where keff = 2gs3Mbeff( ) and λeff = 6gs4ceff( ) with
and beff( ) = -3ceff( ). In this way not only keff 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 → keff( ) and λ = 0 → λeff( ) in a way that all the scale dependence is contained in the leading logs and also achieving keff( ) = λ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* = (k2+M*)1/2 with M* = M-gscwhile the Fermi momentum is given by k2F= (µ-gv0,c)2-M*2. For the vector field one gets
where ρB = (γkF3)/(6π2) is the baryonic density whereas for the scalar field the result is
represents the scalar density and
To get the energy density, e, one can use the relation ε = -P+µρB obtaining
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 kF = 1.42 fm-1 (280.20 MeV). These values are chosen just in order to compare with the original Walecka Model (QHD-I) . The meson masses are ms = 512 MeV and mv = 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, 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  that M*sat/M ≈ 0.74 to 0.82 which can be viewed as approximately describing the Landau effective mass . 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*sat/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 . For this range of K and according to table II the MRHA predicts 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*sat values, M*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 unique range M*sat/M ≈ 0.76 to 0.83 for /M ≈ 0.920 to 0.977 with keff > 0 and λeff < 0.
As a last remark we would like to point out that if one chooses /M = 0.9805, (gs/ms)2 = 9.468 fm2 and (gv/mv)2 = 4.879 fm2 the LHA approach reproduces the same saturation properties as performed by the so-called GM2 parameter set according to : K = 300 MeV, M*sat/M = 0.78, BE = -16.3 MeV, kF = 1.313 fm-1 and ρ0 = 0.153 fm-3. The resulting effective couplings are: beff = 0.005986 and ceff = -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.
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 → beff = 3/π2 ln( /M) and c = 0 → - beff/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*sat as shown by our numerical results which predict M*sat/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 M*sat/M ≈ 0.70 to 0.80 and K ≈ 200 MeV to 300 MeV. To achieve these K values the MRHA predicts either M*sat/M ≈ 0.80 to 0.85 or M*sat/M ≈ 0.65 to 0.69 in the two possible energy scale ranges. Recalling that at /M = 1 the (RHA) results are M*sat/M = 0.718 and 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 keff 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*sat 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 Δ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.  and . 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 . 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.
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.
 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 ]
 J. P. Blaizot, D. Gogny and B. Grammiticos, Nucl. Phys. A265, 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 ]
 A Bielajew and B. Serot, Ann. Phys. (N.Y.) 156, 215 (1984). [ Links ]
 D. P. Menezes, M. B. Pinto and G. Krein, Int. J. of Mod. Phys. E 9, 221 (2000). [ Links ]
 J. Boguta and A. R. Bodmer, Nucl. Phys. A292, 413 (1977). [ Links ]
 E. K. Heide and S. Rudaz, Phys. Lett. B262, 375 (1991). [ Links ]
 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 ]
 P.A. Henning and P. Manakos, Nucl. Phys. A 466, 487 (1987). [ Links ]
 B.L. Friman and P.A. Henning, Phys. Lett.B 206, 579 (1988). [ Links ]
 P.A. Henning, Phys. Rep. 253, 235 (1995). [ Links ]
 Wang Zisheng, Ma Zhongyu and Zhuo Yizhong, Phys. Rev. C 55, 3159 (1997). [ Links ]
 L. Brown, "Quantum Field Theory" (CUP, 1994). [ Links ]
 N. K. Glendenning and S.A. Moszkowski, Phys. Rev. Lett. 67, 2414 (1991). [ Links ]
(Received on 25 August, 2010)
* Electronic address: firstname.lastname@example.org
Electronic address: email@example.com
1 The relation between both scales is given by a constant term, =Λ, where γE = -0.5772....
2 See Ref.  for a RG investigation of the Walecka model.
3 As discussed in Ref.  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,
In this case, the effective potential gives a first (momentum independent) correction to the effective scalar mass in the vacuum, ms,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 ms (as well as gs) must run with the energy scale. However, this more cumbersome approach has almost no effect in our best results for K and M*sat as table III shows indicating the adequacy of the LHA simple prescription previously adopted.