Thomas P. Robitaille

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 circumstellar A _V. 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.

The NDCube 2 API is a Python application programming interface (API) for storing and manipulating N-dimensional coordinate-aware astronomical data. While there are Python packages for handling astronomical data and coordinate transformations separately and for handling specific combinations of dimensions and transformations, none provide a unified and agnostic way of handling them simultaneously. This leads to a proliferation of different APIs for conducting the same analysis tasks on similar types of observations and introduces technical barriers between multi-instrument studies and cross-community collaboration. In this paper, we outline how the NDCube 2 API and its implementation in the open-source, community-developed ndcube package, together with the AstroPy WCS API, help to solve this problem. We discuss the guiding principles underpinning the API design and provide examples of how it is already being used to serve broad sections of the astronomy community, including agency-funded missions. The aim of this paper is to help users better understand the purpose and potential of the NDCube 2 API and ndcube package and hence how to more effectively deploy them in scientific analyses and software development.

The Astropy Project supports and fosters the development of open-source and openly developed Python packages that provide commonly needed functionality to the astronomical community. A key element of the Astropy Project is the core package astropy, which serves as the foundation for more specialized projects and packages. In this article, we summarize key features in the core package as of the recent major release, version 5.0, and provide major updates on the Project. We then discuss supporting a broader ecosystem of interoperable packages, including connections with several astronomical observatories and missions. We also revisit the future outlook of the Astropy Project and the current status of Learn Astropy. We conclude by raising and discussing the current and future challenges facing the Project.

We present powderday (available at https://github.com/dnarayanan/powderday), a flexible, fast, open-source dust radiative transfer package designed to interface with both idealized and cosmological galaxy formation simulations. powderday builds on fsps stellar population synthesis models, and hyperion dust radiative transfer, and employs yt to interface between different software packages. We include our stellar population synthesis modeling on the fly, allowing significant flexibility in the assumed stellar physics and nebular line emission. The dust content follows either simple observationally motivated prescriptions (i.e., constant dust-to-metals ratios, or dust-to-gas ratios that vary with metallicity), direct modeling from galaxy formation simulations that include dust physics, as well as a novel approach that includes the dust content via learning-based algorithms from the simba cosmological galaxy formation simulation. Active galactic nuclei (AGNs) can additionally be included via a range of prescriptions. The output of these models are broadband (912 Å–1 mm) spectral energy distributions (SEDs), as well as filter-convolved monochromatic images. powderday is designed to eliminate last-mile efforts by researchers that employ different hydrodynamic galaxy formation models and seamlessly interfaces with gizmo, arepo, gasoline, changa, and enzo. We demonstrate the capabilities of the code via three applications: a model for the star formation rate–infrared luminosity relation in galaxies (including the impact of AGNs), the impact of circumstellar dust around AGB stars on the mid-infrared emission from galaxy SEDs, and the impact of galaxy inclination angle on dust attenuation laws.

We present a detailed analysis of the disk and circumstellar environment of the forming O-type star AFGL 4176 mm1, placing results from the Atacama Large Millimeter/submillimeter Array (ALMA) into context with multiwavelength data. With ALMA, we detect seventeen 1.2 mm continuum sources within 5″ (21,000 au) of AFGL 4176 mm1. We find that mm1 has a spectral index of 3.4 ± 0.2 across the ALMA band, with >87% of its 1.2 mm continuum emission from dust. The source mm2, projected 4200 au from mm1, may be a companion or a blueshifted knot in a jet. We also explore the morphological differences between the molecular lines detected with ALMA, finding 203 lines from 25 molecules, which we categorize into several morphological types. Our results show that AFGL 4176 mm1 provides an example of a forming O-star with a large and chemically complex disk, which is mainly traced by nitrogen-bearing molecules. Lines that show strong emission on the blueshifted side of the disk are predominantly oxygen-bearing, which we suggest are tracing a disk accretion shock. The molecules C³⁴S, H₂CS, and CH₃CCN trace a slow wide-angle wind or dense structures in the outflow cavity walls. With the Australia Telescope Compact Array (ATCA), we detect a compact continuum source (96% is from ionized gas. The ATCA NH₃(1, 1) and (2, 2) emission traces a large-scale (r ∼ 0.5 pc) rotating toroid with the disk source mm1 in the blueshifted part of this structure offset to the northwest.

We present high-resolution (30 mas or 130 au at 4.2 kpc) Atacama Large Millimeter/submillimeter Array observations at 1.2 mm of the disc around the forming O-type star AFGL 4176 mm1. The disc (AFGL 4176 mm1-main) has a radius of ∼1000 au and contains significant structure, most notably a spiral arm on its redshifted side. We fitted the observed spiral with logarithmic and Archimedean spiral models. We find that both models can describe its structure, but the Archimedean spiral with a varying pitch angle fits its morphology marginally better. As well as signatures of rotation across the disc, we observe gas arcs in CH₃CN that connect to other millimetre continuum sources in the field, supporting the picture of interactions within a small cluster around AFGL 4176 mm1-main. Using local thermodynamic equilibrium modelling of the CH₃CN K-ladder, we determine the temperature and velocity field across the disc, and thus produce a map of the Toomre stability parameter. Our results indicate that the outer disc is gravitationally unstable and has already fragmented or is likely to fragment in the future, possibly producing further companions. These observations provide evidence that disc fragmentation is one possible pathway towards explaining the high fraction of multiple systems around high-mass stars.

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.

We study a very young star-forming region in the outer Galaxy that is the most concentrated source of outflows in the Spitzer Space Telescope GLIMPSE360 survey. This region, dubbed CMa–l224, is located in the Canis Major OB1 association. CMa–l224 is relatively faint in the mid-infrared, but it shines brightly at the far-infrared wavelengths as revealed by the Herschel Space Observatory data from the Hi-GAL survey. Using the 3.6 and 4.5 μm data from the Spitzer/GLIMPSE360 survey, combined with the JHK _s Two Micron All Sky Survey (2MASS) and the 70–500 μm Herschel/Hi-GAL data, we develop young stellar object (YSO) selection criteria based on color–color cuts and fitting of the YSO candidates’ spectral energy distributions with YSO 2D radiative transfer models. We identify 293 YSO candidates and estimate physical parameters for 210 sources well fit with YSO models. We select an additional 47 sources with GLIMPSE360-only photometry as “possible YSO candidates.” The vast majority of these sources are associated with high H₂ column density regions and are good targets for follow-up studies. The distribution of YSO candidates at different evolutionary stages with respect to Herschel filaments supports the idea that stars are formed in the filaments and become more dispersed with time. Both the supernova-induced and spontaneous star formation scenarios are plausible in the environmental context of CMa–l224. However, our results indicate that a spontaneous gravitational collapse of filaments is a more likely scenario. The methods developed for CMa–l224 can be used for larger regions in the Galactic plane where the same set of photometry is available.