The termination of the slow neutron capture process (s-process) in the Pb--Bi mass region plays a key role in shaping the abundances of the heaviest stable nuclei produced in stellar environments. In this work, we revisit the classical description of the s-process termination developed by Clayton and co-authors, incorporating modern experimental neutron capture cross sections and exploring the sensitivity of the results to variations in neutron irradiation parameters. The reaction network describing the Pb--Bi--Po termination cycle is reformulated as a system of Bateman equations and solved numerically using the burn-up matrix formalism combined with the Pad\'e approximation for the matrix exponential. We investigate the influence of updated $(n,\gamma)$ cross sections for Pb and Bi isotopes, as well as variations in the shape and scale parameters of the neutron exposure distribution, on the isotopic abundances, the timescale required to reach equilibrium, and the behavior at maximum neutron exposure. Our results confirm the robustness of the quasi-stationary equilibrium condition $\sigma_A N_A = \mathrm{const}$ under modern nuclear data, while demonstrating that updated cross sections lead to a noticeable redistribution of isotopic abundances and modify the timescale for reaching the equilibrium plateau. The contribution of the s-process termination to the total abundances of Pb and Bi isotopes is further examined using different neutron exposure distributions. We find that the resulting s-process abundances are largely insensitive to the neutron capture cross section of $^{209}\mathrm{Bi}$ over the investigated range, confirming that $^{208}\mathrm{Pb}$ effectively acts as the termination point of the s-process. Among the considered exposure models, the exponential neutron exposure distribution provides the best overall description of the final s-process abundances. These results are in good agreement with observational abundance constraints and demonstrate that the main qualitative features of the classical Clayton scheme remain valid, while quantitative predictions are significantly refined by modern nuclear data and irradiation models.

We examine cyclic nuclear reactions in the lead–bismuth (Pb–Bi) system near the s-process termination point. We present a numerical investigation of the isotopic evolution and decay heat generation in an extended 60-isotope nuclear reaction network under continuous 30 keV neutron irradiation (${10}^{14}$–${10}^{28}$ n cm⁻² s⁻¹) using the Chebyshev Rational Approximation Method (CRAM). The network accounts for 88 transitions, utilising a hybrid data approach that combines neutron capture cross-sections from EAF-2010 and TALYS with fundamental decay properties from the ENDF/B-VIII.0. Our simulations reveal two distinct evolutionary regimes. At moderate fluxes (${10}^{14}$–${10}^{20}$ n cm⁻² s⁻¹), the system establishes a steady cyclic loop driven by the $\alpha$-decay of ${}^{210}\mathrm{Po}$, successfully reproducing the s-process termination isotopic distribution (${}^{208}\mathrm{Pb}{\gg }^{209}\mathrm{Bi}{>}^{207}\mathrm{Pb}{>}^{206}\mathrm{Pb}$), characteristic of low-metallicity AGB stars. As the flux exceeds ${10}^{22}$ n cm⁻² s⁻¹, the classical balance breaks down, propelling mass flow toward heavier trace isotopes and suggesting a potential transition into the intermediate neutron capture (i-process) regime. Heat density calculations demonstrate that while the energy release of the core cycle plateaus near ${10}^5$ W cm⁻³, the extended chain drives an energy surge to over ${10}^8$ W cm⁻³ at ${10}^{24}$ n cm⁻² s⁻¹ before the system enters an unstable transient state at extreme fluxes.