Print version ISSN 0103-9733
Braz. J. Phys. vol.37 no.1a São Paulo Mar. 2007
Brazilian relativistic O(q4) two-pion exchange nucleon-nucleon potential: parametrized version
Carlos A. da RochaI; Renato HigaII; Manoel R. RobilottaIII
INúcleo de Pesquisa em Bioengenharia, Universidade São Judas Tadeu, Rua Taquari, 546, 03166-000, São Paulo, São Paulo, Brazil. Email address email@example.com
IIThomas Jefferson National Accelerator Facility, 12000 Jefferson Avenue, Newport News, VA 23606, USA. Email address: firstname.lastname@example.org
IIIInstituto de Física, Universidade de São Paulo, Caixa Postal 66318, 05315-970, São Paulo, São Paulo, Brazil. Email address: email@example.com
In our recent works we derived a chiral O(q4) two-pion exchange nucleon-nucleon potential (TPEP) formulated in a relativistic baryon (RB) framework, expressed in terms of the so called low energy constants (LECs) and functions representing covariant loop integrations. In order to facilitate the use of the potential in nuclear applications, we present a parametrized version of our configuration space TPEP.
Keywords: Two-pion exchange; Nucleon-nucleon potential; Relativistic baryon; Chiral symmetry; Chiral perturbation theory
The one-pion exchange NN potential (OPEP) is simple, has been well established long ago, and dominates completely partial waves with orbital momentum L > 5. The two-pion exchange potential (TPEP), on the other hand, is rather complex and has become free of important ambiguities only in the 1990s, after the systematic use of chiral symmetry in its theoretical description [1-11].
Chiral perturbation theory is based on the existence of a characteristic scale q, set by both pion four-momenta and nucleon three-momenta, such that q < 1 GeV. Due to this technique, nowadays one understands rather well the internal hierarchies of the NN potential in terms of chiral layers. Leading terms of the chiral TPEP are of order O(q2) and expansions which go up to O(q4) are already available. One of them was produced recently by our group [10, 11]. We departed from a relativistic Lagrangian and evaluated the relevant Feynman diagrams covariantly, without resorting to heavy baryon approximations. The so obtained T-matrix was then transformed into a potential, expressed in terms of covariant loop integrals and observable parameters. Without loss of generality, one may choose these parameters to be either the subthreshold coefficients extracted from pN scattering  or the low-energy constants (LECs) present in the effective Lagrangian.
The research on the TPEP performed in the last decade has set its conceptual foundations on a rather solid basis, comparable to that of the OPEP in the late sixties. On the other hand, the TPEP depends on several LECs, which must be extracted from either pN scattering data or direct fits of NN phase shifts. In general, this last kind of procedure tends to be computationally heavy, for theoretical results are usually given as cumbersome expressions. In order to make applications easier, in this work we present a parametrized version of our O(q4) relativistic configuration space TPEP, which is numerically accurate for distances larger than 1 fm. It is based on the theoretical expressions derived in Ref. .
II. THEORETICAL POTENTIAL
The O(q4) relativistic expansion of the TPEP produced in refs.[10, 11] was based on the evaluation of three families of diagrams given in Fig. 1. The first of them involves only pion and nucleon degrees of freedom into single loops and corresponds to the minimal realization of chiral symmetry. It includes the subtraction of the iterated OPEP and yields the terms in the profile functions given below which are proportional to just or 1/. Terms proportional to 1/, on the other hand, come from two-loop processes, either in the form of t-channel contributions from the second family or s and u-channel terms embodied in the subthreshold coefficients of the third family. Finally, the third group of diagrams includes chiral corrections associated with other degrees of freedom, hidden within the LECs ci and di, and gives rise to contributions which are proportional to either (LEC)/ or (LEC)2/.
The configuration space potential has the isospin structure
and WLS = L·(s(1)+s(2))/2 , WT = 3s(1)· s(2)· - s(1)· s(2), WSS = s(1)· s(2). The form of the operator WQ in configuration space is highly non-local and can be found in ref..
The radial components of the potential are expressed in terms of the following dimensionless profile functions with I = C, LS, T, SS:
where t+ = 3 and t- = 2.
Therefore, the configuration space potential is written in terms of numerical coefficients which multiply dimensionless functions arising form the Fourier transforms of Feynman loop integrals. The former are combinations of external parameters representing the pion and nucleon masses, µ and m, respectively, the pion decay constant fp, the axial coupling constant gA, and the LECs ci and di. The latter, denoted by , depend on just µ, m and x º m r.
III. PARAMETRIZED POTENTIAL
We keep the external quantities as free and parametrize the dimensionless Feynman loop integral functions and their derivatives in terms of Zi º (Gi, Hi, Ii). They have the following structure
which is a combination of an exponential function times a polynomial with coefficients gi. This parametrization is more than 1% accurate in the range 0.8 fm < r < 10 fm.
Using the definition a º m/m, the profile functions are written as
The fitted functions Zi º (Gi, Hi, Ii) given by eq. (7) can be re-written as
where º = ågi xn is the adjusted polynomial, which is given by
For instance, the parametrized function I8 is given by:
Nowadays, the asymptotic expressions for the chiral TPEP have the status of theorems and are written as sums of chiral layers, with little model dependence. These asymptotic functions were our guide to guess what were the correct choice of the parametrized functions to fit the potential profile functions.
The parametrized profile functions given in the preceding section depend explicitly on four well known quantities, namely m, µ, gA, fp, and on the less known LECs ci and di. Therefore, the latter may be extracted from fits to data. When doing this, however, one has to bear in mind that, as discussed in Ref. , the influence of the LECs over the profile functions is rather uneven. Indeed, their influence over , and is rather strong, but barely perceptible in and .
We conclude saying that the various channels of the potential are clearly dominated by either nucleonic (minimal model) or non-nucleonic (LECs) degrees of freedom, mainly because we found the role of intermediate pp scattering (two-loops) to be small.
The work of R. H. was supported by DOE Contract No. DE-AC05-06OR23177, under which JSA operates the Thomas Jefferson National Accelerator Facility, and M. R. R. was supported by FAPESP (Fundação de Amparo à Pesquisa do Estado de São Paulo).
 L.S. Celenza, A. Pantziris, and C.M. Shakin, Phys. Rev. C 46, 2213 (1992); [ Links ]J.L. Friar and S.A. Coon, Phys. Rev. C 49, 1272 (1994); [ Links ]M.C. Birse, Phys. Rev. C 49, 2212 (1994). [ Links ]
 J-L. Ballot, M. R. Robilotta, and C. A. da Rocha, Phys. Rev. C 57, 1574 (1998). [ Links ]
 N. Kaiser, R. Brockman, and W. Weise, Nucl. Phys. A 625, 758 (1997). [ Links ]
 M. C. M. Rentmeester, R. G. E. Timmermans, J. L. Friar, and J. J. de Swart, Phys. Rev. Lett. 82, 4992 (1999); [ Links ]M. C. M. Rentmeester, R. G. E. Timmermans, and J. J. de Swart, Phys. Rev. C 67, 044001 (2003). [ Links ]
 D. R. Entem and R. Machleidt, Phys. Rev. C 66, 014002 (2002). [ Links ]
 R. Higa and M. R. Robilotta, Phys. Rev. C 68, 024004 (2003). [ Links ]
 R. Higa, M. R. Robilotta, and C. A. da Rocha, Phys. Rev. C 69, 034009 (2004).
 G. Höhler, group I, vol. 9, subvol.b, part 2 of Landölt-Bornstein Numerical data and Functional Relationships in Science and Technology, ed. H. Schopper, 1983; [ Links ]G. Höhler, H. P. Jacob and R. Strauss, Nucl. Phys. B 39, 273 (1972). [ Links ]
 N. Hoshizaki and S. Machida, Progr. Theor. Phys. 24, 1325 (1960). [ Links ]
Received on 29 september, 2006