skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: An Updated Modular Set of Synthetic Spectral Energy Distributions for Young Stellar Objects
Abstract Measured properties of young stellar objects (YSOs) are key tools for research into pre-main-sequence stellar evolution. YSO properties are commonly measured by comparing observed radiation to existing grids of template YSO spectral energy distributions (SEDs) calculated by radiative transfer. These grids are often sampled and constructed using simple models of mass assembly/accretion over time. However, because we do not yet have a complete theory of star formation, the choice of model sets the tracked parameters and range of allowed values. By construction, then, the assumed model limits the measurements that can be made using the grid. Radiative transfer models not constrained by specific accretion histories would enable assessment of a wider range of theories. We present an updated version of the Robitaille set (2017) of YSO SEDs, a collection of models with no assumed evolutionary theory. We outline our newly calculated properties: envelope mass, weighted-average dust temperature, disk stability, and circumstellarAV. We also convolve the SEDs with new filters, including JWST, and provide users the ability to perform additional convolutions. We find a correlation between the average temperature and millimeter-wavelength brightness of optically thin dust in our models and discuss its ramifications for mass measurements of pre- and protostellar cores. We also compare the positions of YSOs of different observational classes and evolutionary stages in IR color space and use our models to quantify the extent to which class and stage may be confused due to observational effects. Our updated models are released to the public.  more » « less
Award ID(s):
2008101 2142300
PAR ID:
10487597
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
961
Issue:
2
ISSN:
0004-637X
Format(s):
Medium: X Size: Article No. 188
Size(s):
Article No. 188
Sponsoring Org:
National Science Foundation
More Like this
  1. Context. The Central Molecular Zone (CMZ), a ∼200 pc sized region around the Galactic Centre, is peculiar in that it shows a star formation rate (SFR) that is suppressed with respect to the available dense gas. To study the SFR in the CMZ, young stellar objects (YSOs) can be investigated. Here we present radio observations of 334 2.2 μm infrared sources that have been identified as YSO candidates. Aims: Our goal is to investigate the presence of centimetre wavelength radio continuum counterparts to this sample of YSO candidates which we use to constrain the current SFR in the CMZ. Methods: As part of the GLObal view on STAR formation (GLOSTAR) survey, D-configuration Very Large Array data were obtained for the Galactic Centre, covering −2° < l < 2° and −1° < b < 1° with a frequency coverage of 4-8 GHz. We matched YSOs with radio continuum sources based on selection criteria and classified these radio sources as potential H II regions and determined their physical properties. Results: Of the 334 YSO candidates, we found 35 with radio continuum counterparts. We find that 94 YSOs are associated with dense dust condensations identified in the 870 μm ATLASGAL survey, of which 14 have a GLOSTAR counterpart. Of the 35 YSOs with radio counterparts, 11 are confirmed as H II regions based on their spectral indices and the literature. We estimated their Lyman continuum photon flux in order to estimate the mass of the ionising star. Combining these with known sources, the present-day SFR in the CMZ is calculated to be ∼0.068 M⊙ yr−1, which is ∼6.8% of the Galactic SFR. Candidate YSOs that lack radio counterparts may not have yet evolved to the stage of exhibiting an H II region or, conversely, are older and have dispersed their natal clouds. Since many lack dust emission, the latter is more likely. Our SFR estimate in the CMZ is in agreement with previous estimates in the literature. 
    more » « less
  2. ABSTRACT Forward-modeling observables from galaxy simulations enables direct comparisons between theory and observations. To generate synthetic spectral energy distributions (SEDs) that include dust absorption, re-emission, and scattering, Monte Carlo radiative transfer is often used in post-processing on a galaxy-by-galaxy basis. However, this is computationally expensive, especially if one wants to make predictions for suites of many cosmological simulations. To alleviate this computational burden, we have developed a radiative transfer emulator using an artificial neural network (ANN), ANNgelina, that can reliably predict SEDs of simulated galaxies using a small number of integrated properties of the simulated galaxies: star formation rate, stellar and dust masses, and mass-weighted metallicities of all star particles and of only star particles with age <10 Myr. Here, we present the methodology and quantify the accuracy of the predictions. We train the ANN on SEDs computed for galaxies from the IllustrisTNG project’s TNG50 cosmological magnetohydrodynamical simulation. ANNgelina is able to predict the SEDs of TNG50 galaxies in the ultraviolet (UV) to millimetre regime with a typical median absolute error of ∼7 per cent. The prediction error is the greatest in the UV, possibly due to the viewing-angle dependence being greatest in this wavelength regime. Our results demonstrate that our ANN-based emulator is a promising computationally inexpensive alternative for forward-modeling galaxy SEDs from cosmological simulations. 
    more » « less
  3. Abstract The potential for planet formation of a circumstellar disk depends on the dust and gas reservoirs, which evolve as a function of the disk age. The Atacama Large Millimeter/submillimeter Array AGE-PRO Large Program has measured several disk properties across three star-forming regions of different ages, and in this study, we compare the observational results to dust evolution simulations. UsingDustPyfor the dust evolution, andRADMC-3Dfor the radiative transfer, we ran a large grid of models spanning stellar masses of 0.25, 0.50, 0.75, and 1.0M, with different initial conditions, including: disk sizes, disk gas masses, and dust-to-gas ratio, and viscosity. Our models are performed assuming smooth, weakly, or strongly substructured disks, aiming to investigate if any observational trend can favor or exclude the presence of dust traps. The observed gas masses in the disks of the AGE-PRO sample are not reproducible with our models, which only consider viscous evolution with constantα, suggesting that additional physical mechanisms play a role in the evolution of the gas mass of disks. When comparing the dust continuum emission fluxes and sizes at 1.3 mm, we find that most of the disks in the AGE-PRO sample are consistent with simulations that have either weak or strong dust traps. The evolution of spectral index in the AGE-PRO sample is also suggestive of an unresolved population of dust traps. Future observations at high angular resolution are still needed to test several hypotheses that result from comparing the observations to our simulations, including that more massive disks in gas mass have the potential to form dust traps at larger disk radii. 
    more » « less
  4. Abstract We present ∼8–40μm SOFIA-FORCAST images of seven regions of “clustered” star formation as part of the SOFIA Massive Star Formation Survey. We identify a total of 34 protostar candidates and build their spectral energy distributions (SEDs). We fit these SEDs with a grid of radiative transfer models based on the turbulent core accretion (TCA) theory to derive key protostellar properties, including initial core mass,Mc, clump environment mass surface density, Σcl, and current protostellar mass,m*. We also carry out empirical graybody (GB) estimation of Σcl, which allows a case of restricted SED fitting within the TCA model grid. We also release version 2.0 of the open-source Python packagesedcreator, which is designed to automate the aperture photometry and SED building and fitting process for sources in clustered environments, where flux contamination from close neighbors typically complicates the process. Using these updated methods, SED fitting yields values ofMc∼ 30–200M, Σcl,SED∼ 0.1–3 g cm−2, andm*∼ 4–50M. The GB fitting yields smaller values of Σcl,GB≲ 1 g cm−2. From these results, we do not find evidence for a critical Σclneeded to form massive (≳8M) stars. However, we do find tentative evidence for a dearth of the most massive (m*≳ 30M) protostars in the clustered regions, suggesting a potential impact of environment on the stellar initial mass function. 
    more » « less
  5. null (Ed.)
    Context. FU Orionis is the archetypal FUor star, a subclass of young stellar objects (YSOs) that undergo rapid brightening events, often gaining between four and six magnitudes on timescales of days. This brightening is often associated with a massive increase in accretion, which is one of the most ubiquitous processes in astrophysics for bodies ranging from planets and stars to super-massive black holes. We present multi-band interferometric observations of the FU Ori circumstellar environment, including the first J -band interferometric observations of a YSO. Aims. We investigate the morphology and temperature gradient of the innermost regions of the accretion disk around FU Orionis. We aim to characterise the heating mechanisms of the disk and comment on potential outburst-triggering processes. Methods. Recent upgrades to the MIRC-X instrument at the CHARA array have allowed for the first dual-band J and H observations of YSOs. Using baselines up to 331 m, we present high-angular-resolution data of a YSO covering the near-infrared bands J , H , and K . The unprecedented spectral range of the data allowed us to apply temperature gradient models to the innermost regions of FU Ori. Results. We spatially resolved the innermost astronomical unit of the disk and determine the exponent of the temperature gradient of the inner disk to T ∝ r −0.74 ± 0.02 . This agrees with theoretical works that predict T ∝ r −0.75 for actively accreting, steady-state disks, which is a value only obtainable through viscous heating within the disk. We found a disk that extends down to the stellar surface at 0.015 ± 0.007 au, where the temperature is found to be 5800 ± 700 K. We found a disk inclined at 32 ± 4° with a minor-axis position angle of 34 ± 11°. Conclusions. We demonstrate that J -band interferometric observations of YSOs are feasible with the MIRC-X instrument at CHARA. The temperature gradient power-law derived for the inner disk is consistent with theoretical predictions for steady-state, optically thick, viciously heated accretion disks. 
    more » « less