Many body effects in nuclear matter QCD sum rules
Abstract
We calculate the singleparticle nucleon characteristics in symmetric nuclear matter with inclusion of the and interactions. We calculated the contribution of the interactions earlier, now we add that of the ones. The contribution of the forces to nucleon self energies is expressed in terms of the nonlocal scalar condensate () and of the configurations of the vectorscalar and the scalarscalar quark condensates () in which two diquark operators act on two different nucleons of the matter.These fourquark condensates are obtained in the modelindependent way. The density dependence of the nucleon effective mass, of the vector self energy and of the singleparticle potential energy are obtained. We traced the dependence of the nucleon characteristics on the actual value of the pionnucleon sigma term. We obtained also the nucleon characteristics in terms of the quasifree nucleons, with the noninteracting nucleons surrounded by their pion clouds as the starting point. This approach leads to strict hierarchy of the many body forces.
1 Introduction
The idea of QCD sum rules in vacuum is to express the characteristics of the observed hadrons in terms of the vacuum expectation values of the QCD operators, often referred to as the condensates. Initially the method was suggested for calculation of the meson characteristics [1]. Later it was used for nucleons [2]. It is described in details in [3]. Later the approach was extended for the case of finite baryon density (see the review [4] for references). The idea of QCD sum rules in nuclear matter is to express the self energies of the probe proton in terms of the inmedium condensates.
We gave a short review of the approach in Introduction to our recent paper [5], and there is no need to repeat it here. The self energies of the probe proton are expressed in terms of the inmedium condensate. We just remind the readers that the sum rules are based on the dispersion relations for the function which carries the proton quantum numbers. The dispersion relations are considered at large and negative values of . Here the function can be calculated as a power series of with QCD condensates as the coefficients of the expansion (the Operator Product Expansion, OPE). The imaginary part of on the right hand side of the dispersion relations is formed by the real physical states with the quantum numbers of the proton. Here we separate the lowest lying pole corresponding to the probe proton. In nuclear matter we introduce the 4vector with being the vacuum value of the nucleon mass (we neglect the neutronproton mass splitting). There are three structures of the function proportional to the matrices , and to the unit matrix . The Borel transformed sum rules take the form
(1) 
Here the left hand sides are the Borel transformed OPE terms (the Borel transform is labeled by the tilde sign). They can be presented as
(2) 
Here stands for the dimension of the condensates. The right hand sides of Eq.(1) contain the effective mass of the probe proton and its vector self energy . They are multiplied by the factor with while is the residue of the nucleon pole. The position of the pole is expressed in terms of the characteristics of the probe nucleon and .
The condensates of the lowest dimension are the vector and scalar quark expectation values and Here is the ground state of nuclear matter, the sums are carried out over the quark flavors, is the nuclear matter density. The gluon condensate and the nonlocal contributions to the vector condensate have dimension . They are included in the calculations, providing the contributions which are numerically small. Thus we do not discuss them here. The contributions of dimension are provided mostly by the four quark condensates . The vector condensate , where is the number of the valence quarks in nucleon, is exactly proportional to the nucleon density. The density behavior of the other condensates depends on the approximations used for description of the nuclear matter.
It is reasonable to start the analysis treating the nuclear matter as a relativistic Fermi gas of free noninteracting nucleons. We call this a gas approximation. In such approximation each condensate can be presented as the sum of the vacuum value and nucleon density times expectation value of the same operator in the free nucleon. Neglecting the Fermi motion of the nucleons we can write
(3) 
with the nucleon matrix element for the nucleon at rest
(4) 
It can be expressed through the observable pionnucleon sigma term [6]
(5) 
with the masses of light quarks. The sigma term can be expressed as
(6) 
Here is the pion mass while is the energy of free nucleon at rest.
For the nucleon with three momentum one can find the corrections of the order to the right hand side of Eq.(4), replacing by on the right hand side of Eq.(6). After integration over the Fermi sphere of the radius one finds [7]
(7) 
Recall that . If the kinematical corrections are neglected, i.e. the scalar condensate is expressed by Eq.(3) and is linear in .
The gas approximation corresponds to interaction of the probe proton with each of the nucleons of the matter separately. Hence Eq.(1) provides the nucleon characteristics with inclusion of the forces only [4].
Contribution of the forces to the self energies of the probe proton was found in [5] by inclusion of the inmedium QCD condensates beyond the gas approximation and taking into account only the forces between the nucleons of the matter. Neglecting the kinematical corrections we can write
(8) 
with the energy of interactions between the nucleons of the matter with inclusion of only the forces. Define similar to Eq.(6)
(9) 
and
(10) 
The scalar condensate is thus
(11) 
Here we neglected corrections caused by the nucleon kinetic energy. Note that the latter are not as simple as in the gas approximation (Eq.(7)). Indeed, the nucleons of the matter can be viewed as moving in certain vector and scalar fields and with the energy . Thus expansion in powers of should be rather that in powers of . At saturation value of nuclear density [8], and thus can not be neglected. Note also that depends on .
There is also a contribution to the forces acting on the probe proton, originated by configurations of the fourquark condensates in which two pairs of quarks act on two different nucleons of the matter contribute to the forces.
Now we want to calculate the proton self energies with inclusion of the forces. The nucleon energy is
(12) 
with the contribution of the forces to the energy of interactions between the nucleons of the matter. Introducing similar to Eq.(9)
(13) 
we define
The quark condensate becomes
(14) 
One can view as the value of the operator averaged over the nucleon bound in nuclear matter with the and forces included.
In the meson exchange picture the nucleons of the matter exchange by the mesons . The contributions with should be summed over the mesons . The contribution of each meson to is proportional to the expectation value which is (under the proper normalization) just the total number of quarks and antiquarks in the meson. All the mesons, but the pions contain mainly the valence quarks. Thus the simple quark counting provides . The pions contain the sea of the quarkantiquark pairs, and one can calculate . Thus we expect the pion exchanges to give the leading contribution to . Only these contributions will be included.
The configurations in which two pairs of quark operators act on two different free nucleons contribute to the forces between the nucleons of the matter. Thus they contribute to the forces acting on the probe proton. The contributions of the vector scalar () and of the scalarscalar () fourquark condensates which can be obtained in the modelindependent way [5] depend on the nucleon matrix element . Including the lowest inmedium corrections to we find the contribution of the forces to the fourquark condensates leading to the interactions between the probe nucleon an nuclear matter. The condensate is proportional to . To find the contribution of the interactions to the parameters of the probe nucleon one should change by in the expression for the condensate [5]. In the same way the condensate containers two factors . To obtain the contribution of the interactions one should change one of the matrix elements by multiplying the result by the factor reflecting the permutations.
We can establish connection between our approach and the picture based on the nucleonnucleon interactions. The lowest dimension OPE terms are determined by exchange of weakly correlated quarks with nuclear matter. These terms determine the probe proton self energies and . On the other hand these parameters can be considered as caused by the exchanges of the systems of strongly correlated quarks (the mesons) between the probe nucleon and nuclear matter. Thus the exchange by systems of strongly correlated quarks between the probe nucleon and the nucleons of the matter is expressed through those of weakly correlated quarks with the same quantum numbers between the three quark system and the matter. The Feynman diagrams of nucleonnucleon interactions corresponding to inclusion of and forces were presented in [5]. Additional diagrams for the forces are given in Fig.1.
To analyze the role of various contributions it is convenient to present
(15) 
These equations follow immediately from Eq.(1). The two equations become independent if we put with the vacuum value of the continuum threshold. Under this assumption Eqs.(15) express the nucleon parameters in terms of condensates. In the total solution the unknown ties the equations.
The contributions of the forces to the vector self energy manifest themselves in terms of the fourquark condensates with two pairs of quark operators acting on two different nucleons of the matter located in the same space time point. The probability of such configuration is small due to strong repulsion of the two nucleons, known as the ”core” in traditional nuclear physics. It is known to be due mostly to the Pauli principle for the quarks in the sixquark system [9]. The result was supported by the lattice QCD stimulations [10], [11]. In our calculations the form of the current insures the needed antisymmetrization and leads to numerically small value of the contribution. Still smaller contribution comes from forces since is much smaller than the linear contribution to the scalar condensate (see, e.g. [4]). Thus the many body contributions to exhibit the standard hierarchy with domination of the gas approximation term while the forces provide a small correction and inclusion of forces lead to a still smaller correction.
The contribution of the and forces to the effective mass manifest themselves mainly through the nonlinear scalar condensates (NSC) and correspondingly. These condensates have been obtained in [12] in framework of the Chiral Perturbation Theory (CHPT). Still the NSC is much smaller than the linear part of the condensate and the contribution of and forces is much smaller than that of forces. However it follows from the equations presented in [12] that due to some cancelations the contribution is several times smaller than (). Thus the contribution of forces to the effective mass appears to be larger than that of forces.
The contributions of the and forces to the vector self energy manifest themselves in terms of the fourquark condensates. Due to small value of the contribution of the forces to is much smaller than that of the forces.
We employ the results of [12] for NSC with minor modifications which we discuss in the next Section. We found the density dependence of nucleon parameters with inclusion of the many body effects. The main effect of many body forces is the change of the value of nucleon effective mass . At ( is the phenomenological value of saturation density) and MeV the forces diminish by MeV. The forces increase by MeV. Recall that the forces diminish by 37 MeV. The forces change by several MeV.
Note that the calculations of [12] include the contributions provided by both the direct (Hartree) and the exchange (Fock) interactions between the nucleons of the matter. We do not include the exchanges between our threequark system described by the function and the matter since such terms correspond to the higher order terms of the operator expansion (see, e.g. [4]). Thus in our approach the exchange terms of interaction between the probe nucleon as a whole and the matter are expected to be small and are not included.
The linear part of the condensate and the NSC, as well as of the scalarscalar and vectorscalar fourquark condensates depend on the value of nucleon sigma term . While the conventional value is MeV [6], one can meet the larger values of in literature (see, e.g. [13]). The value MeV [14] is the largest one. We traced dependence of nuclear parameters on the value of .
Note that definition of many body effects is ambiguous. We treated the single pion exchange as a twonucleon interaction. On the other hand it can be expressed as the nucleon self energy in the matter, which can be treated as a single particle characteristics of a quasifree nucleon. The same refers to iterated one pion exchange [15]. The nucleon energy with inclusion of and iterated exchanges can be written as with corresponding to inclusion of and interactions. In the gas approximation for the quasifree nucleons we can write similar to Eq.(6)
(16) 
Beyond the gas approximation and with inclusion of the and forces we can write for the energy of nucleon interaction in the quasifree nucleon presentation
(17) 
and define
(18) 
We can solve the sum rules in the quasifree nucleon presentation employing the equations presented above in which and are replaced by and correspondingly.
In Sec.2 we calculate the contribution of the nonlinear scalar condensate. In Sec.2 we find the condensates contributing to the forces. We obtain the density dependence of nucleon parameters with inclusion of the and forces in Sec.3. In Sec. 4 we solve the sum rules in the quasifree nucleon presentation. We summarize in Sec.5.
2 Condensates contributing beyond the gas approximation
2.1 Nonlinear scalar condensate
The scalar condensate contributes only to equation for in Eq.(1). It provides [4]
(19) 
Note that contributions of the higher real states to the right hand side of the dispersion relation are approximated by continuum with the effective threshold and are accounted for by the function with . We do not write down the lowest order radiative corrections. However they are included in the calculations.
To calculate the condensate we need to know the contributions and expressed through the terms and by Eqs.(10) and (13). The contribution was obtained in [5]. Now we obtain .
It includes the iterated onepion exchange terms involving three inmedium nucleons, expressed by Eqs.(8) and (9) of [12]. We must also include threebody contributions to the twopion exchange terms expressed by Eqs.(12) and (13) of [12]. There are specific for the CHPT vertices proportional to . The threebody terms proportional to are given by Eqs.(22) and (23) of [12]. These terms provide the contribution
(20) 
to the effective sigma term. Here is the sum of the right hand sides of Eqs.(8), (9), (12), (13),(22) and (23) of [12].
There is additional contribution. Note that in the lowest order of CHPT
with standing for the energy of the inmedium nucleon. On the other hand, the effective nucleon sigma term
Thus in the lowest order of CHPT . The higher order CHPT terms are included by Eq.(19) of [12] which can be written as
(21) 
with the high momentum cutoff while is proportional to and does not depend on . The explicit form of the term can be easily deduced from Eq.(19) of [12].
2.2 Fourquark condensates
Going beyond the gas approximation we must include the contribution to the four quark condensate in which two pairs of quarks act on two different nucleons of the matter, which are at the same spacetime point. As we said in Introduction, the strong short range repulsion of the two nucleons is due mostly to the Pauli principle for the quarks in the sixquark system [9][11]. In the present calculations, as well as it was in those carried out in [5] the form of the function insures the needed antisymmetrization, and the latter effectively takes into account the main part of the nucleon short range interaction.
In [5] the contribution to was expressed in terms of the expectation values in the free nucleons. The contribution of the vectorscalar condensate was thus proportional to expectation value determined by Eq.(4). The contribution of the scalarscalar condensate was proportional to . To include the forces in the matter we must replace by in the vectorscalar condensate and by in the scalar condensate. The term would correspond to the forces in nuclear matter and should be omitted.
The vectorscalar condensate is determined by the expectation value of operator. Employing Eq.(36) of [5] we find that it contributes only to the scalar structure with
(24) 
The expectation value of operator provides the scalarscalar condensate contributing to . Employing Eq.(33) of [5] we find
(25) 
3 Contribution of and forces to the values of nucleon parameters
3.1 Contribution of forces
Contribution of forces was studied in [5]. Now we study what do the forces add.
Putting in Eq.(15) we see that the forces caused by the nonlinear scalar condensate (NSC) and the scalarscalar () four quark condensate contribute only to the effective mass. Since the vectorscalar () four quark condensate contributes to the structure , the corresponding forces influence both the effective mass and the vector self energy . In the total solution the unknown ties the equations. Variations of the NSC and SS condensates provide also minor changes in the value of .
Now we solve Eqs. (1) including the forces. The terms corresponding to and forces are given in [4] and [5]. We add the contributions corresponding the forces described in previous Section. The contribution of forces to the effective mass is determined by the contribution to the NSC. We found that MeV at while MeV at this point [5]. Thus the contribution of the forces to the effective mass is several times larger that of the forces. Note that is dominated by the Hartree term of the iterated onepion exchange (Eq.(8) of [14]). The contribution determined by Eq.(23), which makes MeV at is not very important numerically but is needed for the consistency of the procedure. Thus at the forces increase the nucleon effective mass by MeV in agreement with estimations made in [5].
The density dependence of nuclear parameters with and without inclusion of the forces is presented in Figs.2,3. More detailed results in the vicinity of the density are given in Table 1.
,MeV  , MeV  , MeV  ,GeV  , GeV  
0.90  640  196  93  1.43  1.75 
618  166  –145  1.44  1.79  
675  173  80  1.62  1.88  
0.95  620  209  100  1.39  1.73 
596  176  157  1.41  1.79  
661  183  85  1.61  1.88  
1.00  599  222  108  1.36  1.72 
574  185  169  1.39  1.79  
647  192  89  1.59  1.88  
1.05  577  236  115  1.32  1.71 
551  195  183  1.37  1.79  
632  202  93  1.58  1.88  
1.10  555  250  123  1.29  1.70 
528  204  196  1.36  1.80  
618  212  98  1.57  1.88  

The contribution of the forces can be seen by comparing the solid and dashed lines in Figs.2,3 and of the numbers in the third and second lines of Table 1.
As one can see from Table 1, the forces also increase the vector self energy by MeV. Note, however that the result for is determined by the shift of the continuum threshold value in the selfconsistent solution of both equations composing Eq.(15). Indeed putting we find that changes only due to the change of the condensate given by Eq.(25). Direct estimation shows that this makes less than MeV. Thus the result for the shift of depends on the model chosen for the description of the higher states in QCD sum rules. For example, even staying in framework of the ”pole+ continuum” model, we would find another value for , taking different threshold values for different channels. Thus our physical result is that at the forces increase the nucleon effective mass by MeV, while their contribution to is consistent with zero.
These are the results for the conventional value of the cutoff MeV. Since the CHPT parameter and the value determined by Eq.(22) depend on , the contribution of forces also depends on . At MeV the interactions increase the value of the effective mass by MeV at . At MeV they increase by MeV. The value of changes by about MeV in these cases.
Note that our main aim was to calculate the nucleon self energies. We do not claim that our approach is precise enough to reproduce the values of the potential energy with since it is obtained as a result of subtraction of two large positive values and .
In Table 2 we demonstrate how the values of the nucleon self energies at change after separate inclusion of the nonlinear scalar condensate (NSC) and of the fourquark condensates. Note that the separate contributions are not additive, mainly because the self energies are connected with the condensates by the nonlinear expressions given by Eq.(15).
.
,MeV  , MeV  

NSC  74  8 
VS  2  0.1 
SS  3  1.2 
3.2 Total contribution of many body forces
Now we compare the values of nucleon parameters with inclusion of the and forces with those of the gas approximation. The two cases are shown by the solid and dotted lines in Figs.2,3 and in the third and second lines of Table 1.
We find that the many body interactions increase the value of the effective mass by 48 MeV and decrease that of the vector self energy by 30 MeV. Note that there is an uncertainty of about 15 MeV in the value of due to possible variation of the high momentum cutoff . It is caused mainly by the contribution to the NSC of two pion exchanges with two nucleons in the intermediate state. This uncertainty is about MeV for the vector self energy .
We trace also the dependence of nucleon parameters on the value of . The results are shown in Table 3.
.
, MeV  ,MeV  , MeV  MeV 

35  746  210  27 
40  697  201  31 
45  647  192  89 
50  597  185  147 
55  548  178  202 
60  500  171  257 
The sharp dependence of nucleon effective mass on the value of manifests itself in the gas approximation since the scalar condensate proportional to determines the main part of the shift . The mass changes by MeV while changes from MeV to MeV. The many body forces make the dependence somewhat stronger mainly due to the contribution proportional to . These two contributions exhibit linear dependence on . The and four quark condensates also depend on . In the latter case the dependence is not linear. However these contributions are numerically less important.
We obtained reasonable values for the nucleon potential energy for MeV. At smaller values of the potential energy runs positive. At larger values of the value of is too large.
4 Many body effects in the system of quasifree nucleons
.
Until now we treated the single pion exchange as a twonucleon interaction. However another approach is possible. The physical nucleon in vacuum includes the self energy contribution caused by radiation and absorption of the pion. The nucleon propagator included the sum over all possible intermediate states. For the inmedium nucleon one should exclude the occupied nucleon states with momenta smaller than the Fermi momentum see Fig.4a. These are just the Fock terms of the nucleon interaction with onepion exchange which can be treated as a part of the interaction. On the other hand this contribution can be expressed as the nucleon self energy in the matter, which can be viewed as a single particle characteristics of a quasifree nucleon.
The same refers to iterated one pion exchange [15]. In the lowest order they can be expressed in terms of the fourth order selfenergy diagrams for the inmedium nucleon see Fig. 4b. We define the sigma term for quasifree nucleon by Eq.(16). We present
(26) 
with the sum of the right hand sides of Eqs.(4), (6), (7), (8) and (9) of [12]. Also the matrix element of the operator for the quasifree nucleon is
We find MeV at . This provides MeV which is much smaller than in the gas approximation for free nucleons.
Now we include the interactions between the quasifree nucleons of the matter. They determine the interactions between the probe nucleon and the matter in the quasifree nucleon presentation. We calculate the contribution to the NSC
(27) 
with the sum of the right hand sides of Eqs.(11), (14), (17), (18) and (20) of [12]. We must include the contribution of the fourquark condensates found in [5] with changed to in Eqs. (33) and (36). We find MeV at . The main contribution comes from twopion exchanges with virtual excitations.
To include the interactions between the quasifree nucleons of the matter which determine the interactions between the probe nucleon and the matter in the quasifree nucleon presentation, we find the corresponding NSC
(28) 
Here is the sum of the right hand sides of Eqs.(12), (13), (22) and (23) of [12], while is determined by Eq. (23) with replaced by . We obtain MeV at . The leading contributions are determined by the twopion exchanges with virtual excitations and by the term .
Since we obtain the strict hierarchy of contributions of the many body forces to the vector self energy , to the effective mass and to the potential energy .
Density dependence of the nucleon characteristics with inclusion of the , and interactions between the probe nucleon and the matter in the quasifree nucleon presentation is shown in Fig.5. The density dependence of the residue and of the continuum threshold are shown in Fig.6. More detailed data in the vicinity of the point are presented in Table 4.
Since the role of the fourquark condensates increases in the quasifree nucleons presentation. The contribution of the forces to the vector self energy making 9 MeV at is mainly due to the SS fourquark condensate.
,MeV  , MeV  , MeV  ,GeV  , GeV  
0.90  509  182  237  1.14  1.60 
646  167  116  1.55  1.83  
664  173  92  1.57  1.85  
0.95  472  195  262  1.10  1.60 
627  175  127  1.53  1.83  
648  182  98  1.54  1.85  
1.00  433  207  288  1.07  1.59 
607  183  138  1.52  1.83  
631  192  105  1.52  1.84  
1.05  392  221  316  1.04  1.59 
587  191  150  1.50  1.83  
613  202  112  1.50  1.84  
1.10  347  235  347  1.02  1.59 
567  199  162  1.49  1.83  
595  212  129  1.48  1.85  

Note that the results for the total contribution of the and interactions are close to those obtained in the previous Section but are not identical to them. The main reason is that now we replace by in expression for . This changes the contribution of the right hand side of Eq.(20) of [12] to the function (Eq.(26)) and of those of the right hand sides of Eqs.(22,23) to the function . We replaced also by in Eq.(22) for . Thus the values of at differ by 16 MeV.
5 Summary
We calculated the contribution of many body forces to the nucleon characteristics in symmetric nuclear matter. We found the contribution of the interactions earlier [5]. Now we add that of the forces.
The main contribution of the forces to the effective mass is determined by the nonlinear scalar condensate with inclusion the interactions between the nucleons of the matter. The contribution to the vector self energy of the nucleon is dominated by the fourquark condensates in the scalarscalar and vectorscalar channels with two diquark operators acting on two different nucleons of the matter. The values of these condensates are determined by the nonlinear scalar condensate with inclusion of the interactions between the nucleons of the matter. Connection between our approach and the picture based on the nucleonnucleon interactions is illustrated by Fig.1.
The results of calculations at MeV are presented in Figs.2,3 and in Table 1. We found that at the forces increase the effective mass by 73 MeV. Earlier we obtained that inclusion of forces diminishes the value of by 26 MeV [5]. Thus the total contribution of the and forces increase the gas approximation value by MeV. The contribution of forces is larger than that of the forces. The forces change the value of the vector self energy by several MeV while the forces diminished it by MeV. Thus the contribution of the forces to the vector self energy is much smaller than that of the forces.
We traced the dependence of the nucleon parameters on the actual value of . We found reasonable values for the nucleon potential energy for MeV. At larger values of the bound of the nucleons in the matter looks to be too strong.
Calculations of the nucleon energy in the chiral effective field theory reviewed in [16] demonstrated that the forces give smaller contribution than the forces. Our results can not be directly compared with those of [16]. In [16] both Hartree and Fock terms in the scalar channel have been obtained. Our analysis corresponds rather to the mean field approach. The Fiertz transformed Fock terms of [16] contribute to both scalar and vector channels in our approach. In calculations carried out in framework of traditional nuclear physics the authors of [17] found that the contribution of the forces to the binding energy of the nucleon in nuclear matter is smaller than that of the forces. On the other hand, a number of other calculations showed the and correlations to be of the same order [18].
We developed also another approach in which the onepion exchanges and the iterated onepion exchanges which can be expressed in terms of the nucleon self energy are treated as the singlenucleon characteristics. In this approach the gas approximation corresponds to system of noninteracting nucleons with their pion clouds (quasifree nucleons). Parameters of the probe nucleon calculated in the gas approximation for the quasifree nucleons corresponds to inclusion of the forces. Density dependence of the nucleon characteristics corresponding to inclusion of the , and interactions between the probe nucleon and the nucleons of the matter in the quasifree presentation is given in Fig.5,6 and in Table 4. The many body interactions between the probe nucleon and the nucleons of the matter respect the strict hierarchy. The contribution of interactions to the effective mass , to the vector selfenergy and to the potential energy are much smaller than those of interactions. The contributions of forces are much smaller than those of the forces.
References
 [1] M.A. Shifman, A.I. Vainshtein and V.I. Zakharov, Nucl. Phys. B 147 (1979) 385; 448; 519.
 [2] B.L. Ioffe, Nucl. Phys. B 188 (1981) 317; B 191 (1981) 591(E).
 [3] B.L. Ioffe, L.N. Lipatov and V.S. Fadin, Quantum Chromodynamics, Cambridge Univ. Press, Cambridge, 2010.
 [4] E. G. Drukarev, M. G. Ryskin and V. A. Sadovnikova, Phys. Atomic Nuclei 75 (2012) 334. (arXiv:nuclth/1012.0394).
 [5] E. G. Drukarev, M. G. Ryskin and V. A. Sadovnikova, Nucl. Phys. A 959 (2017) 129.
 [6] J. Gasser, H. Leutwyler and M. E. Sainio, Phys. Lett. B 253 (1991) 252.
 [7] N. Kaiser, S. Fritsch and W. Weise, Nucl. Phys. A 697 (2002) 255.
 [8] L. S. Celenza and C. M. Shakin, Relativistic Nuclear Physics, World Scientific, Philadelphia 1988.
 [9] V. G. Neudatchin, Yu. F. Smirnov and R. Tamagaki, Prog. Theor. Phys. 58 (1977) 1072.
 [10] N. Ishii, S. Aoki and T. Hatsuda, Phys. Rev. Lett. 99 (2007) 022001.
 [11] S. Aoki, T. Doi, T. Hatsuda, Y. Ikeda, T. Inoue, N. Ishii, K. Murano, H. Nemuda, K. Sasaki, Prog. Theor. Exper. Phys. 2012 (2012) 01A105.
 [12] N. Kaiser, P. de Homont and W. Weise, Phys.Rev. C 77 (2008) 025204 .
 [13] L. AlvarezRuso, T. Ledwig, J. M. Camalich and M. J. VicenteVacas, Phys. Rev. D 88 (2013) 054507.
 [14] G. E. Hite, W. B. Kaufmann, and R. J. Jacob, Phys. Rev. C 71 (2005) 065201.
 [15] N. Kaiser, R. Brockmann and W. Weise, Nucl. Phys. A 625 (1997) 758.
 [16] E. Epelbaum, H. W. Hammer, U. G. Meissner, Rev. Mod. Phys. 81 (2009) 1773.
 [17] D. W. E. Blatt and B. H. J. McKellar, Phys. Rev. C 11 (1975) 2040.
 [18] H. Bethe, Ann. Rev. Nucl. Sci. 21 (1971) 93.
Appendix A Figure captions
Fig. 1. Feynman diagrams of nucleonnucleon interactions corresponding to inclusion of the forces. Solid lines are for the probe nucleon, bold lines denote the nucleons of the matter. The shaded block denotes the pion field. The dashed line denotes the scalar meson. The dasheddotted line denotes the scalar or vector meson.
Fig. 2. Density dependence of the effective nucleon mass , of the vector self energy and of the singleparticle potential energy . The horizontal axis corresponds to the density related to its saturation value . The vertical axis is for , and . The dotted linesonly the forces between the nucleons of the matter are included. Dashed linesthe and forces are included. Solid linesthe , and forces are included.
Fig. 3. Density dependence of the nucleon residue and of the continuum threshold .
The notations are the same as in Fig. 2.
Fig. 4. The Feynman self energy diagrams for the nucleon in the matter, corresponding to the onepion exchange (Fig.a) and to the iterated onepion exchange (Fig.b). Solid lines are for nucleons, wavy lines are for the pions. The cross on the line marks that the nucleon is on the mass shell.
Fig. 5. Density dependence of the effective nucleon mass , of the vector self energy and of the singleparticle potential energy
in the quasufree nucleons presentation.
The notations are the same as in Fig.2.
Fig. 6 Density dependence of the nucleon residue and of the continuum threshold in the quasifree nucleons presentation.
The notations are the same as in Fig. 2.