Soot Planet Code

Abstract

(a) List of the types and formats of files uploaded The repository includes three text files, one pdf document and two rich text files. Readme.pdf Planet single layer.rtf Planet multi layer.rft (b) Summary of the codes created We created two python codes to calculate mass-radius relations of model planets, either as a single layer of a uniform composition or as multiple layers of distinct chemical compositions. The "Planet single layer" code calculates the radius of a single-layer planet with specified mass and composition. The "Planet multi layer" code calculates the radius of a multi-layer planet with specified mass and composition. (c) Explanations of the uses for uploaded codes To use the "Planet single layer" code, specify an array of masses (in multiples of Earth mass) and a composition (in mass fractions of iron, silicate, soot, and water), and the code will calculate the corresponding array of radii (in multiples of Earth radius) of a single-layer planet with the specified composition. To use the "Planet multi layer" code, specify an array of masses (in multiples of Earth mass) and a composition (in mass fractions of iron, silicate, soot, and water), and the code will calculate the corresponding array of radii (in multiples of Earth radius) of a multi-layer planet consisting of four distinct layers – an iron core, a silicate mantle, a soot layer, and an outer water layer – with the specified composition. We executed the codes in Python (version 3.13.5) using Jupyter Notebook. All computations were performed on a 2022 MacBook Air M2 equipped with 16 GB of memory and running macOS Sonoma 14.7.6. The codes are used in Li, Bergin, Hirschmann, Blake, Ciesla, Kempton, Soot Planets instead of Water World, Astrophysical Journal Letters (under review), 2025.