Long-Term Lithium Abundance Signatures following Planetary Engulfment [Dataset]

MESA inlist files and chemical network used for the published article Sevilla et. al (2022).

The inlist files given are for a 1 solar mass model.  To generate the other models in the paper with a different stellar mass, the initial mass in the pre-ZAMS inlist was changed.

We first used prezams.inlist to generate the model star, with mixing processes besides convection.  This model saved after this run would be used for all other model runs for the same stellar mass.

For a model run without planetary engulfment, we used one inlist to evolve the star through the main-sequence, and then obtained the desired data from the resulting profile and history files.

For model runs with planetary engulfment, we first ran the model using an accretion inlist to simulate the engulfment process by turning on accretion, and then we used a main-sequence (MS) inlist to evolve the star through the main-sequence to produce the desired data.

The corresponding inlists for each model run are as follows:

10 Earth mass engulfment with included all mixing processes considered in our paper, which were convection, overshoot, thermohaline mixing, elemental diffusion, and a min_D_mix coefficient: 10earth_accretion.inlist for accretion, and 10earth_postaccretion.inlist for MS.

10 Earth mass engulfment with all mixing processes except thermohaline mixing: 10earth_accretion_nothermohaline.inlist for accretion, and 10earth_postaccretion_nothermohaline.inlist for MS.

10 Earth mass engulfment with all mixing processes except elemental diffusion: 10earth_accretion_nodiffusion.inlist for accretion, and 10earth_postaccretion_nodiffusion.inlist for MS.

1 Earth mass engulfment with all mixing processes: 1earth_accretion.inlist for accretion, and 1earth_postaccretion.inlist for MS.

100 Earth mass engulfment with all mixing processes: 100earth_accretion.inlist for accretion, and 100earth_postaccretion.inlist for MS.

10 Earth mass engulfment with all mixing processes occurring 1 gyr after ZAMS: We first ran the model with evolve_1gyr.inlist to save a model of the star 1 gyr after ZAMS.  Then, we used 1gyr_accretion.inlist for accretion, and 1gyr_postaccretion.inlist for MS.

No planetary engulfment: No accretion inlist was used; we immediately ran no_accrete.inlist as the MS inlist.

We also uploaded the chemical network used for our model as the file accrete.net.