Issue |
A&A
Volume 518, July-August 2010
Herschel: the first science highlights
|
|
---|---|---|
Article Number | L119 | |
Number of page(s) | 4 | |
Section | Letters | |
DOI | https://doi.org/10.1051/0004-6361/201014599 | |
Published online | 16 July 2010 |
Herschel: the first science highlights
LETTER TO THE EDITOR
Predicted colors and flux densities of protostars in the Herschel
PACS and SPIRE filters
B. Ali1 - J. J. Tobin2 - W. J. Fischer3
- C. A. Poteet3 - S. T. Megeath3 - L. Allen4 - L. Hartmann2 - N. Calvet2 - E. Furlan5,
- M. Osorio6
1 - NHSC/IPAC, California Institute of Technology,
Pasadena, CA USA
2 - University of Michigan, Ann Arbor, Michigan USA
3 - University of Toledo, Toledo, OH USA
4 - NOAO, Tucson, AZ USA
5 - JPL, Caltech, Pasadena, CA, USA
6 - Instituto de Astrofisica de Andalucia, Granada, Spain
Received 31 March 2010 / Accepted 27 April 2010
Abstract
Upcoming surveys with the Herschel Space Observatory will
yield far-IR photometry of large samples of young stellar objects,
which will require careful interpretation. We investigate the color and
luminosity diagnostics based on Herschel broad-band filters
to identify and discern the properties of low-mass protostars. We
compute a grid of 2016 protostars in various physical
congurations, present the expected flux densities and flux density
ratios for this grid of protostars, and compare Herschel
observations of three protostars to the model results. These provide
useful constraints on the range of colors and fluxes of protostar in
the Herschel filters. We find that Herschel data alone is likely a useful diagnostic of the envelope properties of young stars.
Key words: stars: protostars - radiative transfer - submillimeter: stars
1 Introduction
A comprehensive theory of star formation is necessary for understanding processes spanning the entire range of cosmic evolution, from the formation and evolution of galaxies to the formation of planets. Thus, observing the formation of stars is one of the prime science objectives of the Herschel Space Observatory Pilbrat et al. 2010. With its large (3.5 m) primary mirror, Herschel is able to survey a substantial fraction of the current star formation in the Galaxy with unprecedented efficiency and sensitivity. The result will be an explosion of far-IR and sub-mm data on star formation regions. Herschel is particularly suited for carrying out these surveys because the spectral energy distributions (SEDs) of protostars peaks in the far-IR and sub-mm region. A large number of guaranteed- and open- time key (legacy) Herschel programs are, thus, focused on star formation regions. These include PACS (Photodetector Array Camera and Spectrometer, Poglitsch et al. 2010) and SPIRE (Spectral and Imaging Reciever, Griffin et al. 2010) imaging surveys of the Gould Belt (André et al. 2010), and the Herschel Orion Protostar Survey (HOPS) and DIGIT programs, among others.
These surveys will detect hundreds of thousands of protostars. For our
purposes, protostars are low-mass stars in the earliest stages of their
evolution. They are less than a few 105 yrs
old and are still surrounded by an infalling envelope. The SEDs of
these protostars will be complex mixture of their intrinsic emission
coupled with reprocessing and emission from their immediate
environment. Radiative transfer models are an important tool in
understanding and interpretting the observed flux densities of
protostars. Model SEDs can provide crucial constraints on the most
important physical parameters of a protostellar system. Because
protostars can encompass a wide rage of parameter space it is often
necessary to compare observations against a large number of different
models to derive a reasonable solution and to understand the
degeneracies involved. Robitaille et al. (2006)
constructed a grid of 20 000 individual models for SED
fitting; however, these models used an ism model for the envelope,
which is not accurate at far-IR wavelengths. This is because
reprocessing by dust is especially important at Herschel's
wavelengths as this reprocessed light forms the bulk of the protostars'
SED. In this contribution, we have generated a grid of
20 160 SEDs from 2016 individual protostar
models focused primarily on the parameter space appropriate for young stars, and for the Herschel broad-band filters.
2 The model grid and flux predictions
We use the Monte Carlo radiative transfer code developed by Whitney et al. (2003), which is adapted to include the dust mixture properties from Tobin et al.(2008) to calculate the model SEDs. For convenience, we briefly summarize the main components here. The code comprises four physical components: (i) a central object; (ii) a flared disk; (iii) a rotating collapsing envelope; and (iv) a bipolar cavity. These components are illustrated in Fig. 1. The star+disk system emits energy packets (``photons'') that are absorbed or scattered by the disk and/or envelope, subsequently determining the emergent SED of the system. Our adopted dust model uses the methodology described by D'Alessio et al. (2001) with additional optical constants for graphite from Draine & Lee (1984). The dust grains used here are a mixture of graphite, silicate and water ice and their number density has a power-law distribution with grain size to the power -3.5 and include large (1 micron radius) grains; the implied gas-to-dust ratio is 133.
The code was used to produce a grid of 20 160 model SEDs using
photons
for the Monte Carlo computations, and by using various configurations
for the main components listed above. We restricted the variable
parameter set to: the infall rate (
), the maximum envelope size (
), the cavity opening angle (
), the cavity shape exponent (b), the accretion rate (
), the size of the disk (
),
and the inclination angle of the system to the observer's line of
sight (i). These are the critical parameters that predominantly
determine the flux density prediction for Herschel PACS and SPIRE bands (56-610
m). The remaining parameters have fixed values that are typical for low-mass protostars (Tobin et al. 2008).
One important characteristic of the code is that the disk accretion
rate is used in the model grid to vary the luminosity of the
protostellar systems; the central star is fixed at 1 .
This is reasonable since we do not know the relative contribution of
photospheric and accretion luminosities in protostars. However, the
spectrum of the central object's emission does not largely affect the
emergent spectrum of the protostar in the far-IR since nearly all
emission is reprocessed. Thus, dust opacities are a more critical and
relevant factor for our analysis.
Table 1 summarizes the relevant model parameters and their values in our grid. We list only select parameters (see Tobin et al. 2008
for a complete list). The properties of the central object are perhaps
the most crucial of the fixed parameters in that changing the
parameters for central object may lead to degeneracies in the observed
colors and fluxes. We discuss this further in Sect. 3.5.
Other relevant limitations of the model are that dust opacities, which
are known to vary across regions, are assumed to be constant, and that
only the water ice opacities for the 3.6 m feature are included in addition to dust.
2.1 Synthetic photometry
The code computes model SEDs for a given spatial extent (aperture) of the system. We use aperture radius of 20 000 AU. This collects the total integrated flux from the envelopes which are all smaller than the aperture. We convolved the resulting SEDs with the filter+detector transmission curves of PACS and SPIRE (kindly provided by their respective instrument teams). We present the results as flux density values at the central wavelengths of the filters. The model fluxes are calculated at the distance of Orion, adopted to be 420 pc (Menten et al. 2007).
![]() |
Figure 1: Illustration of some of the major parameters used to model the protostar SEDs. |
Open with DEXTER |
Table 1: Model parameters.
3 Results and discussion
Figures 2-6 show the results for the 20 160 model SEDs calculated for our grid for various combinations of flux densities and flux density ratios. The shaded dots shows the full range of predicted flux densities for the entire set of 2016 protostar systems as seen with 10 different inclination angles each. The parameter values used here sample a broad range of physical conditions in the protostars and, thus, provide a reasonable comparative range of expected flux densities and flux density ratios for protostars.
In protostellar systems, model grids, such as the one discussed here, are most accurate and useful when multi-wavelength data are available to provide as many constraints as possible (see e.g. Fischer et al. 2010). Nonetheless, PACS and SPIRE sample the peak of the protostar's SED and we expect Herschel data alone to provide a relatively good estimate of the envelope properties in particular.
We further investigate the effect of each of the variable parameters on
the observed fluxes. To illustrate how a particular model parameter
affects the flux densities, we first selected a typical protostar
system, here defined to have the following values:
= 10-7
yr-1,
= 100 AU,
= 10-5
yr-1,
= 10 000 AU,
= 25 degrees, b=1.5, and i=63.3 degrees.
Then, we show how individual parameters affect the results as they vary
over the range listed in Table 1,
while holding all other parameters constant. The results are
color-coded (as listed in the legend boxes in the upper left hand part
of Figs. 2-6)
and shown as thick dashed lines. We show and discuss below results from
a representative set of correlations investigated from the full grid of
model parameter range. We omit discussion on
,
as this parameter and its use in the code is already discussed above (see Sect. 2).
![]() |
Figure 2:
Predicted results for the PACS 70 |
Open with DEXTER |
![]() |
Figure 3: Same as Fig. 2 but this time for the 3 parameters that appear to produce little or no effect on the observed flux densities. This figure is a magnified portion of the range of values shown in Fig. 2. The effect of these parameters will not be visible if shown over the full range displayed in Fig. 2. |
Open with DEXTER |
![]() |
Figure 4:
The ratio of the PACS flux densities, 70 |
Open with DEXTER |
![]() |
Figure 5:
The SPIRE 250 |
Open with DEXTER |
![]() |
Figure 6:
The PACS (70 |
Open with DEXTER |
3.1 Envelope infall rate
The mass infall rate strongly affects the far-IR flux and not
surprisingly, the predicted flux density and their ratios typically
vary by factors 2 as a function of
.
This is straight forward as increasing the amount of infall simply
increases the amount of cold material and thus the far-IR fluxes. We
expect to be able to constrain the envelope infall rates with
relatively good precision with Herschel data alone, although, some knowledge about the distance to the source (or its luminosity) is required, or must be assumed.
3.2 Envelope size
The size of the envelope directly affects the amount of the coldest
dust present in the protostar. Thus, we note significantly more
variation in the longer wavelength SPIRE filters (see Figs. 5-6)
than in the PACS filters. Increasing the outer radius of the envelope
preferentially creates a comparatively redder protostellar system.
However, the affect of the
is still appreciably lower than
.
That's because the latter more directly controls the envelope density
and total amount of emitting material present in the envelope.
3.3 Cavity opening angle
The PACS fluxes appears to be most strongly affected by this parameter (see Figs. 2 and 4).
This is likely because the width of the cavity preferentially carves
away the part of the envelope responsible for the 70 m flux. By contrast, there is almost no variation as a function of
in the SPIRE filters.
3.4 Inclination angle
The optical and near-IR emission from protostellar systems is dominated by the source itself and the scattered light from the circumstellar envelope. Thus, face-on models are typically brighter and bluer compared to those in which the central star is either wholly or partially obscured by the disk. At optical and near-IR wavelengths, these objects are quite often extremely faint or undetected even in unbiased photometry surveys due to high obscurataion. Figure 3 shows that in the PACS filters this effect becomes largely unimportant as the main emission is from the envelope. We find similar results for the other combinations of flux and flux density diagnostics (not shown).
3.5 The central object
We have concentrated on a single central star; however, the models are applicable to all low mass stars. The observed fluxes and SEDs simply scale with the mass of the central object as determined by the envelope density profile and the accretion luminosity equations. When the mass of the central star is changed, the envelope density remains the same but the infall rate changes. And, the mass accretion rate for a given luminosity scales inversely with the central object mass. Thus, considerable overlap does exist in the observed fluxes for different mass protostars. This result further demonstrates that PACS and SPIRE are best suited for constraining the envelope and disk properties and provide limited constraints on the central object itself.
3.6 The HOPS protostars sample
The HOPS program (Fischer et al. 2010) imaged 4 protostars in Orion L 1641 region for which PACS 70 and 160 m photometry is available for 3 sources. These stars are shown with the diamond symbol in Fig. 2.
Based on this diagnostic alone, we can conclude that all 3 HOPS
sources are much more luminous than our typical protostar of 1
.
Only 3 sources with two-band photometry are in the HOPS field;
nonetheless, the sources are not clustered but spread over a broad
range of color and flux density values in Fig. 2,
implying that a broad range of values must also exist for the other
parameters. These conclusions are indeed confirmed by detailed model
SED fits over the 1-870
m range for the HOPS sources (Fischer et al. 2010).
4 Conclusions
The model grid presented here provides the expected flux densities and flux density ratios for 20 160 low-mass protostellar SEDs, useful for determining the phase space occupied by protostar in the Herschel filters. Disentangling the contributions of the various physical components that make up a protostellar system is observationally challenging and requires multi-wavelength data to properly constrain the properties of the system. Nonetheless, preliminary results show that the fluxes in the Herschel PACS and SPIRE filters are dominated by reprocessed light and emission from the protostellar envelope. Herschel data alone is, therefore, expected to be able to discern envelope properties with some precision provided that a reasonable estimate for the object's distance or its luminosity is available. The source of the emission (combination of central star and disk accretion) is relatively unimportant at Herschel far-IR and sub-mm wavelengths. Thus, Herschel data is unlikely to constrain the mass of the central object, but sampling the peak of the protostars SED allows tight constrains on its total luminosity.
AcknowledgementsWe thank B. Whitney for making her code available publicly. B. Ali acknowledges support from NASA grant IPAC.ALI-OTKP-1-JPL.000094. J. J. Tobin acknowledges support from NASA grant HST-GO-11548.04-A.
References
- André, Ph., et al. 2010, A&A, 518, L102 D'Alessio, P., Calvet, N., & Hartmann, L. 2001, ApJ, 553, 321 [NASA ADS] [CrossRef] [Google Scholar]
- Drain, B. T., & Lee, H. M. 1984, ApJ, 285, 89 [Google Scholar]
- Fischer, W., J. et al. 2010, A&A, 518, L122 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Griffin, M. J., et al. 2010, A&A, 518, L3 [Google Scholar]
- Menten, K. M., Reid, M. J., Forbrich, J., & Brunthaler, A. 2007, A&A, 474, 515 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Pilbratt, G. L., et al. 2010, A&A, 518, L1 [CrossRef] [EDP Sciences] [Google Scholar]
- Poglitsch, A., et al. 2010, A&A, 518, L2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Robitaille, T. P., Whitney, B. A., Indebetouw, R., Wood, K., & Denzmore, P. 2006, ApJS, 167, 256 [NASA ADS] [CrossRef] [Google Scholar]
- Whitney, B. A., Wood, K., Bjorkman, J. E., & Cohen, M. 2003, ApJ, 598 [Google Scholar]
- Tobin, J. J., Hartmann, L., Calvet, N., & D'Alessio, P. 2008, ApJ, 679, 1364 [NASA ADS] [CrossRef] [Google Scholar]
- Whitney, B. A., Wood, K., Bjorkman, J. E., & Wolff, M. J. 2003 ApJ, 591, 1049 [Google Scholar]
Footnotes
- ...Herschel
- Herschel is an ESA space observatory with science instruments provided by Principal Investigator consortia. It is open for proposals for observing time from the worldwide astronomical community.
- ...
- Spitzer Fellow.
- ... protostar
- There are 10 different inclination angles for each of 2016 protostars.
All Tables
Table 1: Model parameters.
All Figures
![]() |
Figure 1: Illustration of some of the major parameters used to model the protostar SEDs. |
Open with DEXTER | |
In the text |
![]() |
Figure 2:
Predicted results for the PACS 70 |
Open with DEXTER | |
In the text |
![]() |
Figure 3: Same as Fig. 2 but this time for the 3 parameters that appear to produce little or no effect on the observed flux densities. This figure is a magnified portion of the range of values shown in Fig. 2. The effect of these parameters will not be visible if shown over the full range displayed in Fig. 2. |
Open with DEXTER | |
In the text |
![]() |
Figure 4:
The ratio of the PACS flux densities, 70 |
Open with DEXTER | |
In the text |
![]() |
Figure 5:
The SPIRE 250 |
Open with DEXTER | |
In the text |
![]() |
Figure 6:
The PACS (70 |
Open with DEXTER | |
In the text |
Copyright ESO 2010
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.