iADD-MESA-Binary 0.1.0

Library for reading and plotting data from iADD's MESA Binary tool

This README is just a copy of what is found on iADD_MESA_Binary's GitHub page: https://github.com/njhurst02/iADD-Mesa-Binary-Package.


Background


Package for reading out and plotting data from MESA binary simulations, designed in tandem with MESA Binary from iadd.astro.illinois.edu. Users are intended to run MESA binary simulations and then create plots as desired with this Package. Please see this page for more information and the MESA binary simulator.

Most of the functions produce animations (at the time of writing, only the kippenhahn plot is a static plot), giving users the ability to watch binary stellar evolution through multiple properties. Users can select the desired FPS of the movie and for specific plots (like abundances), threshold values must be specified. Movies are designed to be formatted such that they are accessible to all; colorblind friendly, large font sizes, etc. However, if you wish to see anything added, please contact us at the information found at the bottom of this page.



Installation & Usage


To install this module, use pip:

pip install iadd-mesa-binary

iADD MESA Binary does not currently provide general methods for reading out data. This may change in the future. Until then, please check out the methods in TULIPS or mesa_reader (MesaData). These packages are designed for single star data but with some creativity, they can work for data from binaries.


Additional Dependencies


Users will need to have the following installed for some or all of the plots. Clicking on the hyperlink takes you to their installation pages, most of which are simply a pip install command.

mesaPlot
mesa_reader
TULIPS (needed for abundances plot)
Latex (needed for text formatting, optional, not a pip install)



Functions


Behind all of these functions are auxiliary helper functions not listed here. They are however included in the download above, if you desire to look at them.

fps is the desired frames per second of the generated movie. Generally, 8-16 fps makes for good movies. Must be an int. PM is a switch for point mass evolution. True means, where applicable, only plot the evolution of the donor. Must be a boolean.

Modeling the system with a point mass companion results in some strange behavior. MESA treats them as a point mass with mass M and records no other properties, like radius or luminosity. It is thus impossible to plot certain graphs normally or at all. In the Roche Lobe geometry plot, the companion star will not display as the plotted circle has 0 radius. I've just set the value as 0 wherever relevant to avoid unnecessarily reading out blank data. The abundances and kippenhahn plots will not work at all for the companion but will behave as normal for the donor. The mass transfer and hertzsprung russel plots will function as normal, just without the companion star.

Additonally, when plotting as a point mass (PM=True), the binary_logs2 file is not used.

plot_Mass_Transfer(fps, PM):
plot_Roche_Lobe(fps, PM):
plot_Hertzsprung_Russel(fps, observers, PM):
Observers is a boolean; set to True, it will slightly modify the HR diagram to include spectral classes and a gradiented colorbar at the bottom. Because this has to be generated for every frame, making an observers HR diagram can take a few minutes.

plot_Abundances(fps, threshold, PM):
If an element is below the given threshold value, it is removed from the plot. 0.001 is typically very reasonable. Must be given as an int.

plot_Kippenhahn(PM):
Unlike the rest of the above functions, plot_Kippenhahn produces a picture, not a movie.

get_Binary_Data(PM):
Returns basic stellar parameters with flexibility for point mass evolution. If PM == False, get_Binary_Data returns (Age, a, M1, M2, Total_M, R1, R2, T1, T2, Lum1, Lum2, period) and if PM == True, get_Binary_Data returns (Age, a, M1, M2, Total_M, R1, T1, Lum1, period).

get_Data(parameter):
Given a parameter, returns data. If the requested data is in a history or binary_history file, data will be returned as a numpy array. If the requested data is in a profile file, user will be asked if they want data from a specific profile or all profiles. For one specific profile, data will be returned as a numpy array. For all profiles, data will be returned as a list of lists (due to the varying length of data from profiles).



Personal Use


If you wish to use this package for simulations not run with MESA_binary, the simplest way to do so is to ensure that: the logs for star 1 are called binary_logs1, the logs for star 2 are called binary_logs2, and the binary history is called binary_history.data. The get_Data function is only programmed to work with the possible outputs found in MESA but could be easily modified. If there are requests for more flexibility with the data reading, please contact me (information below).

Given the source code, the methods used here could easily be modified for single star evolution. If you wish for that to be a feature by default, please contact me (again, information below).

To streamline the process of aligning timesteps, this code assumes that history_interval and profile_interval are equal for star 1, star 2, and the binary data. On iADD's MESA Binary, this is done by default. If you want to use this code for simulations run on your own machine, make sure history_interval and profile_interval are equal for both stars and the binary data.

For the kippenhahn plot, you must edit profile_columns.list to include gradr (radiative gradient) and grada (adiabatic gradient). This is done by simply removing the ! before their entries. For the abundance plot to work, you must edit profiles_columns.list to include the elements/isotopes you want. Again, this is done by simply removing the ! before their entries.

For questions, comments, and concerns about this code, please contact me at [e-mail].