REBOUNDxPaper

Supplemental resources to accompany our implementation paper for new features in REBOUNDx versions 3.0.5 and 3.1.0. Includes mesa2txt tutorial, with scripts to import, isolate, and export data from MESA (Modules for Experiments in Stellar Astrophysics), for use with "Parameter Interpolation" in REBOUNDx.

Prerequisites

• NumPy
• Matplotlib
• REBOUND
• REBOUNDx
• py_mesa_reader (for mesa2txt.ipynb)

Usage

Figures from Baronett, et al. (2022)

This repository contains Jupyter Notebooks (e.g., fig1.ipynb) demonstrating new REBOUNDx features and code to generate some of the figures seen in our implementation paper, Baronett et al. (2022). We also include in subdirectories (e.g., /fig4/) sample REBOUNDx Python scripts used to generate the data for their respective figures.

MESA-REBOUNDx Parameter Interpolation

For those interested in incorporating stellar evolution data from MESA into REBOUND simulations, we recommend following the methodology laid out in mesa2txt.ipynb, to easily extract desired stellar properties from MESA output logs, for use with REBOUNDx's "Paramter Interpolation".

MESA Project Files

We've made available our MESA project files, including custom inlists, etc., under the mesa subfolder. See the subfolder's README for more details.

Acknowledgement

If you find any code here useful in your research, in particular the mesa2txt methodology, we would greatly appreciate a citation of our immplementation paper, Baronett et al. (2022), in your work.

Issues

If you have any questions, suggestions, or encounter any problems with the code in this repository, please let us know by opening an issue. We appreciate any feedback.