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.Comment: Accepted for publication in A&

A Water Maser and Ammonia Survey of GLIMPSE Extended Green Objects (EGOs)

We present the results of a Nobeyama 45-m water maser and ammonia survey of all 94 northern GLIMPSE Extended Green Objects (EGOs), a sample of massive young stellar objects (MYSOs) identified based on their extended 4.5 micron emission. We observed the ammonia (1,1), (2,2), and (3,3) inversion lines, and detect emission towards 97%, 63%, and 46% of our sample, respectively (median rms ~50 mK). The water maser detection rate is 68% (median rms ~0.11 Jy). The derived water maser and clump-scale gas properties are consistent with the identification of EGOs as young MYSOs. To explore the degree of variation among EGOs, we analyze subsamples defined based on MIR properties or maser associations. Water masers and warm dense gas, as indicated by emission in the higher-excitation ammonia transitions, are most frequently detected towards EGOs also associated with both Class I and II methanol masers. 95% (81%) of such EGOs are detected in water (ammonia(3,3)), compared to only 33% (7%) of EGOs without either methanol maser type. As populations, EGOs associated with Class I and/or II methanol masers have significantly higher ammonia linewidths, column densities, and kinetic temperatures than EGOs undetected in methanol maser surveys. However, we find no evidence for statistically significant differences in water maser properties (such as maser luminosity) among any EGO subsamples. Combining our data with the 1.1 mm continuum Bolocam Galactic Plane Survey, we find no correlation between isotropic water maser luminosity and clump number density. Water maser luminosity is weakly correlated with clump (gas) temperature and clump mass.Comment: Astrophysical Journal, accepted. Emulateapj, 24 pages including 24 figures, plus 9 tables (including full content of online-only tables

The RMS Survey: The Bolometric Fluxes and Luminosity Distributions of Young Massive Stars

Context: The Red MSX Source (RMS) survey is returning a large sample of massive young stellar objects (MYSOs) and ultra-compact (UC) \HII{} regions using follow-up observations of colour-selected candidates from the MSX point source catalogue. Aims: To obtain the bolometric fluxes and, using kinematic distance information, the luminosities for young RMS sources with far-infrared fluxes. Methods: We use a model spectral energy distribution (SED) fitter to obtain the bolometric flux for our sources, given flux data from our work and the literature. The inputs to the model fitter were optimised by a series of investigations designed to reveal the effect varying these inputs had on the resulting bolometric flux. Kinematic distances derived from molecular line observations were then used to calculate the luminosity of each source. Results: Bolometric fluxes are obtained for 1173 young RMS sources, of which 1069 have uniquely constrained kinematic distances and good SED fits. A comparison of the bolometric fluxes obtained using SED fitting with trapezium rule integration and two component greybody fits was also undertaken, and showed that both produce considerable scatter compared to the method used here. Conclusions: The bolometric flux results allowed us to obtain the luminosity distributions of YSOs and UC\HII{} regions in the RMS sample, which we find to be different. We also find that there are few MYSOs with L $\geq$ 10$^{5}$\lsol{}, despite finding many MYSOs with 10$^{4}$\lsol{} $\geq$ L $\geq$ 10$^{5}$\lsol{}.Comment: 12 pages, 12 figures, 3 tables, accepted to A&A. The full versions of tables 1 and 2 will be available via the CDS upon publicatio

Infrared point source variability between the Spitzer and MSX surveys of the Galactic mid-plane

We present a list of 552 sources with suspected variability, based on a comparison of mid-infrared photometry from the GLIMPSE I and MSX surveys, which were carried out nearly a decade apart. We were careful to address issues such as the difference in resolution and sensitivity between the two surveys, as well as the differences in the spectral responses of the instruments. We selected only sources where the IRAC 8.0 and MSX 8.28 micron fluxes differ by more than a factor of two, in order to minimize contamination from sources where the difference in fluxes at 8 micron is due to a strong 10 micron silicate feature. We present a subset of 40 sources for which additional evidence suggests variability, using 2MASS and MIPSGAL data. Based on a comparison with the variability flags in the IRAS and MSX Point-Source Catalogs we estimate that at least a quarter of the 552 sources, and at least half of the 40 sources are truly variable. In addition, we tentatively confirm the variability of one source using multi-epoch IRAS LRS spectra. We suggest that most of the sources in our list are likely to be Asymptotic Giant Branch stars.Comment: 47 pages, 12 Figures, 3 Tables, accepted for publication in A

The morphology of minor axis gaseous outflows in edge-on Seyfert galaxies

(Abridged) Spiral galaxies often have extended outflows that permeate beyond the region of the disk. Such outflows have been seen both in starburst galaxies, actively star forming galaxies and galaxies with an AGN. In the latter galaxies it is unknown whether the large-scale outflows are driven by star formation activity or purely by the active nucleus. The aim of our investigation is to study the frequency of extended minor-axis outflows in edge-on Seyfert galaxies to investigate the role of the AGN, the circumnuclear environment and star formation activity within the disk regions, and their importance for IGM enrichment on large scales. Narrowband imaging in two different ionizational stages (H-alpha and [OIII]) was performed to attempt a discrimination between processes associated with the active nucleus and those connected to star forming activity within the disk. The H-alpha morphology of the Seyfert galaxies is usually complex, but only in three out of 14 galaxies did we find evidence for minor axis disk outflows. At the sensitivity of our observations [OIII] emission is generally detected only in the nuclear region. Overall, our results show that extraplanar emission of similar brightness and extent as in the previously known cases of NGC3079 and NGC4388 is not common in Seyfert galaxies of otherwise similar properties. Comparison with our previous results shows that for nearby edge-on spiral galaxies star formation may be a more powerful mechanism for producing DIG than AGN activity. While in general AGN activity undoubtedly plays some role in driving minor-axis outflows, this probably requires higher AGN luminosities than are encountered in our small distance-limited sample.Comment: accepted for publication in Astronomy & Astrophysics, 36 pages, 25 figures, 5 tables (Figs. 6-24 and Tables 2 + 5 in the online part), Article with Full resolution images is available at http://www.stsci.edu/~marel/psgzdir/seyfert_outflows.ps.g

Sky maps without anisotropies in the cosmic microwave background are a better fit to WMAP's uncalibrated time ordered data than the official sky maps

The purpose of this reanalysis of the WMAP uncalibrated time ordered data (TOD) was two fold. The first was to reassess the reliability of the detection of the anisotropies in the official WMAP sky maps of the cosmic microwave background (CMB). The second was to assess the performance of a proposed criterion in avoiding systematic error in detecting a signal of interest. The criterion was implemented by testing the null hypothesis that the uncalibrated TOD was consistent with no anisotropies when WMAP's hourly calibration parameters were allowed to vary. It was shown independently for all 20 WMAP channels that sky maps with no anisotropies were a better fit to the TOD than those from the official analysis. The recently launched Planck satellite should help sort out this perplexing result.Comment: 11 pages with 1 figure and 2 tables. Extensively rewritten to explain the research bette

Duration of Star Formation in Galactic Giant Molecular Clouds. I. The Great Nebula in Carina

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

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.Comment: 29 pages, 9 figures, accepted to Ap