A modular set of synthetic spectral energy distributions for young stellar objects

In this paper, I present a new set of synthetic spectral energy distributions (SEDs) for young stellar objects (YSOs) spanning a wide range of evolutionary stages, from the youngest deeply embedded protostars to pre-main-sequence stars with few or no disks. These models include significant improvements on the previous generation of published models: in particular, the new models cover a much wider and more uniform region of parameter space, do not include highly model-dependent parameters, and include a number of improvements that make them more suited to modeling far-infrared and sub-mm observations of forming stars. Rather than all being part of a single monolithic set of models, the new models are split up into sets of varying complexity. The aim of the new set of models is not to provide the most physically realistic models for young stars, but rather to provide deliberately simplified models for initial modeling, which allows a wide range of parameter space to be explored. I present the design of the model set, and show examples of fitting these models to real observations to show how the new grid design can help us better understand what can be determined from limited unresolved observations. The models, as well as a Python-based fitting tool are publicly available to the community.Comment: Accepted for publication in Astronomy and Astrophysics. The models are available at http://doi.org/10.5281/zenodo.16673

Interpreting Spectral Energy Distributions from Young Stellar Objects. I. A grid of 200,000 YSO model SEDs

We present a grid of radiation transfer models of axisymmetric young stellar objects (YSOs), covering a wide range of stellar masses (from 0.1Msun to 50Msun) and evolutionary stages (from the early envelope infall stage to the late disk-only stage). The grid consists of 20,000 YSO models, with spectral energy distributions (SEDs) and polarization spectra computed at ten viewing angles for each model, resulting in a total of 200,000 SEDs. [...]. These models are publicly available on a dedicated WWW server: http://www.astro.wisc.edu/protostars/ . In this paper we summarize the main features of our models, as well as the range of parameters explored. [...]. We examine the dependence of the spectral indices of the model SEDs on envelope accretion rate and disk mass. In addition, we show variations of spectral indices with stellar temperature, disk inner radius, and disk flaring power for a subset of disk-only models. We also examine how changing the wavelength range of data used to calculate spectral indices affects their values. We show sample color-color plots of the entire grid as well as simulated clusters at various distances with typical {\it Spitzer Space Telescope} sensitivities. We find that young embedded sources generally occupy a large region of color-color space due to inclination and stellar temperature effects. Disk sources occupy a smaller region of color-color space, but overlap substantially with the region occupied by embedded sources, especially in the near- and mid-IR. We identify regions in color-color space where our models indicate that only sources at a given evolutionary stage should lie. [...].Comment: 69 pages, 28 figures, Accepted for publication in ApJS. Preprint with full resolution figures available at http://www.astro.wisc.edu/protostars

The 3-Dimensional Structure of NGC 891 and M51

We investigate the three-dimensional structure of the nearby edge-on spiral galaxy NGC 891 using 3D Monte Carlo radiative transfer models, with realistic spiral structure and fractally clumped dust. Using the spiral and clumpiness parameters found from recently completed scattered light models we produce lower resolution SED models which reproduce the global UV-to-FIR SED of NGC 891. Our models contain a color gradient across the major axis of the galaxy - similar to what is seen in images of the NGC 891. With minor adjustment our SED models are able to match the majority of M51's SED, a similar galaxy at a near face-on different inclination.Comment: 3 pages, 1 figure. To appear in the Conference Proceedings of IAU Symposium No. 284: The Spectral Energy Distribution of Galaxies, R.J. Tuffs & C.C. Popescu, ed

Duration of Star Formation in Galactic Giant Molecular Clouds. I. The Great Nebula in Carina

We present a novel infrared spectral energy distribution (SED) modeling methodology that uses likelihood-based weighting of the model fitting results to construct probabilistic Hertzsprung–Russell diagrams (pHRD) for X-ray-identified, intermediate-mass (2–8 M⊙), pre-main-sequence young stellar populations. This methodology is designed specifically for application to young stellar populations suffering strong, differential extinction (ΔA_V > 10 mag), typical of Galactic massive star-forming regions. We pilot this technique in the Carina Nebula Complex (CNC) by modeling the 1–8 μm SEDs of 2269 likely stellar members that exhibit no excess emission from circumstellar dust disks at 4.5 μm or shorter wavelengths. A subset of ~100 intermediate-mass stars in the lightly obscured Trumpler 14 and 16 clusters have available spectroscopic T_(eff), measured from the Gaia-ESO survey. We correctly identify the stellar temperature in 85% of cases, and the aggregate pHRD for all sources returns the same peak in the stellar age distribution as obtained using the spectroscopic T_(eff). The SED model parameter distributions of stellar mass and evolutionary age reveal significant variation in the duration of star formation among four large-scale stellar overdensities within the CNC and a large distributed stellar population. Star formation began ~10 Myr ago and continues to the present day, with the star formation rate peaking ≾3 Myr ago when the massive Trumpler 14 and 16 clusters formed. We make public the set of 100,000 SED models generated from standard pre-main-sequence evolutionary tracks and our custom software package for generating pHRDs and mass–age distributions from the SED fitting results