Data for 'Partial Tidal Disruption Events: The Elixir of Life'

This folder includes all the data for the paper titled: Partial Tidal Disruption Events: The Elixir of Life

Simulations were run with publicly available phantom code (https://github.com/danieljprice/phantom; https://phantomsph.readthedocs.io/en/latest/) and Kepler code (https://2sn.org/kepler/doc/Setup.html)

We used splash (https://github.com/danieljprice/splash.git) for making Figures 2, 7, 8, 9, 10, 12, 14, B2, F1.

We have three folders called s1_kepler, s3_kepler and s10_kepler. They contain the KEPLER models .wnd and .cnv files along with the dumps for ZAMS, MAMS and TAMS

We have saved the PHANTOM relaxed models in star_(mass_of_star)_(evolution_stage) along with mapping back evolution in KEPLER files (link.dat, .wnd, .cnv). We have also saved Makefile, star.in, star.setup, tde.comp files, linkg.

We have saved the PHANTOM simulation data in folders named star(mass_of_star)_(evolution_stage)_(beta_val). These include the .in, .tdeparams, ptok file, tde dump file used for analysis, KEPLER mapping back files (link.dat, .wnd, .cnv), .ev files. 

We have saved the KEPLER models that had the same mass as the remnants in star(mass_of_star)_(evolution_stage)_(beta_val)_kepler folders. We provide the generator file, .wnd, .cnv files. 

We have saved the KEPLER models that were generated by replacing the composition in the ptok files with the original star composition in the star_(mass_of_star)_(evolution_stage)_(beta_val)_strip

We have saved the KEPLER models that were run until the He4 in the central bin was same as remnant He4 in centre, in the star(mass_of_star)_(evolution_stage)_(beta_val)_he4_same

We have saved the PHANTOM simulation data for zero energy orbits run in star(mass_of_star)_(evolution_stage)_(beta_val)_zero folders. 

We have saved the geodesic files in the geodesic folder. It includes the .py files used to run jobs in parallel, and interpolation files. We have also saved the rp values used in the simulations run. 

We have saved the scripts used to run the KEPLER models (mapping back and stripping) in the python_code folder. We have also saved an EXCEL file which includes all the PHANTOM simulations data, along with the correct beta values shown in the paper. 

-----------------------------------------------------------------------------------------------------

For the figures from the publication, 

Figure 1 uses the data from s1_kepler, s3_kepler and s10_kepler folders. 

Figure 2 uses data files present in the star(mass_of_star)_(evolution_stage)_(beta_val) folders. The folders for 3 solar MID have sub-folders starting with fig2 which contain the PHANTOM dump files.  

Figure 3 uses the data from solar3_mid_beta(value) folders. We have included the .ev files which were used to make this plot. 

Figure 4 uses the mass of the remnant as function of time. We have included the dump_info files that can be used to make this plot. 

Figure 5 and 6 use the mass of the remnant. You can access this information from Table 1 and 2 of the paper. 

Figure 7 uses the remnant files from the star1, star3 and star10 MID models. 

Figure 8 uses the ptok files and tde dump files. These are all included in the folders. 

Figure 9 uses dump files which are included in a sub-folder in the star3_mid_beta18 folder. 

Figure 10 requires the dump files along with files including the binned rotation profile and the break-up velocities. All of this is present in the folders where .exact1 and .exact2 are the binned and break-up profiles.

Figure 11 can be generated by using the ptok files as they include binned rotation profile and radius.

Figure 12 uses data files from the star3_mid_beta18 folder. We have included the remnant files. The same folder includes vphi files for Figure 13. 

Figure 14 requires the dump files along with tde.comp files. All tde.comp files are present in star_(mass_of_star)_(evolution_stage) folders. Dumps are in star_(mass_of_star)_(evolution_stage)_beta(value)folders. 

Figure 15 uses the files provided in the folders ending with _he4_same. 

Figure 16 and 17 can be generated by using the .wnd files provided in the star_(mass_of_star)_(evolution_stage) and star_(mass_of_star)_(evolution_stage)_beta(value) folders. The same files can be used for Figure 18 and Figure 19. 

Figure 20 can be generated by using the files present in folders ending with _kepler. 

Figure 21 uses files from folders ending with _strip. 

Figure 22 uses ptok files which have been provided for all star10_tams_beta(value) folders

Figures A1, A2 can be generated using the ptok files and the original model files. 

Figure B2 uses files from sub-folder called cut_plot in the star3_mid_beta15 folder.

Figure E1 uses .wnd files present in star(mass_of_star)_(evolution_stage) folders. 

Figure E2 uses .wnd files from the star3_tams_beta25 (sub-folder, time_run).

Figure F1 uses files from zero and hyperbolic folders. 

-----------------------------------------------------------------------------------------------------

To load the .wnd and .cnv, use windata and conv modules from python/source code written by Alexander Heger (KEPLER).
To load files starting with tde_0* use splash or sarracen (https://sarracen.readthedocs.io/en/latest/quick-start.html). 
To load files starting with ptok (ASCII files), you can use splash or python. 
To load .ev (ASCII files), you can use splash or python. 
To load rem_* files or exact* files (ASCII) splash or python can be used. 
To load link.dat file (ASCII) one can use splash or python. 
To load KEPLER files containing @ one can use python. These are ASCII files. To convert any KEPLER files containing #, one would need to use Le-data filename to get it into ASCII form. 

-----------------------------------------------------------------------------------------------------
orbit.tdeparams files contain the initial position and velocities required for zero energy orbit. 
tde.tdeparams files contain the beta, pericentre and starting distance along with mass of SMBH (in code units). 

star_00000.tmp files are the relaxed stellar profiles that were disrupted by the SMBH. 
dump_info file is an ASCII file with information about the dumps such as the mass of the remnant, escape velocity, and maximum radius.