1. Introduction
Heavy ion double charge exchange (DCE) reactions are unique as a new tool for investigations of the rather unexplored sector of higher-order nuclear dynamics. DCE research is of generic interest for nuclear reaction and nuclear structure physics because of its large potential for high-precision investigations of nuclear modes, which otherwise are almost impossible to access. A central topic of this article is to show that DCE physics is going significantly beyond the standard approach to peripheral heavy ion reactions as dominated by mean-field dynamics. DCE research is located at the intersection of nuclear and hadron physics, thus broadening the view on the dynamics of nuclear many-body systems.
In a previous paper [
1], the emergence of an effective isotensor interaction and the role of ion–ion elastic interactions in second-order double single charge exchange (DSCE) reactions were investigated. DSCE reactions proceed by acting twice with the nucleon–nucleon (NN) isovector T matrix, where each of the actions generates a single charge exchange (SCE) transition. It was shown that by proper transformations of the operator structures, defined by central spin–scalar, spin–vector and rank-2 spin tensor interactions, effective operators are obtained, acting as rank-2 isotensor operators intrinsically in each nucleus. In addition, in [
1], the role of initial state (ISI) and final state (FSI) ion–ion elastic interactions was investigated. The DSCE investigations led to three significant and far-reaching results:
ISI and FSI interactions lead to distortion coefficients, which act as quenching factors. As a result, the DSCE reaction amplitude and consequently the observed DSCE nuclear matrix elements are strongly suppressed by orders of magnitudes compared to the results expected without ISI/FSI.
The relative motion degree of freedom induces in DSCE reactions in each nucleus a correlation between the pair of SCE vertices, where the correlation length is determined by the kinematical conditions of the reaction.
The pair of NN T matrices can be recast into a set of spin–scalar and spin–vector rank-2 isotensor interactions, acting in each nucleus as effective two-body interactions and forming together a four-body ion–ion interaction.
In this work, we investigate the competing
Majorana DCE (MDCE) scenario. While DSCE theory is the second-order extension of the conventional direct reaction single charge exchange (SCE) theory [
2,
3], MDCE theory takes a completely different view by describing a heavy ion DCE reaction as a combination of pion–nucleon DCE reactions in projectile and target nucleus. A consequence of such an approach is that the second-order aspects inherent to a DCE reaction are treated on the level of isovector pion–nucleon scattering, giving rise to dynamically created effective rank-2 isotensor interactions in the projectile and target nucleus.
By definition, a heavy ion DCE reaction relies finally on an interaction of rank-2 isotensor character. Hitherto, searches for such a kind of nuclear interaction of generic character have been unsuccessful. To date, the existence of neither elementary isotensor mesons [
4,
5,
6] nor signatures of interactions of that kind in single, isolated nuclei [
7] could be confirmed with convincing certainty. Most likely, rank-2 isotensor interactions do not exist as an elementary mode of their own right. The conditions, however, might change if two nuclei are in close contact as in a peripheral ion–ion collision. In such a situation, an effective isotensor interaction can be generated dynamically as a transient phenomenon. MDCE reactions proceed by virtual pion–nucleon double charge exchange scattering, involving sequences of
and
pion–nucleon SCE reactions. Their proper combination leads finally to virtual (
) pion–nucleon DCE reactions in the reacting nuclei. Under nuclear structure aspects, a DCE reaction is determined by excitations of
and
two particle–two hole configurations in the interacting nuclei.
In the past, pion beams were used extensively for DCE research on nuclei at the Los Alamos Meson Physics Facility (LAMPF) [
8]. LAMPF was shut down a long time ago, but the physics issues studied there have become of renewed interest for heavy ion DCE research. The theoretical understanding achieved at that time for pion–nucleon isovector dynamics [
9,
10,
11,
12,
13,
14] and DCE nuclear structure theory [
15,
16,
17,
18,
19] are worth being rediscovered because they are of high value for research on the MDCE mechanism of heavy DCE reactions.
In [
1,
20], the similarity of DSCE and two–neutrino DBD was emphasized. A special aspect of the pionic MDCE scenario is the striking similarity to the heavily discussed neutrinoless
Majorana DBD (MDBD). That similarity is illustrated in
Figure 1 on the elementary level of virtual weak
gauge bosons and highly virtual strong quark–antiquark
modes, the former materializing into a lepton pair on the mass shell, the latter into a pair of mesons off the mass shell. MDBD is searched for as a possible signature for
Beyond the Standard Model (BSM) physics because MDBD relies on the still hypothetical Majorana neutrinos with the claimed property
, see e.g., [
21,
22]. MDBD would lead to the spontaneous creation of matter in the form of lepton pairs, thus violating lepton number conservation. As will be seen, the MDCE mechanism is described by graphs resembling those of neutrinoless DBD. However, MDCE is determined finally by strong hadronic interactions of a quite different range and strength. Spectroscopically, the same nuclear states as in MDBD are involved, and the transitions are induced by the same kind of isovector multipole operators, exciting spin–scalar and spin–vector modes.
In the forthcoming sections, we present a concise, unified picture of the physics of MDCE dynamics and the relation to neutrinoless Majorana double beta decay (MDBD). The theoretical foundations and methods are discussed much beyond the level presented in previous publications [
3,
23]. The overall aspects, the essential features, and theoretical principles of MDCE reaction physics are presented in
Section 2. As mentioned before, the MDCE reaction amplitude is formally given by one-step distorted wave matrix element. The quenching of reaction yields caused by the strongly absorptive ion–ion optical potential discussed in [
1] for the DSCE amplitude is less pronounced but still a highly important effect of significant strength. Therefore, the role of ISI and FSI is elucidated in
Section 3. A different view on ISI and FSI as vertex renormalization is presented in
Section 4, where we point to the formal similarity of ISI/FSI with the treatment of short-range correlations in nuclear structure calculations, especially also used in DBD theory. The MDCE transition form factors and nuclear matrix elements are investigated in
Section 5. There, we also address in some detail the essential features of the box diagram, introduce the closure approximation, which allows to define second-order pion potentials as effective two-body DCE interactions. Pion–nucleon scattering and the construction of the pion–nucleon T matrix, used to describe the excitation of
and
states, are the subjects of
Section 6. Illustrating numerical results are presented in
Section 7. The connections of DCE reactions to DBD are discussed in
Section 8. A summary and an outlook are found in
Section 9. Additional material on distortion amplitudes, details of the box diagram, the pion–nucleon T matrix, and more on the theoretical background of the pion potentials is presented in several appendices.
2. Theory of Heavy Ion MDCE Reactions
The MDCE interaction process for a reaction
is illustrated graphically in
Figure 2. Formally, the MDCE scenario is described in box diagrams, where the dynamical key elements are pion–nucleon isovector interactions. The reaction is described by a first-order distorted wave (DW) reaction amplitude
. The differential cross section for an unpolarized beam and target nuclei is defined as
The cross section is averaged over the initial nuclear spin states (
) and summed over the final nuclear spin states (
), respectively. Reduced masses in the incident and exit channels, respectively, are denoted by
.
and
are the (Lorentz-invariant) momenta in the incident and exit channels, respectively.
Formally, the reaction amplitude has the structure of a first-order distorted wave matrix element:
The charged pions, described by the propagator
, connect the nuclear transition matrix elements (TMEs). Since they describe the intranuclear DCE transitions, they are the key elements for spectroscopic investigations. We introduce the pion and nucleon isospin operators
and
and rewrite the T matrices as
Since in a DCE transition only the ladder parts
are relevant and pion and nucleon operators commute, we rearrange the ladder operators to pion and nucleon rank-2 isotensor operators of complementary charge-lowering and charge-raising properties:
In r–space formulation, the MDCE amplitude is given by
The distorted waves
with asymptotically outgoing and incoming spherical waves, respectively, depend on the invariant channel momenta
and the channel coordinates
, the latter describing the relative distance between the initial nuclei
and the final nuclei
, respectively. The reaction is described in the ion–ion rest frame.
The distorted waves are of central importance for the quantitative description of direct nuclear reactions like heavy ion SCE and DCE scattering. They account for diffractive and dispersive initial state and final state elastic ion–ion interactions. In direct reaction (DR) theory, they are described globally by complex-valued optical model potentials including the long-range Coulomb potential and real and imaginary nuclear potentials of ranges which are defined by the sizes of the density distribution of the colliding nuclei. A key role is played by the strong imaginary parts. They describe the absorption of the probability flux by the coupling to the multitude of non-elastic channels and as such are essential for a realistic description of the magnitudes of SCE and DCE cross sections and shapes of the related angular distributions; see [
24]. As stated in [
25], a proper treatment of distortion effects was badly missed in the theory of pion–DCE, leaving many open questions for a realistic description of pion–DCE data.
The spectroscopic content of the DCE reaction is contained in the transition kernel
. An instructive and successful approach is to use the momentum representation:
As a remarkable first achievement, we have succeeded in separating the nuclear and relative motion degrees of freedom. The latter are represented by the plane waves depending on the relative ion–ion coordinates in the incident and the exit channels, denoted by
and
, respectively. The MDCE transition form factor
defined by the diagram of
Figure 2 contains as key elements the nuclear transition matrix elements (TMEs)
The TMEs are of central interest for DCE research because they account for the spectroscopy of the reaction. For example, the transition
is induced by two consecutive actions of the pion–nucleon isovector T matrices
, each giving rise to a SCE transition,
and
, respectively. The vertices are connected by the Green’s function
, describing the s-channel propagation, i.e., in the direction of the left and right vertical branches of
Figure 2, of the intermediate
system. The transition
follows the same rules. The TMEs will be investigated in more detail in a later section.
The DCE process is driven by the t-channel exchange of charged pions between the projectile and the target nucleus as indicated by the lower and upper horizontal branches in
Figure 2. In lowest order, the exchange is described by the the symmetrized product propagator
Possible pion–pion and pion–matter interactions are neglected.
As discussed in
Appendix B, in the ion–ion rest frame, the four-momentum
,
is purely space-like, while
includes the reaction Q value. For
, we may safely neglect the Q-value dependence and describe the exchange of both mesons by static pion propagators
In the ion-ion rest frame and at the energies relevant for heavy ion MDCE reactions, the isovector pion–nucleon T-matrix
is described adequately by the operator structure [
26,
27]
Nucleon spin degrees of freedom are involved via the spin operators
. The form factors
depend on the invariant pion–nucleon energy
. They are playing the role of energy-dependent coupling constants.
3. Initial State and Final State Interactions
Before further investigating MDCE interactions and form factors, we must understand first the contributions of ISI and FSI to the reaction process. In momentum representation, the MDCE reaction amplitude attains an intriguing form:
The ion–ion ISI/FSI parts are contained in the distortion coefficients
. The distortion coefficients are 3D Fourier transforms of the incoming and outgoing distorted waves
In
Appendix A, the properties of distorted wave, derived from an optical model wave equation, are investigated in detail. On general theoretical grounds, two important results are obtained, namely, that formally the distorted waves are factorizable into plane waves and residual amplitudes
, which are determined essentially by the half off-shell optical model elastic scattering amplitudes. As the central result, the distortion coefficients are derived in closed form.
Anticipating the results of
Appendix A, we write
In the second equation, the well-known relation
is exploited; see [
24,
28,
29].
By defining the 3D Fourier transforms
the distortion amplitudes become
For vanishing elastic interactions, also the residual amplitudes vanish, and the distortion coefficients approach the plane distribution
For realistic optical potentials, accurately describing ion–ion elastic angular distributions and total reaction cross sections, the residual amplitudes attain values of order unity , resulting in . These results explain the pronounced quenching of the cross sections of heavy ion reactions by orders of magnitudes compared to the yields observed in reactions with particles not suffering from the strong absorption of the incoming probability flux.
As implied by Equation (
14), the MDCE reaction amplitude is determined by the product of the initial and final state distortion coefficients. Together, they form the reaction kernel
From Equation (
23), we find that the total kernel is a superposition of two kernels of the diagonal products of plane wave (PW) and DW distributions and two mixed PW/DW kernels. Combining the latter two into a single term, the MDCE kernel becomes a sum of three distinct terms
The product of plane wave coefficients defines the reaction kernel
By exploiting the properties of the Dirac delta distributions, we find the on-shell relations
where
and
denote the average channel three momentum and the three-momentum transfer of the reaction, respectively.
Thus, in the plane wave limit, the momenta
are fixed unambiguously by the (invariant) momenta in the initial and the final channels as derived in
Appendix B.
The ISI/FSI contributions are contained in the remaining two terms, which are determined by the amplitudes of Equation (
19). Two types of ISI/FSI distortion kernels are found:
The kernel
describes the distortion effects exerted on the reaction by one of channels, while the other channel is in the PW mode, i.e., ISI and FSI act separately.
accounts for the combined action of ISI and FSI. In the momentum space approach, the MDCE reaction amplitude is understood as a superposition of essentially three interfering contributions of different origin and structure but of comparable magnitude:
The PW contribution reflects the bare nuclear transition matrix element before ISI/FSI renormalization. The contributions in the second line introduce ISI in the initial channel while the exit channel is in PW mode, and FSI in the exit channel while the initial channel remains in PW mode. In the term of the last line, ISI and FSI act in both channels simultaneously.
4. A Different View: ISI and FSI as Vertex Renormalizations
A standard problem of nuclear many-body theory is to incorporate interactions from outside of the model space into the operators acting between the states in the limited model space. Formally, the projection techniques going back to Feshbach [
30] provide first insight into the problem of induced interactions. Over the years, nuclear many-body theory has developed powerful techniques on how to incorporate induced interactions as consistently as possible into all parts of the theory. Examples are many-body shell model studies of double beta decay as, for example, in [
31,
32,
33,
34,
35,
36,
37] and under slightly different aspects also in [
38,
39,
40], regarding even neutrino effective masses by induced interaction from the coupling to axions [
41]. For example, a widely used approach, introduced into DBD theory by Šimkovic et al., is the Jastrow method which implements short-range correlations into matrix elements by a function, acting repulsively at small distances.
The considerations which led to Equation (
31) are in fact following the same theoretical rationality as in nuclear structure theory, however, as will be seen, in a complementary manner. In order to recognize the relationship, we recall that the model space of MDCE reactions includes the incoming and outgoing channel configurations, where the incoming nuclei are assumed to be in their ground states and the outgoing nuclei are assumed to be again in their ground states or in a well-identified excited state. In addition, the spectrum of intermediate SCE configurations will contribute. However, the intermediate states are acting mainly as a reservoir of unresolved spectroscopic strength, being responsible in the first place for generating the effective two-body interactions for transitions from the incoming nuclei to the emerging ejectiles. Thus, the explicitly treated model space contains only an extremely small subset of states of the total
configuration space. In nuclear reaction theory, the respective optical potentials account for the induced interaction as far as they affect elastic scattering. Hence, to a large extent, ISI and FSI correspond to induced interactions from the vast background of non-elastic channels. As known from nuclear many-body theory, once effective interactions are important in one sector, they also affect all other sectors of the theory. In particular, transition operators have to be renormalized in accordance with the renormalization scheme. In the above cited works, the proper implementation of renormalization into all parts of the theory is a topic of central importance.
Reconsidering under these aspects the MDCE reaction amplitude, we arrive at the conclusion that in Equation (
31), the distortion amplitudes
are playing exactly that role, namely, to renormalize the SCE vertices in agreement and consistently with the induced ion–ion initial and final state interactions. This is performed in a systematic manner starting from the bare matrix element, represented by the PW amplitude, then renormalizing one of the vertices but retaining the second vertex as a bare vertex, and finally renormalizing both vertices simultaneously. Hence, ISI and FSI account for the proper renormalization of the DCE–nuclear matrix element (NME) under the conditions of a heavy ion nuclear reaction.
While in the nuclear structure context, renormalizations typically refer to short-range effects, ISI/FSI renormalization, however, accounts for scales defined by the ion–ion self-energies, subsumed in the respective optical potentials. A decisive role is played by elastic scattering amplitude as discussed in
Appendix A. The most relevant observable, however, is the total reaction cross section as the measure for the amount of probability flux leaving the elastic channel. The redirected flux is absorbed into channels ranging from transfer channels, which are dominated by mean-field dynamics, and channels where the nuclei are excited inelastically by soft vibrational excitations and giant resonances, eventually leading to fission or fusion, to hard central collisions, possibly upending in the complete fragmentation of the incoming nuclei. Thus, renormalization by optical model interaction is of a genuine character by covering a broad range of nuclear modes and interactions from the soft to the hard scale. That mechanism is not specific for first-order DW reactions as considered here. As discussed in [
1], a similar renormalization scheme is also present in the second-order reactions double single charge exchange (DSCE) reaction. In DSCE reactions, the matrix elements, however, are renormalized by second-order distortion amplitudes.
7. Numerical Studies
7.1. Pion–Nucleon Partial Wave Cross Sections
Representative results illustrating the quality of the description for P- and S-wave total cross sections are shown in
Figure 5. The reference data from explicit coupled channels calculations are surprisingly well described, especially in view of the simplicity of the potential approach. In detail, the Delta and the Roper resonances are well reproduced as is the case for the
and
S-wave sector. In the S-wave spectra, the highly disputed
resonance is most prominently visible as a rather narrow structure on a non-resonant background. Interestingly, the
peak is largely the result of interferences of a virtual s-channel state with the smooth t-channel background. A long tome ago, the same explanation was already obtained in coupled channels calculations [
54], and more recent studies have come to similar conclusions. The present potentials model results may be taken as an interesting independent confirmation of the earlier CC results. Overall, the agreement of the present results with the CC-generated reference data is surprisingly good in view of the extremely simplified model. Larger deviations occur in the S-wave spectra. Close to the threshold, the S-wave cross sections show some deficiencies, and deviations are seen in the
channel also towards the highest considered energies. They are, however, of minor importance for the present use in MDCE studies because
is dominated by P-wave interactions.
Partial wave total cross sections are defined by the imaginary parts of the scattering amplitudes,
. Thus, a first important test of the reliability of the model calculations is to compare real and imaginary parts of scattering amplitudes. The agreement between OP and CC scattering amplitudes is very satisfying. An example is shown in
Figure 6, where the
and the
partial waves scattering amplitudes are compared to the corresponding CC amplitudes.
The P-wave cross sections and scattering amplitudes are slightly better reproduced than the corresponding S-wave quantities. The CC calculations show that the S-wave components, which are generally located at higher energies, are strongly affected by coupled channel dynamics. Physically, an important source of CC effects are multimeson decay channels, either by direct decay, possibly passing through intermediate heavy mesons, or sequentially by decay chains passing through lower lying resonances, e.g., . Such details, of course, have not been resolved in the present approach but are taken into account globally by the dispersive parts of the partial wave potentials.
7.2. Construction of the Pion–Nucleon T Matrix
In order to construct the pion–nucleon T matrix, Equation (
13), we need to determine the three vertex form factors
. That goal is achieved by considering the partial wave structure of the T matrix and collecting terms of the proper multipolarity and dependencies on the nucleon spin. That task is well documented in the literature, e.g., [
55] and reviewed briefly in
Appendix D.
In the energy region of our interest, the vertex form factors are obtained with sufficient accuracy by the two S-wave amplitudes
and
and the three P-wave contributions,
,
, and
, respectively. Within this basis, the form factors are
where
is the invariant pion–nucleon three momentum.
The partial wave-scattering amplitudes are normalized to units of 1/MeV. By means of the kinematical factor , the T-matrix amplitudes are normalized to units of 1/MeV. is the pion–nucleon reduced mass and denotes the invariant relative pion–nucleon momentum. For the numerical results displayed below, we follow, however, the widely used practice to present the form factors as function of the pion kinetic energy in the laboratory frame, which is obtained by .
Although each of the (complex-valued) partial wave-scattering amplitudes varies considerably with energy as seen in
Figure 6, their superpositions are much smoother functions as
Figure 7 and
Figure 8 confirm. By multiplication with
, the units may be changed to MeVfm
, which is a typical unit for volume integrals and momentum space form factors of NN interactions.
7.3. Extrapolation into the Subthreshold Region
The most important advantage of the OP approach for MDCE theory, however, is to have at hand a method which allows to extrapolate reliably and easily into the subthreshold region. As illustrated in
Figure 7 for the
partial wave, three different sheets are covered kinematically. The sheets are distinguished by the values of the invariant relative pion–nucleon momentum
:
In the physical region, and the invariant momentum and are positive.
In the interval , one finds and .
If also , positive values of are recovered but remains negative.
In
Figure 7, it is seen that the T matrix changes in a characteristic manner: real and imaginary parts are non-vanishing in the physical region while in the first subthreshold sheet, the imaginary parts vanish but recover as soon as the second subthreshold sheet is entered.
Figure 7.
The T-matrix in the kinematical regions relevant for MDCE reactions. Real and imaginary parts from the potential model (OP) are shown for energies above threshold, , , and the two subthreshold regions , and , .
Figure 7.
The T-matrix in the kinematical regions relevant for MDCE reactions. Real and imaginary parts from the potential model (OP) are shown for energies above threshold, , , and the two subthreshold regions , and , .
Because of the intrinsic momentum spread introduced by ISI and FSI, in a heavy ion MDCE reaction, in principle, all three kinematical sheets will be visited while propagating through the intermediate s-channel pion–nucleon systems. In other words, ISI and FSI lead effectively to a sampling over the distribution of MDCE box diagrams of different kinematical and dynamical content.
The vertex form factors
are shown as functions of the pion energy in the laboratory system in
Figure 8. When traversing the boundaries between the kinematical sheets, the amplitudes develop cusps. In the
amplitude, defined by the S-wave scattering amplitudes, the cusps are most pronounced, while they are washed out in the P-wave amplitudes
. A closer inspection shows that the P-wave amplitudes are, in magnitude, about a factor of 1.5 to 2 times larger than
. That difference will be enhanced further in matrix elements by the fact that in
, the P-wave terms scale by
for SCE transitions and even by
in DCE transitions. Thus, already from these considerations, we expect a prevalence of the momentum-dependent P-wave terms in a DCE reaction.
Figure 8.
The pion-nucleon vertex form factors (left), (center), and (right) are shown as functions of the pion energy in the laboratory frame. The imaginary parts of vanish in the physically inaccessible region, where the invariant Mandelstam energy as demanded by the analytic properties of the T matrix.
Figure 8.
The pion-nucleon vertex form factors (left), (center), and (right) are shown as functions of the pion energy in the laboratory frame. The imaginary parts of vanish in the physically inaccessible region, where the invariant Mandelstam energy as demanded by the analytic properties of the T matrix.
7.4. Form Factors of the Pion Potentials
Since the pion–nucleon T matrix, Equation (
13), consists of three terms, the pion potential
, Equation (
42) or Equation (
44), respectively, is in general, in either version, a superposition of nine terms
,
, which depend on the three-momenta
and
. Likewise, because of
, we may choose one of the momenta and the three-momentum transfer
of the reaction and momentum variables.
A simplification is obtained for vanishing total momentum transfer
, which implies the collinearity of the momenta,
and
. Then, the number of elements reduces to six independent scalar form factors
,
. Under these conditions, we find the diagonal potentials
and three non-diagonal potentials
For simplicity, the potentials are evaluated numerically for the special case that and are collinear as well, and implies .
By expressing the sine and cosine functions in terms of Legendre polynomials or Legendre functions, respectively, the angle integrations can be performed in closed form. The momentum integrals are regularized by dipole form factors with cut-off
MeV/c. The resulting k integrals, given by products of ordinary or spherical Bessel functions and Legendre functions of the second kind, all combined with powers of
k, have to be evaluated numerically. The full propagator, Equation (
33), is used. Excitation energies, however, are neglected, which is justified in view of the rather weak dependence on energies well below the pion rest mass.
Typical results for the pion potentials
for the reaction
O
Ca at
MeV are shown in
Figure 9 and in
Figure 10, respectively. As discussed above, ISI and FSI favor momenta
which are centered around the on-shell momenta of the entrance and exit channels,
MeV/c. Accordingly, the potentials are displayed at
and
. In magnitude, the potentials increase with momentum, which seems to be especially pronounced for the P-wave parts
. However, as a look to Equation (
50) and Equation (
53), respectively, reveals, the enhancement is largely due to the explicit dependence of the P-wave potentials on powers of
p. Compared to that dependence, the S-wave form factors
remain in small-to-moderate magnitude. For
, the enhancement effect decreases, and the S-wave potentials become relatively more important. Comparing the oxygen and calcium potentials, one observes a rather mild dependence on the nuclear system as is expected for a short-range phenomenon.
The s-channel
exchange establishes in fact a rather tight two–nucleon correlation. Overall, the range of the potentials rarely reaches 40% of the range of pion exchange
fm. Hence, the MDCE process is of a pronounced short-range character. The correlated pair of SCE vertices acts as a virtual, polarized pion dipole source. Comparisons of the data of the DCE reaction induced by
O
Ca at
MeV can be found elsewhere [
23].
7.5. Transition Matrix Elements
In closure approximation and with the pion potential formalism, the TMEs are obtained in the condensed form
Thus, the transition
is described by a sum of nine partial TMEs
which are determined by the transition potentials
. They are defined and studied in detail in
Appendix E. There, it is also shown that the useful and successful approach is to express momentum and spin operators in the basis of spherical unit vectors. In that basis, one finds that
are dyadic tensor forms. The x-dependence is given by Yukawa-type form factors of a rather short range of less than half of the range of a (static) pion-exchange potential. Hence, using contact interactions might be a meaningful approximation which, however, will not be considered further here.
The two-body operator connecting, in Equation (
57), the initial and final states is in fact separable into one-body operators. That property is evident for the plane wave factor, considering that
, and also the potentials
are given by products of one-body operators. In practical calculations, the plane waves are expanded into partial waves in
and
, and by the formalism introduced in
Appendix E the potentials can also be treated accordingly. At the end, Equation (
57) reduces to a (finite) sum of a number of multipole components which are determined by the angular momentum and parity selection rules of the DCE transition
.
With the bi-spherical harmonics
we find
The multipole TMEs are given by rank-2 isotensor two-body multipole operators
For
, the MDCE transition operators are of a spin–scalar character, and the matrix elements describe non-spinflip double-Fermi (FF) excitation. The FF modes are described by spin–scalar one-body operators which are given by Riccati–Bessel functions
:
For , spin–vector transition operators are encountered which give rise to double excitations of Gamow–Teller (GG) modes, which include a spin–vector transition of natural and unnatural parity. If but or and , we encounter two-body operators of mixed spin–scalar/spin–vector structure, leading to mixed FG and GF excitation by combination of Fermi and Gamow–Teller modes.
The GG and mixed FG/GF modes are described by spin–vector one-body operators. Their derivation and especially proper implementation into the theory requires a remarkable amount of angular momentum recoupling. The spin–vector formalism for DCE reactions was studied in detail in [
1] and will not be considered further here. As was shown also in [
1], the GG operators support total spin transfers
, to be combined with the total orbital angular momentum transfer
to total angular momentum transfer
. That leads to a rich spectrum of transitions, e.g., a DCE reaction with
may proceed by
and
partial contributions.
7.6. Transition Matrix Elements in Collinear Approximation
For arbitrary values of
and
, the evaluation and the practical handling of the TMEs are theoretically and numerically formidable tasks. The efforts, however, are substantially reduced for collinear external momenta, i.e.,
and also
and
are collinear. Further simplifications are obtained by imposing, in addition, the stronger constraint
, which implies
and
. Then, Equation (
57) simplifies to
where
contains the plane wave factor. As an example, we consider the
component of
double-Fermi transitions. Hence, only the spin–scalar S-wave (
), P-wave (
) and the mixed S/P-wave parts (
) and (
) are considered. In
Appendix F, the spin–scalar collinear transition potentials
are derived, and their multipole structure is investigated. For
transitions, the complexity of the potentials is reduced further. For that case, explicit expression are found also in
Appendix F.
Following [
56], we assume that the states in the DCE daughter nucleus
B are obtained by acting with appropriate many-body operators on the ground state of the parent nucleus
A:
where higher-order quasiparticle configurations may contribute but will not be reached in leading order by the DCE transition operators. The same set of operators and the underlying basis of single article wave functions are used to express the transition potentials in second quantization. In practice, nuclear ground-state properties are described with Hartree–Fock–Bogolyubov (HFB) theory, and Quasiparticle Random Phase Approximation (QRPA) is used for excited SCE-type states, see [
2].
Without going further into the details of the nuclear structure approach, the essence of the approach is that the TMEs, Equation (
62), are given by nuclear transition form factors
The TMEs are obtained by the scheme developed in
Appendix F. As a recipe, we have to replace in the expression derived in the appendix the operators
by the Fourier–Bessel form factors
and finally perform the momentum integrals.
In
Figure 11,
Figure 12 and
Figure 13, partial TMEs, Equation (
60), of
transitions in
O
Ne and
Ca
Ar are shown. The TMEs are relevant for the DCE reaction
Ca
O
Ne
Ar studied in [
23,
57]. State-independent average transition densities are used, which are averaged over the spectral distributions and normalized to the respective non-energy weighted multipole sum rule, corresponding to the
unit strength form factors introduced in [
56]. Hence, the results are representative of monopole FF modes in
Ne and in
Ar, relative to the respective parent nuclei.
Comparing the results, the most outstanding feature are the differences between the S-wave and the P-wave TMEs. The double S-wave TME,
Figure 11, contribute only at small momenta close to the threshold. The TME involving P-wave amplitudes,
Figure 12 and
Figure 13, increase strongly with momentum, exceeding the strength of the S-wave TME by large factors. The P-wave enhancement is largely an effect of the additional polynomial momentum dependencies up to order
, see
Appendix F. However, it has to be remembered that the shown results are the bare TMEs before ISI/FSI renormalization. After renormalization, i.e., in a full distorted wave calculation, the high-momentum regions will especially be quenched in addition to the overall reduction by about two to three orders of magnitude, thus considerably damping the apparent enhancement.
For arbitrary total angular momentum
, various combinations of partial contributions of angular momentum
and
are allowed, constrained, however, by parity,
, and otherwise limited only by the shell structure and other related properties of the nucleus under consideration. For the
case, this means that in principle, all pairs of transition densities of equal angular momentum
may contribute. For the
case, this property of the TME is illustrated in the figures by showing the partial TME, Equation (
60), for
. In magnitude and shape, the partial TMEs are rather similar. Thus, we conclude that the MDCE operators support a large spectrum of multipolarities as is typical for short-range dynamics.
An eye-catching feature visible in all plots is the kinks. They appear at the momenta where the intermediate channels cross from below the on-shell boundary, which produces a pole in the propagator. That happens at MeV/c and MeV/c for F and K, respectively. Another feature is the crossing of the on-shell boundary of the subsystems at the slightly smaller momenta MeV/c for A = 18 and MeV/c for A = 40. The location of these thresholds depends however, on the modeling of in-medium pion dynamics, which here is not considered, as mentioned before. Above these momenta, the potentials develop imaginary parts of moderate strength which are not shown here.
9. Summary
A generic feature of heavy ion DCE reactions is the versatility of reaction mechanisms by which the transition from the initial to the final channel can proceed. Occasionally, the related ambiguities are considered a severe disadvantage of research with heavy ion beams. That point of view is much too pessimistic because in reality, it is of advantage to be able to investigate all facets of a physical system under the same, well-defined experimental conditions and describe the results consistently by the theoretical apparatus of nuclear many-body theory. The theoretical task and challenge is to overcome the traditional separation of nuclear reaction and nuclear structure physics. Heavy ion DCE physics demands a combined approach as indispensable for any research on quantum mechanical many-body systems. In the NUMEN project, this decisive aspect is realized by the multimethod approach as discussed in [
23].
In this work, we investigated the theory of the Majorana (MDCE) mechanism which is an especially interesting part of heavy ion DCE reactions. As it was emphasized repeatedly, MDCE theory requires to go much beyond traditional concepts of nuclear reaction and structure theory. First of all, as a hitherto never considered aspect, the MDCE scenario relies on pion–nucleon dynamics which—from the beginning, most likely unexpected—suddenly involves subnuclear degrees of freedom as nucleon resonances into a low-energy nuclear process. Already, that aspect makes it worth the effort of investigating DCE reactions.
The MDCE process relies on a hitherto unknown mechanism, namely, a dynamically induced rank-2 isotensor interaction. One of the central results was to introduce the MDCE closure approximation, which allowed to derive pion potentials and two-body nuclear matrix elements connecting directly the entrance and the DCE exit channels. The pion potentials include combinations of spin and momentum scalar parts, spin–scalar longitudinal and spin–vector transversal momentum–vector components, all attached to a rank-2 isotensor operator. This rich operator structure allows widespread spectroscopic studies, allowing a detailed tomography of the nuclear wave functions. However, experimentally and theoretically such studies are highly demanding because they require to observe, analyze, and interpreted energy–momentum distributions over large ranges.
An especially appealing aspect of heavy ion DCE physics is the conceptional closeness to double beta decay research. That relationship was elucidated in some detail by considering the deeper levels underlying weak and strong DCE processes. They meet at the level of QCD and electro-weak physics. Clearly, neither DBD nor DCE reactions proceed at those fundamental levels. Rather, both types of process are determined by low-energy realizations of the two fundamental theories of the current standard model of physics. However, the comparison of weak and strong DCE processes at the fundamental level is helpful to understand that nuclear DBD and nuclear DCE phenomena are finally nothing but two realizations of the same kind of fundamental processes. The differences in dynamics and strengths of DBD and hadronic DCE are due to the breaking of the fundamental symmetries in our physical low-energy environment.
As an interesting outlook to future work, the closeness of low-energy DBD and DCE physics was elucidated further by pointing to another competing reaction mechanism in heavy ion DCE reactions. Nothing forbids MDCE reactions from proceeding by the exchange of leptons. Leptonic MDCE proceeds by electro-weak dynamics but relies on diagrams of the same topology as investigated in this paper in detail for hadronic MDCE. It is left for future work to understand the dynamics and physics of lepton MDCE in detail and explore the competition of the two seemingly very different but interfering types of weak and strong MDCE reaction mechanisms.
As a closing remark, we emphasize again that in MDCE, reactions are not governed by NN interactions as is the case for DSCE reactions. MDCE reactions are determined by pion–nucleon interactions, which provide the required isospin operator structures for an effective rank-2 isotensor interaction. In MDCE reactions, the colliding ions dynamically generate their own and specific isotensor interactions. Charge and baryon number conservation and isospin symmetry require that a transition in one nucleus must be accompanied by a transition in the other nucleus. Obviously, all of the involved transitions are allowed and possible by strong nuclear interactions. Hence, hadronic DCE is not suppressed or even forbidden by violating fundamental laws of the standard model as required for decay. While MDBD is constrained trivially to appear on the mass shell, MDCE reactions take advantage of the presence of another nucleus which gives access to a broad spectrum of off-shell processes and new research opportunities. Hence, it depends on our theoretical and experimental skills to identify and prepare the proper conditions under which rare hadronic or even leptonic MDCE events will become observable.