Breakup reactions and their ambiguities

We review the ambiguities in the nuclear information extracted from breakup reactions, focusing on those originating from the description of the reaction mechanism and the overall ambiguity inherent to their interpretation in terms of shell occupancies. We present the current discussion about nucleon knockout reactions and how the understanding of the reaction mechanism would help reducing uncertainties. For the former, we consider the case of 11\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{11}$$\end{document}Li, due to the existing large data set. For the latter, we recall the paradigmatic example of the electro-dissociation of the deuteron to address the question of the scale and scheme dependence from the theoretical framework used for the interpretation.


Introduction
Few years after the discovery and first studies of halo nuclei [1][2][3][4][5], Bertulani and Hussein investigated the dissociation of neutron-rich nuclei from secondary beams impinging on various targets [6]. They concluded, at that time, that the interpretation of almost all recent experimental studies with secondary radioactive beams is ambiguous. We take this statement as a starting point to address the following question thirty years later: how ambiguous is the interpretation of breakup (dissociation) reactions?
To investigate the sources of ambiguity, we review the most recent experimental work and conclusions about the structure of 11 Li to quantify how much various interpretations differ on the same nuclear system as discussed in Ref. [6]. We also overview the sources of ambiguities in the nuclearbreakup-reaction mechanism as they are discussed today. To illustrate the model-dependencies of orbital occupancies, we remind the historical case of the D-state probability of the deuteron and the "simplest" breakup reaction: the electrodisintegration of the deuteron. a e-mail: aobertelli@ikp.tu-darmstadt.de (corresponding author)

General description of break-up reactions
In order to present the sources of ambiguity let us introduce a schematic expression for the T-matrix corresponding to the removal of a nucleon N from the projectile P leaving a core C (P = N + C) via the interaction with a structureless target T , for simplicity: where R is the coordinate between projectile and target and K their asymptotic initial relative momentum, r is the coordinate between nucleon and core and ξ corresponds to all internal coordinates of the core. Φ P is the wave function of the projectile and Ψ is the exact final wavefunction. Usual assumptions taken in the analysis of nucleon removal reactions are to neglect exchange processes, where the emitted nucleon is not the one that interacts with the target, and to prevent the core's internal degrees of freedom from being excited during the reaction, which results in an effective V CT (R, r) core target interaction that does not depend on the internal degrees of freedom of the core. This results in the wave function for the final channel Ψ (R, r, ξ) being factorizable Ψ (R, r, ξ) = ψ(R, r)Φ C (ξ ) in a 3-body wave function ψ(R, r) and the wave function of the residual core Φ C (ξ ): We note that the integral over ξ only involves Φ P (r, ξ) and Φ C (ξ ): Moreover, we can express this integral through second quantization: where a(r) is the annihilation operator for the particle of interest in a position r, obtained from a Hamiltonian that should be consistent to that used to extract Φ C (ξ ) and Φ P (r, ξ). This allows us to recognize the definition of the spectroscopic factor It is usual to define "single-particle spectroscopic factors" based on some effective Hamiltonian (usually a shell model) through where α corresponds to a certain single-particle state. This allows to expand ψ NC (r) as a sum of normalized singleparticle wave functions φ α NC times With this we can express the original T-matrix as a sum for each single-particle configuration which translates to: Note that this expression is general for any removal process in which the core can be assumed to be a spectator and shows the general procedure of separating the ingredients from nuclear structure (S F α ) and nuclear reaction (σ α ). It also shows that uncertainties in one of the terms directly affects the other when comparing to experiment. It also shows the dangers of inconsistency, as both S F α and φ α NC (included in T α ) originate from ψ NC (r) and should be evaluated consistently. However, it is habitual to extract both terms from different prescriptions (shell model and the eigenstate of a Woods-Saxon potential, respectively), while the reaction mechanism is also derived from another Hamiltonian. Equation 4 also allows us to consider the non-observability of the spectroscopic factors [7][8][9]. Indeed a unitary transformation U (λ) produces a change in ψ NC (r): which is reflected in its norm, the spectroscopic factor. It is specially interesting to remark that the evolution of operator a α [10]: which would provide an unchanging ψ NC , clashes directly with the original assumption that the reaction process does not alter the nucleons in the core, thus establishing a connection between the ambiguities in the structure of the nucleus and the assumptions in the reaction mechanism. This poses a formidable challenge to the description of the reactions, as the mechanisms involving the excitation and de-excitation of the nucleons of the core are not well understood, and the few cases that can be tackled via ab initio calculations are very computationally expensive. It is however a general consensus that the description of the nuclear structure through standard shell-model calculations is compatible with the simple spectator-core picture assumed above, with a wide range of reactions described under this assumption. It is however questionable how quantitative these interpretations can be and if, in more extreme cases such as the removal of deeply bound nucleons, the qualitative picture still holds.

The structure of 11 Li
An illustrative example of the ambiguities appearing in the extraction of nuclear structure information from experimental data may be found through a brief and curated, but of course sorely incomplete, story of the study of the 11 Li nucleus. 11 Li has been the subject of intensive study for more than 30 years, since the seminal work of Tanihata and collaborators [2] hinted at its halo structure through its unusually large interaction cross section. Given the large data set on 11 Li, which have been analyzed through multiple models, generally using the factorization presented in the introduction; the current uncertainties in the description of 11 Li serve as an illustration of the magnitude of ambiguities found in the standard analysis of nuclear experiments.
The first experimental result we will be focusing on was published only three years later [11], and corresponds to the production cross section and the momentum distribution of 9 Li after the removal of two nucleons from 11 Li via collision with carbon and lead targets. The momentum distribution, obtained with the carbon target target was interpreted as the sum of two gaussians, one narrow (σ =23 MeV/c) and one wide (σ =95 MeV/c). The narrow component would originate from the removal of the two neutrons from the halo, while the wide component would be produced from neutron removal from the 9 Li core. The above-mentioned seminal paper by Bertulani and Hussein [6] showed that the momentum distribution, together with the dissociation cross section [12], could be explained qualitatively through two different models, one that assumed clustering of 11 Li into a core of 9 Li and a dineutron and the other that viewed the dissociation as the decay of a collective dipolar excitation in 11 Li. The fact that such different reaction mechanisms are both able to describe the considered experimental data shows the importance of pinning down the reaction mechanism before trying to extract nuclear information from the data.
Through exploration of the ( 11 Li, 9 Li) reaction with multiple targets at multiple energies [13,14], it is now well accepted that the reaction takes place via direct dissociation of the two halo neutrons due to the nuclear or Coulomb potential depending on the target [15]. Multiple theoretical studies with this description of the reaction mechanism have successfully reproduced the momentum distributions and other observables [16][17][18][19][20]. In Fig. 1, we compare the data from [11] to an eikonal calculation which follows closely the procedure from [21], in which the formalism from Hussein and McVoy [22,23] is employed to compute the stripping of the dineutron cluster, whose wavefunction with respect to the 9 Li core is obtained via integration of the 3-body wavefunction for 11 Li over the n − n coordinate. The 3-body wavefunction is computed using the P1I parametrization from [24,25] while the density for 9 Li and target have been taken from a Hartree-Fock calculation using the SkX effective interaction [26] and the density for the dineutron has been extracted from the 3-body wavefunction. As can be seen in Fig. 1, the calculation manages to reproduce rather well the experimental data, although it undershoots the wide component for large momenta. This is not surprising, as the used 3-body model does not allow the core to be excited and its neutrons to be emitted, while in [6] it is suggested that the wide component corresponds to the removal of neutrons from the core, though other works point to 3-body effects also influencing the wide component [16][17][18]. Given the rather approximate way in which the interaction between the target and the halo neutrons is treated, a better agreement perhaps should not be expected.
Even within a specific theoretical framework and when the reaction mechanism is believed to be well understood, there can be other sources of ambiguity in the analysis of experimental data. Continuing with 11 Li, we focus now on the analysis of the 11 Li( p, pn) 10 Li * reaction. The unbound nucleus 10 Li has been extensively studied during the years, as can  9 Li for the reaction ( 11 Li, 9 Li) on a carbon target at 780 MeV/nucleon. Experimental data are from [11] (see details in text) be attested from Table 9.3 from [27] but there are still open questions about its structure. In [28], the 11 Li( p, pn) 10 Li * reaction was measured at 280 MeV/nucleon and the energy distribution of the residual 10 Li was obtained. For the analysis of the spectrum, the authors opted for an R-matrix fit, from which they obtained a p-wave resonance at ∼0.5 MeV and a virtual state with a scattering length of -22 fm. Because the R-matrix fit assigns arbitrary weights to each component, they could not extract the contribution of each state to the ground state of 11 Li, though a later analysis [29] obtained relative weights between s and p-wave by studying the profile function. However, as pointed out in [25], the energy of the different states of 10 Li will be correlated to their contribution to the ground state of 11 Li, so an independent fit of the resonance energies and their magnitudes may lead to unphysical values for them. In fact, in [25], the authors presented a method where they extracted the energy distribution for different states of 10 Li * using a three-body model for 11 Li [24] and the Transfer to the Continuum formalism [30] to describe the reaction process to each state of 10 Li, assuming the 10 Li acted as an spectator during the reaction. Figure 2 is an adaptation of their results. P3 corresponds to a description of 11 Li which reproduces the energy of the pwave resonance and the scattering length of the virtual state from [28] with central and spin-orbit potentials between 9 Li and the two neutrons. As can be seen in the top panels in the figure, this description of 11 Li fails to reproduce the experimental results, because the contribution of the s-wave virtual state is too large. Given the strong relation between the position of the virtual state and its contribution to the 11 Li ground state a better agreement with data requires a modification in the model for 11 Li, not only in its parameters. In this case, the interaction P1I considered the splitting of p-wave resonance and s-wave virtual state through the coupling between the spin of 9 Li and the spin of the neutron in 10 Li, while   10 Li after the 11 Li( p, pn) reaction at 280 MeV/nucleon. The panels to the left correspond to calculations without folding with experimental resolution while those to the right include it. Experimental data are from [28]. Figure adapted from [25]. P3 correspond to a simple model reproducing the established energy of the p-wave resonance and the scattering length of the virtual state, while models P1I and P5 are extensions of this model. P1I considers splitting due to the spin-spin interaction between 9 Li and neutron, while P5 considers a low-energy d-wave resonance as suggested in [31] P5 included a d-wave resonance at 1.5 MeV, as suggested in [31]. Both models manage to reproduce data remarkably well, by draining strength from the virtual state. However, this observable is not able to distinguish between them, due to the experimental resolution washing out the d-wave resonance. The same models have managed to reproduce the angular distribution of the 11 Li( p, d) reaction [24], while an analysis of the 9 Li(d, p) reaction seems to favour the model P1I with spin splitting over model P5 with a d-wave resonance [32]. This example points to possible ways to disentangle ambiguities in the analysis of experimental results: more precise measurements of the same observable and different observables that highlight different properties of the nucleus. Turning back to the original discussion of the 9 Li momentum distribution, a better treatment of the three-body projectile, beyond the dineutron clustering presented in Fig. 1 would require a description of the separate interactions between the target and the three components of the projectile, considering both the absorption of one of or the two halo neutrons. Hussein published recently some works to this respect [33]. This relates to the question on whether the reaction happens sequentially (one neutron is removed or absorbed and the remaining 10 Li decays) or simultaneously (both neutrons are removed in the reaction). To this point, it is valuable to mention the work in [34], where the 2n removal from 14 Be was analyzed through Dalitz plots for Pb and C targets at low energies and the signature of simultaneous breakup was found for the Pb target, while sequential breakup was found for the C target. Given the similarities between the Borromean 14 Be and 11 Li it is expectable that these results are applicable to 11 Li. This also leads to the interesting question of the prevalence of the dineutron (spatially-correlated pair of neutrons) in the 11 Li halo. The study of the dineutron has not been bereft of discussions and controversies, with some works [15,[35][36][37] finding evidence against dineutron in 11 Li while mounting evidence during the years [38][39][40] points to 11 Li indeed presenting the spatial neutron-neutron correlations characterizing the dineutron. It should be remarked that in [38] the source of the discrepancies between different experiments was pointed out to originate from the different efficiencies of the experimental setups at low 11 Li energies, which shows the utter importance of adequate characterization of experimental setups when extracting information from experimental data. A very recent work to this subject has been published [41] analyzing neutron-neutron angular correlations through the 11 Li( p, pn) reaction, where they have found the angular correlation to be dependent on the missing momentum of the reaction, with the neutron-neutron correlations being strongest for the neutrons in the surface of 11 Li. In this work, the composition of the ground state of 11 Li has been compiled from different experiments and theoretical calculations. Their compilation is presented here as Table 1. The spread of results in terms of orbital occupancies of the two neutrons of the halo, with typical absolute variations of ±10% from study to study, for s, p and d orbitals, and uncertainties varying from 1 to 10%, shows that despite the great advances in the understanding of 11 Li and the study of the breakup reaction mechanism, significant uncertainties still remain. It is not surprising to find discrepancies for the occupancies, since, as seen in the previous section, the extraction of these quantities from experimental data depends on the model used in the theoretical description.

Open questions on nucleon knockout reaction mechanism
Direct reactions that knock out or transfer a single nucleon from a nucleus are a unique and powerful tool to probe singleparticle dominated states and nucleon correlations of exotic  [50,51]. The discrepancy between experimental and predicted cross sections has been quantified by a so-called reduction factor R s , defined as the experimentalto-theoretical cross section ratio, where the theoretical model includes both a reaction mechanism part (single-particle cross sections) and nuclear structure part (calculated SFs) as described in Sect. 2. The obtained R s from knockout reactions are around 0.6-0.7 for stable nuclei, approximately consistent with those obtained with benchmark (e, e N ) reactions. For the loosely bound states, eikonal-model predictions together with shell-model spectroscopic factors are close to experimental cross sections, resulting in R s around unity. While for the removal of a deeply bound nucleon, the model overestimates the experimental cross section by a factor 3-5, which has been first interpreted as a possible indication of missing correlations for deeply-bound nucleons in shell model calculations [50]. On the other hand, only a weak or even no dependence on proton-neutron asymmetry has been observed in transfer reactions [52][53][54][55]. Although the conclusions are based on few experimental data, still this apparent inconsistent dependence of SFs on proton-neutron asymmetry calls for a deeper understanding on the reaction mechanism and correlations in nuclei. A rather exhaustive discussion can be found in a recent dedicated review [56]. This discrepancy in the behaviour of R s has been attributed to failings in the applied reaction models. The adopted eikonal reaction model in the light-ion induced knockout reaction mainly has two contributions in the calculated single-particle cross section: the stripping/inelastic breakup and the diffractive/elastic breakup. The target nucleus, whose internal structure is not treated explicitly in the calcula-tion, is excited in the former and remains in the ground state in the latter. These two nucleon-removal mechanisms have been identified experimentally from their different kinematics [57]. Depending on the separation energy of the removed nucleon [58], the elastic breakup could contribute 15-35% to the one-nucleon removal cross section. In general, the eikonal model is formulated based on three approximations [22,23,[59][60][61]: the sudden/adiabatic approximation, which ignores the internal motion of nucleons during the reaction, the eikonal approximation, which assumes the scattering nucleons follow a straight line, and the spectator-core approximation, which treats the knockout residue as a spectator that can interact at most elastically with the target. The spectator-core model was initially elucidated by Hussein and McVoy [22] to discuss the stripping process. The validation of the spectator-core assumption has been checked by Al-Khalili when the core itself is loosely bound [62], where less than 10% correction for the stripping cross sections was found. However, in non-direct reaction processes, such as multiple scattering inside of the projectile, the excitation and decay of the spectator core could influence the deeply bound nucleon-removal cross section with a composite target [63]. In particular, recent theoretical and experimental works suggest that there may exist core excitation contributions that could significantly reduce the deeply bound nucleon-removal cross sections [63,64]. In fact, intranuclear cascade calculations, which include these effects but no microscopic information on the detailed nuclear structure, support the important role of evaporation in the cross sections for deeply bound nucleon removal from very asymmetric nuclei, i.e. leading to a fragile residual nucleus of low separation energy [63]. Still, a fully quantum mechanical and microscopic description is missing.
In addition to the knockout cross sections, the momentum of the knockout residue is often measured experimentally and is an observable which offers insights into the underlying reaction mechanism. Concerning parallel momentum distri-  [65] butions (PMDs) from the knockout involving deeply-bound nucleons, asymmetric PMDs have been observed by several experiments [66][67][68], i.e., a long low-momentum tail [66][67][68] and an abrupt high-momentum cutoff [68], in contrast to the symmetric PMDs predicted by the eikonal model. The origin of the observed asymmetric PMDs is still not clear. The low-momentum tail in some limited-asymmetry-tail cases can be fairly well described by the fully quantum mechanical coupled discretized continuum channel (CDCC) calculations including high-order effects for the elastic breakup component [60,69]. Moreover, the advanced transfer to continuum method developed by Bonaccorso and collaborators [70,71], and based on a target-nucleon optical potential, can reasonably describe the recent very asymmetric PMD of ( 14 O, 13 O) [68].
Proton induced quasi-free knockout ( p, pN ) is a complementary tool to intermediate-energy heavy-ion induced one-nucleon removal reactions. The ( p, pN ) reaction mechanism is simpler as the proton can be treated as a structureless probe, while it is expected that it probes a different part of the radial wave function compared to heavy-ioninduced nucleon-removal reactions due to a different absorption. Notably, the recently extracted SFs for the oxygen isotopic chain from quasi-free knockout reactions also show a very weak dependence on proton-neutron asymmetry [72,73] as opposed to the findings in heavy-ion induced knockout reactions. Later analysis with different theoretical methods show a consistent weak dependence [74,75]. Similarly to the heavy-ion induced knockout reactions, the momentum distributions from ( p, pN ) reactions are also sensitive to the orbital angular momentum of the ejected nucleon [65,[76][77][78]. A recent DWIA (distorted-wave impulse approximation) calculation on the 14 O( p, pN ) reaction also predicted large asymmetric PMDs [65], as shown in Fig. 3. The asymmetric PMDs were investigated in detail by the authors. The lowmomentum tail is caused by the attractive potential between the outgoing nucleons and the residue while the cutoff at high momentum side is due to the phase volume effects in the energy and momentum conservation [65], as in the case with heavy-ion targets. Experimental verification of PMD from ( p, pN ) reactions with good resolution is still necessary which will provide a test for a reliable reaction model for SF study. In addition, exclusive measurements of the knockedout proton or neutron and the recoil proton will also provide important information to clarify the relevant processes in the quasi-free knockout mechanism [79,80]. Given the non-observable nature of the spectroscopic and R s factors, the question of model dependence arises. In the next section we present the simpler deuteron electro-disintegration reaction as a proxy to analyse the effects of model dependence.

The electro-disintegration of the deuteron
We discussed above the ambiguities in the modeling of the breakup reaction mechanism and illustrated the typical discrepancies between data and predictions. Beyond the approximations to model the reaction mechanism, the interpretation of cross sections in terms of non-observable quantities such as spectroscopic factors or shell occupancies depends on the theoretical framework. We focus here on the historical example of the electro-disintegration of the deuteron to highlight the essence of ambiguities related to the interpretation of direct-reaction data.
The deuteron has been used as a benchmark for NN interactions. Its wave function can be reduced to a one-dimension problem when expressed in the center of mass coordinates with r the proton-neutron distance, u(r ) and w(r ) are the S-wave and D-wave radial wave functions, Y m are spherical harmonics and χ 1 the triplet spin wave function, as can be found in nuclear physics textbooks. Since the D-state admixture is small, the (observable) quadrupole moment of the deuteron is mostly proportional to the cross term uwr 2 dr, while the D-state probability is proportional to r 2 w 2 dr and is therefore sensitive to the interior of the wave function, as the centrifugal barrier constricts the D-wave component to the nuclear interior as compared to the S-wave component. It was soon understood. that the D-state probability is not an observable and depends quite strongly on the framework used to describe the deuteron and the cutoff scale in that framework (e.g. see Refs. [81,82]). To summarize in few words, unitary transformations, which leave observable quantities unchanged, modify the wave function of the deuteron in terms of S and D wave components, making them model dependent. The same conclusions apply in principle to the above-discussed orbital occupancies of the halo neutrons in 11 Li and any other nucleus, as well as other observables such as momentum distributions [83,84]. Such non-observable quantities are nevertheless important in our description of nuclear structure, but the model dependence has to be kept in mind when results of different studies are being compared. Such ambiguities also apply to the interpretation of reaction observables. The electro-disintegration of the deuteron e + d → e + p + n can be seen as the simplest breakup reaction sensitive to the internal structure of a nucleus: the involved bound nuclei are reduced to the deuteron itself in the entrance channel and the electro-disintegration process relies on the electromagnetic interaction.
First studies of electro-disintegration of the deuteron were performed at rather low incident energy and low momentum transfer, below 300 MeV/c [85][86][87]. At such low momentum transfer, the interpretation of the data required to take into account final state interactions (FSI) [88,89] and mesonexchange currents (MEC). Later measurements at Saclay [90,91], with a momentum transfer between 300 and 700 MeV/c, and at CEBAF [92] allowed to minimize these effects of FSI and EMC as the momentum transfer Q 2 of the virtual photon was increased [93]. Our understanding of the structure of the deuteron, and the high momentum component of the wave function, started from assuming a well controlled reaction mechanism and the plane wave approximation. The additional precision data in different kinematic regions (incident energy, momentum and energy transfer) shifted the focus of research from structure only to both the structure of the deuteron and the electron-scattering reaction mechanism (see [94] for a brief overview).
More recently, the electrodisintegration of the deuteron was used as a controlled laboratory to study how the factorization between reaction and structure components of the cross section is modified through renormalization-group transformations [95]. In particular, the scale dependence of the factorisation between the reaction and structure components of the cross section is well illustrated by a first-order analytical calculation. Assuming a unitary transform operator U (λ), as discussed in Section 2, any operator is evolved as where U is the smooth part of the unitary transform, residual of the identity operator I and defined as U = I + U . With the above, the evolved matrix element that enters the cross section is given by From Eq. 14, it is clear that the unitary transform does not change the observable, while the evolution of the operator modifies the reaction amplitude (first term of the right hand side) by δÔ (second term) which is reshuffled into a modification of the wave functions from modified initial and final state interactions, δ|Ψ λ i and δ|Ψ λ f (third and fourth terms), respectively.
While the results further demonstrate that scale dependence needs to be taken into account in low-energy nuclear physics, the RG methods show promise to be applicable to breakup reactions in the mid-term future. With the development of modern radioactive-ion beam facilities worldwide, allowing the production of intense beams at various energies and the impressive development of nuclear structure theory over the past few decades, it is foreseeable that the interest in the details of nuclear reaction mechanisms as a tool to extract nuclear structure will grow and develop more and more as a topic of study, as was the case for electron-induced deuteron studies over more than 50 years.

Conclusion
In this article, we have reviewed the sources of ambiguities in the interpretation of breakup reaction cross sections through two extensively-studied nuclear systems, the twoneutron halo 11 Li and the stable deuteron. These ambiguities are associated to the description of the reaction mechanisms and to the inherently non-observable character of certain quantities, shell occupancies in the above examples, used to interpret the data.
The exciting case of 11 Li was used to illustrate ambiguities related to the reaction mechanism. Typical variations of 10% on the s, p and d wave components for the two halo neutrons are encountered.
The electro-disintegration of the deuteron, despite its apparent simplicity, has been used as a benchmark to illustrate how that the factorization of the reaction mechanism and structure depends on the scheme and momentum scale of the theoretical framework used for the description of the data [95]. Our interpretation of reaction observables in terms of nuclear-structure non-observable quantities, such as the deuteron d-state probability, is therefore scale and scheme dependent as well. For the deuteron d-state probability this ambiguity in the description results in a range of values from 4 to 8% depending on the chosen NN interaction and its highmomentum behaviour. Similar scale-dependence on shelloccupancies and spectroscopic factors exist on nucleonremoval reactions with stable and unstable nuclei, and can lead to ambiguities if the framework and momentum scale of the model are not established.
The selection of the scale of the model is a formidable challenge, since so far the description of the reaction observ-ables has implied different descriptions with different frameworks and momentum scales for the reaction mechanism and the structure of the involved nuclei. Despite these inconsistencies, the standard approach to the description of nuclear breakup reactions (as described in Sect. 2) has provided remarkable advancements in our understanding of nuclear structure and reactions. There are still discrepancies between the experiment and predictions, in particular in the last decades with radioactive nuclei. These discrepancies, e.g. the isospin dependence of the nucleon removal cross sections, are being investigated to reduce the ambiguities, i.e., here a potential systematic bias in the reaction modeling, in the interpretation of breakup reaction observables.
Finally, we conclude on the words of C. Bertulani and M. Hussein that triggered this contribution: due to the limited experimental information about very exotic nuclei, the interpretation of recent experimental studies often remain ambiguous. The ambiguities can be reduced by confronting nuclear-structure interpretations from different measurements, i.e., different observables. The modeling of nuclear breakup reactions still leads to uncertainties in the interpretation of reaction observables in terms of nuclear structure. In the case of deeply-bound nucleon removal, the interpretation in terms of spectroscopic factors can vary up to a factor of 4 depending on the reaction and models used. Still, another ambiguity arises from the scale and scheme of the theoretical framework when not well defined. A consistent treatment of reaction and structure is therefore a high-priority target of upcoming reaction-theory developments in order to eliminate this ambiguity.
Acknowledgements The authors are thankful to J. Casal for providing the integrated 11 Li wavefunction used for Fig. 1 and to Y. Kubota for providing Table 1. Support from the Alexander von Humboldt foundation is acknowledged. This work has been performed in the context of the Grant SFB1245 from the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation).
Funding Information Open Access funding enabled and organized by Projekt DEAL.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: This work partly relies on published data. No new data are presented.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/.