Following closely the methodology of introducing higher-order
derivatives of zero-range pseudopotentials [11], we obtain
the general form of the finite-range pseudopotential by coupling
derivative operators with spin operators in the spherical-tensor
formalism. The complete list of the building blocks of the tensors
can be found in Ref. [9]. We obtain the
finite-range version by simply substituting the zero-range Dirac
delta with a regularized (smeared) delta,
In such a formalism, the general expression for the antisymmetrized finite-range
pseudopotential reads,
For the sake of completeness, we recall here definitions of
tensors that appear in Eq. (3). Each term of the
pseudopotential is a scalar with a composite internal structure and is
built of higher-rank spherical tensors in the space, spin, and isospin
coordinates. Higher-order tensor derivatives
are built of
relative momenta
, with spherical
components
In the spin space, for particles numbered by 1 or 2, we have as building
blocks the spin unity matrix of rank 0, denoted as
, and
the Pauli matrices
, which
define the spin spherical tensor of rank 1,
The locality delta functions together with the regularized delta must be understood as distributions, whereupon the derivatives in and act. The locality deltas ensure the locality of the energy density only in the zero-range case, whereas it is in principle nonlocal due to the presence of the so called Majorana operator in the exchange part of the pseudopotential.
In Eq. (3), index , taking either value 0 or 1, defines the dependence of the pseudopotential on the isospin. This dependence is in addition to that stemming from the exchange term, and is a new element as compared to the case of the zero-range pseudopotential [11].
In principle, one could write down an equivalent version of the antisymmetrized finite-range pseudopotential, where the new dependence on the isospin exchange operator would have been substituted by the Majorana operator exchanging space coordinates, . This would be possible in virtue of the relation , valid when acting on a fermionic wave function. However, the presence of the Majorana operator in the pseudopotential would break the correspondence between the exchange part of the antisymmetrized pseudopotential and the nonlocal sector of the EDF, according to the common understanding of the nonlocal part of the EDF as resulting from the exchange part of the antisymmetrized interaction. While in a zero-range pseudopotential this operator reduces to a phase depending on the order of the relative-momenta tensors [23,11], in a finite-range interaction this operator acts explicitly, and switches the space coordinates of the two interacting nucleons. Therefore the finite-range character of the interaction, complemented with the action of the Majorana operator, would have given rise to the nonlocal form of the functional stemming from the direct interaction. Although such a construction would be perfectly correct and exactly equivalent to that presented below, it would also be counterintuitive, and thus in our study it is not further pursued.
The constructed pseudopotential possesses all symmetries of the nuclear interaction. Referring the reader to Ref. [11] for a comprehensive discussion, here we only list these symmetries, which are rotational symmetry, time-reversal, parity, and Galilean invariance, supplemented with the hermiticity of the pseudopotential operator and indistinguishability principle operating as invariance under exchange of particles 1 and 2. The pseudopotential (2) is indeed invariant under all the transformations corresponding to the symmetries listed above. This guarantees that the resulting EDF is also invariant with respect to the same symmetries.
All terms of the pseudopotential, derived at zero, second, fourth, and sixth order, , 2, 4, and 6, respectively, are collected in the supplemental material. The supplemental material contains all detailed results of the present study. They were derived by means of symbolic programing, and are presented in readable format as well as in the form directly usable in computer programming. In Table 1 we show numbers of different terms of pseudopotential (2), where we also distinguish between central (), SO (), and tensor () terms, corresponding to different ranks in the coupling of the relative momenta tensors with the spin operator . At each order, the overall numbers of terms equal 4, 14, 30, and 52, giving the total number of 100 terms up to NLO. At each order, these numbers are twice larger than the corresponding numbers for the zero-range pseudopotential, which reflects the addition of the new quantum number .
We note, that here we classified tensor terms as those corresponding to ; however, in the pseudopotential (3) there also appear terms that have purely space-tensor character, e.g., those for and . However, we know from the generalized Cayley-Hamilton (GCH) theorem [24], that any scalar function of two vectors and must have form of a function of three elementary scalars , , and . Therefore, the central term () can always be recoupled to the form where terms with do not appear. An explicit construction of such a form is presented in the Cartesian coordinates in Sec. 4.
Order | Total | ||||
0 | 4 | 0 | 0 | 4 | |
2 | 8 | 2 | 4 | 14 | |
4 | 16 | 4 | 10 | 30 | |
6 | 24 | 8 | 20 | 52 | |
NLO | 52 | 14 | 34 | 100 |
Jacek Dobaczewski 2014-12-07