B(Z,N) = Z Δ(1H) + N Δ(n) − Δ(Z,N) from the atomic mass excesses Δ.
S2n(Z,N) = B(Z,N) − B(Z,N−2), S2p(Z,N) = B(Z,N) − B(Z−2,N).
Δ2n(Z,N) = S2n(Z,N) − S2n(Z,N+2): the drop in S2n just after N. At a magic N it is large; between closures in heavy nuclei it is 0.3 to 1.1 MeV (tin, lead), and in light nuclei it is larger throughout. Δ2p likewise along Z.
Δ(3)n(N) = ½(−1)N+1[B(N−1) − 2B(N) + B(N+1)], the three-point odd-even staggering. It is quoted as a pairing indicator at odd N only: there it is dominated by pairing, though it also contains smooth mass curvature, polarization and, near N = Z, neutron-proton correlations; at even N it also carries the spacing between the last filled and first empty orbital (Satula, Dobaczewski, Nazarewicz, PRL 81, 3599 (1998)), so it spikes at closed shells. The five-point difference smooths further but averages the two parities together, which is why it is not used here.
Along N = Z the Wigner energy adds extra binding, so Δ2n peaks there too (24Mg, 32S, 40Ca): an N = Z ridge is not by itself a shell gap.
Uncertainties: for Sn, Sp, S2n, S2p the AME2020 tabulated ones, which include the correlations between masses; for Δ2n, Δ2p, Δ(3) the mass uncertainties added in quadrature, treating the masses as independent. # marks a value that rests on at least one mass estimated from systematics rather than measured; on the chart those columns are drawn narrow and grey.