Eﬀective ﬁeld theory analysis of the Coulomb breakup of the one-neutron halo nucleus 19 C

We analyse the Coulomb breakup of 19 C measured at 67 A MeV at RIKEN. We use the Coulomb-Corrected Eikonal (CCE) approximation to model the reaction and describe the one-neutron halo nucleus 19 C within Halo Eﬀective Field Theory (EFT). At leading order we obtain a fair reproduction of the measured cross section as a function of energy and angle. The description is insensitive to the choice of optical potential, as long as it accurately represents the size of 18 C. It is also insensitive to the interior of the 19 C wave function. Comparison between theory and experiment thus enables us to infer asymptotic properties of the ground state of 19 C: these data put constraints on the one-neutron separation energy of this nucleus and, for a given binding energy, can be used to extract an asymptotic normalisation coeﬃcient (ANC). These results are conﬁrmed by CCE calculations employing next-to-leading order Halo EFT descriptions of 19 C: at this order the results for the Coulomb breakup cross section are completely insensitive to the choice of the regulator. Accordingly, this reaction can be used to constrain the one-neutron separation energy and ANC of 19 C.


Introduction
Measurements of nuclear reactions along several isotopic chains show that the neutron distribution becomes extended as the neutron dripline is approached [1,2].This has led to the identification of "neutron halos": situations where a significant fraction of the neutron probability distribution resides in the classically forbidden region [3].Up to Z = 6 we already have examples of fourneutron halos, e.g., 8 He, two-neutron halos, e.g., 22 C, 19 B, and one-neutron halos, e.g., 11 Be, 19 C.This last nucleus demonstrates the striking features of a one-neutron halo.Following the dissociation of the halo neutron from the 18 C core, the momentum distribution of either of these fragments is narrow [4][5][6][7], as one would expect from a spatially extended system.Moreover, the breakup cross section of this fragile structure is large [8][9][10].This is particularly true on a heavy target such as Pb, for which the reaction is strongly Coulomb dominated.In that case an enhanced E1 strength is observed at low core-neutron relative energy, which is sometimes called the "pygmy dipole resonance".It is perhaps counterintuitive that properties of the neutron distribution can be probed through an electromagnetic observable, but this significant low-energy E1 strength is a consequence of the extended neutron distribution dragging the center-ofmass of the halo away from the center-of-charge.It is thus related to the halo physics that yields a significant isotope shift in these systems-and this Effective field theory analysis of the Coulomb breakup of the one-neutron halo nucleus relation can be formalised through the non-energy-weighted sum rule.For swave one-neutron halos this physics is "universal" in the sense that it depends only on the one-neutron separation energy and the Asymptotic Normalisation Coefficient (ANC) of the ground-state wave function.
The Coulomb dissociation of 19 C was measured by Nakamura et al. at RIKEN at 67A MeV already in the last millenium [8,9].The large E1 strength below 1 MeV outgoing relative energy of the 18 C-neutron system indicates the presence of a neutron halo.Comparison with models of the reaction implied that this is an s-wave halo, with a one-neutron separation energy S n = 530 ± 130 keV.
Halo Effective Field Theory (Halo EFT) provides a systematic way to analyse the Coulomb dissociation of one-neutron halos.(For a general introduction to Halo EFT and a review of the method's status as it stood in 2017, see Ref. [11].)Halo EFT expands the amplitude for the nuclear reaction in powers of the expansion parameter R core /R halo , where, in this case, R core is the size of 18 C, which amounts to approximately 2.5 fm, and R halo is the size of the neutron halo in 19 C, estimated to be about 6.5 fm.The calculation of Coulomb dissociation in Halo EFT confirms that the amplitude is universal at leading order, depending only on S n and the charge-to-mass ratio of the target [11][12][13].At next-to-leading order the asymptotic normalisation coefficient of the halo affects the amplitude.But, once S n and the ANC are fixed the amplitude is predicted-at least for s-wave halos-up to errors of order Rcore R halo 3 in the Halo EFT expansion.In Ref. [13] the photodissociation of 19 C was computed in Halo EFT and the equivalent photon approximation was used to convert the photodissociation cross section of 19 C into a Coulomb-breakup cross section.Acharya and Phillips extracted the value S n = 575 ± 55(stat.) ± 20(EFT) MeV from the low-energy (E < 1 MeV) and small-angle (θ < 2.2 • ) portion of the data from Ref. [8].
The reaction-theory employed in Ref. [13] was quite rudimentary.In this work we couple an EFT description of the 19 C bound state to a more advanced treatment of the reaction on the 208 Pb target that uses the Coulomb-Corrected Eikonal approximation (CCE) [14][15][16].This approximation corrects the erroneous treatment of the Coulomb interaction within the usual eikonal description of breakup reactions.It enables a computation of breakup cross sections at intermediate beam energies on both light and heavy targets that attains excellent agreement with fully dynamical reaction models while also retaining the simplicity and numerical efficiency of the usual eikonal approximation [15].Within this implementation of the CCE, the 19 C bound state, and the 18 C-neutron continuum, are described within Halo EFT, viz.using a set of 18 C-neutron potentials of Gaussian shape.In a leading-order calculation the depth of the Gaussian is adjusted to reproduce a particular S n .In a nextto-leading-order calculation, an additional term is added to the potential, and its parameter is adjusted to produce a specific ANC.Performing the calculation for a range of Gaussian widths checks whether the breakup cross section is insensitive to details of the potential.This imitates the strategy successfully employed for 11 Be and 15 C reactions on various targets in Refs.[17][18][19][20].Coupling a reliable model of the reaction to a Halo-EFT description of the nucleus provides a detailed account of the reaction mechanism, while enabling us to study very systematically the influence that the halo nucleus' structure has on the reaction cross sections.
The calculations described in this paper were initially performed as part of a week-long set of exercises at the TALENT school "Effective Field Theories in Light Nuclei: from Structure to Reactions" that took place at the Mainz Institute for Theoretical Physics in July-August 2022 [21].Students at the school (the majority of the authors in this paper) tuned the Gaussian potentials to reproduce specific scattering and bound-state parameters for the 19 C system.They then ran CCE calculations, predicted the Coulomb dissociation cross sections, and compared the result with data.The following sections describe their work and its outcomes, as follows.In Sec.2.2 we provide a brief summary of the reaction model and its implementation within the CCE.Section 3 lays out the leading-order (LO) calculation, presenting results for the 19 C system for Gaussians of widths ranging from 0.5 to 2.5 fm.These potentials are then used, together with the CCE, to predict the Coulomb breakup cross section.We find that the cross section scales with the square of the 19 C ANC, demonstrating that the reaction is almost exclusively peripheral.In Sec. 4 we explore the sensitivity of the results to the optical potentials chosen for the 18 C-208 Pb and neutron-208 Pb systems.In Sec. 5 we confirm that the cross section is insensitive to the interior of the 18 C-neutron wave function by performing a NLO calculation and finding (almost) the same result irrespective of the width of the Gaussian employed.In Sec. 6 we return to the LO potentials and vary the binding energy, in order to check the confidence interval given by Acharya and Phillips in Ref. [13].Finally, in Sec.7 we offer some conclusions and point out some interesting aspects of the EFT description of these Coulomb-dissociation data that, we believe, can motivate further theoretical and experimental studies of 19 C.

Reaction model 2.1 Three-body model of Coulomb breakup
To describe the breakup of 19 C on 208 Pb, we consider the usual three-body model of the reaction [16].The projectile P is seen as a two-body structure: a halo neutron (n of mass m n and charge nil) loosely-bound to a 18 C core assumed to be in its 0 + ground state (c of mass m c and charge Z c e).This two-body structure is described by the effective Hamiltonian where r is the c-n relative coordinate, µ = m c m n /(m c + m n ) is the c-n reduced mass, and V cn is an effective potential that describes the c-n interaction.As Effective field theory analysis of the Coulomb breakup of the one-neutron halo nucleus discussed in Secs.3.1 and 5, we consider Halo-EFT interactions up to NLO [11,22].The eigenstates φ of H 0 describe the different states of the projectile.The negative-energy eigenstates correspond to the c-n bound states.They are discrete and, in addition to quantum numbers of the c-n orbital angular momentum l, the total angular momentum j, and its projection m, they are identified by the number of nodes in their radial wave function n r .Asymptotically, the radial part of these bound-state wave functions behaves as where k l is a modified spherical Bessel function of the second kind, κ nrlj is related to the eigenenergy of the state E nr lj = − 2 κ 2 nrlj /2µ, and C nrlj is the asymptotic normalisation coefficient (ANC) associated with that bound state.The positive-energy eigenstates describe the c-n continuum part of the projectile spectrum, viz. the broken up projectile.As such they are identified by their c-n relative energy E, in addition to the quantum number defining the partial wave l, j, and m.
The lead target T is assumed to be a structureless cluster of mass m T and charge Z T e.Its interaction with the projectile components c and n is described by the optical potentials V cT and V nT , respectively.These potentials are found in the literature as explained in Sec. 4.
Within this three-body model of the collision, studying the P -T collision reduces to solving the following Schrödinger equation with the three-body Hamiltonian where R is the coordinate of the projectile center of mass relative to the target, µ P T = m P m T /(m P + m T ) is the P -T reduced mass-with m P = m c + m nand R cT , resp.R nT , are the c-T , resp.n-T , relative coordinates.The total energy E T in Eq. ( 3) is related to the initial P -T kinetic energy and the eigenenergy of the projectile in its initial ground state φ nr0l0j0m0 through where K 0 is the wave vector of the incoming P -T relative motion; that direction defines the Z axis of the system of coordinates.The Schrödinger equation (3) has to be solved with the incoming condition that the projectile, in its ground state, is impinging on the target.Accordingly, the three-body wave function behaves as Various numerical techniques, based on different approximations, have been developed to solve this equation, see Ref. [16] for a recent review.For this study, we have considered the CCE [14][15][16], which is very efficient at the intermediate beam energy considered here.

Coulomb-Corrected Eikonal approximation
At sufficiently high energy, the eikonal approximation is quite reliable to describe the P -T collision [16,23].Within that approximation, the three-body wave function after the collision reads where the eikonal phase is given by with v = K 0 /µ P T the P -T initial velocity.
Being based on an adiabatic description of the reaction, the usual eikonal approximation is valid only for reactions that take place over a short time, viz. that are dominated by the short-ranged nuclear interaction.When the Coulomb interaction is non-negligible, such as for the lead target considered in this study, the breakup cross section inferred from the expression (7) diverges [14,15].Margueron et al. have developed a correction, that efficiently solves that divergence [14].The main idea is to use the first order of the perturbation theory, which accounts for the projectile dynamics, to correct for the erroneous treatment of the Coulomb interaction at the eikonal approximation.In the CCE, the Coulomb contribution to the eikonal phase is replaced, at first order, by its corresponding perturbative estimate [15]: where χ N and χ C are, respectively, the nuclear and Coulomb contributions to the eikonal phase χ (8), and where the first-order phase reads with η = Z c Z T e 2 /(4πǫ 0 v) the P -T Sommerfeld parameter and ω = E − E nr0l0j0 the energy difference between the final continuum state and the initial bound state of the projectile.Effective field theory analysis of the Coulomb breakup of the one-neutron halo nucleus By accounting for the projectile dynamics in the first-order treatment of the Coulomb interaction, this correction solves the aforementioned divergence issue.Moreover the expression (9) enables us to account also for the nuclear part of the P -T interaction at all orders, its interference with the Coulomb force, and, although only in an approximate way, for higher-order Coulomb effects.This CCE leads to breakup cross sections in excellent agreement with fully dynamical models [15].It is thus well suited to describe breakup reactions at intermediate energies on both light and heavy targets, while exhibiting the simplicity and numerical efficiency of a usual eikonal code.In this study we consider the CCE to compute the breakup cross section of 19 C impinging on 208 Pb at 67A MeV and compare these theoretical results with the data of Ref. [8].
3 Leading-order calculation 3.1 Leading-order description of 19

C
The one-neutron halo nucleus 19 C has a 1 2 + ground state that lies slightly more than half an MeV below the one-neutron separation threshold (S n = 0.58 ± 0.09 MeV [24]).Various experiments have confirmed the one-neutron halo structure of that state [4][5][6][7][8][9][10].Therefore it is usually described as a 18 C in its 0 + ground state to which a valence neutron is loosely bound in the s 1/2 partial wave.As mentioned earlier, we consider in the present study a Halo-EFT description of 19 C, assuming the halo neutron sits in a 0s 1/2 bound state, i.e., with n r = 0 nodes in the radial wave function.
At leading order the 18 C-n interaction is described by a Gaussian potential: where the standard deviation of the Gaussian, σ, acts as a regulator.In the limit σ → 0 this becomes a three-dimensional δ-function.We consider σ = 0.5, 1, 1.5, 2, 2.5 fm.For each σ, the potential strength C 0 is adjusted to produce a 0s 1/2 19 C state that is bound by 0.58 MeV with respect to the 18 C-neutron threshold.The corresponding values of the C 0 s and the ANCs C 0s1/2 predicted for 19 C are given in Table 1.
Figure 1 shows the reduced radial wave functions obtained within this LO Halo-EFT model of 19 C (a) normalised to unity and (b) divided by their ANC.Panel (b) shows that the wave functions have the same asymptotic behaviour (thin black dashed line), up to a multiplicative constant-as should be the case given the way in which they were constructed.It also shows they differ markedly at short range, viz.for r ∼ < 3 fm.This provides a straightforward way to test if the reaction is peripheral: if it is, the reaction cross section will scale as the square of the ANC.If the reaction is sensitive to the short-range piece of the 18 C-n wave functions then that scaling will break down.

Coulomb breakup cross sections at LO
We now execute the CCE code with the LO 18 C-n potentials of Sec.3.1 1 .As explained in Sec. 2, the interactions between the projectile constituents and the target are simulated by optical potentials selected from the literature.The reasons for this selection, and its effect on our calculations, will be discussed in Sec. 4. Figure 2 gives the direct CCE results-viz.without data-for the five values of the Gaussian range σ considered in Sec.3.1.The breakup cross section plotted as a function of the 18 C-n relative energy after dissociation is shown in Fig. 2(a), whereas Fig. 2(b) displays it as a function of the scattering angle of the 18 C-n centre of mass for a continuum energy 0 ≤ E ≤ 0.5 MeV.For σ = 1.5 fm (blue dash-dotted lines), the contributions to the cross section from s, p, and d waves in the 18 C-n continuum are shown separately.It is immediately clear that the reaction is dominated by an E1 transition from the s ground state to the p continuum, as expected for the part of the cross section mediated by a single E1 photon exchange between the 208 Pb nucleus and the 19 C projectile.For our LO calculation we take the p-wave phase shifts to be 1 An open-access version of the fortran code "Chaconne" has been developed for the TALENT school "Effective Field Theories in Light Nuclei: from Structure to Reactions" [21].The program, a user's manual, test input file, and the corresponding output, can be downloaded from the school website, see the documents attached to the lecture on nuclear reaction theory in the third week of the school https://indico.mitp.uni-mainz.de/event/279/timetable/#20220808.Effective field theory analysis of the Coulomb breakup of the one-neutron halo nucleus 0, because there is no known state with negative parity at low energy.This makes the overall result rather simple, cf.Eq. ( 15) of Ref. [13].However, the presence of nuclear interactions between 18 C and 208 Pb and between the neutron and 208 Pb, as well as the possibility of multiple photon exchanges, produce noticeable contributions to the cross section from s and d waves in the 18 C-neutron continuum.Both s-and d-wave contributions become a larger fraction of the breakup cross section as the angle increases, although the d-wave piece stays a factor of a few below the s-wave one throughout the angular range of interest here.The s-wave effect is more important at lower relative energy, with the d-wave one growing as energy increases.This is a significant finding, because the CCE is nearly as simple mathematically as the first-order E1 treatment carried out in Ref. [13], but it allows us to quantify the nuclear contribution to breakup, its interference with the Coulomb force, and other quantal interferences seen in the oscillatory pattern of the angular distribution [15].
While the way that the reaction mechanism populates different partial waves in the continuum is interesting, the key finding from Fig. 2 is that population of anything other than the continuum p is small enough that the total cross section scales (nearly) perfectly with C 2 0s1/2 .This shows that reaction is purely peripheral, since it demonstrates that the breakup does not probe the short-range physics of the projectile.
Because the cross sections scale with the C 2 0s1/2 and because the cross section exhibits little sensitivity to the choice of the nuclear part of the optical potential (see Sec. 4), we can infer an ANC by fitting the calculations to the data.To avoid the regions where the nuclear interaction plays a role and where the d waves, which are not well constrained, might affect the calculation, we focus on the forward-angle region-viz.θ < 2 • -of the angular distribution, which is restricted to small 18 C-n relative energies-viz.0 ≤ E ≤ 0.5 MeV.As seen in Fig. 2, that region is dominated by the p-wave contribution.
In order to extract a reliable value for C 0s1/2 from data it is necessary to account for the experimental resolution.This is done by folding the theoretical cross sections with the resolution provided in the experimental paper [8].After folding, we scale the calculations to the data.Minimizing the χ 2 with respect to the scaling factor enables us to infer the ANC: This value, and its uncertainty, are independent of the value of σ chosen for the 18 C-n potential (11), confirming the independence of the calculations to the short-range physics, and hence the accuracy of the method.

(In)Sensitivity of the calculations to the nuclear optical potentials
To test the sensitivity of the calculations to the choice of optical potentials, they have been repeated with different interactions found in the literature.
The results of the previous section followed Typel and Shyam in Ref. [25] and chose for the 18 C-208 Pb interaction a potential developed by Buenerd et al. to reproduce the elastic scattering of 13 C off 208 Pb at 390 MeV (30A MeV).The fact that this is a rather different energy than the one employed in Refs.[8,9] is ignored, but the radii are scaled to the actual size of the core of the projectile.
For the n-208 Pb interaction, in the previous section we also followed Typel and Shyam and use the Becchetti and Greenlees global optical potential (BG) [26].
For a second 18 C-208 Pb optical potential, we use the one considered by Typel and Shyam in Ref. [25] to simulate the interaction between 11 Be and 208 Pb at 70A MeV.That potential is based on an α-208 Pb potential developed by Bonin et al. to reproduce that elastic scattering at 288 MeV (72A MeV), from which we rescale the radii to account for the size of the nucleus.As a second n-208 Pb potential choice, we opt for the Koning-Delaroche global optical potential (KD) [27].
The results of these different calculations obtained with the LO 18 C-n potential with σ = 1.5 fm are shown in Fig. 3.It is clear that these choices have very limited influence on this Coulomb-dominated reaction.We note that if we strictly follow Typel and Shyam's procedure from Ref. [25] and do not adjust the radius of the projectile carbon nucleus then we obtain a higher cross section than is seen here.The size of the core is an important parameter in these calculations, even if the results are not sensitive to the functional form of the potential's radial dependence.

NLO calculations
At NLO the 18 C-neutron potential takes the form: Effective field theory analysis of the Coulomb breakup of the one-neutron halo nucleus where the parameter C0 is not necessarily-indeed not usually-the same as the parameter C 0 .This time we consider potentials with different σs and, in each case adjust them to produce S n = 0.58 MeV and C 0s1/2 = 0.81 fm −1/2 .We achieve this for σ = 1.0, 1.5, and 2.5 fm.The resulting cross sections predicted by the CCE are now completely independent of σ, see Fig. 4, where the theoretical cross sections have been folded with the experimental resolution [8].Moreover, despite the fact that we fit the ANC only to the forward-angle region of the angular distribution limited to E ≤ 0.5 MeV, we find that all calculations match the experimental energy distribution over nearly the entire experimental energy range, viz.out to E = 4 MeV.The excellent agreement with experiment, and the insensitivity of the NLO results to the regulator σ, confirm the value of the ANC we inferred from the data using our LO calculations.This also shows that it is not necessary to go beyond NLO to explain the main features of the data.
Note that no NLO potential could be found for σ = 0.5 fm, i.e, we could not find parameters to fit simultaneously the binding energy and the ANC inferred from the data.This is a realisation of the Wigner bound [28][29][30] for this system: for any C 2 0s1/2 larger than 2(2µS n / 2 ) 1/2 the integral of the asymptotic wave function from zero to infinity is larger than one.It follows that for small enough σ it is simply impossible to produce a normalisable wave function with this ANC 2 .
Significant discrepancies between theory and experiment appear at about E ≈ 1.3 MeV and 2.8 MeV in the energy distribution.At those energies the data seem to be notably larger than our calculations.The large error bars on the experimental data leave open the possibility that these are statistical fluctuations and not due to an effect of final-state interactions 2 .However In both cases the calculations have been folded by the experimental resolution [8].
these deviations could hint at the presence of resonances in the 19 C system at these energies.Refs.[7,10] suggest the existence of a 5 2 + resonance at either E = 1.42 (10) MeV [7] or E = 1.46 (10) MeV [10].This state might have a dominant single-particle structure with a 18 C core in its 0 + ground state and neutron in a d 5/2 resonance and could significantly affect the breakup cross section [31][32][33].Within the usual Halo-EFT power counting, it would therefore enter beyond NLO.Following what has been done in Refs.[17,33] an extension of this work could study this possibility.
At large angles in the angular distribution, we also observe that the calculations slightly overestimate the data.This is a region where the nuclear interaction plays a more significant role, see Fig. 3(b), and hence is subject to caution because this difference might be related to the choice of optical potentials.It could also come from the 18 C-n final-state interaction in the d wave, which is not constrained at NLO; see Fig. 2(b).
The good agreement with experiment suggests that, in absence of more precise measurements, a Halo-EFT description at NLO is both necessary and sufficient to describe most of the breakup data.

Sensitivity to the binding energy
The binding energy quoted in the most recent atomic mass database [24] exhibits a rather large uncertainty: S n = 0.58 ± 0.09 MeV.To gauge the influence of this observable on the calculations, we repeat breakup calculations using LO Halo-EFT 18 C-n potentials fitted to the lower (S n = 0.49 MeV) and upper (S n = 0.67 MeV) end of this 68% confidence interval.We consider σ = 1.5 fm for this test.
Fig. 5 displays the results folded with the experimental resolution and fitted to the data by rescaling the CCE calculation.The ANCs hence obtained differ significantly from the one quoted above: C 0s1/2 (S n = 0.49 MeV) = 0.62 ± 0.02 fm −1/2 and C 0s1/2 (S n = 0.67 MeV) = 1.02 ± 0.03 fm −1/2 .This shows that Effective field theory analysis of the Coulomb breakup of the one-neutron halo nucleus the ANC and binding energy are strongly correlated, as one would expect from the LO Halo EFT relation [11,34] We note that, while the strict scaling fo the ANC-squared at LO is with √ S n , the higher-order terms indicated in Eq. ( 14) are ultimately quite important in the case of 19 C. The ANC-squareds inferred for different binding energies scale markedly more strongly with S n than √ S n .Although the prediction with the lowest binding energy seems to better reproduce the angular distribution throughout the entire experimental angular range [see the green dashed line in Fig. 5(b)], the corresponding energy distribution does not fit the data at E > 0.5 MeV [see Fig. 5(a)].Using a higher binding energy leads to less good agreement with the data in both observables (red solid lines in Fig. 5).This suggests that the actual binding energy is probably close to the central value we have considered up to Sec. 5, viz.S n = 0.58 MeV.It also indicates that, with more thorough uncertainty quantification, including an estimate of the impact of higher-order effects in the EFT [35][36][37], and the uncertainty due to the choice of optical potential, analysis of these data could yield a new, more precise, value, for the 19 C one-neutron separation energy.

Conclusion and Needed Future Work
Many experiments have shown that 19 C exhibits a clear one-neutron halo structure in its 1 2 + ground state [4][5][6][7][8][9][10].However, in contrast to the well-studied cases of 11 Be and 15 C, there is still much to learn about this nucleus, including its one-neutron separation energy S n .In this paper, we present a new analysis of the Coulomb breakup of 19 C on 208 Pb at 67A MeV, which has been measured at RIKEN [8].To this aim, we have used a Halo-EFT description of the projectile within the Coulomb Corrected Eikonal approximation (CCE), which has shown to provide reliable cross sections for this kind of reaction [15], while exhibiting a small numerical cost.
As expected these cross sections are strongly dominated by an E1 transition from the 0s 1/2 ground state of the nucleus towards its 18 C-n continuum.Being Coulomb dominated, they exhibit a minor dependence to the optical potentials used to simulate the nuclear interaction between the projectile constituents ( 18 C and n) and the 208 Pb target.
Using a LO description of 19 C, we have found out that the calculated cross sections are nearly proportional to the square of the ANC of the radial 18 C-n wave function C 0s1/2 .This clearly shows that the reaction is purely peripheral in the sense that it probes only the tail of the ground state wave function.It also indicates that an ANC for the actual nucleus can be inferred from the data.To reduce the uncertainty related to the choice of the optical potentials as well as to avoid the influence of the d-wave continuum, we select forwardangle data at low 18 C-n energy to scale our calculations to the experiment.The value of the ANC hence obtained, C 0s1/2 = 0.81 ± 0.02 fm −1/2 , is independent of the Halo-EFT regulator σ.NLO descriptions of 19 C fitted to reproduce both S n and that value of C 0s1/2 provide an excellent agreement with the data on nearly their entire energy and angular ranges, independently of the value of σ.
Additional tests have shown a strong dependence of the calculations to the binding energy of the nucleus, which, unfortunately is not well known experimentally.However, our tests show that a systematic analysis of Coulombbreakup data, e.g., through Bayesian methods [35][36][37], could provide a significant constraint on that structure observable.Thanks to its small computational cost and accurate description of the reaction process, the CCE would be the ideal reaction-dynamics treatment for such a future statistical analysis.
This theoretical study extends a series of analyses of reactions involving one-neutron halo nuclei, in which a Halo-EFT description of the exotic nucleus is coupled to realistic models of reactions [17][18][19][20].Our work confirms the validity of this approach for the Coulomb breakup of 19 C, and shows that crucial nuclear-structure information can be inferred from such a study.Unfortunately, the experimental uncertainty of the RIKEN data considered in this work [8] is too large to draw reliable conclusions on these structure observables.Accordingly, we advocate for new experiments with smaller uncertainties to pin down these values.Similar breakup data would help us constrain both the binding energy of 19 C and its ANC.Breakup data on a light target, viz. 12C or 9 Be, could help investigate the possible presence of single-neutron resonances in Effective field theory analysis of the Coulomb breakup of the one-neutron halo nucleus the continuum.Transfer measurements, such as 18 C(d,p) in inverse kinematics could help constrain the ANC of the ground state, especially if they are measured at low beam energy and forward angles [18].Knockout measurements with improved uncertainty compared to existing data [4][5][6] would also improve our understanding of this exotic nucleus [20].

Fig. 1
Fig. 1 Reduced radial wave functions of the 0s 1/2 18 C-n bound state (a) normalised to unity; (b) divided by their ANC C 0s1/2 for different values of σ, as indicated in the legend.

Fig. 2
Fig. 2 Breakup cross section of 19 C on 208 Pb at 67A MeV (a) plotted as a function of the 18 C-n relative energy E after dissociation, and (b) plotted as a function of the scattering angle of the 18 C-n centre of mass for energies 0 ≤ E ≤ 0.5 MeV.In both cases, s, p, and d components are shown separately for the σ = 1.5 fm case (blue dash-dotted lines).

Fig. 3
Fig. 3 Influence of the optical potential choice on the breakup cross section of 19 C on 208 Pb at 67A MeV; (a) energy distribution; (b) angular distribution.

Fig. 4
Fig.4NLO calculations of the breakup of19 C on 208 Pb at 67A MeV compared to the data of Ref.[8].NLO Halo-EFT18 C-n potentials are fitted to reproduce the binding energy and the ANC inferred from the comparison of the LO calculations to the data (angular distribution restricted to forward angles); (a) energy distribution and (b) angular distribution.In both cases the calculations have been folded by the experimental resolution[8].

Fig. 5
Fig.5Sensitivity of breakup calculation to the 18 C-n binding energy for19 C impinging on 208 Pb at 67A MeV.LO Halo-EFT 18 C-n potentials are fitted to reproduce three binding energies at the centre (0.58 MeV, blue dash-dotted line), lower bound (0.49 MeV, green dashed line) and higher bound (0.67 MeV, red solid line) of the experimental uncertainty range[24].The calculations have been scaled to the data as explained in Sec.5; (a) energy distribution; (b) angular distribution.In both cases the calculations have been folded with the experimental resolution[8].

Table 1
Strengths of the LO 18 C-n potentials for the different regulators σ considered in this study [see Eq. (11)].They have been fitted to reproduce the ground state energy at E 0s1/2 = −0.58MeV.The corresponding ANCs C 0s1/2 are listed as well.