The neutron drip line — the termination of the isotopic chain where neutron separation energies turn negative — is known experimentally only up to the lightest elements, and calcium (Z = 20) is the benchmark region where microscopic theory most sharply disagrees with itself. An exact pairing-plus-continuum calculation (Dassie and Id Betan, 2020) using the modified Richardson equations places the calcium two-neutron drip line at 60Ca with an isospin-independent parametrization but at 66Ca with an isospin-dependent one, and obtains a one-neutron drip line at 57Ca that conflicts with the experimentally discovered — and hence particle-bound — 59Ca; the authors attribute the conflict to a possible overestimate of the pairing gap in the hybrid odd-mass treatment, omit deformation, and call for mass measurements of 60Ca and heavier isotopes to discriminate the models. As of 2026, masses are measured only through 57Ca, 60Ca remains the heaviest identified calcium isotope (a single candidate 61Ca event exists), and theoretical predictions of the two-neutron drip line span 60Ca to 71Ca. The problem is to determine the last particle-bound calcium isotopes — the isotopes where the one- and two-neutron separation energies change sign — either by direct mass measurements of 58Ca and heavier isotopes or by a theoretical determination that models continuum coupling, pairing (including odd-mass staggering), and deformation and is consistent with all existing data.

Background

The neutron drip line is the boundary of nuclear existence: past it, the energy gained by adding a neutron no longer binds it, and the one-neutron separation energy SnS_n (energy to remove one neutron) or the two-neutron separation energy S2nS_{2n} becomes negative. Empirically the drip line is established only up to about neon (Z = 10); for every heavier element its location is a prediction, not a measurement. The calcium chain (Z = 20) is the best-studied frontier: precision mass spectrometry first pinned down the neutron-rich calcium masses and their implications for three-nucleon forces (Wienholtz et al., Nature 2013), later measurements reached 57Ca (Michimasa et al., PRL 2018, the heaviest calcium masses known to date), subshell closures were established at N = 32 and N = 34, and the isotopes 59Ca and 60Ca were discovered by projectile fragmentation at RIKEN (Tarasov et al., PRL 2018). Event-by-event in-flight identification implies 59Ca and 60Ca are particle-bound, but their masses have never been measured.

Pairing — the correlated addition of neutrons in time-reversed orbits — is the dominant correlation near the drip line, and it couples strongly to the continuum of unbound single-particle states. The source paper (Dassie and Id Betan, Phys. Rev. C 102, 064301 (2020)) addresses this by solving the pairing Hamiltonian exactly: the reduced (BCS-type) pairing problem is integrable through Richardson's equations, a set of nonlinear equations for collective pair energies, which the authors modify so that the continuum enters through the continuum single-particle level density. From the exact solution they compute binding energies, separation energies, and occupation probabilities for the calcium chain. Their results define the problem addressed here: two isospin-independent parametrizations (with 40Ca and 48Ca inert cores) both place the two-neutron drip line at 60Ca, while allowing the mean field and pairing strength to depend on isospin extends it to 66Ca; the one-neutron drip line, obtained from a hybrid treatment in which odd-isotope separation energies are assembled from the exact even-N solution combined with Bardeen-Cooper-Schrieffer (BCS) relations for the pairing gap Δ\Delta and Fermi level λ\lambda, falls at 57Ca. That last result disagrees with the observed particle stability of 59Ca, which the authors ascribe to a possible overestimation of the pairing gap — i.e. to the approximate odd-mass treatment itself. They also note that deformation of very neutron-rich calcium isotopes, expected to be important for the single-particle level ordering and continuum coupling, is absent from their spherical framework, and they state that mass measurements of 60Ca and heavier isotopes (e.g. at the Facility for Rare Isotope Beams, FRIB) are needed to discriminate between the competing predictions.

The wider theoretical landscape shows that the model dependence is not an artifact of one method. Ab initio coupled-cluster calculations — a wave-function method expanding the many-body state in particle-hole excitations of a reference state — with chiral two- and three-nucleon forces found 60Ca bound (Hagen et al., PRL 2012), and interactions derived from chiral effective field theory (EFT) have generally failed to bind calcium beyond neutron number N = 40. The Gamow shell model — a configuration-interaction shell model whose basis includes resonant and scattering (continuum) states — predicts 70Ca as the dripline nucleus and, like the source, 57Ca as the heaviest bound odd isotope (Li et al., PRC 2020). Bayesian model averaging over global mass models gives 68Ca about a 76% posterior probability of being bound and leaves 61Ca near even odds (Neufcourt et al., PRL 2019). Nuclear density functionals typically place the drip line near 70Ca, and deformed relativistic Hartree-Bogoliubov calculations in continuum — a self-consistent mean-field method treating pairing, deformation, and scattering states simultaneously — show that switching on deformation shifts predicted drip-line nuclei in the oxygen-to-calcium region (In et al., IJMPE 2021). Most recently, the chiral interaction N3LO_Texas — a combination of two- and three-nucleon potentials fitted at fourth and third chiral order — binds the chain out to a predicted two-neutron drip line at 71Ca (Hu et al., arXiv:2512.11723), widening rather than narrowing the spread. Any credible answer must also respect the measured spectra of the neutron-rich calcium isotopes, including spectroscopy of 56,58Ca that casts doubt on a doubly magic 60Ca (Chen et al., PLB 2023).

Problem Statement

Determine the location of the neutron drip line of the calcium isotopes. Concretely: (i) identify the heaviest calcium isotope with positive one-neutron separation energy SnS_n — currently known only to lie at or beyond 60^{60}Ca, with the binding of 61^{61}Ca undecided; and (ii) identify the heaviest calcium isotope with positive two-neutron separation energy S2nS_{2n} — currently known only to lie at or beyond 60^{60}Ca, with theoretical predictions spanning 60^{60}Ca to 71^{71}Ca. A complete answer fixes the sign of SnS_n and S2nS_{2n} for every calcium isotope from 58^{58}Ca past the drip line, thereby identifying the heaviest isotope with positive SnS_n and the heaviest with positive S2nS_{2n}, either by direct mass (separation-energy) measurements of 58^{58}Ca and heavier calcium isotopes, or by a theoretical determination that reproduces the existing data — the measured masses through 57^{57}Ca, the particle stability of 59,60^{59,60}Ca, and the measured spectra of 54^{54}-58^{58}Ca — and that carries a quantified uncertainty on its predicted SnS_n and S2nS_{2n}, so that the drip-line location it reports is stated with the range its uncertainties support. The single scientific objective uniting these components is the empirical or constrained-theoretical termination of the calcium isotopic chain.

References

  1. A. C. Dassie and R. M. Id Betan, 'Estimate of the location of the neutron drip line for calcium isotopes from an exact Hamiltonian with continuum pair correlations', Phys. Rev. C 102, 064301 (2020). DOI: 10.1103/PhysRevC.102.064301; arXiv:2012.00680, https://arxiv.org/abs/2012.00680
  2. O. B. Tarasov, D. S. Ahn, D. Bazin, et al., 'Discovery of 60Ca and Implications for the Stability of 70Ca', Phys. Rev. Lett. 121, 022501 (2018). DOI: 10.1103/PhysRevLett.121.022501, https://doi.org/10.1103/PhysRevLett.121.022501
  3. S. Michimasa, M. Kobayashi, Y. Kiyokawa, S. Ota, et al., 'Magic nature of neutrons in 54Ca: First mass measurements of 55-57Ca', Phys. Rev. Lett. 121, 022506 (2018). DOI: 10.1103/PhysRevLett.121.022506, https://doi.org/10.1103/PhysRevLett.121.022506
  4. F. Wienholtz, D. Beck, K. Blaum, et al., 'Masses of exotic calcium isotopes pin down nuclear forces', Nature 498, 346-349 (2013). DOI: 10.1038/nature12226, https://doi.org/10.1038/nature12226
  5. G. Hagen, M. Hjorth-Jensen, G. R. Jansen, R. Machleidt, T. Papenbrock, 'Evolution of shell structure in neutron-rich calcium isotopes', Phys. Rev. Lett. 109, 032502 (2012). DOI: 10.1103/PhysRevLett.109.032502, https://doi.org/10.1103/PhysRevLett.109.032502
  6. L. Neufcourt, Y. Cao, W. Nazarewicz, E. Olsen, F. Viens, 'Neutron drip line in the Ca region from Bayesian model averaging', Phys. Rev. Lett. 122, 062502 (2019). DOI: 10.1103/PhysRevLett.122.062502; arXiv:1901.07632, https://arxiv.org/abs/1901.07632
  7. J. G. Li, B. S. Hu, Q. Wu, Y. Gao, S. J. Dai, F. R. Xu, 'Neutron-rich calcium isotopes within realistic Gamow shell model calculations with continuum coupling', Phys. Rev. C 102, 034302 (2020). DOI: 10.1103/PhysRevC.102.034302; arXiv:2008.11948, https://arxiv.org/abs/2008.11948
  8. A. Magilligan, B. A. Brown, S. R. Stroberg, 'Data-driven configuration-interaction Hamiltonian extrapolation to 60Ca', Phys. Rev. C 104, L051302 (2021). DOI: 10.1103/PhysRevC.104.L051302, https://doi.org/10.1103/PhysRevC.104.L051302
  9. S. Chen, F. Browne, P. Doornenbal, et al., 'Level structures of 56,58Ca cast doubt on a doubly magic 60Ca', Phys. Lett. B 843, 138025 (2023). DOI: 10.1016/j.physletb.2023.138025; arXiv:2307.07077, https://arxiv.org/abs/2307.07077
  10. T. Koiwai, K. Wimmer, P. Doornenbal, et al., 'A first glimpse at the shell structure beyond 54Ca: Spectroscopy of 55K, 55Ca, and 57Ca', Phys. Lett. B 827, 136953 (2022). DOI: 10.1016/j.physletb.2022.136953, https://doi.org/10.1016/j.physletb.2022.136953
  11. E. J. In, P. Papakonstantinou, Y. Kim, S.-W. Hong, 'Neutron drip line in the deformed relativistic Hartree-Bogoliubov theory in continuum: Oxygen to Calcium', Int. J. Mod. Phys. E 30, 2150009 (2021). DOI: 10.1142/S0218301321500099; arXiv:2007.00893, https://arxiv.org/abs/2007.00893
  12. B. S. Hu, Z. H. Sun, G. Hagen, G. R. Jansen, T. Papenbrock, 'Ab initio computations from 78Ni towards 70Ca along neutron number N = 50', Phys. Lett. B 858, 139010 (2024). DOI: 10.1016/j.physletb.2024.139010; arXiv:2408.07856, https://arxiv.org/abs/2408.07856
  13. B. S. Hu, A. Ekström, C. Forssén, G. Hagen, W. G. Jiang, T. Miyagi, T. Papenbrock, 'The neutron dripline in calcium isotopes from a chiral interaction', arXiv:2512.11723 (2025, preprint), https://arxiv.org/abs/2512.11723
  14. O. B. Tarasov, B. M. Sherrill, A. C. Dombos, et al., 'Discovery of new isotopes in the fragmentation of 82Se and insights into their production', Phys. Rev. C 112, 034604 (2025). DOI: 10.1103/573p-7fjp, https://doi.org/10.1103/573p-7fjp
  15. W.-J. Ong, Z. Y. Xu, R. Grzywacz, et al., 'Beta-decay half-lives beyond 54Ca: A systematic survey of decay properties approaching the neutron dripline', Phys. Rev. Lett. 136, 092502 (2026). DOI: 10.1103/5yrq-s3g5, https://doi.org/10.1103/5yrq-s3g5
  16. F. Nowacki, A. Obertelli, A. Poves, 'The neutron-rich edge of the nuclear landscape: Experiment and theory', Prog. Part. Nucl. Phys. 120, 103866 (2021). DOI: 10.1016/j.ppnp.2021.103866; arXiv:2104.06238, https://arxiv.org/abs/2104.06238
  17. B. A. Brown, 'The nuclear shell model towards the drip lines', Physics 4, 525-547 (2022). DOI: 10.3390/physics4020035; arXiv:2204.06088, https://arxiv.org/abs/2204.06088