1. Introduction
The study of higher-order nuclear processes is a very demanding field of research, especially when they are driven by hadronic interactions. The complexities are introduced by the equal importance of nuclear many-body aspects and the effective nature of in-medium low-energy nuclear interactions. Since both sectors are intimately intertwined, a clear separation of effects is hardly possible, thus inhibiting a straightforward perturbative approach. Still, over the years, nuclear reactions and structure theory have developed a tool box of methods allowing now systematic investigations of rare nuclear processes and their spectroscopy. Phenomenological approaches and the meanwhile rather successful many-body approaches based on effective nuclear field theory have reached a level of accuracy that fine details of nuclear spectroscopy are now accessible—and predictable—by theory. Our recent work on a first-time quantitative description of heavy ion double charge exchange (DCE) data on microscopic grounds [
1] is a prominent example for that kind of achievement.
Although there is a general consensus on the importance of studying higher-order nuclear processes, detailed studies are rare. Recent examples are the multi-phonon description of the extremely rare nuclear double-gamma emission [
2] and, using similar nuclear structure methods, the study of the quenching of low-energy Gamow–Teller strength [
3]. On the experimental side of nuclear reactions, the NUMEN project [
4] is the trend-setting case of a research project fully devoted to a higher-order nuclear reaction, namely to investigate nuclear DCE reactions with heavy ion beams, aiming to make an independent probe for the nuclear matrix elements of nuclear double beta decay (DBD) available.
Higher-order reactions are of interest for nuclear reaction physics. They have a high potential to reveal rare reaction mechanisms hitherto undiscovered because they are not present or suppressed in first-order processes. Heavy ion second-order reactions have rarely to never been used as spectroscopic tools. An exception is statistical multi-step reactions in the pre-equilibrium region of nuclear spectra, which have been studied in the past for light and heavy ion reactions using Multi-step Direct Reaction (MSDR) theory [
5,
6,
7]. In [
8,
9], the MSDR scheme was used to study neutron-induced single charge exchange (SCE) reactions in the continuum region of the spectra in a two-step approach. In the present context, we are interested in a much more selective case, namely on
reactions leading to the discrete part of the spectrum. These reactions take place in a complementary manner in the projectile and the target system. Hence, for a complete description, both nuclear systems must be described and understood simultaneously. From a general point of view, however, that apparent complication can be considered an advantage because it allows us to probe two DBD processes simultaneously in a single reaction, namely a
transition in one nucleus by a complementary
transition in the other nucleus. Thus, DCE reactions pose a double challenge to nuclear theory.
The reaction mechanism of a DCE reaction with composite nuclei is by no means obvious. In principle, DCE reactions can proceed either by mutual nucleon transfer processes or by acting twice with the isovector nucleon–nucleon (NN) interaction. A first detailed discussion on that important issue is found in our recent review [
10]. Historically, after first heavy ion DCE data were measured, a pair transfer scenario was favored, by which DCE reactions are assumed to proceed as a simultaneous mutual exchange of a proton pair in one direction and of a neutron pair in the other direction [
11,
12]. The transfer mechanism is a soft process driven by mean-field dynamics. The minimal scenario is a sequence of two pair transfer reactions, e.g.,
, interfering with a second reaction path where the proton pair is exchanged first. Hence, in leading order, the pair transfer scenario is at least of fourth order in the nucleon binding potentials. Single nucleon exchange processes are of even higher order. Transfer reaction mechanisms are most important in general at low incident energies where the kinematical conditions are favorable for probing mean-field dynamics. We will not consider further transfer DCEs which, in fact, have been found to be negligible for the reactions considered here, as confirmed by recent experimental and theoretical investigations [
13,
14].
For a long time, pair transfer was thought to be the dominant heavy ion DCE reaction mechanism. After the impressive successful use of heavy ion SCE reactions for spectroscopic studies (see, e.g., [
10]), a first attempt to use heavy ion DCE reactions for spectroscopic purposes was made by Blomgren et al. [
15], intending to measure the excitation of the spin-flip double Gamow–Teller resonance (DGTR). However, at that time, the results were disappointing, which led the authors to rather pessimistic conclusions on the usefulness of heavy ion reaction for DCE studies. About a decade later, the situation changed when the feasibility of DCE reactions and their potential for spectroscopic investigations was shown for the reaction
at
270 MeV by Cappuzzello et al. [
16] in an experiment at LNS Catania. That experiment was important for narrowing down the conditions under which nuclear structure information can be extracted from data, as is now the central goal of the NUMEN project [
4].
The observed angular distributions of Reference [
16] in fact show a puzzling similarity to SCE reactions in shape and, to a lesser extent, in magnitude. Moreover, the data cover a surprisingly large range of linear momentum transfer, extending up to about 500 MeV/c over the measured angular range. Thus, these properties demand a reaction mechanism different from low momentum-centered mean-field dynamics. The appropriate candidate is charge exchange by hard collisional interactions, as provided by the mesonic DCE scenario, introduced for the first time in [
17]. Actually, as pointed out in [
10,
17], there are two competing mesonic DCE reaction mechanisms. Here, we consider specifically the double single charge exchange (DSCE) scenario. As illustrated in
Figure 1, the DSCE reaction mechanism is given by two consecutive SCE events, both occurring half off-shell. As was discussed in detail by Bellone et al. [
1], the DSCE process is of second order in the isovector NN T-matrix. The measured angular distribution was described close to perfection in magnitude and very satisfactorily in shape by distorted wave theory, free space NN T-matrices, and microscopic nuclear structure input. Hartree–Fock–Bogoliubov ground state densities were used for the optical potentials and response functions, and transition form factors were obtained by QRPA theory, as discussed before in Reference [
18]. Alternative approaches to the nuclear structure aspects of DCE reactions (as for DBD theory) are of course highly desirable. The interacting boson model (IBM) belongs to the frequently used approaches in DBD theory. Using both eikonal and closure approximation, the aforementioned DCE reaction was analyzed in terms of an IBM–NME by Santopinto et al. [
19].
In [
1], the usual
t-channel approach for calculating form factors was used, according to the scheme displayed in
Figure 1. That formulation is perfectly suited for the proper description of DSCE cross sections if the interest is focused on the reproduction or prediction of cross sections. However, that approach is not suitable for investigations and/or extraction of DBD nuclear matrix elements (NME) from cross sections. The latter are connecting the two SCE-type vertices within the same nucleus, while the standard reaction theoretical approach is directed towards the description of the pair of vertices excited in the projectile and target in the first or the second steps of the DCE reaction. Thus, a change from the conventional
t-channel formulation to an appropriate
s-channel formulation is required, not to the least as a necessary prerequisite for establishing the connection to the NME entering DBD theory.
Keeping this goal in mind, the program of this paper is a purely theoretical one, namely to recast the second-order DSCE reaction amplitude into an
s-channel representation. As seen below, this requires a demanding amount of recoupling of various kinds of angular momenta, including the spectroscopic ones intrinsic to the involved nuclei and those describing the multipolarities acting on the relative motion degrees of freedom. Moreover, the total number of form factors to be considered increases to the fourth power (or stronger) by the number of elementary NN-interaction vertices. Thus, the full account of rank-0 central, rank-1 spin orbit, and rank-2 tensor NN-vertices results in general in a total of at least
form factors, distributed half by half in the projectile and target. In order to keep the presentation at a manageable level, we therefore restrict the discussion to the vertices of the rank-0 central interactions. As was discussed already in References [
1,
18], they involve already the complete set of relevant fundamental isovector vertices, describing non-spin flip
Fermi-type and spin-flip
Gamow–Teller-type nuclear transitions of any multipolarity. The algebraic rules developed below can be used in the same manner also for more extended sets of NN-vertices. Practical calculations, e.g., those in Reference [
1], account of course for the complete set of interactions.
The paper is organized as follows: The reaction theoretical aspects of DSCE reactions are presented in
Section 2, adding additional theoretical background to our recent work [
1]. The DSCE reaction amplitude is derived and discussed in
Section 3, leading to a factorized form that separates ion–ion interactions and nuclear matrix elements. A key element of DSCE theory is discussed in
Section 4. The reaction amplitude and the transition form factors are transformed from the
t-channel to the
s-channel representation, thus recasting the theory into a form compatible with and comparable to the formulations used in DBD theory. The investigations lead also to the result that a rich spectrum of multipoles contributes to a DSCE reaction, confirming our previous numerical results on theoretical grounds [
1]. The physics content of the form factors and accordingly of DCE cross sections is investigated in
Section 5 by considering a few limiting cases. Form factors are derived in closure approximation. A reduction scheme that allows for a first-time derivation of DSCE multipole unit cross sections, which account for the reaction dynamical aspects and may serve in the future to extract DSCE–NME directly from data, is presented. A few representative examples of unit cross sections are shown. The work is summarized and an outlook is given in
Section 6. Certain coefficients resulting from the recoupling of angular momenta are found in
Appendix A.
2. Theory of Sequential Double Charge Exchange Reactions
As depicted schematically in
Figure 1, the double single charge exchange reactions are a sequence of two consecutive single charge exchange processes. After the first SCE event, the system propagates in a combination of
configurations, concluded by a follow-up second charge exchange process. Each of the single charge exchange processes is induced by the two-body NN–isovector interaction
. The T-matrix is used in a form given by one-body operators acting in the projectile and the target nucleus, respectively. For a reaction
, the reaction amplitude is written down readily as a quantum mechanical second-order reaction matrix element [
1]:
Initial (ISI) and final state (FSI) interactions are taken into account by the distorted waves , depending on the center-of-mass (c.m.) momenta and obeying outgoing and incoming spherical wave boundary conditions, respectively. The available c.m. energy is , where .
As discussed in [
18], we use an (anti-symmetrized) isovector NN T-matrix of the form
In non-relativistic notation, the rank-0 central and rank-2 tensor interactions are (In our notation, the form factor of the rank-2 tensor interaction includes an additional factor
.)
As indicated by the dot product, the spherical harmonics has to be contracted with the rank-2 spin tensor to a total scalar. The spin operators act in the projectile and the target. Summation over all target and projectile nucleons is implicit. The form factors and are given by a superposition of meson-exchange propagators connecting projectile and target nucleons. The strength factors are given by complex-valued coupling functionals in general depending on the energy in the NN center-of-momentum frame and the nuclear densities. Other operator structures such as two-body spin-orbit interactions will not be considered but are included easily.
The intermediate propagator for the evolution of the intrinsic nuclear states and relative motion,
is given by the nuclear Hamiltonians
and the relative motion Hamiltonian
. The latter is described by an optical model Hamiltonian,
. With the set of intermediate SCE-type states
and
in the projectile and target, respectively, we obtain
where the channel propagator is
is the total c.m. energies of the intermediate nuclei in states
c and
C, respectively. As discussed in [
1], by means of the bi-orthogonal set of distorted waves
, the reaction amplitude is finally obtained as
where
is the S-matrix element from the dual states
, being solutions of
.
denotes the kinetic energy related to the (off-shell) momentum
. The half off-shell SCE amplitudes are of the form
with the transition form factor
.
The DSCE differential cross section (for unpolarized ions) is given as
averaged over the initial nuclear spin states (
and
) and summed over the final nuclear spin states (
, respectively). Reduced masses in the incident and exit channel, respectively, are denoted by
.
3. The Heavy Ion DSCE Reaction Amplitude
In momentum space, the operator structure of the isovector part of the NN T-matrix is determined in all tensorial parts by the operators [
18]
With the nuclear transition form factors
, the half off-shell SCE amplitudes become
and accordingly for
. The above—on first sight, unusual—form was chosen in virtue of displaying the factorization of the SCE reaction amplitude into a distortion coefficient
, containing the elastic ion–ion interactions and—in brackets—the nuclear transition form factors, describing intranuclear SCE dynamics. The full information on elastic ion–ion interactions is contained in the distortion coefficients:
They can be considered an extension of the S-matrix concept into the off-shell region.
By contour integration, the propagator Equation (
5) is separated into the intrinsic nuclear and the relative motion propagators,
by which a formal separation of the relative motion and intrinsic nuclear evolution is achieved. The integration path
extends over the upper half of the complex
–plane. Applying the momentum representation, Equation (
12), the DSCE reaction amplitude becomes
where contributions of higher rank tensor operators have been left out for the reasons discussed in the Introduction. In numerical calculations, the full spectrum of tensor operators is, of course, taken into account.
The projectile and target sequential SCE responses are now contained in the nuclear polarization tensor:
combining, however, both projectile and target transitions. As indicated by the dot products, a total spin-scalar tensor is obtained.
4. Multipole Structure of the Transition Form Factors and Nuclear Matrix Elements
The result of Equation (
16) is in fact perfectly well suited for DCE reaction calculations, as in [
1]. The focus of this section is to clarify the relation of a DSCE reaction to nuclear matrix elements (NMEs) of the projectile and target. Hence, we develop a formalism by which the contributions of the two nuclei to the combined ion–ion NME can be separated. Such a program requires decomposing and rearranging the nuclear tensor of Equation (
16) in an appropriate manner. The momentum representation provides the suitable formalism.
As seen by the results of the previous section, the polarization tensor is given by products of Fourier–Bessel transforms of transition densities. The momentum structure of the transition densities is probed by the operators
(Equation (
11)). By expanding the plane waves into multipoles, we find the spin-scalar (
) Fermi-like and spin-vector (
) Gamow–Teller-like isovector (
) one-body operators:
where
denotes a spherical Bessel function of order
ℓ. The nuclear SCE form factors become
where the Wigner–Eckardt theorem was used to derive the reduced matrix elements:
We use the notation .
In order to clarify the physics content, we emphasize that the matrix elements (Equation (
19)) are in fact momentum-dependent transition form factors. As such, they do not follow the rules known from beta-decay on the enhancement or suppression of multipolarities already by the operator structure alone. Only for
the transition operators of Equation (
17) approach the long–wave length limit underlying the weak and the electromagnetic operators commonly used in nuclear structure theory. For sufficiently large
p—as easily realized in a heavy ion reaction—essentially all multipole operators are of the same magnitude as dictated by the asymptotics of the Bessel–Riccati functions.
The standard ordering of transitions in a DSCE reaction follows the scheme indicated in
Figure 1, as imposed by meson exchange. However, in order to comply with the goal to access projectile and target DSCE nuclear matrix elements, a regrouping and recoupling of terms and correspondingly of angular momenta is required in order to follow the evolution of the intrinsic nuclear states instead of focusing on meson exchange. In other words, a change in representation from the
t-channel to the
s-channel has to be performed.
The complexities of the second-order process are reflected of course in a correspondingly involved formalism. As a rule of thumb, momenta and quantities such as operators and quantum numbers of spins and angular momenta related to the first and the second SCE vertexes will be denoted by indices 1 and 2, respectively. Quantities related to processes in the projectile or target nuclei are usually denoted by the index a and A, respectively, occasionally complemented by indices and if states in the projectile-like and the target-like intermediate and final nuclei have to be distinguished.
Considering only the central spin–scalar and spin–vector interactions by the reasons discussed in the Introduction, the result is -4.6cm0cm
thus now being in an order appropriate for the separation of projectile and target response functions. A second contour integration is used to separate completely projectile and target NMEs:
where care has been taken in maintaining a total spin-scalar result. The target tensor, for example, is
and the projectile tensor is defined accordingly. Replacing
and re-interpreting
as the lepton energy, a striking similarity to the NME of
decay, e.g., [
20], is immediately identified.
By several steps of angular momentum recoupling, the nominator of Equation (
22) is finally obtained as a superposition of irreducible multipole components:
where we introduced the bi-spherical harmonics:
The reduced form factors themselves are defined by a superposition of multipole contributions:
resulting from the coupling of transferred spin and orbital angular momenta to the total angular momentum transfers
in the first and second SCE interactions, respectively. The recoupling coefficients are defined in
Appendix A.
We introduce the multipole polarization propagators
with the reduced polarization propagators
by which Equation (
22) becomes
Inserting Equation (
28) into Equation (
21) and the corresponding expressions for the complementary sequence
, the summation over the magnetic spin quantum numbers
can be performed (see
Appendix A) and we obtain the intermediate result
The coupling indicated in the last line of the above formula is finally exploited to recouple the two bi-spherical harmonics into a single one, as shown in
Appendix A. As a consequence, the angular dependencies are stripped off the nuclear tensors, and we find
with the transition form factor of total multipolarity
where
is found in
Appendix A. This allows us to define the reduced reaction amplitudes
By comparison to Equation (
16), the essence of the exercise is that we have achieved a reduction in the
nuclear transition tensor to a form displaying explicitly the target and projectile response functions and their multipole structure. In addition, the
and
angular momentum coupling coefficients have been split off such that, for the angular distribution (Equation (
10)), the summations over the magnetic quantum numbers can be performed, resulting in
As a closing remark to this section, we emphasize that the formulation has been kept very general, intending to cover for future use the full multipole spectrum. For special cases, especially for transitions, the situation simplifies considerably. The total angular momentum transfer is constrained to , which for the total orbital and spin angular momentum transfer implies the two combinations and , respectively. In the first case, the intermediate channels are restricted to sequential excitations of Fermi modes or Gamow–Teller modes where and , i.e., the same multipolarity is excited in each of the two SCE steps. The case is accessible only by sequential Gamow–Teller-type transitions of the same total multipolarity, . For a reaction, the combination is forbidden by parity conservation.
4.1. Nuclear Structure Aspects
In order to evaluate the polarization tensors, nuclear wave functions are required for the involved states. Nuclear ground states are described by Hartree–Fock–Bogoliubov (HFB) theory, as discussed in [
21,
22,
23]. The SCE excited states are obtained by QRPA calculations; see, e.g., [
1,
18] for recent results.
The DCE parent and daughter nuclei are connected by an isospin rotation perpendicular to the
-axis. The rotation is such that the isospin in each nucleus is changed by two units but the total isospin as defined by the incident projectile–target system is conserved, of course. Since isospin is a conserved symmetry in strong interactions, we are eligible to expect that the ground states of the parent and the daughter nuclei are related in leading order by a rather transformation, changing by a rotation in quasiparticle space, e.g., a pair of protons into a pair of neutrons (or vice versa). Hence, in the HFB mean-field picture, the final state reached by a
transition is dominantly given by an
-configuration in the valence shells. Thus, assuming for
a
ground state, states in the
daughter nucleus will be considered as a 4-quasi particle configuration with respect to
. Thus, we use
where
is a single quasiparticle operator. The proper normalization is taken care of by the constant
NB. Typically, such two particle–two hole states are rather stable against perturbations. Thus, good approximation admixtures of higher-order configurations, caused by residual interactions inducing core polarization, can be neglected.
In DBD theory, the ground state of the nucleus B is usually treated as the quasiparticle vacuum state of the daughter nucleus. However, in a DCE reaction, that point of view does not match the sequential character of the transition. We emphasize that a DSCE reaction probes the 4 quasiparticle (QP) content of the states in B with respect to the parent nucleus. A clear advantage of that picture is that the whole spectrum of final states is accessible by the same theoretical methods.
The
intermediate states are of a more complex structure. Starting from an even–even ground state—as is common practice—the odd–odd character of the SCE states and collectivity have to be taken into account. Thus, residual interactions have to be included into the theoretical description for which QRPA theory is an appropriate approach. Hence, the intermediate states
of energy
are described by acting with a two-quasiparticle (2QP) QRPA operator
onto the ground state:
with the 2QP operators
and the time-reversed state
. The 4QP states in
B are grouped and coupled accordingly:
where the coefficient
C accounts for recoupling and normalization. Hence, even in the simplest case of a 4QP configuration given by
pairs of protons and neutrons, a rich spectrum of multipolarities is encountered when transformed to the particle–hole representation.
In second quantization, the transition operators (Equation (
17)) become in the
-channel
while in the
-channel, we find
where scattering terms
have been neglected. We use the same angular momentum coupling scheme for the
and the
cases and exploit the reduced matrix elements obeying the relation
. The Bogoliubov–Valatin QP amplitudes are denoted by
and
, respectively. Within the 2QP–representation, the evaluation of the two sequential SCEs transition is a comparatively easy task, especially for a
reference state. For
transitions of
character, the NME of Equation (
19) is
Since the
transitions start from an already excited nucleus, the form factors are of a quite different structure: the form factors of the second
transition are superpositions of contributions given by one 2QP NME times an overlap amplitude of the second 2QP pair with the previous SCE excitation.
The quasiparticle rescattering contributions, neglected here, would lead to form factors involving a quasiparticle from the intermediate phonon and a quasiparticle from the final configuration. Different from the -type form factors, the scattering terms are of the order and , respectively. Thus, these transitions proceed by decomposing the state vectors of the intermediate configurations into their single quasiparticle components, thereby destroying the coherence of the transition.
4.2. Brief on Spectral Properties of DCE Transitions
As an important message from the above results, we notice that the structure of the final B-configurations plays an essential role in selecting the admissible intermediate SCE states. This, of course, affects also the reaction mechanism because the structure of the final DCE state B determines the path through the pool of intermediate SCE states by constraining the accessible multipolarities.
As an example, we consider more closely the DCE reaction
, which was observed a few years ago [
16] and studied theoretically recently in [
1]. The incident channel involves only (
)-shell nuclei. In the exit channel, the (
)-nucleus
and the (
)–hell ejectile
are present. The HFB results predict that
is given with respect to
in good approximation by two hole states in the
-proton shell and two particle states in the
-neutron shell. Hence, the recoupling leads to 2QP proton–neutron states of negative parity, implying a clear preference for negative parity intermediate states in
. On the projectile side,
may be considered in leading order as a
relative to
, as predicted by our HFB calculation. Since only positive parity 2QP-states are involved, this implies a selectivity for a route through positive parity states in
(Spectral distributions for
and
are found in [
18]).