We want to understand the kinematic and thermal properties of young massive gas clumps prior to and at the earliest evolutionary stages of high-mass star formation. Do we find signatures of gravitational collapse? Do we find temperature gradients in
the vicinity or absence of infrared emission sources? Do we find coherent velocity structures toward the center of the dense and cold gas clumps? To determine kinematics and gas temperatures, we used ammonia, because it is known to be a good tracer and thermometer of dense gas. We observed the NH$_3$(1,1) and (2,2) lines within seven very young high-mass star-forming regions with the VLA and the Effelsberg 100m telescope. This allows us to study velocity structures, linewidths, and gas temperatures at high spatial resolution of 3-5$$, corresponding to $sim$0.05 pc. We find on average cold gas clumps with temperatures in the range between 10 K and 30 K. The observations do not reveal a clear correlation between infrared emission peaks and ammonia temperature peaks. We report an upper limit for the linewidth of $sim$1.3 km s$^{-1}$, at the spectral resolution limit of our VLA observation. This indicates a relatively low level of turbulence on the scale of the observations. Velocity gradients are present in almost all regions with typical velocity differences of 1 to 2 km s$^{-1}$ and gradients of 5 to 10 km s$^{-1}$ pc$^{-1}$. These velocity gradients are smooth in most cases, but there is one exceptional source (ISOSS23053), for which we find several velocity components with a steep velocity gradient toward the clump centers that is larger than 30 km s$^{-1}$ pc$^{-1}$. This steep velocity gradient is consistent with recent models of cloud collapse. Furthermore, we report a spatial correlation of ammonia and cold dust, but we also find decreasing ammonia emission close to infrared emission sources.
The chemical evolution in high-mass star-forming regions is still poorly constrained. Studying the evolution of deuterated molecules allows to differentiate between subsequent stages of high-mass star formation regions due to the strong temperature d
ependence of deuterium isotopic fractionation. We observed a sample of 59 sources including 19 infrared dark clouds, 20 high-mass protostellar objects, 11 hot molecular cores and 9 ultra-compact HII regions in the (3-2) transitions of the four deuterated molecules, DCN, DNC, DCO+ and N2D+ as well as their non-deuterated counterpart. The overall detection fraction of DCN, DNC and DCO+ is high and exceeds 50% for most of the stages. N2D+ was only detected in a few infrared dark clouds and high-mass protostellar objects. It can be related to problems in the bandpass at the frequency of the transition and to low abundances in the more evolved, warmer stages. We find median D/H ratios of ~0.02 for DCN, ~0.005 for DNC, ~0.0025 for DCO+ and ~0.02 for N2D+. While the D/H ratios of DNC, DCO+ and N2D+ decrease with time, DCN/HCN peaks at the hot molecular core stage. We only found weak correlations of the D/H ratios for N2D+ with the luminosity of the central source and the FWHM of the line, and no correlation with the H2 column density. In combination with a previously observed set of 14 other molecules (Paper I) we fitted the calculated column densities with an elaborate 1D physico-chemical model with time-dependent D-chemistry including ortho- and para-H2 states. Good overall fits to the observed data have been obtained the model. It is one of the first times that observations and modeling have been combined to derive chemically based best-fit models for the evolution of high-mass star formation including deuteration.
Context: How do molecular clouds form out of the atomic phase? And what are the relative fractions of carbon in the ionized, atomic and molecular phase? These are questions at the heart of cloud and star formation. Methods: Using multiple observatori
es from Herschel and SOFIA to APEX and the IRAM 30m telescope, we mapped the ionized, atomic and molecular carbon ([CII]@1900GHz, [CI]@492GHz and C18O(2-1)@220GHz) at high spatial resolution (12-25) in four young massive infrared dark clouds (IRDCs). Results: The three carbon phases were successfully mapped in all four regions, only in one source the [CII] line remained a non-detection. Both the molecular and atomic phases trace the dense structures well, with [CI] also tracing material at lower column densities. [CII] exhibits diverse morphologies in our sample, from compact to diffuse structures probing the cloud environment. In at least two out of the four regions, we find kinematic signatures strongly indicating that the dense gas filaments have formed out of a dynamically active and turbulent atomic/molecular cloud, potentially from converging gas flows. The atomic-to-molecular carbon gas mass ratios are low between 7% and 12% with the lowest values found toward the most quiescent region. In the three regions where [CII] is detected, its mass is always higher by a factor of a few than that of the atomic carbon. The ionized carbon emission depends as well on the radiation field, however, we also find strong [CII] emission in a region without significant external sources, indicating that other processes, e.g., energetic gas flows can contribute to the [CII] excitation as well.
The massive infrared dark cloud G0.253+0.016 projected 45pc from the Galactic centre contains ~10^5Msun of dense gas whilst being mostly devoid of observed star-formation tracers. To scrutinise the physical properties, dynamics and structure of this
cloud with reference to its star-forming potential, we have carried out a concerted SMA and IRAM 30m study of this cloud in dust continuum, CO isotopologues, shock tracing molecules, as well as H$_2$CO to trace the gas temperature. We detect and characterise the dust cores within G0.253+0.016 at ~1.3 mm and find that the kinetic temperature of the gas is >320K on size-scales of ~0.15 pc. Analysis of the position-velocity diagrams of our observed lines show broad linewidths and strong shock emission in the south of the cloud, indicating that G0.253+0.016 is colliding with another cloud at v(LSR)~70 km/s. We confirm via an analysis of the observed dynamics in the CMZ that it is an elongated structure, orientated with Sgr B2 closer to the Sun than Sgr A*, however our results suggest that the actual geometry may be more complex than an elliptical ring. We find that the column density PDF of G0.253+0.016 is log-normal with no discernible power-law tail, consistent with little star formation, and that its width can be explained in the framework of theory predicting the density structure of clouds created by supersonic, magnetised turbulence. We also present the delta-variance spectrum of this region, and show it is consistent with that expected for clouds with no star formation. Using G0.253+0.016 as a test-bed of the conditions required for star formation in a different physical environment to that of nearby clouds, we also conclude that there is not one column density threshold for star formation, but instead this value is dependant on the local physical conditions. [Abbrv.]
The mid- and far-infrared view on high-mass star formation, in particular with the results from the Herschel space observatory, has shed light on many aspects of massive star formation. However, these continuum studies lack kinematic information. W
e study the kinematics of the molecular gas in high-mass star-forming regions. We complemented the PACS and SPIRE far-infrared data of 16 high-mass star-forming regions from the Herschel key project EPoS with N2H+ molecular line data from the MOPRA and Nobeyama 45m telescope. Using the full N2H+ hyperfine structure, we produced column density, velocity, and linewidth maps. These were correlated with PACS 70micron images and PACS point sources. In addition, we searched for velocity gradients. For several regions, the data suggest that the linewidth on the scale of clumps is dominated by outflows or unresolved velocity gradients. IRDC18454 and G11.11 show two velocity components along several lines of sight. We find that all regions with a diameter larger than 1pc show either velocity gradients or fragment into independent structures with distinct velocities. The velocity profiles of three regions with a smooth gradient are consistent with gas flows along the filament, suggesting accretion flows onto the densest regions. We show that the kinematics of several regions have a significant and complex velocity structure. For three filaments, we suggest that gas flows toward the more massive clumps are present.
Understanding the chemical evolution of young (high-mass) star-forming regions is a central topic in star formation research. Chemistry is employed as a unique tool 1) to investigate the underlying physical processes and 2) to characterize the evolut
ion of the chemical composition. We observed a sample of 59 high-mass star-forming regions at different evolutionary stages varying from the early starless phase of infrared dark clouds to high-mass protostellar objects to hot molecular cores and, finally, ultra-compact HII regions at 1mm and 3mm with the IRAM 30m telescope. We determined their large-scale chemical abundances and found that the chemical composition evolves along with the evolutionary stages. On average, the molecular abundances increase with time. We modeled the chemical evolution, using a 1D physical model where density and temperature vary from stage to stage coupled with an advanced gas-grain chemical model and derived the best-fit chi^2 values of all relevant parameters. A satisfying overall agreement between observed and modeled column densities for most of the molecules was obtained. With the best-fit model we also derived a chemical age for each stage, which gives the timescales for the transformation between two consecutive stages. The best-fit chemical ages are ~10,000 years for the IRDC stage, ~60,000 years for the HMPO stage, ~40,000 years for the HMC stage, and ~10,000 years for the UCHII stage. The total chemical timescale for the entire evolutionary sequence of the high-mass star formation process is on the order of 10^5 years, which is consistent with theoretical estimates. Furthermore, based on the approach of a multiple-line survey of unresolved data, we were able to constrain an intuitive and reasonable physical and chemical model. The results of this study can be used as chemical templates for the different evolutionary stages in high-mass star formation.
Aims: We study the fragmentation and dynamical properties of a massive starless gas clump at the onset of high-mass star formation. Methods: Based on Herschel continuum data we identify a massive gas clump that remains far-infrared dark up to 100mum
wavelengths. The fragmentation and dynamical properties are investigated by means of Plateau de Bure Interferometer and Nobeyama 45m single-dish spectral line and continuum observations. Results: The massive gas reservoir fragments at spatial scales of ~18000AU in four cores. Comparing the spatial extent of this high-mass region with intermediate- to low-mass starless cores from the literature, we find that linear sizes do not vary significantly over the whole mass regime. However, the high-mass regions squeeze much more gas into these similar volumes and hence have orders of magnitude larger densities. The fragmentation properties of the presented low-to high-mass regions are consistent with gravitational instable Jeans fragmentation. Furthermore, we find multiple velocity components associated with the resolved cores. Recent radiative transfer hydrodynamic simulations of the dynamic collapse of massive gas clumps also result in multiple velocity components along the line of sight because of the clumpy structure of the regions. This result is supported by a ratio between viral and total gas mass for the whole region <1. Conclusions: This apparently still starless high-mass gas clump exhibits clear signatures of early fragmentation and dynamic collapse prior to the formation of an embedded heating source. A comparison with regions of lower mass reveals that the linear size of star-forming regions does not necessarily have to vary much for different masses, however, the mass reservoirs and gas densities are orders of magnitude enhanced for high-mass regions compared to their lower-mass siblings.
Despite increasing research in massive star formation, little is known about its earliest stages. Infrared Dark Clouds (IRDCs) are cold, dense and massive enough to harbour the sites of future high-mass star formation. But up to now, mainly small sam
ples have been observed and analysed. To understand the physical conditions during the early stages of high-mass star formation, it is necessary to learn more about the physical conditions and stability in relatively unevolved IRDCs. Thus, for characterising IRDCs studies of large samples are needed. We investigate a complete sample of 218 northern hemisphere high-contrast IRDCs using the ammonia (1,1)- and (2,2)-inversion transitions. We detected ammonia (1,1)-inversion transition lines in 109 of our IRDC candidates. Using the data we were able to study the physical conditions within the star-forming regions statistically. We compared them with the conditions in more evolved regions which have been observed in the same fashion as our sample sources. Our results show that IRDCs have, on average, rotation temperatures of 15 K, are turbulent (with line width FWHMs around 2 km s$^{-1}$), have ammonia column densities on the order of $10^{14}$ cm$^{-2}$ and molecular hydrogen column densities on the order of $10^{22}$ cm$^{-2}$. Their virial masses are between 100 and a few 1000 M$_odot$. The comparison of bulk kinetic and potential energies indicate that the sources are close to virial equilibrium. IRDCs are on average cooler and less turbulent than a comparison sample of high-mass protostellar objects, and have lower ammonia column densities. Virial parameters indicate that the majority of IRDCs are currently stable, but are expected to collapse in the future.
(Context) Many physical parameters change with time in star forming regions. Here we attempt to correlate changes in infall and outflow motions in high mass star forming regions with evolutionary stage using JCMT observations. (Aims) From a sample of
45 high mass star forming regions in three phases of evolution, we investigate the presence of established infall and outflow tracers to determine whether there are any trends attributable to the age of the source. (Methods) We obtained JCMT observations of HCO+/H13CO+ J=4-3 to trace large scale infall, and SiO J=8-7 to trace recent outflow activity. We compare the infall and outflow detections to the evolutionary stage of the host source (high mass protostellar objects, hypercompact HII regions and ultracompact HII regions). We also note that the integrated intensity of SiO varies with the full width at half maximum of the H13CO+. (Results) We find a surprising lack of SiO detections in the middle stage (Hypercompact HII regions), which may be due to an observational bias. When SiO is detected, we find that the integrated intensity of the line increases with evolutionary stage. We also note that all of the sources with infall signatures onto Ultracompact HII regions have corresponding outflow signatures as well.
The intermediate-mass star-forming core UYSO 1 has previously been found to exhibit intriguing features. While deeply embedded and previously only identified by means of its (sub-)millimeter emission, it drives two powerful, dynamically young, molecu
lar outflows. Although the process of star formation has obviously started, the chemical composition is still pristine. We present Herschel PACS and SPIRE continuum data of this presumably very young region. The now complete coverage of the spectral energy peak allows us to precisely constrain the elevated temperature of 26 - 28 K for the main bulge of gas associated with UYSO1, which is located at the interface between the hot HII region Sh 2-297 and the cold dark nebula LDN 1657A. Furthermore, the data identify cooler compact far-infrared sources of just a few solar masses, hidden in this neighbouring dark cloud.
