Leaky Dust Traps

This repository contains the DustPy setups and a notebook producing the figures of the publication Leaky dust traps: How fragmentation impacts dust filtering by planets.

Installation

The models contain Fortran modules that need to be compiled first. To do so clone the repository and run the setup script.
NumPy and a Fotran compiler is required.

git clone git@github.com:stammler/leaky_dust_traps.git
cd leaky_dust_traps
pip install .

Running Models

To run a model change to the respective directory and run the model script

cd models/toy_saturn/
python start.py

Please note, that some models will take several days to run.
For more details, please have a look at the DustPy documentation.

Models

The following model setups are available in the repository

Name	Planet mass	Description
toy_saturn	$M_\mathrm{sat}$	Toy model, Saturn mass planet
toy_noPlanet	–	Toy model, without planet
toy_saturn_noFrag	$M_\mathrm{sat}$	Toy model, Saturn mass planet, no fragmentation
toy_saturn_bouncing	$M_\mathrm{sat}$	Toy model, Saturn mass planet, bouncing
toy_saturn_vFrag100	$M_\mathrm{sat}$	Toy model, Saturn mass planet, $v_\mathrm{frag}=1$ m/s
full_noPlanet	–	Full model, without planet
full_30Me	$30M_\oplus$	Full model, 30 Earth mass planet
full_50Me	$50M_\oplus$	Full model, 50 Earth mass planet
full_saturn_deltar2	$M_\mathrm{sat}$	Full model, Saturn Earth mass planet, $\delta_r = 10^{-2}$
full_saturn	$M_\mathrm{sat}$	Full model, without planet
full_saturn_deltar4	$M_\mathrm{sat}$	Full model, Saturn Earth mass planet, $\delta_r = 10^{-4}$
full_saturn_deltar5	$M_\mathrm{sat}$	Full model, Saturn Earth mass planet, $\delta_r = 10^{-5}$
full_200Me	$200M_\oplus$	Full model, 200 Earth mass planet
full_jupiter	$M_\mathrm{jup}$	Full model, Jupiter mass planet
full_MtSlow	$M\left(t\right)$	Full model, slowly growing planet
full_MtSlow_deltar4	$M\left(t\right)$	Full model, slowly growing planet, $\delta_r = 10^{-4}$
full_MtSlow_deltar5	$M\left(t\right)$	Full model, slowly growing planet, $\delta_r = 10^{-5}$
full_MtSlow_bouncing	$M\left(t\right)$	Full model, slowly growing planet, bouncing
full_MtSlow_bouncing_deltas5	$M\left(t\right)$	Full model, slowly growing planet, bouncing, $\delta_i = 10^{-5}$
full_MtFast_bouncing_deltas5	$M\left(t\right)$	Full model, rapidly growing planet, bouncing, $\delta_i = 10^{-5}$

Notebook

The notebook producing the figures of the publication can be found here.

Citation

Please cite the following publication

@ARTICLE{2023A&A...670L...5S,
       author = {{Stammler}, Sebastian Markus and {Lichtenberg}, Tim and {Dr{\k{a}}{\.z}kowska}, Joanna and {Birnstiel}, Tilman},
        title = "{Leaky dust traps: How fragmentation impacts dust filtering by planets}",
      journal = {\aap},
     keywords = {meteorites, meteors, meteoroids, methods: numerical, protoplanetary disks, planets and satellites: formation, planets and satellites: composition, Astrophysics - Earth and Planetary Astrophysics},
         year = 2023,
        month = feb,
       volume = {670},
          eid = {L5},
        pages = {L5},
          doi = {10.1051/0004-6361/202245512},
archivePrefix = {arXiv},
       eprint = {2301.05505},
 primaryClass = {astro-ph.EP},
       adsurl = {https://ui.adsabs.harvard.edu/abs/2023A&A...670L...5S},
      adsnote = {Provided by the SAO/NASA Astrophysics Data System}
}

Acknowledgements

This project has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme under grant agreement No 714769. This project has received funding by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) through grants FOR 2634/1 and 361140270. This research was supported by the Munich Institute for Astro-, Particle and BioPhysics (MIAPbP) which is funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy – EXC-2094 – 390783311. JD was funded by the European Union under the European Union’s Horizon Europe Research & Innovation Programme 101040037 (PLANETOIDS). Views and opinions expressed are however those of the authors only and do not necessarily reflect those of the European Union or the European Research Council. Neither the European Union nor the granting authority can be held responsible for them. TL was supported by grants from the Simons Foundation (SCOL Award No. 611576) and the Branco Weiss Foundation.