Acknowledgements: MUSCLES: If you make use of the MUSCLES spectral data products, please cite Papers I, II, and III from the top of this page and specify the version number of the spectra you used in your work. The version number is in the filenames. For example, "v22" denotes version 2.2.
Mega-MUSCLES: If the spectral data products are from the Mega-MUSCLES survey (which have versions of "v23" and are denoted by a double asterisk ** in the Data Access table), please cite Mega-MUSCLE Paper II from the top of this page.
Synthetic Models:If you use the semi-empirical synthetic models, please cite the relevant paper(s) in the "Semi-Empirical Synthetic Models" section at the top of the page.
MUSCLES is a spectral survey of 11 low-mass, planet-hosting stars, 7 M and 4 K dwarfs. The follow-on Mega-MUSCLES survey observed 12 more M dwarfs, extending the spectral type coverage down to M8. The spectra cover wavelengths from 5 Å to 5.5 µm (MUSCLES) / 1000 µm (Mega-MUSCLES), with emphasis on high-energy radiation. More information on the MUSCLES survey and collaboration can be found at the MUSCLES team page. Data sources for the various regions of the spectra are:
The data products from all stars are organized into subdirectories based on the star name. Several files are available for each target. We provide a brief summary below, for full details, consult the README file.
Primary Data Products
*const-res-sed.fits = Panchromatic SED binned to a constant 1 Å resolution. One file per star.
*var-res-sed.fits = Panchromatic SED that retains the native instrument resolutions where possible. One file per star.
*adapt-const-res-sed.fits = Panchromatic SED binned to a constant 1 Å resolution, downsampled in low signal-to-noise regions to avoid negative fluxes. One file per star.
*adapt-var-res-sed.fits = Panchromatic SED that retains the native instrument resolutions where possible, downsampled in low signal-to-noise regions to avoid negative fluxes. One file per star.
Ancillary Data Products
*hst_cos*component-spec.fits = Coadded HST COS spectrum. One file per observed grating (three for each star, except V-EPS-ERI, which has two).
*hst_stis*component-spec.fits = Coadded HST STIS spectrum. One file per observed grating (three for each star, except V-EPS-ERI, which has four).
*euv-scaling*component-spec.fits = Extreme UV spectrum based on empirical scaling relation using Lyα flux. One file per star.
*lya-reconstruction*component-spec.fits = Reconstructed spectrum from model fits to Lyα line profiles accounting for the stellar flux and the interstellar medium. One file per star.
*phx*component-spec.fits = Synthetic spectrum in the visible and infrared from PHOENIX models. One file per star.
Proxima Centauri (GJ 551)
In response to popular demand after the announcement of the discovery of a habitable-zone planet orbiting Prox Cen, we added Prox Cen to the MUSCLES sample and created an SED from archival data using methodology consistent with the other MUSCLES SEDs. We highly encourage users of the Prox Cen spectra to read these notes on the reduction of the Prox Cen spectra, as they address an important issue regarding the stellar effective temperature and bolometric flux.
Also available for download is a synthetic spectrum of the MUSCLES M2 dwarf GJ 832 from Fontenla et al. 2016. The synthetic spectrum was created from a one-dimensional model of the stellar atmosphere that incorporates non-LTE radiative transfer techniques and many molecular lines. The synthetic spectrum covers 1 Å - 1 mm at a resolving power R = 100,000 (denoted by "r1e5" in the file names). A single resolution element equals two wavelength samples. For consistency with the MUSCLES data products, the synthetic spectrum's wavelength is in Angstroms and the flux is reported as would be observed from Earth (erg/cm2/s/Å) under the assumption that the stellar radius = 0.499 R_solar (Houdebine 2010) and distance = 4.95 pc (van Leeuwen 2007). Note that Houdebine et al. 2016 determine GJ 832's stellar radius to be 0.458 ± 0.039 R_solar. The atomic and molecular transitions appear at their laboratory rest wavelengths (vacuum). Semi-empirical models of the other 10 MUSCLES stars as well as Proxima Centauri and TRAPPIST-1 are coming soon. If you utilize the GJ 832 synthetic spectrum, please cite Fontenla et al. 2016, ApJ, 830, 154.
Click on a star name to visit the Simbad page for that target. Download a .tar.gz containing all FITS files for a target using the (tar) link. Click on "FITS" to directly download the FITS file. Click on "PL" to use the MAST interactive plotter.
Target's with a double asterisk (**) denote Mega-MUSCLES targets and should cite the Mega-MUSCLES papers when used in publications.
Note: The STIS observations used a few different modes for a given central wavelength, e.g., the _140_ column contains spectra observed with either the G140L, G140M, or E140M modes. The EUV column might include _euv-scaling_ and _dem_ spectra. In all cases, consult the FITS header and/or file name to get the full observation modes used for a given target.
NOTICE: Until 20 Feb. 2023, the STIS G140M files for GJ1132, GJ163, LHS-2686 contained incorrect spectra. If you downloaded these spectra before that date then please re-download them, as any analysis based on the older spectra will be incorrect. The rest of the files and SED are unaffected for these targets.
†NOTE: The X-ray spectrum for HD 97658 does not come from XMM or Chandra observations. It is scaled from X-ray data of HD 85512, which shows similar lines in the UV. For more details, consult the README or primary publication.
‡NOTE: The links in the top row of the "_230_" cell are for the E230M data, while the links in the bottom row are for the E230H data.
✠NOTE: The X-ray spectrum for L980-5 is based on the upper limit flux from a Chandra non-detection.
The example file used here, hlsp_muscles_multi_multi_gj436_broadband_v10_adapt-const-res-sed.fits, can be downloaded by clicking the "FITS" link in the GJ 436 / Adaptive Res row under the Panchromatic / Const column.
from astropy.io import fits
import matplotlib.pyplot as plt
spec = fits.getdata('hlsp_muscles_multi_multi_gj436_broadband_v10_adapt-const-res-sed.fits',1)
plt.ylabel('Flux Density (erg/cm2/s/Ang)')