Initial phases of massive star formation in high infrared extinction clouds. II. Infall and onset of star formation

The onset of massive star formation is not well understood because of observational and theoretical difficulties. To find the dense and cold clumps where massive star formation can take place, we compiled a sample of high infrared extinction clouds, which were observed previously by us in the 1.2 mm continuum emission and ammonia. We try to understand the star-formation stages of the clumps in these high extinction clouds by studying the infall and outflow properties, the presence of a young stellar object (YSO), and the level of the CO depletion through a molecular line survey with the IRAM 30m and APEX 12m telescopes. Moreover, we want to know if the cloud morphology, quantified through the column density contrast between the clump and the clouds, has an impact on the star formation occurring inside it. We find that the HCO+(1-0) line is the most sensitive for detecting infalling motions. SiO, an outflow tracer, was mostly detected toward sources with infall, indicating that infall is accompanied by collimated outflows. The presence of YSOs within a clump depends mostly on its column density; no signs of YSOs were found below 4E22 cm-2. Star formation is on the verge of beginning in clouds that have a low column density contrast; infall is not yet present in the majority of the clumps. The first signs of ongoing star formation are broadly observed in clouds where the column density contrast between the clump and the cloud is higher than two; most clumps show infall and outflow. Finally, the most evolved clumps are in clouds that have a column density contrast higher than three; almost all clumps have a YSO, and in many clumps, the infall has already halted. Hence, the cloud morphology, based on the column density contrast between the cloud and the clumps, seems to have a direct connection with the evolutionary stage of the objects forming inside

A 10-M⊙M_{\odot} YSO with a Keplerian disk and a nonthermal radio jet

We previously observed the star-forming region G16.59$-$0.05 through interferometric observations of both thermal and maser lines, and identified a high-mass young stellar object (YSO) which is surrounded by an accretion disk and drives a nonthermal radio jet. We performed high-angular-resolution (beam FWHM ~0.15") 1.2-mm continuum and line observations towards G16.59$-$0.05 with the Atacama Large Millimeter Array (ALMA). The main dust clump, with size ~10$^4$ au, is resolved into four relatively compact (diameter ~2000 au) millimeter (mm) sources. The source harboring the high-mass YSO is the most prominent in molecular emission. By fitting the emission profiles of several unblended and optically thin transitions of CH$_3$OCH$_3$ and CH$_3$OH, we derived gas temperatures inside the mm-sources in the range 42--131 K, and calculated masses of 1--5 $M_{\odot}$. A well-defined Local Standard of Rest velocity (Vlsr) gradient is detected in most of the high-density molecular tracers at the position of the high-mass YSO, pinpointed by compact 22-GHz free-free emission. This gradient is oriented along a direction forming a large (~70 degree) angle with the radio jet, traced by elongated 13-GHz continuum emission. The butterfly-like shapes of the P-V plots and the linear pattern of the emission peaks of the molecular lines at high velocity confirm that this Vlsr gradient is due to rotation of the gas in the disk surrounding the high-mass YSO. The disk radius is ~500 au, and the Vlsr distribution along the major axis of the disk is well reproduced by a Keplerian profile around a central mass of 10$\pm$2 $M_{\odot}$. The position of the YSO is offset by >~ 0.1" from the axis of the radio jet and the dust emission peak. To explain this displacement we argue that the high-mass YSO could have moved from the center of the parental mm source owing to dynamical interaction with one or more companions.Comment: 16 pages, 12 figures, accepted by Astronomy & Astrophysics, Main Journa

Trigonometric Parallaxes of High Mass Star Forming Regions: the Structure and Kinematics of the Milky Way

Over 100 trigonometric parallaxes and proper motions for masers associated with young, high-mass stars have been measured with the BeSSeL Survey, a VLBA key science project, the EVN, and the Japanese VERA project. These measurements provide strong evidence for the existence of spiral arms in the Milky Way, accurately locating many arm segments and yielding spiral pitch angles ranging from 7 to 20 degrees. The widths of spiral arms increase with distance from the Galactic center. Fitting axially symmetric models of the Milky Way with the 3-D position and velocity information and conservative priors for the solar and average source peculiar motions, we estimate the distance to the Galactic center, Ro, to be 8.34 +/- 0.16 kpc, a circular rotation speed at the Sun, To, to be 240 +/- 8 km/s, and a rotation curve that is nearly flat (a slope of -0.2 +/- 0.4 km/s/kpc) between Galactocentric radii of 5 and 16 kpc. Assuming a "universal" spiral galaxy form for the rotation curve, we estimate the thin disk scale length to be 2.44 +/- 0.16 kpc. The parameters Ro and To are not highly correlated and are relatively insensitive to different forms of the rotation curve. Adopting a theoretically motivated prior that high-mass star forming regions are in nearly circular Galactic orbits, we estimate a global solar motion component in the direction of Galactic rotation, Vsun = 14.6 +/- 5.0 km/s. While To and Vsun are significantly correlated, the sum of these parameters is well constrained, To + Vsun = 255.2 +/- 5.1 km/s, as is the angular speed of the Sun in its orbit about the Galactic center, (To + Vsun)/Ro = 30.57 +/- 0.43 km/s/kpc. These parameters improve the accuracy of estimates of the accelerations of the Sun and the Hulse-Taylor binary pulsar in their Galactic orbits, significantly reducing the uncertainty in tests of gravitational radiation predicted by general relativity.Comment: 38 pages, 6 tables, 6 figures; v2 fixed typos and updated pulsar section; v3 replaced fig 2 (wrong file

Techniques for Accurate Parallax Measurements for 6.7-GHz Methanol Masers

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

Italian Science Case for ALMA Band 2+3

The Premiale Project "Science and Technology in Italy for the upgraded ALMA Observatory - iALMA" has the goal of strengthening the scientific, technological and industrial Italian contribution to the Atacama Large Millimeter/submillimeter Array (ALMA), the largest ground based international infrastructure for the study of the Universe in the microwave. One of the main objectives of the Science Working Group (SWG) inside iALMA, the Work Package 1, is to develop the Italian contribution to the Science Case for the ALMA Band 2 or Band 2+3 receiver. ALMA Band 2 receiver spans from ~67 GHz (bounded by an opaque line complex of ozone lines) up to 90 GHz which overlaps with the lower frequency end of ALMA Band 3. Receiver technology has advanced since the original definition of the ALMA frequency bands. It is now feasible to produce a single receiver which could cover the whole frequency range from 67 GHz to 116 GHz, encompassing Band 2 and Band 3 in a single receiver cartridge, a so called Band 2+3 system. In addition, upgrades of the ALMA system are now foreseen that should double the bandwidth to 16 GHz. The science drivers discussed below therefore also discuss the advantages of these two enhancements over the originally foreseen Band 2 system.Comment: 43 pages, 21 figure

Clustered star formation and outflows in AFGL 2591

We report on a detailed study of the water maser kinematics and radio continuum emission toward the most massive and young object in the star-forming region AFGL 2591. Our analysis shows at least two spatial scales of multiple star formation, one projected across 0.1 pc on the sky and another one at about 2000 AU from a ZAMS star of about 38 Solar masses. This young stellar object drives a powerful jet- and wind-driven outflow system with the water masers associated to the outflow walls, previously detected as a limb-brightened cavity in the NIR band. At about 1300 AU to the north of this object a younger protostar drives two bow shocks, outlined by arc-like water maser emission, at 200 AU either side of the source. We have traced the velocity profile of the gas that expands along these arc-like maser structures and compared it with the jet-driven outflow model. This analysis suggests that the ambient medium around the northern protostar is swept up by a jet-driven shock (>66 km/s) and perhaps a lower-velocity (~10 km/s) wind with an opening angle of about 20 degrees from the jet axis.Comment: 21 pages, 4 figures, 2 tables, accepted by The Astrophysical Journa

What determines the density structure of molecular clouds? A case study of Orion B with <i>Herschel</i>

A key parameter to the description of all star formation processes is the density structure of the gas. In this Letter, we make use of probability distribution functions (PDFs) of Herschel column density maps of Orion B, Aquila, and Polaris, obtained with the Herschel Gould Belt survey (HGBS). We aim to understand which physical processes influence the PDF shape, and with which signatures. The PDFs of Orion B (Aquila) show a lognormal distribution for low column densities until AV ~ 3 (6), and a power-law tail for high column densities, consistent with a ρα r-2 profile for the equivalent spherical density distribution. The PDF of Orion B is broadened by external compression due to the nearby OB stellar aggregates. The PDF of a quiescent subregion of the non-star-forming Polaris cloud is nearly lognormal, indicating that supersonic turbulence governs the density distribution. But we also observe a deviation from the lognormal shape at AV > 1 for a subregion in Polaris that includes a prominent filament. We conclude that (1) the point where the PDF deviates from the lognormal form does not trace a universal AV -threshold for star formation, (2) statistical density fluctuations, intermittency, and magnetic fields can cause excess from the lognormal PDF at an early cloud formation stage, (3) core formation and/or global collapse of filaments and a non-isothermal gas distribution lead to a power-law tail, and (4) external compression broadens the column density PDF, consistent with numerical simulations

Anomalous peculiar motions of high-mass young stars in the Scutum spiral arm

We present trigonometric parallax and proper motion measurements toward 22 GHz water and 6.7 GHz methanol masers in 16 high-mass star-forming regions. These sources are all located in the Scutum spiral arm of the Milky Way. The observations were conducted as part of the Bar and Spiral Structure Legacy (BeSSeL) survey. A combination of 14 sources from a forthcoming study and 14 sources from the literature, we now have a sample of 44 sources in the Scutum spiral arm, covering a Galactic longitude range from 0$^\circ$ to 33$^\circ$. A group of 16 sources shows large peculiar motions of which 13 are oriented toward the inner Galaxy. A likely explanation for these high peculiar motions is the combined gravitational potential of the spiral arm and the Galactic bar.Comment: 27 pages, 52 figures, 4 tables, accepted for publication in A&

Characterizing the structure of diffuse emission in Hi-GAL maps

We present a study of the structure of the Galactic interstellar medium through the Delta-variance technique, related to the power spectrum and the fractal properties of infrared/sub-mm maps. Through this method, it is possible to provide quantitative parameters which are useful to characterize different morphological and physical conditions, and to better constrain the theoretical models. In this respect, the Herschel Infrared Galactic Plane Survey carried out at five photometric bands from 70 to 500 \mu m constitutes an unique database for applying statistical tools to a variety of regions across the Milky Way. In this paper, we derive a robust estimate of the power-law portion of the power spectrum of four contiguous 2{\deg}x2{\deg} Hi-GAL tiles located in the third Galactic quadrant (217{\deg} < l < 225{\deg}, -2{\deg} < b < 0{\deg}). The low level of confusion along the line of sight testified by CO observations makes this region an ideal case. We find very different values of the power spectrum slope from tile to tile but also from wavelength to wavelength (2 < \beta < 3), with similarities between fields attributable to components located at the same distance. Thanks to the comparison with models of turbulence, an explanation of the determined slopes in terms of the fractal geometry is also provided, and possible relations with the underlying physics are investigated. In particular, an anti-correlation between ISM fractal dimension and star formation efficiency is found for the two main distance components observed in these fields. A possible link between the fractal properties of the diffuse emission and the resulting clump mass function is discussed.Comment: Accepted by Ap