PySSED: an automated method of collating and fitting stellar spectral energy distributions

Iain McDonald, Albert A. Zijlstra, Nick L. J. Cox, Emma L. Alexander, Alexander Csukai, Ria Ramkumar, Alexander Hollings
2024-02-20
Abstract:Stellar atmosphere modelling predicts the luminosity and temperature of a star, together with parameters such as the effective gravity and the metallicity, by reproducing the observed spectral energy distribution. Most observational data comes from photometric surveys, using a variety of passbands. We herein present the Python Stellar Spectral Energy Distribution (PySSED) routine, designed to combine photometry from disparate catalogues, fit the luminosity and temperature of stars, and determine departures from stellar atmosphere models such as infrared or ultraviolet excess. We detail the routine's operation, and present use cases on both individual stars, stellar populations, and wider regions of the sky. PySSED benefits from fully automated processing, allowing fitting of arbitrarily large datasets at the rate of a few seconds per star.
Instrumentation and Methods for Astrophysics
What problem does this paper attempt to address?
### Problems the Paper Aims to Solve This paper introduces a new method called PySSED (Python Stellar Spectral Energy Distribution), which aims to automate the collection and fitting of stellar spectral energy distributions (SEDs). The main objectives are as follows: 1. **Integrate Multi-band Data**: Automatically collect photometric data from different catalogs and combine them into a complete SED. 2. **Stellar Parameter Fitting**: Quickly and accurately fit stellar temperature and luminosity, and identify deviations from standard stellar atmosphere models, such as infrared or ultraviolet excess. 3. **Automated Processing**: Achieve a fully automated processing workflow, enabling the handling of datasets of any size at a speed of several seconds per star. ### Research Background and Motivation Traditional stellar atmosphere modeling typically relies on spectral observations to calculate parameters such as stellar temperature, surface gravity, and metallicity. However, this method has some limitations: - Requires high-quality spectral data, which is difficult to obtain for fainter stars. - The precision of surface gravity is usually not high. - Computationally expensive. In contrast, by analyzing stellar spectral energy distributions (SEDs), multi-band observational data can be used to better fit the basic parameters of stars. This not only improves data utilization but also better identifies objects that do not conform to simple stellar models, such as stars with binary companions, gas disks, or dust winds. ### Main Contributions The main features of PySSED include: - Automated Processing: Acquiring data from multiple databases and performing photometric matching and correction. - Outlier Exclusion: Automatically identifying and excluding bad data points during the fitting process. - Distance Estimation: Prioritizing the use of probabilistically corrected distance estimates to reduce Lutz-Kelker bias. - Extinction Correction: Using 3D extinction maps provided by Gaia data for correction. These features enable PySSED to efficiently handle large-scale astronomical data, whether for individual stars or stellar populations.