Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

# The 29F nucleus as a lighthouse on the coast of the island of inversion

## Abstract

The exotic, neutron-rich and weakly-bound isotope 29F stands out as a waymarker on the southern shore of the island of inversion, a portion of the nuclear chart where the effects of nuclear forces lead to a reshuffling of the single particle levels and to a reorganization of the nuclear structure far from stability. This nucleus has become very popular, as new measurements allow to refine theoretical models. We review the latest developments and suggest how to further assess the structure by proposing predictions on electromagnetic transitions that new experiments of Relativistic Coulomb Excitation should soon become able to measure.

## Introduction

Weakly bound nuclear systems severely test our ability to disentangle the mysteries and oddities of the nuclear interactions and how nuclear systems achieve stability. One of the latest conundrums in this respect is the structure of the exotic 29-fluorine isotope. With nine protons and 20 neutrons, it is located almost on the edge of the stability valley of the nuclide chart, very close to the neutron drip-line, i.e., the dividing line (Sn = 0, null separation energy) between bound and unbound neutron-rich nuclei. Pioneering researches on this isotope have recently skyrocketed, as more and more theoretical and experimental papers are being published. The main interest lies in understanding if it belongs to the so-called island of inversion, a portion of the nuclear chart where the standard list of single-particle energy levels (that are filled by nucleons, very similarly to the atomic physics counterpart, in a sort of Aufbauprinzip, a.k.a. the building-up principle, i.e., the rule that states how electronic orbitals are filled up in the atomic shell model) shows an inversion between orbitals belonging to the sd- and pf-shells, see Fig. 1. This is crucial to ascertain the presence of a neutron halo (a diffused tail of nuclear matter that spreads around the central core) and its extent.

This problem relates also to the disappearance of the N = 20 neutron magic number when moving away from doubly magic nuclei. This shell evolution has profound connections with tensor nuclear interactions1,2 and with the deformation of the nuclear surface, occurring typically at mid-shell. This system, due to its mass A = 29, is still a hard nut to crack for ab initio approaches exploiting gargantuan numerical calculations. At present, simpler schematic models based on a few degrees of freedom are the preferred way to discern any pattern in the limited set of spectroscopic data.

## Theoretical approaches and new experiments

In one of the earliest theoretical three-body models for this nucleus3, the interactions between the 27-fluorine core and each of the two neutrons were tuned to reproduce the scarcely known properties of the binary subsystem 28-fluorine. Four different possible scenarios have been envisaged, called A, B, C, D, that we might rename as standard, intruder, degenerate, and inverted scenarios, referring to the order of single-particle levels and resonances appearing in the phase-shifts. The consequences of these choices have been analyzed and the relationships between the matter radius and the two-neutron separation energy with the percentage of p-wave content of the ground-state wave function have been explored3. These results show a mild to moderate sensitivity of these properties on the composition of the ground-state configuration and hint at the presence of a moderate halo structure connected with the p-orbital. The main reason to setup such a model is that 28F is neutron unbound with ground state at 220(50) keV and an unassigned resonance at 810 keV according to ref. 4,5, where the g.s. value has now been remeasured in ref. 6 to 199(6) keV, and the correlations between the two neutrons bring the nuclear system 29F back to stability, in a way similar to 6He and 11Li. When the binary subsystems are all unbound, as in all these cases mentioned, we speak of a Borromean nucleus, the archetypal system for studying two-neutron correlations. Another recent theoretical work7, based on the Gamow shell model that treats bound and continuum states on the same footing, was applied to neutron-rich fluorine isotopes up to 31F. The comparative plots of one-body density distributions suggest that 31F has a halo structure with neutrons occupying an orbital with mixed configurations: 2p3/2, 1d3/2, and 1f7/2, while 27F and 29F have instead a smaller radius, less compatible with a halo. These calculations produce a 29F-wave function where the p3/2 components amount to only a few percent8, more akin to the standard scenario A above. A novel mechanism to form the halo structure in 31F was discussed in ref. 9, in which the effect of pairing on the energy gaps between these orbits are essential to determine the configuration mixing. This puzzle is not to be resolved by theory alone, without further input from experiments.

Now, a new proton/neutron removal experiment on exotic beams6, performed at the Radioactive Isotope Beam Factory (RIBF) of the RIKEN Nishina Center in Japan, indicates that 28F exists as a resonant state at about 199 keV above the 27F + n threshold with a dominant  = 1 orbital angular momentum component (79%). This measurement is consistent with the position of the p-wave resonance in models C and D above, i.e., 200 keV (cfr. Table 1 in ref. 3, using some prescriptions from ref. 10). The analysis of this measurement is in stark contrast with the conclusions of ref. 4 where it is stated that “…pf shell intruder components play only a small role in the ground-state structure of 28F…”. The new data are interesting because they also show the presence of a resonance at about 0.9 MeV with a strong  = 2 component (72%). This is located half-way between scenarios C and D.

Previous measurements11 found a resonance at 1080(18) keV in 29F, to which the Jπ = 1/2+ quantum numbers were assigned. The analysis of shell-model calculations with effective interactions, contained in the same paper, suggested that the ground state is probably a 5/2+ state, arising from the unpaired proton in the 1d5/2 orbital, while the 1/2+ state arises mainly from the coupling of this proton with the 2+ excitation of 28O. The latter is still unmeasured and requires further confirmation. The fact that the two-neutron separation energy was measured by Gaudefroy et al.12 at 1443 (436) keV, implies that this excited state is also bound, a singular coincidence for such a neutron-rich system that awaits further confirmation through the observation of gamma-ray transitions to the 29F ground state. The neutron–neutron correlations at N = 20 are probably so strong as to make room not just for one, but for two bound states! These data are collected in Fig. 2 for the purposes of the present perspective, but we refer the reader to the original papers for more details.

When the older data11,12 on 29F were published, a theoretical study13 interpreted them as arising from the Nilsson deformed harmonic oscillator sd-states (N = 2 shell) within the particle-rotor model with a moderate deformation. Now, the newly measured dominant  = 1 component is incompatible with this view, apparently the model space should be enlarged to include intruder configurations from the N = 3 shell. Early indications that a large mixing with intruder 2p–2h and 4p–4h configuration is indeed necessary to help stabilizing the heavier fluorine isotopes came from Monte Carlo shell-model calculations with the SDPF-M interaction14.

### New results

In view of the recent experimental results, we have adjusted our model for 28F to the new data, producing the scenario D, intermediate between C and D of Fig. 1 of ref. 3. By simply adjusting the position of the d-wave resonant state, we get also an f-wave resonance at about the same energy of one of the states indicated in Fig. 2 of ref. 6. The phase-shifts are shown in Fig. 3 and the energy at which they cross 90 is the position of the resonant state. Using this updated model, we have performed new three-body (27F + n + n) calculations for 29F to investigate its ground-state properties (as in ref. 3) and, for the first time, the effect of configuration mixing on its dipole response. We obtain, as expected due to parity inversion, a ground state dominated by a $${(2{p}_{3/2})}^{2}$$ component of 57.5%, followed by 28.1% of $${(1{d}_{3/2})}^{2}$$ and a smaller 6% of $${(1{f}_{7/2})}^{2}$$. The dominant p-wave configuration leads to a non-negligible increase of the matter radius with respect to the 27F core, with ΔR = 0.192 fm, which may be an indication of moderate halo formation in 29F. Intruder components seem to favor also the dineutron configuration, which becomes more pronounced in the wave-function density as a direct consequence of the mixing. This is shown in Fig. 4a.

From the ground state, we have considered electric dipole transitions into the continuum. It is known that two-neutron halo nuclei, and in general weakly bound nuclei at the driplines, exhibit large E1 transition probabilities (and therefore large cross-sections in Coulomb excitation experiments)15,16,17. Since the halo structure is favored by p-wave intruder components, a measurement of the E1 strength could confirm whether 29F belongs to the island of inversion. Another way to pin down the content of p-waves is through neutron knock-out reactions. Our prediction using the three-body model constrained by the new data on 28F is shown in Fig. 4b. With the present model, the B(E1) distribution is strongly concentrated at low 27F + n + n energies. The integrated E1 strength up to 6 MeV is 1.59 e2fm2, which is a large fraction of the total cluster sum rule, that amounts to 1.92 e2fm2. Note the presence of a prominent peak at about 0.8 MeV above the three-body threshold, which exhausts about 50% of the total strength. This peak arises from a resonant state, which has been identified as an eigenstate of a resonance operator using the method proposed in ref. 18. The composition of this state yields 73% of (2p3/2)(1d3/2) components, so that the computed E1 distributions arises mainly from $${(2{p}_{3/2})}^{2}\to (2{p}_{3/2})(1{d}_{3/2})$$ and $${(1{d}_{3/2})}^{2}\to (2{p}_{3/2})(1{d}_{3/2})$$ transitions.

It is worth noting that the E1 probability is highly sensitive to the ground-state radius and configuration mixing. In particular, the same calculations using the scenario A (standard shell-model order) of ref. 3 yields a reduction of the total B(E1) by  ~40%. Such a difference should lead to very different experimental curves, making this observable a good probe for the configuration mixing and parity inversion in this region of the nuclear chart.

### Estimate of Relativistic Coulomb Excitation cross-section

As a final note, we can estimate the cross-section for Relativistic Coulomb Excitation (RCE) to the continuum, which in the semiclassical approach it is just proportional to the distribution of electromagnetic strength19,20, i.e.,

$$\frac{d\sigma }{d\varepsilon }={(\alpha {Z}_{{\rm{2}}})}^{2}\sum _{\pi \lambda \mu }{\left(\frac{\Delta E}{\hslash c}\right)}^{2\lambda -2}\frac{dB}{d\varepsilon }(\pi \lambda ;{\rm{gs}}\to {\rm{Cont.}})| {G}_{\pi \lambda \mu }({\beta }^{-1}){| }^{2}{g}_{\mu }(\xi )$$

where α is the fine structure constant, ΔE = S2n + ϵ, Gπλμ are the relativistic Winther–Adler functions and g is a function of the adiabaticity parameter ξ. Note that the proportionality factor, despite having no direct dependence on energy when λ = 1, is not a constant, because the function g has a decreasing dependence on energy through ξ (at this bombarding energy it amounts to a factor of 2 over the first 10 MeV’s). This modulates slightly the dipole distribution, by flattening its tail (see Fig. 4b). Choosing the same kinematic parameters of the experiment at RIBF6, i.e., a beam of 29F (Z1 = 9) at an energy of EB = 235 MeV/nucleon, impinging on a target of 208Pb (Z2 = 82), we get an estimate for the total cross-section of about 600 mb, integrated up to 15 MeV (with more than 90% of the cross-section concentrated below 5–6 MeV). As discussed in ref. 10, the nuclear contribution is not negligible, about half of the Coulomb contribution, in case of 31Ne. With a simple geometric model, the nuclear contribution can roughly be estimated as $$2\pi \sqrt{\frac{5}{3}}$$(R(27F) + R(208Pb))ΔR, resulting in  ≈200 mb. The semiclassical theory allows also to estimate the forward-peaked angular range for this process to about Δθ ~ 2Z1Z2e2/(RA1EB) ~ 0.2219. Our predictions may guide the interpretation of future experiments focused on the structure of 29F and the extent of the island of inversion. Let us also mention that a new method, alternative to the Equivalent Photon Method, has recently been proposed21 to extract B(E1) strengths, disentangling nuclear and coulomb contributions.

Slightly after our Perspective paper has been submitted, another paper was published in Physical Review Letters22 where new measurements of reaction cross-sections are reported for 27,29F on C at RIKEN. From the large value of the cross-section, out of the systematic trend for fluorine isotopes, the authors conclude that 29F is the heaviest Borromean two-neutron halo nucleus to date and explain their findings with a large occupancy of the p3/2 orbital, supported by shell-model calculations and compared with coupled-cluster calculations. This finding is fully consistent with our discussion in Table III of ref. 3, despite the suggested value of ΔR = 0.35 fm being somewhat larger than our predicted value for ΔR, thus pointing to a larger halo structure. The comparison of our three-body calculations with the shell-model results presented in ref. 22 shows a reasonable agreement.

## Outlook

In conclusion, we believe that new experiments by Revel and collaborators6 and by Bagchi and collaborators22 are clearing the field of any doubt on the theoretical interpretation on the structure of the interesting 29F nucleus. Isotopes sitting near the southern shore of the island of inversion are crucial for our understanding of nuclear stability in general and for the peculiar interplay of nuclear forces that determines the insurgence of a small p-wave halo, the shell-evolution and the disappearance of the N = 20 magic number in this region of the nuclear chart. We propose that Relativistic Coulomb Excitation experiments might lead to additional confirmations of these new findings.

## Data availability

There are no data associated with the present paper. The theoretical calculations plotted in Figs. 3 and 4 are available from the corresponding author upon request.

## References

1. Otsuka, T. et al. Magic numbers in exotic nuclei and spin-isospin properties of the NN interaction. Phys. Rev. Lett. 87, 082502 (2001).

2. Otsuka, T., Suzuki, T., Fujimoto, R., Grawe, R. H. & Akaishi, Y. Evolution of nuclear shells due to the tensor Force. Phys. Rev. Lett. 95, 232502 (2005).

3. Singh, J., Casal, J., Horiuchi, W., Fortunato, L. & Vitturi, A. Exploring two-neutron halo formation in the ground state of 29F within a three-body model. Phys. Rev. C. 101, 024310 (2020).

4. Christian, G. et al. Exploring the Low-Z Shore of the Island of Inversion at N = 19. Phys. Rev. Lett. 108, 032501 (2012).

5. Christian, G. et al. Spectroscopy of neutron-unbound 27,28F. Phys. Rev. C. 85, 034327 (2012).

6. Revel, A. et al. Extending the southern shore of the island of inversion to 28F. Phys. Rev. Lett. 124, 152502 (2020).

7. Michel, N., Li, J. G., Xu, F. R. & Zuo, W. Two-neutron halo structure of 31F. Phys. Rev. C 101, 031301 (2020).

8. Michel, N. private communication (2020).

9. Masui, H., Horiuchi, W. & Kimura, K. Two-neutron halo structure of 31F and a novel pairing antihalo effect. Phys. Rev. C 101, 041303 (2020).

10. Horiuchi, W., Suzuki, Y., Capel, P. & Baye, D. Probing the weakly-bound neutron orbit of 31Ne with total reaction and one-neutron removal cross sections. Phys. Rev. C. 81, 024606 (2010).

11. Doornenbal, P. et al. Low-Z shore of the “island of inversion” and the reduced neutron magicity toward 28O. Phys. Rev. C 95, 041301 (2017)..

12. Gaudefroy, L. et al. Direct mass measurements of 19B, 22C, 29F, 31Ne, 34Na and other light exotic nuclei. Phys. Rev. Lett. 109, 202503 (2012).

13. Macchiavelli, A. O. et al. Structure of 29F in the rotation-aligned coupling scheme of the particle-rotor model. Phys. Lett. B 775, 160–162 (2017).

14. Fallon, P. et al. Two-proton knockout from 32Mg: intruder amplitudes in 30Ne and implications for the binding of 29,31F. Phys. Rev. C 81, 041302 (2010).

15. Aumann, T. Low-energy dipole response of exotic nuclei. Eur. Phys. J. A 55, 234 (2019).

16. Nakamura, T. et al. Observation of strong low-lying E1 strength in the two-neutron halo nucleus 11Li. Phys. Rev. Lett. 96, 252502 (2006).

17. Aumann, T. & Nakamura, T. The electric dipole response of exotic nuclei. Phys. Scr. 2013, 014012 (2013).

18. Casal, J. & Gómez-Camacho, J. Identifying structures in the continuum: application to 16Be. Phys. Rev. C. 99, 014604 (2019).

19. Winther, A. & Alder, K. Relativistic Coulomb excitation. Nucl. Phys. A 319, 518–532 (1979).

20. Bertulani, C. & Baur, G. Electromagnetic processes in relativistic heavy ion collisions. Phys. Rep. 163, 299–408 (1988).

21. Moro, A. M., Lay, J. A. & Gómez-Camacho, J. Determining B(E1) distributions of weakly bound nuclei from breakup cross sections using Continuum Discretized Coupled Channels calculations. Application to 11Be. Preprintat https://arxiv.org/abs/2004.14612 (2020).

22. Bagchi, S. et al. Two-neutron halo is unveiled in 29F. Phys. Rev. Lett. 124, 222504 (2020).

## Acknowledgements

This work was supported by SID funds 2019 (Universitá degli Studi di Padova, Italy) under project No. CASA_SID19_01, by JSPS KAKENHI Grants No. 18K03635, No. 18H04569, and No. 19H05140 and by the collaborative research programs 2020, information initiative center, Hokkaido University.

## Author information

Authors

### Contributions

L.F. drafted this paper and calculated the Relativistic Coulomb Excitation cross-section; J.C. calculated the phase-shifts and made the three-body calculations. J.S., W.H., and A.V. contributed with ideas, discussions and revisions. W.H. also provided estimates for the nuclear contribution.

### Corresponding author

Correspondence to L. Fortunato.

## Ethics declarations

### Competing interests

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Rights and permissions

Reprints and Permissions

Fortunato, L., Casal, J., Horiuchi, W. et al. The 29F nucleus as a lighthouse on the coast of the island of inversion. Commun Phys 3, 132 (2020). https://doi.org/10.1038/s42005-020-00402-5

• Accepted:

• Published:

• DOI: https://doi.org/10.1038/s42005-020-00402-5