The DESK is an SED-fitting python scripts for fitting data from evolved stars


Keywords
Dusty-Evolved-Star-Kit
License
BSD-3-Clause
Install
pip install desk==1.6.25

Documentation

Dusty-Evolved-Star-Kit

Build Documentation Status Codacy Badge codecov pypi Steve Goldman

The DESK is an SED-fitting python package for fitting data from evolved stars (photometry or spectra) with radiative transfer model grids. The package is currently in development and all contributions are welcomed. For current progress, see the Issues tab at the top of the page. The package is ideal for fitting small samples of dusty evolved stars. It will soon utilize a bayesian-fitting strategy with mass-loss rate and luminosity distributions as inputs (priors), and will provide a better fit to these broader sample properties.

Input: A csv file with the first column as wavelength in um and second column as flux in Jy. To fit multiple csv files, put them in a directory, and use the directory name as the input.

Output: Two results files including the best fit model and corresponding stellar parameters, as well as an optional figure of the fit SED.

Available model grids: Several grids are already available upon installation. Descriptions of the model grids can be found in the Documentation. You will soon be able to specify the state-of-the-art dust growth models by Nanni et al. (2019) and the 2D GRAMS model grid based on the 2DUST code, which are automatically downloaded and used when selected.

A module for creating your own DUSTY grid is under development, but for now, please email me (Dr. Steven Goldman) directly for grid requests or for help with the package.

Documentation

The documentation can be found on readthedocs.

Install Using Python

1). Install the package with the command pip install desk.

Using the DESK

2). Go to the directory where your target csv file (or target directory of files) is.

3). Use the command (without starting python)

desk fit --source='target_name.csv'

or if you have a folder of csv files

desk fit --source='folder_of_csvs'

additional options are:

desk fit --source='target_name.csv' --distance=50 --grid='Oss-Orich-bb'

These additional arguments are the distance (in kpc) and the grid of models you would like to use (options listed below). You can select 'oxygen' or 'carbon' to use the default models. To see other available grids use:

desk grids

To create a figure showing all of the fits of the SED, use the following command in the same directory.

desk sed

This is an example of the output_sed.png file fitting three massive oxygen-rich AGB stars from the LMC.

To produce individual figures for each SED instead use the command:

desk sed_indiv

Retrieve model

To retrieve a model from the DUSTY model grids or to interpolate a model in the grid parameter space, use the save_model module with the following options:

desk save_model grid_name=Oss-Orich-bb luminosity=10000 teff=2700 tinner=1000 tau=0.4 distance_in_kpc=50

Attribution

The method used is similar to that of Goldman et al. 2017; a more in-depth publication is in prep. If used please add the following to your acknowledgements:

This research has made use of the Dusty Evolved Star Kit (DESK; https://github.com/s-goldman/Dusty-Evolved-Star-Kit).

Please also specify the options selected and make the data publicly available for reproducibility.

License

This project is Copyright (c) Dr. Steven Goldman and licensed under the terms of the BSD 3-Clause license.