This directory contains supplementary materials for the following publication:
Robitaille, 2017, A modular set of synthetic spectral energy distributions for young stellar objects
The models themselves are available on Zenodo in a form usable by the SED fitting code. You can also find the raw Hyperion input/output files here.
To get started, take a look at this Jupyter notebook to find out more about how to use the models.
The packages required to run the scripts in this directory are:
- Python 2.7 or 3.3 or later (http://www.python.org)
- Numpy 1.10 or later (http://www.numpy.org)
- Matplotlib 1.4 or later (http://matplotlib.sourceforge.net)
- Astropy 1.0 or later (http://www.astropy.org)
- Jupyter notebook (https://jupyter-notebook.readthedocs.io)
- sedfitter 1.0 or later (http://sedfitter.readthedocs.io/en/stable/)
For any questions or issues with using the models, please open an issue in this repository. For issues with using the sedfitter Python package, open an issue in the repository for that package.
What is the value of r_0 in the paper?
There was an oversight in the paper which meant that the value of r_0 was not given. It is set to 1000 au.
How can I find out what the dust sublimation radius is for each model?
There is currently no table of r_sub values for each model, but you can find this out by downloading the Hyperion input/output files for the models you are interested in from here and then use code along the lines of:
>>> from hyperion.model import ModelOutput
>>> from hyperion.util.constants import au
>>> mo = ModelOutput('grids-1.1/s---s-i/output/a3/A3kQmQtj.rtout')
>>> grid = mo.get_quantities()
>>> grid.r[1] / au
0.0023819780734823503
to find out the r_sub radius in au for each model.
What is mu in Equations (4) and (6) of paper?
The mu variable is defined as cos(theta).
This repository will evolve over time. In particular, I may:
- Add the scripts to reproduce the results from the paper
- Add the scripts used to construct, run, and post-process the model grid
To get announcements when these are made available, you can sign up to the following mailing list:
https://groups.google.com/forum/#!forum/protostars
Note that you do not need a Google account, and can sign up by sending an empty email to:
protostars+subscribe@googlegroups.com
For any questions, please open an issue in this repository.