Low-energy dipole strength in 8He

In this work, we present new ab initio coupled-cluster calculations of dipole-excited state properties of 8He based on the chiral effective field theory interaction 1.8/2.0 (EM). We focus on the dipole polarizability, and compare the results to our previous study [Phys. Rev. C 105, 034313 (2022)] and subsequent theoretical work. With the aim of connecting the presence of low-lying dipole strength to structure properties of 8He, we compute the point-neutron radius, finding excellent agreement with available experimental data, and investigate its correlation with the dipole polarizability.


Introduction
At the boundaries of the nuclear chart, the force binding protons and neutrons in the nucleus gives rise to intriguing phenomena, such as the emergence of halo structures, with weakly-bound nucleons orbiting a compact core at a distance.Halo nuclei are characterized by small separation energies and extended matter radii, which do not follow the typical A 1/3 scaling found in the valley of stability.
Among halo systems, 8 He represents an interesting case study.It has the most extreme neutron-to-proton ratio (N/Z = 3) in the nuclear chart and up to now, it is the only known four-neutron halo.This unique configuration has been leveraged in a recent groundbreaking experiment, performed at RIKEN, where a 8 He knockout reaction at large momentum transfer led to the observation of a correlated four-neutron state [1].While from the theory point of view the existence of a bound tetraneutron is disputed (see e.g.Ref. [2] and references therein), new experiments have either been made [3] or they are being planned for [4].The debate around the interpretation of the RIKEN 4n signal has triggered new investigations of the neutron halo distribution and the possible correlations among the loosely-bound neutrons in 8 He [5,6].
In halo nuclei, the extended size of the neutron cloud is often mentioned together with a strong enhancement of the electric dipole response at low excitation energies.Revealed by early measurements with radioactive ion beams [7], this so-called "soft E1 excitation" ranks among the main features identifying halo nuclei [8].However, the presence of a soft dipole mode in 8 He is still a controversial issue from both the experimental and theoretical point of view.
Coulomb excitation [9] and nuclear fragmentation [10] experiments first supported the existence of a 1 − resonance at an excitation energy of around 4 MeV in the spectrum of 8 He.Later, (t, p) transfer reaction data [11] confirmed the presence of low-lying dipole strength, but at a lower energy of 3 MeV.In contrast, the Coulomb excitation experiment of Ref. [12] attributed a relatively small fraction of the total energy-weighted dipole sum rule (less than 3%) to a potential soft dipole mode.Also, a measurement of breakup of 8 He on carbon [13] found the spin-parity assignment of the excited state at 4 MeV to 1 − highly uncertain.These results were validated by a recent inelastic proton scattering experiment [14], where the measured angular distribution was found to be incompatible with a low-lying dipole resonance.Unpublished, highstatistics Coulomb excitation data obtained at RIKEN by the SAMURAI collaboration are expected to shed light on this controversy [15].
From the theory point of view, our recent work [16] was the first to tackle dipole excited-state properties of 8 He in the framework of ab initio calculations based on Chiral Effective Field Theory (χEFT) interactions [17][18][19].Employing Coupled-Cluster (CC) theory [20] merged with the Lorentz Integral Transform (LIT) technique [21] in the so-called LIT-CC approach [22,23], we observed low-energy dipole strength at around 5 MeV.Subsequent works based on the cluster orbital shell model (COSM) [24,25] and the equation-of-motion multiphonon approach (EMPM) [26] support the presence of a soft dipole mode.In contrast, a random-phase approximation (RPA) calculation within a density functional theory (DFT) framework [27] disfavours such scenario, attributing it to spurious center-of-mass contaminations.
Motivated by this ongoing debate in both the theoretical and experimental communities, in this work we complement our previous analysis of 8 He in two ways.First, we provide new LIT-CC predictions of dipole excited-state properties of 8 He employing the χEFT interaction 1.8/2.0(EM) [28].This interaction contains SRG-evolved two-nucleon forces at next to next to next to leading order (N3LO) and three-body forces at next to next to leading order (NNLO).The nuclear force models used in Ref. [16] (namely, NNLO sat [29], ∆NNLO GO (450) and ∆NNLO GO (394) [30]) are instead derived at NNLO in the chiral expansion.Second, using these four χEFT interactions, we investigate the halo structure of this nucleus by calculating the point-neutron radius and explore its correlation with the dipole polarizability.
The structure of this paper is the following.In Section 2, we give an overview of the LIT-CC method employed in this work.In Section 3, we present our new results, devoting one subsection to the discretized response function and dipole polarizability and one subsection to the analysis of correlations between dipole polarizability and point-neutron radius.In Section 4, we draw our conclusions.

Review of the method
Electric dipole excitations of nuclei can be studied introducing the dipole response function, which is defined as: where |Ψ µ ⟩ are the excited states obtained from the ground state |Ψ 0 ⟩ applying the dipole operator Θ, and ω is the photon energy.
To rewrite Eq. ( 1) in the CC formalism [20], we will briefly go through how |Ψ 0 ⟩ and |Ψ µ ⟩ can be calculated in this approach.Given the nuclear Hamiltonian H, in CC theory one starts from a reference state |Φ 0 ⟩ expanded on the harmonic oscillator basis, and constructs the correlated many-body wavefunction via an exponential ansatz The cluster operator T can be written in terms of a sum of n-particle n-hole excitations, 2), the Schrödinger equation becomes where H = e −T He T is the similarity-transformed Hamiltonian.As the latter is non-Hermitian, left and right eigenstates of H are different.While the right eigenstate of H corresponds to the reference state |Φ 0 ⟩, the left eigenstate can be determined via the relation |Ψ 0,L ⟩ = (1 + Λ) |Φ 0 ⟩, where Λ is a de-excitation operator having a structure analogous to T .Excited states are computed employing the equation-of-motion coupled-cluster (EOM-CC) method [31].The latter is based on the following ansatz for right and left excited states: where the EOM operators R µ and L µ can also be expanded in terms of particle-hole excitations.The EOM-CC ansatz allows to reformulate the computation of excited states as a right (left) matrix eigenvalue problem, with the excitation energies E exc,µ = E µ − E 0 as eigenvalues and the amplitudes of R µ (L µ ) as eigenvectors.
In CC theory, the particle-hole expansion of T and Λ for the ground state and of the EOM operators R µ and L µ for the excited states is truncated due to computational limitations.The most frequently employed approximation is coupled-cluster singles and doubles (CCSD), where we include up to 2p-2h excitations.Higher accuracy can be obtained within the CCSDT-1 framework, where leading order 3p-3h configurations are also present.In this work, we choose to adopt the same approximation scheme for the ground state and the excited state computations, either CCSD or CCSDT-1.Substituting Eq. ( 2), (4) in Eq. ( 1), we can express R(ω) as where Θ = e −T Θe T is the similarity-transformed dipole operator.We point out that the sum over µ in Eq. ( 1) encompasses both bound and continuum excited states of the nucleus.This makes the calculation of response functions particularly challenging, as all the possible unbound configurations arising from the break-up of the nucleus at a given energy are involved.To avoid the issue of explicitly computing nuclear continuum wave functions, one can resort to the Lorentz Integral Transform (LIT) method.The LIT technique is based on the calculation of an integral transform with Lorentzian kernel of the response: Within the so-called LIT-CC method [22,23], the LIT is connected to the solution of a EOM-CC equation with a source term.Once L(σ, Γ) is determined, one can recover R(ω) via a numerical inversion procedure [21].
Starting from the LIT, one can access the moments of the response function, i.e., the moments of R(ω) interpreted as a distribution function, which can be written as: where n is an integer.As the Lorentzian kernel tends to a Dirac delta function in the limit Γ → 0, we get: Taking advantage of Eq. ( 8), we can directly compute the moments from the LIT as: It is worth pointing out that this strategy does not require the inversion of L(σ, Γ), which represents an additional source of uncertainty, and it has been proved to be equivalent to the integration of the response [32].
In this work we will focus in particular on the inverse-energy-weighted dipole sum rule m −1 , which is proportional to the electric dipole polarizability α D : 3 Results

Discretized response function and dipole polarizability
In this Section, we revisit our previous analysis of the dipole excited-state spectrum of 8 He, by considering new calculations performed with the chiral force 1.8/2.0(EM) [28].We first consider the discretized dipole response function, that in our approach corresponds to the limit of the LIT for Γ → 0. Fig. 1 shows the LIT for Γ = 10 −4 MeV for the 1.8/2.0(EM) interaction in comparison to the results obtained with the χEFT potentials used in our previous work [16].We observe that for the 1.8/2.0(EM) interaction the response is characterized by low-energy strength emerging at slightly lower energies than predicted by the other χEFT interactions employed.This reinforces the consistency of our coupled-cluster predictions with the experiments of Ref. [9][10][11], as well as with the COSM [24,25] and EMPM [26] theoretical calculations.
Starting from the discretized response function and employing Eq. ( 10), we can study the electric dipole polarizability α D .In Table 1, the theoretical value of α D for the 1.8/2.0(EM) interaction is reported together with the results given by the set of chiral EFT Hamiltonians used in our previous work 1 .RPA [27] and EMPM [26] predictions of α D are also included for comparison.
Table 1: Theoretical predictions for the dipole polarizability of 8 He with the CCSDT-1 approximation for the four chiral EFT interactions under analysis.RPA predictions based on three relativistic mean-field (RMF) energy density functionals [27] and the EMPM result [26], obtained with the NNLO sat interaction, are also reported.The uncertainty bands indicated in Table 1 for the coupled-cluster results take into account (i) the convergence with respect to the model space N max and (ii) the effect of the many-body coupled cluster truncation.(i) is estimated considering the residual ℏΩ dependence at the highest value of N max available, corresponding to N max = 14.(ii) is instead calculated comparing the two available CC approximation schemes, namely CCSD and CCSDT-1, and taking half of their difference as truncation uncertainty [33].(i) and (ii) are then summed in quadrature, leading to an uncertainty ranging between 6 and 8% of the central value for the different interactions.In Ref. [34] we discussed the uncertainty stemming from the truncation of the χEFT expansion on α D for this nucleus.Focusing on the ∆NNLO GO (450) model, NLO and NNLO results came out to agree with errorbars.
Looking at Table 1, we notice that while α D values from RPA [27] and EMPM [26] range between 0.20 and 0.26 fm 3 , CC predictions are appreciably larger.Moreover, it is worth to point out the significant model dependence shown by our CC results for the polarizability of 8 He.The NNLO sat interaction, which performs well in reproducing the dipole polarizability in the medium-mass region [35], yields the lowest value for α D in 8 He.The 1.8/2.0(EM) interaction delivers instead the largest prediction, which is more than 20% higher than the NNLO sat result.This is a consequence of the location of the low-lying dipole strength in 8 He, which appears at slightly lower energies in the 1.8/2.0(EM) case with respect to the other interactions, as shown in Fig. 1.Following 1 The value quoted for ∆NNLOGO(394) in our previous work was 0.39(2) fm 3 , and it was obtained varying ℏΩ between 12 and 16 MeV with the largest model space size available (Nmax = 14).We update here our prediction on the basis of new calculations including ℏΩ = 10 MeV.We point out that the two values are compatible between errorbars.
the line of previous works as Ref. [36], a global sensitivity analysis of α D in 8 He would help in identifying which low-energy constants in the interaction are driving the substantial spread of the results obtained with different chiral forces.

Neutron radius and correlation with dipole polarizability
To assess our understanding of the halo structure of 8 He and its impact on low-energy dipole excitations, we calculate the point-neutron radius R n with the previously employed four χEFT interactions and investigate correlations with the dipole polarizability.
Fig. 2 shows coupled-cluster predictions for α D as a function of R n , in comparison to three experimental determinations of R n [37][38][39] and the theoretical estimates obtained by the RPA-DFT [27] and EMPM [26]  Fig. 2: Electric dipole polarizability of 8 He as a function of the corresponding point-neutron radius for the χEFT interactions 1.8/2.0(EM), ∆NNLO GO (450), ∆NNLO GO (394) and NNLO sat .CC predictions for α D and R n have been obtained in the CCSDT-1 approximation scheme.Comparison with the experimental determinations of Ref. [37][38][39] as well as with the RPA-DFT [27] and EMPM [26] results is shown.tions of R n have been obtained on the basis of interaction cross section measurements in inverse kinematics on carbon [37] and hydrogen targets [38,39], respectively.A Glauber analysis of the experimental data allows for the extraction of the nuclear matter radius R m , which in combination with the point-proton radius R p can be employed to estimate R n according to the relation R We find that our theoretical predictions for R n are in excellent agreement with all three experimental determinations.We estimated our theoretical uncertainty following the procedure outlined in Sect.3.1.The residual dependence on convergence parameters and the truncation of the many-body CC expansion lead to a combined errorband varying between 1 and 2% for the different interactions.The EMPM value for R n is also in good accordance with experiment, while DFT results overestimate the data.
Looking at our predictions for different chiral forces, we notice the presence of a correlation between α D and R n , which emerges also from the RPA-DFT calculations of Ref. [27].This may reflect a link between the low-lying dipole strength, which mainly determines α D , and an excitation of the weakly-bound excess neutrons.We can arrive to similar conclusions also for neutron-rich nuclei in the medium-mass region, as 48 Ca [33] and 68 Ni [40].In fact, Refs.[33,40] highlight the presence of a correlation between α D and the point-proton radius R p for these two nuclei.At the same time, the authors also observed that increasing values of R p were associated to a corresponding rise in R n , leading as a consequence to the emergence of a correlation between α D and R n .

Conclusions
Motivated by recent theoretical and experimental interest towards the possible presence of a soft dipole mode in 8 He, we reconsider the analysis of our previous work [16] in the light of new coupled-cluster calculations of the discretized response function and dipole polarizability with the chiral force 1.8/2.0(EM).As a difference with respect to the Hamiltonian models employed in Ref. [16], derived at NNLO in the chiral expansion, 1.8/2.0(EM) contains also the SRG-evolved two-nucleon forces at N3LO.
The 1.8/2.0(EM) interaction confirms the presence of low-lying dipole strength in 8 He at around 5 MeV, and predicts the largest value for the polarizability, corresponding to 0.48(3) fm 3 .The presence of low-energy dipole strength in 8 He still remains an open problem which calls for further theoretical and experimental investigation.The new Coulomb excitation data from RIKEN [15] are expected to sort out the soft dipole mode controversy and serve as a benchmark for theoretical predictions of α D .
Moreover, employing the 1.8/2.0(EM) interaction and the set of chiral forces from our previous study, we computed the point-neutron radius, obtaining excellent agreement with experiment.We pointed out the presence of a correlation between dipole polarizability and point-neutron radius, which indicates the sensitivity of the low-energy dipole strength to the extent of the pont-neutron distribution in such an exotic halo nucleus. of the Department of Energy under Contract No. DE-AC05-00OR22725, and of the supercomputer Mogon at Johannes Gutenberg Universität Mainz.