An all-sky Support Vector Machine selection of WISE YSO Candidates

We explored the AllWISE catalogue of the Wide-field Infrared Survey Explorer mission and identified Young Stellar Object candidates. Reliable 2MASS and WISE photometric data combined with Planck dust opacity values were used to build our dataset and to find the best classification scheme. A sophisticated statistical method, the Support Vector Machine (SVM) is used to analyse the multi-dimensional data space and to remove source types identified as contaminants (extragalactic sources, main sequence stars, evolved stars and sources related to the interstellar medium). Objects listed in the SIMBAD database are used to identify the already known sources and to train our method. A new all-sky selection of 133,980 Class I/II YSO candidates is presented. The estimated contamination was found to be well below 1% based on comparison with our SIMBAD training set. We also compare our results to that of existing methods and catalogues. The SVM selection process successfully identified >90% of the Class I/II YSOs based on comparison with photometric and spectroscopic YSO catalogues. Our conclusion is that by using the SVM, our classification is able to identify more known YSOs of the training sample than other methods based on colour-colour and magnitude-colour selection. The distribution of the YSO candidates well correlates with that of the Planck Galactic Cold Clumps in the Taurus--Auriga--Perseus--California region.Comment: 27 pages, 15 figures, 15 table

A Broadband Study of Galactic Dust Emission

We have combined infrared data with HI, H2 and HII surveys in order to spatially decompose the observed dust emission into components associated with different phases of the gas. An inversion technique is applied. For the decomposition, we use the IRAS 60 and 100 micron bands, the DIRBE 140 and 240 micron bands, as well as Archeops 850 and 2096 micron wavelengths. In addition, we apply the decomposition to all five WMAP bands. We obtain longitude and latitude profiles for each wavelength and for each gas component in carefully selected Galactic radius bins.We also derive emissivity coefficients for dust in atomic, molecular and ionized gas in each of the bins.The HI emissivity appears to decrease with increasing Galactic radius indicating that dust associated with atomic gas is heated by the ambient interstellar radiation field (ISRF). By contrast, we find evidence that dust mixed with molecular clouds is significantly heated by O/B stars still embedded in their progenitor clouds. By assuming a modified black-body with emissivity law lambda^(-1.5), we also derive the radial distribution of temperature for each phase of the gas. All of the WMAP bands except W appear to be dominated by emission from something other than normal dust, most likely a mixture of thermal bremstrahlung from diffuse ionized gas, synchrotron emission and spinning dust. Furthermore, we find indications of an emissivity excess at long wavelengths (lambda > 850 micron) in the outer Galaxy (R > 8.9 kpc). This suggests either the existence of a very cold dust component in the outer Galaxy or a temperature dependence of the spectral emissivity index. Finally, it is shown that ~ 80% of the total FIR luminosity is produced by dust associated with atomic hydrogen, in agreement with earlier findings by Sodroski et al. (1997).Comment: accepted for publication by A&

Grain size limits derived from 3.6 {\mu}m and 4.5 {\mu}m coreshine

Recently discovered scattered light from molecular cloud cores in the wavelength range 3-5 {\mu}m (called "coreshine") seems to indicate the presence of grains with sizes above 0.5 {\mu}m. We aim to analyze 3.6 and 4.5 {\mu}m coreshine from molecular cloud cores to probe the largest grains in the size distribution. We analyzed dedicated deep Cycle 9 Spitzer IRAC observations in the 3.6 and 4.5 {\mu}m bands for a sample of 10 low-mass cores. We used a new modeling approach based on a combination of ratios of the two background- and foreground-subtracted surface brightnesses and observed limits of the optical depth. The dust grains were modeled as ice-coated silicate and carbonaceous spheres. We discuss the impact of local radiation fields with a spectral slope differing from what is seen in the DIRBE allsky maps. For the cores L260, ecc806, L1262, L1517A, L1512, and L1544, the model reproduces the data with maximum grain sizes around 0.9, 0.5, 0.65, 1.5, 0.6, and > 1.5 {\mu}m, respectively. The maximum coreshine intensities of L1506C, L1439, and L1498 in the individual bands require smaller maximum grain sizes than derived from the observed distribution of band ratios. Additional isotropic local radiation fields with a spectral shape differing from the DIRBE map shape do not remove this discrepancy. In the case of Rho Oph 9, we were unable to reliably disentangle the coreshine emission from background variations and the strong local PAH emission. Considering surface brightness ratios in the 3.6 and 4.5 {\mu}m bands across a molecular cloud core is an effective method of disentangling the complex interplay of structure and opacities when used in combination with observed limits of the optical depth.Comment: 23 pages, 18 figures, accepted for publication in A&

Modeling and predicting the shape of the far-infrared to submillimeter emission in ultra-compact HII regions and cold clumps

Dust properties are very likely affected by the environment in which dust grains evolve. For instance, some analyses of cold clumps (7 K- 17 K) indicate that the aggregation process is favored in dense environments. However, studying warm (30 K-40 K) dust emission at long wavelength ($\lambda$$>$300 $\mu$m) has been limited because it is difficult to combine far infared-to-millimeter (FIR-to-mm) spectral coverage and high angular resolution for observations of warm dust grains. Using Herschel data from 70 to 500 $\mu$m, which are part of the Herschel infrared Galactic (Hi-GAL) survey combined with 1.1 mm data from the Bolocam Galactic Plane Survey (BGPS), we compared emission in two types of environments: ultra-compact HII (UCHII) regions, and cold molecular clumps (denoted as cold clumps). With this comparison we tested dust emission models in the FIR-to-mm domain that reproduce emission in the diffuse medium, in these two environments (UCHII regions and cold clumps). We also investigated their ability to predict the dust emission in our Galaxy. We determined the emission spectra in twelve UCHII regions and twelve cold clumps, and derived the dust temperature (T) using the recent two-level system (TLS) model with three sets of parameters and the so-called T-$\beta$ (temperature-dust emissvity index) phenomenological models, with $\beta$ set to 1.5, 2 and 2.5. We tested the applicability of the TLS model in warm regions for the first time. This analysis indicates distinct trends in the dust emission between cold and warm environments that are visible through changes in the dust emissivity index. However, with the use of standard parameters, the TLS model is able to reproduce the spectral behavior observed in cold and warm regions, from the change of the dust temperature alone, whereas a T-$\beta$ model requires $\beta$ to be known.Comment: Accepted for publication in A&A. 19 pages, 8 figures, 7 table

The Luminous Blue Variable RMC127 as seen with ALMA and ATCA

We present ALMA and ATCA observations of the luminous blue variable \rmc. The radio maps show for the first time the core of the nebula and evidence that the nebula is strongly asymmetric with a Z-pattern shape. Hints of this morphology are also visible in the archival \emph{HST} $\rm H\alpha$ image, which overall resembles the radio emission. The emission mechanism in the outer nebula is optically thin free-free in the radio. At high frequencies, a component of point-source emission appears at the position of the star, up to the ALMA frequencies. The rising flux density distribution ($S_{\nu}\sim \nu^{0.78\pm0.05}$) of this object suggests thermal emission from the ionized stellar wind and indicates a departure from spherical symmetry with $n_{e}(r)\propto r^{-2}$. We examine different scenarios to explain this excess of thermal emission from the wind and show that this can arise from a bipolar outflow, supporting the suggestion by other authors that the stellar wind of \rmc is aspherical. We fit the data with two collimated ionized wind models and we find that the mass-loss rate can be a factor of two or more smaller than in the spherical case. We also fit the photometry obtained by IR space telescopes and deduce that the mid- to far-IR emission must arise from extended, cool ($\sim80\,\rm K$) dust within the outer ionized nebula. Finally we discuss two possible scenarios for the nebular morphology: the canonical single star expanding shell geometry, and a precessing jet model assuming presence of a companion star.Comment: Accepted for publication in ApJ (minor revision included

The spiral structure of our Milky Way Galaxy

The spiral structure of our Milky Way Galaxy is not yet known. HII regions and giant molecular clouds are the most prominent spiral tracers. We collected the spiral tracer data of our Milky Way from the literature, namely, HII regions and giant molecular clouds (GMCs). With weighting factors based on the excitation parameters of HII regions or the masses of GMCs, we fitted the distribution of these tracers with models of two, three, four spiral-arms or polynomial spiral arms. The distances of tracers, if not available from stellar or direct measurements, were estimated kinetically from the standard rotation curve of Brand & Blitz (1993) with $R_0$=8.5 kpc, and $\Theta_0$=220 km s$^{-1}$ or the newly fitted rotation curves with $R_0$=8.0 kpc and $\Theta_0$=220 km s$^{-1}$ or $R_0$=8.4 kpc and $\Theta_0$=254 km s$^{-1}$. We found that the two-arm logarithmic model cannot fit the data in many regions. The three- and the four-arm logarithmic models are able to connect most tracers. However, at least two observed tangential directions cannot be matched by the three- or four-arm model. We composed a polynomial spiral arm model, which can not only fit the tracer distribution but also match observed tangential directions. Using new rotation curves with $R_0$=8.0 kpc and $\Theta_0$=220 km s$^{-1}$ and $R_0$=8.4 kpc and $\Theta_0$=254 km s$^{-1}$ for the estimation of kinematic distances, we found that the distribution of HII regions and GMCs can fit the models well, although the results do not change significantly compared to the parameters with the standard $R_0$ and $\Theta_0$.Comment: 34 Pages, 10 Figures, 5 Tables. Accepted for publication in A&A. Edited

Unexpectedly large mass loss during the thermal pulse cycle of the red giant R Sculptoris!

The asymptotic giant branch star R Sculptoris is surrounded by a detached shell of dust and gas. The shell originates from a thermal pulse during which the star undergoes a brief period of increased mass loss. It has hitherto been impossible to constrain observationally the timescales and mass-loss properties during and after a thermal pulse - parameters that determine the lifetime on the asymptotic giant branch and the amount of elements returned by the star. Here we report observations of CO emission from the circumstellar envelope and shell around R Sculptoris with an angular resolution of 1.3 arcsec. What was hitherto thought to be only a thin, spherical shell with a clumpy structure, is revealed to contain a spiral structure. Spiral structures associated with circumstellar envelopes have been seen previously, from which it was concluded that the systems must be binaries. Using the data, combined with hydrodynamic simulations, we conclude that R Sculptoris is a binary system that underwent a thermal pulse approximately 1800 years ago, lasting approximately 200 years. About 0.003 Msun of mass was ejected at a velocity of 14.3 km s-1 and at a rate approximately 30 times higher than the prepulse mass-loss rate. This shows that approximately 3 times more mass is returned to the interstellar medium during and immediately after a pulse than previously thought.Comment: Accepted by Natur

Large-scale environments of binary AGB stars probed by Herschel. II: Two companions interacting with the wind of pi1 Gruis

Context. The Mass loss of Evolved StarS (MESS) sample observed with PACS on board the Herschel Space Observatory revealed that several asymptotic giant branch (AGB) stars are surrounded by an asymmetric circumstellar envelope (CSE) whose morphology is most likely caused by the interaction with a stellar companion. The evolution of AGB stars in binary systems plays a crucial role in understanding the formation of asymmetries in planetary nebul{\ae} (PNe), but at present, only a handful of cases are known where the interaction of a companion with the stellar AGB wind is observed. Aims. We probe the environment of the very evolved AGB star $\pi^1$ Gruis on large and small scales to identify the triggers of the observed asymmetries. Methods. Observations made with Herschel/PACS at 70 $\mu$m and 160 $\mu$m picture the large-scale environment of $\pi^1$ Gru. The close surroundings of the star are probed by interferometric observations from the VLTI/AMBER archive. An analysis of the proper motion data of Hipparcos and Tycho-2 together with the Hipparcos Intermediate Astrometric Data help identify the possible cause for the observed asymmetry. Results. The Herschel/PACS images of $\pi^1$ Gru show an elliptical CSE whose properties agree with those derived from a CO map published in the literature. In addition, an arc east of the star is visible at a distance of $38^{\prime\prime}$ from the primary. This arc is most likely part of an Archimedean spiral caused by an already known G0V companion that is orbiting the primary at a projected distance of 460 au with a period of more than 6200 yr. However, the presence of the elliptical CSE, proper motion variations, and geometric modelling of the VLTI/AMBER observations point towards a third component in the system, with an orbital period shorter than 10 yr, orbiting much closer to the primary than the G0V star.Comment: 13 pages, 11 figures, accepted for publication in Astronomy & Astrophysic

The wonderful complexity of the Mira AB system

We have mapped the CO(3-2) line emission around the Mira AB system at 0.5 resolution using the Atacama Large Millimeter/submillimeter Array (ALMA). The CO map shows amazing complexity. The circumstellar gas has been shaped by different dynamical actors during the evolution of the system and several morphological components can be identified. The companion is marginally resolved in continuum emission and is currently at 0.487$\pm$0.006 separation. In the main line component, centered on the stellar velocity, spiral arcs around Mira A are found. The spiral appears to be relatively flat and oriented in the orbital plane. An accretion wake behind the companion is clearly visible and the projected arc separation is of order 5''. In the blue wing of the line emission, offset from the main line, several large ($\sim$5-10''), opposing arcs are found. We tentatively suggest that this structure is created by the wind of Mira B blowing a bubble in the expanding envelope of Mira A.Comment: Letter accepted in A&