Source author record

Thomas P. Robitaille

Thomas P. Robitaille appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

24works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

24 published item(s)

preprint2022arXiv

The Astropy Project: Sustaining and Growing a Community-oriented Open-source Project and the Latest Major Release (v5.0) of the Core Package

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 $\texttt{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 for 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.

preprint2016arXiv

Rapid Circumstellar Disk Evolution and an Accelerating Star Formation Rate in the Infrared Dark Cloud M17 SWex

We present a catalog of 840 X-ray sources and first results from a 100 ks Chandra X-ray Observatory imaging study of the filamentary infrared dark cloud G014.225$-$00.506, which forms the central regions of a larger cloud complex known as the M17 southwest extension (M17 SWex). In addition to the rich population of protostars and young stellar objects with dusty circumstellar disks revealed by Spitzer Space Telescope archival data, we discover a population of X-ray-emitting, intermediate-mass pre--main-sequence stars (IMPS) that lack infrared excess emission from circumstellar disks. We model the infrared spectral energy distributions of this source population to measure its mass function and place new constraints on the inner dust disk destruction timescales for 2-8 $M_{\odot}$ stars. We also place a lower limit on the star formation rate (SFR) and find that it is quite high ($\dot{M}\ge 0.007~M_{\odot}$ yr$^{-1}$), equivalent to several Orion Nebula Clusters in G14.225$-$0.506 alone, and likely accelerating. The cloud complex has not produced a population of massive, O-type stars commensurate with its SFR. This absence of very massive (${\ge}20~M_{\odot}$) stars suggests that either (1) M17 SWex is an example of a distributed mode of star formation that will produce a large OB association dominated by intermediate-mass stars but relatively few massive clusters, or (2) the massive cores are still in the process of accreting sufficient mass to form massive clusters hosting O stars.

preprint2015arXiv

A Keplerian-like disk around the forming O-type star AFGL 4176

We present Atacama Large Millimeter/submillimeter Array (ALMA) line and continuum observations at 1.2mm with ~0.3" resolution that uncover a Keplerian-like disk around the forming O-type star AFGL 4176. The continuum emission from the disk at 1.21 mm (source mm1) has a deconvolved size of 870+/-110 AU x 330+/-300 AU and arises from a structure ~8 M_sun in mass, calculated assuming a dust temperature of 190 K. The first-moment maps, pixel-to-pixel line modeling, assuming local thermodynamic equilibrium (LTE), and position-velocity diagrams of the CH3CN J=13-12 K-line emission all show a velocity gradient along the major axis of the source, coupled with an increase in velocity at small radii, consistent with Keplerian-like rotation. The LTE line modeling shows that where CH3CN J=13-12 is excited, the temperatures in the disk range from ~70 to at least 300 K and that the H2 column density peaks at 2.8x10^24 cm^-2. In addition, we present Atacama Pathfinder Experiment (APEX) 12CO observations which show a large-scale outflow from AFGL 4176 perpendicular to the major axis of mm1, supporting the disk interpretation. Finally, we present a radiative transfer model of a Keplerian disk surrounding an O7 star, with a disk mass and radius of 12 M_sun and 2000 AU, that reproduces the line and continuum data, further supporting our conclusion that our observations have uncovered a Keplerian disk around an O-type star.

preprint2014arXiv

Main-sequence stars masquerading as Young Stellar Objects in the central molecular zone

In contrast to most other galaxies, star-formation rates in the Milky Way can be estimated directly from Young Stellar Objects (YSOs). In the Central Molecular Zone (CMZ) the star-formation rate calculated from the number of YSOs with 24 microns emission is up to order of magnitude higher than the value estimated from methods based on diffuse emission (such as free-free emission). Whether this effect is real or whether it indicates problems with either or both star formation rate measures is not currently known. In this paper, we investigate whether estimates based on YSOs could be heavily contaminated by more evolved objects such as main-sequence stars. We present radiative transfer models of YSOs and of main-sequence stars in a constant ambient medium which show that the main-sequence objects can indeed mimic YSOs at 24 microns. However, we show that in some cases the main-sequence models can be marginally resolved at 24 microns, whereas the YSO models are always unresolved. Based on the fraction of resolved MIPS 24 microns sources in the sample of YSOs previously used to compute the star formation rate, we estimate the fraction of misclassified "YSOs" to be at least 63%, which suggests that the star-formation rate previously determined from YSOs is likely to be at least a factor of three too high.

preprint2014arXiv

Spitzer View of Massive Star Formation in the Tidally Stripped Magellanic Bridge

The Magellanic Bridge is the nearest low-metallicity, tidally stripped environment, offering a unique high-resolution view of physical conditions in merging and forming galaxies. In this paper we present analysis of candidate massive young stellar objects (YSOs), i.e., {\it in situ, current} massive star formation (MSF) in the Bridge using {\it Spitzer} mid-IR and complementary optical and near-IR photometry. While we definitely find YSOs in the Bridge, the most massive are $\sim10 M_\odot$, $\ll45 M_\odot$ found in the Large Magellanic Cloud (LMC). The intensity of MSF in the Bridge also appears decreasing, as the most massive YSOs are less massive than those formed in the past. To investigate environmental effects on MSF, we have compared properties of massive YSOs in the Bridge to those in the LMC. First, YSOs in the Bridge are apparently less embedded than in the LMC: 81% of Bridge YSOs show optical counterparts, compared to only 56% of LMC sources with the same range of mass, circumstellar dust mass, and line-of-sight extinction. Circumstellar envelopes are evidently more porous or clumpy in the Bridge's low-metallicity environment. Second, we have used whole samples of YSOs in the LMC and the Bridge to estimate the probability of finding YSOs at a given \hi\ column density, N(HI). We found that the LMC has $\sim3\times$ higher probability than the Bridge for N(HI) $>10\times10^{20}$ cm$^{-2}$, but the trend reverses at lower N(HI). Investigating whether this lower efficiency relative to HI is due to less efficient molecular cloud formation, or less efficient cloud collapse, or both, will require sensitive molecular gas observations.

preprint2013arXiv

Astropy: A Community Python Package for Astronomy

We present the first public version (v0.2) of the open-source and community-developed Python package, Astropy. This package provides core astronomy-related functionality to the community, including support for domain-specific file formats such as Flexible Image Transport System (FITS) files, Virtual Observatory (VO) tables, and common ASCII table formats, unit and physical quantity conversions, physical constants specific to astronomy, celestial coordinate and time transformations, world coordinate system (WCS) support, generalized containers for representing gridded as well as tabular data, and a framework for cosmological transformations and conversions. Significant functionality is under active development, such as a model fitting framework, VO client and server tools, and aperture and point spread function (PSF) photometry tools. The core development team is actively making additions and enhancements to the current code base, and we encourage anyone interested to participate in the development of future Astropy versions.

preprint2012arXiv

A self-consistent model of Galactic stellar and dust infrared emission and the abundance of polycyclic aromatic hydrocarbons

We present a self-consistent three-dimensional Monte-Carlo radiative transfer model of the stellar and dust emission in the Milky-Way, and have computed synthetic observations of the 3.6 to 100 microns emission in the Galactic mid-plane. In order to compare the model to observations, we use the GLIMPSE, MIPSGAL, and IRAS surveys to construct total emission spectra, as well as longitude and latitude profiles for the emission. The distribution of stars and dust is taken from the SKY model, and the dust emissivities includes an approximation of the emission from polycyclic aromatic hydrocarbons in addition to thermal emission. The model emission is in broad agreement with the observations, but a few modifications are needed to obtain a good fit. Firstly, by adjusting the model to include two major and two minor spiral arms rather than four equal spiral arms, the fit to the longitude profiles for |l|>30 degrees can be improved. Secondly, introducing a deficit in the dust distribution in the inner Galaxy results in a better fit to the shape of the IRAS longitude profiles at 60 and 100 microns. With these modifications, the model fits the observed profiles well, although it systematically under-estimates the 5.8 and 8.0 microns fluxes. One way to resolve this discrepancy is to increase the abundance of PAH molecules by 50% compared to the original model, although we note that changes to the dust distribution or radiation field may provide alternative solutions. Finally, we use the model to quantify which stellar populations contribute the most to the heating of different dust types, and which stellar populations and dust types contribute the most to the emission at different wavelengths.

preprint2012arXiv

Observing Simulated Protostars with Outflows: How Accurate are Protostellar Properties Inferred from SEDs?

The properties of unresolved protostars and their local environment are frequently inferred from spectral energy distributions (SEDs) using radiative transfer modeling. We perform synthetic observations of realistic star formation simulations to evaluate the accuracy of properties inferred from fitting model SEDs to observations. We use ORION, an adaptive mesh refinement (AMR) three-dimensional gravito-radiation-hydrodynamics code, to simulate low-mass star formation in a turbulent molecular cloud including the effects of protostellar outflows. To obtain the dust temperature distribution and SEDs of the forming protostars, we post-process the simulations using HYPERION, a state-of-the-art Monte-Carlo radiative transfer code. We find that the ORION and HYPERION dust temperatures typically agree within a factor of two. We compare synthetic SEDs of embedded protostars for a range of evolutionary times, simulation resolutions, aperture sizes, and viewing angles. We demonstrate that complex, asymmetric gas morphology leads to a variety of classifications for individual objects as a function of viewing angle. We derive best-fit source parameters for each SED through comparison with a pre-computed grid of radiative transfer models. While the SED models correctly identify the evolutionary stage of the synthetic sources as embedded protostars, we show that the disk and stellar parameters can be very discrepant from the simulated values. Parameters such as the stellar accretion rate, stellar mass, and disk mass show better agreement, but can still deviate significantly, and the agreement may in some cases be artificially good due to the limited range of parameters in the set of model SEDs. Lack of correlation between the model and simulation properties in many individual instances cautions against over-interpreting properties inferred from SEDs for unresolved protostellar sources. (Abridged)

preprint2012arXiv

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.

preprint2012arXiv

The Milky Way Project: A statistical study of massive star formation associated with infrared bubbles

The Milky Way Project citizen science initiative recently increased the number of known infrared bubbles in the inner Galactic plane by an order of magnitude compared to previous studies. We present a detailed statistical analysis of this dataset with the Red MSX Source catalog of massive young stellar sources to investigate the association of these bubbles with massive star formation. We particularly address the question of massive triggered star formation near infrared bubbles. We find a strong positional correlation of massive young stellar objects (MYSOs) and H II regions with Milky Way Project bubbles at separations of < 2 bubble radii. As bubble sizes increase, a statistically significant overdensity of massive young sources emerges in the region of the bubble rims, possibly indicating the occurrence of triggered star formation. Based on numbers of bubble-associated RMS sources we find that 67+/-3% of MYSOs and (ultra)compact H II regions appear associated with a bubble. We estimate that approximately 22+/-2% of massive young stars may have formed as a result of feedback from expanding H II regions. Using MYSO-bubble correlations, we serendipitously recovered the location of the recently discovered massive cluster Mercer 81, suggesting the potential of such analyses for discovery of heavily extincted distant clusters.

preprint2011arXiv

A Mid-Infrared Census of Star Formation Activity in Bolocam Galactic Plane Survey Sources

We present the results of a search for mid-infrared signs of star formation activity in the 1.1 mm sources in the Bolocam Galactic Plane Survey (BGPS). We have correlated the BGPS catalog with available mid-IR Galactic plane catalogs based on the Spitzer Space Telescope GLIMPSE legacy survey and the Midcourse Space Experiment (MSX) Galactic plane survey. We find that 44% (3,712 of 8,358) of the BGPS sources contain at least one mid-IR source, including 2,457 of 5,067 (49%) within the area where all surveys overlap (10 deg < l < 65 deg). Accounting for chance alignments between the BGPS and mid-IR sources, we conservatively estimate that 20% of the BPGS sources within the area where all surveys overlap show signs of active star formation. We separate the BGPS sources into four groups based on their probability of star formation activity. Extended Green Objects (EGOs) and Red MSX Sources (RMS) make up the highest probability group, while the lowest probability group is comprised of "starless" BGPS sources which were not matched to any mid-IR sources. The mean 1.1 mm flux of each group increases with increasing probability of active star formation. We also find that the "starless" BGPS sources are the most compact, while the sources with the highest probability of star formation activity are on average more extended with large skirts of emission. A subsample of 280 BGPS sources with known distances demonstrates that mass and mean H_2 column density also increase with probability of star formation activity.

preprint2011arXiv

A Pan-Carina YSO Catalog: Intermediate-Mass Young Stellar Objects in the Carina Nebula Identified Via Mid-Infrared Excess Emission

We present a catalog of 1439 young stellar objects (YSOs) spanning the 1.42 deg^2 field surveyed by the Chandra Carina Complex Project (CCCP), which includes the major ionizing clusters and the most active sites of ongoing star formation within the Great Nebula in Carina. Candidate YSOs were identified via infrared (IR) excess emission from dusty circumstellar disks and envelopes, using data from the Spitzer Space Telescope Vela--Carina survey and the Two-Micron All Sky Survey. We model the 1--24 /mu m IR spectral energy distributions of the YSOs to constrain physical properties. Our Pan-Carina YSO Catalog (PCYC) is dominated by intermediate-mass (2 Msun < m < 10 Msun) objects with disks, including Herbig Ae/Be stars and their less evolved progenitors. The PCYC provides a valuable complementary dataset to the CCCP X-ray source catalogs, identifying 1029 YSOs in Carina with no X-ray detection. We also catalog 410 YSOs with X-ray counterparts, including 62 candidate protostars. Candidate protostars with X-ray detections tend to be more evolved than those without. In most cases, X-ray emission apparently originating from intermediate-mass, disk-dominated YSOs is consistent with the presence of low-mass companions, but we also find that X-ray emission correlates with cooler stellar photospheres and higher disk masses. We suggest that intermediate-mass YSOs produce X-rays during their early pre-main sequence evolution, perhaps driven by magnetic dynamo activity during the convective atmosphere phase, but this emission dies off as the stars approach the main sequence. Extrapolating over the stellar initial mass function scaled to the PCYC population, we predict a total population of >2x10^4 YSOs and a present-day star formation rate (SFR) of >0.008 Msun/yr. The global SFR in the Carina Nebula, averaged over the past ~5 Myr, has been approximately constant.

preprint2011arXiv

Candidate X-ray-Emitting OB Stars in the Carina Nebula Identified Via Infrared Spectral Energy Distributions

We report the results of a new survey of massive, OB stars throughout the Carina Nebula using the X-ray point source catalog provided by the Chandra Carina Complex Project (CCCP) in conjunction with infrared (IR) photometry from the Two Micron All-Sky Survey and the Spitzer Space Telescope Vela--Carina survey. Mid-IR photometry is relatively unaffected by extinction, hence it provides strong constraints on the luminosities of OB stars, assuming that their association with the Carina Nebula, and hence their distance, is confirmed. We fit model stellar atmospheres to the optical (UBV) and IR spectral energy distributions (SEDs) of 182 OB stars with known spectral types and measure the bolometric luminosity and extinction for each star. We find that the extinction law measured toward the OB stars has two components: Av=1--1.5 mag produced by foreground dust with a ratio of total-to-selective absorption Rv=3.1 plus a contribution from local dust with Rv>4.0 in the Carina molecular clouds that increases as Av increases. Using X-ray emission as a strong indicator of association with Carina, we identify 94 candidate OB stars with Lbol\geq10^4 Lsun by fitting their IR SEDs. If the candidate OB stars are eventually confirmed by follow-up spectroscopic observations, the number of cataloged OB stars in the Carina Nebula will increase by ~50%. Correcting for incompleteness due to OB stars falling below the Lbol cutoff or the CCCP detection limit, these results potentially double the size of the young massive stellar population.

preprint2011arXiv

HYPERION: An open-source parallelized three-dimensional dust continuum radiative transfer code

HYPERION is a new three-dimensional dust continuum Monte-Carlo radiative transfer code that is designed to be as generic as possible, allowing radiative transfer to be computed through a variety of three-dimensional grids. The main part of the code is problem-independent, and only requires an arbitrary three-dimensional density structure, dust properties, the position and properties of the illuminating sources, and parameters controlling the running and output of the code. HYPERION is parallelized, and is shown to scale well to thousands of processes. Two common benchmark models for protoplanetary disks were computed, and the results are found to be in excellent agreement with those from other codes. Finally, to demonstrate the capabilities of the code, dust temperatures, SEDs, and synthetic multi-wavelength images were computed for a dynamical simulation of a low-mass star formation region. HYPERION is being actively developed to include new features, and is publicly available (http://www.hyperion-rt.org).

preprint2011arXiv

Massive Young Stellar Objects in the Galactic Center. I. Spectroscopic Identification from Spitzer/IRS Observations

We present results from our spectroscopic study, using the Infrared Spectrograph (IRS) onboard the Spitzer Space Telescope, designed to identify massive young stellar objects (YSOs) in the Galactic Center (GC). Our sample of 107 YSO candidates was selected based on IRAC colors from the high spatial resolution, high sensitivity Spitzer/IRAC images in the Central Molecular Zone (CMZ), which spans the central ~300 pc region of the Milky Way Galaxy. We obtained IRS spectra over 5um to 35um using both high- and low-resolution IRS modules. We spectroscopically identify massive YSOs by the presence of a 15.4um shoulder on the absorption profile of 15um CO2 ice, suggestive of CO2 ice mixed with CH3OH ice on grains. This 15.4um shoulder is clearly observed in 16 sources and possibly observed in an additional 19 sources. We show that 9 massive YSOs also reveal molecular gas-phase absorption from CO2, C2H2, and/or HCN, which traces warm and dense gas in YSOs. Our results provide the first spectroscopic census of the massive YSO population in the GC. We fit YSO models to the observed spectral energy distributions and find YSO masses of 8 - 23 Msun, which generally agree with the masses derived from observed radio continuum emission. We find that about 50% of photometrically identified YSOs are confirmed with our spectroscopic study. This implies a preliminary star formation rate of ~0.07 Msun/yr at the GC.

preprint2011arXiv

The standard model of low-mass star formation applied to massive stars: multi-wavelength modelling of IRAS 20126+4104

In order to investigate whether massive stars form similarly to their low-mass counterparts, we have used the standard envelope plus disc geometry successfully applied to low-mass protostars to model the near-IR to sub-millimetre SED and several mid-IR images of the embedded massive star IRAS20126+4104. We have used a Monte Carlo radiative transfer dust code to model the continuum absorption, emission and scattering through two azimuthally symmetric dust geometries, the first consisting of a rotationally flattened envelope with outflow cavities, and the second which also includes a flared accretion disc. Our results show that the envelope plus disc model reproduces the observed SED and images more accurately than the model without a disc, although the latter model more closely reproduces the morphology of the mid-IR emission within a radius of 1.1" or ~1800au. We have put forward several possible causes of this discontinuity, including inner truncation of the disc due to stellar irradiation, or precession of the outflow cavity. Our best fitting envelope plus disc model has a disc radius of 9200 au. We find that it is unlikely that the outer regions of such a disc would be in hydrostatic or centrifugal equilibrium, however we calculate that the temperatures within the disc would keep it stable to fragmentation.

preprint2010arXiv

On the Modified Random Walk for Monte-Carlo Radiation Transfer

Min et al. (2009) presented two complementary techniques that use the diffusion approximation to allow efficient Monte-Carlo radiation transfer in very optically thick regions: a modified random walk and a partial diffusion approximation. In this note, I show that the calculations required for the modified random walk method can be significantly simplified. In particular, the diffusion coefficient and the mass absorption coefficients required for the modified random walk are in fact the same as the standard diffusion coefficient and the Planck mean mass absorption coefficient.

preprint2010arXiv

Spitzer Space Telescope observations of the Carina Nebula: The steady march of feedback-driven star formation

We report the first results of imaging the Carina Nebula with Spitzer/IRAC, providing a catalog of point sources and YSOs based on SED fits. We discuss several aspects of the extended emission, including dust pillars that result when a clumpy molecular cloud is shredded by massive star feedback. There are few "extended green objects" (EGOs) normally taken as signposts of outflow activity, and none of the HH jets detected optically are seen as EGOs. A population of "extended red objects" tends to be found around OB stars, some with clear bow-shocks. These are dusty shocks where stellar winds collide with flows off nearby clouds. Finally, the relative distributions of O stars and subclusters of YSOs as compared to dust pillars shows that while some YSOs are located within pillars, many more stars and YSOs reside just outside pillar heads. We suggest that pillars are transient phenomena, part of a continuous outwardly propagating wave of star formation driven by massive star feedback. As pillars are destroyed, they leave newly formed stars in their wake, which are then subsumed into the young OB association. Altogether, the current generation of YSOs shows no strong deviation from a normal IMF. The number of YSOs suggests a roughly constant star-formation rate over the past 3Myr, implying that star formation in pillars constitutes an important mechanism to construct unbound OB associations. Accelerated pillars may give birth to O-type stars that, after several Myr, could appear to have formed in isolation.

preprint2010arXiv

The present-day star formation rate of the Milky-Way determined from Spitzer detected young stellar objects

We present initial results from a population synthesis model aimed at determining the star formation rate of the Milky-Way. We find that a total star formation rate of 0.68 to 1.45 Msun/yr is able to reproduce the observed number of young stellar objects in the Spitzer/IRAC GLIMPSE survey of the Galactic plane, assuming simple prescriptions for the 3D Galactic distributions of YSOs and interstellar dust, and using model SEDs to predict the brightness and color of the synthetic YSOs at different wavelengths. This is the first Galaxy-wide measurement derived from pre-main-sequence objects themselves, rather than global observables such as the total radio continuum, Halpha, or FIR flux. The value obtained is slightly lower than, but generally consistent with previously determined values. We will extend this method in the future to fit the brightness, color, and angular distribution of YSOs, and simultaneously make use of multiple surveys, to place constraints on the input assumptions, and reduce uncertainties in the star formation rate estimate. Ultimately, this will be one of the most accurate methods for determining the Galactic star formation rate, as it makes use of stars of all masses (limited only by sensitivity) rather than solely massive stars or indirect tracers of massive stars.

preprint2009arXiv

First Spectroscopic Identification of Massive Young Stellar Objects in the Galactic Center

We report the detection of several molecular gas-phase and ice absorption features in three photometrically-selected young stellar object (YSO) candidates in the central 280 pc of the Milky Way. Our spectra, obtained with the Infrared Spectrograph (IRS) onboard the Spitzer Space Telescope, reveal gas-phase absorption from CO2 (15.0um), C2H2 (13.7um) and HCN (14.0um). We attribute this absorption to warm, dense gas in massive YSOs. We also detect strong and broad 15um CO2 ice absorption features, with a remarkable double-peaked structure. The prominent long-wavelength peak is due to CH3OH-rich ice grains, and is similar to those found in other known massive YSOs. Our IRS observations demonstrate the youth of these objects, and provide the first spectroscopic identification of massive YSOs in the Galactic Center.

preprint2009arXiv

IRS Characterization of a Debris Disk around an M-type star in NGC2547

We present 5 to 15 micron Spitzer Infrared Spectrograph (IRS) low resolution spectral data of a candidate debris disk around an M4.5 star identified as a likely member of the ~40 Myr old cluster NGC2547. The IRS spectrum shows a silicate emission feature, indicating the presence of warm, small, (sub)micron-sized dust grains in the disk. Of the fifteen previously known candidate debris disks around M-type stars, the one we discuss in this paper is the first to have an observed mid-infrared spectrum and is also the first to have measured silicate emission. We combined the IRS data with ancillary data (optical, JHKs, and Spitzer InfraRed Array Camera and 24 micron data) to build the spectral energy distribution (SED) of the source. Monte Carlo radiation transfer modeling of the SED characterized the dust disk as being very flat (h100=2AU) and extending inward within at least 0.13AU of the central star. Our analysis shows that the disk is collisionally dominated and is likely a debris disk.

preprint2009arXiv

The Extended Environment of M17: A Star Formation History

M17 is one of the youngest and most massive nearby star-formation regions in the Galaxy. It features a bright H II region erupting as a blister from the side of a giant molecular cloud (GMC). Combining photometry from the Spitzer GLIMPSE survey with complementary infrared (IR) surveys, we identify candidate young stellar objects (YSOs) throughout a 1.5 deg x 1 deg field that includes the M17 complex. The long sightline through the Galaxy behind M17 creates significant contamination in our YSO sample from unassociated sources with similar IR colors. Removing contaminants, we produce a highly-reliable catalog of 96 candidate YSOs with a high probability of association with the M17 complex. We fit model spectral energy distributions to these sources and constrain their physical properties. Extrapolating the mass function of 62 intermediate-mass YSOs (M >3 Msun), we estimate that >1000 stars are in the process of forming in the extended outer regions of M17. From IR survey images from IRAS and GLIMPSE, we find that M17 lies on the rim of a large shell structure ~0.5 deg in diameter (~20 pc at 2.1 kpc). We present new maps of CO and 13CO (J=2-1) emission, which show that the shell is a coherent, kinematic structure associated with M17 at v = 19 km/s. The shell is an extended bubble outlining the photodissociation region of a faint, diffuse H II region several Myr old. We provide evidence that massive star formation has been triggered by the expansion of the bubble. The formation of the massive cluster ionizing the M17 H II region itself may have been similarly triggered. We conclude that the star formation history in the extended environment of M17 has been punctuated by successive waves of massive star formation propagating through a GMC complex.

preprint2006arXiv

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. [...].

preprint2006arXiv

Interpreting Spectral Energy Distributions from Young Stellar Objects. II. Fitting observed SEDs using a large grid of pre-computed models

We present a method to analyze the spectral energy distributions (SEDs) of young stellar objects (YSOs). Our approach is to fit data with pre-computed 2-D radiation transfer models spanning a large region of parameter space. This allows us to determine not only a single set of physical parameter values but the entire range of values consistent with the multi-wavelength observations of a given source. In this way we hope to avoid any over-interpretation when modeling a set of data. We have constructed spectral energy distributions from optical to sub-mm wavelengths, including new Spitzer IRAC and MIPS photometry, for 30 young and spatially resolved sources in the Taurus-Auriga star-forming region. We demonstrate fitting model SEDs to these sources, and find that we correctly identify the evolutionary stage and physical parameters found from previous independent studies, such as disk mass, disk accretion rate, and stellar temperature. We also explore how fluxes at various wavelengths help to constrain physical parameters, and show examples of degeneracies that can occur when fitting SEDs. A web-based version of this tool is available to the community at http://www.astro.wisc.edu/protostars .