Nuclear drip line

Chart of nuclides for carbon to fluorine. Decay modes:

   proton emission

   positron emission or electron capture

   stable isotope

   beta decay

   neutron emission

The nuclear drip line is the boundary beyond which atomic nuclei are unbound with respect to the emission of a proton or neutron.

An arbitrary combination of protons and neutrons does not necessarily yield a stable nucleus. One can think of moving up or to the right across the table of nuclides by adding a proton or a neutron, respectively, to a given nucleus. However, adding nucleons one at a time to a given nucleus will eventually lead to a newly formed nucleus that immediately decays by emitting a proton (or neutron). Colloquially speaking, the nucleon has leaked or dripped out of the nucleus, hence giving rise to the term drip line.

Drip lines are defined for protons and neutrons at the extreme of the proton-to-neutron ratio; at p:n ratios at or beyond the drip lines, no bound nuclei can exist. While the location of the proton drip line is well known for many elements, the location of the neutron drip line is only known for elements up to neon. ^[1]

General description

Nuclear stability is limited to those combinations of protons and neutrons described by the chart of the nuclides, also called the valley of stability. The boundaries of this valley are the neutron drip line on the neutron-rich side, and the proton drip line on the proton-rich side. ^[2] These limits exist because of particle decay, whereby an exothermic nuclear transition can occur by the emission of one or more nucleons (not to be confused with particle decay in particle physics). As such, the drip line may be defined as the boundary beyond which proton or neutron separation energy becomes negative, favoring the emission of a particle from a newly formed unbound system. ^[2]

Allowed transitions

When considering whether a specific nuclear transmutation, a reaction or a decay, is energetically allowed, one only needs to sum the masses of the initial nucleus/nuclei and subtract from that value the sum of the masses of the product particles. If the result, or Q-value, is positive, then the transmutation is allowed, or exothermic because it releases energy, and if the Q-value is a negative quantity, then it is endothermic as at least that much energy must be added to the system before the transmutation may proceed. For example, to determine if ¹²C, the most common isotope of carbon, can undergo proton emission to ¹¹B, one finds that about 16 MeV must be added to the system for this process to be allowed. ^[3] While Q-values can be used to describe any nuclear transmutation, for particle decay, the particle separation energy quantity S, is also used, and it is equivalent to the negative of the Q-value. In other words, the proton separation energy S_p indicates how much energy must be added to a given nucleus to remove a single proton. Thus, the particle drip lines defined the boundaries where the particle separation energy is less than or equal to zero, for which the spontaneous emission of that particle is energetically allowed. ^[4]

Although the location of the drip lines is well defined as the boundary beyond which particle separation energy becomes negative, the definition of what constitutes a nucleus or an unbound resonance is unclear. ^[2] Some known nuclei of light elements beyond the drip lines decay with lifetimes on the order of 10⁻²² seconds; this is sometimes defined to be a limit of nuclear existence because several fundamental nuclear processes (such as vibration and rotation) occur on this timescale. ^[4] For more massive nuclei, particle emission half-lives may be significantly longer due to a stronger Coulomb barrier and enable other transitions such as alpha and beta decay to instead occur. This renders unambiguous determination of the drip lines difficult, as nuclei with lifetimes long enough to be observed exist far longer than the timescale of particle emission and are most probably bound. ^[2] Consequently, particle-unbound nuclei are difficult to observe directly, and are instead identified through their decay energy. ^[4]

Nuclear structure origin of the drip lines

The energy of a nucleon in a nucleus is its rest mass energy minus a binding energy. In addition to this, there is an energy due to degeneracy: for instance, a nucleon with energy E₁ will be forced to a higher energy E₂ if all the lower energy states are filled. This is because nucleons are fermions and obey Fermi–Dirac statistics. The work done in putting this nucleon to a higher energy level results in a pressure, which is the degeneracy pressure. When the effective binding energy, or Fermi energy, reaches zero, ^[5] adding a nucleon of the same isospin to the nucleus is not possible, as the new nucleon would have a negative effective binding energy — i.e. it is more energetically favourable (system will have lowest overall energy) for the nucleon to be created outside the nucleus. This defines the particle drip point for that species.

One- and two-particle drip lines

In many cases, nuclides along the drip lines are not contiguous, but rather are separated by so-called one-particle and two-particle drip lines. This is a consequence of even and odd nucleon numbers affecting binding energy, as nuclides with even numbers of nucleons generally have a higher binding energy, and hence greater stability, than adjacent odd nuclei. These energy differences result in the one-particle drip line in an odd-Z or odd-N nuclide, for which prompt proton or neutron emission is energetically favorable in that nuclide and all other odd nuclides further outside the drip line. ^[5] However, the next even nuclide outside the one-particle drip line may still be particle stable if its two-particle separation energy is non-negative. This is possible because the two-particle separation energy is always greater than the one-particle separation energy, and a transition to a less stable odd nuclide is energetically forbidden. The two-particle drip line is thus defined where the two-particle separation energy becomes negative, and denotes the outermost boundary for particle stability of a species. ^[5]

The one- and two-neutron drip lines have been experimentally determined up to neon, though unbound odd-N isotopes are known or deduced through non-observance for every element up to magnesium. ^[2] For example, the last bound odd-N fluorine isotope is ²⁶F, ^[6] though the last bound even-N isotope is ³¹F. ^[1]

Nuclei near the drip lines are uncommon on Earth

Of the three types of naturally occurring radioactivities (α, β, and γ), only alpha decay is a type of decay resulting from the nuclear strong force. The other proton and neutron decays occurred much earlier in the life of the atomic species and before the earth was formed. Thus, alpha decay can be considered either a form of particle decay or, less frequently, as a special case of nuclear fission. The timescale for the nuclear strong force is much faster than that of the nuclear weak force or the electromagnetic force, so the lifetime of nuclei past the drip lines are typically on the order of nanoseconds or less. For alpha decay, the timescale can be much longer than for proton or neutron emission owing to the high Coulomb barrier seen by an alpha-cluster in a nucleus (the alpha particle must tunnel through the barrier). As a consequence, there are no naturally-occurring nuclei on Earth that undergo proton or neutron emission; however, such nuclei can be created, for example, in the laboratory with accelerators or naturally in stars. ^[7] The Facility for Rare Isotope Beams (FRIB) at Michigan State University came online in mid-2022 and is slated to create novel radioisotopes, which will be extracted in a beam and used for study. It uses a process of running a beam of relatively stable isotopes through a medium, which disrupts the nuclei and creating numerous novel nuclei, which are then extracted. ^[8]

Nucleosynthesis

Explosive astrophysical environments often have very large fluxes of high-energy nucleons that can be captured on seed nuclei. In these environments, radiative proton or neutron capture will occur much faster than beta decays, and as astrophysical environments with both large neutron fluxes and high-energy protons are unknown at present, the reaction flow will proceed away from beta-stability towards or up to either the neutron or proton drip lines, respectively. However, once a nucleus reaches a drip line, as we have seen, no more nucleons of that species can be added to the particular nucleus, and the nucleus must first undergo a beta decay before further nucleon captures can occur.

Photodisintegration

While the drip lines impose the ultimate boundaries for nucleosynthesis, in high-energy environments the burning pathway may be limited before the drip lines are reached by photodisintegration, where a high-energy gamma ray knocks a nucleon out of a nucleus. The same nucleus is subject both to a flux of nucleons and photons, so an equilibrium between neutron capture and photodisintegration is reached for nuclides with a sufficiently low neutron separation energy, particularly those near waiting points. ^[9]

As the photon bath will typically be described by a Planckian distribution, higher energy photons will be less abundant, and so photodisintegration will not become significant until the nucleon separation energy begins to approach zero towards the drip lines, where photodisintegration may be induced by lower energy gamma rays. At 10⁹ kelvin, the photon distribution is energetic enough to knock nucleons out of any nuclei that have particle separation energies less than 3 MeV, ^[10] but to know which nuclei exist in what abundances one must also consider the competing radiative captures.

As neutron captures can proceed in any energy regime, neutron photodisintegration is unimportant except at higher energies. However, as proton captures are inhibited by the Coulomb barrier, the cross sections for those charged-particle reactions at lower energies are greatly suppressed, and in the higher energy regimes where proton captures have a large probability to occur, there is often a competition between the proton capture and the photodisintegration that occurs in explosive hydrogen burning; but because the proton drip line is relatively much closer to the valley of beta-stability than is the neutron drip line, nucleosynthesis in some environments may proceed as far as either nucleon drip line.^{[ citation needed ]}

Waiting points and time scales

Once radiative capture can no longer proceed on a given nucleus, either from photodisintegration or the drip lines, further nuclear processing to higher mass must either bypass this nucleus by undergoing a reaction with a heavier nucleus such as ⁴He, or more often wait for the beta decay. Nuclear species where a significant fraction of the mass builds up during a particular nucleosynthesis episode are considered nuclear waiting points, since further processing by fast radiative captures is delayed.

As has been emphasized, the beta decays are the slowest processes occurring in explosive nucleosynthesis. From the nuclear physics side, explosive nucleosynthesis time scales are set simply by summing the beta decay half-lives involved, ^[11] since the time scale for other nuclear processes is negligible in comparison, although practically speaking this time scale is typically dominated by the sum of a handful of waiting point nuclear half lives.

The r-process

The rapid neutron capture process is believed to operate very close to the neutron drip line, though the astrophysical site of the r-process, while widely believed to take place in core-collapse supernovae, is unknown. While the neutron drip line is very poorly determined experimentally, and the exact reaction flow is not precisely known, various models predict that nuclei along the r-process path have a two-neutron separation energy (S_2n) of approximately 2 MeV. Beyond this point, stability is thought to rapidly decrease in the vicinity of the drip line, with beta decay occurring before further neutron capture. ^[12] In fact, the nuclear physics of extremely neutron-rich matter is a fairly new subject, and already has led to the discovery of the island of inversion and halo nuclei such as ¹¹Li, which has a very diffuse neutron skin leading to a total radius comparable to that of ²⁰⁸Pb.^{[ clarification needed ]} Thus, although the neutron drip line and the r-process are linked very closely in research, it is an unknown frontier awaiting future research, both from theory and experiment.

The rp-process

The rapid proton capture process in X-ray bursts runs at the proton drip line except near some photodisintegration waiting points. This includes the nuclei ²¹Mg, ³⁰S, ³⁴Ar, ³⁸Ca, ⁵⁶Ni, ⁶⁰Zn, ⁶⁴Ge, ⁶⁸Se, ⁷²Kr, ⁷⁶Sr, and ⁸⁰Zr. ^{[13] [14]}

One clear nuclear structure pattern that emerges is the importance of pairing, as one notices all the waiting points above are at nuclei with an even number of protons, and all but ²¹Mg also have an even number of neutrons. However, the waiting points will depend on the assumptions of the X-ray burst model, such as metallicity, accretion rate, and the hydrodynamics, along with the nuclear uncertainties, and as mentioned above, the exact definition of the waiting point may not be consistent from one study to the next. Although there are nuclear uncertainties, compared to other explosive nucleosynthesis processes, the rp-process is quite well experimentally constrained, as, for example, all the above waiting point nuclei have at the least been observed in the laboratory. Thus as the nuclear physics inputs can be found in the literature or data compilations, the Computational Infrastructure for Nuclear Astrophysics allows one to do post-processing calculations on various X-ray burst models, and define for oneself the criteria for the waiting point, as well as alter any nuclear parameters.

While the rp-process in X-ray bursts may have difficulty bypassing the ⁶⁴Ge waiting point, ^[14] certainly in X-ray pulsars where the rp-process is stable, instability toward alpha decay places an upper limit near A = 100 on the mass that can be reached through continuous burning. ^[15] The exact limit is a matter presently under investigation; ^104–109Te are known to undergo alpha decay whereas ¹⁰³Sb is proton-unbound. ^[6] Even before the limit near A = 100 is reached, the proton flux is thought to considerably decrease and thus slow down the rp-process, before low capture rate and a cycle of transmutations between isotopes of tin, antimony, and tellurium upon further proton capture terminate it altogether. ^[16] However, it has been shown that if there are episodes of cooling or mixing of previous ashes into the burning zone, material as heavy as ¹²⁶Xe can be created. ^[17]

Neutron stars

In neutron stars, neutron heavy nuclei are found as relativistic electrons penetrate the nuclei and produce inverse beta decay, wherein the electron combines with a proton in the nucleus to make a neutron and an electron-neutrino:

p

+

e−

→

n

+

ν e

As more and more neutrons are created in nuclei the energy levels for neutrons get filled up to an energy level equal to the rest mass of a neutron. At this point any electron penetrating a nucleus will create a neutron, which will "drip" out of the nucleus. At this point we have:^{[ citation needed ]}

${\displaystyle E_{\text{F}}^{n}=m_{n}c^{2}\,}$

And from this point onwards the equation

${\displaystyle E_{\text{F}}^{n}={\sqrt {\left(p_{\text{F}}^{n}\right)^{2}c^{2}+m_{n}^{2}c^{4}}}\,}$

applies, where p_Fⁿ is the Fermi momentum of the neutron. As we go deeper into the neutron star the free neutron density increases, and as the Fermi momentum increases with increasing density, the Fermi energy increases, so that energy levels lower than the top level reach neutron drip and more and more neutrons drip out of nuclei so that we get nuclei in a neutron fluid. Eventually all the neutrons drip out of nuclei and we have reached the neutron fluid interior of the neutron star.

Known values

Neutron drip line

The values of the neutron drip line are only known for the first ten elements, hydrogen to neon. ^[18] For oxygen (Z = 8), the maximal number of bound neutrons is 16, rendering ²⁴O the heaviest particle-bound oxygen isotope. ^[19] For neon (Z = 10), the maximal number of bound neutrons increases to 24 in the heaviest particle-stable isotope ³⁴Ne. The location of the neutron drip line for fluorine and neon was determined in 2017 by the non-observation of isotopes immediately beyond the drip line. The same experiment found that the heaviest bound isotope of the next element, sodium, is at least ³⁹Na. ^{[20] [21]} These were the first new discoveries along the neutron drip line in over twenty years. ^[1]

The neutron drip line is expected to diverge from the line of beta stability after calcium with an average neutron-to-proton ratio of 2.4. ^[2] Hence, is predicted that the neutron drip line will fall out of reach for elements beyond zinc (where the drip line is estimated around N = 60) or possibly zirconium (estimated N = 88), as no known experimental techniques are theoretically capable of creating the necessary imbalance of protons and neutrons in drip line isotopes of heavier elements. ^[2] Indeed, neutron-rich isotopes such as ⁴⁹S, ⁵²Cl, and ⁵³Ar that were calculated to lie beyond the drip line have been reported as bound in 2017–2019, indicating that the neutron drip line may lie even farther away from the beta-stability line than predicted. ^[22]

The table below lists the heaviest particle-bound isotope of the first ten elements. ^[23]

Z	Isotope
01	03H
02	08He
03	011Li
04	014Be
05	017B
06	022C
07	023N
08	024O
09	031F
10	034Ne

Not all lighter isotopes are bound. For example, ³⁹Na is bound, but ³⁸Na is unbound. ^[1] As another example, although ⁶He and ⁸He are bound, ⁵He and ⁷He are not.

Proton drip line

The general location of the proton drip line is well established. For all elements occurring naturally on earth and having an odd number of protons, at least one species with a proton separation energy less than zero has been experimentally observed. Up to germanium, the location of the drip line for many elements with an even number of protons is known, but none past that point are listed in the evaluated nuclear data. There are a few exceptional cases where, due to nuclear pairing, there are some particle-bound species outside the drip line, such as ⁸B and ¹⁷⁸Au.^{[ verification needed ]} One may also note that nearing the magic numbers, the drip line is less understood. A compilation of the first unbound nuclei known to lie beyond the proton drip line is given below, with the number of protons, Z and the corresponding isotopes, taken from the National Nuclear Data Center. ^[24]

Z	Isotope
02	02He
03	05Li
04	06Be
05	07B, 09B
06	08C
07	11N
08	12O
09	16F
10	16Ne
11	19Na
12	19Mg
13	21Al
15	25P
17	30Cl
18	30Ar [25]
19	34K
21	39Sc
22	38Ti [26]
23	42V
25	45Mn
27	50Co
29	55Cu
30	54Zn [27]
31	59Ga
32	58Ge
33	65As
35	69Br
37	73Rb
39	77Y
41	81Nb
43	85Tc
45	89Rh
47	93Ag
49	97In
51	105Sb
53	110I
55	115Cs
57	119La
59	123Pr
61	128Pm
63	134Eu
65	139Tb
67	145Ho
69	149Tm
71	155Lu
73	159Ta
75	165Re
77	171Ir
79	175Au, 177Au
81	181Tl
83	189Bi
85	195At
87	201Fr
89	207Ac
91	214Pa
93	219Np [28]

See also

• Extended periodic table
• Table of nuclides
• Radioactive decay

Further reading

• Blumenfeld, Yorick (16 November 2022). "Probing the Limits of Nuclear Existence". Physics Magazine. 15. American Physical Society (APS): 177. Retrieved 1 March 2024.

References

1. Tarasov, O.B. (2017). "Production of very neutron rich isotopes: What should we know?".
2. Thoennessen, M. (2004). "Reaching the limits of nuclear stability" (PDF). Reports on Progress in Physics. 67 (7): 1187–1232. Bibcode:2004RPPh...67.1187T. doi:10.1088/0034-4885/67/7/R04. S2CID   250790169.
3. Wang, M.; Audi, G.; Kondev, F. G.; Huang, W. J.; Naimi, S.; Xu, X. (2017). "The AME2016 atomic mass evaluation (II). Tables, graphs, and references" (PDF). Chinese Physics C. 41 (3): 030003-1–030003-442. doi:10.1088/1674-1137/41/3/030003.
4. Thoennessen, M. (2016). The Discovery of Isotopes: A Complete Compilation. Springer. pp. 275–292. doi:10.1007/978-3-319-31763-2. ISBN   978-3-319-31761-8. LCCN   2016935977.
5. Smolańczuk, R.; Dobaczewski, J. (1993). "Particle-drip lines from the Hartree-Fock-Bogoliubov theory with Skyrme interaction". Physical Review C . 48 (5): R2166–R2169. arXiv: nucl-th/9307023v1 . Bibcode:1993PhRvC..48.2166S. doi:10.1103/PhysRevC.48.R2166. PMID   9969127. S2CID   12117057.
6. Audi, G.; Kondev, F. G.; Wang, M.; Huang, W. J.; Naimi, S. (2017). "The NUBASE2016 evaluation of nuclear properties" (PDF). Chinese Physics C. 41 (3): 030001. Bibcode:2017ChPhC..41c0001A. doi:10.1088/1674-1137/41/3/030001.
7. Mueller, Alex C.; Sherrill, Bradley M. (1993). "Nuclei at the Limits of Particle Stability" (PDF). Annual Review of Nuclear and Particle Science. 43: 529–583. doi:10.1146/annurev.ns.43.120193.002525.
8. "Biggest expansion of known chemical universe targeted by FRIB nuclear facility". 16 February 2021.
9. Zhao, B.; Zhang, S. Q. (2019). "The r-process with the Newly Developed High-precision Mass Model WS4". The Astrophysical Journal. 874 (5): 5. doi: 10.3847/1538-4357/ab0702 .
10. Thielemann, Friedrich-Karl; Kratz, Karl-Ludwig; Pfeiffer, Bernd; Rauscher, Thomas; et al. (1994). "Astrophysics and nuclei far from stability". Nuclear Physics A. 570 (1–2): 329. Bibcode:1994NuPhA.570..329T. doi:10.1016/0375-9474(94)90299-2.
11. van Wormer, L.; Goerres, J.; Iliadis, C.; Wiescher, M.; et al. (1994). "Reaction rates and reaction sequences in the rp-process". The Astrophysical Journal. 432: 326. Bibcode:1994ApJ...432..326V. doi: 10.1086/174572 .
12. Wang, R.; Chen, L.W. (2015). "Positioning the neutron drip line and the r-process paths in the nuclear landscape". Physical Review C. 92 (3): 031303–1–031303–5. arXiv: 1410.2498 . Bibcode:2015PhRvC..92c1303W. doi:10.1103/PhysRevC.92.031303. S2CID   59020556.
13. Koike, O.; Hashimoto, M.; Arai, K.; Wanajo, S. (1999). "Rapid proton capture on accreting neutron stars – effects of uncertainty in the nuclear process". Astronomy and Astrophysics. 342: 464. Bibcode:1999A&A...342..464K.
14. Fisker, Jacob Lund; Schatz, Hendrik; Thielemann, Friedrich-Karl (2008). "Explosive Hydrogen Burning during Type I X-Ray Bursts". The Astrophysical Journal Supplement Series. 174 (1): 261. arXiv: astro-ph/0703311 . Bibcode:2008ApJS..174..261F. doi:10.1086/521104. S2CID   119342620.
15. Schatz, H.; A. Aprahamian; V. Barnard; L. Bildsten; et al. (April 2001). "End Point of the rp Process on Accreting Neutron Stars" (subscription required). Physical Review Letters. 86 (16): 3471–3474. arXiv: astro-ph/0102418 . Bibcode:2001PhRvL..86.3471S. doi:10.1103/PhysRevLett.86.3471. PMID   11328001. S2CID   46148449 . Retrieved 2006-08-24.
16. Lahiri, S.; Gangopadhyay, G. (2012). "Endpoint of rp process using relativistic mean field approach and a new mass formula". International Journal of Modern Physics E. 21 (8): 1250074. arXiv: 1207.2924 . Bibcode:2012IJMPE..2150074L. doi:10.1142/S0218301312500747. S2CID   119259433.
17. Koike, Osamu; Hashimoto, Masa-aki; Kuromizu, Reiko; Fujimoto, Shin-ichirou (2004). "Final Products of the rp-Process on Accreting Neutron Stars". The Astrophysical Journal. 603 (1): 242–251. Bibcode:2004ApJ...603..242K. doi:10.1086/381354. S2CID   121805380.
18. "Three First-ever Atomic Nuclei Created; New Super-heavy Aluminum Isotopes May Exist". Sciencedaily.com. 2007-10-27. Retrieved 2010-04-06.
19. "Nuclear Physicists Examine Oxygen's Limits". Sciencedaily.com. 2007-09-18. Retrieved 2010-04-06.
20. Ahn, D.S.; et al. (2018). New isotope of ³⁹Na and the neutron dripline of neon isotopes using a 345 MeV/nucleon ⁴⁸Ca beam (Report). RIKEN Accelerator Progress Reports. Vol. 51. p. 82.
21. Ahn, D. S.; Amano, J.; Baba, H.; Fukuda, N.; Geissel, H.; Inabe, N.; Ishikawa, S.; Iwasa, N.; Komatsubara, T.; Kubo, T.; Kusaka, K.; Morrissey, D. J.; Nakamura, T.; Ohtake, M.; Otsu, H. (2022-11-14). "Discovery of ³⁹Na". Physical Review Letters. 129 (21): 212502. Bibcode:2022PhRvL.129u2502A. doi: 10.1103/PhysRevLett.129.212502 . PMID   36461972.
22. Neufcourt, L.; Cao, Y.; Nazarewicz, W.; Olsen, E.; Viens, F. (2019). "Neutron drip line in the Ca region from Bayesian model averaging". Physical Review Letters. 122 (6): 062502–1–062502–6. arXiv: 1901.07632 . Bibcode:2019PhRvL.122f2502N. doi:10.1103/PhysRevLett.122.062502. PMID   30822058. S2CID   73508148.
23. "Interactive Chart of Nuclei". www.nndc.bnl.gov. Archived from the original on 2005-12-23.
24. "National Nuclear Data Center" . Retrieved 2010-04-13.
25. Mukha, I.; et al. (2018). "Deep excursion beyond the proton dripline. I. Argon and chlorine isotope chains". Physical Review C. 98 (6): 064308–1–064308–13. arXiv: 1803.10951 . Bibcode:2018PhRvC..98f4308M. doi:10.1103/PhysRevC.98.064308. S2CID   119384311.
26. Meierfrankenfeld, D.; Bury, A.; Thoennessen, M. (2011). "Discovery of scandium, titanium, mercury, and einsteinium isotopes". Atomic Data and Nuclear Data Tables . 97 (2): 134–151. arXiv: 1003.5128 . Bibcode:2011ADNDT..97..134M. doi:10.1016/j.adt.2010.11.001. S2CID   97263560.
27. Gross, J.L.; Claes, J.; Kathawa, J.; Thoennessen, M. (2012). "Discovery of zinc, selenium, bromine, and neodymium isotopes". Atomic Data and Nuclear Data Tables . 98 (2): 75–94. arXiv: 1012.2027 . Bibcode:2012ADNDT..98...75G. doi:10.1016/j.adt.2011.12.001. S2CID   67813061.
28. Zhang, Z. Y.; Gan, Z. G.; Yang, H. B.; et al. (2019). "New isotope ²²⁰Np: Probing the robustness of the N = 126 shell closure in neptunium". Physical Review Letters . 122 (19): 192503. Bibcode:2019PhRvL.122s2503Z. doi:10.1103/PhysRevLett.122.192503. PMID   31144958. S2CID   169038981.