Effects of stellar density on the photoevaporation of circumstellar discs [2021MNRAS.501.1782C]

Francisca Concha-Ramírez, Martijn J. C. Wilhelm, Simon Portegies Zwart, Sierk E. van Terwisga, Alvaro Hacar

(Code to run the simulations and generate figures of the paper: Effects of stellar density on the photoevaporation of circumstellar discs.)

For questions please contact Francisca Concha-Ramírez, fconcha at strw.leidenuniv.nl

Getting Started

Requisites

• Python 2.7. Should work fine with Python 3 but it has not been tested.
• AMUSE: https://amusecode.github.io
• VADER: https://bitbucket.org/krumholz/vader/src
• scipy

How to set up

You need to install AMUSE with the developer option so that you can access the source code. Then, download VADER and put it inside the folder: /amuse/src/amuse/community. In the /vader folder in this repository there are several files related to VADER. They should go in the following directories:

• interface.cc and interface.py should go on /amuse/src/amuse/community/vader/
• userFunc_pedisk.c and userFunc_none.c should go on amuse/src/amuse/community/vader/src/prob (these 2 files are redundant but it was a way to go around some compilation issues)
• Makefile_interface should go on /amuse/src/amuse/community/vader/ and renamed Makefile
• Makefile_source should go on /amuse/src/amuse/community/vader/src/ and renamed Makefile
• Compile VADER from the main AMUSE directory with make vader.code
• Add this line to amuse/src/amuse/lab.py:

from amuse.community.vader.interface import vader, vaderInterface

You should now be ready to run the simulation scripts.

Running the simulations

The main file to run the simulations is vader_cluster_parallel.py located in the src folder. You can run an individual simulation by using the AMUSE script directly from the home directory:

amuse.sh src/vader_cluster_parallel.py

The script has extensive options which can be passed through the command line, for number of stars, radius of the cluster, etc. For a list of these options run:

amuse.sh src/vader_cluster_parallel.py --help

Creating the plots

There are several different scripts to create the Figures of the paper. All of them are located in the plots directoy:

• Figure 1: tracks.py
• Figure 2: disc_fractions.py
• Figures 3 and 4: binned.py

Each script can be run with

amuse.sh plots/<scriptname>.py

A list of options is available for each script, including the path to the files that you want to use for the plots. To see the list of options add --help or -h to the line above. The results used on the paper are available in the results/ directory.

License

This project is licensed under the GPL v3 License - see the LICENSE.md file for details