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
A&A 518, L119 (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 $\times$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 $2 \times 10^7$ 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 ( $\dot{M}_{\rm env}$), the maximum envelope size ( $R_{\rm env,\max}$), the cavity opening angle ( $\theta_{C}$), the cavity shape exponent (b), the accretion rate ( $\dot{M}_{\rm disk}$), the size of the disk ( $R_{\rm disk,\max}$), 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 $\mu $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 $L_{\odot}$. 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 $\mu $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).

\begin{figure}
\par\includegraphics[width=8cm,clip]{14599fg1.eps}
\end{figure} 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: $\dot{M}_{\rm disk}$ = 10-7 $M_{\odot}$ yr-1, $R_{\rm disk,\max}$ = 100 AU, $\dot{M}_{\rm env}$ = 10-5 $M_{\odot}$ yr-1, $R_{\rm env,\max}$ = 10 000 AU, $\theta_{C}$ = 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 $\dot{M}_{\rm disk}$, as this parameter and its use in the code is already discussed above (see Sect. 2).

\begin{figure}
\par\includegraphics[width=8.3cm,clip]{14599fg2.eps}
\end{figure} Figure 2:

Predicted results for the PACS 70 $\mu $m filter vs. the ratio of 70 to 160 $\mu $m flux density. The fluxes are reported for a protostar at the distance of Orion. The shaded region shows the full range in the 20 160 SEDs in the grid. The color-coded (inset) thick dashed lines illustrate the magnitude and direction of the parameter when all other parameters are held constant. We do this by choosing typical values (see text) and plotting the variation from only one parameter for that protostar. The diamond symbols show the locations of the 3 HOPS protostars (see Sect. 3.6).

Open with DEXTER

\begin{figure}
\par\includegraphics[width=8cm,clip]{14599fg3.eps}
\end{figure} 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

\begin{figure}
\par\includegraphics[width=8cm,clip]{14599fg4.eps}
\end{figure} Figure 4:

The ratio of the PACS flux densities, 70 $\mu $m/100 $\mu $m vs. 100 $\mu $m/160 $\mu $m. The symbols and lines are as shown in Fig. 2.

Open with DEXTER

\begin{figure}
\par\includegraphics[width=8.5cm,clip]{14599fg5.eps}
\end{figure} Figure 5:

The SPIRE 250 $\mu $m flux density vs. the ratio 250 $\mu $m/350 $\mu $m. The symbols and lines are as in Fig. 2.

Open with DEXTER

\begin{figure}
\par\includegraphics[width=8.2cm,clip]{14599fg6.eps}
\end{figure} Figure 6:

The PACS (70 $\mu $m/160 $\mu $m) ratio vs. the SPIRE (250 $\mu $m/350 $\mu $m) ratio.

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 $\sim$2 as a function of $\dot{M}_{\rm env}$. 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 $R_{\rm env,\max}$ is still appreciably lower than $\dot{M}_{\rm env}$. 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 $\mu $m flux. By contrast, there is almost no variation as a function of $\theta_{C}$ 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 $\mu $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 $L_{\odot}$. 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 $\mu $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.

Acknowledgements
We 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

  1. 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]
  2. Drain, B. T., & Lee, H. M. 1984, ApJ, 285, 89 [Google Scholar]
  3. Fischer, W., J. et al. 2010, A&A, 518, L122 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Griffin, M. J., et al. 2010, A&A, 518, L3 [Google Scholar]
  5. Menten, K. M., Reid, M. J., Forbrich, J., & Brunthaler, A. 2007, A&A, 474, 515 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Pilbratt, G. L., et al. 2010, A&A, 518, L1 [CrossRef] [EDP Sciences] [Google Scholar]
  7. Poglitsch, A., et al. 2010, A&A, 518, L2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Robitaille, T. P., Whitney, B. A., Indebetouw, R., Wood, K., & Denzmore, P. 2006, ApJS, 167, 256 [NASA ADS] [CrossRef] [Google Scholar]
  9. Whitney, B. A., Wood, K., Bjorkman, J. E., & Cohen, M. 2003, ApJ, 598 [Google Scholar]
  10. Tobin, J. J., Hartmann, L., Calvet, N., & D'Alessio, P. 2008, ApJ, 679, 1364 [NASA ADS] [CrossRef] [Google Scholar]
  11. 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

  \begin{figure}
\par\includegraphics[width=8cm,clip]{14599fg1.eps}
\end{figure} Figure 1:

Illustration of some of the major parameters used to model the protostar SEDs.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=8.3cm,clip]{14599fg2.eps}
\end{figure} Figure 2:

Predicted results for the PACS 70 $\mu $m filter vs. the ratio of 70 to 160 $\mu $m flux density. The fluxes are reported for a protostar at the distance of Orion. The shaded region shows the full range in the 20 160 SEDs in the grid. The color-coded (inset) thick dashed lines illustrate the magnitude and direction of the parameter when all other parameters are held constant. We do this by choosing typical values (see text) and plotting the variation from only one parameter for that protostar. The diamond symbols show the locations of the 3 HOPS protostars (see Sect. 3.6).

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=8cm,clip]{14599fg3.eps}
\end{figure} 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

  \begin{figure}
\par\includegraphics[width=8cm,clip]{14599fg4.eps}
\end{figure} Figure 4:

The ratio of the PACS flux densities, 70 $\mu $m/100 $\mu $m vs. 100 $\mu $m/160 $\mu $m. The symbols and lines are as shown in Fig. 2.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=8.5cm,clip]{14599fg5.eps}
\end{figure} Figure 5:

The SPIRE 250 $\mu $m flux density vs. the ratio 250 $\mu $m/350 $\mu $m. The symbols and lines are as in Fig. 2.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=8.2cm,clip]{14599fg6.eps}
\end{figure} Figure 6:

The PACS (70 $\mu $m/160 $\mu $m) ratio vs. the SPIRE (250 $\mu $m/350 $\mu $m) ratio.

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.