SUPPLEMENTARY INFORMATION (ONLINE RESOURCES)

Article title: The Trap Number T: A Dimensionless Criterion for Boulder-Boulder
  Gravitational Dynamics on Near-Earth Rubble-Pile Asteroid Surfaces
Journal: [target journal - to be finalised]
Author: [Author Name]
Affiliation: [Department, Institution, City, Country]
Corresponding author e-mail: [e-mail address]

CONTENTS
--------
Data/
  Supplementary_Data_TrapNumber.xlsx   Parameters, Monte Carlo, Hessian, and d_eff
                                       sensitivity results, with live formulas.
  boulder_catalogue_Bennu.csv          Synthetic boulder catalogue (Bennu).
  boulder_catalogue_Ryugu.csv          Synthetic boulder catalogue (Ryugu).
  boulder_catalogue_Itokawa.csv        Synthetic boulder catalogue (Itokawa).
  asteroid_params_Bennu.csv            Adopted physical parameters (Bennu).
  asteroid_params_Ryugu.csv            Adopted physical parameters (Ryugu).
  asteroid_params_Itokawa.csv          Adopted physical parameters (Itokawa).

Code/
  code_trap_number.py        Core trap-number computation (surface map).
  code_monte_carlo.py        Monte Carlo statistics.
  code_hessian_stability.py  Hessian saddle-to-minimum analysis.
  code_deff_sensitivity.py   Sensitivity to the interaction scale d_eff.
  code_generate_figures.py   Regenerates Figures 1-6.
  TrapNumber_Derivations.wl  Symbolic verification of the derivations (Mathematica).

Derivations/
  Supplementary_Derivations.pdf  Step-by-step derivation of all equations, with diagrams.

REPRODUCIBILITY
---------------
Run the Python files with numpy and matplotlib installed. code_monte_carlo.py,
code_hessian_stability.py, and code_deff_sensitivity.py import from code_trap_number.py.
code_generate_figures.py reproduces all six manuscript figures.

REAL-DATA ANALYSIS (added)
--------------------------
Code/code_real_ryugu_trap_number.py   Computes the trap number directly from the
  real Ryugu boulder catalogue of Ray et al. (2026). Download that catalogue from
  the figshare deposit cited in the paper (see Boulder_Data_Sources.md), place
  Ryugu_Boulder_Attributes.csv beside the script, and run it. It reproduces the
  real-data result T_max ~= 0.5 (0.34-0.67) at Otohime. The Ray et al. catalogue is
  the authors' dataset and is cited, not redistributed here.
Code/run_trap_number_real_catalogue.py  Generic version for any lon/lat/size CSV.
Boulder_Data_Sources.md               Sources and links for the real catalogues.
