arXiv:2604.15931 Code and Data

This repository contains the computational material accompanying the manuscript arXiv:2604.15931. The calculations investigate the stability of the Starobinsky inflationary trajectory and the evolution of cosmological perturbations in the two-field scalar-tensor model studied in the manuscript. The principal computational resource is `starobinsky_modification_numerics-0.5.0.zip`. This archive contains the Python implementation used to generate and verify the numerical results reported in the current version of the manuscript. It calculates the homogeneous background trajectories, finite-time attractor diagnostics, complete scalar and tensor perturbation stability matrices, curvature and entropy power spectra, Bunch-Davies initialisation convergence tests, and entropy-transfer diagnostics. It also contains the scripts used to generate the manuscript figures and tables. The Python archive includes pinned software dependencies, debug and production configurations, restartable numerical workflows, automated validation tests, validation records, a documented HDF5 output format, and detailed instructions for reproducing the calculations. Readers should begin with README.md; the complete production workflow is described in docs/reproduce_manuscript.md. The repository also retains the Wolfram Mathematica scripts RS_Run_*.wls, the curvature-perturbation datasets R_S_Run_*_Data_Q.dat, the entropy-perturbation datasets R_S_Run_*_Data_s.dat, and the notebook R_S_Analysis_Equall_Masses.nb. These files belong to an earlier numerical workflow and are preserved for provenance and comparison. The Python archive is the recommended implementation for reproducing and independently checking the results presented in the current manuscript.