# The partonic structure of the nucleon from generalized transverse momentum-dependent parton distributions

###### Abstract

We discuss the general formalism for the calculation in light-front quark models of the fully unintegrated, off-diagonal quark-quark correlator of the nucleon, parametrized in terms of generalized transverse momentum dependent parton distributions (GTMDs). By taking specific limits or projections, these GTMDs yield various transverse-momentum dependent and generalized parton distributions, thus providing a unified framework to simultaneously model different observables. The corresponding distributions in impact-parameter space are the Wigner functions which provide multidimensional images of the quark distributions in phase space. We present results within a light-front constituent quark model, discussing some of the complementary aspects encoded in the different distributions and the relation to the quark orbital angular momentum of the proton.

## 1 Introduction

Parton distributions entering many hard and exclusive processes play a key role to describe the nonperturbative structure of hadrons. The most complete information is contained in the generalized transverse momentum dependent parton distributions (GTMDs) Meissner:2009ww (); Meissner:2008ay (); Lorce:2011dv () which parametrize the unintegrated off-diagonal quark-quark correlator, depending on the quark longitudinal and transverse momentum, and , respectively, and on the 4-momentum which is transferred by the probe to the hadron; for a classification see Refs. Meissner:2009ww (); Meissner:2008ay (). The GTMDs give the full one-quark density matrix in the momentum space and reduce to different parton distributions and form factors as is shown in Fig. 1.

The different arrows in this
figure represent particular projections in the hadron and quark momentum space, and give
the links between the matrix elements of different reduced density matrices.
Such matrix elements can in turn be parametrized in terms of generalized parton distributions
(GPDs), transverse-momentum dependent parton distributions (TMDs) and
form factors (FFs). These are the quantities which enter the description of various
exclusive (GPDs), semi-inclusive (TMDs), and inclusive (PDFs) deep inelastic scattering
processes, or parameterize elastic scattering processes (FFs). At leading twist, there are
sixteen complex GTMDs, which are defined in terms of the independent polarization states
of quarks and hadron. In the forward limit they reduce to eight TMDs which depend
on the longitudinal momentum fraction and transverse momentum of quarks,
and therefore give access to the three-dimensional picture of the hadrons in momentum
space. On the other hand, the integration over of the GTMDs leads to eight GPDs
which are probability amplitudes related to the off-diagonal matrix elements of the parton
density matrix in the longitudinal momentum space. The common limit of TMDs and GPDs is
given by the standard parton distribution functions (PDFs), related to the diagonal matrix
elements of the longitudinal-momentum density matrix for different polarization states of
quarks and hadron. The integration over of the GTMDs leads to a bilocal operator restricted to the
plane transverse to the light-front direction and brings to the lower plane of the box in Fig. 1.
The off-forward matrix elements of this operator can be parametrized in terms of so-called
transverse-momentum dependent form factors (TMFFs). Starting from the TMFFs, we
can follow the same path as in the case of the GTMDs, and at each vertex of the basis of
the box of Fig. 1 we find the restricted version of the operator defining the distributions
in the upper plane. Therefore, integrating out the dependence on the quark transverse
momentum, we encounter matrix elements parametrized in terms of form factors (FFs),
while the forward limit of TMFFs leads to transverse-momentum dependent spin densities
(TMSD). Both FFs and TMSDs have the charges as common limit.

After
appropriate Fourier transform, the GTMDs can be interpreted as Wigner or phase-space distributions
Ji:2003ak (); Belitsky:2003nz (); Belitsky:2005qn (); Lorce:2011kd (); Lorce:2011ni (),
giving access to the correlations between quark momentum and transverse position.
The Wigner distributions reduce to the Fourier transform of the GPDs in impact-parameter space (or impact-parameter dependent distributions) after integration over the quark transverse momentum, and, upon further integration over the longitudinal quark momentum,
to the charge densities in the transverse coordinate plane.

Although a variety of models has been employed to explore separately the different
observables related to GTMDs, a unifying formalism for modeling the GTMDs has been presented only recently Lorce:2011dv ().
In the following, we will review some of the results discussed in Ref. Lorce:2011dv (), using the language of light-front wave functions
(LFWFs) and focusing on the three-quark (3Q) contribution.
In Sect. 2 we present the formal derivation of
the LFWF overlap representation of the quark contribution to GTMDs, specializing the
results to two light-front quark models, namely the light-front chiral quark-soliton model (LFQSM) and
the light-front constituent quark model (LFCQM).
In Sect. 3 we introduce the Wigner distributions, discussing the case of unpolarized quarks in the longitudinally polarized nucleon
and its relation to the quark OAM.
Then, in Sects. 4 and 5, we discuss
some of the complementary aspects encoded in the GPDs and the TMDs, in particular with regards to the information on the quark OAM. Concluding remarks are given in the final section.

## 2 Quark-quark Correlator

The fully-unintegrated quark-quark correlator for a spin- hadron is defined as Meissner:2009ww (); Meissner:2008ay ()

(1) |

This correlator is a function of the initial and final hadron light-front helicities and , the average hadron and quark four-momenta and , and the four-momentum transfer to the hadron (see Fig. 2 for the kinematics).

The superscript stands for any element of the basis in Dirac space. A Wilson line ensures the color gauge invariance of the correlator, connecting the points and via the intermediary points and by straight lines. This induces a dependence of the Wilson line on the light-front direction . Since any rescaled four-vector with some positive parameter could be used to specify the Wilson line, the correlator actually only depends on the four-vector , where is the hadron mass. The parameter gives the sign of the zeroth component of , i.e. indicates whether the Wilson line is future-pointing () or past-pointing ().

Since the parton light-front energy is particularly difficult to access in high-energy experiments, the relevant correlators are actually obtained from the integrated version of Eq. (1), setting all the fields at the same light-front time :

(2) |

where we used for a generic four-vector the light-front components and the transverse components , and is the average fraction of longitudinal momentum carried by the quark. A complete parametrization of this object in terms of GTMDs has been achieved in Meissner:2009ww (); Meissner:2008ay ().

### 2.1 Overlap Representation

Following
the lines of Diehl:2000xz (); Brodsky:2000xy (), we obtain in the light-front gauge an overlap representation for the correlator (2) at the twist-two level restricted to the 3Q Fock sector^{1}^{1}1Quark flavor and color indices have been omitted for clarity. In the processes considered here the flavor and color of a given quark remain unchanged.

(3) |

where the integration measures are defined as

(4) |

Furthermore, in Eq. (3) the function selects the active quark average momentum (we choose to label the active quark with and the spectator quarks with ). The 3Q LFWF depends on the momentum coordinates of the quarks relative to the hadron momentum (collectively indicated by ), and the index which stands for the set of the quark light-front helicities . The transition from the initial quark light-front helicity to the final one is described by a complex-valued matrix . In particular, we have for the spectator quarks . For the active quark, the matrix depends on the twist-two Dirac structure used in the correlator, see e.g. Boffi:2002yy (); Boffi:2003yj (); Pasquini:2005dk (). We choose to work in an infinite momentum frame such that is large, and . The four-momenta involved are then

(5) | ||||||

Note that the form used for is not the most general one, but leads to an appropriate definition of TMDs for semi-inclusive deep inelastic and Drell-Yan processes. For the active and spectator quarks the initial and final momentum coordinates are then

(6) | ||||||

So far, the exact 3Q LFWF derived directly from the QCD Lagrangian is not known. Nevertheless, we can try to reproduce the gross features of hadron structure
at low scales using constituent quark models. Many models exist on the market based on the concept of constituent quarks. However only a few incorporate consistently relativistic effects. We focus here on two such models, the light-front constituent quark model (LFCQM) Boffi:2002yy (); Boffi:2003yj (); Pasquini:2005dk () and the light-front chiral quark-soliton model (LFQSM) Petrov:2002jr (); Diakonov:2005ib (); Lorce:2007as (); Lorce:2007fa ().
However, the formalism can be easily generalized to other quark models as explained in Refs. Lorce:2011zta (); Pasquini:2012jm ().

The LFWFs used in LFCQM and in QSM have a very similar structure, given by

(7) |

where is a global symmetric momentum wave function, is the spin-flavor wave function, and is an matrix connecting light-front helicity and canonical spin

(8) |

The explicit expressions for the momentum wave function in Eq. (7) and the vector in Eq. (8) in LFCQM read

(9) |

where is a normalization factor, is the free invariant mass, is the free energy of quark , is the constituent quark mass, and are model parameters fitted to reproduce the anomalous magnetic moments of the nucleon Pasquini:2007iz (). On the other hand, within the QSM one has

(10) |

where is the soliton mass, is the energy of the discrete level in the spectrum, and are the upper and lower components of the Dirac spinor describing this discrete level.

For further convenience we introduce the tensor correlator

(11) |

where and with the Pauli matrices. We now use the LFWF given by Eq. (7) and write the overlap representation of the correlator tensor as

(12) |

where stands for

(13) |

In Eq. (13), and the matrix is given by

(14) |

The tensor correlator in Eq. (12) has two indices. The index refers to the transition in terms of hadron light-front helicity, while the index refers to the transition in terms of the active quark light-front helicity. For example, the components and correspond to the matrix elements of the and operators in the case of an unpolarized hadron, respectively. Equation (12) gives the explicit expression for the tensor correlator in terms of the overlap of initial and final symmetric (instant-form) momentum wave functions with the tensor for a fixed average momentum of the active quark . The tensor contains the spin-flavor structure derived from the overlap of the three initial and final quarks. Taking into account the possible couplings of the helicities of the active and spectator quarks to give the hadron helicity, the coefficient and in Eq. (11) for SU(6) spin-flavor wave functions are

(15) |

Furthermore, the matrix in Eq. (12) describes the overlap of the initial and final quark state. The columns are labeled by the index which indicates the type of transition in terms of quark light-front helicity. The rows are labeled by the index which indicates the type of transition in terms of quark canonical spin. This matrix reduces to for the spectator quarks, since in this case the light-front helicity is conserved.

## 3 Wigner distributions

By performing a Fourier transform of the GTMDs to the impact-parameter space, we obtain quark distributions which are naturally interpreted as Wigner distributions Ji:2003ak (); Belitsky:2003nz (); Lorce:2011kd ()

(16) |

Although the GTMDs are in general complex-valued functions, their two-dimensional Fourier transforms are always real-valued functions, in accordance with their interpretation as phase-space distributions. We note that, like in the usual quantum-mechanical Wigner distributions, and are not Fourier conjugate variables. However, they are subjected to Heisenberg’s uncertainty principle because the corresponding quantum-mechanical operators do not commute . As a consequence, the Wigner functions can not have a strict probabilistic interpretation. There are in total 16 Wigner functions at twist-two level, corresponding to all the 16 possible configurations of nucleon and quark polarizations. Here we will discuss only one particular case, namely the distortion in the distribution of unpolarized quarks due to the longitudinal polarization of the nucleon which has a close connection with the quark orbital angular momentum (OAM). Other configurations for the quark and nucleon polarizations can be found in Ref. Lorce:2011kd () .

In Fig. 4, the upper panels show the distortions in impact-parameter space for (left panel) and (right panel) quarks with fixed transverse momentum and GeV. We observe a clear dipole structure in these distributions which indicates that the (quasi-)probability for finding the quark orbiting clockwise is not the same as for the quark orbiting anticlockwise, leading in average to a nonvanishing OAM. The lower panels of Fig. 4 describes the distribution in impact-parameter space of the quark average transverse momentum in a longitudinally polarized nucleon

(17) |

We observe that the average transverse momentum is always orthogonal to the impact-parameter vector . This is not surprising since a nonvanishing radial component of the average transverse momentum would indicate that the proton size and/or shape are changing. We can also clearly notice that quarks tend to orbit anticlockwise inside the nucleon, corresponding to positive OAM aligned with the nucleon spin which is pointing out of the figure. For the quarks, we see two regions. In the central region of the nucleon, fm, the quarks tend to orbit anticlockwise like the quarks, while in the peripheral region, fm, the quarks tend to orbit clockwise, with a flip of the local net quark OAM. Note that such information about the OAM can not be accessed through GPDs and TMDs since none of them describe at leading twist the distortion in the distribution of unpolarized quarks due to the longitudinal polarization of the nucleon. This is because one needs the correlation between and which is lost by integrating over or .

The Wigner distributions were originally constructed as the quantum mechanical analogue of the classical density operator in the phase space. In particular, any matrix element of a quark operator can be rewritten as a phase-space integral of the corresponding classical quantity weighted by the Wigner distribution. It is therefore natural to define the quark OAM as follows Lorce:2011kd ()

(18) |

Since the Wigner distribution involves in its definition a gauge link, it inherits a path dependence Lorce:2012ce (); Lorce:2012rr () . The simplest choice is a straight gauge link. In this case, Eq. (18) gives the kinetic OAM associated with the quark OAM operator appearing in the Ji decomposition Ji:1996ek (); Ji:2012sj () , where is the usual covariant derivative. According to the Ji’s relation Ji:1996ek (), this kinetic quark OAM can be extracted from the GPDs

(19) |

with

(20) | |||||

(21) |

In order to connect the Wigner distributions to the TMDs, it is more natural to consider instead a staple-like gauge link consisting of two longitudinal straight lines connected at by a transverse straight line. In this case, Eq. (18) gives the canonical OAM associated with the quark OAM operator that appears in the Jaffe-Manohar decomposition in the gauge Jaffe:1989jz (); Lorce:2011ni (); Hatta:2011ku (), i.e. .

## 4 GPDs in impact-parameter space

In this section, we discuss a few examples of spin densities parametrized in terms of GPDs in impact-parameter space. As outlined in the introduction, they can be obtained from the Wigner distributions after integration over the quark transverse momentum and can be interpreted as probability densities of quarks with longitudinal momentum fraction and transverse location with respect to the nucleon center of momentum Burkardt:2005td (). In Fig. 4 we show the results within the LFCQM Boffi:2007yc (); Pasquini:2007xz () in the case of unpolarized quarks in a transversely polarized nucleon. This spin density is given by the sum of a nucleon spin-independent contribution related to the GPD and a nucleon spin-dependent contribution from the GPD , corresponding to monopole and dipole distributions in impact-parameter space, respectively. The dipole contribution introduces a large distortion perpendicular to both the nucleon spin and the momentum of the proton, with opposite sign for and quarks. Such a distortion reflects the large value of the anomalous magnetic moments .

With the present model, and , to be compared with the values and derived from data. This effect can serve as a dynamical explanation of a non-vanishing Sivers function which measures the correlation between the intrinsic quark transverse momentum and the transverse nucleon spin Burkardt:2002ks (). This connection between the GPD and the Sivers function has recently been exploited in Ref. Bacchetta:2011gx () to determine the total quark angular momentum from the Ji’s relation (20), reconstructing the GPD in the collinear limit from the available experimental information on in SIDIS Hermes05a (); Alekseev:2008aa (). Though this estimate is based on a model-dependent relation, the consistency with constraints on the angular momentum arising from DVCS measurements Airapetian:2008aa (); Mazouz:2007aa () is encouraging.

## 5 TMDs in momentum space

The eight leading-twist TMDs are a
natural extension of standard parton distribution from one to three dimensions in
momentum space, being function of both the longitudinal quark momentum fraction and the transverse momentum .
The knowledge of TMDs allow us to build tomographic images of the inner structure of the nucleon
in momentum space, complementary to the impact-parameter space tomography that can be achieved by studying
GPDs.

The LFWF overlap representation of the TMDs has been explicitly derived in Refs. Pasquini:2008ax (); Pasquini:2010af (); Brodsky:2010vs () and can be also obtained using the results of Sect. 2 for the GTMDs in the forward limit .
This representation is
well suited to illustrate the relevance of the different orbital angular momentum components
of the nucleon wave function, and provide an intuitive picture for the physical meaning of
the quark TMDs. Moreover, they can be regarded as initial input for phenomenological
studies for the semi-inclusive processes where quark TMDs play a very important role Boffi:2009sh (); Pasquini:2011tk ().
Most of the TMDs would simply vanish
in absence of quark OAM.
Recently, it has been suggested, on the basis of some
quark-model calculations She:2009jq (); Avakian:2010br (), that the TMD
may be related to the quark OAM:

(22) |

However, Eq. (22) is not a rigorous expression and holds only in a specific class of quark models. For a detailed discussion on the the physical origin of this relation and the underlying model assumptions for its validity we refer to Lorce:2011kn (). The TMD describes the distortion due to the transverse polarizations in perpendicular directions of the quark and the nucleon Miller:2003sa (). In this case, the nucleon helicity flips in the direction opposite to the quark helicity, with a mismatch of two units for the orbital angular momentum of the initial and final LFWFs. The corresponding quadrupole structure in the momentum space for both and quarks is shown in Fig. 5, as obtained from the model of Ref. Pasquini:2008ax ().

Finally, in Table 1, we summarize the results from the LFCQM and the LFQSM for the quark OAM obtained from the Ji relation (Eq. (19)), the Wigner distributions (Eq. (18)) and the TMD (Eq. (22)).

Model | LFCQM | LFQSM | ||||
---|---|---|---|---|---|---|

Total | Total | |||||

As expected in a pure quark model, all the definitions give the same value for the total quark OAM, with nearly twice more net quark OAM in the LFCQM than in the LFQSM. The difference between the various definitions appears in the separate quark-flavour contributions. Note in particular that unlike the LFCQM, the LFQSM predicts a negative sign for the -quark OAM in agreement with lattice calculations Hagler:2007xi (). It is surprising that since it is generally believed that the Jaffe-Manohar and Ji’s OAM should coincide in absence of gauge degrees of freedom. Note that a similar observation has also been made in the instant-form version of the QSM Wakamatsu:2005vk (). On the other hand, the individual quark contributions to the OAM obtained from Eq. (22) do not correspond to the intrinsic quark orbital angular momentum, and therefore do not coincide with the results for . The two calculations agree only for the total OAM, since in the sum over the individual quark contributions the spurious terms due to the transverse centre of momentum cancel out.

## 6 Conclusions

In this work we presented a study of GTMDs, which parametrize the fully-unintegrated quark-quark correlators with the quark fields are taken at the same light-front time. By taking specific limits or projections of these GTMDs, they yield PDFs, TMDs, GPDs, FFs, and charges, accessible in various inclusive, semi-inclusive, exclusive, and elastic scattering processes. The GTMDs therefore provide a unified framework to simultaneously model these different observables. We discussed a first step in this modeling, by considering a light-front wave function (LFWF) overlap representation of the GTMDs and by restricting ourselves to the 3Q Fock components in the nucleon LFWF. At twist-two level, we studied the most general transition which the active quark light-front helicity can undergo in a polarized nucleon, corresponding to the general helicity amplitudes of the quark-nucleon system. We develop a formalism which is quite general and can be applied to many quark models as long as the nucleon state can be represented in terms of 3Q without mutual interactions. By Fourier transform in the transverse space of the GTMDs we obtain the Wigner distributions which provide the multidimensional images of the quark distributions in the phase space. In particular, we discussed results for the Wigner distributions of unpolarized quarks in a longitudinally polarized nucleon that allow us to calculate the phase-space average of the quark OAM. Other ways to access information about the quark OAM from GPDs and TMDs have been also discussed, comparing the corresponding results obtained within a light-front constituent quark model and the light-front chiral quark-soliton model.

## Acknowledgments

This work was supported in part by the Research Infrastructure Integrating Activity ÒStudy of Strongly Interacting MatterÓ (acronym HadronPhysic3, Grant Agreement n. 283286) under the Seventh Framework Programme of the European Community, by the Italian MIUR through the PRIN 2008EKLACK ÒStructure of the nucleon: transverse momentum, transverse spin and orbital angular momentumÓ, and by the P2I (ÒPhysique des deux InfinisÓ) network.

## References

- (1) S. Meissner, A. Metz and M. Schlegel, JHEP 0908, (2009) 056.
- (2) S. Meissner, A. Metz, M. Schlegel and K. Goeke, JHEP 0808, (2008) 038.
- (3) C. Lorcé, B. Pasquini and M. Vanderhaeghen, JHEP 1105, (2011) 041.
- (4) X. d. Ji, Phys. Rev. Lett. 91, (2003) 062001.
- (5) A. V. Belitsky, X. d. Ji and F. Yuan, Phys. Rev. D 69, (2004) 074014.
- (6) A. V. Belitsky and A. V. Radyushkin, Phys. Rept. 418, (2005) 1.
- (7) C. Lorcé and B. Pasquini, Phys. Rev. D 84, (2011) 014015.
- (8) C. Lorcé, B. Pasquini, X. Xiong and F. Yuan, Phys. Rev. D 85, (2012) 114006.
- (9) M. Diehl, T. Feldmann, R. Jakob and P. Kroll, Nucl. Phys. B 596, (2001) 33 [Erratum-ibid. B 605, (2001) 647].
- (10) S. J. Brodsky, M. Diehl and D. S. Hwang, Nucl. Phys. B 596, (2001) 99.
- (11) S. Boffi, B. Pasquini and M. Traini, Nucl. Phys. B 649, (2003) 243.
- (12) S. Boffi, B. Pasquini and M. Traini, Nucl. Phys. B 680, (2004) 147.
- (13) B. Pasquini, M. Pincetti and S. Boffi, Phys. Rev. D 72, (2005) 094029; Phys. Rev. D 76, (2007) 034020.
- (14) V. Y. Petrov and M. V. Polyakov, arXiv:hep-ph/0307077.
- (15) D. Diakonov and V. Petrov, Phys. Rev. D 72, (2005) 074009.
- (16) C. Lorcé, Phys. Rev. D 78, (2008) 034001.
- (17) C. Lorcé, Phys. Rev. D 79, (2009) 074027.
- (18) C. Lorcé and B. Pasquini, Phys. Rev. D 84, (2011) 034039.
- (19) B. Pasquini and C. Lorcé, arXiv:1203.5006 [hep-ph].
- (20) B. Pasquini and S. Boffi, Phys. Rev. D 76, (2007) 074011.
- (21) C. Lorcé, arXiv:1210.2581 [hep-ph].
- (22) C. Lorcé, arXiv:1205.6483 [hep-ph].
- (23) X. D. Ji, Phys. Rev. Lett. 78, (1997) 610.
- (24) X. Ji, X. Xiong and F. Yuan, Phys. Rev. Lett. 109, (2012) 152005.
- (25) R. L. Jaffe and A. Manohar, Nucl. Phys. B 337, (1990) 509.
- (26) Y. Hatta, Phys. Lett. B 708, (2012) 186
- (27) M. Burkardt, Int. J. Mod. Phys. A 21, (2006) 926; Int. J. Mod. Phys. A 18, (2003) 173; Phys. Rev. D 62, (2000) 071503 [Erratum-ibid. 66, (2002) 119903].
- (28) S. Boffi and B. Pasquini, Riv. Nuovo Cim. 30, (2007) 387.
- (29) B. Pasquini and S. Boffi, Phys. Lett. B 653, (2007) 23.
- (30) M. Burkardt, Phys. Rev. D 66, (2002) 114005; M. Burkardt and D. S. Hwang, Phys. Rev. D 69, (2004) 074032.
- (31) A. Bacchetta and M. Radici, Phys. Rev. Lett. 107, (2011) 212001.
- (32) A. Airapetian et al. (Hermes Collaboration), Phys. Rev. Lett. 94, (2005) 012002.
- (33) M. Alekseev et al. [COMPASS Collaboration], Phys. Lett. B 673, (2009) 127.
- (34) A. Airapetian et al. [HERMES Collaboration], JHEP 0806, (2008) 066.
- (35) M. Mazouz et al. [Jefferson Lab Hall A Collaboration], Phys. Rev. Lett. 99, (2007) 242501.
- (36) B. Pasquini, S. Cazzaniga and S. Boffi, Phys. Rev. D 78, (2008) 034025.
- (37) B. Pasquini and F. Yuan, Phys. Rev. D 81, (2010) 114013.
- (38) S. J. Brodsky, B. Pasquini, B. -W. Xiao and F. Yuan, Phys. Lett. B 687 (2010) 327.
- (39) S. Boffi, A. V. Efremov, B. Pasquini and P. Schweitzer, Phys. Rev. D 79, (2009) 094012.
- (40) B. Pasquini and P. Schweitzer, Phys. Rev. D 83, (2011) 114044.
- (41) J. She, J. Zhu and B. -Q. Ma, Phys. Rev. D 79, (2009) 054008.
- (42) H. Avakian, A. V. Efremov, P. Schweitzer and F. Yuan, Phys. Rev. D 81, (2010) 074035.
- (43) C. Lorcé and B. Pasquini, Phys. Lett. B 710, (2012) 486
- (44) G. A. Miller, Phys. Rev. C 68, (2003) 022201.
- (45) Ph. Hägler et al. [LHPC Collaborations], Phys. Rev. D 77, (2008) 094502.
- (46) M. Wakamatsu, H. Tsujimoto, Phys. Rev. D 71, (2005) 074001; M. Wakamatsu, Eur. Phys. J. A 44, (2010) 297.