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

Emission from Very Small Grains and PAH Molecules in Monte Carlo Radiation Transfer Codes: Application to the Edge-On Disk of Gomez's Hamburger

We have modeled optical to far infrared images, photometry, and spectroscopy of the object known as Gomez's Hamburger. We reproduce the images and spectrum with an edge-on disk of mass 0.3M_sun and radius 1600AU, surrounding an A0 III star at a distance of 280pc. Our mass estimate is in excellent agreement with recent CO observations. However, our distance determination is more than an order of magnitude smaller than previous analyses which inaccurately interpreted the optical spectrum. To accurately model the infrared spectrum we have extended our Monte Carlo radiation transfer codes to include emission from polycyclic aromatic hydrocarbon (PAH) molecules and very small grains (VSG). We do this using pre-computed PAH/VSG emissivity files for a wide range of values of the mean intensity of the exciting radiation field. When Monte Carlo energy packets are absorbed by PAHs/VSGs we reprocess them to other wavelengths by sampling from the emissivity files, thus simulating the absorption and re-emission process without reproducing lengthy computations of statistical equilibrium, excitation and de-excitation in the complex many level molecules. Using emissivity lookup tables in our Monte Carlo codes gives the flexibility to use the latest grain physics calculations of PAH/VSG emissivity and opacity that are being continually updated in the light of higher resolution infrared spectra. We find our approach gives a good representation of the observed PAH spectrum from the disk of Gomez's Hamburger. Our models also indicate the PAHs/VSGs in the disk have a larger scaleheight than larger radiative equilibrium grains, providing evidence for dust coagulation and settling to the midplane.Comment: ApJ accepte

Calibration of evolutionary diagnostics in high-mass star formation

The evolutionary classification of massive clumps that are candidate progenitors of high-mass young stars and clusters relies on a variety of independent diagnostics based on observables from the near-infrared to the radio. A promising evolutionary indicator for massive and dense cluster-progenitor clumps is the L/M ratio between the bolometric luminosity and the mass of the clumps. With the aim of providing a quantitative calibration for this indicator we used SEPIA/APEX to obtain CH3C2H(12-11) observations, that is an excellent thermometer molecule probing densities > 10^5 cm^-3 , toward 51 dense clumps with M>1000 solar masses, and uniformly spanning -2 < Log(L/M) < 2.3. We identify three distinct ranges of L/M that can be associated to three distinct phases of star formation in massive clumps. For L/M <1 no clump is detected in CH3C2H , suggesting an inner envelope temperature below 30K. For 1< L/M < 10 we detect 58% of the clumps, with a temperature between 30 and 35 K independently from the exact value of L/M; such clumps are building up luminosity due to the formation of stars, but no star is yet able to significantly heat the inner clump regions. For L/M> 10 we detect all the clumps, with a gas temperature rising with Log(L/M), marking the appearance of a qualitatively different heating source within the clumps; such values are found towards clumps with UCHII counterparts, suggesting that the quantitative difference in T - L/M behaviour above L/M >10 is due to the first appearance of ZAMS stars in the clumps.Comment: Astrophysical Journal Letters, Accepte

Far-infrared signatures and inner hole sizes of protoplanetary discs undergoing inside-out dust dispersal

By means of radiative transfer simulation we study the evolution of the far-infrared colours of protoplanetary discs undergoing inside-out dispersal, often referred to as transition discs. We show that a brightening of the mid and far-infrared emission from these objects is a natural consequence of the removal of the inner disc. Our results can fully explain recent observations of transition discs in the Chamaleon and Lupus star forming regions from the Herschel Gould Belt Survey, which show a higher median for the 70?um (Herschel PACS 1) band of known transition objects compared with primordial discs. Our theoretical results hence support the suggestion that the 70?um band may be a powerful diagnostic for the identification of transition discs from photometry data, provided that the inner hole is larger than tens of AU, depending on spectral type. Furthermore we show that a comparison of photometry in the K , 12?um and 7u0?m bands to model tracks can provide a rough, but quick estimate of the inner hole size of these objects, provided their inclination is below ?85 degrees and the inner hole size is again larger than tens of AU.Comment: 8 pages, figure 4, accepted for MNRA

2-D and 3-D Radiation Transfer Models of High-Mass Star Formation

2-D and 3-D radiation transfer models of forming stars generally produce bluer 1-10 micron colors than 1-D models of the same evolutionary state and envelope mass. Therefore, 1-D models of the shortwave radiation will generally estimate a lower envelope mass and later evolutionary state than multidimensional models. 1-D models are probably reasonable for very young sources, or longwave analysis (wavelengths > 100 microns). In our 3-D models of high-mass stars in clumpy molecular clouds, we find no correlation between the depth of the 10 micron silicate feature and the longwave (> 100 micron) SED (which sets the envelope mass), even when the average optical extinction of the envelope is >100 magnitudes. This is in agreement with the observations of Faison et al. (1998) of several UltraCompact HII (UCHII) regions, suggesting that many of these sources are more evolved than embedded protostars. We have calculated a large grid of 2-D models and find substantial overlap between different evolutionary states in the mid-IR color-color diagrams. We have developed a model fitter to work in conjunction with the grid to analyze large datasets. This grid and fitter will be expanded and tested in 2005 and released to the public in 2006.Comment: 10 pages, 8 figures, to appear in the proceedings of IAU Symp 227, Massive Star Birth: A Crossroads of Astrophysics, (Cesaroni R., Churchwell E., Felli M., Walmsley C. editors