Finite-size effects in heavy halo nuclei from effective field theory

Halo/Cluster Effective Field Theory describes halo/cluster nuclei in an expansion in the small ratio of the size of the core(s) to the size of the system. Even in the pointparticle limit, neutron-halo nuclei have a finite charge radius, because their center of mass does not coincide with their center of charge. This point-particle contribution decreases as 1/Ac, where Ac is the mass number of the core, and diminishes in importance compared to other effects, e.g., the size of the core to which the neutrons are bound. Here we propose that for heavy cores the EFT expansion should account for the small factors of 1/Ac. As a specific example, we discuss the implications of this organizational scheme for the inclusion of finite-size effects in expressions for the charge radii of halo nuclei. We show in particular that a short-range operator could be the dominant effect in the charge radius of one-neutron halos bound by a P-wave interaction. The pointparticle contribution remains the leading piece of the charge radius for one-proton halos, and so Halo EFT has more predictive power in that case.


Introduction
Theoretical models of many-body systems usually treat the constituent particles as having no internal structure. This point-particle approximation is also used in cluster models, e.g., for the description of halo systems [1][2][3][4], even though one might encounter situations for which the core is rather large. Finite-size effects are included a posteriori, but can become significant for certain observables. As an example, the total charge radius is usually calculated by simply adding in quadrature [5] the charge radii of the constituents and the a e-mail: christian.forssen@chalmers.se calculated point-particle radius; see e.g. the calculation of charge radii for neutron-rich helium isotopes in the Gamow Shell Model [6]. Instead, it would be useful to construct a framework in which finite-size effects can be included systematically.
Constituent-size effects can be accounted for in effective field theories (EFTs), where they appear through derivative interactions. For example, the nucleon charge radius (and, more generally, nucleon form factors) can be calculated [7] in Chiral Perturbation Theory (ChPT) in an expansion in powers of k π /M QCD , where k π ∼ 150 MeV is a momentum scale associated with the lightest carrier of the nuclear force, the pion, and M QCD ∼ 1 GeV is the characteristic mass scale of QCD. The relevant pion parameters are its mass and decay constant. Chiral EFT [8] is a generalization of ChPT to a typical nucleus, for which the binding energy per nucleon is B/A ∼ k 2 π /M QCD and the radius R ∼ A 1/3 /k π . The nuclear charge radius includes the sum of the nucleons' radii plus many-body effects generated by internucleon interactions and currents [9]. We would like to have a similar framework for clusterized systems.
Clusterized systems, with much smaller energies and larger radii, are additionally characterized by scales beyond the pion scales. These nuclei can be viewed as a collection of valence nucleons orbiting around either no core (fewnucleon systems), one core (halo nuclei) or many cores (cluster nuclei). The cores themselves frequently-but not always-have properties of typical nuclei. The generic existence of such systems can be understood as a consequence of a fine-tuning in QCD, which introduces a lighter momentum scale ℵ ∼ 30 MeV [10,11]. For such loosely bound systems we can devise EFTs that exploit the separation of scales without involving pions explicitly. In these EFTs one considers processes with typical momenta k lo , such that k lo k hi , 0123456789().: V,-vol where k hi k π is a high-momentum scale. One then develops an expansion for observables in powers of k lo /k hi . The very lightest nuclei are dilute systems with no core, where the dominant (two-and three-nucleon) interactions are S-wave. The corresponding EFT is Pionless EFT, for which power counting is relatively well understood [8].
Halo/Cluster EFT, here labeled HEFT, was proposed as an EFT for systems with one [12,13] or more [14] cores and valence nucleons [15,16]. (See Ref. [17] for a recent review.) HEFT power counting is a generalization of the power counting for Pionless EFT allowing for dominant interactions in waves with non-vanishing angular momentum and for a breakdown scale k hi estimated from the first excitation of the core and/or its size. In the first cases considered, 5,6 He [12,13,[18][19][20], there is an alpha-particle core, and the neutron-alpha (n-α) interaction is mostly of P-wave nature, generating a near-threshold 5 He resonance. The α-α interaction, in turn, is obtained [14] from α-α scattering and the lowest 8 Be state. HEFT has since been extended to heavier cores and to proton halo systems [17]. In most of these cases the high scale k hi in HEFT is associated with the size of the core, i.e. k hi ∼ 1/R c . This means HEFT is an expansion in powers of R c /R halo , where R halo is the unnaturally large size of the halo system.
The different contributions to the charge radius of a halo nucleus can then be organized in the HEFT expansion and the size of the effect due to the finite size of the core (and of the nucleon) estimated. We do that and thereby derive-for both S-and P-wave one-neutron halos-the charge-radius formula that is frequently used in nuclear theory. However, we also point out that there is, in principle, another expansion parameter present when HEFT is applied to systems with a relatively large number, A c , of core nucleons. To leading order in 1/A c the core is static, its recoil being small compared to nucleon recoil. Consequently the center of mass of a neutron halo coincides with its center of charge. Thus, whereas for light halos (e.g. 6 He) the difference between these two generates an important contribution to the charge radius [21], for heavier systems the corresponding effect goes to zero. Correctly assessing the impact of the finite size of the constituents on the charge radius requires keeping track of factors of 1/A c .
Moreover, the charge radii of halo nuclei are affected by a short-range operator, which is subleading in R c /R halo but leading in 1/A c . We show that for one-neutron halos bound by a P-wave interaction (e.g. the excited state of 11 Be) this effect may be as important as the long-distance contributions to the halo's charge radius that have been previously computed in HEFT [22]. Similar considerations also apply to the form factors of two-neutron halos such as those discussed in Ref. [23]. They are not, however, as pressing for proton halos, where a finite charge radius will be generated by the photon coupling to the valence proton(s) even if the core is infinitely heavy.
This exemplifies the importance of keeping track-to the extent possible-of factors of 1/A c in observables, rather than just counting powers of k lo . This is quite similar to the need to distinguish between relativistic corrections that are suppressed by powers of the inverse nucleon mass, 1/m N , and other corrections that only carry powers of 1/k hi [24][25][26]. The significant difference between these two scales produces a hierarchy between effects that scale with the same power of k lo . Only once that hierarchy has been identified can the power counting be formulated in an efficient manner.
We isolate the A c dependence that enters the charge radius through kinematic effects, i.e., because the nucleon-core mass ratio, m N /m c , is small. In contrast, we assume that all Lagrangian coefficients (LECs) scale with a power of k hi that is solely determined by the naive engineering dimension of the operator they multiply, i.e., we use naive dimensional analysis with respect to k hi and do not attach any additional A c dependence to the LECs. It is true that k hi is also A c dependent, because k hi will generically be of order the inverse core size, 1/R c , and R c can be taken to be ∝ A 1/3 c . But any additional accounting of the A c dependence of short-distance physics in HEFT beyond this would require a more microscopic understanding of the A c dependence of all the cofficients in the EFT. This could be achieved by matching HEFT to a microscopic calculation, but it could be argued that such matching goes beyond the EFT philosophy of writing down a theory that is independent of the short-distance physics. In contrast, the kinematic effects we identify here are universal, in the sense that they occur irrespective of the nature of the forces between the halo nucleons and the core(s).
Here we focus on the charge radius of single-neutron halo nuclei, where the point contribution is suppressed by 1/A c for the reason described above. But the presence of the heavy-core propagator is ubiquitous, so similar effects will affect other observables as well. For example, one expects the Born-Oppenheimer approximation to emerge in systems with multiple heavy cores and/or valence nucleons.
Our paper is structured as follows: In Sect. 2 we discuss the interplay between the various scales that are involved in an EFT for a halo system with a heavy core. We also discuss the low-energy scattering parameters for the nucleoncore system and introduce the charge radius in terms of the momentum expansion of the low-energy charge form factor. In Sect. 3 we derive the observable charge radius for Sand P-wave one-neutron-halo states. The power counting is exemplified by considering the charge radius for selected halo states. We summarize our findings in Sect. 4. An appendix discusses the corresponding results for proton halos, where considering factors of 1/A c does not lead to any change in the hierarchy of the various physical mechanisms that contribute to the charge radius.

Power counting
Once the relevant degrees of freedom are chosen, a model consists of a specific set of interactions among them. In contrast, with an EFT one considers the most general dynamics consistent with the known symmetries. It is crucial to organize the corresponding infinity of contributions to any observable according to their size ("power counting").
We are interested in a clusterized system where the size R halo ∼ 1/k lo of the system is sufficiently large that the constituents can be taken as point-like in a first approximation. This system might be probed with particles (photons, electrons, neutrinos, nucleons) of wavelength ∼ 1/k lo that cannot resolve the inner structure of the constituents. For simplicity we consider a few valence nucleons orbiting around a single core of radius R c consisting of A c 1 nucleons. The arguments below can be generalized straightforwardly to multiple-core systems.
The power counting of HEFT is based, like that of other EFTs, on the ratio of momentum scales, k lo /k hi 1. The high-momentum scale k hi is determined by physics not accounted for explicitly in HEFT. Since nuclei are bigger than nucleons we must have k hi 1/R N ∼ k π , with R N the size of a nucleon, which is generically set by pionic dynamics described by ChPT. But a more restrictive condition on k hi arises from the requirement that details of the core are not resolved: Adopting the standard rule for the scaling of the nuclear size with A c we then have k hi ∼ A −1/3 c k π , although we note that several of the cores considered up until now in HEFT are somewhat larger than this valley-of-stability lore indicates. If there exist low-lying excited states of the core, corresponding to a lower momentum scale than the inverse size of the core, the high-momentum scale needs to be adjusted accordingly, or else explicit degrees of freedom must be introduced for the low-lying excited states. One example of this is in the 7 Be( p, γ ) 8 B reaction, where the low-lying excited state of 7 Be must be included as an explicit degree of freedom in the HEFT if Eq. (1) is to prevail [27,28]. k −1 hi can therefore be expected to increase for heavier systems, rendering finite-size effects more important. Meanwhile, 1/A c decreases and will generically be smaller than k lo /k hi . It introduces an additional expansion parameter. Explicit factors of 1/A c 1 enter through the core mass, where the average nucleon mass m N 940 MeV and we have neglected the effect of nuclear binding for the kinematic purposes we have in mind here.
In our non-relativistic theory approximate Galilean invariance guarantees that the interactions depend on the mass of the particles only in trivial ways that can be scaled out of the theory. Therefore all explicit occurrences of 1/A c are associated with the kinematics of the two-particle system, and once again, Galilean invariance means that they must be encoded in the halo's reduced mass, where is the total nucleon-core mass, with f 1/(A c + 1) ∼ 1/A c . Specifically, m R has a fractional difference from m N of ≈ 1/A c , reflecting the extent to which the core is still dynamical in the (effective) one-body problem.
The large mismatch in masses evident in Eq.
(2) means that the pertinent energy scale for the valence nucleon is the one-nucleon separation energy which is much smaller than the binding energy of the core but much larger than the recoil energy of the core For energies of the order of the typical nucleon energy, E ∼ B s , nucleon recoil is a leading-order effect. Beyond leading order the ratio k lo /m N occurs only in (small) relativistic corrections. In contrast, core recoil is suppressed by a factor of 1/A c . Thus at leading order (LO) the core propagator is static, that is, Thus, for low-order calculations one can simplify the procedure by using a static core, and including recoil effects perturbatively as higher-order corrections. Relativistic corrections that scale like k lo /m c will be even smaller.
In addition to these kinematic factors of 1/A c , there may be dependence on A c coming through the interaction coefficients, or "low-energy constants" (LECs). As a trivial example, electromagnetic interactions add up constructively for protons and the corresponding LECs in general depend on the core charge Z c = A c − N c . It is not clear how to deal with this quantity a priori. In neutron halos, Z c can be significantly smaller than A c /2, but this is not necessarily so for proton halos. We will keep factors of Z c explicit and deal with them on a case-by-case basis.
Likewise, the LECs for strong interactions might in specific cases represent some constructive or destructive interference in the interactions of the valence nucleon with the core nucleons. One way to determine the A c dependence of these LECs is by matching HEFT to the ab initio solution of the same system with a more fundamental EFT [27][28][29], in a region where both EFTs are valid. Another way is to look at systematic trends in LECs fitted to data for different cores. In either case a manifestation of strong A c dependence would be a particularly large or small LEC value with respect to the expected power of k hi . Since there is no clear case at this point, below we limit ourselves to the kinematical factors arising from the core mass, although the counting of factors of 1/A c could be improved later if needed.

Nucleon-core scattering
EFTs incorporate from the start the coupling to the continuum, so that most calculations of halo structure, including form factors, are intimately connected with nucleon-core scattering. A discussion of nucleon-core interactions is therefore necessary for the calculation of form factors, and we briefly review previous work on the subject here.
First we consider a halo system where the dominant nucleon-core interaction is S-wave. In this case, an EFT where all forces are of contact type reduces to the effective range expansion (ERE) [8,30]. One can think of the scattering length a 0 as what characterizes the size of the halo system, and the effective range r 0 (and higher ERE parameters) as reflecting the breakdown scale, namely the size of the core. The large size of the halo system is manifest in a large scattering length, while higher effective-range parameters are assumed to have sizes set by 1/k hi : For an S-wave nucleon bound to the core with separation energy B s0 > 0, the nucleon-core T matrix has a pole at k = iγ 0 with The power counting for this system is almost identical to that of Pionless EFT for an S-wave bound state [8,30], which was adopted, for example, in the form-factor calculation of Ref. [22]. Note, however, that the suppression of the core recoil by Eq. (6) means that in LO the nucleon-core reduced mass that enters scattering is m N ; see Eq. (3).
Similar considerations apply to higher partial waves, but differences arise from the different renormalization: the more singular character of the interactions requires more LECs at any given order. For P-wave nucleon-core scattering, for example, both the scattering "length" a 1 and the effective "range" r 1 appear at LO [12,13] 1 . The mildest assumption [13] is that the effective range, just as for S-waves, is not fine tuned and directly reflects the breakdown scale, In this case r 1 k 2 is larger than the unitarity term ik 3 , and Smatrix poles of non-zero energy require a single fine tuning, Assuming the higher ERE parameters still scale with k hi , they are all subleading, and at LO there are two poles: depending on the sign of a 1 r 1 , a resonance on the real axis or a real/virtual bound-state pair with binding momentum in terms of the separation energy B s1 > 0. Thus, again, 1/k lo is related to the size of the halo system. At NLO the unitarity term needs to be included. If treated exactly, a third pole appears with momentum ∼ k hi , that is, outside the EFT. In the unlikely case where there are three low-energy poles in the low-energy region, the k hi in Eqs. (10) and (11) should be replaced by k lo [12]. Just as for S-waves, the assumption that no further powers of 1/A c appear in the ERE parameters implies that the only change in power counting when treating 1/A c as small is the extra expansion (3).

Charge form factor
The sizes of the halo system and its components are manifest not only in the ERE parameters but also in the charge form factor. The charge form factor is obtained as the matrix element of the zeroth component of the electromagnetic fourcurrent, J μ , according to where Z h is the proton number of the halo nucleus and Q is the momentum transfer. A measure of the size of the nucleus is the charge radius r ch . We now look at the power counting for the observable charge form factor of one-nucleon halo systems. The discussion here follows Ref. [22], but makes explicit the factors of 1/A c that were not incorporated into the power counting there. We discuss contributions to the charge form factor in the point-like limit, due to the finite size of the constituents, and from additional two-body (shortrange) contributions, as displayed in Fig. 1. More details and explicit examples will be presented in Sect. 3. We start by discussing the point-like part, r pt , which comes from the photon coupling to the charge of the constituents. The corresponding operators in the Lagrangian have the general form ψ † A 0 ψ where ψ denotes either the core or the nucleon field and A 0 is the zeroth component of the photon four-vector field. The point-like contribution is kinemat- Fig. 1 Charge form-factor diagrams for a nucleon-core system. A solid/dashed/double/wavy line denotes a nucleon/core/dicluster/photon field. An unmarked (solid square) photon vertex is due to minimal (non-minimal) coupling, which is independent of (quadratic in) the photon momentum. The diagrams (a-c) give contributions to the point-like part of the charge radius. Diagrams (d-f) enter with the finite-size contribution of the core (e) and the nucleon (f), and through a short-range contribution (d) ically generated by nucleon-core one-loop diagrams, where the photon couples to the core (in the neutron-halo case) or to both the core and the proton (for the proton halo)see Fig. 1b, c. For a neutron halo, momentum dependence requires recoil of the core, and thus the charge radius, which involves two powers of momenta, is proportional to f 2 , where f is defined in Eq. (4). Indeed, the loop for neutronhalo systems was calculated by Hammer and Phillips [22] with the leading-order (LO) results: The scalings are then given by assuming Eqs. (9), (10), and (12). However, it should be noted that typically r 1 < 0 so, if |r 1 | is close to 3γ 1 , then |r 1 + 3γ 1 | can in practice be ∼ k lo rather than the formally correct assignment ∼ k hi we have used here. The point-like contribution to the charge radius is interesting because it can be calculated from known properties of the constituents, but it exists against a backdrop of other, less well-known contributions. One type of these other contributions comes from the finite sizes of the constituent core and nucleon, which enter through the same loop diagrams as the point contributions; see Fig. 1e, f. The finite-size contributions correspond to operators of the form ψ † (∇ 2 A 0 )ψ. The two extra powers of the small momentum compared to the point-like terms coming from minimal substitution means that, by naive dimensional analysis, this operator carries a suppression of k −2 hi . Such terms produce a direct contribu-tion to the charge radius that is not suppressed by A −2 c . The nucleon finite-size contribution to r 2 ch will be denoted by ρ 2 N and should be proportional to ρ 2 N ∼ R 2 N , i.e., to the proton or neutron charge radius squared, respectively, ρ 2 p = 0.766 fm 2 [31] and ρ 2 n = −0.116 fm 2 [32]. This contribution scales as Meanwhile, the core-size contribution to the charge radius squared, ρ 2 c , will scale as ρ 2 c ∼ R 2 c ∼ 1/k 2 hi , S-wave or P-wave halo.
Taking the ratio of Eqs. (16) and (17) shows that in general the nucleon finite-size contribution to the charge radius squared is smaller than the core-finite-size contribution by both (R N /R c ) 2 and a factor of the total charge Z h of the system. For the canonical estimate There is another contribution to the charge radius, but this time it is not determined by data from other processes. Shortrange contributions to the charge density are encoded in a contact operator of the form Ψ † (∇ 2 A 0 )Ψ , where Ψ denotes the dicluster field for either the S-or the P-wave systemsee Fig. 1d. Because of the two derivatives, this operator is suppressed by a factor of k −2 hi with respect to Ψ † A 0 Ψ , which originates in the minimal substitution of the dicluster kinetic term- Fig. 1a. The minimal substitution term ensures that the charge comes out correct; the term with two additional derivatives comes with short-range parameters, which we denote in S-and P-wave halos by, respectively, ρ 2 σ,π ∼ k −2 hi . In the case of an S-wave halo, the dicluster kinetic term is itself a relative k −1 hi effect (it gives rise to the range, Eq. (8)), for an overall k −3 hi suppression. For a P-wave system, the dicluster kinetic term leads to the "range" which scales with k hi ; see Eq. (10), and there is no extra suppression [22]. We expect short-range contributions to the charge radius that scale as again shown explicitly in Sect. 3. Thus, for S-wave onenucleon halos this short-range operator enters one order after the core-size contribution. For P-wave halos both contribute at the same order.
In summary, these power-counting arguments make explicit that the point-like contribution for one-nucleon halos involves a kinematical suppression factor 1/A 2 c for neutron halos. But this has the consequence that, for P-wave one-neutron halos, a short-range operator enters at the same order as the finitesize core contribution. The existence of such additional shortrange operators will have negative influence on the predictive power of LO calculations.

The charge-radius formula
In this section we will derive charge-radius formulas in HEFT with the heavy-core power counting. In the process we will critique some results from Ref. [22] where finite-size effects were not treated explicitly. For example, the charge radius formula used in Ref. [22] for a one-neutron-halo system is where ρ 2 c is the charge radius squared of the core. In principle, one should also add the neutron charge-radius contribution ρ 2 n /Z c , where Z c is the core charge, but this term is usually neglected since |ρ 2 n | is tiny. The key point is that Eq. (20) has not been derived within the field theory: finite-size contributions were instead added to the point-like result in a rather ad hoc fashion. In what follows we will show that ρ 2 c (and, for that matter, ρ 2 n /Z c ) indeed add to the charge radius squared, but in principle other contributions of similar size can appear. We also argue that it is important to keep track of the large suppression in neutron halos of the point-like radius for A c 1, when the main contribution to the charge radius of the system comes from the finite size of the constituents.
This derivation, carried out in the next subsections, starts from the HEFT Lagrangian. We will consider explicitly only dominant S-or P-wave interactions, and a spin-0 coregeneralizations are straightforward but cumbersome to write. The Lagrangian for a spin-1/2 nucleon N s , where s = −1/2, 1/2, and a spin-0 core c with S-and P-wave interactions is given by The field σ s is the spin-1/2 dicluster field, which we introduce for convenience. Its kinetic term has a sign η 0 and it has a residual mass Δ 0 , which is a parameter to be fixed. The most important S-wave coupling of nucleon and core has strength g 0 . For the P-wave interaction, corresponding to the last two terms, we have, for simplicity, included only the J = 1/2 channel, through a spin-1/2 dicluster field π s and the Clebsch-Gordan coefficient C s si = 1 2 s1i 1 2 1 1 2 s [22]. The indices take values according to s, s = −1/2, 1/2 and i = −1, 0, 1. As above, η 1 is a sign, and Δ 1 and g 1 are parameters to be fixed, while f was defined in Eq. (4). The field A 0 is the zeroth component of the photon four-vector field. Here τ 3 is the third isospin Pauli matrix, and we have defined the charge number of the core Z c . Note that Z h is the charge (Z c or Z c + 1) of the nucleon-core system. The charge radius of the nucleon (core) field is ρ N (ρ c ), while ρ σ and ρ π are additional short-range parameters with sizes given by Eq. (19), which will be discussed below.
This Lagrangian includes all operators that will contribute to the charge radius of the halo system up to NLO. Higherorder terms, such as the one that leads to the shape parameter in the ERE and terms that do not contribute to the charge form factor, are denoted by the ellipsis. The coefficients of the ∇ 2 A 0 terms in Eq. (21) encode the finite size and composite nature of the core. We assign the same scaling to them as Hammer and Phillips, who used naive dimensional analysis to argue that ρ 2 c , ρ 2 σ , ρ 2 π are all ∼ R 2 c 2 . However, whereas Ref. [22] argued that this scaling rendered these effects of higher order, here we contend that they can provide the main contribution to the charge radius for heavy-core systems. + + ... Fig. 2 Diagrams for the charge form factor of the core. The first diagram has a vertex ieZ c and the second diagram gives the charge radius of the core, through the vertex −ieZ c ρ 2 c Q 2 /6 The finite constituent sizes are reflected in the non-trivial momentum-transfer dependence of the core and nucleon form factors. The form factor of the core is given by the part of the Lagrangian where a photon couples to the core field, that is, terms of the form c † A 0 c (and derivatives of A 0 ). The resulting charge form factor of the core is thus given by the diagrams shown in Fig. 2 as and as such the charge radius of the core is given by ρ c . Note that since Fig. 2 contains no loops the Lagrangian parameter ρ 2 c appears directly here without being affected by renormalization-at this order. The nucleon charge form factors and charge radii can be considered in a similar fashion, with the exception that the electric charge of the neutron is zero: where Z N = 0 (1) for the neutron (proton). We now examine the effect of these terms in the charge radii of S-and P-wave neutron halos, while the case of proton halos is discussed in the Appendix.

S-wave neutron halos
Here we compute the expectation value of the zeroth component of the electromagnetic current, J 0 , which appears in Eq. (13), for an S-wave one-neutron halo. The longdistance contributions to this quantity are well known, cf.
Refs. [22,25,30], where the diagrams in Fig. 1a, b were evaluated (although only for A c = 1 in Refs. [25,30]). Here we include diagrams Fig. 1d-f as well, and so account for finitesize effects and the leading short-range, two-body operator. The charge form factor of an S-wave one-neutron halo can be computed from the amputated correlator of the σ s field with one insertion of all possible couplings to an A 0 photon. The diagrams that contribute up to O((k lo /k hi ) 3 ) are shown in Fig. 1. We must also apply a wavefunction renormalization factor Z 0 , which defines the overlap of the field σ s with the physical one-neutron-halo state. The contributions from tree and loop diagrams can then be separated, viz.
The wavefunction renormalization factor is [22,30] Z 0 = 2πγ 0 where we kept some higher-order terms as well by not expanding the (1 − γ 0 r 0 ) −1 ratio. Note that Z 0 is finite to this order. At Q = 0 finite-size effects cannot play a role as the photon only "sees" the entire charge of the system. The leading contribution to the form factor at Q = 0 is then from the loop diagram in Fig. 1b, where the virtual photon is attached to the core via its charge. This diagram also gives rise to subleading corrections to F ch ( Q): it generates the point contribution to the form factor, but away from Q = 0, this is suppressed by 1/A 2 c and not as large as other effects once A c 1. The most important such other effect is due to the loop diagram Fig. 1e, i.e., the coupling of the A 0 photon to the finite size of the core inside the loop. The contribution of this graph can be combined with the corresponding coupling for the neutron, Fig. 1f, to obtain the contribution stemming from the constituent form factors, Eqs. (22) and (23). The result can be expressed as a coordinate-space integral, We expand the integral (26) up to order Q 2 to arrive at At O(k lo /k hi ), F ch ( Q) also receives a contribution from the tree-level diagram, Fig. 1a. Considering also the O((k lo / k hi ) 3 ) tree diagram in Fig. 1d, which represents the shortrange contribution to the halo charge radius, we arrive at The first term here is a consequence of charge conservation and ensures that as Q → 0 we have F ch (0) = 1, that is, the charge form factor is correctly normalized at zero momentum transfer. Moving away from Q → 0 we insert Eqs. (28) and (27) in Eq. (24), and compare with the term quadratic in momentum in Eq. (13), to obtain the charge-radius formula for S-wave neutron halos, 123 Table 1 Orders of the various contributions to the charge radius of S-wave neutron halos listed in Eq. (30). In each column effects of a particular order in the usual HEFT expansion parameter k lo /k hi appear.
Meanwhile the rows organize contributions due to additional small factors: inverse powers of the number of core nucleons (A c ) and protons (Z c ) where the "· · · " represent higher-order contributions. Note that at this order the Lagrangian parameter ρ 2 σ appears directly here, without renormalization. This is because the LO loop that gives the point-charge contribution is finite. The point-charge contribution to r 2 ch was computed in Ref. [22] as (1 − γ 0 r 0 ) −1 r 2 pt,LO , with r 2 pt,LO given by Eq. (14). Expanding in γ 0 r 0 , where the orders of various contributions are summarized in Table 1, assuming that R N /R c A −1/3 c . If A c ∼ 1, the most important terms are given by the point-radius r 2 pt,LO . For larger A c , these terms rapidly lose importance, as do contributions from the neutron size. Unless we are dealing with light halos, we expect the dominant contribution to the difference in charge radii between halo and core to be given by γ 0 r 0 (ρ 2 c −ρ 2 σ ). This is an example of a term that is missed if one simply adds the core radius by hand, rather than including it in the EFT as any other operator.
Unfortunately, while ρ c can be extracted from the core form factor (22), ρ σ is a short-range term that cannot easily be extracted from a quantity other than the halo form factor itself. Since this is a short-range effect it is possible that it can be efficiently extracted from ab initio calculations of the charge radius. In such a calculation the difference between ρ σ and ρ c can be viewed as originating in two effects: 1. A change in the size of the core when it is placed in the bound state with the neutron. 2. Pieces of the ab initio wave function not in the core + neutron piece of the Hilbert space, e.g., those due to excited states of the core.
These effects cannot, however, be separated in a modelindependent way, and only their combination enters through ρ 2 c − ρ 2 σ .
As a concrete example we consider the S-wave ground state of the one-neutron halo 11 Be, whose form factor and photodisintegration were investigated in Ref. [22]. The neutron separation energy is B s0 = 0.502 MeV [33], corresponding to k lo ∼ γ 0 30 MeV through Eq. (9). Using the charge radius of the 10 Be core, ρ c = 2.357(18) fm [34], as an estimate for its size, the breakdown scale is k hi ∼ 1/R c 80 MeV. This is also the momentum √ 2m R E ex 80 MeV ∼ k hi corresponding to the first excitation of the core at E ex = 3.368 MeV [35], so there is no need to include a field for this state. These scales then give us the expansion parameter k lo /k hi ≈ 0.4. This means that r 2 pt,LO is numerically of the same size as O((k lo /k hi ) 3 ) corrections. Since Z c ∼ A c k lo /k hi , the neutron-radius contributions are suppressed by more than five powers of k lo /k hi compared to ρ c . At LO there is a charge-radius prediction, but it is trivial since it is just the charge radius of the 10 Be core, ρ 2 c 5.56(4) fm 2 . This does, though, explain most of the measured value of r 2 11 Be 6.07(8) fm 2 (r11 Be = 2.463 (16) fm [34]). Estimating γ 0 r 0 from the EFT expansion parameter ∼ 0.4, we find that the point-charge contribution to r 2 ch , i.e., the first term of Eq. (29), is r 2 pt,LO / (1 − γ 0 r 0 ) 0.3 fm 2 . This explains more than half of the difference r 2 11 Be − r 2 10 Be . The rest must come from the short-distance effect ρ 2 σ : the experimental value for r11 Be can be reproduced if the shortrange parameter is given by ρ 2 σ ≈ 5.1 fm 2 , which is of the expected order of magnitude, 1/k hi ∼ 2.5 fm. This supports the power counting presented here. We thus see that ρ 2 σ must be a little smaller than ρ 2 c in order to explain the data, although the errors on the atomic measurements of the 10 Be and 11 Be radii make it difficult to extract a precise value for ρ 2 c − ρ 2 σ .

P-wave neutron halos
An important aspect of the S-wave halo system is that all the charge form-factor diagrams we considered are finite. For P waves this is not the case. The increased singularity of the P-wave interaction can be seen already in the need for the effective-range term (10) at LO to allow proper renormalization of nucleon-core scattering. As before we will consider operators up to second order in the photon momentum Q and will show that, if the charge-radius contributions of the constituents are to be considered explicitly, we need an additional short-range operator to renormalize the halo charge radius.
Since the cancelation of divergences will be critical to what follows we recapitulate the formulas for neutron-core scattering derived in Refs. [12,22]. The power counting discussed in Sect. 2.1 indicates that neutron-core scattering proceeds through the bare dicluster propagator at LO, and by an insertion of one nucleon-core bubble at NLO. Up to this order, elastic scattering with a P-wave interaction gives the matching where the L n = d p p n−1 are divergent integrals in the irreducible dicluster self-energy, It is evident from Eqs. (31) and (32) that two parameters, Δ 1 and g 1 , are needed to renormalize scattering up to NLO. The P-wave wavefunction renormalization is given by Note that, contrary to the S-wave wavefunction renormalization (25), Z 1 is not finite to this order. The P-wave charge form-factor diagrams are similar to those for the S-wave interaction, Fig. 1. The tree diagrams amount to while the one-loop diagrams give Expanding in powers of the momentum transfer Q, The only difference with respect to the S-wave, apart from Z 1 , is that the P-wave bound-state wavefunction is 1 + 1/(γ 1 r ) exp (−γ 1 r ), which is irregular at the origin. As a consequence, the integral that appears already in the momentum-independent term is divergent, and related to one of the divergent integrals in Σ(E), the L 1 of Eq. (32). The divergence in the momentum-independent contribution cancels between Eqs. (35) and (37), and we obtain a properly normalized form factor, F ch (0) = 1 [22]. The terms quadratic in momentum give the charge radius, where the LO point-charge contribution, defined in Eq. (14), agrees with Ref. [22], and the (finite) short-range contribution ρ 2 π is related to the counterterm ρ 2 π bȳ An interesting point here is that the finite contribution of ρ 2 c to r 2 ch from Eq. (37) is suppressed by an additional factor γ 1 /r 1 ∼ k lo /k hi with respect to the estimate (17). The appearance of the full ρ 2 c in (38) is a consequence of the particular renormalization condition (39).
In this renormalization scheme the effect beyond the "standard" charge-radius formula depends on the extent to which the dicluster counterterm differs from the core radius. In contrast to the S-wave case, here the difference ρ 2 π − ρ 2 c − ρ 2 n /Z c must go to zero as the regulator is taken to infinity, in order to yield a finiteρ 2 π . It is important to consider what would happen if we were to include the finite-size contributions, but not the ρ 2 π short-range operator. In Eq. (37) we see that the constituent charge radii enter with a prefactor that corresponds to a divergent integral. Since the parameters ρ 2 c and ρ 2 n are observables-these are the charge radii of the core and the neutron-they cannot absorb this divergence. The only parameter available for this purpose is the ρ 2 π . As such it is not possible to add the finite-size contributions without also including the short-range operator. Formally, the scaling of ρ 2 π is ρ 2 π ∼ρ 2 π ∼ R 2 c for renormalization scales such that L 1 ∼ k hi . However, the crucial difference between ρ 2 π and ρ 2 π is that the latter is an observable, while the former absorbs a divergence and so is scheme-and regulator-dependent.
As an explicit example, let us consider the P-wave excited state of 11 Be with neutron separation energy B s1 = 0.182 MeV [33]. The breakdown scale for this EFT was argued in Sect. 3.1 to be k hi ∼ 80 MeV. These scales then give us the expansion parameter k lo /k hi ≈ 0.2 for the P-wave system. The corresponding charge-radius formula is simply organized as r 2 11 Be * = ρ 2 c +ρ 2 π 1/k 2 hi + · · · + r 2 pt,LO + · · · 1/(k lo k hi ) × 1/A 2 c + · · · In this case, the LO result is given by the combination of the charge radius of 10 Be and an undetermined short-range parameter. The dots refer to corrections due to non-included interactions and the finite neutron size. We show the pointcharge contribution explicitly to emphasize that it appears at N 2 LO in the heavy-core power counting. This means that the charge radius for the P-wave state in 11 Be cannot be predicted in HEFT using the heavy-core power counting, unless the short-range parameter can be fixed to some other observable.

Conclusion
HEFT offers a systematic approach to make modelindependent predictions of low-energy observables. In this paper we have discussed a new power-counting scheme for systems with a heavy-core nucleus, and we have derived the finite-size contributions to charge radii of one-nucleon halos. HEFT in general is restricted by appearances of short-range operators at rather low orders. With the heavy-core power counting, these restrictions are even enhanced for some systems and observables. For one-neutron halos where the core is much heavier than the neutron, the point-particle result for the charge radius is demoted from leading to subleading order since the core recoil due to the photon interaction is very small. In contrast, in the case of an S-wave system, the LO charge radius is given by the finite-size contributions of the constituents. For a P-wave one-neutron halo the heavycore version of HEFT is non-predictive at LO, since the LO charge radius includes an undetermined short-range operator 3 . Note, however, that not all systems are made less predictive in the heavy-core power counting. For proton halos there are no issues for the charge-radius results due to the core being heavy (as shown in the appendix). This is due to the fact that the photon also couples to the proton field, which has a larger recoil than the core field. Furthermore, the expectation for the future is that more cluster data will become available and that this data can then be used to fix the parameters of the corresponding HEFT.
While we considered in detail the case of one-nucleon halo charge radii, the suppression of some contributions by factors of the inverse of the number of core nucleons is not restricted to this class of observables. The suppression for radii can be traced to the small recoil of the core or, equivalently, to the fact that the heavy-core propagator is static at leading order. Similar effects will in principle be present in any calculation at the loop level, where the propagator appears, for example the structure (energies, form factors, etc.) of two-nucleon Data Availability Statement This manuscript has no associated data or the data will not be deposited. [Authors' comment: The manuscript has no associated data, as all calculations are fully explained and justified in the text.] 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/.

A Single-proton halos
For proton halos one needs additionally to account for Coulomb effects. Denoting by Z c the charge of the core and by α = e 2 /(4π) the fine-structure constant, the strength of the Coulomb interaction is characterized by the momentum scale k C ≡ Z c αm R . At energy E the relative importance of Coulomb is given by the Sommerfeld parameter For moderate Z c , as in light nuclei, we expect k C k lo . As Z c increases the Coulomb force experienced by the halo proton increases. In that case, it might be appropriate to consider the limit k lo /k C 1. The effects of Coulomb in proton halo systems have been examined in Refs. [28,[37][38][39][40][41][42].
The charge form factor for proton halos involves the coupling of the photon to both the core and the proton. The pointlike contribution to the charge radius is therefore not kinematically suppressed by 1/A 2 c and the scalings are naively given by 1/k 2 lo , S-wave proton halo, 1/(k lo k hi ) , P-wave proton halo. (A.1) These naive scalings are only valid if the Coulomb momentum is a low-momentum scale, that is, k C k lo . If instead we are in the strong Coulomb regime, k C k lo , the predictive power of LO calculations is reduced, which has been analyzed for S-wave one-proton halo states in Ref. [41,42]. For example, the point-like contribution for S-waves scales as r pt ∼ 1/k C if k C k lo , which implies that the pointlike contribution becomes suppressed by the strong Coulomb repulsion.

A.1 S-wave proton halos
The procedure for deriving the charge-radius formula for Swave proton halos is similar to the neutron case; see Refs. [38,41]. However, there are four main differences: (i) The total charge of the system is Z c + 1.
(ii) The Coulomb interaction enters proton-core scattering and the wavefunction renormalization is given by , NLO, with ψ being the polygamma function. For k C γ 0 , As before, some higher-order terms are kept at NLO in Eq. (A.2). (iii) The photon couples also to the proton in the proton-core loop of Fig. 1c, according to Eq. (23). (iv) Coulomb interactions enter the loops in Fig. 1. The bound-state wavefunction is the Whittaker W -function W −k C /γ 0 ,1/2 (2γ 0 r ) instead of the exponential exp (−γ 0 r ).
Taking these differences into account, the resulting chargeradius formula for an S-wave proton halo system is The leading-order point-charge contribution without the effective-range correction is given by where the loop diagram is given in Ref. [41] as where j 0 is a spherical Bessel function of the first kind. The 1/2 + excited state of 17 F was considered in HEFT by Ryberg et al. [38,41]. The 1/2 + excited state is located at 0.105 MeV below threshold [43], which then defines the low-momentum scale k lo ∼ 14 MeV. The first excitation of the 16 O core is at about 6 MeV [43], so the size of the core defines the breakdown scale k hi of about 1/R c ∼ 60-70 MeV, giving an expansion parameter k lo /k hi ∼ 0.2. However, for the 16 O-proton system the Coulomb momentum scale k C = 51.2 MeV is much larger than the low-momentum scale and 3k C r 0 is very close to unity [41]. This makes the effective-range prefactor in Eq. (A.5) very large, so this proton halo state cannot be well described without the inclusion of effective-range corrections. In practice, the ρ 2 σ counterterm then enters at the same order as the finite-size contributions. Furthermore, in this strong Coulomb regime the LO pointlike charge-radius contribution, r pt,LO , scales with 1/k C , as was discussed in Ref. [41]. Therefore, organizing the chargeradius formula for the S-wave proton halo at hand, we have since 3k C r 0 ∼ 1 and k C γ 0 for 17 F * . The point-like contribution to the charge radius of 17 F * , including the finite-range correction, was evaluated by Ryberg et al. [41] to r 2 pt,LO / (1−3k C r 0 ) = (2.20 ± 0.11 fm) 2 . Here, the assigned uncertainty originates from the asymptotic normalization coefficient (ANC) obtained from a fit to proton radiative-capture on 16 O. One should note that corrections from finite-size contributions and the short-range operator can be as large as k C /k hi ≈ 70%, unless at least parts of those corrections can be resummed. But Eq. (A.8) means that-at least formally-one may not add the finitesize contributions to the charge-radius result of 17 F * without also including the unknown short-range parameter ρ 2 σ .
A.2 P-wave proton halos As for the S-wave proton halos, many of the details in the derivation of charge-radius formula for a P-wave proton halo are the same as for the P-wave neutron halo. Again, the differences are that the photon couples also to the proton, the bound-state wavefunctions are Whittaker W-functions, the total charge is Z c + 1 and the wavefunction renormalization is now given by with h 1 (η) = p 2 (1+η 2 )h 0 (η). For details, see Refs. [28,39]. The resulting charge-radius formula for a P-wave proton halo system is whereρ 2 π is the renormalized short-range parameter. Here we emphasize the contribution r 2 pt,LO is not kinematically suppressed for systems where both constituents are charged, since the photon couples through minimal substitution also to the proton. Therefore, ifρ 2 π is subleading to r 2 pt,LO , it is possible to give a low-order charge-radius prediction of P-wave proton halos without including the short-range parameterρ 2 π . The point-charge contribution r 2 pt,LO was derived in Ref. [39] and it is given by The one-proton separation energy of 8 B is B s1 0.138 MeV [35], giving a low-momentum scale k lo ∼ 15 MeV. The 7 Be core has two low-lying states that were both included explicitly into the field theory: the 3/2 − ground state and the 1/2 − excited state at 0.429 MeV [44]. It can be argued that the breakdown scale for the EFT is given by the alphaparticle threshold at a momentum scale of k α 51 MeV [44] and thus the expansion parameter is about k lo /k hi ∼ 0.3, or even as large as k C /k hi ∼ 0.5 with k C = 23.8 MeV.
The leading-order contribution to the charge radius of 8 B was evaluated by Ryberg et al. [39] to r 2 pt,LO = (2.56 ± 0.08 fm) 2 . Again, the uncertainty estimate comes from the relevant ANCs, which in this case were adopted from a microscopic ab initio computation by Nollett and Wiringa [45]. Alternatively, one can obtain the ANCs from a fit to proton radiative-capture 7 Be( p, γ ) 8 B. Such a fit was performed by Zhang et al. [28] and the ANCs are very consistent with the computed ones. However, the large expansion parameter suggests the finite-size and (unknown) short-range contributions to the charge radius of 8 B can be significant.