Neutron Economy and Freeze-out Dynamics in the Late-Time Cold r-Process

Standard r-process models often treat the freeze-out phase as a passive, uniform smoothing of the global abundance pattern. We demonstrate, however, that the final pattern is instead actively sculpted by a dynamic "neutron economy" that operates well after the free neutron source is exhausted. By analyzing cold neutron star merger outflows, we identify a distinct neutron exporter-importer relationship: the second-peak region ($A\sim 130$) transitions early into a $\beta$-decay dominated regime, acting as a net exporter, while the third-peak ($A\sim 195$) and Rare Earth ($A\sim 164$) regions remain capture-dominated neutron importers. Targeted simulations reveal that the final peak structures emerge from a competition between two opposing forces: a global "push" driven by neutron capture shifting mass toward heavier nuclei, and a local "pull" from beta delayed neutron emission returning material back toward lighter masses. We identify ${}^{136}Pd$, ${}^{135}Ag$, and ${}^{130}Ru$ as the primary exporters driving this economy, while third-peak fine structure is governed by $N=126$ isotones (Sm to Tb). Physically, this transport reflects a fundamental drive toward shell closure, where nuclei that overshoot the $N=82$ closed shell shed mass to migrate back toward this stable configuration, while nuclei in the third peak capture these redistributed neutrons to anchor themselves at the $N=126$ shell. Crucially, our models overproduce the odd-even staggering in the third peak compared to solar data. This discrepancy suggests current models likely overestimate the $\beta$-delayed neutron emission probabilities ($P_n$) for these $N=126$ nuclei, highlighting them as high-priority targets for theoretical refinement and future FRIB experiments.