Initial conditions for the article "Survival of satellites during the migration of a Hot Jupiter"

This repository contains the initial conditions to reproduce the results of the article: "Survival of satellites during the migration of a Hot Jupiter".

These simulations were done with Posidonius (https://github.com/marblestation/posidonius/), version commit b3c5c71 (20/02/2022). To reproduce the results corresponding to the Equilibrium tide section (4.1), the two lines 20 & 21 in src/effects/tides/constant_time_lag.rs have to be commented and the code should be re-compiled before execution. By default, the dynamical tide is considered in both star and planet. 

The data file for the planet (evolution of the structure, following Leconte & Chabrier, 2013 [1]), including the dynamical tide dissipation as described in Section 2.3.3 is called Jupiter.dat. 

The initial conditions cover the following cases:

a disk density of $\Sigma_0 = $1000 g/cm$^2$

disk lifetimes of 5'000, 7'000, 10'000, 12'500, 20'000, 35'000, 50'000 yr

initial semi-major axis of the satellite of 0.3, 0.4, 0.6, 0.8, 1.0 $a_\mathrm{Io}$ (where $a_\mathrm{Io}$ is the distance between Jupiter and Io)

masses of the satellite of 1 $M_\mathrm{Io}$ and 10 $M_\mathrm{Io}$ (where $M_\mathrm{Io}$ is the mass of Io)

dissipations in the satellite of: $0$, $10^{-4}$, $10^{-2}$, $10^{-1}$, $1.0\times \Delta t_{\oplus}$ (where $\Delta t_\oplus$ is the time lag of the Earth)

initial stellar rotation periods of 1~day, 3~days and 8~days

The names of the initial conditions files contained in the Initial_conditions folder are as such: 

diskd_1000.0_gcm2_t50000_10.0_mIO_0.8_aIO_dissipFactor_1.0_rots_3.0.json

For that specific example: 

the disk density $\Sigma_0$ is 1000 g/cm$^2$

the disk lifetime is 50'000 yr

the mass of the satellite is 10 $M_\mathrm{Io}$

the initial satellite distance is 0.8 $a_\mathrm{Io}$

the dissipation factor (time lag) $\Delta t = 1.0 \Delta t_\oplus$ 

the initial rotation of the star is 3 days

References: [1] Leconte, J. and Chabrier, G.~2013, Layered convection as the origin of Saturn's luminosity anomaly, Nature Geoscience, 6, 347, doi: 10.1038/ngeo1791 (https://ui.adsabs.harvard.edu/abs/2013NatGe...6..347L)