Evolutionary Models of Binary, Triple and Quadruple Star Systems

mse.zip contains the code used for this project. Refer to the documentation in doc/doc.pdf for details about the code and logs. The data is contained in the remaining 12 zip files. Each data set of binary (1+1), triple (2+1) and quaruple (2+2 and 3+1) are contained in the named zipped files. The slow kick models have a Maxwelling balck hole natal kick velocity peaked at 10 km/s The fast kick models have a Maxwelling balck hole natal kick velocity peaked at 50 km/s The filename describes if the calculations include the effects of stellar flybys or not. Each data zip file contains a file name hdf5 with the initial conditions and a .txt text file with the indices of the completed systems. The data for each system is stored in a .pkl file in the subdirectories, the log contained details of the system at evolutionary points of interest (stellar type change, mass transfer, mergers, supernova etc.). Refer to the documentation for details of the logs. In addition, for fixed time steps additional information is saved for the slow kick models. The individual system .pkl files require mse.py require some of the tools from mse. This is most easily managed by 1)Download mse.zip. 2)Unzip and create a subdirectory named data 3)Download any of the dataset zip files and store them in data. 4)Each file can be opened in python with a snippet such as: import pickle from mse import Tools with open("test1.pkl",'rb') as file: data = pickle.load(file) The following quantities can be extracted from each individual system .pkl file: Log - data['log'] For the slow kick models the following is also available: Number of orbits - data['N_orbits_status']Number of bodies - data['N_bodies_status']Time - data['t_print']Mass - data['m_print']Core mass - data['mc_print']Stellar type - data['k_print']Radius - data['R_print']Core Radius - data['Rc_print']X location - data['X_print']Y location - data['Y_print']Z location - data['Z_print']Semi-major axis - data['a_print']Eccentricity - data['e_print']Effective Temparature - data['T_eff_print']Luminosity - data['L_print']Relative Incluination - data['rel_INCL_print']Spin - spin_frequency_print = data['spin_frequency_print']