Data of the manuscript: Formation of Warm sub-Saturns through Structure-Coupled von Zeipel–Lidov–Kozai Migration

# Article Data Share This directory contains the data, plotting notebooks, generated figures, and representative simulation code used for the article Formation of Warm sub-Saturns through Structure-Coupled von Zeipel–Lidov–Kozai Migration. Authors: Jiapeng Gao and Gongjie Li ## Directory Layout ```text article_data_share/ ├── code/ ├── figures/ │ ├── article_figs/ │ ├── plot_data/ │ └── fig_01_...ipynb through fig_09_...ipynb └── raw_data/ ``` ## Environment The simulation code in `code/` requires a working MESA-r11532 installation and the corresponding MESA SDK (x86_64-linux-20190315) or equivalent toolchain. The plotting notebooks do not require rerunning MESA unless you want to regenerate raw simulation histories. ## Typical Use To rerun secular code without coupling, first use article_data_share/code/no_inflation_toi1842b_fast/planet_parameter_set.py to set up the initial parameters. Then run: ```bash cd article_data_share/code/no_inflation_toi1842b_fast ./run.sh ``` To rerun coupled code, you can set up initial parameters in planet_parameter_set.py and initialize.py. Then run the initialize.py and make and run. ```bash cd article_data_share/code/toi1842_10ME_surface_heating_1e-3_final python initialize.py ./mk ./rn ``` The population_10ME_surface_heating_1e-3_final code enables simulation with different stellar masses for population study. You can run the code with the following commands: ```bash cd article_data_share/code/population_10ME_surface_heating_1e-3_final python initialize.py ./mk ./rn ``` ## Data Layout `raw_data/` contains larger source products from the simulation and post-processing pipeline. - Per-run folders such as `toi1842_10ME_surface_heating_1e-3_tvp004_stellar_final_` contain saved `read_and_plot.py` outputs for selected final times. - `fig_06_coremass_histories/` contains the MESA `history.data` files used by `figures/fig_06_coremass_plot.ipynb` to regenerate `X_heatingrate.npy`, `Y_coremass.npy`, and `radius.npy`. `figures/plot_data/` contains compact, figure-ready inputs used directly by the notebooks. - `fig_01_Surface_10ME_time/`, `fig_02_final_scatter_with_marginals_two_panels_v2/`, etc. hold per-figure data subsets. - `fig_03_sub_saturns_plot_data.csv` stores the sampled sub-Saturn plotting data, including obliquity columns. - `PS_2026.03.25_20.56.30.csv` is the planet catalog input used by the sub-Saturn figures. The helper script `copy_read_and_plot_raw_data.py` copies selected `read_and_plot.py` outputs into `raw_data/`.