## Abstract

The exotic, neutron-rich and weakly-bound isotope ^{29}F 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 (*S*_{n} = 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 interactions^{1,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 nucleus^{3}, 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 explored^{3}. 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 ^{28}F 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 ^{29}F back to stability, in a way similar to ^{6}He and ^{11}Li. 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 work^{7}, 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 ^{31}F. The comparative plots of one-body density distributions suggest that ^{31}F has a halo structure with neutrons occupying an orbital with mixed configurations: 2*p*_{3/2}, 1*d*_{3/2}, and 1*f*_{7/2}, while ^{27}F and ^{29}F have instead a smaller radius, less compatible with a halo. These calculations produce a ^{29}F-wave function where the *p*_{3/2} components amount to only a few percent^{8}, more akin to the standard scenario A above. A novel mechanism to form the halo structure in ^{31}F 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 beams^{6}, performed at the Radioactive Isotope Beam Factory (RIBF) of the RIKEN Nishina Center in Japan, indicates that ^{28}F exists as a resonant state at about 199 keV above the ^{27}F + 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 “…*p**f* shell intruder components play only a small role in the ground-state structure of ^{28}F…”. 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 measurements^{11} found a resonance at 1080(18) keV in ^{29}F, 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 1*d*_{5/2} orbital, while the 1/2^{+} state arises mainly from the coupling of this proton with the 2^{+} excitation of ^{28}O. 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 ^{29}F 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 data^{11,12} on ^{29}F were published, a theoretical study^{13} 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 interaction^{14}.

### New results

In view of the recent experimental results, we have adjusted our model for ^{28}F 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 (^{27}F + n + n) calculations for ^{29}F 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 ^{27}F core, with Δ*R* = 0.192 fm, which may be an indication of moderate halo formation in ^{29}F. 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 *E*1 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 *E*1 strength could confirm whether ^{29}F 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 ^{28}F is shown in Fig. 4b. With the present model, the *B*(*E*1) distribution is strongly concentrated at low ^{27}F + *n* + *n* energies. The integrated *E*1 strength up to 6 MeV is 1.59 e^{2}fm^{2}, which is a large fraction of the total cluster sum rule, that amounts to 1.92 e^{2}fm^{2}. 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 (2*p*_{3/2})(1*d*_{3/2}) components, so that the computed *E*1 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 *E*1 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*(*E*1) 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 strength^{19,20}, i.e.,

where *α* is the fine structure constant, Δ*E* = *S*_{2n} + *ϵ*, *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 RIBF^{6}, i.e., a beam of ^{29}F (*Z*_{1} = 9) at an energy of *E*_{B} = 235 MeV/nucleon, impinging on a target of ^{208}Pb (*Z*_{2} = 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 ^{31}Ne. With a simple geometric model, the nuclear contribution can roughly be estimated as \(2\pi \sqrt{\frac{5}{3}}\)(*R*(^{27}F) + *R*(^{208}Pb))Δ*R*, resulting in ≈200 mb. The semiclassical theory allows also to estimate the forward-peaked angular range for this process to about Δ*θ* ~ 2*Z*_{1}*Z*_{2}*e*^{2}/(*R**A*_{1}*E*_{B}) ~ 0.22^{∘}^{19}. Our predictions may guide the interpretation of future experiments focused on the structure of ^{29}F 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 proposed^{21} to extract *B*(*E*1) strengths, disentangling nuclear and coulomb contributions.

### Note added during revision

Slightly after our Perspective paper has been submitted, another paper was published in Physical Review Letters^{22} where new measurements of reaction cross-sections are reported for ^{27,29}F on C at RIKEN. From the large value of the cross-section, out of the systematic trend for fluorine isotopes, the authors conclude that ^{29}F is the heaviest Borromean two-neutron halo nucleus to date and explain their findings with a large occupancy of the *p*_{3/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 collaborators^{6} and by Bagchi and collaborators^{22} are clearing the field of any doubt on the theoretical interpretation on the structure of the interesting ^{29}F 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.

## References

Otsuka, T. et al. Magic numbers in exotic nuclei and spin-isospin properties of the NN interaction.

*Phys. Rev. Lett.***87**, 082502 (2001).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).Singh, J., Casal, J., Horiuchi, W., Fortunato, L. & Vitturi, A. Exploring two-neutron halo formation in the ground state of

^{29}F within a three-body model.*Phys. Rev. C.***101**, 024310 (2020).Christian, G. et al. Exploring the Low-

*Z*Shore of the Island of Inversion at*N*= 19.*Phys. Rev. Lett.***108**, 032501 (2012).Christian, G. et al. Spectroscopy of neutron-unbound

^{27,28}F.*Phys. Rev. C.***85**, 034327 (2012).Revel, A. et al. Extending the southern shore of the island of inversion to

^{28}F.*Phys. Rev. Lett.***124**, 152502 (2020).Michel, N., Li, J. G., Xu, F. R. & Zuo, W. Two-neutron halo structure of

^{31}F.*Phys. Rev. C***101**, 031301 (2020).Michel, N. private communication (2020).

Masui, H., Horiuchi, W. & Kimura, K. Two-neutron halo structure of

^{31}F and a novel pairing antihalo effect.*Phys. Rev. C***101**, 041303 (2020).Horiuchi, W., Suzuki, Y., Capel, P. & Baye, D. Probing the weakly-bound neutron orbit of

^{31}Ne with total reaction and one-neutron removal cross sections.*Phys. Rev. C.***81**, 024606 (2010).Doornenbal, P. et al. Low-Z shore of the “island of inversion” and the reduced neutron magicity toward

^{28}O.*Phys. Rev. C***95**, 041301 (2017)..Gaudefroy, L. et al. Direct mass measurements of

^{19}B,^{22}C,^{29}F,^{31}Ne,^{34}Na and other light exotic nuclei.*Phys. Rev. Lett.***109**, 202503 (2012).Macchiavelli, A. O. et al. Structure of

^{29}F in the rotation-aligned coupling scheme of the particle-rotor model.*Phys. Lett. B***775**, 160–162 (2017).Fallon, P. et al. Two-proton knockout from

^{32}Mg: intruder amplitudes in^{30}Ne and implications for the binding of^{29,31}F.*Phys. Rev. C***81**, 041302 (2010).Aumann, T. Low-energy dipole response of exotic nuclei.

*Eur. Phys. J. A***55**, 234 (2019).Nakamura, T. et al. Observation of strong low-lying E1 strength in the two-neutron halo nucleus

^{11}Li.*Phys. Rev. Lett.***96**, 252502 (2006).Aumann, T. & Nakamura, T. The electric dipole response of exotic nuclei.

*Phys. Scr.***2013**, 014012 (2013).Casal, J. & Gómez-Camacho, J. Identifying structures in the continuum: application to

^{16}Be.*Phys. Rev. C.***99**, 014604 (2019).Winther, A. & Alder, K. Relativistic Coulomb excitation.

*Nucl. Phys. A***319**, 518–532 (1979).Bertulani, C. & Baur, G. Electromagnetic processes in relativistic heavy ion collisions.

*Phys. Rep.***163**, 299–408 (1988).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).

Bagchi, S. et al. Two-neutron halo is unveiled in

^{29}F.*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 and Affiliations

### 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

## Ethics declarations

### Competing interests

The authors declare no competing interests.

## Additional information

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

## Rights and permissions

**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 license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license 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 license, visit http://creativecommons.org/licenses/by/4.0/.

## About this article

### Cite this article

Fortunato, L., Casal, J., Horiuchi, W. *et al.* The ^{29}F 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

Received:

Accepted:

Published:

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

## Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.