• ### New detections of (sub)millimeter hydrogen radio recombination lines towards high-mass star-forming clumps(1804.06999)

April 19, 2018 astro-ph.GA, astro-ph.SR
Previous radio recombination line (RRL) observations of dust clumps identified in the APEX Telescope Large Area Survey of the Galaxy (ATLASGAL) have led to the detection of a large number of RRLs in the 3mm range. Here, we aim to study their excitation with shorter wavelength (sub)millimeter radio recombination line (submm-RRL) observations. We made observations of submm-RRLs with low principal quantum numbers ($n$ $\leq$ 30) using the APEX 12 m telescope, toward 104 HII regions associated with massive dust clumps from ATLASGAL. The observations covered the H25$\alpha$, H28$\alpha$, and H35$\beta$ transitions. Toward a small subsample the H26$\alpha$, H27$\alpha$, H29$\alpha$, and H30$\alpha$ lines were observed to avoid contamination by molecular lines at adjacent frequencies. We have detected submm-RRLs (signal-to-noise $\geq$ 3 $\sigma$) from compact HII regions embedded within 93 clumps. The submm-RRLs are approximately a factor of two brighter than the mm-RRLs and consistent with optically thin emission in local thermodynamic equilibrium (LTE). The average ratio (0.31) of the measured H35$\beta$/H28$\alpha$ fluxes is close to the LTE value of 0.28. No indication of RRL maser emission has been found. The Lyman photon flux, bolometric, and submm-RRL luminosities toward the submm-RRL detected sources present significant correlations. The trends of dust temperature and the ratio of bolometric luminosity to clump mass, $L_{\rm bol}/M_{\rm clump}$, indicate that the HII regions are related to the most massive and luminous clumps. By estimating the production rate of ionizing photons, $Q$, from the submm-RRL flux, we find that the $Q$(H28$\alpha$) measurements provide estimates of the Lyman continuum photon flux consistent with those determined from 5 GHz radio continuum emission.
• ### The search for high-mass protostars with ALMA revealed up to kilo-parsec scales (SPARKS): I. Indication for a centrifugal barrier in the environment of a single high-mass envelope(1804.06482)

April 17, 2018 astro-ph.GA, astro-ph.SR
The formation of the most massive O-type stars is poorly understood. We present a case study of a young massive clump from the ATLASGAL survey, G328.2551-0.5321. It exhibits a bolometric luminosity of 1.3$\times$10$^4$ L$_{\odot}$ corresponding to a current protostellar mass of $\sim$11 and 16 M$_{\odot}$. We analyze high angular-resolution observations with ALMA at $\sim$0.17" corresponding a physical scale of $\sim$400 au in dust continuum and molecular lines. The dust continuum emission reveals a single high-mass protostellar envelope and shows evidence for a marginally resolved continuum source. We detect a rotational line of CH$_3$OH within its $v_{\rm t}$=1 torsionally excited state revealing two bright peaks of emission spatially offset from the dust continuum peak, and exhibiting a distinct velocity component $\pm$4.5 km s$^{-1}$ offset compared to the source $v_{\rm lsr}$. Local thermodynamic equilibrium analysis suggests N(CH$_3$OH)=1.2$-$2$\times$10$^{19}$ cm$^{-2}$, and kinetic temperatures of 160$-$200 K at the position of these peaks. Their velocity shifts correspond well to the expected Keplerian velocity around a central object with 15M$_{\odot}$ consistent with the mass estimate based on the source's bolometric luminosity. We propose a picture where the CH$_3$OH emission peaks trace the accretion shocks around the centrifugal barrier, pinpointing the interaction region between the collapsing envelope and an accretion disk. Because the HC$_3$N $v_{\rm 7}$=1e ($J$=38-37) line shows compact emission, and a velocity pattern consistent with models of Keplerian rotation, we suggest that this could be a new tracer for compact accretion disks around high-mass protostars. The estimated physical properties of the accretion disk suggest a specific angular momentum several times larger than typically observed towards low-mass protostars.
• ### Weak-lensing mass calibration of the Sunyaev-Zel'dovich effect using APEX-SZ galaxy clusters(1804.03671)

April 10, 2018 astro-ph.CO
The use of galaxy clusters as precision cosmological probes relies on an accurate determination of their masses. However, inferring the relationship between cluster mass and observables from direct observations is difficult and prone to sample selection biases. In this work, we use weak lensing as the best possible proxy for cluster mass to calibrate the Sunyaev-Zel'dovich (SZ) effect measurements from the APEX-SZ experiment. For a well-defined (ROSAT) X-ray complete cluster sample, we calibrate the integrated Comptonization parameter, $Y_{\rm SZ}$, to the weak-lensing derived total cluster mass, $M_{500}$. We employ a novel Bayesian approach to account for the selection effects by jointly fitting both the SZ Comptonization, $Y_{\rm SZ}\text{--}M_{500}$, and the X-ray luminosity, $L_{\rm x}\text{--}M_{500}$, scaling relations. We also account for a possible correlation between the intrinsic (log-normal) scatter of $L_{\rm x}$ and $Y_{\rm SZ}$ at fixed mass. We find the corresponding correlation coefficient to be $r= 0.47_{-0.35}^{+0.24}$, and at the current precision level our constraints on the scaling relations are consistent with previous works. For our APEX-SZ sample, we find that ignoring the covariance between the SZ and X-ray observables biases the normalization of the $Y_{\rm SZ}\text{--}M_{500}$ scaling high by $1\text{--}2\sigma$ and the slope low by $\sim 1\sigma$, even when the SZ effect plays no role in the sample selection. We conclude that for higher-precision data and larger cluster samples, as anticipated from on-going and near-future cluster cosmology experiments, similar biases (due to intrinsic covariances of cluster observables) in the scaling relations will dominate the cosmological error budget if not accounted for correctly.
• ### OH absorption in the first quadrant of the Milky Way as seen by THOR(1803.04794)

March 13, 2018 astro-ph.GA
The hydroxyl radical (OH) is present in the diffuse molecular and partially atomic phases of the interstellar medium (ISM), but its abundance relative to hydrogen is not clear. We aim to evaluate the abundance of OH with respect to molecular hydrogen using OH absorption against cm-continuum sources over the first Galactic quadrant. This OH study is part of the HI/OH/Recombination line survey (THOR). THOR is a Karl G. Jansky Very Large Array large program of atomic, molecular and ionized gas in the range 15{\deg}$\leq$l$\leq$67{\deg} and |b|$\leq$1{\deg}. It is the highest-resolution unbiased OH absorption survey to date towards this region. We combine the derived optical depths with literature 13CO(1-0) and HI observations to determine the OH abundance. We detect absorption in the 1665 and 1667 MHz transitions for continuum sources stronger than $F_{\rm cont}\geq$0.1 Jy/beam. OH absorption is found against $\sim$15% of these continuum sources with increasing fractions for stronger sources. Most of the absorption is associated with Galactic HII regions. We find OH and 13CO gas to have similar kinematic properties. The OH abundance decreases with increasing hydrogen column density. The OH abundance with respect to the total hydrogen nuclei column density (atomic and molecular phase) is in agreement with a constant abundance for $A_V$ < 10-20. Towards the lowest column densities, we find sources that exhibit OH absorption but no 13CO emission, indicating that OH is a well suited tracer of the low column density molecular gas. We present spatially resolved OH absorption towards W43. The unbiased nature of the THOR survey opens a new window onto the gas properties of the ISM. The characterization of the OH abundance over a large range of hydrogen gas column densities contributes to the understanding of OH as a molecular gas tracer and provides a starting point for future investigations.
• ### A revised distance to IRAS 16293$-$2422 from VLBA astrometry of associated water masers(1802.03234)

Feb. 9, 2018 astro-ph.GA, astro-ph.SR
IRAS 16293-2422 is a very well studied young stellar system seen in projection towards the L1689N cloud in the Ophiuchus complex. However, its distance is still uncertain with a range of values from 120 pc to 180 pc. Our goal is to measure the trigonometric parallax of this young star by means of H$_2$O maser emission. We use archival data from 15 epochs of VLBA observations of the 22.2 GHz water maser line. By modeling the displacement on the sky of the H$_2$O maser spots, we derived a trigonometric parallax of $7.1\pm1.3$ mas, corresponding to a distance of $141_{-21}^{+30}$ pc. This new distance is in good agreement with recent values obtained for other magnetically active young stars in the L1689 cloud. We relate the kinematics of these masers with the outflows and the recent ejections powered by source A in the system.
• ### ATLASGAL-selected massive clumps in the inner Galaxy: VI. Kinetic temperature and spatial density measured with formaldehyde(1711.10012)

Nov. 27, 2017 astro-ph.GA, astro-ph.SR
We aim to directly determine the kinetic temperature and spatial density with formaldehyde for the $\sim$100 brightest ATLASGAL-selected clumps at 870 $\mu$m representing various evolutionary stages of high-mass star formation. Ten transitions ($J$ = 3-2 and 4-3) of ortho- and para-H$_2$CO near 211, 218, 225, and 291 GHz were observed with the APEX 12 m telescope. Using non-LTE models with RADEX, we derive the gas kinetic temperature and spatial density using the measured p-H$_2$CO 3$_{21}$-2$_{20}$/3$_{03}$-2$_{02}$, 4$_{22}$-3$_{21}$/4$_{04}$-3$_{03}$, and 4$_{04}$-3$_{03}$/3$_{03}$-2$_{02}$ ratios. The gas kinetic temperatures derived from the p-H$_2$CO 3$_{21}$-2$_{20}$/3$_{03}$-2$_{02}$ and 4$_{22}$-3$_{21}$/4$_{04}$-3$_{03}$ line ratios are high, ranging from 43 to $>$300 K with an unweighted average of 91 $\pm$ 4 K. Deduced $T_{\rm kin}$ values from the $J$ = 3-2 and 4-3 transitions are similar. Spatial densities of the gas derived from the p-H$_2$CO 4$_{04}$-3$_{03}$/3$_{03}$-2$_{02}$ line ratios yield 0.6-8.3 $\times$ 10$^6$ cm$^{-3}$ with an unweighted average of 1.5 ($\pm$0.1) $\times$ 10$^6$ cm$^{-3}$. A comparison of kinetic temperatures derived from p-H$_2$CO, NH$_3$, and the dust emission indicates that p-H$_2$CO traces a distinctly higher temperature than the NH$_3$ (2,2)/(1,1) transitions and the dust, tracing heated gas more directly associated with the star formation process. The H$_2$CO linewidths are found to be correlated with bolometric luminosity and increase with the evolutionary stage of the clumps, which suggests that higher luminosities tend to be associated with a more turbulent molecular medium. It seems that the spatial densities measured with H$_2$CO do not vary significantly with the evolutionary stage of the clumps. However, averaged gas kinetic temperatures derived from H$_2$CO increase with time through the evolution of the clumps.
• ### Galactocentric variation of the gas-to-dust ratio and its relation with metallicity(1710.05721)

Oct. 16, 2017 astro-ph.GA
(Abridged) Context: The assumption of a gas-to-dust mass ratio (\gamma) is a common approach to estimate the basic properties of molecular clouds, such as total mass and column density of molecular hydrogen, from (sub)mm continuum observations of the dust. In the Milky Way a single value is used at all galactocentric radii, independently of the observed metallicity gradients. Both models and extragalactic observations suggest that this quantity increases for decreasing metallicity Z, typical of the outer regions in disks, where fewer heavy elements are available to form dust grains. Aims: We aim to investigate the variation of the gas-to-dust ratio as a function of galactocentric radius and metallicity, to allow a more accurate characterisation of the quantity of molecular gas across the galactic disk, as derived from observations of the dust. Methods: Observations of the optically thin C\$^{18}\$O (2-1) transition were obtained with the APEX telescope for a sample of 23 massive and dense star-forming regions in the far outer Galaxy (galactocentric distance greater than 14 kpc). From the modelling of this line and of the spectral energy distribution of the selected clumps we computed the gas-to-dust ratio and compared it to that of well-studied sources from the ATLASGAL TOP100 sample in the inner galactic disk. Results: The gradient in gas-to-dust ratio is found to be 0.087 dex/kpc (or equivalently \gamma\ proportional to Z\$^{-1.4}\$). The dust-to-metal ratio, decreases with galactocentric radius, which is the most common situation also for external late-type galaxies. This suggests that grain growth dominates over destruction. The predicted gas-to-dust ratio is in excellent agreement with the estimates in Magellanic clouds, for the appropriate value of Z.
• ### Circumstellar ammonia in oxygen-rich evolved stars(1710.01027)

Oct. 3, 2017 astro-ph.GA, astro-ph.SR
The circumstellar ammonia (NH$_3$) chemistry in evolved stars is poorly understood. Previous observations and modelling showed that NH$_3$ abundance in oxygen-rich stars is several orders of magnitude above that predicted by equilibrium chemistry. In this article, we characterise the spatial distribution and excitation of NH$_3$ in the O-rich circumstellar envelopes (CSEs) of four diverse targets: IK Tau, VY CMa, OH 231.8+4.2, and IRC +10420 with multi-wavelength observations. We observed the 1.3-cm inversion line emission with the Very Large Array (VLA) and submillimetre rotational line emission with the Heterodyne Instrument for the Far-Infrared (HIFI) aboard Herschel from all four targets. For IK Tau and VY CMa, we observed the rovibrational absorption lines in the $\nu_2$ band near 10.5 $\mu$m with the Texas Echelon Cross Echelle Spectrograph (TEXES) at the NASA Infrared Telescope Facility (IRTF). We also attempted to search for the rotational transition within the $v_2=1$ state near 2 mm with the IRAM 30m Telescope towards IK Tau. Non-LTE radiative transfer modelling, including radiative pumping to the vibrational state, was carried out to derive the radial distribution of NH$_3$ in these CSEs. Our modelling shows that the NH$_3$ abundance relative to molecular hydrogen is generally of the order of $10^{-7}$, which is a few times lower than previous estimates that were made without considering radiative pumping and is at least 10 times higher than that in the C-rich CSE of IRC +10216. Incidentally, we also derived a new period of IK Tau from its $V$-band light curve. NH$_3$ is again detected in very high abundance in O-rich CSEs. Its emission mainly arises from localised spatial-kinematic structures that are probably denser than the ambient gas. Circumstellar shocks in the accelerated wind may contribute to the production of NH$_3$. (Abridged abstract)
• ### Kinetic temperature of massive star-forming molecular clumps measured with formaldehyde. III. The Orion Molecular Cloud 1(1709.07694)

Sept. 22, 2017 astro-ph.GA, astro-ph.SR
We mapped the kinetic temperature structure of the Orion molecular cloud 1 with para-H2CO(303-202, 322-221, and 321-220) using the APEX 12m telescope. This is compared with the temperatures derived from the ratio of the NH3(2,2)/(1,1) inversion lines and the dust emission. Using the RADEX non-LTE model, we derive the gas kinetic temperature modeling the measured averaged line ratios of para-H2CO 322-221/303-202 and 321-220/303-202. The gas kinetic temperatures derived from the para-H2CO line ratios are warm, ranging from 30 to >200 K with an average of 62 K at a spatial density of 10$^5$ cm$^{-3}$. These temperatures are higher than those obtained from NH3(2,2)/(1,1) and CH3CCH(6-5) in the OMC-1 region. The gas kinetic temperatures derived from para-H2CO agree with those obtained from warm dust components measured in the mid infrared (MIR), which indicates that the para-H2CO(3-2) ratios trace dense and warm gas. The cold dust components measured in the far infrared (FIR) are consistent with those measured with NH3(2,2)/(1,1) and the CH3CCH(6-5) line series. With dust at MIR wavelengths and para-H2CO(3-2) on one side and dust at FIR wavelengths, NH3(2,2)/(1,1), and CH3CCH(6-5) on the other, dust and gas temperatures appear to be equivalent in the dense gas of the OMC-1 region, but provide a bimodal distribution, one more directly related to star formation than the other. The non-thermal velocity dispersions of para-H2CO are positively correlated with the gas kinetic temperatures in regions of strong non-thermal motion (Mach number >2.5) of the OMC-1, implying that the higher temperature traced by para-H2CO is related to turbulence on a 0.06 pc scale. Combining the temperature measurements with para-H2CO and NH3(2,2)/(1,1) line ratios, we find direct evidence for the dense gas along the northern part of the OMC-1 10 km s$^{-1}$ filament heated by radiation from the central Orion nebula.
• ### Testing the variability of the proton-to-electron mass ratio from observations of methanol in the dark cloud core L1498(1709.03103)

Sept. 10, 2017 astro-ph.CO
The dependence of the proton-to-electron mass ratio, mu, on the local matter density was investigated using methanol emission in the dense dark cloud core L1498. Towards two different positions in L1498, five methanol transitions were detected and an extra line was tentatively detected at a lower confidence level in one of the positions. The observed centroid frequencies were then compared with their rest frame frequencies derived from least-squares fitting to a large data set. Systematic effects, as the underlying methanol hyperfine structure and the Doppler tracking of the telescope, were investigated and their effects were included in the total error budget. The comparison between the observations and the rest frame frequencies constrains potential mu variation at the level of Dmu/mu < 6 x 10^(-8), at a 3 sigma confidence level. For the dark cloud we determine a total CH3OH (A+E) beam averaged column density of 3-4 x 10^(12) cm(-2) (within roughly a factor of two), an E- to A-type methanol column density ratio of N(A-CH3OH)/N(E-CH3OH) = 1.00 +/- 0.15, a density of n(H2) = 3 x 10^5 cm^(-3) (again within a factor of two), and a kinetic temperature of Tkin = 6 +/- 1 K. In a kinetic model including the line intensities observed for the methanol lines, the n(H2) density is higher and the temperature is lower than that derived in previous studies based on different molecular species; the intensity of the 1_0 --> 1_-1 E line strength is not well reproduced.
• ### ATLASGAL - Ammonia observations towards the southern Galactic Plane(1708.07839)

Aug. 25, 2017 astro-ph.GA, astro-ph.SR
The initial conditions of molecular clumps in which high-mass stars form are poorly understood. In particular, a more detailed study of the earliest evolutionary phases is needed. The APEX Telescope Large Area Survey of the whole inner Galactic disk at 870 micron, ATLASGAL, has been conducted to discover high-mass star-forming regions at different evolutionary phases. Using the Parkes telescope, we observed the NH3 (1,1) to (3,3) inversion transitions towards 354 ATLASGAL clumps in the fourth quadrant. For a subsample of 289 sources, the N2H+ (1-0) line was measured with the Mopra telescope. We measured a median NH3(1,1) line width of about 2 km/s and rotational temperatures from 12 to 28 K with a mean of 18 K. For a subsample with detected NH3 (2,2) hyperfine components, we found that the commonly used method to compute the (2,2) optical depth from the (1,1) optical depth and the (2,2) to (1,1) main beam brightness temperature ratio leads to an underestimation of the rotational temperature and column density. A larger median virial parameter of about 1 is determined using the broader N2H+ line width than is estimated from the NH3 line width of about 0.5 with a general trend of a decreasing virial parameter with increasing gas mass. We found a warmer surrounding of ATLASGAL clumps than the surrounding of low-mass cores and smaller velocity dispersions in low-mass than high-mass star-forming regions. The NH3 (1,1) inversion transition of 49% of the sources shows hyperfine structure anomalies. The intensity ratio of the outer hyperfine structure lines with a median of 1.27+/-0.03 and a standard deviation of 0.45 is significantly higher than 1, while the intensity ratios of the inner satellites with a median of 0.9+/-0.02 and standard deviation of 0.3 and the sum of the inner and outer hyperfine components with a median of 1.06+/-0.02 and standard deviation of 0.37 are closer to 1.
• ### Complex organic molecules in diffuse clouds along the line of sight to Sgr B2(1708.07292)

Aug. 24, 2017 astro-ph.GA
Up to now, mostly relatively simple molecules have been detected in interstellar diffuse molecular clouds in our galaxy, but more complex species have been reported in the diffuse/translucent medium of a z = 0.89 spiral galaxy. We aim at searching for complex organic molecules (COMs) in diffuse molecular clouds along the line of sight to Sgr B2(N), taking advantage of the high sensitivity and angular resolution of the Atacama Large Millimeter/submillimeter Array (ALMA). We use data acquired as part of the EMoCA survey performed with ALMA. To analyse the absorption features of the molecules detected towards the ultracompact HII region K4 in Sgr B2(N), we calculate synthetic spectra for these molecules and fit their column densities, line widths, centroid velocities, and excitation temperatures. We report the detection of CH3OH, CH3CN, CH3CHO, HC3N, and NH2CHO in Galactic center (GC) diffuse clouds and CH3OH and CH3CN in a diffuse cloud in the Scutum arm. The chemical composition of one of the diffuse GC clouds is found to be similar to the one of the diffuse/translucent medium of the z=0.89 spiral galaxy. The chemical processes leading to chemical complexity in the diffuse molecular ISM appear to have remained similar since z=0.89. As proposed in previous studies, the presence of COMs in diffuse molecular clouds may result from a cyclical interstellar process of cloud contraction and expansion between diffuse and dense states.
• ### Organic molecules, ions, and rare isotopologues in the remnant of the stellar-merger candidate, CK Vulpeculae (Nova 1670)(1708.02261)

Aug. 18, 2017 astro-ph.SR
CK Vul is a star whose outburst was observed in 1670-72. A stellar-merger event was proposed to explain its ancient eruption. Aims: We aim to investigate the composition of the molecular gas recently discovered in the remnant of CK Vul. Methods: We observed millimeter and submillimeter-wave spectra of CK Vul using the IRAM 30m and APEX telescopes. Radiative-transfer modeling of the observed molecular features was performed to yield isotopic ratios for various elements. Results: The spectra of CK Vul reveal a very rich molecular environment of low excitation ($T_{\rm ex} \lesssim$12 K). Atomic carbon and twenty seven different molecules, including two ions, were identified. They range from simple diatomic to complex polyatomic species of up to 7 atoms large. The chemical composition of the molecular gas is indicative of carbon and nitrogen-driven chemistry but oxides are also present. Additionally, the abundance of F may be enhanced. The spectra are rich in isotopologues that are very rare in most known sources. All stable isotopes of C, N, O, Si, and S are observed and their isotopic ratios are derived. Conclusions: The composition of the remnant's molecular gas is most peculiar and gives rise to a very unique millimeter and submillimeter spectrum. The observation of ions and complex molecules suggests the presence of a photoionizing source but its nature (a central star or shocks) remains unknown. The elemental and isotopic composition of the gas cannot be easily reconciled with standard nucleosynthesis but processing in hot CNO cycles and partial He burning can explain most of the chemical peculiarities. The isotopic ratios of CK Vul are remarkably close to those of presolar nova grains but the link of Nova 1670 to objects responsible for these grains is unclear.
• ### HIFI Spectroscopy of ${\rm H_2O}$ submm Lines in Nuclei of Actively Star Forming Galaxies(1707.04914)

July 16, 2017 astro-ph.GA
We present a systematic survey of multiple velocity-resolved H$_2$O spectra using Herschel/HIFI towards nine nearby actively star forming galaxies. The ground-state and low-excitation lines (E$_{\rm up}\,\le 130\,{\rm K}$) show profiles with emission and absorption blended together, while absorption-free medium-excitation lines ($130\,{\rm K}\, \le\, E_{\rm up}\,\le\,350\,{\rm K}$) typically display line shapes similar to CO. We analyze the HIFI observation together with archival SPIRE/PACS H$_2$O data using a state-of-the-art 3D radiative transfer code which includes the interaction between continuum and line emission. The water excitation models are combined with information on the dust- and CO spectral line energy distribution to determine the physical structure of the interstellar medium (ISM). We identify two ISM components that are common to all galaxies: A warm ($T_{\rm dust}\,\sim\,40-70\,{\rm K}$), dense ($n({\rm H})\,\sim\,10^5-10^6\,{\rm cm^{-3}}$) phase which dominates the emission of medium-excitation H$_2$O lines. This gas phase also dominates the FIR emission and the CO intensities for $J_{\rm up} > 8$. In addition a cold ($T_{\rm dust}\,\sim\,20-30\,{\rm K}$), dense ($n({\rm H})\sim\,10^4- 10^5\,{\rm cm^{-3}}$) more extended phase is present. It outputs the emission in the low-excitation H$_2$O lines and typically also produces the prominent line absorption features. For the two ULIRGs in our sample (Arp 220 and Mrk 231) an even hotter and more compact (R$_s\,\le\,100$ pc) region is present which is possibly linked to AGN activity. We find that collisions dominate the water excitation in the cold gas and for lines with $E_{\rm up}\le300\,{\rm K}$ and $E_{\rm up}\le800\,{\rm K}$ in the warm and hot component, respectively. Higher energy levels are mainly excited by IR pumping.
• ### Techniques for Accurate Parallax Measurements for 6.7-GHz Methanol Masers(1706.03128)

June 9, 2017 astro-ph.IM
The BeSSeL Survey is mapping the spiral structure of the Milky Way by measuring trigonometric parallaxes of hundreds of maser sources associated with high-mass star formation. While parallax techniques for water masers at high frequency (22 GHz) have been well documented, recent observations of methanol masers at lower frequency (6.7 GHz) have revealed astrometric issues associated with signal propagation through the ionosphere that could significantly limit parallax accuracy. These problems displayed as a "parallax gradient" on the sky when measured against different background quasars. We present an analysis method in which we generate position data relative to an "artificial quasar" at the target maser position at each epoch. Fitting parallax to these data can significantly mitigate the problems and improve parallax accuracy.
• ### SiS in the circumstellar envelope of IRC +10126: maser and quasi-thermal emission(1706.02446)

June 8, 2017 astro-ph.GA, astro-ph.SR
We present new Effelsberg-100 m, ATCA, and VLA observations of rotational SiS transitions in the circumstellar envelope (CSE) of IRC +10216. Thanks to the high angular resolution achieved by the ATCA observations, we unambiguously confirm that the molecule's J=1-0 transition exhibits maser action in this CSE, as first suggested more than thirty years ago. The maser emission's radial velocity peaking at a local standard of rest velocity of -39.862$\pm$0.065 km/s indicates that it arises from an almost fully accelerated shell. Monitoring observations show time variability of the SiS (1-0) maser. The two lowest-$J$ SiS quasi-thermal emission lines trace a much more extended emitting region than previous high-J SiS observations. Their distributions show that the SiS quasi-thermal emission consists of two components: one is very compact (radius<1.5", corresponding to <3$\times 10^{15}$ cm), and the other extends out to a radius >11". An incomplete shell-like structure is found in the north-east, which is indicative of existing SiS shells. Clumpy structures are also revealed in this CSE. The gain of the SiS (1-0) maser (optical depths of about -5 at the blue-shifted side and, assuming inversion throughout the entire line's velocity range, about -2 at the red-shifted side) suggests that it is unsaturated. The SiS (1-0) maser can be explained in terms of ro-vibrational excitation caused by infrared pumping, and we propose that infrared continuum emission is the main pumping source.
• ### Radio continuum of galaxies with H$_{2}$O megamaser disks: 33 GHz VLA data(1706.02699)

June 8, 2017 astro-ph.GA
We investigate the nuclear environment of galaxies with observed 22 GHz water megamaser in their subparsec edge-on accretion disks, using 33 GHz (9mm) radio continuum data from VLA, with a resolution of ~ 0.2-0.5 arcsecs, and relate the maser and host galaxy properties to those of its radio continuum emission. Eighty-seven percent (21 out of 24) galaxies in our sample show 33 GHz radio continuum emission at levels of 4.5-240 $\sigma$. Five sources show extended emission, including one source with two main components and one with three main components. The remaining detected 16 sources exhibit compact cores within the sensitivity limits. Little evidence is found for extended jets (>300 pc) in most sources. Either they do not exist, or our chosen frequency of 33 GHz is too high for a detection of these supposedly steep spectrum features. In only one source among those with known maser disk orientation, NGC4388, we found an extended jet-like feature that appears to be oriented perpendicular to the water megamaser disk. Smaller 100-300 pc sized jets might also be present, as is suggested by the beam-deconvolved morphology of our sources. Whenever possible, central positions with accuracies of 20-280 mas are provided. A correlation analysis shows that the 33 GHz luminosity weakly correlates with the infrared luminosity. The 33 GHz luminosity is anticorrelated with the circular velocity of the galaxy. The black hole masses show stronger correlations with water maser luminosity than with 1.4 GHz, 33 GHz, or hard X-ray luminosities. Furthermore, the inner radii of the disks show stronger correlations with 1.4 GHz, 33 GHz, and hard X-ray luminosities than their outer radii, suggesting that the outer radii may be affected by disk warping, star formation, or peculiar density distributions.
• ### ISM properties of a Massive Dusty Star-Forming Galaxy discovered at z ~ 7(1705.07912)

May 25, 2017 astro-ph.GA
We report the discovery and constrain the physical conditions of the interstellar medium of the highest-redshift millimeter-selected dusty star-forming galaxy (DSFG) to date, SPT-S J031132-5823.4 (hereafter SPT0311-58), at $z=6.900 +/- 0.002$. SPT0311-58 was discovered via its 1.4mm thermal dust continuum emission in the South Pole Telescope (SPT)-SZ survey. The spectroscopic redshift was determined through an ALMA 3mm frequency scan that detected CO(6-5), CO(7-6) and [CI](2-1), and subsequently confirmed by detections of CO(3-2) with ATCA and [CII] with APEX. We constrain the properties of the ISM in SPT0311-58 with a radiative transfer analysis of the dust continuum photometry and the CO and [CI] line emission. This allows us to determine the gas content without ad hoc assumptions about gas mass scaling factors. SPT0311-58 is extremely massive, with an intrinsic gas mass of $M_{\rm gas} = 3.3 \pm 1.9 \times10^{11}\,M_{\odot}$. Its large mass and intense star formation is very rare for a source well into the Epoch of Reionization.
• ### Deep submillimeter and radio observations in the SSA22 field. I. Powering sources and Ly{\alpha} escape fraction of Ly{\alpha} blobs(1704.05101)

April 17, 2017 astro-ph.GA
We study the heating mechansims and Ly{\alpha} escape fractions of 35 Ly{\alpha} blobs (LABs) at z = 3.1 in the SSA22 field, using the submillimeter (submm) data from the ALMA and the SCUBA-2 on the JCMT, and the radio data from the VLA. With ALMA and JCMT/SCUBA-2, we have identified dust continuum sources in 11 out of 35 LABs in the SSA22 field, with star formation rates (SFRs) of above 100Msun/yr. Using VLA, we have detected likely radio counterparts in 9 out of 29 LABs in the SSA22 field. We find that a radio excess in the submm/radio detected LABs is common, suggesting that active galactic nuclei (AGN) may often be present in these systems. Most radio sources without X-ray counterparts are located in the centers of the LABs. However, it is surprising that all X-ray counterparts avoid the central regions of their LABs. Possibly these X-ray sources are associated with merging systems, which re-distribute the interstellar meidium (ISM) and circumgalactic medium (CGM), leading to highly inhomogeneous gas. Alternatively, AGN and its feedback may also be responsible for inhomogeneous Ly{\alpha} morphology in these LABs. We obtain escape fractions of Ly{\alpha} photons for the LABs based on the submm and Ly{\alpha} data, and find that the Ly{\alpha} escape fractions in the LABs are significantly lower than those found for Ly{\alpha} emitters (LAEs) in other studies. We suspect that this large difference is due to the SCUBA-2 detection threshold, the high dust attenuation indicated by the large SFRs in our sample, as well as the dense large-scale environment in SSA22.
• ### Exploring molecular complexity with ALMA (EMoCA): Detection of three new hot cores in Sagittarius B2(N)(1703.09544)

March 28, 2017 astro-ph.GA
The SgrB2 molecular cloud contains several sites forming high-mass stars. SgrB2(N) is one of its main centers of activity. It hosts several compact and UCHII regions, as well as two known hot molecular cores (SgrB2(N1) and SgrB2(N2)), where complex organic molecules are detected. Our goal is to use the high sensitivity of ALMA to characterize the hot core population in SgrB2(N) and shed a new light on the star formation process. We use a complete 3 mm spectral line survey conducted with ALMA to search for faint hot cores in SgrB2(N). We report the discovery of three new hot cores that we call SgrB2(N3), SgrB2(N4), and SgrB2(N5). The three sources are associated with class II methanol masers, well known tracers of high-mass star formation, and SgrB2(N5) also with a UCHII region. The chemical composition of the sources and the column densities are derived by modelling the whole spectra under the assumption of LTE. The H2 column densities are computed from ALMA and SMA continuum emission maps. The H2 column densities of these new hot cores are found to be 16 up to 36 times lower than the one of the main hot core Sgr B2(N1). Their spectra have spectral line densities of 11 up to 31 emission lines per GHz, assigned to 22-25 molecules. We derive rotational temperatures around 140-180 K for the three new hot cores and mean source sizes of 0.4 for SgrB2(N3) and 1.0 for SgrB2(N4) and SgrB2(N5). SgrB2(N3) and SgrB2(N5) show high velocity wing emission in typical outflow tracers, with a bipolar morphology in their integrated intensity maps suggesting the presence of an outflow, like in SgrB2(N1). The associations of the hot cores with class II methanol masers, outflows, and/or UCHII regions tentatively suggest the following age sequence: SgrB2(N4), SgrB2(N3), SgrB2(N5), SgrB2(N1). The status of SgrB2(N2) is unclear. It may contain two distinct sources, a UCHII region and a very young hot core.
• ### ALMA survey of massive cluster progenitors from ATLASGAL. Limited fragmentation at the early evolutionary stage of massive clumps(1703.03273)

March 9, 2017 astro-ph.GA
The early evolution of massive cluster progenitors is poorly understood. We investigate the fragmentation properties from 0.3 pc to 0.06 pc scales of a homogenous sample of infrared-quiet massive clumps within 4.5 kpc selected from the ATLASGAL survey. Using the ALMA 7m array we detect compact dust continuum emission towards all targets, and find that fragmentation, at these scales, is limited. The mass distribution of the fragments uncovers a large fraction of cores above 40 $M_\odot$, corresponding to massive dense cores (MDCs) with masses up to ~400 $M_\odot$. 77 % of the clumps contain at most 3 MDCs per clump, and we also reveal single clumps/MDCs. The most massive cores are formed within the more massive clumps, and a high concentration of mass on small scales reveals a high core formation efficiency. The mass of MDCs highly exceeds the local thermal Jeans-mass, and observational evidence is lacking for a sufficiently high level of turbulence or strong enough magnetic fields to keep the most massive MDCs in equilibrium. If already collapsing, the observed fragmentation properties with a high core formation efficiency are consistent with the collapse setting in at parsec scales.
• ### Hydrogen in diffuse molecular clouds in the Milky Way: Atomic column densities and molecular fraction along prominent lines of sight(1702.05956)

Feb. 20, 2017 astro-ph.GA
Recent submillimeter and far-infrared wavelength observations of absorption in the rotational ground-state lines of various simple molecules against distant Galactic continuum sources have opened the possibility of studying the chemistry of diffuse molecular clouds throughout the Milky Way. In order to calculate abundances, the column densities of molecular and atomic hydrogen, HI, must be known. We aim at determining the atomic hydrogen column densities for diffuse clouds located on the sight lines toward a sample of prominent high-mass star-forming regions that were intensely studied with the HIFI instrument onboard Herschel. Based on Jansky Very Large Array data, we employ the 21 cm HI absorption-line technique to construct profiles of the HI opacity versus radial velocity toward our target sources. These profiles are combined with lower resolution archival data of extended HI emission to calculate the HI column densities of the individual clouds along the sight lines. We employ Bayesian inference to estimate the uncertainties of the derived quantities. Our study delivers reliable estimates of the atomic hydrogen column density for a large number of diffuse molecular clouds at various Galactocentric distances. Together with column densities of molecular hydrogen derived from its surrogates observed with HIFI, the measurements can be used to characterize the clouds and investigate the dependence of their chemistry on the molecular fraction, for example.
• ### ATLASGAL-selected massive clumps in the inner Galaxy IV. Millimeter hydrogen recombination lines from associated HII regions(1702.02062)

Feb. 7, 2017 astro-ph.GA, astro-ph.SR
Observations of millimeter wavelength radio recombination lines (mm-RRLs) are used to search for HII regions in an unbiased way that is complementary to many of the more traditional methods previously used. We carried out targeted mm-RRL observations (39$\leq$ principal quantum number ($n$) $\leq$65 and $\Delta n$ = 1, 2, 3 and 4, i.e. H$n\alpha$, H$n\beta$, H$n\gamma$ and H$n\delta$) using the IRAM 30m and Mopra 22m telescopes. We observed 976 compact dust clumps selected from a catalog of $\sim$10,000 ATLASGAL sources. We detected H$n\alpha$ mm-RRL emission toward 178 clumps; H$n\beta$, H$n\gamma$ and H$n\delta$ were also detected toward 65, 23, and 22 clumps, respectively. This is the largest sample of mm-RRLs detections published to date. Comparing the positions of these clumps with radio continuum surveys we identified compact radio counterparts for 134 clumps. The nature of the other 44 detections is unclear, but 8 detections are thought to be potentially new HII regions. Broad linewidths are seen toward nine clumps (linewidth > 40 km s$^{-1}$) revealing significant turbulent motions within the ionized gas. We find that the linewidth of the H$^{13}$CO$^{+}$ (1$-$0) emission is significantly wider than those without mm-RRL detection. We also find a correlation between the integrated fluxes of the mm-RRLs and the 6 cm continuum flux densities of their radio counterparts. We find that the mm-RRL emission is associated with HII regions with $n_{e}<10^{5}$ cm$^{-3}$ and HII region diameter >0.03 pc. We detected mm-RRLs toward 178 clumps and identified eight new HII region candidates. The broad mm-RRL from nine clumps may indicate that they arise in very young hyper-compact HII regions. The mm-RRLs trace the radio continuum sources detected by high-resolution observations and their line parameters show associations with the embedded radio sources and their parental molecular clumps.
• ### Exploring molecular complexity with ALMA (EMoCA): Simulations of branched carbon-chain chemistry in Sgr B2(N)(1701.07160)

Jan. 25, 2017 astro-ph.GA
Using mm-wavelength data from ALMA, the EMoCA spectral line survey revealed the presence of both the straight-chain and branched forms of propyl cyanide (C$_3$H$_7$CN) toward the Galactic Center star-forming source Sgr B2(N2). This was the first interstellar detection of a branched aliphatic molecule. Through computational methods, we seek to explain the observed $i$:$n$ ratio for propyl cyanide, and to predict the abundances of the four different forms of the homologous nitrile, butyl cyanide (C$_4$H$_9$CN). We also investigate whether other molecules will show a similar degree of branching, by modeling alkanes up to pentane (C$_5$H$_{12}$). We use a coupled three-phase chemical kinetics model to simulate the chemistry of Sgr B2(N2), using an updated chemical network that includes grain-surface/ice-mantle formation routes for branched nitriles and alkanes. We use the EMoCA survey data to search for the straight-chain form of butyl cyanide toward Sgr B2(N2). The observed $i$:$n$ ratio for propyl cyanide is reproduced by the models. Butyl cyanide is predicted to show similar abundances to propyl cyanide, and to exhibit strong branching, with the $sec$ form clearly dominant over all others. The addition of CN to acetylene and ethene is found to be important to the production of vinyl, ethyl, propyl, and butyl cyanide. We report a non-detection of $n$-C$_4$H$_9$CN toward Sgr B2(N2), with an abundance at least 1.7 times lower than that of $n$-C$_3$H$_7$CN. This value is within the range predicted by the chemical models. The models indicate that the degree of branching rises with increasing molecular size. The efficiency of CN addition to unsaturated hydrocarbons boosts the abundances of nitriles in the model, and enhances the ratio of straight-to-branched molecule production. The predicted abundance of $s$-C$_4$H$_9$CN makes it a good candidate for future detection toward Sgr B2(N2).
• ### Rotational spectroscopy, tentative interstellar detection, and chemical modelling of N-methylformamide(1701.04640)

Jan. 17, 2017 astro-ph.GA
N-methylformamide, CH3NHCHO, may be an important molecule for interstellar pre-biotic chemistry because it contains a peptide bond. The rotational spectrum of the most stable trans conformer of CH3NHCHO is complicated by strong torsion-rotation interaction due to the low barrier of the methyl torsion. We use two absorption spectrometers in Kharkiv and Lille to measure the rotational spectra over 45--630 GHz. The analysis is carried out using the Rho-axis method and the RAM36 code. We search for N-methylformamide toward the hot molecular core Sgr B2(N2) using a spectral line survey carried out with ALMA. The astronomical results are put into a broader astrochemical context with the help of a gas-grain chemical kinetics model. The laboratory data set for the trans conformer of CH3NHCHO consists of 9469 line frequencies with J <= 62, including the first assignment of the rotational spectra of the first and second excited torsional states. All these lines are fitted within experimental accuracy. We report the tentative detection of CH3NHCHO towards Sgr B2(N2). We find CH3NHCHO to be more than one order of magnitude less abundant than NH2CHO, a factor of two less abundant than CH3NCO, but only slightly less abundant than CH3CONH2. The chemical models indicate that the efficient formation of HNCO via NH + CO on grains is a necessary step in the achievement of the observed gas-phase abundance of CH3NCO. Production of CH3NHCHO may plausibly occur on grains either through the direct addition of functional-group radicals or through the hydrogenation of CH3NCO. Provided the detection of CH3NHCHO is confirmed, the only slight underabundance of this molecule compared to its more stable structural isomer acetamide and the sensitivity of the model abundances to the chemical kinetics parameters suggest that the formation of these two molecules is controlled by kinetics rather than thermal equilibrium.