What do calcium and tin isotopes reveal about the invisible matter inside neutron stars?

A dense star that cannot be squeezed in a lab

A neutron star packs more mass than the Sun into a sphere only tens of kilometres wide. That density is far beyond anything a laboratory can create, so astronomers read it from the stars themselves. One such star, the heaviest known pulsar J0740+6620, now enters that reading at 2.08 plus or minus 0.07 solar masses. Because a star of that weight must not collapse into a black hole, the densest matter inside it has to be stiff enough to bear the load. Other measurements, including four NICER radius readings and a merger signal, keep that stiffness from becoming arbitrary.

A Bayesian calibration turns those observations into probability distributions rather than one fixed answer for nuclear matter. The pulsar mass acts as a ceiling on how soft the high-density part can be, while the radii and merger deformation act as checks on how much space such a star may occupy. Together they narrow the family of acceptable models, but they do not close it into a single curve. Each surviving model still carries quantified uncertainty, so the result is a range of possible dense-star behaviour with stated odds.

That probabilistic form matters because the same calculation must also stay connected to ordinary nuclear matter near saturation density. A free high-density guess could support heavy stars while contradicting what we know about nuclei in the laboratory. The protocol therefore passes the low-density information forward, so the dense core and the outer layers remain part of one physical description. In that way, the heaviest pulsar does not simply add a number. It anchors an otherwise unreachable region of matter to evidence we can actually measure.

Open-shell calcium and tin pull on the nuclear model

A probability-based nuclear calibration now uses open-shell calcium and tin nuclei to test how it handles unpaired particles. Those nuclei matter because their measured masses and sizes respond to pairing and to how nuclear energy changes with density. The calibration added four calcium isotopes, 42, 44, 46 and 50, plus three tin isotopes, 112, 120 and 124. Their binding energies and charge radii give the model a wider set of low-density anchors than closed-shell nuclei alone. A neutron-pairing signal from tin-120 also helps fix an added pairing strength in the calculation. That pairing term matters because open-shell nuclei can have unpaired neutrons, and their masses reveal how strongly those neutrons bind. The update pushed the binding energy at ordinary nuclear density to lower values after these isotopes were included. It also moved the symmetry-energy parameters upward only modestly. Here the symmetry energy is the extra cost of having more neutrons than protons in a nucleus. Its slope describes how quickly that neutron-excess penalty rises as matter gets denser, and the inferred slope stays low compared with many other estimates. Yet values above 50 MeV are no longer strongly excluded once the open-shell data enter the fit. The result is therefore a probability distribution over possible nuclear parameters, not one settled answer for dense matter. That spread matters because binding energies pull toward one region while neutron-rich observables pull differently from another. So the calcium and tin measurements tighten some parts of the model without erasing the tension between them.

The crust gains heavier ions while its bulk strength stays put

An updated statistical inference built from measurements of atomic nuclei pushes the probability range for ion charge in the neutron-star crust toward higher values. That shift reflects a larger value of the symmetry-energy term at low density, which measures the energy cost of an unequal neutron-proton mix. So the crust’s ions move toward heavier charges, not because the lattice is rebuilt from scratch but because the allowed proton fraction rises. The same inference leaves the crust’s bulk mechanical properties almost unchanged, including pressure and shear modulus.

The cell size that sets how those ions sit in the solid also stays close to its earlier value. Central density, gas density, and radius show no material change when the new nuclear inputs are added. Because the ions remain embedded in a lattice with nearly the same spacing, the crust’s resistance to shear does not follow the heavier charge upward. Even finite-size corrections to the ions leave that shear modulus essentially where it was.

The effect is therefore local in composition rather than global in structure. The equation of state and mass-radius relation remain close to the earlier neutron-star analysis and still fit NICER radius measurements. That separation matters, because a crust can carry different ions without forcing the star’s total size or pressure profile to move.

A softer low-density symmetry energy, with a stiffer high-density tail

A Bayesian calibration of a Skyrme model leaves a probability distribution for the slope of the symmetry energy that is relatively low near ordinary nuclear density. A Skyrme model is a mathematical recipe for nuclear energy, and this calibration draws on data from open-shell calcium and tin nuclei, other finite-nucleus observables, direct calculations of pure neutron matter, X-ray size measurements from NICER, and tidal limits from gravitational waves. The useful part is that correlations are carried forward instead of being reset at each step. A compact statistical summary lets other groups use those linked uncertainties directly, rather than assuming that every nuclear parameter varies on its own.

This inference turns a difficult many-parameter problem into a shared starting point. Open-shell calcium and tin isotopes add binding energies, charge radii, and a pairing-gap constraint, so the result does not rest on one nuclear family or one poorly described mass. The neutron-star step then links those nuclear changes to the star’s crust and core, while keeping the overall mass-radius curve compatible with NICER radius measurements.

The limitation is that finite-nucleus data constrain the model most tightly near ordinary nuclear density, while core densities are supplied by a flexible high-density extension. Within this calibration, the allowed high-density behavior tends to stiffen above ordinary nuclear density if the model is to remain compatible with NICER radius measurements and with a pulsar mass of about 2.08 solar masses used as the maximum-mass constraint. Gravitational-wave tidal constraints then limit how stiff that high-density region can become. The compact Gaussian summary also misses skewness and multimodality, especially for a surface parameter whose probability distribution is broad and unevenly shaped, so other groups using the summary need to remember that it is not the full shape of uncertainty.

Because the inference remains probabilistic, future radius and tidal measurements will test how much high-density nuclear matter must stiffen while keeping the low-density symmetry energy relatively soft.