Halbarath/Jupiter_Fuzzy_Meier_2025: Zenodo release

This repository contains the data and scripts to create all figures (except those that are SPH simulation snapshots) in the publication Meier et. al., 2025. Each subfolder contains the necessary files to create one or two of the figures. Figure_01: Analytic calculations for the definition of embryo sizes, no data files Figure_02_and_03: Statistics of collisions in the N-body simulations. 2080 data files split into 4 groups for the 3 different inclinations and 2 different separations. The data structure of these files is defined by the python classes Load_MyLogcolFile and LogCollision in Figure_02_and_03.py Figure_05: Time series of mixing metrics for the collision simulations of the 3 cases at different resolutions. 9 data files. File names correspond to the naming convention used in the paper (e.g. file 512_timeseries.txt corresponds to run 5.1.2). The rows are time, M_{mix,\alpha}, M_{mix,\beta}, M_{mix,\gamma}, M_{mix,\delta,0.2}, M_{mix,\delta,0.4}, M_{mix,\delta,0.6}, M_{mix,\beta\gamma}, M_{mix,\epsilon}, A_{I} as described in Appendix A of the paper. Figure_08: Density profile of the initial model for Jupiter and the 4 chosen simulation results as described in the paper, as well as 3 profiles from interior modeling literature. 7 data files. target.model, CMS19+HG23.csv and out10967_layers_mrt_with_header.txt have headers describing their structure. The remaining files have names corresponding to the naming convention used in the paper (e.g. file 511.00100.profile.txt corresponds to run 5.1.1 timestep 100). It contains the following collumns: Radius in Earth radii, enclosed mass in Earth masses, density in cgs, water mass fraction, rock mass fraction, hydrogen-helium mass fraction, temperature in Kelvin. The line for Militzer et al. 2022 is contained directly in the python script, as it is read from the plot in their paper. Figure_09_and_10: Heavy element fraction plot of the profiles before and after the thermal evolution simulations. 18 data files. All files contain 3 columns: normalized radius, normalized mass and heavy element fraction. File names are descriptive. Figure_11: Plot of M_{mix,\beta} normalized by the total core mass as a function of the core-mantle-boundary density contrast. Data is contained in the python script as (x,x,..;y,y,...) arrays for all cases. Variable names correspond to the nomenclature in the paper (e.g. Values711_4 are the values of runs 7.1.1 to 7.1.4). Figure_12: Resolution study of the mixing metrics. 56 data files in directories named for the resolution, file names are descriptive. Each directory also contains the python file to create the tipsy input files. The data files contain these entries as single numbers time (always 0), M_{mix,\alpha}, M_{mix,\beta}, M_{mix,\gamma}, M_{mix,\delta,0.2}, M_{mix,\delta,0.4}, M_{mix,\delta,0.6}, M_{mix,\beta\gamma}, M_{mix,\epsilon}, A_{I} as described in Appendix A of the paper. Figure_13: Plot of M_{mix,\beta}, M_{mix,\delta,0.2} and M_{mix,\beta\gamma} normalized by the total core mass as a function of the resolution for the pure material simulations. Data is contained in the python script (arrays of 3 metrics with 4 entries for the different resolutions). Figure_14: Same as Figure_13, but for the simulations with 2 material bodies. Figure_16: Plot of M_{mix,\beta}, M_{mix,\delta,0.2} and M_{mix,\beta\gamma} normalized by the total core mass as a function of time for the demixing simulations with Jupiter-like bodies. 3 data files. File names correspond to the nomenclature in the paper (e.g. 541_timeseries.txt corresponds to run 5.4.1). The rows are time, M_{mix,\alpha}, M_{mix,\beta}, M_{mix,\gamma}, M_{mix,\delta,0.2}, M_{mix,\delta,0.4}, M_{mix,\delta,0.6}, M_{mix,\beta\gamma}, M_{mix,\epsilon}, A_{I} as described in Appendix A of the paper. Figure_17: Plot of M_{mix,\beta} normalized by the total core mass as a function of time for the demixing simulations with rocky bodies. 12 data files. File names correspond to the nomenclature in the paper (e.g. 771_timeseries.txt corresponds to run 7.7.1). The rows are time, M_{mix,\alpha}, M_{mix,\beta}, M_{mix,\gamma}, M_{mix,\delta,0.2}, M_{mix,\delta,0.4}, M_{mix,\delta,0.6}, M_{mix,\beta\gamma}, M_{mix,\epsilon}, A_{I} as described in Appendix A of the paper.