Nuclear electric dipole moment of three-body system

Nuclear electric dipole moment of three-body system


Nuclear electric dipole moments of and are calculated using Time Reversal Invariance Violating (TRIV) potentials based on the meson exchange theory, as well as the ones derived by using pionless and pionful effective field theories, with nuclear wave functions obtained by solving Faddeev equations in configuration space for the complete Hamiltonians comprising both TRIV and realistic strong interactions. The obtained results are compared with the previous calculations of EDM and with time reversal invariance violating effects in neutron-deuteron scattering.

24.80.+y, 11.30.Er, 21.10.Ky

I Introduction

The electric dipole moment (EDM) of particles is a very important parameter in searching for Time Reversal Invariance Violation (TRIV) and for the possible manifestation of new physics. The discovery of non-zero value of the EDM would be a clear evidence of TRIV Landau (1957), therefore, it has been a subject for intense experimental and theoretical investigations for more than 50 years. The search for TRIV also has fundamental importance for the explanation of the baryon asymmetry Beringer et al. (2012) of the Universe which requires a source of CP-violation Sakharov (1967) beyond that entering the Cabibbo-Kobayashi- Maskawa matrix of the Standard Model. Any observation of EDM in the near future will be a direct indication of new physics beyond the Standard Model. However, theoretical estimates for the values of particle EDMs are extremely small which results in many difficulties in experimental search for neutron and electron EDMs. Therefore, it is desirable to consider more complex systems, where EDMs or another TRIV parameters could be enhanced. This also would provide assurance that there would be enough observations to avoid a possible “accidental” cancelation of T-violating effects due to unknown structural factors related to strong interactions. The study of the EDM is particularly important for the simplest few-nucleon systems, the result of which may lead to a better understanding of TRIV effects in heavier nuclei. Moreover, few-nucleon systems meet requirements for a number of proposals to measure EDMs of light nuclei in storage rings Khriplovich (1998); Farley et al. (2004); Semertzidis (2009); Lehrach et al. (2012).

In this paper, we calculate nuclear EDMs of and using TRIV potential in meson exchange model, as well as in pionless and pionful EFT. Weak TRIV potentials are used in conjunction with realistic strong interaction Hamiltonians. Several realistic nucleon-nucleon potentials have been tested to represent the strong interaction: the Argonne v18(AV18), the Reid soft core(Reid93), the Nijmegen(NijmII), the INOY, as well as the AV18 in conjunction with the three-nucleon Urbanna IX(UIX) potential. Three-nucleon wave functions have been obtained by solving Faddeev equations in the configuration space for the complete Hamiltonians, comprising both TRIV and strong interactions.

Ii Time reversal violating potentials

The most general form of time reversal violating and parity violating part of nucleon-nucleon Hamiltonian in the first order of relative nucleon momentum can be written as the sum of momentum independent and momentum dependent parts, Herczeg (1966),


where the exact form of depends on the details of the particular theory. Because of the additional factor, the non-static potential contributions are suppressed by a factor , therefore, we consider here only static TRIV interactions which could be obtained within three different approaches: in a meson exchange model, pionless EFT, and pionful EFT.

The TRIV meson exchange potential in general involves exchanges of pions (, MeV), -mesons(, MeV), and - and -mesons (, MeV). To derive this potential, one can use strong and TRIV Lagrangians, which can be written as Herczeg (1987); Liu and Timmermans (2004)


where , and are iso-vector and scalar magnetic moments of a nucleon ( and ), and are TRIV meson-nucleon coupling constants.

Then, a TRIV potential obtained from these Lagrangians can be written as

where , , ,.

Comparing eq. (1) with this potential, one can see that functions in a meson exchange model can be identified as

The TRIV potentials in pionless EFT (without explicit pion contributions) contain only point-like nucleon-nucleon interactions which are proportional to the local delta functions. Thus, one can write the corresponding functions as


where low energy constants (LECs) of pionless EFT have the dimension of . Here, for numerical calculations we approximate the singular delta functions by the Yukawa type functions as , as it was done in Liu and Timmermans (2004), with a natural scale of the parameter .

In the pionful EFT, the long range terms of the potential are due to the one pion exchange whereas the short range terms are similar to the ones obtained within the pionless EFT. Then, by ignoring the contribution of the two pion exchange at the middle range scale, as well as other higher order corrections, one can write functions for the pionful EFT as


For this potential, the cutoff scale is larger than pion mass, because pion is an explicit degree of freedom of the theory.1 The following identities might be useful in order to compare this potential with the one used in reference de Vries et al. (2011):


However, the parameters and were not included at the leading order potential in de Vries et al. (2011) because they were considered as higher order terms with additional assumptions related to the source of TRIV interactions. In this paper, we calculate contributions from all the operators without making any assumption about the possible value of the coefficients for each operator.

It is important that all these three potentials which come from different approaches have exactly the same operator structure. Thus, the only difference between them is related to the difference in corresponding scalar functions which, in turn, differ only by the values of the characteristic masses: , , , and . Therefore, to unify the notation, it is convenient to define the new constants (having dimension of ) together with the scalar function (having dimension of ) as


where the expressions of and are provided in eqs. (II), (7), and (II).

Since the non-static TRIV potential terms, with , do not appear either in a meson exchange model or in the lowest order EFTs, they can be considered as a higher order correction to the lowest order EFT or be related to heavy meson or multi-meson contributions in the meson exchange model.

Iii and EDMs

The value of nuclear EDM is defined as


where is a nuclear wave function with a total spin and its projection equal to . The EDM operator contains direct contributions from the intrinsic nucleon EDMs (current operators)


and contributions from the nuclear EDM polarization operator


which describe polarization of the nuclei due to TRIV potentials. Here, and are neutron and proton EDMs, and and are charge and position of -th nucleon. In this work, we do not consider possible TRIV meson exchange current contributions. Therefore, the value of (or ) EDM can be expressed as


where and represent the time reversal invariant and TRIV parts of the nuclear wave function.

First, we will analyze the contribution of the intrinsic nucleon EDM momenta to the nuclear EDM which are defined by the action of operator defined in eq.(12) on the nuclear wave function. For the deuteron (the two-body system), the contributions of the intrinsic neutron and proton EDMs simply add , i.e. this value does not depend on the nuclear wave function and, consequently, on the choice of the particular strong interaction model. The situation is different in the three-body system, or nuclei, where the EDM contributions from intrinsic nucleon EDMs become wave function dependent. One can see that nucleonic contributions to the nuclear EDMs are in rather good agreement for all strong potentials with local interactions: the AV18, the Reid93 and the Nijm II. Nevertheless, these values differ for the models that includ non-locality or for the ones where a three-nucleon force is added (see Table 1 where EDM calculations from the references Stetcu et al. (2008); de Vries et al. (2011) are also listed). Softer two-nucleon interaction models2 have a tendency to provide the nuclear EDM values closer to the ones of the unpaired nucleon (i.e. neutronic EDM for case, and protonic EDM for ). The addition of the three-nucleon force provides an effect similar to the one of hardening the interaction. This effect is clearly related to the strength of the tensor force, which permits one to reduce the pairing of the nucleons in the three-nucleon system. Let us mention that our results for this single-nucleonic operator are in excellent agreement with those from references Stetcu et al. (2008); de Vries et al. (2011).

CD-BONNStetcu et al. (2008); de Vries et al. (2011) -
AV18Stetcu et al. (2008); de Vries et al. (2011) -
EFT NNStetcu et al. (2008); de Vries et al. (2011) -
EFT NN+NNNStetcu et al. (2008); de Vries et al. (2011) -
Table 1: The nucleon electric dipole moment contributions to nuclear EDMs calculated for different strong interaction potentials.

In order to calculate polarization contributions to the nuclear EDM, we solve Faddeev equations in a configuration space Faddeev (1961) by including TRIV potentials. We consider neutrons and protons as isospin-degenerate states of the same particle nucleon whose mass is fixed to MeVfm. By using the isospin formalism, the three Faddeev equations become formally identical, which for pairwise interactions reads


where are particle indexes, is the kinetic energy operator, is a two body force between particles , and , and is the so-called Faddeev component. In the last equation, the potential formally contains both a strong interaction (TRI conserving) part () and a TRIV (parity violating) part (), i.e.: . Due to the presence of TRIV potential, the system’s wave function does not have a definite parity and contains both positive and negative parity components. As a consequence, the Faddeev components of the total wave function can be split into the sum of positive- and negative-parity parts:


Three-nucleon bound state wave function has a strongly predominant positive-parity component. The TRIV interaction is weak (). Then, by neglecting second-order terms in TRIV potential, one obtains a system of two differential equations:


One can see that the first equation (17) defines only the positive-parity part of the wave function. This equation contains only a strong nuclear potential and corresponds to the standard three-nucleon problem: a bound state of helium or triton. The solution of the second differential equation (18), which contains an inhomogeneous term , gives us the negative-parity components of the wave functions.

To solve these equations numerically, we use our standard procedure described in detail in Lazauskas (2003). Using a set of Jacobi coordinates, defined by and , we expand each Faddeev component of the wave function in bipolar harmonic basis:


where index represents all allowed combinations of the quantum numbers presented in the brackets, and are the partial angular momenta associated with respective Jacobi coordinates, and and are spins and isospins of the individual particles. Functions are called partial Faddeev amplitudes. In the expansion (19), we consider both possible total isospin channels and , regardless of the fact that positive-parity components have predominant contribution of state.

Equations (17) and (18) must be supplemented the appropriate boundary conditions for Faddeev partial amplitudes : partial Faddeev amplitudes are regular at the origin


and the system’s wave function vanishes exponentially as either or becomes large. This condition is imposed by setting Faddeev amplitudes to vanish at the borders of a chosen grid, i.e.:


This formalism can be easily generalized to accommodate three-nucleon forces, as is described in paper Lazauskas (2009).

In Table 2, we summarize the calculation of the matrix elements ; values for each TRIV operator from Eq.(II) obtained for a different choice of the strong interaction are tabulated. In this table, operators of TRIV potential (II) are calculated when combined with a unified scalar function


and calculated when taking various values of the parameter , which was chosen to coincide with the masses of , , and mesons. Therefore, this table can be used to analyze the TRIV potentials in the meson exchange model, whereas once multiplied by an additional cutoff factor, it also can be used for the TRIV potentials in pionless EFT or pionful EFT.

operator AV18 Reid93 NijmII AV18UIX INOY
Table 2: Contribution of the different TRIV operators (Eq.(II) to the expectation value of . Calculations has been performed for several different strong potentials and for () nucleus; values are given in e-fm units.

Similar to the intrinsic nucleon EDM contribution, presented in table 1, dynamical nuclear EDMs are also rather insensitive to the choice of the local strong interaction potential. The addition of the three-nucleon force affects the results by 10-20%, whereas the presence of the strong non-locality in two-nucleon interaction (as in the INOY model) has a very large impact. This fact is also confirmed by the cutoff dependence behavior of the relative deviations of the matrix element for operators 1 and 5 calculated for different strong interaction potentials in relation to the matrix element calculated with AV18 potential (see Fig. (1)), where the first plot corresponds to the operator 1 and second one to the operator 5.

Figure 1: The relative deviations of the value from the one obtained for AV18 potential . Results are presented for the operators 1(upper) and 5(lower) and as a function of the cutoff parameter.

Using meson exchange TRIV potential from Eq.(II), we can present the results of the calculations of nuclear EDMs as the sum of contributions of the different TRIV potential terms to the matrix element, as is shown in Table 3.

Couplings AV18 Reid93 NijmII AV18UIX INOY AV18Stetcu et al. (2008)
77.2(-76.9) 79.5(-79.3) 80.0(-79.8) 71.9(-71.6) 134(-134) 157
141(144) 143(145) 145(148) 138(141) 142(145) 288
88.3(-91.8) 88.7(-91.6) 88.3(-91.3) 87.1(-90.1) 98.5(-102) 444
-0.803(0.814) -1.06(1.08) -0.964(0.974) -0.742(0.751) -4.40(4.41) -1.65
1.20(1.21) 1.34(1.35) 1.49(1.50) 1.09(1.09) 2.36(2.37) 2.48
-0.836(0.879) -0.824(0.858) -0.811(0.845) -0.846(0.885) -1.77(1.85) 4.13
1.84(-1.85) 2.05(-2.06) 1.91(-1.91) 1.54(-1.54) 6.88(-6.91) 4.13
-4.33(-4.46) -4.74(-4.86) -5.19(-5.32) -4.27(-4.38) -8.49(-8.71) -9.08
-1.28(1.28) -1.36(1.36) -1.31(1.31) -1.07(1.07) -3.43(3.45) -
2.40(2.41) 2.57(2.59) 2.75(2.77) 2.18(2.18) 3.48(3.50) -
Table 3: Contributions to for EDMs from different terms of meson exchange TRIV potential in e-fm units. We use the following values for strong couplings constants: , , , . A similar table can be inferred for the case of pionless and pionful EFT from Table 2.

The last column in Table 3 shows the results for EDM obtained in reference Stetcu et al. (2008). Comparing results of Stetcu et al. (2008) with our calculations for AV18 potential, one can see that there is a systematic discrepancy for all the values of the matrix elements. For these calculations, we have used the same strong and TRIV potential as in Stetcu et al. (2008). It, therefore, points at the possible systematic error in one of the calculations. We use solutions of the Faddeev equations, while the calculations of the wave functions in Stetcu et al. (2008) have been done in a no-core shell model framework using perturbative expansion for the negative parity states. Based on the presented results, it is impossible to figure out if the discrepancy is the result of a numerical error in one of the algorithms or if there is an intrinsic limitation of the perturbative expansion used in reference Stetcu et al. (2008) with no-core shell model approach3.

On the other hand, our calculations for the deuteron (two-body system) using the same formalism and employing AV18 strong interaction gives


which is in excellent agreement with the result of reference Liu and Timmermans (2004)


Iv Conclusions

Using the data from Table 3, one can obtain the polarization parts of and EDMs for the choice of AV18UIX strong potential. Then, the expressions for and EDMs can be written as




It should be noted that in general neutron and proton EDMs can not be related to the meson-nucleon TRIV constants from TRIV potential. However, it is convenient to present expressions for these EDMs, obtained in the chiral limit Crewther et al. (1979) with the assumption that nucleon EDM are resulted from TRIV potential


which could be used for some models of CP-violation.

Finally, one can compare the obtained expressions for nuclear EDMs with TRIV effects in neutron deuteron elastic scattering related to the correlation, where is the neutron spin, is the target spin, and is the neutron momentum, which can be observed in the transmission of polarized neutrons through a target with polarized nuclei. This correlation leads to the difference Stodolsky (1982) between the total neutron cross sections for parallel and anti-parallel to