Saint-Maurice, commander of the Theban legion
This collects a few simple functions or scripts that can be useful for studying accreting (forming) planets.
-
RpTeff_from_Mdot_Mp_Popsynthfit.gnu,Rp_from_Mdot_Mp_Popsynthfit.{pro,py}: Planet radius and Teff as a function of accretion rate and mass, from the population syntheses of Mordasini et al. (2012), as fit by Aoyama, Marleau, Mordasini & Ikoma (2020). The files are organised by language (gnuplot(4.2/5),idl(8.7),python(3.7); other versions might/should work too).
New (13.02.2021): The Teff fit and test functions have been added to all three languages. -
accretionflowproperties: Computes approximate properties of the accretion flow onto a forming gas giant, including, at the shock, the density, temperature, velocity as a function of radius (distance from the planet) and time since the beginning of the fall. Assumes spherical symmetry. Based on the (semi-)analytical and simulation results of Marleau et al. (2019b). This submodule can be found under https://github.com/gabrielastro/accretionflowproperties. -
Upcoming, when time permits: the relations between accretion luminosity Lacc and line luminosity Lline from our Marleau & Aoyama (2022) Research Note "Planetary Line-to-accretion Luminosity Scaling Relations: Extrapolating to Higher-order Hydrogen Lines", complementing the Lacc–Lline scalings in Aoyama et al. (2021, ApJL). The explicit relations given there are already available in
pythonin Tomas Stolker'sspeciestoolkit; see the tutorial inspeciesanalysing the Paschen beta emission from GQ Lup B as an example.
Author: Gabriel-Dominique Marleau (c) 2023
Collaborator: Yuhiko Aoyama
Contact: uni-tuebingen.de with gabriel.marleau and @ in front
Note that the free, open-source Gnu Data Language (GDL) and also the closed-sourced but free Fawlty Language FL exist as drop-in replacements for IDL!
Comments and requests welcome.